gen_data_ecgsyn_N.m

close all;

sfecg = 1000; % ECG sampling frequency
N_HBeats = 450; % approximate number of heart beats
Anoise = 0.0; % Additive uniformly distributed measurement noise
hrmean = 60; % Mean heart rate [bpm]
hrstd = 2; % Standard deviation of heart rate (default : 1)
lfhfratio = 0.5; % LF/HF ratio (default : 0.5)
sfint = 1000; % Internal sampling frequency

% Order of extrema: [P Q R S T]
ti = [-70 -15 0 15 100];
ai = [1.2 -5 30 -7.5 0.75];
bi = [0.25 0.1 0.1 0.1 0.4];

addpath('./ecgsyn/');

Fs = sfecg;
Lx = Fs*30;

N = 10;

% ecgsyn_N = zeros(Lx, N);
% 
% for n=1:N
%     [s_syn_init, ipeaks] = ecgsyn(sfecg,N_HBeats,Anoise,hrmean,hrstd,lfhfratio,sfint,ti,ai,bi);
%     ecgsyn_N(:, n) = s_syn_init(1:Lx);
% end

save("data_ecg_syn_N.mat", 'ecgsyn_N', 'Fs', 'hrmean', 'hrstd');