迫零均衡代码

很好用的zf算法 经过我的改正。希望对大家有帮助。

x1=randsrc（1101）;%生成长度为101，1，-1等概出现的随机输入序列
x2=[0.1 -0.25 1 -0.25 0.1];  %信道的系统函数
x=conv（x1x2）; %输入信号通过信道的ISI波形
subplot（211）
stem（x1‘bo‘）
hold on
stem（x‘gv‘）
axis（[1101-1.51.5]）
grid on
hold on                    
title（‘信道接收信号波形‘）;

X=zeros（33）;
for m=1:3
    for n=1:3
        X（mn）=x2（（length（x2）+1）/2+（m-1）-（n-1））;
    end
end
H=zeros（13）‘;
H（2）=1;
C=（X\H）‘

Y=conv（xC）;
subplot（212）
stem（x1‘bo‘）
hold on
stem（Y‘gv‘）
axis（[1101-1.51.5]）
grid on
hold on                    
title（‘均衡后接收信号波形‘）;

eyediagram（x2）;         %均衡前眼图
title（‘均衡前眼图‘）;
hold on;
eyediagram（Y2）;         %均衡后眼图
title（‘均衡后眼图‘）;
hold on;

for i=1:length（Y）
    if Y（i）>=0
        Z（i）=1;
    else Z（i）=-1;
    end
end

e=0;
for i=1:length（X）
    if Z（i）~=X（i）
        e=e+1;
    end
end

Pe=e./length（x）

        
        
    
    
    

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

     文件      22713  2009-11-24 10:49  迫零均衡\signal.png

     文件      20935  2009-11-24 10:50  迫零均衡\signal1eyediagram.png

     文件      15770  2009-11-24 10:50  迫零均衡\signal2eyediagram.png

     文件        962  2009-11-25 23:46  迫零均衡\ZeroForcingEqualization.m

     文件      99874  2009-11-25 23:49  迫零均衡\迫零均衡.docx

     目录          0  2010-09-08 15:19  迫零均衡

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

               160254                    6