lms_eq_changemu.m

自己编的基于最抖下降法

data_len=500;
repeat=200;
lamta=1;
var_v=0.001;
delta=0.004;
rou=2.9;
p=11;
tao=7;
mu=[0.025 0.075 0.0075];

for dif=1:3

E=zeros(repeat,488);
for repeat_n=1:repeat

data=2*binornd(1,0.5*ones(1,data_len))-1;
l=1;
for n=1:3
    h(n)=0.5*(1+cos(2*pi*(n-2)/rou));
end
x=conv(data,h);
xi=x((length(h)+1):length(x)); %x[4]~x[502]
xn=xi+sqrt(var_v)*randn(length(xi),1)';
d=data((tao+1):data_len);
w=(zeros(1,p))';




for n=1:(length(xn)-p)
    x1=xn(n:n+p-1)';
    y(n)=w'*x1;
    e(n)=d(n)-y(n);
    w=w+mu(dif)*e(n)*x1;
end
E(repeat_n,:)=e;
end


Q(1:488,dif)=mean(E.^2)';


semilogy(Q);
set(gca,'xlim',[1 500]);
title('LMS算法学习曲线');
xlabel('迭代次数n');
ylabel('集平均平方误差');
end