IIR滤波器MATLAB设计.docVIP

IIR滤波器的MATLAB设计 姓名：王照义 学号：12010245259 班级：2010级通信工程 日期：2012年12月12日 摘 要在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号分析都是基于滤波器而进行的，(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，精度高，而且能够保留一些模拟滤波器的优良特性，因此应用很广。M软件以矩阵运算为基础，把计算可视化及程序设计有机融合到交互式工作环境中，并且为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。尤其是Matlab中的信号处理工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。本文MATLAB函数语言编程，信号处理图形界面FDATool来设计滤波器Sptool界面设计法，并用FDATool模拟IIR数字滤波器处理信号介绍最优化设计IIR 滤波器FDATool Sptool Simulink MATLAB是矩阵实验室（Matrix Laboratory）之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多.在新的版本中也加入了对C,FORTRAN,c++ ,JAVA的支持.可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用，非常的方便。lp21p,lp2hp,lp2bp,lp2bs。 运用脉冲响应不变法或双线性变法把模滤波器转数字滤波器，调用的函数是impinvar和bilinear。脉冲响应不变法适用于采样频率大于4倍截止频率的锐截止低通带通滤波器，而双线性变换法适合于相位特性要求不高的各型滤波器。 根据输出的分子分母系数，调用函数buttord计算N和wc，有系数向量可以写出数字滤波器系统函数Z,再用freqz函数验证设计结果。 IIR数字滤波器的典型设计法 设计巴特沃斯数字低通滤波器和椭圆数字低通滤波器，要求通带边界频fp=2.1kHZ，通带最大衰减Rp=0.5dB；阻带边界频率fs=8kHZ，阻带最小衰减Rs=30dB,采样频率为Fs=20kHZ。 低通巴特沃斯滤波器设计步骤如下： 确定阶数N。 取N=4 求极点 ，，，， 归一化低通原型系统函数为 由N=4直接查表得到： 极点： 归一化低通滤波器系统函数为 式中，0.0000,0.0999,0.1914,0.0252 将去归一化最终得到 通过计算可以总结出过程太麻烦，而且容易出错，结果不直观。 下面用M程序设计来实现 用脉冲响应不变法设计的巴特沃斯数字低通滤波器的M程序如下： fp=2100; fs=8000; Fs=20000; Rp=0.5; Rs=30; T=1/Fs; %设计指标 W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率 [N,Wn]=buttord(W1p,W1s,Rp,Rs,s); %确定butterworth的最小介数N和频率参数Wn [z,p,k]=buttap(N); %设计模拟低通原型的零极点增益参数 [bp,ap]=zp2tf(z,p,k); %将零极点增益转换成分子分母参数 [bs,as]=lp2lp(bp,ap,Wn*pi*Fs);%将低通原型转换为模拟低通 [bz,az]=impinvar(bs,as,Fs); %用脉冲响应不变法进行模数变换 sys=tf(bz,az,T); %给出传输函数H(Z) [H,W]=freqz(bz,az,512,Fs); %生成频率响应参数 subplot(2,1,1); plot(W,20*log10(abs(H))); %绘制幅频响应 grid on; %加坐标网格 xlabel(频率/Hz); ylabel(振幅/dB); subplot(2,1,2); plot(W,abs(H)); grid on; xlabel(频率/Hz); ylabel(振幅/H); 运行后的波形如下： 运行结果： N =4 bz