数字信号处理第七章演示文稿.pptVIP

* 大连民族学院 * 八 频率抽样法的优缺点 优点：设计方法简单，在频域直接进行设计，适合优化设计，可利用FFT算法，通过改变N和设置过渡带样值，能获得满意结果； 缺点：不容易控制边缘频率点的频响值。 * 大连民族学院 * 九 IIR与FIR数字滤波器的比较（P221） 相同的技术指标下，IIR由于存在反馈，因此比FIR滤波器的阶数低，存储单元、运算次数少，经济。 FIR容易实现线性相移；IIR一般不易，除非增大阶数。 FIR主要采用非递归结构和算法，因而在有限精度的运算中总是稳定的，由于无反馈，所以无累积误差；IIR必须采用递归结构，由于四舍五入和有限精度运算，可能会引起寄生振荡。 FIR的h(n)为有限长，可用FFT算法，IIR一般不能用。 从设计角度，IIR可以利用模拟滤波器的结果转换；FIR则无现成公式可用，现有方法也不易控制增益、衰减和边界频率。 IIR主要实现规格化的、频率特性为分段常数的标准低通、高通、带通、带阻、全通等；FIR比较灵活，可以设计出理想的正交变换器、理想微分器等，应用较广。 * 大连民族学院 * 汉宁窗的幅度特性 * 大连民族学院 * 海明（Hamming）窗（改进的升余弦） 主瓣宽8π/N，旁瓣比Hanning窗更小。 * 大连民族学院 * 布莱克曼（Blackman）窗（二阶升余弦） 主瓣宽12π/N。 * 大连民族学院 * 常用的窗函数 * 大连民族学院 * 常用窗函数的幅度特性 (a)矩形窗； (b)巴特利特窗(三角形窗)； (c)汉宁窗； (d)哈明窗； (e)布莱克曼窗 * 大连民族学院 * 理想低通加窗后的幅度特性(N=51,ωc=0.5π) (a)矩形窗； (b)巴特利特窗(三角形窗)； (c)汉宁窗； (d)哈明窗； (e)布莱克曼窗 * 大连民族学院 * 凯泽（Kaiser）窗 其中I0(x)是第一类修正零阶贝塞尔函数，其幂级数展开式为： ?称形状参数，用以调整主瓣宽度和旁瓣电平。选择?可使Kaiser窗变为其它窗。 * 大连民族学院 * 表7.2 凯塞窗参数对滤波器的性能影响 * 大连民族学院 * 几种窗函数的参数 窗名称 3dB 带宽 主瓣 过零点 宽度 最大旁瓣峰值A 旁瓣谱峰渐近衰减速度D 阻带 最小衰减 矩形窗 0.89 ?? 4?/N －13dB -6dB/OCT －21dB Bartlett （三角） 1.28 ?? 8?/N －27dB -12dB/OCT －25dB Hanning （汉宁） 1.44 ?? 8?/N －32dB -18dB/OCT －44dB Hamming （海明） 1.3 ?? 8?/N －43dB -6dB/OCT －53dB Blackman （布莱克曼） 1.68 ?? 12?/N －58dB -18dB/OCT －74dB * 大连民族学院 * 三 窗函数法的设计步骤 选定理想的频率响应函数Hd(ejω)； 由Hd(ejω)求出理想单位样值响应hd(n)； 根据允许的过渡带宽与阻带最小衰减的要求，选择窗函数w(n)的形式，并估计窗口长度N； 求出h(n)=hd(n)w(n)，若要求线性相位，则hd(n)和w(n)均对(N-1)/2对称； 求出H(ejω)＝DTFT[h(n)]，检验是否满足设计要求，若不满足，则需重新设计。 * 大连民族学院 * 例题1 设计一个线性相位FIR低通数字滤波器，所希望的频率响应Hd(ej?)在0 ≤|?|≤0.25?之间为1，在0.25?≤|?|≤?之间为0，分别取N＝11,21,41，观察其幅频响应的特点。 解：由题意要求，得 * 大连民族学院 * 例题1续 当N＝11时，求得　 h(0) = h(N-1) = h(10) = -0.045 h(1) = h(N-2) = h(9) = 0 h(2) = h(N-3) = h(8) = 0.075 h(3) = h(N-4) = h(7) = 0.1592 h(4) = h(N-5) = h(6) = 0.2251 h(5) = 0.25 0 0.1 0.2 0.3 0.4 0.5 H(ej?)| N=11 N=21 N=41 吉布斯(Gibbs)现象 右图分别给出了N＝11,21及41时的|H(ej?)|的幅频特性曲线。 * 大连民族学院 * 小结 窗函数设计方法是选取某种好的窗函数对hd(n)予以截短来得到h(n)的，所以称之为“窗函数法”。 FIR－DF设计的窗函数法不像IIR－DF的设计那样能精确指定通带、阻带的边缘频率?p、?s，也不能精确地给定在这两个带内的衰减?p、?s，而是仅给出通带截止频率?p（或?c），其它几个参数是靠h