LMS自适应信道均衡程序以及学习曲线绘制_matlab_LMS_信道均衡

%%LMS自适应均衡算法学习曲线，信噪比为20dB clc; clear all; M=1000; N=5; L=11; a=5*1e-3; d=6; H=[-0.557 -0.5 1.0 0.033 0.35]; ee=zeros(1,M-d); for nn=1:1:200; x=rand(1,M); x=x-0.5; for n=1:1:length(x) if x(n)>=0 x(n)=1; else x(n)=-1; end end y=filter(H,1,x); y=y+randn(1,M)/5; w=zeros(1,L); %w=(rand(1,L)-0.5)*2; E1=1; E2=0; P=0; d=6; %while abs(E2-E1)>1e-5 E2=E1; E1=0; for m=d+1:1:L+1 sum=0; for n=1:1:m-1 sum=sum+y(m-n)*w(n); end e=x(m-d)-sum; for n=1:1:m-1 w(n)=w(n)+2*a*e*y(m-n); end E1=E1+e*e; ee(m-d)=ee(m-d)+e.^2; end for m=L+2:1:M sum=0; for n=m-1:-1:m-L sum=sum+y(n)*w(m-n); end e=x(m-d)-sum; for n=m-1:-1:m-L w(m-n)=w(m-n)+2*a*e*y(n); end E1=E1+e*e; ee(m-d)=ee(m-d)+e.^2; end E1=E1/(M-d); clear x y; end %end ee=ee/200; M=10000; x=rand(1,M); x=x-0.5; for n=1:1:length(x) if x(n)>=0 x(n)=1; else x(n)=-1; end end y=filter(H,1,x); y=y+randn(1,M)/10; %下面进行测试 P=0; for m=d+1:1:L+1 sum=0; for n=1:1:m-1 sum=sum+y(m-n)*w(n); end if sum>=0 sum=1; else sum=-1; end if sum~=x(m-d) P=P+1; end end for m=L+2:1:M sum=0; for n=m-1:-1:m-L sum=sum+y(n)*w(m-n); end if sum>=0 sum=1; else sum=-1; end if sum~=x(m-d) P=P+1; end end str=sprintf('错误率是：'); disp(str); disp(100*(double(P)/(M-d))); hold on; plot(ee);