example.asv

电子血压计的显示程序,用MATLAB编写

datafile='62-33.dcm';
fs=133;
n=8;
wn=1.2/(fs/2);
fid=fopen(datafile,'r');
data=fread(fid,5001,'int16');
for i=1:1250
    data1(i)=data(2*i);
    data2(i)=data(2*i+1);
end
[b,a]=butter(n,wn);
for i=1:max(length(a),length(b))-1
    z(i)=0
end
[data3,t]=filter(b,a,data2(1:1250),z);
for i=1250:2500
    data4(i)=data(2*i-1);
    data2(i)=data4(i);
end
[b,a]=butter(n,wn)
[data5,g]=filter(b,a,data4(1250:2500),t);
for i=1:1250
    data4(i)=data3(i);
    data4(i+1250)=data5(i);
end
%n=5;
%p=2;
%[c,l]=wavedec(data2,5,'db4');
%nc=wthcoef2('t',c,l,n,p,'h');
%data3=waverec(nc,1,'db4');
subplot(411),plot(data2(1:1250));
subplot(412),plot(data3);
subplot(413),plot(data2(1250:2500));
subplot(414),plot(data4);
%subplot(313),plot(t);