function data = rrgen(tmax, NoEctopic, NoNoise, seed, filename); % data = rrgen(tmax, NoEctopic, NoNoise, seed, 'filename'); % generates artificial rr time series over a time tmax seconds % and dumps the output to 'filename'. If filename is not specified, % the output is returned into the data array. If it is specified, then the % data array is zero-valued and of length(tmax). % If NoEctopic == 1 (YES) then no ectopics are added to the time series. % If NoNoise == 1 (YES) then no artifact is added to the time series. % Defaults : tmax=3600*24; NoEctopic=1; NoNoise=1; seed=1; % % Requires rrgen binary - compilation of rrgen.c on your system: % gcc -Wall rrgen.c -lm -o rrgen % % (C) P.E. McSharry & G.D. Clifford 2002 under the % GNU Public License. Contact P. McSharry (patrick@mcsharry.net) % or G. Clifford (gari@mit.edu) % Also see - McSharry P.E., Clifford G.D., Tarassenko L., % Smith L.: Method for generating an artificial RR tachogram of a typical % healthy human over 24-hours, Computers in Cardiology, 29:225-228, IEEE % Computer Society Press, September 2002. if nargin < 5 filename = 'tmp.dat'; end if nargin < 4 seed=1; end if nargin < 3 NoNoise=1; end if nargin < 2 NoEctopic=1; end if nargin < 1 tmax=3600*24; end if tmax<115 fprintf('Warning, tmax minimum value is 115, using this as default\n') end tmax=round(tmax); eval(['!rrgen ' int2str(seed) ' ' int2str(tmax) ' ' int2str(NoEctopic) ' ' int2str(NoNoise) ' > ' filename]); fid = fopen(filename,'r'); data = zeros(1,tmax); for(i=1:tmax) fscanf(fid,'%f\n',data(i)); end fclose(fid); % if not file name read in from tmp.dat and then delete it. if nargin < 5 load tmp.dat; data = tmp; if(ispc==1) system('del tmp.dat'); else system('rm tmp.dat'); end end