% ecg.m - ECG generator.
%
% x = ecg(L) = column vector
%
% generates piecewise linear ECG signal of length L
% must post-smooth it with an N-point smoother:
% y = sgfilt(d, N, x), usually with d=0, and N=3,5,9, etc.

function x = ecg(L)

a0 = [0,1,40,1,0,-34,118,-99,0,2,21,2,0,0,0];               % template
d0 = [0,27,59,91,131,141,163,185,195,275,307,339,357,390,440];
a = a0 / max(a0);
d = round(d0 * L / d0(15));            % scale them to fit in length L
d(15)=L;

for i=1:14,
       m = d(i) : d(i+1) - 1;
       slope = (a(i+1) - a(i)) / (d(i+1) - d(i));
       x(m+1) = a(i) + slope * (m - d(i));
end