基于matlab的IIR数字滤波器的设计毕业设计.docVIP

基于matlab的IIR数字滤波器的设计 摘 要：IIR数字滤波器在MATLAB环境下的设计方法和实现方法，在无限脉冲响应（IIR）数字滤波器设计中，先进行模拟滤波器的设计，然后进行模拟—数字滤波器转换，即采用脉冲响应不变法及双线性Z变化法设计数字滤波器，最后进行滤波器的频带转换。 关键词：IIR数字滤波器；matlab；频带转换； 引言，它的平方幅度函数为： 式中,为滤波器频率，为3dB截止频率，N表示滤波器的阶次。 MATLAB信号处理工具箱为低通模拟巴特沃斯滤波器的产生提供了函数buttap,其调用的格式为：[z,p,k]=buttap(N)，其中，z表示零点,p表示极点,k表示增益，N表示阶次。 %巴特沃斯低通模拟滤波器 clear all; n=0:0.01:2; for i=1:4 switch i case 1 N=2; case 2 N=5; case 3 N=10; case 4 N=20; end [z,p,k]=buttap(N); [b,a]=zp2tf(z,p,k); [H,w]=freqs(b,a,n); magH2=(abs(H)).^2; hold on plot(w,magH2) axis([0 2 0 1]); end xlabel(w/wc); ylabel(|H(jw)|^2); title(巴特沃斯低通模拟滤波器); grid on 由图2-1可清楚地分析出，巴特沃斯滤波器拥有平滑的频率响应，在截止频率以外，频率响应单调下降。其过渡带的陡峭程度正比于滤波器的阶数，高阶巴特沃斯滤波器的频率响应近似于理想低通滤波器。由图2-1所示，对N=2阶滤波器和N=20阶滤波器的平方幅度响应进行比较后，证明了高阶巴特沃斯滤波器有着更好的幅度特性，更接近理想低通滤波器。 2.1.2 切比雪夫低通滤波器设计 巴特沃斯滤波器的频率特性在通带和阻带内部都是随频率单调变化的，如果在通带的边缘能够满足指标，那么在通带的内部会超过设计的指标要求，造成滤波器的阶次N比较高。因此，有效地设计方法是将指标的精确度均匀的分布在整个通带内，或者均匀分布在整个阻带内，或者同时分布在两者之内，这样就可用阶数较低的系统来满足要求。这可通过选择具有等波纹特性的逼近函数来达到。 切比雪夫滤波器的振幅特性就具有这种等波纹特性。它有两种形式：振幅特性在通带内是等波纹的、在阻带内是单调的切比雪夫I型滤波器；振幅特性在通带内是单调的、在阻带内是等波纹的切比雪夫II型滤波器。采用何种形式的切比雪夫滤波器取决于实际用途。 切比雪夫I型(Chebyshev-I)滤波器平方幅度响应函数为： (2-2) 其中是一个小于1的正数，它与通带波纹有关，越大，波纹也越大；是的切比雪夫多项式，它的定义为： 而切比雪夫II型(Chebyshev-II)滤波器平方幅度响应函数为： (2-4) 其中是一个小于1的正数，表示波纹变化情况；为截止频率；N为滤波器的阶次，也是的阶次。 图2-2和图2-3分别画出了理想时的切比雪夫I型与切比雪夫II型滤波器阶次N为奇数与偶数时的幅频特性。而通过MATLAB信号处理工具箱中的函数cheb1ap及cheb2ap，可以实现切比雪夫滤波器设计，其调用格式为： [z,p,k]=cheb1ap(N,Rp) [z,p,k]=cheb2ap(N,Rs) 其中，z表示零点,p表示极点,k表示增益，N表示阶次，Rp为通带波纹(dB), Rs为阻带波纹(dB)。 %切比雪夫I型低通模拟滤波器 clear all; n=0:0.01:2; for i=1:2 switch i case 1 N=7; case 2 N=8; end Rp=1; [z,p,k]=cheb1ap(N,Rp); [b,a]=zp2tf(z,p,k); [H,w]=freqs(b,a,n); magH2=(abs(H)).^2; %绘制图形 posplot=[12 num2str(i)]; subplot(posplot) plot(w,magH2) axis([0 2 0 1]); xlab