原创 低通滤波的效果不理想

http://www.chinavib.com/forum/thread-71992-1-41.html

 

处理之前的信号如下：
放大图

附件

以前对信号处理没什么接触，也不知道有没有更合适的方法来解决这个问题，只是自己的一个尝试，所以如果各位有什么想法敬请直言——是不是参数设置不合理，或者就根本不应该用低通滤波，而应该用其他更合适的方法。
稍微补充一点，就是如下图中的红线就是想要的结果——那个是我用鼠标画上去的。

--------------------------------- 版主的问题不是很容易解决的，因为实际上并没有低频信号存在，所以不是能用低频滤波器能过滤出来的。版主所要的波形幅值在0.985~1.015之间，我把“粗线”中的一部分取出，并作了谱分析（版主没有给采样频率，我假设为1000Hz），从图中看出这“粗线”中的主要频率在250Hz（见下图）。
我在处理中对signal信号进行了截幅，只取幅值在0.985~1.015之间的波形（见图），并通过250Hz的带陷滤波器，滤波器的输出（见下图）比较接近版主的要求，我想再做些平滑处理就能达到要求了。
程序为：
load signal.mat;
fs=1000;
x=signal';
N=length(x);
for k=1 : N   %截幅
    if x(k)>1.015
        x(k)=1.015;
    elseif x(k)<0.985
        x(k)=0.985;
    end
end
x=x-mean(x);
plot(x);        
t=(0:N-1)/fs;
%带陷滤波
x1=cos(2*pi*250*t/fs);
x2=sin(2*pi*250*t/fs);
u=0.05; w1=0.1; w2=0.1;
e=zeros(1, N);
for i=1: N
  y(i)=w1 * x1(i)+ w2 * x2(i);
  e(i) =x(i)-y(i);
  w1=w1+u * e(i) * x1(i);
  w2=w2+u * e(i) * x2(i);
end
figure
plot(y); ylim([-0.015 0.015]);
title('带陷滤波器输出(在范围-0.015~0.015之间)');

-------------------------------