《数字滤波器实现介绍(共13页).doc》由会员分享，可在线阅读，更多相关《数字滤波器实现介绍(共13页).doc（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上实验八 数字滤波器实现实验室名称： 信息学院2204 实验时间：2015年11月26日姓 名 一、实验目的1、通过逼近给定的频率响应指标而得到数字传输函数G（z），并通过软件实现。2、利用MATLAB来研究冲激响应或有限冲激响应数字滤波器结构的软件实现。3、并联型实现、并联II型实现、直接II型仿真滤波器的验证仿真。二、实验内容Q8.1程序P8.1设计了什么类型的滤波器？其指标是什么？滤波器的阶数是多少？为了检验仿真，需要计算多少个冲激响应样本？仿真是正确的吗？Q8.2用转置的直接II型仿真滤波器，修改程序P8.1并运行修改后的程序。仿真是正确的吗？Q8.3生成Q8.
2、1中产生的传输函数的一个级联实现，并编写一个程序来仿真它。其中每个单独的部分用一个直接型实现。验证仿真。Q8.5 生成Q8.1中传输函数的一个并联型实现，并编写一个程序来仿真它。每一个单独的部分用直接型实现。验证仿真。Q8.6 生成Q8.1中传输函数的一个并联II型实现，并编写一个程序来仿真它。每一个单独的部分用直接型实现。验证仿真。Q8.7运行程序P8.2并产生被仿真的级联格型结构的传输函数。该传输函数与式（8.35）中的相同吗？Q8.9程序设计了什么类型的滤波器？其指标是什么？滤波器的阶数是多少？形成输入的正弦序列的频率是多少？Q8.10 运行P8.3并产生两个图形。哪种输入成分会在滤波器
3、输出出现？为什么输出序列的开始部分不是一种理想的正弦曲线？修改P8.3，以便纸过滤序列X2n。产生的输出序列和预料的是一样的吗？证明你的答案。Q8.14 程序P8.4设计了什么类型的滤波器？其指标是什么？滤波器的阶数是多少？为了验证仿真，需要多少个冲激响应样本？仿真是正确的吗？Q8.15修改程序P8.4，用直接II型仿真该滤波器，运行修改后的程序。仿真是正确的吗？Q8.16使用程序P8.4生成在习题Q8.4中产生的传输函数的一个级联实现，写出对其仿真的一个程序。其中每一个部分以直接II型实现。验证仿真。三、实验器材及软件1. 微型计算机1台2. MATLAB 7.0软件四、实验原理1、可调数字
4、滤波器的设计基础是复变量的z谱变换，通过z替换为F(z)，将具有给定截止频率的数字滤波器的实现调整为具有不同截止频率的另一个实现。2、利用全通结构代替每一个延时块，可能会导致实现Gnew(z)时具有不能被实现的无延时环结构。3、设计可调有限冲激响应滤波器的一种简单方法是基于加窗的傅里叶级数方法。4、通常情况下，复数乘法和复数加法总数量是N（log2N），通使用Wn的对称性质，可进一步简化计算得复杂性。五、实验步骤1、 打开计算机中的MATLAB2、 在MATLAB中新建一个M-file，在其中输入实验要求的相关代码3、 输入代码后调试代码，代码无误后可以得到相应的实验结果4、 观察实验结果回答
5、相关问题5、记录实验结果，将图形等截图粘贴，完成实验报告六、实验记录（数据、图表、波形、程序等）Q8.1% Program P8_1clear all;clc;Wp = 0.4 0.5; Ws = 0.1 0.8; Rp = 1; Rs = 30;N1, Wn1 = buttord(Wp, Ws, Rp, Rs)num,den = butter(N1,Wn1);disp(分子系数是 );disp(num);disp(分母系数是 );disp(den);impres = direct2(num,den,1 zeros(1,4*N1);p,d = strucver(impres,2*N1);dis
6、p(实际分子系数是 ); disp(p);disp(实际分母系数是 ); disp(d);运行结果：N1 = 2 Wn1 = 0.3608 0.5421分子系数是 0.0571 0 -0.1143 0 0.0571分母系数是 1.0000 -0.5099 1.2862 -0.3350 0.4479实际分子系数是 0.0571 -0.0000 -0.1143 0.0000 0.0571实际分母系数是 1.0000 -0.5099 1.2862 -0.3350 0.4479Q8.2% Program P8_1clear all;clc;Wp = 0.4 0.5; Ws = 0.1 0.8; Rp
7、= 1; Rs = 30;N1, Wn1 = buttord(Wp, Ws, Rp, Rs)num,den = butter(N1,Wn1);disp(分子系数是 );disp(num);disp(分母系数是 );disp(den);impres = filter(num,den,1 zeros(1,4*N1);p,d = strucver(impres,2*N1);disp(实际分子系数是 ); disp(p);disp(实际分母系数是 ); disp(d);运行结果N1 = 2Wn1 = 0.3608 0.5421分子系数是 0.0571 0 -0.1143 0 0.0571分母系数是 1
8、.0000 -0.5099 1.2862 -0.3350 0.4479实际分子系数是 0.0571 0.0000 -0.1143 0 0.0571实际分母系数是 1.0000 -0.5099 1.2862 -0.3350 0.4479Q8.3% Program P8clear all;clc;Wp = 0.4 0.5; Ws = 0.1 0.8; Rp = 1; Rs = 30;N, Wn = buttord(Wp, Ws, Rp, Rs)num,den = butter(N,Wn);z,p,k = tf2zp(num,den)sos = zp2sos(z,p,k)num1=sos(1,1)
9、sos(1,2) sos(1,3);den1=sos(1,4) sos(1,5) sos(1,6);disp(第一级的分子系数是 );disp(num1);disp(第一级的分母系数是 );disp(den1);num2=sos(2,1) sos(2,2) sos(2,3);den2=sos(2,4) sos(2,5) sos(2,6);disp(第二级的分子系数是 );disp(num2);disp(第二级的分母系数是 );disp(den2);impres1 = direct2(num1,den1,1 zeros(1,2*N);p1,d1 = strucver(impres1,N);dis
10、p(第一级的实际分子系数是 ); disp(p1);disp(第一级的实际分母系数是 ); disp(d1);impres2 = direct2(num2,den2,1 zeros(1,2*N);p2,d2 = strucver(impres2,N);disp(第二级的实际分子系数是 ); disp(p2);disp(第二级的实际分母系数是 ); disp(d2);运行结果N = 2 Wn = 0.3608 0.5421sos = 0.0571 0.1143 0.0571 1.0000 0.0877 0.6601 1.0000 -2.0000 1.0000 1.0000 -0.5976 0.6
11、785第一级的分子系数是 0.0571 0.1143 0.0571第一级的分母系数是 1.0000 0.0877 0.6601第二级的分子系数是 1.0000 -2.0000 1.0000第二级的分母系数是 1.0000 -0.5976 0.6785第一级的实际分子系数是 0.0571 0.1143 0.0571第一级的实际分母系数是 1.0000 0.0877 0.6601第二级的实际分子系数是 1.0000 -2.0000 1.0000第二级的实际分母系数是 1.0000 -0.5976 0.6785Q8.5clear all;clc;Wp = 0.4 0.5; Ws = 0.1 0.8;
12、 Rp = 1; Rs = 30;N, Wn = buttord(Wp, Ws, Rp, Rs)num,den = butter(N,Wn);r1,p1,k1 = residuez(num,den);disp(并联I型)disp(留数是);disp(r1);disp(极点是);disp(p1);disp(常数);disp(k1);numI1 = 2*real(r1(1) -2*(real(r1(1)*real(p1(1)+imag(r1(1)*imag(p1(1);denI1 = 1 -2*real(p1(1) real(p1(1)2+imag(p1(1)2;numI2 = 2*real(r1
13、(3) -2*(real(r1(3)*real(p1(3)+imag(r1(3)*imag(p1(3);denI2 = 1 -2*real(p1(3) real(p1(3)2+imag(p1(3)2;disp(并联I型的常数是);disp(k1);disp(并联I型的第一级的分子系数是 );disp(numI1);disp(并联I型的第一级的分母系数是 );disp(denI1);disp(并联I型的第二级的分子系数是 );disp(numI2);disp(并联I型的第二级的分母系数是 );disp(denI2);impresI1 = direct2(numI1,denI1,1 zeros(1
14、,2*N);pI1,dI1 = strucver(impresI1,N);disp(并联I型的第一级的实际分子系数是 ); disp(pI1);disp(并联I型的第一级的实际分母系数是 ); disp(dI1);impresI2 = direct2(numI2,denI2,1 zeros(1,2*N);pI2,dI2 = strucver(impresI2,N);disp(并联I型的第二级的实际分子系数是 ); disp(pI2);disp(并联I型的第二级的实际分母系数是 ); disp(dI2);运行结果如下N = 2 Wn = 0.3608 0.5421并联I型留数是 -0.0235
15、+ 0.1978i -0.0235 - 0.1978i -0.0117 - 0.2132i -0.0117 + 0.2132i极点是 0.2988 + 0.7676i 0.2988 - 0.7676i -0.0438 + 0.8113i -0.0438 - 0.8113i常数 0.1276并联I型的常数是 0.1276并联I型的第一级的分子系数是 -0.0470 -0.2897并联I型的第一级的分母系数是 1.0000 -0.5976 0.6785并联I型的第二级的分子系数是 -0.0234 0.3448并联I型的第二级的分母系数是 1.0000 0.0877 0.6601并联I型的第一级的实
16、际分子系数是 -0.0470 -0.2897 -0.0000并联I型的第一级的实际分母系数是 1.0000 -0.5976 0.6785并联I型的第二级的实际分子系数是 -0.0234 0.3448 0.0000并联I型的第二级的实际分母系数是 1.0000 0.0877 0.6601Q8.6% Program P8clear all;clc;Wp = 0.4 0.5; Ws = 0.1 0.8; Rp = 1; Rs = 30;N, Wn = buttord(Wp, Ws, Rp, Rs)num,den = butter(N,Wn);r2,p2,k2 = residue(num,den);d
17、isp(并联II型)disp(留数是);disp(r2);disp(极点是);disp(p2);disp(常数是);disp(k2);numII1 = 0 2*real(r2(1) -2*(real(r2(1)*real(p2(1)+imag(r2(1)*imag(p2(1);denII1 = 1 -2*real(p2(1) real(p2(1)2+imag(p2(1)2;numII2 = 0 2*real(r2(3) -2*(real(r2(3)*real(p2(3)+imag(r2(3)*imag(p2(3);denII2 = 1 -2*real(p2(3) real(p2(3)2+ima
18、g(p2(3)2;disp(并联II型的常数是);disp(k2);disp(并联II型的第一级的分子系数是 );disp(numII1);disp(并联II型的第一级的分母系数是 );disp(denII1);disp(并联II型的第二级的分子系数是 );disp(numII2);disp(并联II型的第二级的分母系数是 );disp(denII2);impresII1 = direct2(numII1,denII1,1 zeros(1,2*N);pII1,dII1 = strucver(impresII1,N);disp(并联II型的第一级的实际分子系数是 ); disp(pII1);di
19、sp(并联II型的第一级的实际分母系数是 ); disp(dII1);impresII2 = direct2(numII2,denII2,1 zeros(1,2*N);pII2,dII2 = strucver(impresII2,N);disp(并联II型的第二级的实际分子系数是 ); disp(pII2);disp(并联II型的第二级的实际分母系数是 ); disp(dII2);运行结果N = 2 Wn = 0.3608 0.5421并联II型留数是 -0.1589 + 0.0411i -0.1589 - 0.0411i 0.1734 - 0.0002i 0.1734 + 0.0002i极点
20、是 0.2988 + 0.7676i 0.2988 - 0.7676i -0.0438 + 0.8113i -0.0438 - 0.8113i常数是 0.0571并联II型的常数是 0.0571并联II型的第一级的分子系数是 0 -0.3178 0.0319并联II型的第一级的分母系数是 1.0000 -0.5976 0.6785并联II型的第二级的分子系数是 0 0.3469 0.0155并联II型的第二级的分母系数是 1.0000 0.0877 0.6601并联II型的第一级的实际分子系数是 0 -0.3178 0.0319并联II型的第一级的实际分母系数是 1.0000 -0.5976
21、0.6785并联II型的第二级的实际分子系数是 0 0.3469 0.0155并联II型的第二级的实际分母系数是 1.0000 0.0877 0.6601Q8.7% Program P8_2clear all;clc;x1 = 1 zeros(1,6); %产生单位冲激序列q3old = 0; q2old = 0; q1old = 0; %初始条件% 输入滤波器系数D1 = 0.; D2 = 0.27083; D3 =-0.2;alpha1 = 0.02; alpha2 = 0.352; alpha3 = 0.; alpha4 = - 0.19016;%计算前七个冲激响应样本for n = 1
22、:7w3 = x1(n) - D3*q3old;w2 = w3 - D2*q2old;q1new = w2 - D1*q1old;q2new = D1*q1new + q1old;q3new = D2*w2 + q2old;y1 = D3*w3 + q3old;yo(n) = alpha1*y1 + alpha2*q3new + alpha3*q2new + alpha4*q1new;q3old = q3new; q2old = q2new; q1old = q1new;endnum,den = strucver(yo,3);disp(分子系数);disp(num);disp(分母系数);di
23、sp(den);运行结果分子系数 -0.0000 0.4400 0.3600 0.0200分母系数1.0000 0.4000 0.1800 -0.2000Q8.9% Program P8_3%一个无限冲激响应滤波器滤波的说明clear all;clc;clf;%产生输入序列k = 0:50;w2 = 0.7*pi;w1 = 0.2*pi;x1 = 1.5*cos(w1*k); x2 = 2*cos(w2*k);x = x1+x2;% 确定滤波器传输函数N, Wn = ellipord(0.25, 0.55, 0.5, 50)num, den = ellip(N,0.5, 50,Wn);% 产生
24、输出序列y = filter(num,den,x);% 画出输入及输出序列subplot(2,2,1);stem(k,x1); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输入序列中的低频成分);subplot(2,2,2);stem(k,x2); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输入序列中的高频成分);subplot(2,2,3);stem(k,x); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);
25、title(输入序列);subplot(2,2,4);stem(k,y); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输出序列);运行结果N = 4 Wn = 0.2500% Program P8_3%一个无限冲激响应滤波器滤波的说明clear all;clc;clf;%产生输入序列clear all;clc;clf;k = 0:50;w2 = 0.7*pi;w1 = 0.2*pi;x1 = 1.5*cos(w1*k); x2 = 2*cos(w2*k);x = x1+x2;% 确定滤波器传输函数N, Wn = ellipor
26、d(0.55, 0.25, 0.5, 50)num, den = ellip(N,0.5, 50,Wn,high);% 产生输出序列y = filter(num,den,x);% 画出输入及输出序列subplot(2,2,1);stem(k,x1); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输入序列中的低频成分);subplot(2,2,2);stem(k,x2); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输入序列中的高频成分);subplot(2,2,3
27、);stem(k,x); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输入序列);subplot(2,2,4);stem(k,y); grid; axis(0 50 -4 4);xlabel(时间序列 n); ylabel(振幅);title(输出序列);运行结果N = 4 Wn = 0.5500Q8.14% Program P8_4clear all;clc;num = remez(9, 0 0.3 0.5 1,1 1 0 0);disp(滤波器系数是);disp(num);impres = filter(num,1 ,1 z
28、eros(1,9);disp(滤波器实际系数是 ); disp(impres);程序运行结果滤波器系数是 -0.0209 -0.0926 0.0064 0.1932 0.3721 0.3721 0.1932 0.0064 -0.0926 -0.0209滤波器实际系数是 -0.0209 -0.0926 0.0064 0.1932 0.3721 0.3721 0.1932 0.0064 -0.0926 -0.0209Q8.15% Program P8clear all;clc;num = remez(9, 0 0.3 0.5 1,1 1 0 0);disp(滤波器系数是);disp(num);im
29、pres = direct2(num,1 ,1 zeros(1,9);disp(滤波器实际系数是 ); disp(impres);运行结果滤波器系数是 -0.0209 -0.0926 0.0064 0.1932 0.3721 0.3721 0.1932 0.0064 -0.0926 -0.0209滤波器实际系数是 -0.0209 -0.0926 0.0064 0.1932 0.3721 0.3721 0.1932 0.0064 -0.0926 -0.0209Q8.16% Program P8clear all;clc;num = remez(9, 0 0.3 0.5 1,1 1 0 0);de
30、n = 1 zeros(1,length(num)-1);z,p,k = tf2zp(num,den);sos = zp2sos(z,p,k)% size(sos,1)为sos的行数，代表级数% size(sos,2)为sos的列数N=length(num)/size(sos,1)for i=1:size(sos,1) num1=sos(i,1) sos(i,2) sos(i,3); disp(第,num2str(i), 级滤波器的系数是 );disp(num1); impres = direct2(num,1,1 zeros(1,N); disp(第,num2str(i),级滤波器的实际系数
31、是 ); disp(impres);end 运行结果sos = -0.0209 -0.0870 0 1.0000 0 0 1.0000 -0.9623 -1.9623 1.0000 0 0 1.0000 -0.2695 -0.1223 1.0000 0 0 1.0000 1.2970 1.0000 1.0000 0 0 1.0000 0.2002 1.0000 1.0000 0 0第1级滤波器的系数是 -0.0209 -0.0870 0第1级滤波器的实际系数是 -0.0209 -0.0870 0第2级滤波器的系数是 1.0000 -0.9623 -1.9623第2级滤波器的实际系数是 1.0000 -0.9623 -1.9623第3级滤波器的系数是 1.0000 -0.2695 -0.1223第3级滤波器的实际系数是 1.0000 -0.2695 -0.1223第4级滤波器的系数是 1.0000 1.2970 1.0000第4级滤波器的实际系数是 1.0000 1.2970 1.0000第5级滤波器的系数是 1.0000 0.2002 1.0000第5级滤波器的实际系数是 1.0000 0.2002 1.0000七、实验思考题及解答八、实验结果分析与总结专心-专注-专业