|
这段时间看了ARMFLY的DSP教程,并自己动手设计一个FIR滤波,我的本意是滤除一个振动干扰波形如下图1:
FIR低通滤波后效果图
请看我的原码:
x = [10019 10019 10019 10019 10019 10020 10020 10020 10017 10017 10018 10018 9899 9899 9899 9827 9827 10696 10696 10696 9811 9811 9857 9857 9834 9834 9834 10187 10187 10266 10266 9790 9790 9790 10020 10020 9999 9999 9999 9926];
blo = fir1(16,0.02);
outlo=filter(blo,1,x);
subplot(2,1,1);
Fs = 60;
t = (0:length(x)-1)/Fs;
plot(t,x)
title('Original Signal')
subplot(2,1,2)
plot(t,outlo)
title('Lowpass Filtered Signal')
xlabel('Time (s)')
我的X是原波形,受到振动干扰的一段数据,现在经过FIR低通滤除之后的效果,请问我这里有没有问题,滤波后的数据和原始数据有一定差别,实际我要得到的数据是10019
请教论坛大神给个方法怎么处理这个信号 |
|