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