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