chase.m

TPC的Chase迭代译码时测试图样的程序

function D=chase(p,r_s,r_h,H)
 q=2^p;
% %---------------------确定y的最不可靠位元-----------------------------------
% %--------------------构造测试图样------------------------------------------
% %-------------找到3个最小数的标号---------------
[E,co]=test_error(p,r_s);
%********可在次将软输入转变为硬输入，减少输入项******
%------------------构造测试图样--------------------------------------------------
T=test_graph(E,co); 
%-----------------------构造测试序列--------------------------------------

for j=1:q
        Z(j,:)=mod((r_h+T(j,:)),2);
    end
% %-----------------------对Z代数译码--------------------------------------
    for j=1:q
        e(j,:)=yi_ma_biao(H,Z(j,:));
        CC(j,:)=mod(e(j,:)+Z(j,:),2);
    end
% %----------------------chase译码的硬输出码字D--------------------------------
n2=length(r_h);
    distance=zeros(1,q);
    for j=1:q
        for k=1:n2
         distance(j)=distance(j)+abs(CC(j,k)-r_s(k));
        end
    end
    for j=1:q-1
        if distance(j)<=distance(j+1)
            temp=j;
        else temp=j+1;
        end
    end
    D=CC(temp,:);