lms.rar_LMS MATLAB_matlab 陷波_notch_信号 陷波_信号陷波

clear;close all; N=40; NS=400; Fs=8000; u=0.002; v=rand(1,NS)-rand(1,NS); s=2*sin(2*pi*[0:(NS-1)]*100/Fs); figure,hold on,grid on x=s+v; %plot(x,'k-.') delayx=[0,x(1:NS-1)]; [w,y,e]=lms(x,delayx,u,N); plot(v,'g-') plot(s,'b-') plot(y,'r-') %legend('加噪信号','噪声信号','理想信号','滤波后信号'); legend('噪声信号','理想信号','滤波后信号'); %m=fft(w); [m,a]=freqz(w); figure %subplot(211),plot(a,20*log10(abs(m))),title('幅频响应'); %subplot(212),plot(a,20*log10(angle(m))),title('相频响应'); subplot(211),plot(a,abs(m)),title('幅频响应');grid on subplot(212),plot(a,angle(m)),title('相频响应');grid on