%例程7－10  基于LMS算法的自適應均衡器
%仿真參數

isi=[0.18,0.3,1,0.18];       %ISI信道參數
order=63;                %濾波器階數
snr=10;                  %AWGN信道信噪比
len=1000;                %訓練序列長度
u=0.02;                  %調整步長

M=(order-1)/2;
N=len+length(isi)-1;
e=zeros(1,N);
error=e;
out=zeros(1,N);            %FIR濾波器輸出
for i=1:1000
    x=sign(rand(1,len)-0.5);         %采用PN碼作為訓練序列
    noise=randn(1,N)/10.^(snr/10);   %AWGN
    y=conv(isi,x)+noise;
    wk=zeros(1,order);
    for n=order:N-M+1            %LMS算法
        y1=y(n+M-1:-1:n-M-1);
        d1=wk*y1';
        e(n)=x(n-2)-d1;
        wk=wk+u*e(n)*y1;
        e(n)=10*log10(abs(e(n)));
    end
    error=error+e;
end
error=error(order:N-M+1)/1000;     %誤差值

t=1:length(error);
plot(t,error);
title('誤差收斂曲線');
ylabel('誤差dB');
xlabel('n');