如何确定巴特沃斯滤波器的阶数与截止频率？滤波失真问题咨询

我正尝试设计最优高通滤波器以去除信号中的毛刺。最初，我对5Hz正弦信号应用巴特沃斯（Butterworth）高通滤波器，具体代码如下：

% Sample rate and time vector
Fs = 1000;            % Sampling frequency (Hz)
t = 0:1/Fs:1;         % Time vector (1 second duration)

% Create a simple signal with a glitch
signal = sin(2*pi*5*t);               % A sine wave at 5 Hz
signal(500:510) = signal(500:510) + 2; % Introduce a glitch

% Define high-pass filter characteristics
cutoff_frequency = 4; % Cutoff frequency (Hz)
order = 8;            % Filter order

% Design and apply the high-pass filter
[b, a] = butter(order, cutoff_frequency / (Fs/2), 'high');
filtered_signal = filter(b, a, (signal));

当阶数为8、截止频率为4Hz时，滤波效果良好。但将截止频率调整至3Hz以下时，信号出现严重失真。我无法理解信号失真的原因，请问这是怎么回事？是否有解决该问题的建议？

失真原因 #

• 目标信号与截止频率过于接近：你用的目标信号是5Hz正弦波，当截止频率降到3Hz以下时，5Hz和截止频率的比值大幅缩小。巴特沃斯高通滤波器的通带边缘（靠近截止频率的区域）存在幅度衰减和相位偏移，高阶滤波器（8阶）的滚降特性虽陡，但信号频率越接近截止频率，这些影响就越显著，最终导致波形失真。
• IIR滤波的非线性相位影响：filter函数实现的是直接IIR滤波，巴特沃斯IIR滤波器本身是非线性相位的，高阶滤波器的相位偏移会更严重。当目标信号频率靠近截止频率时，相位失真被放大，直接表现为信号波形变形。

解决建议 #

• 调整滤波器参数：如果必须使用低截止频率，可降低滤波器阶数（比如降到4阶以内），减小相位偏移带来的影响；或者将截止频率设为目标信号频率的1/3以下（比如1.5Hz），让5Hz信号远离通带边缘，避免衰减和相位偏移的干扰。
• 改用零相位滤波：用filtfilt替代filter函数，它会对信号做正向、反向两次滤波，抵消IIR滤波器的相位失真，同时等效阶数翻倍，既能保留高阶滤波的陡滚降特性，又能避免波形失真。修改后的代码如下：

filtered_signal = filtfilt(b, a, signal);

• 换成FIR滤波器：FIR滤波器天生具有线性相位特性，不会产生相位失真。可以用fir1设计高通FIR滤波器，示例代码：

% 设计100阶FIR高通滤波器，截止频率2Hz
b = fir1(100, 2/(Fs/2), 'high');
filtered_signal = filtfilt(b, 1, signal);

虽然FIR需要更高阶数才能达到和IIR类似的滚降效果，但能保证目标信号的波形不失真。

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