- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
設计数字低通滤波器(用matlab实现)
DSP设计滤波器报告(一)实验目的:IIR数字低通滤波器的设计方法。
(二)实验原理:
1中的有用成分和希望取出的成分各自占有不同的频带。这样,当通过一个线性系统(即滤波器)后可讲欲去除的成分有效的去除。
现代滤波器理论研究的主要内容是从含有噪声的数据记录(又称时间序列)中估计出信号的某些特征或信号本身。
经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器(AF)和数字滤波器(DF)。对数字滤波器,又有IIR滤波器和FIR滤波器。
IIR DF 的转移函数是:
FIR DF的转移函数是:
FIR滤波器可以对给定的频率特性直接进行设计,而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。
2、滤波器的技术要求
低通滤波器::通带截止频率(又称通带上限频率) :阻带下限截止频率
:通带允许的最大衰减 :阻带允许的最小衰减 (,的单位dB)
:通带上限角频率 :阻带下限角频率
(,)即
3、IIR 数字滤波器的设计步骤:
按一定规则将给出的数字滤波器的技术指标转换魏模拟低通滤波器的技术指标。
根据转换后的技术指标设计模拟低通滤波器;
再按一定的规则将转换成。
若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标,然后按上述步骤2)设计出低通,再将转换为所需的。
4.几种不同类型的滤波器的介绍:
因为我们设计的滤波器的冲击响应一般都为实数,所以有
这样,如果我们能由,,,求出,那么就容易得到所需要的。
不同类型的的表达式,代表了几种不同类型的滤波器。
(1)巴特沃思(Butterworth)滤波器:
C为待定常数,N为待定的滤波器阶次。
(2)切比雪夫I型(Chebyshev – I )滤波器:
5.巴特沃思模拟低通滤波器的设计
由于每一个滤波器的频率范围将直接取决于设计者的应用目的,因此必然是千差万别。为了使设计规范化,我们需要将滤波器的频率参数作归一化处理。设所给的实际频率为(或),归一化后的频率为,对低通模拟滤波器,令
=
显然,,。又令归一化复数变量为,,显然
所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。
(1)将实际频率规一化
(2)求C和N
这样C和N可求。
若令=3dB,则C=1,这样巴特沃思滤波器的设计就只剩一个参数N,这时
(3)确定
因为,根据上面公式有
由 解得
,=1,2,···,2N
这样可得
求得后,用代替变量,即得实际需要得。
6.用双线性Z变换法设计IIR数字低通滤波器
s平面到z平面的映射关系
称为双线性Z变换,由此关系求出
及
即
因为设计滤波器时系数会被约掉,所以又有
(三)实验内容
Hz,Hz,dB, dB,抽样频率Hz。
提示:首先应该得到角频率,然后再
将数字滤波器的技术要求转换为模拟滤波器的技术要求。
设计低通滤波器
由=依次求出,,再求出N,可得
然后由转换成
由求
源程序:fp=100;
fs=300;
ap=3; %通带最大衰减
as=20; %阻带最小衰减
Fs1=1000; %抽样频率
wp=2*pi*fp/Fs1;
ws=2*pi*fs/Fs1; %数字角频率
Fs=Fs1/Fs1;
T=1/Fs;
Op=2/T*tan(wp/2);
Os=2/T*tan(ws/2); %角频率畸变
[N Wn]=buttord(Op,Os,ap,as,s);
[z p k]=buttap(N); %G(p)的零点,极点,增益
[b a]=zp2tf(z,p,k); %G(p)的分子、分数系数
[B A]=lp2lp(b,a,Op); %H(s)的分子、分数系数
[Bz Az]=bilinear(B,A,Fs*T/2); %H(z)的分子、分数系数
[H w]=freqz(Bz,Az,256,Fs*Fs1); %H
plot(w,abs(H),r);
title(低通滤波器);
grid on;
实验结果:
MATLAB参考程序和仿真内容
%******************************************************
文档评论(0)