example.m

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

clear
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,z);
for i=1250:2500
    data4(i)=data(2*i-1);
    data2(i)=data4(i);
end
for i=1:1250
    data5(i)=data4(i+1249);
end
[b,a]=butter(n,wn)
[data6,h]=filter(b,a,data5,z);
[b,a]=butter(n,wn)
[data7,g]=filter(b,a,data5,t);
for i=1:1250
    data4(i)=data3(i);
    data4(i+1250)=data7(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');
figure(1)
subplot(511),plot(data2(1:1250));
subplot(512),plot(data3);
subplot(513),plot(data5);
subplot(514),plot(data6);
subplot(515),plot(data4);
figure(2)
subplot(311),plot(data3);
subplot(312),plot(data6);
subplot(313),plot(data4);
hold on
subplot(313),plot(1000,data4(1000),'go');
%subplot(313),plot(t);