Encode--LDPC校验矩阵H的高斯变换
这两天看LDPC得代码,刚开头,看到H矩阵的高斯变换部分,即得到H=[I P],输出P。
按自己的想法写了段matlab代码(注释那句matlab里应该是%):
for k = 1:rows if ( H(k,k) == 0 ) ind = find( H(k+1:rows,k)~=0 ); if ~isempty(ind) x = k+ind(unidrnd(length(ind))); H(k,:) = rem( H(k,:) + H(x,:), 2 ); end end for i = [k+1:rows] if ( H(i,k) == 1) H(i,:) = rem( H(i,:)+H(k,:), 2 ); end end end //查找有无不符合的行需要进行列变换的 for k = 1:rows if H(k,k) == 0 x = rows; cnd = find( H(k,rows+1:cols) ~= 0 ); if ~isempty(cnd) temp = H(:,k); H(:,k) = H(:,x+cnd(1)); H(:,x+cnd(1)) = temp; rearranged_cols(k) = x + cnd(1); end end end
上面实现了上三角阵,继续下面的代码就能完成要求的操作:
//将上三角矩阵变为前半部分单位矩阵 for k = 1:rows x = k; ind = find( H(k,k+1:rows) ~= 0 ); while ( ~isempty(ind) ) x = ind(1)+x; H(k,:) = rem( H(k,:) + H(x,:), 2 ); ind = find( H(k,x+1:rows) ~= 0 ); end end
虽然比较满意能够编出来,因为前几天看了但没理解,但还是网上给的方法编出来更简单些,而且用的循环少些。不管怎么说,是一个小节点的完成吧~~继续继续!!