- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验五 IIR 滤波器实验
一 实验目的
(1) 了解IIR 滤波器的原理及使用方法;
(2) 了解使用Matlab 语言设计IIR 滤波器的方法;
(3) 了解DSP 对IIR 滤波器的设计及编程方法;
(4) 熟悉对IIR 滤波器的调试方法;
二 实验内容
本试验要求设计滤波器采样频率为600hz ,截止频率200hz 的高通滤波器。设计IIR 滤
波器实现上面要求。
输入信号频率为100HZ 和300HZ 的合成信号,目的是通过我们设计的滤波器将100HZ
的信号滤掉,余下300HZ 的信号成分,达到滤波的效果。
三 实验原理
数字滤波器的输入x[k]和输出y[k]之间的关系可以用如下常系数线性差分方程及其z变换描述:
系统的转移函数为:
设N=M ,则传输函数变为:
转换成极零点表示为:
式中,表示零点,表示极点,它具有N个零点和N个极点,如果任何一个极点在Z平面单位圆外,则系统不稳定。如果系数bj全都为0,滤波器成为非递归的FIR滤波器, 这时系统没有极点,因此FIR 滤波器总是稳定的。对于IIR 滤波器,有系数量化敏感的缺点。由于系统 对序列施加的算法,是由加法、延时和常系数乘三种基本运算的组合,所以可以用不同结构的数字滤波器来实现而不影响系统总的传输函数。图4.1 是四阶直接型IIR 滤波器的结构。
图 4.1 四阶直接型IIR 滤波器的结构
四 实验步骤
1.滤波器的Matlab 语言设计
在Matlab 中使用滤波器设计工具箱(FDA )设计滤波器。
图 4.2 打开的FDA 主界面图
上图 4.2 是打开的滤波器设计的主界面图,在图4.2 上可以设计各种满足用户要求的
滤波器,包括滤波器的类型(IIR 或 FIR )、滤波器的阶数、滤波器的种类、滤波器的截至频率、带宽、纹波系数、采样频率等各种和设计滤波器有关的所有参数。
图中,我们设计IIR 型的3 阶切比雪夫I 型的采样频率为600HZ ,截至频率为200HZ 的低通滤波器。图中中间显示的数据为设计好的滤波器参数。
下图 4 .3 是对应频率特性图。
图 4.3 滤波器的频率特性图
下图 4.4 是对应的冲击响应特性图。
图 4 .4 滤波器的冲击响应特性图
图 4 .5 是对应的滤波器阶跃特性图。
图 4 .6 是对应的滤波器极零点示意图。
图 4 .7 是对应的滤波器结构示意图。
图 4 .5 滤波器阶跃特性图。
图 4 .6 滤波器的极零点示意图
图 4 .7 滤波器的结构示意图
2.从 Matlab 语言转换成通用语言
使用Matlab 语言设计出滤波器的参数之后,紧接着需要编写通用的语言,以便向DSP
所需要的汇编语言转换,也直接使用C 语言,然后调用CCS 自带的C 编译器将C 语言转换
成汇编语言,但一般情况下,滤波器对实时要求比较高,而整个滤波器的程序编写也不是很
大,所以建议采用汇编语言编写。
滤波器设计的第二步需要编写通用的语言,在这里,我们仍然使用Matlab 语言编写易
懂的语句。下面程序是编写后的程序。浅显易懂,移值性强。
********************* 滤波器程序 *****************************************
clear all; 清寄存器值
clf; 清屏
N=256; 数据点数
Fs=600; 采样频率
dt=1/fs
for K=1:N
f1=300; 信号频率
f2=100; 信号频率
y(k)=sin(2*pi*f1*k*dt)+sin(2*pi*f2*k*dt)+sin(2*(f1+20)*k*dt); 产生信号
end
lp=200; 截止频率
wn1=2lp/fs; 函数的参数
[zl,pl,kl]=CHEBY1(3,0.5,wnl,’high’); 滤波器的极零点表示
[B,A]= CHEBY1(3,0.5,wnl,’high’); 滤波器的传递函数表示
yy1(1)=0;
yy1(2)=0;
b(1)=0.0580;
文档评论(0)