关于均衡的MATLAB程序合集

%  完成一个自适应均衡器，要求均衡器的阶数为64，信源信号经过一个两径信道，                       
%  在接收端通过自适应均衡器得到期望信号，并比较有无均衡器情况下的误比特率。   
%  文件名为：adaequ.m

len = 20000;
Tlen = 2000;     %  假设用前2000个数据作为均衡器的训练数据
step = 0.001;
N = 64;          %  均衡器的抽头个数
s = zeros（1len）;
s1 = s;
x = zeros（1N）;
w = zeros（1N）;
s = randsrc（1len）;      %  信源
s1（2:len） = s（1:len-1）;  %  第2径
p = 0.9;
SNR = 0:10;
for db = 1:length（SNR）
     s2 = sqrt（p）*s + sqrt（1-p）*s1;
     s3 = awgn（s2db‘measured‘）;
     for i = N:len
           u（1:N） = s3（i:-1:i-N+1）;
           y（i） = u*conj（w.‘）;
           e（i） = u*w‘-conj（s（i））;
           w = w - step*u*（e（i））;
           if  y（i）>0                %  有均衡器的判决输出
                y1（i） = 1;
           else
                y1（i） = -1;
           end
           if  s3（i）>0              %  无均衡器的判决输出
                y2（i） = 1;
           else
               y2（i） = -1;
           end
     end
     errornum1 = sum（y1（Tlen:end） ~= s（Tlen:end））;
     errornum2 = sum（y2（Tlen:end） ~= s（Tlen:end））;
     ber1（db） = errornum1/（len-Tlen）;
     ber2（db） = errornum2/（len-Tlen）;
end
semilogy（SNRber1‘+-‘）;hold on;
semilogy（SNRber2）;
xlabel（‘信噪比‘）; ylabel（‘误比特率‘）;
legend（‘有均衡器‘‘无均衡器‘）;
title（‘均衡器的误比特率曲线‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    4815360  2007-04-25 09:35  zishiyinglvbodebiyesheji.doc

----------- ---------  ---------- -----  ----

              4815578                    2