qc_ldpc矩阵构造

clc
clear all

%%%%%%基于自然序列构造1
A1=eye（73）;
B1=[182764125;2123680150;3164596175;42054112200];
for m=1:4
    for n=1:5;
H1（73*（m-1）+1:73*m73*（n-1）+1:73*n）=circshift（A1B1（mn）2）;
    end
end


%%%%%%%%基于自然序列构造2
A2=eye（127）;
B2=[11111;12345;1491625;182764125];
for m=1:4
    for n=1:5;
H2（127*（m-1）+1:127*m127*（n-1）+1:127*n）=circshift（A2B2（mn）2）;
    end
end

%%%%%%%%%%   EG_LDPC 无4环

C1=repmat（[1:6]21）;
C1 = reshape（C1 1 12）;
R1=[121413242334];
H = full（sparse（R1 C1 1 4 6））;
r=3;
a=eye（r）;   b=z