求教:关于自适应陷波后的频谱图

xiaocheng_2007

感谢songzy41网友提供的程序,用该网友提供的程序,我对我的数据进行了处理,发现效果不错.程序和图如下所示.
load 'jd1.txt'
x=jd1;
N=length(x);
%自适应陷波器
x1=cos(2*pi*50*t);
x2=sin(2*pi*50*t);
u=0.1; 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
%画图
out=e;
plot(t,out,'r',t,x,'g');
hold off;
除此之外,我还有一个不确定的问题.自适应陷波处理后
,得出的结果里有些幅值挺高.我想问一下,这样是否正常?
数据为肌电数据,采样频率5000Hz