LMS均衡.rar_LMS-4QAM_nobodycpp_信号均衡_水声信道_水声信道均衡

以下是基于LMS的程序 clear all; t=0:0.0001:1.9999; N=20000; %生成有用s s=3*sin(2*pi*100*t); %(1); %sub(3,1,1); %plot(t,real(s)); %title('有用信号s'); %xlabel('n'); %ylabel('s'); %生成两个相关的噪声n1,n2 n1=wgn(1,N,3); n2=zeros(1,N); n2(1:18000)=n1(1:18000); n3=wgn(1,N,3); n2(18001:N)=n3(18001:N); %figure(1); %subplot(3,1,2); %plot(t,real(n2)); %title('噪声参考'); %xlabel('n'); %ylabel('n2'); %输入端期望信号d=s+n1 d=s+n1; subplot(3,1,3); plot(t,real(d)); title('含有噪声的信号d'); xlabel('n'); ylabel('d'); %初始化各参量 M=8; %滤波器阶数 w=zeros(M,1); %滤波器权系数 mu=0.01; %更新步长 iter=2000; %迭代次数 x=zeros(M,1); %输入矢量 %LMS算法 for i=1:iter i x(1)=n2(i); y(i)=dot(w,x); e=d(i)-y(i); w=w+mu*e*x; for j=M:-1:2 x(j)=x(j-1); end (i)=e; %误差 cost(i)=e*e; end figure(2); subplot(1,1,1); plot(error); title('除噪后的信号'); xlabel('n'); ylabel('error');