d.m

心电信号滤波

b=[0.09036,0.09036*2,0.09036];
a1=[1,-1.2686,0.7051];
a2=[1,-1.0106,0.3583];
a3=[1,-0.9044,0.2155];
[H1,f]=freqz(b,a1,25,pi);
[H2,f]=freqz(b,a2,25,pi);
[H3,f]=freqz(b,a3,25,pi);
H=H1.*H2.*H3;
plot(f,abs(H),'r');
%title('滤器频响特性曲线')