MATLAB中FIR滤波器阶数确定方法咨询：低通等波纹滤波器设计

在MATLAB中，对应IIR巴特沃斯滤波器的buttord函数，FIR等波纹滤波器（基于Parks-McClellan算法）可以通过firpmord函数计算满足指标的最小阶数和修正后的频率参数，再配合firpm生成滤波器系数，全程无需依赖图形化工具。

核心步骤与代码示例 #

假设你要求的归一化频率ω₁=0.3 rad/s（需转换为MATLAB的归一化频率，即Ω/π，范围0~1），具体实现如下：

1. 定义设计参数

   • 通带边缘频率Wp：需小于阻带边缘，示例设为0.08（对应0.08π rad/s，可按需调整）
   • 阻带边缘频率Ws：将给定的ω₁转换为归一化值，即0.3/pi
   • 通带最大纹波Rp：示例设为0.1 dB（可根据教学需求调整）
   • 阻带最小衰减Rs：按要求设为80 dB

2. 用firpmord计算阶数与修正参数
   该函数会自动计算满足指标的最小滤波器阶数，以及修正后的频率点、幅值目标和权重：

   % 定义核心参数
   Wp = 0.08;       % 归一化通带边缘
   Ws = 0.3/pi;     % 归一化阻带边缘（对应0.3 rad/s）
   Rp = 0.1;        % 通带最大纹波(dB)
   Rs = 80;         % 阻带最小衰减(dB)
   
   % 计算滤波器阶数和设计参数
   [N, Fo, Ao, W] = firpmord([Wp Ws], [1 0], [0.01*(10^(Rp/20)-1) 10^(-Rs/20)]);
   

   注：firpmord的第三个输入是线性幅值容差，由dB指标转换而来：通带容差为0.01*(10^(Rp/20)-1)，阻带容差为10^(-Rs/20)。

3. 用firpm生成滤波器系数
   基于上述参数生成等波纹FIR滤波器的系数：

   b = firpm(N, Fo, Ao, W);
   

4. 验证滤波器性能
   用freqz绘制幅频响应，确认阻带衰减是否达标：

   freqz(b, 1, 1024);
   title('FIR等波纹低通滤波器幅频响应');
   grid on;
   

关键说明 #

• firpm是实现Parks-McClellan算法的专用函数，对波纹和衰减指标的控制精度高于窗函数法
• 若使用窗函数法（如fir1），也可通过firpmord先获取阶数，再搭配凯泽窗等合适窗型，但等波纹法更契合你的指标需求
• 调整Wp时需确保Wp < Ws，否则参数逻辑无效

内容的提问来源于stack exchange，提问作者Boat01