(最新整理)数字信号处理实验报告实验三

1、(完整)数字信号处理实验报告实验三(完整)数字信号处理实验报告实验三 编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)数字信号处理实验报告实验三）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快 业绩进步，以下为(完整)数字信号处理实验报告实验三的全部内容。数字信号处理实验报告实验三：iir滤波器设计实验班 级: 13080241 姓

2、 名： 徐可 学 号：1308024121 实验三 iir滤波器设计实验一、实验目的1。学习模拟数字变换滤波器的设计方法2.掌握双线性变换滤波器的设计方法3.掌握实现数字滤波的具体方法。二、实验要求 1。 用双线性变换法设计一个巴特沃斯低通iir数字滤波器。设计指标参数为：在通带内频率低于0。2时，最大衰减小于1db；在阻带内0。3, 频率区间上，最小衰减大于15db.2。0.02为采样间隔，打印出数字滤波器在频率区间0， /2上的频率响应特性曲线。3。 用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，观察总结滤波作用与效果附：心电图采样序列x（n）人体心电图信号在测量过程中往往受到

3、工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息.下面给出一实际心电图信号采样序列样本x(n),其中存在高频干扰.在实验中以x（n)作为输入序列,滤除其中的干扰成分.x（n)=4，2，0，-4，6，4，-2，-4,-6，6,4,-4,-6，6,2，6,12,8，0，16,38，-60，84，90,66，-32，4,2，4,8,12,12,10，6,6，6，4，0,0,0，0,0，2,4,0，0,0,-2，-2,0,0，2,-2,-2,-2,0三、实验设备计算机、matlab软件四、实验报告要求1。给出详细的滤波器设计说明书；2。整理好经过运行并证明是正确的程序并且加上详

4、细注释；3。用所设计的滤波器对心电信号进行滤波，打印滤波后的波形;五、实验结果及分析iir滤波器设计步骤：1. 将数字滤波器的技术指标转换为巴特沃斯模拟低通滤波器的技术指标.2. 算出巴特沃斯低通滤波器的阶数及3. 确定归一化模拟滤波器模型的零极点形式，确定。4. 用双线性变换法将转换为。根据双线性法设计滤波器的步骤实现。ap=1；as=15； %数字滤波器指标t=0.02pi; %采样频率与间隔fsa=1/t；wp=0。2pi;ws=0。3pi； 转换为数字角频率wp=2/ttan（wp/2);ws=2/ttan（ws/2); %由数字角频率转换为模拟角频率n，wc=buttord(wp,w

5、s，ap,as，s); %获取模拟滤波器的阶数和3db截止频率 z,p，k=buttap（n)； %归一化模拟滤波器模型的零极点形式参数b,a=zp2tf（z，p，k）； %归一化模拟滤波器传递函数的系数bl，al=lp2lp（b,a，wc）； 把模拟滤波器原型转换成截至频率为wc的低通滤波器b,a=bilinear(bl,al，fsa)；%用双线性变换法实现模拟滤波器到数字滤波器的转换 h，w=freqz(b，a)； %获取频率响应l=length（w）/2+1；figure（1)，plot(w（1：l)/pi，abs（h（1：l),grid，xlabel（角频率（pi）),ylabel(频

6、率响应幅度);plot(wfsa/（2pi）,abs（h);grid； %绘制频率响应曲线%xlabel（频率（hz)）;ylabel（频率响应幅度）；2.消除心电图中干扰成分ap=1;as=15; %数字滤波器指标t=0。02pi; %采样频率与间隔fsa=1/t；wp=0。2pi；ws=0.3*pi； %转换为数字角频率wp=2/ttan(wp/2)；ws=2/t*tan(ws/2)； 由数字角频率转换为模拟角频率n，wc=buttord(wp，ws，ap,as,s)； %获取模拟滤波器的阶数和3db截止频率 z,p,k=buttap（n）; %归一化模拟滤波器模型的零极点形式参数b，a=

7、zp2tf（z,p,k）； 归一化模拟滤波器传递函数的系数bl,al=lp2lp（b，a,wc）； %把模拟滤波器原型转换成截至频率为wc的低通滤波器b,a=bilinear（bl，al,fsa)；%用双线性变换法实现模拟滤波器到数字滤波器的转换 h，w=freqz（b，a）; 获取频率响应l=length(w）/2+1；figure(1)，plot(w（1:l）/pi，abs(h(1：l)）,grid,xlabel(角频率（pi），ylabel(频率响应幅度)；plot(wfsa/（2*pi），abs(h）)；grid； %绘制频率响应曲线xlabel(频率(hz）;ylabel(频率响应幅度）;x=4，-2,0,4,-6，4，-2，-4,-6，6,-4，4,-6，-6,2,6，12,8,0，-16,-38，-60,-84，90,66，-32,-4,2，4,8，12,12,10,6,6，6，4,0,0，0，0,0,2，-4，0,0,0，-2，-2，0，0,2,-2，2,2,0；y=filter(b,a，x);figure(2),subplot(2，1，1),plot(x），title（原心电图信号);subplot（2,1，2）,plot（y)，title(滤波后的心电图信号）；六、实验体会1消除心电图中的干扰成分我们从图中看