getg.m

Mackey1随机构造LDPC码的校验矩阵进行编译码的仿真过程（其中译码算法采用置信传播算法）

m=504;n=1008;
H=getH(m,n);
[G,valid]=H2G(H); 
while valid==0                   %valid作为校验矩阵是否为满秩的标志,若为非满秩(valid=0),则返回重新利用1A随机构造
H=getH(m,n);                     %构造校验矩阵H
[G,valid]=H2G(H);                %将校验矩阵H转化为生成矩阵G
end

t=mod(G*H,2);
[i.j]=find(t);

save G