clear all;
close all;
N=8000;%运算循环的次数
x=zeros(3,N);
x(1,:)=randn(1,N);%产生N个正态分布高斯白噪声
x(2,:)=rand(1,N);%产生N个均匀分布的随机信号
x(3,:)=sin(0.01*pi*N);%产生正弦信号
d=zeros(1,N);%真实信号
e=zeros(1,N);%误差值
u=[0.1,0.15,0.2;0.1,0.15,0.2;0.1,0.15,0.2];%设定三种信号的步长
coefficient=[1,-0.5,0.3,-0.2,1,-0.8,0.6,-0.3];%真实系数
w=zeros(N,8);%系数估计值
for n0=1:3%for1开始
for n1=1:length(u(n0,:))%for2开始
    figure((n0-1)*3+n1);
for n=5:N%for3开始
    D=[fliplr(d(1,n-4:n-1)),fliplr(x(n0,n-3:n))];
    d(1,n)=D*coefficient.';
    e(1,n)=d(1,n)-w(n-1,:)*D.';
    w(n,:)=w(n-1,:)+u(n0,n1)*conj(e(1,n))*D;
end%for3结束
for n2=1:8%for4开始
    subplot(2,4,n2);
    plot(1:N,w(:,n2));
    xlabel('n');
    ylabel({['信号',num2str(n0),'系数估计值','(u=',num2str(u(n0,n1)),')']});
end%for4结束
end%for2结束
end%for1结束