码间干扰信道 均衡算法 ISI equalization

通过仿真ISI信道码间干扰信道并使用自适应均衡器来对信道进行均衡，提高ISI信道接收效果

clear;clc;echo off;close all;
N=10000;                             %指定信号序列长度
info=random_binary（N）;                  %产生二进制信号序列
SNR_in_dB=8:1:18;                     %AWGN信道信噪比
for j=1:length（SNR_in_dB）
    [ylen]=channel（infoSNR_in_dB（j））;   %通过既有码间干扰又有白噪声信道
    numoferr=0;                       %初始误码统计数
    for i=len+1:N+len                  %从第len个码元开始为真实信号码元
        if （y（i）<0）                    %判决译码
            decis=-1;          
        else
            decis=1;          
        end;
        if （decis~=info（i-5））             %判断是否误码，统计误码码元个数  
            numoferr=numoferr+1; 
        end;
    end;
    Pe（j）=numoferr/N;                   % 未经均衡器均衡，得到的误码率
end;
semilogy（SNR_in_dBPe‘red*-‘）;           %未经均衡器，误码率结果图
    hold on; 

delta_1=0.11;                           %指定自适应均衡器的步长
delta_2=0.09;                           %指定自适应均衡器的步长
for j=1:length（SNR_in_dB）
    y=channel（infoSNR_in_dB（j））;         %通过信道 
    z=lms_equalizer（yinfodelta_1）;        %通过自适应均衡器，并设置步长为0.11
    numoferr=0;
    for i=1:N
       if （z（i）<0）
          decis=-1;          
       else
          decis=1;          
       end;
       if （decis~=info（i））     
           numoferr=numoferr+1;
       end;
    end;
    Pe（j）=numoferr/N;                  % 经自适应均衡器均衡后，得到的误码率
end;
semilogy（SNR_in_dBPe‘blacko-‘）;        %自适应均衡器均衡之后，误码率结果图
    hold on;

    for j=1:length（SNR_in_dB）
    y=channel（infoSNR_in_dB（j））;        %通过信道 
    z=lms_equalizer（yinfodelta_2）;        %通过自适应均衡器，并设置步长为0.09
    numoferr=0;
    for i=1:N
       if （z（i）<0）
          decis=-1;          
       else
          decis=1;          
       end;
       if （decis~=info（i））     
           numoferr=numoferr+1;
       end;
    end;
    Pe（j）=numoferr/N;                   % 经自适应均衡器均衡后，得到的误码率
end;
semilogy（SNR_in_dBPe‘blue.-‘）;           %自适应均衡器均衡之后，误码率结果图
    hold on;

xlabel（‘SNR in dB‘）; 
ylabel（‘Pe‘）;
title（‘ISI信道自适应均衡系统仿真‘）;
legend（‘未经均衡器均衡‘‘经自适应均衡器均衡，步长detla=0.11‘‘经自适应均衡器均衡，步长detla=0.09‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2331  2012-09-24 13:43  main_plot.m