d.m

来自「心电信号滤波」· M 代码 · 共 11 行

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('滤器频响特性曲线')