mylms.m

自己编写的自适应notch fileter算法

t =0:1/1500:2;
s0 = 0.1*cos(0.5+cos(500*pi*t));
t =2:1/1500:4;
s1 = 0.1*cos(0+cos(500*pi*t));
q=length(t);n=1;
while n<=q
    s(n)=s0(n);
    n=n+1;
end
while n<=q*2
    s(n)=s1(n-q);
    n=n+1;
end
n=1:q*2;
subplot(231)
plot(n,s);
axis([2,6000,-0.3,0.3]);
title('噪声源信号曲线图');