testall.m

来自「自己设计的LT数字喷泉码（信道编码）的解码和测试程序。」· M 代码 · 共 56 行

clear all;
k=1000;ml=1;max=100;indx=1;
rr=[1:0.1:2];
for rate=rr
    q=0;
    for cnt=1:max
        for j=1:k
            A(j,:)=randint(1,ml)';
        end;
        % LT coding
        [X,n,g,c_deg,v_deg,c_node,v_node] = LT_encode(k,ml,A,rate);
        
        % LT decoding
        [B,ret] = LT_decode(k,n,ml,X,g,c_deg,v_deg,c_node,v_node);
        
        check = sum(sum(abs(A-B)));
        fprintf('LT [%d,%d]译码结果%d,ret=%d\n',k,n,check,ret);

        if check ~=0 || ret ~= 0
            q=q+1;
        end;
    end;
    q = q/max;
    eq(indx)=q;
    indx=indx+1;
end;
figure(1);
plot(rr,eq,'b');hold on;grid on;

indx=1;
for rate=rr
    q=0;
    for cnt=1:max
        for j=1:k
            A(j,:)=randint(1,ml)';
        end;
        % LT coding
        [X,n,g,c_deg,v_deg,c_node,v_node] = LT_encodeYang(k,ml,A,rate);
        
        % LT decoding
        [B,ret] = LT_decode(k,n,ml,X,g,c_deg,v_deg,c_node,v_node);
        
        check = sum(sum(abs(A-B)));
        fprintf('LT [%d,%d]译码结果%d,ret=%d\n',k,n,check,ret);

        if check ~=0 || ret ~= 0
            q=q+1;
        end;
    end;
    q = q/max;
    eq(indx)=q;
    indx=indx+1;
end;
figure(1);
plot(rr,eq,'r');hold on;grid on;