数字信号处理实验七

云南大学信息学院数字信号处理实验报告1实验七数字滤波器设计实验室名称：信息学院2204实验时间：2015年11月26日姓名：蒋逸恒学号：20131120038专业：通信工程指导教师：陶大鹏成绩教师签名：年月日一、实验目的1、学会使用matlab设计符合要求的数字滤波器的参数。2、学会使用基本的matlab命令，把获得的数字滤波器参数带到传输函数中写出传输函数，之后利用之前的实验知识，画出滤波器级联框图，实现滤波器。二、实验内容Q7.1用MATLAB确定一个数字无限冲击响应低通滤波器所有四种类型的最低阶数。指标如下：40kHz的抽样率，4kHz的通带边界频率，8kHz的阻带边界频率，0.5dB的通带波纹，40dB的最小阻带衰减。评论你的结果。Q7.2用MATLAB确定一个数字无限冲击响应高通滤波器所有四种类型的最低阶数。指标如下：3500Hz的抽样率，1050Hz的通带边界频率，600Hz的阻带边界频率，1dB的通带波纹，50dB的最小阻带衰减。评论你的结果。Q7.5通过运行程序P7.1来设计巴特沃兹带阻滤波器。写出所产生的传输函数的准确表达式。滤波器的指标是什么？你的设计符合指标吗？使用MATLAB，计算并绘制滤波器的未畸变的相位响应及群延迟响应。Q7.6修改程序P7.1来设计符合习题Q7.1所给指标的切比雪夫1型低通滤波器。写出所产生的传输函数的准确表达式。你的设计符合指标吗？使用MATLAB，计算并绘制滤波器的未畸变的相位响应及群延迟响应。Q7.20使用函数fir1，设计一个线性相位有限冲击响应低通滤波器，使其满足习题Q7.13给出的指标，并画出其增益和相位响应。使用习题Q7.13中用凯泽公式估计出的阶数。用表格形式显示滤波器系数。你的设计满足指标吗？若不满足，调整滤波器阶数直到设计满足指标。满足指标的滤波器阶数是多少？Q7.23用凯泽窗设计一个有限冲击响应低通滤波器。滤波器的指标是：wp=0.31,ws=0.41,As=50dB。注意，函数kaiser需要参数β及阶数N的值，他们必须先用式（7.36）和式（7.37）分别算出。你的设计满足指标吗？Q7.25用fir2设计一个95阶有限冲击响应滤波器，它具有三个不同的常熟幅度级：在频率范围0到0.25中为0.4，在频率范围0.3到0.45中为1.0，在频率范围0.5到1.0中为0.8。画出所设计的滤波器的幅度响应。你的设计满足指标吗？Q7.27用remez设计具有如下指标的有限冲击响应带通滤波器：通带边界为云南大学信息学院数字信号处理实验报告21.8kHz和3.0kHz，阻带边界为1.5kHz和4.2kHz，通带波纹δp=0.1，阻带波纹δs=0.02，抽样率为12kHz。用kaiserord估计滤波器的阶数。你的设计是一个最优有限冲击响应滤波器吗？你的设计满足指标吗？若不满足，增加滤波器阶数在满足指标方面有用吗？指标由一个较低阶数的滤波器来满足而不是由kaiserord得到的来满足吗？在不等过渡带的情形下，用设计的滤波器可能在较大的过度带宽中以增益响应表现不满意的行为。改进该行为的一种方法是：通过移动阻带边界减少过渡带宽，直到使设计在过渡带中以平滑的下降来满足指标。在通带边界保持固定的情况下，尝试这种方法并确定新的指标，它在过渡带中提供平滑的下降。三、实验器材及软件1.微型计算机1台2.MATLAB7.0软件四、实验原理IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。相对的，FIR滤波器是直接型结构，没有递归结构，也就没有反馈环路，FIR滤波器运算结构通常也是由由延时、乘以系数和相加等基本运算组成，可以组合成各种结构形式的滤波器组来实现滤波。IIR滤波器的单位脉冲响应为无限长，网络中有反馈回路。FIR（FiniteImpulseResponse）滤波器的单位脉冲响应是有限长的，一般网络中没有反馈回路。IIR滤波器的系统函数一般是一个有理分式，分母多项式决定滤波器的反馈网络。FIR滤波器的系统函数一般是一个分母为1的多项式，所以没有反馈网络。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。FIR的设计在matlab中其实也同样很简单，其主要思想是通过先设计出无限长的IIR滤波器再通过截尾的方式获得一个FIR滤波器，但是简单的截尾会出现吉布斯现象，所以此时还需要对截尾之后的结果加窗；来减小吉布斯震荡现象。IIR数字滤波器幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上；FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变，这是很好的性质。FIR数字滤波器是有限的单位响应也有利于对数字信号的处理，便于编程。五、实验步骤1、进行本实验，首先必须熟悉matlab的运用，所以第一步是学会使用matlab。2、学习相关基础知识，根据《数字信号处理》课程的学习理解实验内容和目的。3、在充分熟悉基础知识的情况下进行实验，利用matlab完成各种简单的波形产生和云南大学信息学院数字信号处理实验报告3观察，理解各种波形产生的原理和方法。4、从产生的图形中学习新的知识，掌握实验的目的，充分学习数字滤波器的设计和运用。5、最后需要思考各种滤波器的联系和区别，以及如何根据要求设计出最符合要求且最容易实现的滤波器6、小结本实验，通过小结发现问题并解决问题。六、实验记录（数据、图表、波形、程序等）7.1输入：[N1,wn1]=buttord(0.2,0.4,0.5,40)[N2,wn2]=cheb1ord(0.2,0.4,0.5,40)[N3,wn3]=cheb2ord(0.2,0.4,0.5,40)[N4,wn4]=ellipord(0.2,0.4,0.5,40)输出：N1=8wn1=0.2469N2=5wn2=0.2000N3=5wn3=0.4000N4=4wn4=0.20007.5输入：Ws=[0.40.6];Wp=[0.20.8];Rp=0.4;Rs=50;[N1,Wn1]=buttord(Wp,Ws,Rp,Rs);[num,den]=butter(N1,Wn1,'stop');disp('分子系数：');disp(num);disp('分母系数：');disp(den);[g,w]=gain(num,den);plot(w/pi,g);gridaxis([01-605]);xlabel('\omega/\pi');ylabel('增益，dB');title('巴特沃斯带阻滤波器的增益响应');输出：分子系数：0.04930.00000.24650.00000.49300.00000.49300.00000.24650.00000.0493分母系数：1.00000.0000-0.08500.00000.63600.0000-0.02880.00000.05610.0000-0.000800.10.20.30.40.50.60.70.80.91-60-50-40-30-20-100/增益，dB巴特沃斯带阻滤波器的增益响应云南大学信息学院数字信号处理实验报告4相位响应和群延迟响应：w=0:pi/255:pi;Ws=[0.40.6];Wp=[0.20.8];Rp=0.4;Rs=50;[N1,Wn1]=buttord(Wp,Ws,Rp,Rs);[num,den]=butter(N1,Wn1,'stop');h=freqz(num,den,w);subplot(1,2,1);plot(w/pi,unwrap(angle(h)));gridon;k=[diff(angle(h))./diff(w)0];subplot(1,2,2);plot(w/pi,k);axis([0,1,3,10]);gridon;7.6输入：Ws=0.4;Wp=0.2;Rp=0.5;Rs=40;[N1,Wn1]=buttord(Wp,Ws,Rp,Rs);[num,den]=butter(N1,Wn1);disp('分子系数：');disp(num);disp('分母系数：');disp(den);[g,w]=gain(num,den);plot(w/pi,g);gridaxis([01-605]);xlabel('\omega/\pi');ylabel('增益，dB');title('切比雪夫1型低通滤波器的增益响应');输出：分子系数：0.00040.00200.00400.00400.00200.0004分母系数：1.0000-3.82696.2742-5.44642.4915-0.479700.10.20.30.40.50.60.70.80.91-60-50-40-30-20-100/增益，dB切比雪夫1型低通滤波器的增益响应云南大学信息学院数字信号处理实验报告5相位响应和群延迟响应：Ws=0.4;Wp=0.2;Rp=0.5;Rs=40;%估计滤波器阶数[N1,Wn1]=cheb1ord(Wp,Ws,Rp,Rs);%设计滤波器[num,den]=cheby1(N1,Rp,Wn1);h=freqz(num,den,w);subplot(1,2,1);plot(w/pi,unwrap(angle(h)));gridon;k=[diff(angle(h))./diff(w)0];subplot(1,2,2);plot(w/pi,k);axis([0,1,-5,20]);gridon;7.20输入：w=0:pi/255:pi;Fs=2500;Fp=2000;Rp=0.005;Rs=0.005;Ft=10000;N=kaiord(Fp,Fs,Rp,Rs,Ft);b=fir1(N,2*Fp/Ft);[g,w]=gain(b,1);subplot(2,1,1);plot(w/pi,g);grid;axis([01-1255])xlabel('\omega/\pi');ylabel('增益，dB');title('FIR低通滤波器的增益响应');h=freqz(b,1,w);subplot(2,1,2);plot(w/pi,unwrap(angle(h)));gridon;xlabel('\omega/\pi');ylabel('相位');title('FIR低通滤波器的相频响应');00.10.20.30.40.50.60.70.80.91-100-500X:0.4Y:-6.017/增益，dBFIR低通滤波器的增益响应00.10.20.30.40.50.60.70.80.91-40-30-20-100/相位FIR低通滤波器的解缠绕的相位响应云南大学信息学院数字信号处理实验报告6输出：b=-0.00070.00070.0014-0.0000-0.0023-0.00190.00250.0052-0.0000-0.0083-0.00640.00790.0157-0.0000-0.0233-0.01760.02150.0431-0.0000-0.0706-0.06000.09190.30130.39980.30130.0919-0.0600-0.0706-0.00000.04310.0215-0.0176-0.0233-0.00000.01570.0079-0.0064-0.0083-0.00000.00520.0025-0.0019-0.0023-0.00000.00140.0007-0.00077.23输入：w2=0:pi/255:pi;Wp=0.31;Ws=0.41;Wn=Wp