process.m

计算ECG压缩程式

%
clc;
m=10;
BASE=250;
LIMIT=256;
ecgX = ECG_ORG(1, BASE, LIMIT);
ecgY = ECG_ORG(1, BASE+16, LIMIT);
%ecgY = ECG_ORG(2, BASE, LIMIT);
figure(1)
subplot(211)
plot(ecgX(1:length(ecgX), 2), 'r')
hold on
plot(ecgY(1:length(ecgY), 2), 'g')
hold off
axis([1 LIMIT -1.0 1.0])
%
ampX = ecgX(:,2);
ampY = ecgY(:,2);
ampX = ECG_ORG_single(1,LIMIT);
ampY = ECG_ORG_single(1,LIMIT);
subplot(212)
plot(ampX, 'r-')
hold on
plot(ampY, 'g-')
hold off
axis([1 LIMIT -1.0 1.0])
%
ampZ = ampX+i*ampY;
FFTz = fft(ampZ);
REALz = real(ampZ);
IMAGz = imag(ampZ);
figure(2)
subplot(211)
plot(REALz, 'r-');
hold on
plot(IMAGz, 'g-');
hold off
axis([1 LIMIT min(min(REALz), min(IMAGz)) max(max(REALz), max(IMAGz))])
subplot(212)
stem(abs(FFTz), 'k');
axis([1 LIMIT min(abs(FFTz)) max(abs(FFTz))])
%
FFTz_pro = FFTz;
%FFTz_pro(241:1:256) = FFTz_pro(240:-1:225);
IFFTz_pro=ifft(FFTz_pro);
REALz_pro=real(IFFTz_pro);
IMAGz_pro=imag(IFFTz_pro);
figure(3)
subplot(211)
stem(abs(FFTz_pro), 'k')
axis([1 LIMIT min(abs(FFTz_pro)) max(abs(FFTz_pro))])
subplot(212)
plot(REALz_pro, 'r-')
hold on
plot(IMAGz_pro, 'g-')
hold off
axis([1 LIMIT min(min(REALz_pro), min(IMAGz_pro)) max(max(REALz_pro), max(IMAGz_pro))])
%