数字滤波器实现分析Word下载.docx
《数字滤波器实现分析Word下载.docx》由会员分享，可在线阅读，更多相关《数字滤波器实现分析Word下载.docx（15页珍藏版）》请在冰豆网上搜索。
验证仿真。
Q8.5 生成Q8.1中传输函数的一个并联Ⅰ型实现，并编写一个程序来仿真它。
每一个单独的部分用直接Ⅱ型实现。
Q8.6 生成Q8.1中传输函数的一个并联II型实现，并编写一个程序来仿真它。
Q8.7运行程序P8.2并产生被仿真的级联格型结构的传输函数。
该传输函数与式
（8.35）中的相同吗？
Q8.9程序设计了什么类型的滤波器？
形成输入的正弦序列的频率是多少？
Q8.10运行P8.3并产生两个图形。
哪种输入成分会在滤波器输出出现？
为什么输出序列的开始部分不是一种理想的正弦曲线？
修改P8.3，以便纸过滤序列X2[n]。
产生的输出序列和预料的是一样的吗？
证明你的答案。
Q8.14 程序P8.4设计了什么类型的滤波器？
为了验证仿真，需要多少个冲激响应样本？
Q8.15修改程序P8.4，用直接II型仿真该滤波器，运行修改后的程序。
Q8.16使用程序P8.4生成在习题Q8.4中产生的传输函数的一个级联实现，写出对其仿真的一个程序。
其中每一个部分以直接II型实现。
三、实验器材及软件
1.微型计算机1台
2.MATLAB7.0软件
2
四、实验原理
1、可调数字滤波器的设计基础是复变量的z谱变换，通过z替换为F（z），将具有给定截止频率的数字滤波器的实现调整为具有不同截止频率的另一个实现。
2、利用全通结构代替每一个延时块，可能会导致实现Gnew（z）时具有不能被实现的无延时环结构。
3、设计可调有限冲激响应滤波器的一种简单方法是基于加窗的傅里叶级数方法。
4、通常情况下，复数乘法和复数加法总数量是N（log2N），通使用Wn的对称性质，可进一步简化计算得复杂性。
五、实验步骤
1、打开计算机中的MATLAB
2、在MATLAB中新建一个M-file，在其中输入实验要求的相关代码
3、输入代码后调试代码，代码无误后可以得到相应的实验结果
4、观察实验结果回答相关问题
5、记录实验结果，将图形等截图粘贴，完成实验报告
六、实验记录（数据、图表、波形、程序等）
Q8.1
%ProgramP8_1clearall;
clc;
Wp=[0.40.5];
Ws=[0.10.8];
Rp=1;
Rs=30;
[N1,Wn1]=buttord（Wp,Ws,Rp,Rs）[num,den]=butter（N1,Wn1）;
disp（'
分子系数是'
）;
disp（num）;
分母系数是'
disp（den）;
impres=direct2（num,den,[1zeros（1,4*N1）]）;
[p,d]=strucver（impres,2*N1）;
实际分子系数是'
disp（p'
实际分母系数是'
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
实际分母系数是
Q8.2
impres=filter（num,den,[1zeros（1,4*N1）]）;
运行结果
0.0571 0.0000 -0.1143 0 0.0571
Q8.3
%ProgramP8clearall;
[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）sos（1,2）sos（1,3）];
den1=[sos（1,4）sos（1,5）sos（1,6）];
第一级的分子系数是'
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（num2）;
13
第二级的分母系数是'
disp（den2）;
impres1=direct2（num1,den1,[1zeros（1,2*N）]）;
[p1,d1]=strucver（impres1,N）;
第一级的实际分子系数是'
disp（p1'
第一级的实际分母系数是'
disp（d1'
impres2=direct2（num2,den2,[1zeros（1,2*N）]）;
[p2,d2]=strucver（impres2,N）;
第二级的实际分子系数是'
disp（p2'
第二级的实际分母系数是'
disp（d2'
N=
sos=
Wn=
0.3608
0.5421
0.0571
0.1143
1.0000
0.0877
0.6601
-2.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.6785
第一级的实际分子系数是
第一级的实际分母系数是
第二级的实际分子系数是
第二级的实际分母系数是
Q8.5
clearall;
[r1,p1,k1]=residuez（num,den）;
1
并联I型'
）
留数是'
disp（r1）;
极点是'
disp（p1）;
常数'
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（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];
并联I型的常数是'
并联I型的第一级的分子系数是'
disp（numI1）;
并联I型的第一级的分母系数是'
disp（denI1）;
并联I型的第二级的分子系数是'
disp（numI2）;
并联I型的第二级的分母系数是'
disp（denI2）;
impresI1=direct2（numI1,denI1,[1zeros（1,2*N）]）;
[pI1,dI1]=strucver（impresI1,N）;
并联I型的第一级的实际分子系数是'
disp（pI1'
并联I型的第一级的实际分母系数是'
disp（dI1'
impresI2=direct2（numI2,denI2,[1zeros（1,2*N）]）;
[pI2,dI2]=strucver（impresI2,N）;
并联I型的第二级的实际分子系数是'
disp（pI2'
并联I型的第二级的实际分母系数是'
disp（dI2'
运行结果如下
N=2 Wn=0.36080.5421
并联I型留数是
-0.0235+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
升级会员 