qcldpc是非验证邮箱的正则表达式ldpc的吗
点击联系发帖人
时间：2017-03-10 10:31
503 Service Temporarily Unavailable
503 Service Temporarily Unavailable
openresty/1.9.7.4下载费用：20 元（CNY） &
&&&&&&&&&&&&qc-ldpc码编码构造
毕业论文热点推荐
qc-ldpc码编码构造
下载文档到电脑，查找使用更方便
20 元（CNY）&&0人已下载
还剩页未读，继续阅读
资源关键词暂无 资源描述：中国科学技术大学硕士论文 摘要摘要低密度奇偶校验(Low．Density Parity．Check，LDPC)码是一种先进的信道编码技术，它具有逼近香农限的纠错性能。准循环低密度奇偶校验码(Quasi。Cyclic Low—Density Parity—Check，QC．LDPC)码是LDPC码的一个重要子类，它的校验矩阵具有准循环形式，这种结构特征决定了其较低的编解码复杂度。QC．LDPC码的主要构造方法包括有限几何法，均衡不完全区组 相关资源
暂无评论，赶快抢占沙发吧。rs、ldpc错码及其级联的研究.pdf -max上传文档投稿赚钱-文档C2C交易模式-100%分成比例文档分享网
rs、ldpc错码及其级联的研究.pdf
文档名称：rs、ldpc错码及其级联的研究.pdf
格式：pdf&&&大小：4.42MB&&&总页数：63
可免费阅读页数：63页
下载源文档需要：18元人民币
预览与实际下载的一致，文档内容不会超过预览的范围，下载前请务必先预览，自行甄别内容是否完整、是否存在文不对题等情况（本网站为文档分享平台性质），一旦付费下载，本站不支持退款
我已知晓：实际下载内容以预览为准！
文档介绍：摘要摘要本课题来源于与研究所合作项目，日的是实现某无线传输通信系统中的误码纠错。现代通信系统中设计的核心问题是在有随机噪声的信道中如何克服干扰，减小信息传输的差错。里德一索罗蒙RS码是一种具有很强的纠正突发和随机错误的信道编码方式，在深空通信、移动通信、磁盘阵列、数字视频广播D、偈以及军事等系统中都具有广泛的应用。低密度校验LDPC]i马是一类用非常稀疏的校验矩阵或二分图定义的线性分组纠错码，它的设计、构造、译码、快速编码、性能分析以及在数字通信和存储系统领域中的应用成为了当前纠错编码研究的焦点。本文首先对RS码的编解码原理作了详细的分析，并采用现场可编程门阵列gS31，l5码。然后分析和研究了准循环LDPC码的构造以及编解码算法，给出了LDPC编设计并实现了LDPC36,24的编码，并采用复杂度较低的硬判决译码方法进行译码。由于LDPC编解码的过程中涉及了大量的矩阵相乘运算，本文设计了一种循环移位求和电路，有效的降低了LDPC编解码中矩阵相乘运算的复杂度以及存储空间。并对其性能进行分析，该级联方法具有较好的纠错性能，同时又具有较低的复杂度的特点。通过实验表明，该级联码取得了比RS纠错码更好的纠错性能，能够符合该通信系统的要求。关键词：里德．索罗蒙码，低密度校验码，现场可编程门阵列，级联码AbstractTheissuestemsfromtheinstituteandtheaimistofulfillcooperationprojects、析ⅡlerrorinwirelesscoⅡm删cationmaininmodemcorr．；ringsystem．Theproblemcommunicationishowtodisturbancesandtoreducethetransmissionsystemconquererrorrandomnoiseexists．whenisone当前位置： &
& 查看源码
ldpc_decode_spa_m.m
ldpc_decode_spa_m.m：源码内容
function&[code]=ldpc_decode_SPA_M(r_wave,H,sigma_2,M)
%LDPC码的Log_BP(SPA)算法译码算法
%接口参数校验矩阵H,接收到信号波形r_wave
%sigma^2方差
%调制方式&M=2:BPSK，M=4:QPSK，M=16:16QAM,M=64:64QAM
[rows,cols]=size(H);
code(1:cols)=0;%定义硬判决输出
%(1)初始化模块
switch&(M)
&&&&case&2
&&&&&&&&Lq=-2*r_wave/sigma_2;%L(qi)
&&&&otherwise
&&&&&&&&[Lq]=Lq_initialization(r_wave,sigma_2,M);
Lq(cols+1:length(Lq))=[];
%&switch&(M)
%&&&&&case&2
%&&&&&&&&&%BPSK调制方式下的初始化
%&&&&&&&&&Lq=-2*r_wave/sigma_2;%L(qi)
%&&&&&&&&&Lq=-log(exp((r_wave-1).^2*-0.5/sigma_2)./exp((r_wave+1).^2*-0.5/sigma_2));
%&&&&&case&4
%&&&&&&&&&%QPSK调制方式下的初始化
%&&&&&&&&&Lq(1:2:cols)=-log((exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./(exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:2:cols)&&=-log((exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2))./(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&case&16
%&&&&&&&&&Lq(1:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(3:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(4:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&case&16&&&&&&&&&&&&&&&&&&&&&&&&
%&&&&&&&&&Lq(1:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(3:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&
L_q=repmat(Lq,rows,1);&%初始化L(qij)矩阵,为信息节点的值
[hi,hj]=find(H==1);
L_q=L_q.*H;
L_r=zeros(rows,cols);
for&iteration=1:50&%&迭代次数
&&&&iteration
&&&&&%&1)校验节点更新
&&&&for&i=1:rows
&&&&&&&&col_ind=find(hi==i);%找出与i个校验节点相连的信息节点集合
&&&&&&&&col_num=length(col_ind);
&&&&&&&&%temp=zeros(1,col_num);
&&&&&&&&temp=L_q(i,hj(col_ind));
&&&&&&&%&for&j=1:col_num
&&&&&&&&%&&&&temp(j)=L_q(i,hj(col_ind(j)));
&&&&&&&&a=sign(temp);&%&aij=sgn(L(qij))
&&&&&&&&a1=prod(a);&%求积
&&&&&&&&a2=a1*a;%rji
&&&&&&&&B=abs(temp);&%&Bij=abs(Lqij);
&&&&&&&&temp=tanh(B/2);
&&&&&&&&for&temp_i=1:length(temp)
&&&&&&&&&&&&if&temp(temp_i)&10e-324&
&&&&&&&&&&&&&&&&&temp(temp_i)=10e-324;
&&&&&&&&&&&&
&&&&&&&&Q=-log(temp);
&&&&&&&&Q1=sum(Q);
&&&&&&&&Q2=Q1-Q;%rji
&&&&&&&&temp=tanh(Q2/2);
&&&&&&&&for&temp_i=1:length(temp)
&&&&&&&&&&&&if&temp(temp_i)&10e-324&
&&&&&&&&&&&&&&&&&temp(temp_i)=10e-324;
&&&&&&&&&&&&
&&&&&&&&Q_2=-log(temp);
&&&&&&&&L_r_temp=a2.*Q_2;
&&&&&&&&for&j=1:col_num
&&&&&&&&&&&&L_r(i,hj(col_ind(j)))=L_r_temp(j);
&&&&%校验节点更新完毕
&&&&%&2)+3)
&&&&for&j=1:cols
&&&&&&&&row_ind=find(hj==j);
&&&&&&&&row_num=length(row_ind);
&&&&&&&&temp=zeros(1,row_num);
&&&&&&&&for&i=1:row_num
&&&&&&&&&&&&temp(i)=L_r(hi(row_ind(i)),j);
&&&&&&&&L_r_sum=sum(temp);
&&&&&&&&%&3)&尝试译码&
&&&&&&&&L_Q=Lq(j)+L_r_%软判决
%&&&&&&&&&code(j)=(sign(L_Q)+1)/2;%硬判决
&&&&&&&&if&L_Q&0
&&&&&&&&&&&&&&&&code(j)=1;
&&&&&&&&&else
&&&&&&&&&&&&&&&&&code(j)=0;&
&&&&&&&&end
&&&&&&&&%&2)&比特节点(信息节点)更新
&&&&&&&&L_r_temp=L_Q-
&&&&&&&&for&i=1:row_num
&&&&&&&&L_q(hi(row_ind(i)),j)=L_r_temp(i);
&&&&%比特节点更新完毕
&&&&s=mod(code*H',2);%检查是否满足中止条件
&&&&if&s==0
CopyRight & 2008- All Rights reserved. 苏ICP备
号 京公网安备：95}
503 Service Temporarily Unavailable
openresty/1.9.7.4下载费用：20 元（CNY） &
&&&&&&&&&&&&qc-ldpc码编码构造
毕业论文热点推荐
qc-ldpc码编码构造
下载文档到电脑，查找使用更方便
20 元（CNY）&&0人已下载
还剩页未读，继续阅读
资源关键词暂无 资源描述：中国科学技术大学硕士论文 摘要摘要低密度奇偶校验(Low．Density Parity．Check，LDPC)码是一种先进的信道编码技术，它具有逼近香农限的纠错性能。准循环低密度奇偶校验码(Quasi。Cyclic Low—Density Parity—Check，QC．LDPC)码是LDPC码的一个重要子类，它的校验矩阵具有准循环形式，这种结构特征决定了其较低的编解码复杂度。QC．LDPC码的主要构造方法包括有限几何法，均衡不完全区组 相关资源
暂无评论，赶快抢占沙发吧。rs、ldpc错码及其级联的研究.pdf -max上传文档投稿赚钱-文档C2C交易模式-100%分成比例文档分享网
rs、ldpc错码及其级联的研究.pdf
文档名称：rs、ldpc错码及其级联的研究.pdf
格式：pdf&&&大小：4.42MB&&&总页数：63
可免费阅读页数：63页
下载源文档需要：18元人民币
预览与实际下载的一致，文档内容不会超过预览的范围，下载前请务必先预览，自行甄别内容是否完整、是否存在文不对题等情况（本网站为文档分享平台性质），一旦付费下载，本站不支持退款
我已知晓：实际下载内容以预览为准！
文档介绍：摘要摘要本课题来源于与研究所合作项目，日的是实现某无线传输通信系统中的误码纠错。现代通信系统中设计的核心问题是在有随机噪声的信道中如何克服干扰，减小信息传输的差错。里德一索罗蒙RS码是一种具有很强的纠正突发和随机错误的信道编码方式，在深空通信、移动通信、磁盘阵列、数字视频广播D、偈以及军事等系统中都具有广泛的应用。低密度校验LDPC]i马是一类用非常稀疏的校验矩阵或二分图定义的线性分组纠错码，它的设计、构造、译码、快速编码、性能分析以及在数字通信和存储系统领域中的应用成为了当前纠错编码研究的焦点。本文首先对RS码的编解码原理作了详细的分析，并采用现场可编程门阵列gS31，l5码。然后分析和研究了准循环LDPC码的构造以及编解码算法，给出了LDPC编设计并实现了LDPC36,24的编码，并采用复杂度较低的硬判决译码方法进行译码。由于LDPC编解码的过程中涉及了大量的矩阵相乘运算，本文设计了一种循环移位求和电路，有效的降低了LDPC编解码中矩阵相乘运算的复杂度以及存储空间。并对其性能进行分析，该级联方法具有较好的纠错性能，同时又具有较低的复杂度的特点。通过实验表明，该级联码取得了比RS纠错码更好的纠错性能，能够符合该通信系统的要求。关键词：里德．索罗蒙码，低密度校验码，现场可编程门阵列，级联码AbstractTheissuestemsfromtheinstituteandtheaimistofulfillcooperationprojects、析ⅡlerrorinwirelesscoⅡm删cationmaininmodemcorr．；ringsystem．Theproblemcommunicationishowtodisturbancesandtoreducethetransmissionsystemconquererrorrandomnoiseexists．whenisone当前位置： &
& 查看源码
ldpc_decode_spa_m.m
ldpc_decode_spa_m.m：源码内容
function&[code]=ldpc_decode_SPA_M(r_wave,H,sigma_2,M)
%LDPC码的Log_BP(SPA)算法译码算法
%接口参数校验矩阵H,接收到信号波形r_wave
%sigma^2方差
%调制方式&M=2:BPSK，M=4:QPSK，M=16:16QAM,M=64:64QAM
[rows,cols]=size(H);
code(1:cols)=0;%定义硬判决输出
%(1)初始化模块
switch&(M)
&&&&case&2
&&&&&&&&Lq=-2*r_wave/sigma_2;%L(qi)
&&&&otherwise
&&&&&&&&[Lq]=Lq_initialization(r_wave,sigma_2,M);
Lq(cols+1:length(Lq))=[];
%&switch&(M)
%&&&&&case&2
%&&&&&&&&&%BPSK调制方式下的初始化
%&&&&&&&&&Lq=-2*r_wave/sigma_2;%L(qi)
%&&&&&&&&&Lq=-log(exp((r_wave-1).^2*-0.5/sigma_2)./exp((r_wave+1).^2*-0.5/sigma_2));
%&&&&&case&4
%&&&&&&&&&%QPSK调制方式下的初始化
%&&&&&&&&&Lq(1:2:cols)=-log((exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./(exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:2:cols)&&=-log((exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2))./(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&case&16
%&&&&&&&&&Lq(1:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(3:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(4:4:cols)=log((exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)));
%&&&&&case&16&&&&&&&&&&&&&&&&&&&&&&&&
%&&&&&&&&&Lq(1:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(2:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&Lq(3:6:cols)=log((exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+7).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+5).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&(exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)-3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2))./...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)-3).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+7).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+5).^2).*-0.5./sigma_2)+...
%&&&&&&&&&&&&&&&&&&&&&&&&&&&exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+1).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+1).^2).*-0.5./sigma_2)+exp(((r_wave(1,:)+3).^2+(r_wave(2,:)+3).^2).*-0.5./sigma_2)));
%&&&&&&&&&
L_q=repmat(Lq,rows,1);&%初始化L(qij)矩阵,为信息节点的值
[hi,hj]=find(H==1);
L_q=L_q.*H;
L_r=zeros(rows,cols);
for&iteration=1:50&%&迭代次数
&&&&iteration
&&&&&%&1)校验节点更新
&&&&for&i=1:rows
&&&&&&&&col_ind=find(hi==i);%找出与i个校验节点相连的信息节点集合
&&&&&&&&col_num=length(col_ind);
&&&&&&&&%temp=zeros(1,col_num);
&&&&&&&&temp=L_q(i,hj(col_ind));
&&&&&&&%&for&j=1:col_num
&&&&&&&&%&&&&temp(j)=L_q(i,hj(col_ind(j)));
&&&&&&&&a=sign(temp);&%&aij=sgn(L(qij))
&&&&&&&&a1=prod(a);&%求积
&&&&&&&&a2=a1*a;%rji
&&&&&&&&B=abs(temp);&%&Bij=abs(Lqij);
&&&&&&&&temp=tanh(B/2);
&&&&&&&&for&temp_i=1:length(temp)
&&&&&&&&&&&&if&temp(temp_i)&10e-324&
&&&&&&&&&&&&&&&&&temp(temp_i)=10e-324;
&&&&&&&&&&&&
&&&&&&&&Q=-log(temp);
&&&&&&&&Q1=sum(Q);
&&&&&&&&Q2=Q1-Q;%rji
&&&&&&&&temp=tanh(Q2/2);
&&&&&&&&for&temp_i=1:length(temp)
&&&&&&&&&&&&if&temp(temp_i)&10e-324&
&&&&&&&&&&&&&&&&&temp(temp_i)=10e-324;
&&&&&&&&&&&&
&&&&&&&&Q_2=-log(temp);
&&&&&&&&L_r_temp=a2.*Q_2;
&&&&&&&&for&j=1:col_num
&&&&&&&&&&&&L_r(i,hj(col_ind(j)))=L_r_temp(j);
&&&&%校验节点更新完毕
&&&&%&2)+3)
&&&&for&j=1:cols
&&&&&&&&row_ind=find(hj==j);
&&&&&&&&row_num=length(row_ind);
&&&&&&&&temp=zeros(1,row_num);
&&&&&&&&for&i=1:row_num
&&&&&&&&&&&&temp(i)=L_r(hi(row_ind(i)),j);
&&&&&&&&L_r_sum=sum(temp);
&&&&&&&&%&3)&尝试译码&
&&&&&&&&L_Q=Lq(j)+L_r_%软判决
%&&&&&&&&&code(j)=(sign(L_Q)+1)/2;%硬判决
&&&&&&&&if&L_Q&0
&&&&&&&&&&&&&&&&code(j)=1;
&&&&&&&&&else
&&&&&&&&&&&&&&&&&code(j)=0;&
&&&&&&&&end
&&&&&&&&%&2)&比特节点(信息节点)更新
&&&&&&&&L_r_temp=L_Q-
&&&&&&&&for&i=1:row_num
&&&&&&&&L_q(hi(row_ind(i)),j)=L_r_temp(i);
&&&&%比特节点更新完毕
&&&&s=mod(code*H',2);%检查是否满足中止条件
&&&&if&s==0
CopyRight & 2008- All Rights reserved. 苏ICP备
号 京公网安备：95}
我要回帖
更多关于 ip地址的正则表达式 的文章
更多推荐
- ·侵犯肖像权律师函头像来啦，快来挑选你喜欢的😍
- ·怎么想生个男孩有什么方法儿子啊，大家可以推荐嘛？
- ·我公司还欠对方公司欠货款不还怎么办最有效未付,但是对方公司已经注销了.这个要怎么付款,怎么走账?直接打给法人个人卡能行不
- ·一般纳税人海产品税率餐饮企业购水产海鲜取得增值税专用发票可以抵税吗?
- ·单方交通事故过了7天报警了没报保险后6个小时报警会影响吗，交警不给出具事故过了7天报警了没报保险责任证明也不给出具事故过了7天报警了没报保险证明？
- ·曲江二首诗意简短
- ·怎么像画画儿呢？仿写句
- ·怎么给老人/老太太化妆
- ·绿茶绿豆薏米做面膜可以加白糖吗
- ·性骚扰15岁女儿小说
- ·脸部过敏一直不消了，怎么办
- ·手脚皮肤角质厚多怎么办？
- ·廊坊市东晖化妆新东方烹饪学校学费表多少？
- ·美发纳米喷雾机效果怎么样
- ·lol找黑色玫瑰大神带我双排？
- ·加里奥被削重做r是不是只能飞队友
- ·奥拉星超级真诚之心在哪得 怎么得
- ·龙之谷手游地狱犬首通奖励怎么获得 地狱犬地狱模式三星攻略
- ·unity里面怎么设置pbr材质
- ·天天塔防中贵族11可以培养什么塔
- ·qcldpc是非验证邮箱的正则表达式ldpc的吗
- ·西游之路免费辅助小号多开辅助有吗
- ·安卓7.0小米5玩王者荣耀掉帧是不是掉帧啊？有一点点不兼容的样子
- ·最早穿越火线大炮是不是88cf永久大炮
- ·王者荣耀孙尚香铭文6神装是什么装备
- ·卖个csol卖号吧号 50元 有买的吗
- ·复仇者魔化流加点点有问题吗
- ·前期也是相当牛b，后期打魔界裂缝boss怎么打还是换成道士或者战士.虽然前期升级和打初中
- ·部落冲突靠谱无敌版在哪里下了，靠助手里面都下不到。
- ·龙之谷赤兔马怎么样手游羊咩咩和赤兔马哪个好
- ·口袋妖怪复刻技能解放假勇敢技能 怎么用
- ·使命召唤搞笑视频OL新出的这个动作是不是有点搞笑？
- ·天天酷跑仙灵萝拉属性战地犀牛属性
- ·拳风传奇和dnf黑暗禁咒镰刀哪个更值得跨界
- ·三国志13三国志9pk威力加强版版哪些武将有忠勇义烈的
