iir.m

一个简单的IIR滤波器实现编程

%H(z)分子，分母的系数向量n,m
n=[0.0026 0.0155 0.0388 0.0517 0.0388 0.0155 0.0026];
m=[1.0000 -2.3797 2.9104 -2.0551 0.8779 -0.2099 0.0218];
[h,w]=freqz(n,m,512,200);
figure(1);
subplot(221);
plot(w,abs(h));
title('IIR滤波器');
%求输出信号
load BBB.TXT;
subplot(223);
plot(BBB);
%y=linspace(0,0,512);
for i=1:7
 y(i)=0;
end
title('原心电信号图');
N=length(BBB);
for i=7:N
   y(i)=2.3797*y(i-1)-2.9104*y(i-2)+2.0551*y(i-3)-0.8779*y(i-4)+0.2099*y(i-5)-0.0218*y(i-6)+0.0026*BBB(i)+0.0155*BBB(i-1)+0.0388*BBB(i-2)+0.0517*BBB(i-3)+0.0388*BBB(i-4)+0.0155*BBB(i-5)+0.0026*BBB(i-6);
end
subplot(224);
plot(y);
title('滤波后的心电图');