新数字信号处理 教学课件 尹为民 25 6.3FIR DF设计(二).pptVIP

§6.3 FIR数字滤波器设计 引言 窗函数设计法 频率抽样设计法 设计FIR滤波器的最优化方法 IIR与FIR数字滤波器的比较 §6.3 FIR数字滤波器设计 review 三、 频率抽样设计法 1 频率抽样法基本原理 1 频率抽样法基本原理 2 线性相位约束 2 线性相位约束 2 线性相位约束 2 线性相位约束 2 线性相位约束 3 频率抽样法设计线性相位FIR DF的步骤 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 4 过渡带抽样的优化设计 4 过渡带抽样的优化设计 4 过渡带抽样的优化设计 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 利用1个过渡点设计的FIR低通滤波器的增益响应 利用2个过渡点设计的FIR低通滤波器的增益响应 频率抽样法的优、缺点 四、 FIR数字滤波器的优化设计 1 问题的提出 1 问题的提出 1 问题的提出 2 等波纹FIR滤波器设计 2 等波纹FIR滤波器设计 2 等波纹FIR滤波器设计 2 等波纹FIR滤波器设计 2 等波纹FIR滤波器设计 2 等波纹FIR滤波器设计 3 优化设计的MATLAB实现 五、 IIR与FIR数字滤波器的比较 五、 IIR与FIR数字滤波器的比较 重点：了解FIR数字滤波器优化设计的基本思想及基本方法。 问题的提出 等波纹FIR滤波器设计 优化设计的MATLAB实现 矩形窗函数法设计出的FIR滤波器是方均误差最小意义下的最优FIR滤波器。 等于零，即h(n)= hd(n), 0?n ?N-1 最小 方均误差 由Parseval定理 解决方法： 采用最大误差最小化准则优化设计！ 问题： 通带和阻带波动相等，无法单独控制； 通带或阻带内波动不均匀，峰值较大 最大最小准则(minimax criterion)或Chebyshev准则 其中H d(?)为所需系统的幅度函数， H (?)为设计所得系统的幅度函数，W(?)为加权函数，且W(?)?0 。 FIR数字滤波器优化设计的基本思想： 在一定的误差准则下，设计线性相位FIR滤波器H(z)，使得H(z)所对应的幅度函数H(?)和理想滤波器的幅度函数Hd(?)的误差在区间上达到最小。 为使四种类型的线性相位滤波器能有统一的优化算法，将它们的幅度函数统一表示为 设计方法：利用Chebyshev误差准则设计FIR滤波器 等波纹FIR滤波器的设计方法，由 Parks-McClellan 在1972提出，在解决问题过程中使用了数学优化中的Remez交换算法，故称其为PM算法或Remez算法。 设计出的FIR滤波器的特点： 幅度响应在通带和阻带都会呈现等波波纹的波动。 设计准则：确定系数{?(n)}，使得误差 达到最小。 其中 交替定理 最佳逼近函数P(?)存在并惟一的充要条件为在集合A的范围内至少存在r+1个极值点 ?1 ? 2 ? ? r+1，使得 e (?i ) = -e(? i+1 )， 1 ? i ? r, |e (?i )| = e, 1 ? i ? r 如果选择通带的加权为1，阻带的加权为dp/ds，则设计得到的最优滤波器通带和阻带的波动分别为 如果r+1个极值点已找到，由交替定理 共有r+1个方程，从中可解出r+1个系数{?(n), 0 ? n ? r}和误差e 。 等波纹线性相位FIR滤波器的设计步骤： (1) 用Kaiser提出的经验公式估计滤波器的阶数 并由N及类型(I型、II型)确定r。 (2) 选定初始极值点{?i ; 1?n?r+1}。 (3) 计算误差e。 (4) 计算函数P(?)在A中的抽样值。 (5) 寻找新的极值点 ，并求出最大误差值 。 等波纹线性相位FIR滤波器的设计步骤： (7) 解方程获得?(n)。 (8) 由?(n)求出h(n) 。 (6) 如果