- 1、本文档共31页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
《数字信号处理》课程研究性学习报告
姓名
学号
同组成员
指导教师
时间 2014/5/24
IIR和FIR滤波器设计专题研讨
【目的】
(1) 掌握根据滤波器指标设计IIR和FIR数字滤波器的原理和方法。
(2) 熟悉通过IIR和FIR数字滤波器进行实际系统设计的方法。
(3) 培养学生自主学习能力,以及发现问题、分析问题和解决问题的能力。
【研讨题目】
1. 设计一个数字滤波器,在频率低于的范围内,低通幅度特性为常数,并且不低于0.75dB。在频率和之间,阻带衰减至少为20dB。
(1)试求满足这些条件的最低阶Butterworth滤波器。
(2)试求满足这些条件的最低阶Chebyshev I滤波器。
(3)自主选择一段带限信号,通过所设计的(1)、(2)两种滤波器,比较各自的输入和输出信号。讨论两种滤波器在结构和性能上的差异。
【题目分析】
本题讨论模拟滤波器和数字滤波器的设计。结合课本分析各类滤波器的特性。
【IIR模拟滤波器设计的基本方法】
在设计IIR模拟滤波器的时候,
先将待设计的模拟滤波器技术指标转换为原型模拟低通滤波器的指标;
然后设计原型模拟低通滤波器(有Butterworth、Chebyshev和椭圆等设计方法);
再通过频率变换将模拟低通滤波器转换为所需的滤波器。
【仿真程序】
Butterworth滤波器:
Wp=0.2613*pi;
Ws=0.4018*pi;
Ap=0.75;
As=20;
Fs=1;
wp=Wp*Fs;
ws=Ws*Fs;
[N,wc]=buttord(wp,ws,Ap,As,s);
[num_a,den_a]=butter(N,wc,s);
[num_d,den_d]=impinvar(num_a,den_a,Fs);
w=linspace(0,pi,1024);
h=freqz(num_d,den_d,w);
norm=max(abs(h));
num_d=num_d/norm;
plot(w/pi,20*log10(abs(h/norm)));
xlabel(归一化频率);
ylabel(Gain,dB);
title(Butterwort滤波器);
Chebyshev I滤波器:
Wp=0.2613*pi;
Ws=0.4018*pi;
Ap=0.75; As=20;
Fs=1;
wp=Wp*Fs;
ws=Ws*Fs;
[N,wc]=cheb1ord(wp,ws,Ap,As,s);
[num_a,den_a]=cheby1(N,Ap,wc,s);
[num_d,den_d]=impinvar(num_a,den_a,Fs);
w=linspace(0,pi,1024);
h=freqz(num_d,den_d,w);
norm=max(abs(h));
num_d=num_d/norm;
plot(w/pi,20*log10(abs(h/norm)));
xlabel(归一化频率);
ylabel(Gain,dB);
title(Chebyshev Ⅰ滤波器);
滤波:
1)Butterworth滤波器
Wp=0.2613*pi;
Ws=0.4018*pi;
Ap=0.75;
As=20;
Fs=1;
wp=Wp*Fs; ws=Ws*Fs;
[N,wc]=buttord(wp,ws,Ap,As,s);
[num_a,den_a]=butter(N,wc,s);
[num_d,den_d]=impinvar(num_a,den_a,Fs);
w=linspace(0,pi,1024);
h=freqz(num_d,den_d,w);
x=w;
s=h.*w;
plot(w,abs(s));
2)Chebyshev I滤波器
Wp=0.2613*pi;
Ws=0.4018*pi;
Ap=0.75;
As=20;
Fs=1;
wp=Wp*Fs;
ws=Ws*Fs;
[N,wc]=cheb1ord(wp,ws,Ap,As,s);
[num_a,den_a]=cheby1(N,Ap,wc,s);
[num_d,den_d]=impinvar(num_a,den_a,Fs);
w=linspace(0,pi,1024);
h=freqz(num_d,den_d,w);
x=w;
s=h.*w;
plot(w,abs(s));
【结果分析】
Butterworth
文档评论(0)