离散傅里叶变换DF.ppt

* 例7.26 窗函数法设计带通滤波器 使用fir1函数b = fir1(N, wc, window) 编程 参数?c为行向量?c = [?lp/?，?hp/?] 根据阻带最小衰减Rs = 60dB选择窗函数类型和阶次。可以查上面列出的“窗函数设计滤波器时的阶数选择表”。选blackman窗，其滤波器阻带最小衰减可达到74dB，其窗口长度M由过渡带宽度B = 0.15? 决定，Blackman窗设计的滤波器过渡带宽度为12?/M，故M取80。因M = N+1，所以滤波器阶数N = 79。 * 例7.27 用remez函数设计低通滤波器 解：先由题意计算设计参数 f = [1/4,5/16]，m = [1,0]; dev的计算稍复杂一些，由于 所以 有了这几个参数就可以调用remezord和remez函数了. * 例7.27 用remez函数设计低通滤波器 横线为-3dB，两条竖线分别位于频率?/4和5?/16。显然，通带指标稍有富裕，过渡带宽度和阻带最小衰减刚好满足指标要求。 程序输出的幅频特性 * 例7.28 remez函数设计高通滤波器 观察等波纹逼近法中加权系数w(?)及滤波器阶数N的作用和影响。期望逼近的滤波器通带为[3?/4，?]，阻带为[0，23?/32]。 解:在滤波器设计中，技术指标越高，实现滤波器的阶数也就越高。另外，对固定的阶数，通带与阻带指标可以互换，过渡带宽度与通带波纹和阻带衰减指标可以互换。 取f = [0, 3/4, 23/32, 1]，m = [0, 0, 1, 1]。其余参数分三种情况进行设计，①N = 30，w = [1, 1]；②N = 30，w = [1, 5]; ③N = 60，w = [1, 1]。 * 例7.28 remez函数设计高通滤波器 程序运行结果如图 由图可见，w较大的频段逼近精度较高；w较小的频段逼近精度较低。N较大时逼近精度较高,N较小时逼近精度较低 。 * 7.6 IIR数字滤波器设计 IIR数字滤波器设计的主要方法是先设计低通模拟滤波器，进行频率变换，将其转换为相应的（高通、带通等）模拟滤波器，再转换为高通、带通或带阻数字滤波器。对设计的全过程的各个步骤，MATLAB都提供了相应的工具箱函数，使IIR数字滤波器设计变得非常简单。本节主要结合例题介绍这些IIR滤波器设计的工具箱函数。 IIR数字滤波器的设计步骤由以下的流程图来表示。下面以巴特沃斯滤波器设计函数为典型，介绍此流程图中函数的功能和用法。 * IIR数字滤波器设计流程图 模拟低通滤波器原型设计 Buttap,cheb1ap,cheb2ap besselap，ellipap函数 频率变换(变为高通，带通，带阻等) lp2lp，lp2hp，lp2bp，lp2bs 模拟数字变换 bilinear impinvar 合为一步的设计函数 butter，cheb1，cheb2，ellip， besself 求最小阶数N buttord, cheb1ord cheb2ord,ellipord 设计指标 wp, ws, Rp, Rs N，wc 滤波器系数 B，A 滤波器系数 B，A * 巴特沃斯滤波器设计流程 （1）求最小阶数N的函数buttord [N, wc] = buttord (wp, ws, Rp, Rs, ‘s’) 根据滤波器指标wp，ws，Rp，Rs，求出巴特沃斯模拟滤波器的阶数N及频率参数wc，此处wp，ws及wc均以弧度/秒为单位。 （2）得到N后,调用设计函数buttap [z，p，k] = buttap(N) 得到[z, p, k]后，很容易求出滤波器系数B，A。 （3）调用模拟频率变换函数lp2lp [Bt, At] = lp2lp(B, A, wo) （4）调用模拟数字变换函数 ? [Bd, Ad] = bilinear (B, A, Fs) * 集成的数字滤波器设计函数 把（2）、（3）、（4）合为一步的数字滤波器设计函数butter(N, wc, ftype) ? [B, A] = butter (N, wc) 设计低通或带通数字滤波器系数B，A（当为带通滤波器时，第(1)类函数由wp = [wp1, wp2]会自动生成wc = [w1, w2]）。 ? [B, A] = butter (N, wc, high) 设计高通数字滤波器系数B，A。 ? [B, A] = butter (N, wc, stop) 设计带阻数字滤波器系数B，A。 butter(N, wc, ftype)还有零极增益和状态空间形式，读者可用help命令查阅。 * 例7.29 巴特沃斯模拟滤波器设计 设计一个低通巴特沃斯