基于MATLAB对于设计IIR数字滤波器.docVIP

基于MATLAB对于设计IIR数字滤波器巴特沃斯低通滤波器，同时利用matlab软件绘制出相应的图像，这两个程序证明了这两种方法设计巴特沃斯低通滤波器的可行性。 关键词：低通滤波器，脉冲响应不变法，双线性变换法，matlab Design of MATLAB digital filter based on IIR Yin yuchen Abstract:Application of pulse response method and bilinear transformation method, a Butterworth low-pass filter is designed at the same time, at the same time, using MATLAB software to draw the corresponding image, the two procedures proved the feasibility of the two kinds of method for the design of Butterworth low pass filter. Keywords: Low pass filter, impulse response invariant method, bilinear transformation method, matlab 1 IIR数字滤波器 引言：数字滤波器是数字信号处理的重要基础。在对信号的一系列处理中，数字滤波器是极为重要的也是使用最多的线性系统。而IIR滤波器是递归型结构，通常可以写成封闭函数的形式。 正文：数字滤波器的设计一般需要包括： 1）确定性能指标 按照需求，确定滤波器 2）函数的确定 用一个因果稳定的系统函数(或差分方程、脉冲响应h(n))去逼近上述性能要求。 达到从s平面到z平面的转化。 图1 如图1所示，从s平面的阴影部分，即每一横条的左半边都会映射到z平面单位圆以内，而右半边都会映射到z平面单位圆外。同时，s平面上的虚轴都会映射到z平面的圆周上。因此映射关系并不是简单的一一对应，也不是简单的代数关系，而是 ……① 由公式①得，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为 也就是说数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。 所以只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，才可以让模拟滤波器的频率响应在折叠频率内，在数字滤波器的频率响应中得以重现，才能不产生混叠失真，即为 ， 图2脉冲响应不变法的频响混叠现象 当模拟滤波器的频率响应不严格限带，变换后就会产生周期延拓的频谱混叠，即为产生频响混叠。 因此为了减小频响混叠，只需要增加，让采样时间减少即可，则会让系统频率响应各周期延拓之间距离加大。。T=2;采样时间间隔 Fs=1/T;采样频率 wp=0.2*pi;归一化频率 ws=0.3*pi; Rp=1;最小衰减 As=15;最大衰减 [n,wc]=buttord(wp,ws,Rp,As,s);利用buttord函数确定阶数 [b,a]=butter(n,wc,s);利用butter函数完成滤波器 w=linspace(0,pi,400*pi);指定一段频率 [d,c]=impinvar(b,a,Fs);使用脉冲不变法 H=freqz(d,c,w);返回频率值 plot(w/pi,abs(H)/abs(H(1))),grid on生成图像 title(低通滤波器); xlabel(频率); ylabel(级数); 将S1平面上的横条变换到z平面的单位圆周里。这样便建立了s-s1-z平面的一一对应关系，很好的解决了混叠失真现象。 再将s1平面通过 映射到z平面。 因此，我们得到了s平面到z平面的转换关系： ….② ….③ ②③两个式子都为线性函数，所以此变换称作双线性变换。 但是一般我们为了模拟滤波器和数字滤波器有相应的关系，引入了常数c，是上述表达式变成了： …..④ 将代入后，可以得到 同时由上述公式可得双线性变换的频率变换关系，如下图所示 图4频率变换关系 2）常数c的选择 ⒈当模拟滤波器与数字滤波器，在低频存在对应关系时 同时由公式④得： 所以 ⒉当数字滤波器某频率和模拟滤波器的严格对应的时候 所以 3）优缺点 优点：双线性变换法完美的解决了频率响应而引起的混叠现象 缺点：由于s-s1-z的这个变换是通过严重非线性关系得到的，所以当一个线性相位的模拟滤波器经过双线性变换之后，不会再保持原有的线性相位了。 结束语：显然，双线性变换法结局了脉冲响应不变法的混叠现象，但是同时自身也出现了线性相位不再保持问题。 代码： T=2; %采样频率 %归一化频率wp=0.2*pi; ws=0.3*pi;