GALLAGER校验矩阵详细构造代码

matlab代码实现gallager校验矩阵构造，过程详细，已验证可行，适合初学者，大家可以自己先试试编一下在对比一下，不难。

%参数设置
n = 20; % Number of columns
w_c = 3; % Column weight
w_r = 4; % Row weight
k = n*w_c/w_r; % Number of rows
H_sub = zeros（n/w_rn）; % First sub-matrix; there are w_c such sub-matrices.

%产生第一个子矩阵
for i = 1:n/w_r
for j = （i-1）*w_r+1:i*w_r
H_sub（ij） = H_sub（ij） + 1;
end
end

%对后两个矩阵进行列置换
H_pre = H_sub;
for t = 2:w_c
 x = randperm（n）;
 H_sub_perm = H_sub（:x）;
H_pre = [H_pre H_sub_perm];
end
%将生成的三个子矩阵上下排列好，结束
H = zeros（kn）;
for p = 1:w_c
 
H（（p-1）*（n/w_r）+1:（p）*（n/w_r）1:n） = H（（p-1）*（n/w_r）+1:（p）*（n/w_r）1:n） + H_pre（:（p-1）*n+1:p*n）;
end
H

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        626  2018-11-03 20:28  GALLAGER.m

----------- ---------  ---------- -----  ----

                  626                    1