您现在的位置:首页 >> ★免费资源 >> 源码下载 >> 内容

MATLAB代做|FPGA代做|Python代做--LDPC译码算法的matlab程序

时间:2019-1-11 21:40:13 点击:

  核心提示:MATLAB代做|FPGA代做|Python代做--LDPC译码算法的matlab程序...
function x_hat = func_Dec(y,sigma,H,max_iter)

f1    = 1./(1+exp(-2*y/(sigma^2)));         
f0    = 1-f1;

[m,n] = size(H); 
if m>n
   H     = H'; 
   [m,n] = size(H); 
end

%如果不是稀疏矩阵,将H变为稀疏矩阵
if ~issparse(H) 
   [ii,jj,sH] = find(H);
   H = sparse(ii,jj,sH,m,n);
end

%初始化
[ii,jj] = find(H);        
%将非零元素标上
indx    = sub2ind(size(H),ii,jj); 
q0      = H * spdiags(f0(:),0,n,n);
sq0     = full(q0(indx)); 
sff0    = sq0;

q1 = H * spdiags(f1(:),0,n,n); 
sq1 = full(q1(indx));
sff1 = sq1;

k=0;
success = 0;

while ((success == 0) & (k < max_iter)),
   k = k+1;
   %若先验概率一样的话,令sdq为一个很小很小的数
   sdq = sq0 - sq1; 
   sdq(find(sdq==0)) = 1e-20; 
   %用概率差sdq构成m*n的稀疏矩阵
   dq = sparse(ii,jj,sdq,m,n);
   %将dq在对数域求和,即将dq相乘,构成数组
   Pdq_v = full(real(exp(sum(spfun('log',dq),2)))); 
   %更新了rij                                                                                                 
   Pdq = spdiags(Pdq_v(:),0,m,m) * H;
   sPdq = full(Pdq(indx));
   
   sr0 = (1+sPdq./sdq)./2; 
   sr0(find(abs(sr0) < 1e-20)) = 1e-20;%将每个元素的rij都归一化,更新rij。概率译码的第一步完成
   sr1 = (1-sPdq./sdq)./2; 
   sr1(find(abs(sr1) < 1e-20)) = 1e-20;%若比值非常小,则设为1e-20
   %将更新完的rij存入对应的位置,方便后面要用
   r0 = sparse(ii,jj,sr0,m,n);
   r1 = sparse(ii,jj,sr1,m,n);
   
   Pr0_v = full(real(exp(sum(spfun('log',r0),1))));
   Pr0 = H * spdiags(Pr0_v(:),0,n,n);
   sPr0 = full(Pr0(indx));
   Q0 = full(sum(sparse(ii,jj,sPr0.*sff0,m,n),1))';
   sq0 = sPr0.*sff0./sr0;
   
   Pr1_v = full(real(exp(sum(spfun('log',r1),1))));
   Pr1 = H * spdiags(Pr1_v(:),0,n,n);
   sPr1 = full(Pr1(indx)); 
   Q1 = full(sum(sparse(ii,jj,sPr1.*sff1,m,n),1))';
   sq1 = sPr1.*sff1./sr1;
   
   sqq = sq0+sq1;
   sq0 = sq0./sqq;
   sq1 = sq1./sqq;
   
   %伪后验概率
   QQ = Q0+Q1;
   Q0 = Q0./QQ;
   Q1 = Q1./QQ;
   
   %判决输出
    x_hat = (sign(Q1-Q0)+1)/2;
    if rem(H*x_hat,2) == 0, success = 1; end
end

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

Tags:LDPC译码 
作者:LDPC译码 来源:LDPC译码
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做|matlab代写(www.hslogic.com) © 2019 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment