bp_display.asv

来自「电子血压计的显示程序,用MATLAB编写」· ASV 代码 · 共 46 行

%---bp_display

% data-file
DATAFILE='ttt.dcm';         

% lowpass filters parameter
%
N = 256;
fs = 133;
n = 8;
%wn = 10/(fs/2);
wn = 1.2/(fs/2);
rp = 3;
rs = 60;
%

fid=fopen(DATAFILE,'r');
data=fread(fid,5001,'int16');
for i=1:2500
data1(i)=data(2*i);
data2(i)=data(2*i+1);
end
[b,a]=butter(n,wn);
%data10=filter(b,a,data1);
data10=filtfilt(b,a,data1);
subplot(511),plot(data1);
subplot(512),plot(data10);
%data20=filter(b,a,data2);
data20=filtfilt(b,a,data2);
subplot(513),plot(data2);
subplot(514),plot(data20);
t=[1:2500];
[up,down]=envelope(t,data20,'linear');

%[up,down]=envelope(t,data2,'linear');
subplot(515),plot(t,data20,'g-');hold on;
%subplot(515),plot(t,data2,'g-');hold on;
subplot(515),plot(t,up,'r-.');
subplot(515),plot(t,down,'r-.');
t=0;
for i=1:2500
    if data(i)>t
        t=data(i);
    end
end
hold off;