lms1.m

接收端自适应均衡器的LMS算法MATLAB仿真实现

clc;
clear;
N=1000;
v=randn(1,N);
x(1)=v(1);d(1)=v(1);
d(2)=1.6*v(1)+v(2);y(1)=0;x(2)=d(2);
d(3)=1.6*x(2)-0.95*x(1)+v(3);x(3)=d(3);y(2)=x(1);
for n=4:N
     x(n)=1.6*x(n-1)-0.95*x(n-2)+v(n);
     y(n-1)=x(n-2);
     d(n)=x(n);
 end
 u=0.0003;
 w1(1)=0;
 w1(2)=0;
 w2(1)=0;
 w2(2)=0;
 for n=2:length(x)-1
     w1(n+1)=(1.0-2*u*y(n)*y(n))*w1(n)-2*u*y(n)*y(n-1)*w2(n)+2*u*y(n)*d(n);
     w2(n+1)=(1.0-2*u*y(n-1)*y(n-1))*w2(n)-2*u*y(n)*y(n-1)*w1(n)+2*u*y(n-1)*d(n);
 end
 subplot(2,1,1)
 plot(w1);
 ylabel('w1');
 title('w1(n)随n变化的曲线');
 subplot(2,1,2)
 plot(w2);
 ylabel('w2');
 title('w2(n)随n变化的曲线');