7-10.m

matlab数字信号与应用源代码7

%例程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');