LDPC码Gallager构造校验矩阵(MATLAB)

M=9;%H的行数
N=12;%H的列数
j=3;%每列1的个数
k=4;%每行1的个数

H1=zeros(M/j,N);%初始化
for a=1:(M/j)
    H1(a,(a-1)*k+1:a*k) =1;
  
end
H2=[];%初始化
for c=1:j-1
b=randperm(N);
H2=[H2;H1(:,b)];
end
H=[H1;H2];

%消除4环
for m=1:M
for n=m+1:M
   w=and(H(m,:),H(n,:)) ;
d=find(w);
e=length(d);
if e>1
    if  length(find(H(m,:)))>length(find(H(n,:)))
        for f=1:e-1
           H(m,d(f)) =0;
          
        end
    else
        for f=1:e-1
           H(n,d(f)) =0;
        end
    end
end
    
    
end  
end  

    
disp(H)
 

  • 1
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值