我正在使用remez功能编写脚本以使用Octave将高通滤波器应用于ECG信号。基本上,代码是这样的:
fs = 500; % sampling frequency
nyq = fs/2; % Nyquist frecuency
% remez parameters
N = 256; % filter order
f1 = 0.2; % begin of transition band in hz
f2 = 4; % end of transition band in hz
F = [0 f1/nyq f2/nyq 1];
A = [0 0 1 1];
hc1 = remez(N,F,A);
xf1 = filter(hc1,1,data_vector);
如果我绘制结果,我会得到:
所以,我的问题是看起来高通滤波器运行良好,但信号有偏移。我希望看到信号以零为中心(没有直流偏移)。
这种抵消的原因是什么?
PD:在该图中我没有包含前N + 2个值,因为它们是过滤器的初始化值。
谢谢。