数字信号处理第四次实验报告.docx
- 文档编号:27039799
- 上传时间:2023-06-26
- 格式:DOCX
- 页数:9
- 大小:430.95KB
数字信号处理第四次实验报告.docx
《数字信号处理第四次实验报告.docx》由会员分享,可在线阅读,更多相关《数字信号处理第四次实验报告.docx(9页珍藏版)》请在冰豆网上搜索。
数字信号处理第四次实验报告
电子科技大学生命科学与技术学院
标准实验报告
(实验)课程名称数字信号处理
2017-2018-第2学期
电子科技大学教务处制表
电子科技大学
实验报告
学生姓名
学号
指导教师
实验时间
年月日
一、实验室名称:
二、实验名称:
滤波器传输函数的零点和极点对滤波特性的影响
三、实验学时:
2学时
四、实验原理:
五、实验目的:
1.直观地了解滤波器传输函数的零点和极点(的个数和位置)对滤波特性的影响。
2.利用设计的滤波器进行滤波。
六、实验内容:
1、给定某个滤波器的传输函数H(z):
增加其零点和极点的个数,分析新的滤波器的滤波特性(幅度谱)。
2、利用设计滤波器进行滤波,分析滤波器的相位谱特性对滤波结果的影响
九、实验数据及结果分析:
1.b=conv([0,0.0798],[1,0.6566]);
a=conv([1,-0.5335],[1,-(0.625+0.3773i)]);
a=conv(a,[1,-(0.625-0.3773i)]);
zplane(b,a);
[h,w]=freqz(b,a);
figure
(2);plot(w/pi,abs(h))
r=0.5;
w=pi/4;
Rre=r*cos(w);
Rim=r*sin(w);
a1=Rre+i*Rim;
a2=Rre-i*Rim;
b1=0.8*(Rre+i*Rim);
b2=0.8*(Rre-i*Rim);
b=conv(b,[1,b1]);
b=conv(b,[1,b2]);
a=conv(a,[1,a1]);
a=conv(a,[1,a2]);
figure(3)
zplane(b,a);
[h,w]=freqz(b,a);
b=b/max(abs(h));
figure(4);plot(w/pi,abs(h)/max(abs(h)))
clear
x=1:
0.01:
10;
y=sin(5*x).*exp(0.5*x)+sin(5*x);
ffty=fft(y);
amp=abs(ffty);
pha=angle(ffty);
fliterPha=2*rand(1,length(pha));
phaNew=pha+fliterPha;
cc1=amp.*exp(i*phaNew);
yy1=ifft(cc1,'symmetric');
fliterPha2=linspace(0,-200,length(y));
phaNew2=pha+fliterPha2;
cc2=amp.*exp(i*phaNew2);
yy2=ifft(cc2,'symmetric');
plot(1:
length(x),y,'g',1:
length(x),yy1,'b',1:
length(x),yy2,'r')
原系统图像
新增零极点后的图像
幅度谱对比图
n=0:
99;
xs=cos(n*0.1*pi);
xn=sin(n*0.9*pi);
x=xs+xn;
figure(5)
plot(n,xs,'r',n,x,'b');legend('有用信号','加噪声信号');title('信号幅度')
[h4,w4]=freqz(x);
figure(6);
plot(w4/pi,abs(h4));title('滤波前信号频谱')
q=filtfilt(b,a,x);
figure(7)
[h3,w3]=freqz(q);
plot(w/pi,abs(h3));title('滤波后信号频谱');
figure(8);
plot(n,q,'r',n,xs,'b');legend('有用信号','滤波后');title('滤波前后')
十、实验结论:
(详细填写)
1.添加零点可以将幅度谱曲线向下弯曲
2.合理添加零极点可以更好地滤波
十一、总结及心得体会:
(详细填写)
增加了对z变换意义的理解,熟悉了关于z变换的matlab编程方法。
十二、对本实验过程及方法、手段的改进建议:
报告评分:
指导教师签字:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 第四 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)