(2)求Hd(e
jw
)的傅立叶反变换:
hd(n)
1 jwjwn
Hd(e)edw (12) 2
(3)根据技术要求(在通带 p处衰减不大于k1;在阻带 s处衰减不小于k2),确定窗函数形式w n 。并且根据采样周期T,确定相应的数字频率 p p, S ST。
(4)确定滤波器长度N。滤波器长度可以根据Hd(e波器的过度带宽有关。
(5)求所得设计滤波器的单位采样响应h n :
h(n) hd(n)w(n) (13)
jw
)的相位特性来确定,也和滤
(6)考察Hejw的指标:
H(ejw) h(n)e jwn (14)
n 0
N 1
(7)审核技术指标是否已经满足。如不满足,则重新选取较大的N进行(5),(6)计算;如果满足有余,则选取较小的N进行(5),(6)项计算。
设计完成后对已设计的滤波器的频率响应要进行校核,要得到幅频响应特性,运算量也是很大的。平时所要设计的数字滤波器,阶数和类型并不一定是完全给定的,很多时候都是要根据设计要求和滤波效果不断的调整,以达到设计的最优。
3.2.2FIR数字带通滤波器的设计
由于IIR数字滤波器利用模拟滤波器的理论进行设计,保留了模拟滤波器优良的幅度特性,因而设计中只考虑了幅度特性,没有考虑相位特性,所以一般情况IIR数字滤波器是非线性的[7]。而FIR数字滤波器很容易做到严格意义的线性相位特征,这对于信号处理是非常重要的。此外FIR数字滤波器单位冲激响应是有限长的,所以它永远是稳定的。稳定和线性是FIR数字滤波器最突出的两个优点[3]。
MATLAB提供了许多常用的窗函数,其中布莱克曼窗函数的调用形式为:
W=Blackman(N)
其中N是窗函数的长度,返回的变量w是一个长度为N的列向量,给出窗函数在N点的取值。
根据设计任务,用MATLAB设计的FIR数字带通滤波器的具体程序如下:
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库基于MATLAB的数字滤波器设计与仿真(很全--免费下载)(12)在线全文阅读。
相关推荐: