checkh.m

来自「包含五种LDPC码编码生成矩阵的编码算法及其比较」· M 代码 · 共 25 行

function [row_flag,col_flag,row_var,col_var]=checkH(H)
%[row_flag,col_flag,row_var,col_var]=checkH(H)

[rows cols]=size(H);

%计算行重
row_flag(1:rows)=0;
for i=1:rows
   ind=find(H(i,:)==1);
   row_flag(i)=length(ind);
end

%每行1的个数的方差
row_var=var(row_flag);

%计算列重
col_flag(1:cols)=0;
for j=1:cols
   ind=find(H(:,j)==1);
   col_flag(j)=length(ind);
end

%每列1的个数的方差
col_var=var(col_flag);