首页理论教育传统方法FIR数字滤波器设计的MATLAB实现函数及性能分析

传统方法FIR数字滤波器设计的MATLAB实现函数及性能分析

2023-06-23 理论教育版权反馈

【摘要】：1）h=h，且N为奇数。图2-3 奇长度偶对称h设计FIR数字滤波器特性2）h=h，且N为偶数。MATLAB实现如下：式中，b=2hn=1，2，…图2-6 偶长度奇对称h设计FIR滤波器特性综合以上4种条件对应的FIR滤波器设计，给出这种传统方法的MATLAB实现函数[Hg，φω]=traFIR（h，p）。此函数包含两个输入参数即单位冲击响应h及滤波器类型p，函数的输出为系统的幅频特性和相频特性。此外，此方法得到的滤波器在通带、阻带存在幅度较大的波动，总体性能比较低下。

FIR数字滤波器的输入和输出关系是

它所对应的频率响应为

式中，H_g（ω）为H（e^jω）的增益，是一个可正可负的实函数；φ（ω）为H（e ^jω）的相频特性。具有线性相位系统的相位满足：

φ（ω）=β-αω （2-4）

式中，α、β为常数。具有线性相位的系统可保证在通频带内各种频率的信号都具有线性延迟，这样滤波后的信号就不会发生相位失真。这在诸如数据传输、语音综合等很多领域中都是非常希望得到的特性。

把式（2-3）两端用余弦函数展开有

把上面两式对应相除再交叉相乘展开，并利用三角恒等式得

式（2-6）恒成立的两组解为

当N分别取偶数和奇数时，FIR数字滤波器的频率特性有四种不同的形式。

1）h（n）=h（N-1-n），且N为奇数。

式中，

MATLAB实现如下：

运行得到的结果如图2-3所示。

图2-3 奇长度偶对称h（n）设计FIR数字滤波器特性

2）h（n）=h（N-1-n），且N为偶数。

式中，b（n）=2h（N/2-n）n=1，2，…，N/2。

MATLAB实现如下：

运行得到的结果如图2-4所示。

图2-4 偶长度偶对称h（n）设计FIR数字滤波器特性

3）h（n）=-h（N-1-n），且N为奇数。

式中，

MATLAB实现如下：

运行得到的结果如图2-5所示。

图2-5奇长度奇对称h（n）设计FIR数字滤波器特性

4）h（n）=-h（N-1-n），且N为偶数。

式中，

MATLAB实现如下：

运行得到的结果如图2-6所示。

图2-6 偶长度奇对称h（n）设计FIR滤波器特性

综合以上4种条件对应的FIR滤波器设计，给出这种传统方法的MATLAB实现函数[H_g，φω]=traFIR（h，p）。此函数包含两个输入参数即单位冲击响应h及滤波器类型p，函数的输出为系统的幅频特性和相频特性。该函数的MATLAB实现如下：

其中函数chkpara用于检测输入参数的合法性，这里不再赘述。

传统方法设计FIR滤波器需要知道系统的单位冲击响应序列，需要变换到频域才可以了解滤波器是低通或高通，这给设计带来了不方便。此外，此方法得到的滤波器在通带、阻带存在幅度较大的波动，总体性能比较低下。