- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
.
..
XXXX大学XXXX学院
课程名称:
数字信号处理
班级:
姓名:
学号:
实验地点:
日期:
实验名称 IIR数字滤波器的设计
实验目的:
加深理解IIR数字滤波器的时域特性和频域特性,掌握IIR数字滤波器的设计原理与设计方法,以及IIR数字滤波器的应用。
实验内容:
IIR数字滤波器一般为线性移不变的因果离散系统,N阶IIR数字滤波器的系统函数可以表达为的有理多项式,即
式中:系数至少有一个非零。对于因果IIR数据滤波器,应满足。
IIR数字滤波器的设计主要通过成熟的模拟滤波器设计方法来实现。
首先在频域将数字滤波器设计指标转换为模拟滤波器设计指标,然后将任意的模拟滤波器为原型模拟低通滤波器指标,根据模拟滤波器的设计指标来设计出模拟低通滤波器,然后又经过相应的复频域转换得到H(s),最后又H(s)经过脉冲响应不变法或双线性变换法得到所需要的IIIR数字滤波器H(z)。
由此可见,IIR数字滤波器设计的重要环节是模拟滤波器的设计。设计模拟低通滤波器的主要方法有Butterwort、Chebyshev、和椭圆等滤波器设计方法。
实验步骤
Butterwort数字滤波器设计
(1) Butterwort滤波器是通带阻带都单调衰减的滤波器。调用buttord函数可以确定巴特沃斯滤波器的阶数,其格式为:[N,Omegac]=buttord(Omegap,Omegas,Rp,As,’s’)。
其中,输入参数Rp,As分别为通带最大衰减和阻带最小衰减,以dB为单位;Omegap,Omegas分别为通带截止频率和阻带截止频率,‘s’说明所设计的是模拟滤波器。输出参数为滤波器的阶数,Omegac为3dB截止频率。
(2) 调用归一化巴特沃斯模拟原型滤波器的函数,其格式为
[z0,po,k0]=buttap(N)
其中,只要输入巴特沃斯滤波器的阶数,它就可以返回零点和极点数组z0、p0以及增益k0。当需要任意Omegac的未归一化的巴特沃斯滤波器时,就要用Omegac乘以p0或k0来进行归一化。
调用脉冲响应不变法或双线性变换法来设计巴特沃斯数字滤波器的函数,其格式分别如下。
= 1 \* GB3 ① 脉冲响应不变法:[bd,ad]=impinvar(b,a,Fs)。
其中,b为模拟滤波器分子系数向量;a为模拟滤波器分母系数向量;Fs为采样频率;
bd为数字滤波器的分子多项式系数;ad为数字滤波器的分母多项式系数。
= 2 \* GB3 ② 双线性变换法:[bd,ad]=bilinear(b,a,Fs),参数含义与impinvar一致。
Chebyshev = 1 \* ROMAN I型数字滤波器设计
Chebyshev = 1 \* ROMAN I型数字滤波器为通带纹波控制器:在通带呈现纹波特性,在阻带单调衰减。其格式如下:
[N,Omegac]= cheb1ord(Omegap,Omegas,Rp,As,’s’)
[z0,p0,k0]= cheb1ap(N,Rp)
参数含义与buttord buttap中参数一致。
Chebyshev = 2 \* ROMAN II型数字滤波器设计
Chebyshev = 2 \* ROMAN II型数字滤波器为阻带纹波控制器:在阻带呈现纹波特性,在通带单调衰减。其格式如下:
[N,Omegac]= cheb2ord(Omegap,Omegas,Rp,As,’s’)
[z0,p0,k0]= cheb2ap(N,Rp)
实验结果及分析:
例3-1
实验代码
wp=0.4*pi; %数字通带频率
ws=0.6*pi; %数字阻带频率
Rp=0.5; %通带波动(dB)
As=50; %阻带波动(dB)
Fs=1000; %置Fs=1000
OmegaP=wp*Fs; %原型通带频率
OmegaS=ws*Fs; %原型阻带频率
ep=sqrt(10^(Rp/10)-1); %通带波动参数
Ripple=sqrt(1/(1+ep*ep)); %通带波动
Attn=1/(10^(As/20)); %阻带衰减
%模拟巴特沃斯原型滤波器计算
[N,OmegaC]=buttord(OmegaP,OmegaS,Rp,As,s); %原型的阶数和截止频率的计算
[z0,p0,k0]=buttap(N); %归一化巴特沃斯原型设计函数
p=p0*OmegaC;z=z0*OmegaC;
%将零点极点乘以OmegaC, 得到非归一化零极点
k=k0*OmegaC^N; %将k()乘以Omegac^N,得到非归一化k
ba=k*real(poly(z)); %有零点计算分子系数向量
aa=real(poly
您可能关注的文档
- 数组应用,插入字符,简单家庭财务管理程序的设计.doc
- 酸碱平衡紊乱的分类及判断.doc
- 新深基坑围护结构钻孔灌注桩施工方案资料.doc
- 异步电动机与异步电动机的拖动实验.doc
- 临床试验术语解释.doc
- 庆云中心小学开展为树木穿冬衣活动.doc
- 住宅建设中常见质量问题成因分析及对策.doc
- 地方课教案(五下).doc
- 会员OA办公系统二次开发项目建设方案书.doc
- 物业电梯维保招标文件.doc
- 原电池电动势的测定实验报告.pdf
- 与业主、设计、总包、监理和他承包人的配合措施.pdf
- 公司管理流程.pptx
- 2024_2025学年新教材高中地理第1章地球的运动素养综合训练新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第3章大气的运动第1节常见天气系统第1课时锋与天气分层作业新人教版选择性必修1.doc
- 2024_2025学年新教材高中地理第1章地球的运动第2节地球运动的地理意义第4课时正午太阳高度的变化四季更替和五带划分分层作业课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第2章地表形态的塑造第2节构造地貌的形成第1课时地质构造与地貌课件新人教版选择性必修1.pptx
- 2024_2025学年新教材高中地理第1章地球的运动问题研究人类是否需要人造月亮课件新人教版选择性必修1.pptx
- 五片小雪花课件.pdf
- 2024_2025学年新教材高中地理第3章大气的运动第2节气压带和风带第1课时气压带和风带的形成分层作业课件新人教版选择性必修1.pptx
文档评论(0)