数字滤波器实现分析Word格式.docx
- 文档编号:811035
- 上传时间:2023-04-29
- 格式:DOCX
- 页数:15
- 大小:81.82KB
数字滤波器实现分析Word格式.docx
《数字滤波器实现分析Word格式.docx》由会员分享,可在线阅读,更多相关《数字滤波器实现分析Word格式.docx(15页珍藏版)》请在冰点文库上搜索。
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.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型的第一级的实际分子系数是
-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.6601
Q8.6
[r2,p2,k2]=residue(num,den);
并联II型'
disp(r2);
disp(p2);
常数是'
disp(k2);
numII1=[02*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=[02*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+imag(p2(3))^2];
并联II型的常数是'
并联II型的第一级的分子系数是'
disp(numII1);
并联II型的第一级的分母系数是'
disp(denII1);
并联II型的第二级的分子系数是'
disp(numII2);
并联II型的第二级的分母系数是'
disp(denII2);
impresII1=direct2(numII1,denII1,[1zeros(1,2*N)]);
[pII1,dII1]=strucver(impresII1,N);
并联II型的第一级的实际分子系数是'
disp(pII1'
并联II型的第一级的实际分母系数是'
disp(dII1'
impresII2=direct2(numII2,denII2,[1zeros(1,2*N)]);
[pII2,dII2]=strucver(impresII2,N);
并联II型的第二级的实际分子系数是'
disp(pII2'
并联II型的第二级的实际分母系数是'
disp(dII2'
运行结果
并联II型留数是
-0.1589+0.0411i
-0.1589-0.0411i
0.1734-0.0002i
0.1734+0.0002i
常数是
并联II型的常数是
并联II型的第一级的分子系数是
0 -0.3178 0.0319
并联II型的第一级的分母系数是
并联II型的第二级的分子系数是
0 0.3469 0.0155
并联II型的第二级的分母系数是
并联II型的第一级的实际分子系数是
并联II型的第一级的实际分母系数是
并联II型的第二级的实际分子系数是
并联II型的第二级的实际分母系数是
Q8.7
%ProgramP8_2clearall;
x1=[1zeros(1,6)];
%产生单位冲激序列q3old=0;
q2old=0;
q1old=0;
%初始条件
%输入滤波器系数
D1=0.357377;
D2=0.27083;
D3=-0.2;
alpha1=0.02;
alpha2=0.352;
alpha3=0.276533;
alpha4=-0.19016;
%计算前七个冲激响应样本forn=1:
7
w3=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;
end
[num,den]=strucver(yo,3);
分子系数'
disp(num'
分母系数'
disp(den'
运行结果分子系数
-0.0000 0.4400 0.3600 0.0200
分母系数
1.0000 0.4000 0.1800 -0.2000
Q8.9
%ProgramP8_3
%一个无限冲激响应滤波器滤波的说明clearall;
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);
%产生输出序列
y=filter(num,den,x);
%画出输入及输出序列
subplot(2,2,1);
stem(k,x1);
grid;
axis([050-44]);
xlabel('
时间序列n'
ylabel('
振幅'
title('
输入序列中的低频成分'
subplot(2,2,2);
stem(k,x2);
输入序列中的高频成分'
subplot(2,2,3);
stem(k,x);
输入序列'
subplot(2,2,4);
stem(k,y);
输出序列'
N= 4 Wn= 0.2500
□中中中中中中中中中
4
□中中中中中中中中中
2 2
□中
0 0
-2
-4
0 10 20
30 40 50
0 10
20 30
40 50
□中中中n
□中中中
□中中中n
□中中中
-2
0 10 20 30 40 50
□中中中n 中中中中n
%产生输入序列clearall;
k=0:
[N,Wn]=ellipord(0.55,0.25,0.5,50)
[num,den]=ellip(N,0.5,50,Wn,'
high'
N= 4 Wn= 0.5500
Q8.14
%ProgramP8_4clearall;
num=remez(9,[00.30.51],[1100]);
滤波器系数是'
impres=filter(num,1,[1zeros
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字滤波器 实现 分析
![提示](https://static.bingdoc.com/images/bang_tan.gif)