数字信号处理实验报告

1、数字信号处理报告IIR数字滤波器上海理工大学教师：苏湛组员：王世豪 徐骞 刘新 2016.1.4一、实验简介Butterworth和Chebyshev低通滤波器方法：1) 根据性能参数，先设计一个模拟滤波器，按照一定的算法转换为满足预定指标的数字滤波器。利用模拟原型滤波器的逼近算法和特性。2）计算机辅助设计，从统计概念出发，对所要提取的有用信号从时域进行估计，在统计指标最优的意义下，使得估计值最优逼近有用信号，减弱或消除噪声。1）Butterworth低通滤波器1 幅频特性：，其中N为滤波器的阶数，为通带截止频率。在=0处，有最大值；2）在通带截止频率处，不同阶次的幅频量值都相同，即为；3）阶

2、数N增加时，通带幅频特性变平，阻带衰减更快，逐渐趋近于理想滤波器的幅频特性。幅频特性通常用衰减函数描述。分贝（dB）2 极点一共有2N个，并且以圆点为对称中心成对的出现。 k=1,2,N 系统函数：3 通带衰减函数、阻带衰减函数和系统幅频特性的关系： 为通带截止频率 为阻带截止频率4 阶数N 5 通带截止频率 确定了滤波器的阶数N和通带截止频率，就可以求出系统的极点，从而求出系统函数，这样就完成了Butterworth低通滤波器的设计。通常这是在给定技术指标的前提下进行的。例题： 设计一个Butterworth低通滤波器，要求频率小于20rad/s范围内幅频响应衰减不大于2dB，频率大于30r

3、ad/s的幅频响应的衰减不小于10dB。解：该滤波器的技术指标为：=20rad/s, =2dB, =30rad/s, =10dB代入阶数N的计算公式可得：=3.371取满足以上条件的最小整数N=4。2）Chebyshev低通滤波器Chebyshev滤波器最主要的特点是引入了Chebyshev多项式。这是其特殊幅频特性的数学基础。1 Chebyshev多项式：，为信号的模拟角频率。由此多项式可以得出如下特性：1）时，在-1和1之间波动；2）=1时，=1；3）=0时，若N为奇数，则=0；若N为偶数，则等于1或-1；4）时，随单调增大，N越大，得增幅越大。2 幅频特性：3 系统函数：=4 极点分布：

4、2N个极点成对分布在椭圆的圆周上。 ， 5 通带波动函数：6 波纹系数：7 阶数：例题：设计一个Chebyshev低通滤波器，技术指标为：通带波动=1dB，截止频率，阻带衰减函数，阻带边界频率为。解：对频率进行归一化处理 波纹系数 ，代入N的求解公式，求得取N=4。2）概念1.脉冲响应不变法：脉冲响应不变法是一种将模拟滤波器转化为数字滤波器的基本方法。它利用模拟滤波器理论设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n) 模仿模拟滤波器的冲击响应ha(t), 使h(n)正好等于ha

5、(t)的采样值，即h(n)=ha(nT)T为采样周期。3）基本思想2.双线性变换法：在将模拟滤波器H(s)转换为数字滤波器H(z)时，不是直接从s域到z域，而实现将非带限的H(s)映射为代限的H(s),再通过脉冲响应不变法将s域映射到z域，即H(s)->H(s)->H(z).4）对比3.脉冲响应不变法及双线性变换法优缺点脉冲响应不变法的优点：1，模拟频率到数字频率的转换时线性的；2，数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应，因此时域特性逼近好。 缺点：会产生频谱混叠现象，只适合带限滤波器双线性变换法优点：克服多值映射得关系，可以消除频率的混叠 缺点：时域到频

6、域的变换是非线性的，在高频处有较大的失真5）IIR数字滤波器的特点1、封闭函数IIR数字滤波器的系统函数可以写成封闭函数的形式。2、IIR数字滤波器采用递归型结构IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。3、借助成熟的模拟滤波器的成果IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不

7、高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。4、需加相位校准网络IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。6）Matlab运行结果如图二、实验设计1. 实验讨论根据所要设计滤波器的参数去确定一个模拟滤波器的传输函数，然后再根据这个传输函数，通过双线性变换、或脉冲响应不变法来进行数字滤波器的设计。它的设计比较复杂，复杂在于它的模拟滤波器传输函数H（s）的确定。这一点我们可以让软件来实现。然后，我们说一下它的具体实现步骤：首先你要先确定你需要一个什么样的滤波器，巴特沃斯型，切比雪夫

8、型，还是其它什么型的滤波器。当你选定一个型号后，你就可以根据设计参数和这个滤波器的计算公式来确定其阶数、传输函数的表达式。通常这个过程中还存在预扭曲的问题（这只是双线性变换法所需要注意的问题，脉冲响应不变法不存在这种问题）。确定H（S）后，就可以通过双线性变换得到其数字域的差分方程。2. 实验目的在理论学习的基础上，掌握不同IIR滤波器的性质、特点，并通过实验学习如何设计各种常用的IIR滤波器，以便在实际工作中能根据具体情况使用IIR滤波器。3. 设计流程模拟原型低通 模拟带通滤波器 数字带通滤波器4.实验原理根据模拟带通滤波器的上、下截频确定变换式中的参数B和w0B=wpu-wpl=0.31

9、42;w02=wpl*wpu=1.5629;确定原型低通滤波器的通带截频wp、阻带截频wswp1=(wpu2-w02)/(wpu-wpl)*wpu)=1.0000wp2=(wpl2-w02)/(wpu-wpl)*wpl)=-1.0000ws1=(wsu2-w02)/(wpu-wpl)*wsu)=1.8750ws2=(wsl2-w02)/(wpu-wpl)*wsl)=-2.1875原型低通滤波器的通带截频wp=1rad/s,但变换后的阻带截频|ws1|和|ws2|通常不等，只有选择其中较小者作为原型滤波器的阻带截频ws,相应的模拟带通滤波器在阻带衰减才能满足设计，即ws=min|ws1|ws2|

10、=1.8750设计技术指标wp=1rad/s,Ap=3dB,ws=1.8750rad/s,As=10dB的原型低通滤波器HL (s)Nlg(1000.1As-1)/(1000.1Ap-1)/2lg(ws)=1.0936N=2，Butterworth低通滤波器的wc为wc=1/(100.1Ap-1)1/2N=1.0012由N和wc可得所设计的原型低通滤波器函数HL(s)= 1.0024-s2 + 1.416 s + 1.0024利用复频率变换获得模拟带通滤波器的系统函数为HBS(s)= 0.09893 s2 + 1.237e-016 s - 4.437e-017-s4 + 0.4448 s3 +

11、 4.984 s2 + 1.087 s + 5.967利用脉冲响应不变法将模拟滤波器转换为数字滤波器HBS(z)=6.939e-018 z3 + 0.02237 z2 - 0.07855 z + 0.02793-z4 - 0.03108 z3 + 1.562 z2 - 0.01985 z + 0.6409三、程序代码clear; wpu=0.55*pi;wpl=0.45*pi; wsu=0.6*pi;wsl=0.4*pi;B=0.55*pi-0.45*pi;w0=(wpl*wpu)0.5; %根据模拟带通滤波器的上、下截频确定变换式中的参数B和w0.wp1=abs(wpu2-w02)/(wpu

12、-wpl)*wpu)wp2=abs(wpl2-w02)/(wpu-wpl)*wpl)ws1=abs(wsu2-w02)/(wpu-wpl)*wsu)ws2=abs(wsl2-w02)/(wpu-wpl)*wsl)wp=min(wp1,wp2)ws=min(ws1,ws2)%确定原型低通滤波器的通带截频wp、阻带截频ws./通带截频一般相等，阻带截频一般不等。选择其中较小者作为原型低通滤波器的阻带截频。 Ap=3; As=10; Fs=1; n1=buttord(wp,ws,Ap,As,'s') %确定阶数 wn1=wp/(10(0.1*Ap)-1)(1/2/n1); num,d

13、en=butter(n1,wn1,'s'); %确定BW AFw=0:pi/512:pi; h=freqz(numd,dend,w); norm=max(abs(h);%确定归一化DF的幅度响应 numd=numd/norm; g=20*log10(abs(h)/norm); plot(w/pi,g); w1=wpl wpu wsl wsu; h=freqz(numd,dend,w1) fprintf('Ap1=%.4fn',-20*log10(abs(h(1); fprintf('Ap2=%.4fn',-20*log10(abs(h(2); fprintf('As1=%.4fn',-20*log10(abs(h(3); fprintf('As1=%.4fn',-20*log10(abs(h(4);四、实验结果1.实验数据wp1 =1.0000wp2