电子科技大学生命科学与技术学院标准实验报告
（实验）课程名称数字信号处理
2017-2018-第2学期
电子科技大学教务处制表
电子科技大学
实验报告
一、实验室名称：
二、实验名称：
滤波器传输函数的零点和极点对滤波特性的影响
三、实验学时：2学时
四、实验原理：
五、实验目的：
1. 直观地了解滤波器传输函数的零点和极点（的个数和位置）对滤波特性的影响。

2. 利用设计的滤波器进行滤波。

六、实验内容：
H(z):
1、给定某个滤波器的传输函数
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编程方法。

十二、对本实验过程及方法、手段的改进建议：
报告评分：
指导教师签字：。