中衡LDPC码程序【设计明细】.docxVIP

LDPC码编/译码程序设计 要求： 用Matlab编程实现（2016,1008）LDPC码的编码器和译码器，并搭建仿真系统统计误码性能。 设计内容： 一、LDPC码参数 二、编码器 三、噪声信道 四、译码器 和积算法译码器 最小和算法译码器 修正最小和算法译码器五、仿真解析 一、LDPC码参数 本设计采用准循环LDPC码码型（QuasiCyclicLDPC，QC-LDPC），此种码型性能优异， 并有利于硬件实现，现已被各种工业标准所采纳。比方， Wimax标准，DVB-S2标准。详尽 矩阵请见：Matrix(2016,1008)Block56.mat 。拥有如下所示的结构： 本设计采用系统码，码长为 2016比特，1/2码率，信息序列长度为 1008比特。其H矩 阵表示如下： H1,1 H1,2 H1,nb H H2,1 H2,2 H2,nb Hmb,1 Hmb,2 Hmb,nb 其中Hi,j 是大小为z×z的循环移位矩阵，行重为1，它的值表示该矩阵的循环移位偏 移量。比方： Hi,j＝n，则表示此矩阵第 1行的第n列为1，其余列为 0，其余各行均是上 一行的循环移位；如果Hi,j0 ，则表示该矩阵是一个全零矩阵。比方 Hi,j＝4，z＝7，则 对应的矩阵结构为： 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 Hij 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 本设计给出的 QC-LDPC码，z取56。 由于本设计采用系统码， H矩阵被分为两个部分 Hp|Hs ： H1,1 H1,kb H1,kb 1 H1,nb H H2,1 H2,kb H2,kb1 H2,nb Hmb,1 Hmb,kb Hmb,kb1 Hmb,nb 其中Hp对应校验比特部分，大小为mbzmbz；Hs对应信息比特部分，大小为 mbzkbz，其中kb (nbmb)z。Hp拥有如下规则结构： 1 0 0 0 0 a 1 1 0 0 0 0 0 1 1 0 0 Hp 0 1 0 0 0 0 1 1 0 0 0 0 0 1 1 其中Hp右上角a表示的矩阵拥有如下结构： 000 0 01 00 0010 【注】 上文Hp中的子矩‘a’由文档直接出，不利用循移位方法算； 本次作中的LDPC系，其中校位在前，信息在后。 二、LDPC码编码器 LDPC属于性分，可以用性分的方式行， 即利用生成 矩行。但LDPC平时用其稀疏校矩（H矩）行定，把校矩化成 生成矩运算量极大。由于 LDPC一般，校矩大，就更增加了运算担。 并且化的生成矩并不是稀疏的， 就硬件来了困。中一般采用校矩 直接行，利用公式 VHT 0（其中，V生成的字），直接解出字V。里 在不影响性能的前提下， 一般需要特殊构的校矩， 以化解方程运算。 详尽算法如 下： 由于采用的是系，利用入信息比特s s1,s2,,skbz，获得校比特 p1,p2,,pmbz即可完成。分如下两个步行： 步一：利用入信息矢量s，算中果xx1,x2,,xmbz： xsHs 步二：利用x算校比特p，采用如下公式：  T xi ,i 1 p x p 1)zi1 ,1 iz，其中 表示GF(2)中的加法。 i i (m b xi piz ,i z 首先获得第一个校比特p1 x1，然后依次获得pz1 xz1p1， p2z 1x2z1 pz 1，??，p(mb1)z1 x(mb1)z1 p(mb2)z1，p2x2 p(m 1)z1，??， b pmbz xmbz p(mb1)z。 三、噪声信道 本需要在AWGN信道下行任性能仿真，制方式采用BPSK制（0映射 成1，1映射成-1），可能用到的Matlab函数有：wgn也许awgn。 四、LDPC码译码器 本设计要求编程实现三种LDPC码译码器，一种为在系统仿真中较为常用的和积算法 译码器，另一种为利于硬件实现的最小和算法以及其修正算法译码器。详尽算法流程如下： 和积算法程序流程 为了方便说明，定义LDPC码长为N，校验方程数目为M，则该LDPC码的二分图包 含了N个变量节点，M个校验节点以及条边线。在该二分图中，和积算法（sum-product algorithm）的译码过程被分为以下三个主要步骤，如图1所示： 初始化： 每个比特从信道中接收到的信道信息作为初始置信度，赋予相应的变量节点，它们 的置信度uji u0 i，其中u0i 表示第i个比特从信道接收到的置信度 4yi ，其 SNR 中，yi为信道输入译码器的软信息， SNR Es 为符号能量与噪声的比值（注意： N0 (SNR)dB 这里单位不是 dB，而是单纯的比值，即10 10 ）。 2.迭代： a) 每个变量节点将自己的置信度传