- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
滤波器设计步骤实现程序
数字滤波器的设计步骤及程序实现
湖南理工学院信息与通信工程学院
一、IIR脉冲响应不变法设计步骤
1、已知实际数字指标
2、将数字指标化为原型模拟指标,可设T=pi,
3、求原型模拟滤波器的,其中:
4、根据N写出归一化原型系统函数
5、用代入得原型系统函数
6、将化为部分分式展开形式
7、写出的极点,并写出的部分分式展开形式
8、将化为分子分母形式,验证设计结果。
二、IIR双线性变换法设计步骤
1、已知实际数字指标
2、将数字指标化为原型模拟指标,可设T=2,
3、求原型模拟滤波器的,其中:
4、根据N写出归一化原型系统函数
5、用代入得原型系统函数
6、用代入原型系统函数得
8、将整理成分子分母形式,验证设计结果。
三、FIR窗函数法设计步骤
1、已知实际数字指标
2、根据选窗的类型:矩形窗as21dB, A=1.8π,窗函数是 boxcar(N);三角窗as25dB, A=6.1π,窗函数是 bartlett(N);汉宁窗as44dB, A=6.2π,窗函数是 hanning(N);哈明窗as53dB, A=6.6π,窗函数是 hamming(N);布莱克曼窗as74dB, A=11π,窗函数是 blackman(N)。
3、根据过渡带和窗类型求总点数。
4、根据写出理想频响指标
5、根据算出
6、对加窗得设计结果
8、写出,验证设计结果。
四、FIR频率采样法设计步骤
1、已知实际数字指标
2、根据选过渡带点数
3、根据过渡带和过渡带点数求总点数。
4、根据求出,设置过渡值
5、根据约束条件构建理想频响的采样指标
6、对进行IDFT变换得,取实部。
7、写出,验证设计结果,优化过渡值大小、过渡点位置和过渡点多少。
一、IIR滤波器设计:脉冲响应不变法实现程序
%用脉冲响应不变法设计butterworth数字低通滤波器
%技术指标:wp=0.3*pi rad, ap=2dB, ws=0.5*pi rad, as=10dB
clc; clear; close all; format compact;%程序初始化
wp=0.3*pi, ap=2, ws=0.5*pi, as=10,%输入数字指标
T=pi,%假设采样周期,用于设计原型模拟滤波器,不影响H(z)的设计结果
Wp=wp/T, Ap=ap, Ws=ws/T, As=as,%将数字指标转化为原型模拟指标
M=log10( (10 .^ (0.1*Ap) - 1)./(10 .^ (0.1*As) - 1) ) / ...
(2*log10(Wp/Ws)) ,%计算滤波器阶数
N = ceil( M),%滤波器阶数向上取整
Wcp = Wp / ( (10^(.1*Ap) - 1)^(1/(2*N))),%通带边界精确满足的截止频率
Wcs = Ws / ( (10^(.1*As) - 1)^(1/(2*N))),%阻带边界精确满足的截止频率
Wc=Wcp,%截止频率用通带边界精确满足的截止频率
%Wc=(Wcp+Wcs)/2,%通带阻带边界都有余量的截止频率
%Wc=Wcs,%截止频率用阻带边界精确满足的截止频率
[bp,ap]=butter(N,1,s),%求归一化原型滤波器系统函数Ga(p)P157
tf(bp,ap,variable,p),%显示Ga(p)
[bs,as]=lp2lp(bp,ap,Wc),%去归一化得原型滤波器系统函数Ha(s)
tf(bs,as),%显示Ha(s),分子不足前面补0
[Ak,sk]=residue(bs,as),%将Ha(s)按部分分式形式展开
ak=T*Ak,zk=exp(sk*T),%将Ha(s)的部分分式参数转换为H(z)的部分分式参数
[bz,az]=residuez(ak,zk,0),%将H(z)的部分分式形式化为分子分母等阶形式
tf(bz,az,variable,z^-1),%显示系统函数
%[bz1,az1] = impinvar(bs,as,1/T)%调用impinvar函数验证
%tf(bz1,az1,variable,z^-1),%显示验证系统函数,z^-1式的分子不足是后面补0
freqz(bz,az,100),%绘出频率特性曲线,检验设计指标
二、IIR滤波器设计:双线性变换法实现程序
%用双线性变换法设计butterworth数字低通滤波器
%技术指标:wp=0.3*pi rad, ap=2dB, ws=0.5*pi rad, as=10dB
clc; clear; close all; format compact;%程序初始化
wp=0.3*pi, ap=2,
您可能关注的文档
- 深刻理解mai函数的命令行参数.doc
- 深圳七下暑假作答案(五科).doc
- 深圳中考地理知(2013年第1版)点.doc
- 深圳国际园林花博览园1MWp并网光伏电站系统描述及效益2004_10_19.doc
- 深圳市历年中考语试题命题规律回顾.doc
- 深圳市宝荣娱乐2014年消防演练报告.doc
- 深水射流曝气机射流式潜水曝气机价格.docx
- 淳化县苹果种植术普及现状.doc
- 混凝土立方体试劈裂抗拉强度尺寸效应研究.doc
- 混凝土结构设计卷一B卷.doc
- 第十一章 电流和电路专题特训二 实物图与电路图的互画 教学设计 2024-2025学年鲁科版物理九年级上册.docx
- 人教版七年级上册信息技术6.3加工音频素材 教学设计.docx
- 5.1自然地理环境的整体性 说课教案 (1).docx
- 4.1 夯实法治基础 教学设计-2023-2024学年统编版九年级道德与法治上册.docx
- 3.1 光的色彩 颜色 电子教案 2023-2024学年苏科版为了八年级上学期.docx
- 小学体育与健康 四年级下册健康教育 教案.docx
- 2024-2025学年初中数学九年级下册北京课改版(2024)教学设计合集.docx
- 2024-2025学年初中科学七年级下册浙教版(2024)教学设计合集.docx
- 2024-2025学年小学信息技术(信息科技)六年级下册浙摄影版(2013)教学设计合集.docx
- 2024-2025学年小学美术二年级下册人美版(常锐伦、欧京海)教学设计合集.docx
文档评论(0)