资源描述：
《IIR数字滤波器设计及软件实现.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字系统设计与仿真实验报告学院：南湖学院系别：机械与电子工程系班级：11电信——1F姓名：学号：指导老师：实验三：IIR数字滤波器设计及软件实现1．实验目的（1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法。（2）学会根据滤波需求确定滤波器指标参数。（3）学会调用滤波器设计分析工具fdatool设计各种IIR数字滤波器。（4）掌握IIR数字滤波器的MATLAB实现方法。（5）通过观察滤波器输入输出信号的波形，建立数字滤波的概念。2．实验原理设计IIR数字滤波器一般采用间接法，即先设计一个模拟滤波器原型，再利用脉冲响应不变法或双线性变换法变换成数字滤波
2、器，其中，应用最广泛的是双线性变换法。本实验要求能根据滤波需求确定滤波器指标参数，学会按设计原理编程设计和调用滤波器设计分析工具fdatool设计各型IIR数字滤波器。本实验软件实现是调用MATLAB信号处理工具箱函数filter函数对给定输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。3.实验内容及步骤(1)IIR数字滤波器设计根据双线性变换法设计IIR数字滤波器的原理和步骤，设计一个butterworth数字低通滤波器，要求通带临界频率fp=3400Hz,阻带临界频率fs=5000Hz，通带内的最大衰减Ap=2dB，阻带内的最小衰减As=20d
3、b,采样频率Fs=22050Hz。（2)IIR数字滤波器软件实现利用第（1）步设计的数字滤波器，调用filter函数对信号进行滤波，观察滤波前后的信号波形变化，建立数字滤波的概念。4．思考题（1）简述双线性变换法设计IIR数字滤波器的设计步骤。1.学习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR数字滤波器的内容，按照教材中方法，用双线性变换法设计数字滤波器系统函数H（z）。方法一：根据教材中已求出满足本实验要求的数字滤波器系统函数： 方法二：根据设计指标，调用MATLAB信号处理工具箱函数buttord和butter，也可得到H（z）。 2.编写滤
4、波器仿真程序，计算H(z)对心电图信号采样序列x(n)的相应序列y(n)。 3.在通过计算机上运行仿真滤波程序，并调用通用绘图子程序，完成实验内容（2）和（3）。 本实验要用的MATLAB绘图函数参阅教材。（2）为了观察数字滤波的效果，应怎样设计滤波前的信号？可以用MATLAB产生一个由两个正弦波混叠的信号，然后将其输出地数据存储在一个文本文件中，再导入滤波器的测试文件中就可以了：N=256;n=1:256;x=fix(128+(2^7-10)*sin(2*pi*n/N));plot(x)fid=fopen('e:/sin.txt','wt');fprin
5、tf(fid,'%x',x);fclose(fid);5．实验报告要求（1）结合实验内容打印程序清单和信号波形。（2）对实验结果进行简单分析和解释。（3）简要回答思考题。附：IIR数字滤波器的工具设计法工具设计法是利用MATLAB提供的滤波器设计与分析工具（FDATool）进行设计的一种方法。在命令窗口输入FDATool，即可打开如图所示的图形用户界面设计窗口。该窗口分为上下两部分：上面是设计结果显示；下面用来设定所需的技术参数。FDATool需设置的参数主要有响应类型、设计方法、滤波器阶数及选项、频率参数和幅度参数等项目，不同类型和不同方法的滤波器设
6、计参数不尽相同，图中给出的是某设计实例的设计指标。设置好技术指标后，点击DesignFilter按键，即可完成设计。想对设计结果进行分析，可通过分析菜单或相关工具进行分析。如要将设计结果以传输函数的分子分母形式输出，先运行Edit/converttosinglesection命令，再运行File/Export…命令，在弹出的对话框中给变量命名为bz和az，即可在工作空间得到设计结果。%(1)IIR数字滤波器设计clear;clc;closeall;formatcompactfp=3400,Ap=2,fs=5000,As=20,Fs=22050,wp=2*p
7、i*fp/Fs,ws=2*fs/Fs,T=2,FS=1/T,Wp=tan(wp/2),Ws=tan(ws/2),N=ceil(log10((10.^(0.1*Ap)-1./(10.^(0.1*As)-1)/...(2*log10(Wp/Ws))),Wcp=Wp/((10^(.1*Ap)-1)^(1/(2*N))),Wcs=Ws/((10^(.1*As)-1)^(1/(2*N))),Wc=Wcs,[bp,ap]=butter(N,1,'s'),[bs,as]=lp2lp(bp,ap,Wc),[bz,az]=bilinear(bs,as,Fs),sys=tf(
8、bz,az,1/Fs),freqz(bz,az,512,Fs),图
