lms.rar_LMS信号抵消_LMS噪声_自适应 噪声 抵消_自适应法

clear all; clc f=[200;300]; A=[1,2]; fs=1000; t=(0:fs)/fs; s=A*sin(2*pi*f*t); %产生有用周期信号 u=wgn(1,length(t),8);%产生一个m行n列的高斯白噪声的矩阵，强度为8dBW。 m=length(t); d=s+u; d=zeros(1,m); delay=3; x1=zeros(1,m); for i=delay+1:m x1(i)=x(i-delay); end mu=0.02; order=4; w0=zeros(1,4); S=initlms(w0,mu); %初始化S信息结构体 [y,e,S] =adaptalms(x1,d,S); %自适应滤波 plot(t,s) title('周期信号2sin(2*pi*200*t)+sin(2*pi*300*t)时域图'); subplot(412); plot(t,d); title('高斯白噪声+有用信号的时域图'); subplot(413); plot(t,y); title('经过滤波的周期信号时域图'); subplot(414); plot(t,e); title('误差信号时域图'); figure; subplot(411); plot(abs(fft(s))); title('周期信号2sin(2*pi*200*t)+sin(2*pi*300*t)频域图'); subplot(412); plot(abs(fft(d))); title('高斯白噪声+有用信号的频域图'); subplot(413); plot(abs(fft(y))); title('经过滤波的周期信号频域图'); subplot(414); plot(abs(fft(e))); title('误差信号频域图');