《数字信号处理》课程设计基于切比雪夫I型的IIR数字高通滤波器设计Word文件下载.docx
- 文档编号:13415293
- 上传时间:2022-10-10
- 格式:DOCX
- 页数:21
- 大小:478.38KB
《数字信号处理》课程设计基于切比雪夫I型的IIR数字高通滤波器设计Word文件下载.docx
《《数字信号处理》课程设计基于切比雪夫I型的IIR数字高通滤波器设计Word文件下载.docx》由会员分享,可在线阅读,更多相关《《数字信号处理》课程设计基于切比雪夫I型的IIR数字高通滤波器设计Word文件下载.docx(21页珍藏版)》请在冰豆网上搜索。
3.2源程序代码--------------------------------------------------------------14
3.3结果分析------------------------------------------------------------------17
(四)总结和心得体会--------------------------------------------20
(五)参考文献-----------------------------------------------------21
基于切比雪夫I型的IIR数字高通滤波器设计
(一)数字滤波器的概述
1.1数字滤波器的设计方法
数字滤波器的设计方法有多种,如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。
1.2数字滤波器的性能要求
我们在进行滤波器设计时,需要确定其性能指标。
一般来说,滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。
以低通滤波器特性为例,频率响应有通带、过渡带及阻带三个范围。
如图1-1所示:
在通带内:
在阻带中:
其中
为通带截止频率,
为阻带截止频率,
为通带最大衰减,
为阻带最大衰减。
图1-1低通滤波器的幅频特性指标示意图
1.3数字滤波器的技术要求
滤波器技术要求主要包括4个方面。
即:
滤波器的截止频率
低通滤波器的截止频率主要包括通带截止频率(又称通带上线频率)
和阻带下限截止频率
;
高通滤波器的截止频率主要包括通带截止频率(下限频率)
和阻带上限截止频率
带通滤波器的截止频率主要包括通带下限截止频率
,通带上限截止频率
,下限通带截止频率
,以及上阻带截止频率
带阻滤波器的截止频率与带通滤波器一致,也主要包括通带下限截止频率
,下阻带截止频率
。
带通带阻的容限
滤波器中带通带阻的容限
与
的具体技术指标,往往由允许的最大衰减
及阻带应达到的最小衰减
给出。
通带及阻带的衰减
,
分别定义为:
式中均假定
已被归一化为1.例如当
在
处下降为0.707时,
,在
处降到0.01时,
.
(二)基于切比雪夫I型无限脉冲响应IIR数字高通滤波器的设计依据和原理
2.1课设任务
设计一个数字高通滤波器,要求通带截止频率fp=3000Hz,通带衰减不大于1dB(rp=1dB);
阻带截止频率fr=1500Hz,阻带衰减不小于50dB(rs=50dB)。
采样频率fs=15000Hz,采用切比雪夫I型IIR滤波器。
2.2IIR数字滤波器
IIR数字滤波器设计原理
IIR数字滤波器的设计步骤是:
(a)按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标;
(b)根据转换后的技术指标设计模拟低通滤波器
(c)在按一定规则将
转换为
若所设计的数字滤波器是低通的,那么上述设计工作可以结束,若所设计的是高通、带通或者带阻滤波器,那么还有步骤:
(d)将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标,然后按上述步骤(b)设计出模拟低通滤波器
,再由冲击响应不变法或双线性变换将
转换为所需的
s-z映射的方法有:
冲激响应不变法、阶跃响应不变法、双线性变换法等。
下面讨论双线性变换法。
双线性变换法是指首先把s平面压缩变换到某一中介平面s1的一条横带(宽度为
即从
到
),然后再利用
的关系把s1平面上的这条横带变换到整个z平面。
这样s平面与z平面是一一对应关系,消除了多值变换性,也就消除了频谱混叠现象。
s平面到z平面的变换可采用
(2-5)
(2-6)
令
有:
(2-7)
从s1平面到z平面的变换,即
(2-8)
代入上式,得到:
(2-9)
一般来说,为使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系,可引入代定常数c,
(2-10)则
(2-11)
这种s平面与z平面间的单值映射关系就是双线性变换。
有了双线性变换,模拟滤波器的数字化只须用进行置换。
2.3由模拟滤波器设计IIR数字滤波器
理想的滤波器是非因果的,即物理上不可实现的系统。
工程上常用的模拟滤波器都不是理想的滤波器。
但按一定规则构成的实际滤波器的幅频特性可逼近理想滤波器的幅频特性,例如巴特奥兹(Butterworth)、切比雪夫(Chebyshev)滤波等。
切比雪夫滤波器的原理:
特点:
误差值在规定的频段上等幅变化。
巴特沃兹滤波器在通带内幅度特性是单调下降的,如果阶次一定,则在靠近截止频率
处,幅度下降很多,或者说,为了使通常内的衰减足够小,需要的阶次
很高,为了克服这一缺点,采用切比雪夫多项式逼近所希望
切比雪夫滤波器的
在通带范围内是等幅起伏的,所以同样的通带衰减,其阶数较巴特沃兹滤波器要小。
可根据需要对通带内允许的衰减量(波动范围)提出要求,如要求波动范围小于1db。
振幅平方函数为
(3-7)
式中
—有效通带截止频率
—与通带波纹有关的参量,
大,波纹大,
—N阶切比雪夫多项式,定义为
(3-8)
(3-9)
如图3-1,通带内
,变化范围
,随着
(迅速趋于零)
当
时,
(3-10)
N为偶数,
(3-11)
N为奇数,
(3-12)
有关参数的确定:
a.通带截止频率,预先给定;
b.由通带波纹表为
(3-13)
(3-14)
给定通带波纹值分贝数后,可求
(3-15)
c.阶数N—由阻带的边界条件确定。
(
A事先给定)
(3-16)
(3-17)
(3-19)
得
(3-20)
2.4数字高通滤波器的设计(本设计采用双线性变换法)
双线性变换法的原理:
双线性变换法是采用非线性频率压缩的方法,将整个频率轴上的频率范围压缩到
之间,再用
转换到Z平面上。
也就是说,第一步先将整个S平面压缩映射到S1平面的
一条横带里;
第二步再通过标准变换关系
将此横带变换到整个Z平面上去。
这样就使S平面与Z平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象,映射关系如图所示。
图1-2映射关系
为了将S平面的整个虚轴
压缩到S1平面
轴上的
段上,可以通过以下的正切变换实现
由
经过0变化到
由-∞经过0变化到+∞,也即映射了整个
轴。
将式(3-5)写成
将此关系解析延拓到整个S平面和S1平面,令
则得
再将S1平面通过以下标准变换关系映射到Z平面
从而得到S平面和Z平面的单值映射关系为:
(3-6)
(3-7)
式(3-6)与式(3-7)是S平面与Z平面之间的单值映射关系,这种变换都是两个线性函数之比,因此称为双线性变换
式(3-5)与式(3-6)的双线性变换符合映射变换应满足的两点要求。
首先,把
,可得
即S平面的虚轴映射到Z平面的单位圆。
其次,将
代入式(3-8),得
因此
由此看出,当
也就是说,S平面的左半平面映射到Z平面的单位圆内,S平面的右半平面映射到Z平面的单位圆外,S平面的虚轴映射到Z平面的单位圆上。
因此,稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的.
本设计总体设计步骤为:
(1)确定数字高通滤波器的技术指标
、
(2)将数字高通滤波器的技术指标转换成高通模拟滤波器的技术指标
,转换公式为
(3)利用频率变换
将模拟高通滤波器技术指标转换成归一化模拟低通滤波器G(p)的技术指标;
(4)设计模拟低通滤波器G(p),并去归一化得:
(5)采用双线性变换将模拟低通滤波器H(s)转换成数字低通滤波器H(z):
(6)采用频带变换,将数字低通滤波器转换成所需类型的数字高通滤波器。
(三)基于切比雪夫I型无限脉冲响应IIR数字高通滤波器的具体设计过程
3.1计算过程
(1)数字高通滤波器的技术指标:
(2)将数字高通滤波器的技术指标转换成模拟高通滤波器的技术指标(令T=2s):
(3)相应低通滤波器的技术指标:
(4)设置归一化低通模拟滤波器G(p)各项指标:
(5)将极点
和
代入
,求低通模拟滤波器G(p):
可得:
归一化:
(6)利用双线性变换将模拟低通传输函数
转换为数字低通滤波器系统函数
:
(7)用频带变换法将数字低通
转换成数字高通
3.2源程序代码
%%%%chebfilter.m
%%%%highpass
fp=3000;
rp=1;
fr=1500;
rs=50;
Fs=15000;
[N,Wn]=cheb1ord(2*fp/Fs,2*fr/Fs,rp,rs);
[b,a]=cheby1(N,rp,Wn,'
high'
);
[H,W]=freqz(b,a,512,'
whole'
n=1:
1000;
signal=8*sin(n*pi/6
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号处理 数字信号 处理 课程设计 基于 IIR 数字 滤波器 设计
![提示](https://static.bdocx.com/images/bang_tan.gif)