我正在尝试以这种方式实现低通,如Digital signal processing guide
书中所述,但我alfa = 2.660393305517957e-55
的显然是错误的。这是代码:
x = exp(-2*pi*20);
a0 = 1 - x;
b1 = x;
eda2(1) = eda_raw2(1);
for n = 2:length(eda_raw2)
eda2(n) = a0*eda_raw2(n) + b1.*eda2(n - 1);
end
我正在尝试以这种方式实现低通,如Digital signal processing guide
书中所述,但我alfa = 2.660393305517957e-55
的显然是错误的。这是代码:
x = exp(-2*pi*20);
a0 = 1 - x;
b1 = x;
eda2(1) = eda_raw2(1);
for n = 2:length(eda_raw2)
eda2(n) = a0*eda_raw2(n) + b1.*eda2(n - 1);
end
问题是您如何计算x
离散时间滤波器的极点 ( )。您没有考虑采样频率!计算应为x = exp(-2*pi*20/fm);
,其中fm
是以 Hz 为单位的采样频率。
查看在线数学作品页面:低通滤波器设计