毕业设计（论文）-基于MATLAB的LDPC码的仿真与研究

1.2LDPC发展动态 2LDPC码的原理 3LDPC码校验矩阵构造 3.2校验矩阵的结构化构造 4LDPC码的编码 4.2LU分解 4.4LDPC码的有效编码 5LDPC码的译码 5.2降低复杂度的BP译码 6运用MATLAB对LDPC码仿真 6.1不同码长对LDPC码的影响 6.2不同迭代次数对LDPC码的影响 6.3不同列重对LDPC码的影响 参考文献 致谢 淮南师范学院2012届本科毕业论文1基于MATLAB的LDPC码的仿真与研究学生：指导教师：淮南师范学院电气信息工程学院摘要:随着无线通信技术的发展和各种传输方式对可靠性要求的不断提高，信道编码作为抗干扰技术的重要手段之一，在数字通信技术领域和数字传输领域显示出越来越重要的作用。目前3G和4G无线通信系统的信道编码的主流技术是Turbo码，其具有接近香农极限的优异性能，但是在码长较长的情况下，其误码性能还有待于提高。LDPC码是一种线性分组码，和Turbo码同属于复合码类。两者的性能相近，且两者的译码延迟都很长，其适用一些实时性要求不是很高的通信。但是在码长较长的情况下，LDPC码比Turbo码的译码性能更高而且在误码率上，在硬件应用和应用方面LDPC具有更大的优点。本课题主要研究LDPC码在不同条件下的误码率，译码算法和应用上相比其他编解码技术的优点，为B3G和4G无线通信系统的信道编码提供理论依据。本论文所做工作如下:本文阐述了LDPC码的基本原理和分析LDPC码关键技术及影响性能的因素;对LDPC码进行编码在MATLAB软件的环境下进行仿真，在仿真的图像上进行性能的分析;运用仿真的图像对LDPC和Turbo在误码率上进行分析比较;在研究LDPC码基本理论的基础上,利用MATLAB仿真比较不同码长、列重和迭代次数对LDPC码性能的影响，为B3G和4G移动通信提供了有利的参考价值。关键词：LDPC码；MATLAB；仿真ResearchandSimulationofLDPCcodeBasedonMATLAB 基于MATLAB的LDPC码的仿真与研究2Keywords:LDPCcode；MATLAB；simulati 淮南师范学院2012届本科毕业论文321世纪是一个信息化的时代，信息的传输即通信起着支撑的作用。随着科学技术的发展，人们对通信的要求与日俱增，世界各国都在致力于现代通信技术的研究与开发。无线通信是现代通信系统的重要组成部分。经过专家和学者的努力，在过去几十年中，无线通信已经取得了很大程度上的进步。但是现代社会的发展是更加迅速的，这就对无线通信提出了进一步的要求，不仅要高速率的传输数据，而且要高质量的传输数据。为了满足这种高速率和高质量的传输需求，本文提出了一种新的编码调制技术。古时“千里眼，顺风耳”表达了人们对通信特别是远距离通信的向往。19世纪末，马可尼无线通信实验的电报声宣布了一个新时代的到来。二十世纪是无线通信发展的黄金时期，40年代无线通信理论得到了长远的发展。可捷尔尼可夫和香农为无线通信理论的发展做出了卓越的贡献。近20年，无线通信的发展经过了三代。1978年底，美国贝尔实验室研发了高级移动电话系统（AMPS，AdvancedMobilePhoneService建成了模拟蜂窝语音通信系统，诞生了第一代无线通信系统。同时，其他国家也研发出采用频分多址（FDMA，FrequencyDivisionMultipleAccess）方式的模拟蜂窝无线通信系统，包括以美国AMPS，英国TACS（TotalAccessCommunicationSystem北欧NMT450/900等。二十世纪八十年代，进入了数字通信的时代，欧洲率先研发出了全球移动通信系统（GSM，GlobalSystemforMobileCommunication日本和美国也研发了自己的标准。二十世纪九十年代，美国推出了窄带码分多址（CDMA，CodeDivisionMultipleAccess）数字蜂窝无线通信系统[1]。这也就是目前正普遍使用的无线通信系统，即第二代无线通信系统。第二代无线通信系统主要支持语音和低速率的数据业务。当今世界的第二代数字无线标准，包括GSM、D-AMPS、JDC（JapanDigitalCellular）和IS-95CDMA等，仍然是窄带系统。现有的无线通信网络主要以GSM和CDMA为主，采用GSM、GPRS、CDMA的IS-95B技术，速率可达115.2kbit/s，全球移动通信系统（GSM）采用增强型数据速率（EDGE）技术，速率可达384kbit/s。随着人们的物质和文化水平的提高，对无线通信业务的需求日益增多，这种现存的无线通信网难以满足新的需求。为适应新的需求，第三代（3G，3rdGeneration）无线 基于MATLAB的LDPC码的仿真与研究4TelecommunicationUnion）提出。目前，被世界公认的主流标准为：欧洲和日本提出的WCDMA、北美的CDMA2000和中国的TD-SCDMA方案。虽然3G在通信容量和质量上较2G有了较大的提高，但是其核心技术没有发生革命性的变化，所以3G可看作是2G向未来无限通信系统发展的一个过渡。随着时代的发展，人们已经提出了B3G乃至4G无线通信系统，该系统有着不可替代的优点，不仅进一步提高了系统的容量和数据的传输质量，而且可实现数据的高速率和高质量传输。B3G和4G系统涵盖了现有的3G和3G增强型技术以及新的移动接入和本地接入系统。在各种技术中，正交频分复用（OFDM）技术受到了极大的关注。B3G蜂窝移动通信系统需要具备以下能力：（1）支持全IP高速分组数据传输，数据速率可高达数百Mbit/s。（2）支持高传输质量，传输数据的误码率低于106。（3）提供高的频谱利用率和功率效率，发射功率降低10dB以上。（4）支持高终端移动性，移动速度高达几百km/h。（5）能够支持在用户数据速率、用户容量、服务质量和移动速度等方面大范围的动态变化。为满足这些技术上的需求，B3G移动通信系统必须在系统理论、网络结构和关键技术等方面具有突破性的改变，具体为：（1）在网络结构方面，将采用全IP、分布式、自组织和多层的无线广带个人通信新体制和新模式，以对抗2G以上电波传输特性的挑战，并适应未来移动通信以数据业（2）在传输体制方面，传统的单载波时分多址技术和码分多址技术很难直接推广到广带传输，必须采用多载波并行传输体制，而设计高度灵活的多载波传输方案是设计B3G移动通信系统的一个关键。（3）在编码与调制技术方面，将采用新型的自适应编码调制技术，而包括其在内的、高效的自适应链路技术则是B3G移动通信的另一个关键技术。（4）在空中接口方面，将采用分布式的接入方式，多天线环境下多输入多输出（MIMO）无线通信系统的理论将突破传统的无线通信系统理论，成为未来移动与无线通信系统理论的核心。（5）在天线与射频技术方面，将采用新型的多天线和阵列天线技术以及宽高线性 淮南师范学院2012届本科毕业论文5度射频技术。在20世纪60年代Gallager[3]在他的博士论文中提出了低密度分组校验码LDPC（Low经数十年的沉寂，随着计算机能力的增强和相关理论的发展，Mackay和Neal重新发现了它，并证明它在与基于BP（BeliefPropagation）的迭代译码相结合的条件下具有逼近Shannon极限的性能。LDPC的重新发现是继Turbo码后纠错编码领域又一重大进展。LDPC码的特点是：性能优于Turbo码，灵活性大；译码复杂度低于Turbo码，可以完全并行操作，硬件复杂度较低，所以硬件易于实现；描述简单，对严格的理论分析具有可验证性；吞吐量很大，能够高速译码。经过大量实验表明，AWGN信道下，码长足够长时码率为1/2的非正则LDPC码可达到距离香农限0.13dB[4]。LDPC码在数据可靠传输中的良好应用前景引起了学术界和IT业界的高度重视，成为现今信道编码领域最受瞩目的研究热点之一。下一代数字卫星视频广播标准DVB-S2已经采纳了基于LDPC码的编码方案，在第四代通信系统LDPC码将得到广泛的应用。本论文做了如下工作：对数字通信与差错控制编码的基本理论进行了概述。介绍了LDPC码的历史与发展，并对LDPC码的编译码算法的研究现状进行阐述了其在各个方面的应用，介绍了LDPC码的编解码的设计和译码的算法。重点运用matlab软件对其LDPC码码长，列重和迭代次数进行仿真：（1）应用matlab软件针对规则LDPC码码长分别为300、500和1000,列重选择。（2）应用matlab软件针对码长为500，译码的最大迭代次数为20，校验矩阵每列的列重分别为2、3和4情况下进行了仿真实验。（3）应用matlab软件针对码长500的规则LDPC码,列重为2,译码迭代次数分别为10、20和40的情况下进行了仿真实验。(4)对LDPC码和Turbo码在误码的性能上进行仿真并且进行系统的分析2LDPC码的原理LDPC码是一种非常特殊的线性分组纠错码[5]。通过生成矩阵G可以将线性分组纠错码要发送的信息s转换成被传输的码字t，与G对应的是校验矩阵H，满足H×t=0。LDPC 基于MATLAB的LDPC码的仿真与研究6码校验矩阵中非0元素的个数远远小于0元素的个数，是稀疏矩阵。LDPC码可以分为规则（Regular）码和非规则（Irregular）码两种，二者的主要区别在于它们的校验矩阵的行和列权重是否变化。假设LDPC码的校验矩阵H是M×N维，而且满秩，则LDPC码长为N，校验位为M，信息位则为k=M-N，码率r=k/N。H矩阵中每行“1”的个数为行权重，每列中“1”的个数为列权重。H矩阵可以用二部图Node）或比特节点（BitNode上边的M个节点代表M个校验式，称为校验节点（CheckNode）。当矩阵中元素Hmn1时，信息节点n和校验节点m就可以通过对应的边（Edge）连接起来，连接两个节点的边称为这两个节点的相邻边，相互连接的两个节点称为相邻节点，和每个节点连接的边的数量称为该节点的度（Degree）[6]。对于规则LDPC码，校验矩阵H中行权重和列权重都是固定不变的，在二部图上，信息节点和校验节点的度都分别为固定值，我们用（N，j，k）来表示，其中N为码长，j为校验矩阵的列权重，k为行权重。对于非规则LDPC码，其Tanner图中上、下任一个节点的度数都不是固定值，分别占上、下总度数的一定比例i、i。规则码是非规则码的一个特例。3LDPC码校验矩阵构造在介绍LDPC码校验矩阵的构造之前，首先阐述一下什么是girth。图2中，粗线部 淮南师范学院2012届本科毕业论文7图2二部图中的girthGallager基于GF（2）域上定义的（n，j，k）LDPC码[7]，其校验矩阵H的构造如（1）将Gallager码的监督矩阵按行划分成j个部分（每部分包含相同的行数每一部分的每一列中只包含一个“1”。（2）第一部分构造的矩阵中，“1”比特在行中按降幂排列，在第一行中，第1到（3）其余j-1部分的构造是对第一部分进行列的随机重排。法构造的（20，3，4）的LDPC码校验矩阵，码长为20，j=3，k=4。110000000000000000001111000000000000000000111100000000000000000011110000000000000000001111001000100010000000000100010000001000100010000001000100010000001000100010000001000100010001000100000100000100000010001000010000100001000010000010010000100001001000001000010000100001图3Gallager构造的（20，3，4）的LDPC码校验矩阵 基于MATLAB的LDPC码的仿真与研究82Mackay构造法Mackay提出了3种校验矩阵的构造方法[8]，使其二部图中循环的数目最少，得到的矩阵去掉了长度为4的短环，分别如下：（1）构造1A这是一种最基本的构造方法，保证矩阵列重t，通常取t=3，保证行重在每行中均匀分布，而且任意列间不存在长度为4的短环。如图4所示。图41A构造（2）构造2A矩阵中有M/2列列重为2，用2个M/2×M/2的单位矩阵上下摆放，其余的（N-M/2）列按照1A法构造，任意列间的重叠不超过1。如图5所示。3图52A构造（3）构造1B、2B从1A、2A中特意删去二部图中一些出现短环的列，再插入随Davey构造法[9]把每个节点按照它们的度d变成d个节点，然后构造二部图。设H中所有非0元素的个数为T，则有Tniimii，下半部和上半部分别有T个节点（D1,D2Ti的列用i个下节点表示，用其列号表示这些节点，同样把行重为j的行用j个上节点表示。接下来对上节点进行随机排列构造校验式，为确保不出现双边，即一个上节点的两条边参与同一个校验式，在排列时要保证同名的上、下节点不相互对应。4.比特填充及扩展的比特填充方法 淮南师范学院2012届本科毕业论文9m，其中g是偶数，构造一个m×n维校验矩阵H，每列有a个1，每行最多有b个1，girth为g，码率尽量大；二是给定变量节点数n、度为a，校验节点数m、度为b，满足n×a=m×b，使得g较大；三是给定n、a、g，使得m较小，也就是使码率较大。比特填充算法的实现：假设我们已得到一个校验矩阵H，n列、列重a、行重未超过b、girth大于等于g，现在在H上添加第n+1列：用一个初值为空的集合U1表示第n+1列，该集合长度为a，元素为各校验节点，是{l，2，…，m}的子集之一，再假设已添加j个元素到U1（0<j<a按条件再添加第j+1个元素到U1，这样循环往复下去就会得到整个校验矩阵H。过减少girth到g-2这种方式，在列不能增加时，保证列可以继续增加下去。5.PEG方法PEG（ProgressiveEdge-Growth）方法是Xiao-YuHu等提出的，是一种在二部图上以增加girth长度为目的的方法[12]。具体操作如下：对于给定的信息节点数n、校验节没有较大影响，然后选择程序放置，接着继续搜索下一边放置，直至结束。具体算法如下：j)少E，E是比特节点bi的第1条入射边，cj是校验节点在当前图集合Eb0Eb1…Ebi1中度数最低的。 Else，在当前图中从比特节点bi开始采用树形结构扩大到深度为l，直到Ni式，j)少E，E是比特节点bi的第k条入射边，cj是集合Ni中度数最低的校验节点。1有限几何构造法 基于MATLAB的LDPC码的仿真与研究10有限几何构造法是基于有限几何中的线和点来进行的，包括欧式几何（EuclideanGeometry，EG）和投影几何（ProjectiveGeometry，PG）。EG和PG具有n个点和J条线，满足以下条件[13]：每条线有p个点；任意两点间有且只有一条线；每个点只能落在q条线上；两条线平行或者只有一个交点。与有限几何G对应的GF（2）上的J×n维矩阵为H=（hij矩阵的行和列分别对应G的线和点。若第i条线包含第j个点，则hij=1，否则hij=0。矩阵的某一行表示这条线包含的所有点，重量为q；某一列表示穿过这个点的所有线，重量为q。矩阵的行称为G中线的入射矢量，H矩阵为G中行的入射矩阵，矩阵的行对应信息比特，列对应校验方程，可产生LDPC码。2.BIBD组合法数，则（V，B）构成区组设计，而且对任意一对元素（i，j=l，2，…，v，i式j有个区组同时包括它们，称为平衡不完全区组设计，记为（v，k-BIBD。=1时的BIBD称为Steiner系统。在LDPC码中，关联矩阵和校验矩阵有对应关系。如：假设X=（1，2，3，4，5，6，7，8，9A是12个3元素块的集合，各区组92，6，73，4，81，6，82，4，93，5，7）}X，A）就是一个（9，3，1）-BIBD。它的关联矩阵如下：次次次次次次在BIBD中，任两个元素相遇的次数都是1，构造的校验矩阵可以确保不存在girth 淮南师范学院2012届本科毕业论文114LDPC码的编码本小节主要讨论LDPC码随机构造的几种通用的编码方法，包括密集编码方法、LU分解法、高斯消去法和LDPC码的有效编码方法，并对它们的编码复杂度进行了比较。设校验矩阵H的所有行都是线性无关的。根据分组码定义，编码后的码字x必须满足：Hx=0。为了在接收区易于区分信息位和校验位，一般采用系统码。但对于随机构造的校验矩阵H具有非系统码的形式。因此首先将H进行列变换，将H分割成M×M的左方阵I和M×（N-M）的右矩阵P。同样将x分成M个校验比特C和N-M个系统比特S，有得到AC+BS=0（2）因此，得到校验位C'I1BP该方法的计算复杂度表现在计算I1P，约为O（M3但实际中如果采用相同的校验矩阵，其计算复杂度为M（N-M）。LU分解法的思想是：若I是非奇异矩阵，则可将I分解为一个上三角矩阵U和下三角矩阵L的乘积，其中L和U也是M×N维稀疏矩阵。基本步骤如下：（1）对H矩阵进行LU分解，得到重排后的H、B、L、U。（2）计算Z=BS。（3）通过前向消元法解方程LY=Z，得到Y，其中Y是M维列向量。（4）通过后向消元法解方程UC=Y，得到C。LU分解的一个基本算法如下所示：（1）设U和L为全零矩阵。（2）设F=H。找一个F矩阵中第i行第i列的非零元素，或在其后行（或列）的非零元素。 基于MATLAB的LDPC码的仿真与研究12重新排列F和H的行和列，使得这个非零元素放入第i行第i列。把F矩阵第i列中的从顶到第i行的元素拷贝到U矩阵的第i列。End（4）把B矩阵设置为重排后的H矩阵最后N-M列。高斯消去法把校验矩阵化简为图6所示的等价下三角矩阵。图6下三角校验矩阵构造出系统码c=（u，p使用式（3）计算校验比特。nmHnmpj前面介绍的LDPC编码方法都没有办法使校验矩阵始终保持它的稀疏性，为了能够充分保证和利用校验矩阵的稀疏性，提出了一种有效的编码方法，该编码方法复杂度较低，近似线性复杂度，具体步骤如下：为不改变校验矩阵H的稀疏性，仅对H的行和列进行重组，得到近似下三角稀疏矩阵H，如图7所示： 淮南师范学院2012届本科毕业论文13表示如下：H近似下三角的校验矩阵FDEy*其中，A的维数为（M-g）×（N-MB为（M-g）×g，T为（M-g）×（M-gF过程中，使矩阵的间隔g尽量小。对H左乘矩阵ET1Iy*使其成为右下角为0的矩阵次ABTET1AFET1BD0y*检验ET1BD是否是非奇异，如果是奇异，则与前面列重排，直至为非奇异。次s设码字x=1，其中s为信息符号，p1,p2为校验符号。据方程Hx0可得：2y*AsBp1Tp20（5）(ET1BD)p10（6）得到p 基于MATLAB的LDPC码的仿真与研究14p2T1(AsBp1)((ET1AF)s的计算复杂度矩阵运算计算复杂度AsO（N）T1O（N）O（N）O（N）ET1[As]+FsO（N）(ET1AsFs)的计算复杂度矩阵运算计算复杂度AsO（N）Bp1O（N）As+Bp1O（N）T1O（N）该算法总的计算复杂度为O(Ng2)，优于前几种算法。5LDPC码的译码LDPC码的迭代译码方法是LDPC码能够得以迅速发展的主要原因，该译码方法使得LDPC码不仅描述简单，译码复杂度低，而且可以并行操作，便于硬件实现，具有接近Shannon极限的优异性能[14]。Gallager提出的译码方法可以被当作是Tanner图上的置信传播算法（BP，BeliefPropagationAlgorithm亦可Algorithm）。本小节阐述了概率域的BP算法和对数域的LLRBP算法，然后介绍并分析了几种简化的BP译码算法，并对它们的复杂度和性能进行比较。1概率BP算法 淮南师范学院2012届本科毕业论文15和积（BP）算法[15]是由多变量函数通过因子图中的比特节点和校验节点，以及比特节点和校验节点间的更新规则来得到以某个变量为参数的边缘函数，函数值消息沿边传输，在节点更新。定义如下符号：pP(xnb)为比特节点xnb的概率，L(m)={n：Hmn=1}为参与校验集zm所有比特节点的集合，L(m)\n为除了xn外的L(m)，qn是在件下满足校验zm=0的概率。译码采用基于Tanner图的概率迭代算法。译码过程中，比特节点xn向与其相连的所有校验节点zm发送消息qn，通知zm节点xn处于状态b的概率，用于更新节点zm上的消息rn。然后校验节点zm向父节点xn发送己更新的消息rn，通知它满足校验时节点xn应处于的状态。一轮更新后，产生一个估计码字，若满足H=0，则译码成功；否则再迭代上述过程。如果达到迭代的最大值后仍没有满足校验，则译码失败。可以表述（2）校验节点更新（3）比特节点更新qnmnp'rn'，qnmnp'rn'，选择合适的mn，使qn+qn=1。n）col[m]\(n)n）col[m]\(n)（4）更新伪后验概率n，qmnprn，选择合适的n，使q+q=1。m）M(n)m）M(n)n=1；否则n=0。如果使HT0或到达最大迭代次数， 基于MATLAB的LDPC码的仿真与研究162LLRBP算法LLRBP算法是将概率消息用对数似然比来表示，由此就可以用计算量较小的加法运算来代替计算量较大的乘法运算，这样计算时间就可以大大减少[16]。定义：LLRBP算法计算过程如下：PP（1）初始化首先算出初始概率似然比消息L(Pi)，该消息是由信道传向信息节点的，接着对所有信息节点i和与其相邻的校验节点j，给出信息节点传递给校验节点的初始消息L（2）校验节点更新计算所有的校验节点j第l次迭代时变量节点传递给校验节点L（3）信息节点更新计算所有的信息节点i第l次迭代时校验节点传递给变量节点L)LL1。如果使HiT0或达到最大迭代次数，则结束；否则5.2降低复杂度的BP译码BP算法能够得以普遍应用，不仅因为其具有优异的译码性能，而且因为其具有较等人提出了几种简易算法，这些简易算法都是基于LLRBP算法的，仅仅采用加法运算就可以成功地译码，在很大程度上降低了译码的复杂度[17]。下面加以介绍并进行分析1迭代APP算法 淮南师范学院2012届本科毕业论文17算，但算法的实现复杂度还可进一步降低。如果在要传递的变量节点之间引入相关性，这样传递的就是内部消息而不再是外部消息，需要计算和存储的就仅仅是一个变量消息值，从而可以降低算法的复杂度。迭代APP用LLRBP算法中的L(qi)代替L(qij)参与校验消息的迭代，这样就通过传递内部消息而大大降低了系统的复杂度。2UMPBP-Based算法（最小和或最大积）LLRBP算法中校验节点的处理可以进一步推导为2)所以UMP（UniformlyMostPowerful）BP-Based算法[19]又称为最小和（MinSum）算法或最大积（MaxProduct）算法，该算法是用式来处理LLRBP算法中校验节点的消息，此时校验节点的迭代只有比较算法和加法运算，计算的复杂度就大大降低了。对于加性高斯白噪声信道，该算法不需要信道估计。3APP-Based算法迭代APP算法仅对变量消息的计算进行简化，而UMPBP-Based算法仅对校验消息的计算进行简化。如果对校验消息和变量消息的计算都进行简化，运算的复杂度将会进一步降低，UMP（UniformlyMostPowerful）APP-Based算法就将BP-Based算法与迭代APP算法结合在一起。译码步骤如下：（1）初始化首先算出初始概率似然比消息L(Pi)，该消息是由信道传向信息节点 基于MATLAB的LDPC码的仿真与研究18的，接着对每一变量节点i，设定变量节点传递给校验节点的初始消息（2）校验节点更新计算所有的校验节点j第l次迭代时变量节点传递给校验节点L（3）变量节点更新计算所有的信息节点i第l次迭代时校验节点传递给变量节点L1。如果使HiT0或达到最大迭代次数，则结束；否则转到步骤1）。该算法同样也不需要信道估计。4与标准BP算法的比较为了比较LLRBP算法、UMPBP-Based算法和UMPAPP-Based算法的复杂度，对于码率1/2的LDPC码（n，p，2p对以上三种算法进行了比较，每一次迭代计算量[20]表3译码算法计算量比较译码算法加法运算乘法运算除法运算LLRBPUMPBP-Based UMPAPP-Based 由表3可知，UMPBP-Based算法和UMPAPP-Based算法只有加法运算，跟BP算法比较，计算量大大减小了。同时，为了比较了三种译码算法的误码性能，本论文对三种算法进行了仿真比较，仿真采用的是（1532，3，6）LDPC码，QPSK调制，AWGN信道，仿真结果如图8所示。从仿真结果可以看出，采用两种简化算法，虽然计算复杂度有所降低，但误码性能同时也有一定程度的降低。BP-Based算法只有加法运算，SNR 淮南师范学院2012届本科毕业论文19在误码率为104时比LLRBP算法多损失约0.5dB；而APP-Based算法的计算量约为BP-Based算法的1/2，在误码率为104时又多损失了约0.5dB。相比较而言，BP-Based算法更实用些。图8三种译码算法性能比较 基于MATLAB的LDPC码的仿真与研究206运用MATLAB对LDPC码仿真Matlab是一个大家常用的通信仿真平台，本章利用Matlab软件在其平台上通过仿真得出图形来理论分析说明码长、列重和迭代次数对LDPC码性能的影响从而更好的在通信中运用LDPC码[21]。6.1不同码长对LDPC码的影响应用Matlab软件在列重和迭代次数选择一定值时选取三种不同的码长进行仿真比较对上所述选取LDPC码码长分别为300、500和1000，列重选择2，最大迭代次数设置为20,仿真结果如图9所示。图9不同码长对LDPC码性能的影响从图9的仿真图像可以清楚地看出：在信噪比相同的条件下，LDPC码的性能随着码长的增加而不断提高但是在小信噪比区域，码长的增加对误码率的改进不大，但随着信噪比的增大，LDPC码的误码率得到明显的提高。在误码率为很小时，码长为1000的LDPC码比码长为500和300的LDPC码，信噪比分别降低了约0.3dB和0.6dB。但随着码长的增加，LDPC码性能的提高是相对的，当达到一定码长后，性能将会有很小的提高。这是因为一定码长下编码性能有一定的极限，随着码长的增大，编码和译码的复杂度也增加，编码性能就会更接近极限，性能随码长增加改善的就更少。6.2不同迭代次数对LDPC码的影响应用Matlab软件在列重和码长选择一定值时选取三个不同的迭代次数进行仿真比较对上所述选取的迭代次数分别为10,20,40，列重选择2，码长选择为500,仿真结果如下 淮南师范学院2012届本科毕业论文21图10不同迭代次数对LDPC性能的影响从图10所示的仿真结果可以看出：在相同的信噪比下，LDPC码的性能随着迭代次数的增加而逐渐提高。在误码率为l0时，译码迭代次数为40的LDPC码比迭代次数为20和l0的LDPC码，信噪比分别降低了约0．16dB和0．5dB。但是LDPC码的误码率并不能随着迭代次数的增加无限地减小，当迭代次数足够大的时候，再增加LDPC码的迭代次数，只能增加系统的时延和复杂度，而LDPC码的性能不会再有提高。6.3不同列重对LDPC码的影响应用Matlab软件在迭代次数和码长选择一定值时选取三个不同的列重进行仿真比较，对上所述选取的列重分别为2，3和4迭代次数选择20，码长选择为500,仿真结果如下图11所示。图11不同列重对LDPC性能的影响从图11的仿真结果可以看出，在相同的信噪比下，随着列重的增加，LDPC码的误 基于MATLAB的LDPC码的仿真与研究22码率增大。分析其原因，因为仿真过程中所用的码长不够大，校验矩阵不是足够稀疏，增加列重，会在一定程度上降低检验矩阵的稀疏性，在校验矩阵不是足够稀疏的情况下，稀疏性的一定提高会给编码对应的Tanner图带来大量短长度的环。而短长度的环，尤其是长度为4的环的增加将会使BP算法的性能恶化，导致LDPC码的性能的下降。因此随着列重的增加，LDPC码的性能有所下降。这种LDPC码性能上的下降将会随着码长的增加而逐渐减少，当码长足够长，列重的增加对校验矩阵的稀疏性的影响相对较少，且在译码时，列重大的LDPC码比列重小的LDPC码得到更多的校验信息，从而得到更可靠的译码，所以此时，随着列重的增加，LDPC码的性能将会得到改善。但是当列重增加到较大时，因为校验矩阵不具有稀疏性，性能会随着列重的增加严重下降。本文对LDPC码做了深入的研究，然后应用MATLAB软件对码长、列重和迭代次数对LDPC码性能的影响进行了仿真实验，通过仿真和理论分析得到如下结论：LDPC码长码的误码性能优于短码的误码性能，但当码长达到一定值后，再增加码长，LDPC码的误码率降低的幅度将不大。当码长较小时，增加列重，LDPC码的性能将变差；但当码长足够大时，增加列重，LDPC码的性能将得到改善，但当列重达到一定值时，随着列重增加LDPC码的性能将变差。增加译码迭代次数，LDPC码的性能将得到改善；但当迭代次数足够大时，再增加迭代次数，LDPC码的误码率将不会再降低。因此，为了达到一定的性能，要综合考虑LDPC码的码长、列重和迭代次数。随着科学技术的发展，人们对无线数据和多媒体业务的需求越来越高，Beyond3G以及4G移动通信系统的目标是实现高质量的移动多媒体传输。上世纪90年代末，信道编码领域的一个重大发现是低密度奇偶校验码的再度利用。受到了Turbo码的启示，人们在Gallager六十年代提出的LDPC码的基础上，构造出新型的低密度校验（LDPC）码。LDPC码是一种校验矩阵为稀疏矩阵的线性分组纠错码，采用迭代译码算法时表现出接近Shannon限的优异性能，且译码算法可以并行操作，能够实现数据的高质量传输。本文通过对LDPC码的发展研究状况做了详细的阐述。并且在LDPC码方面，在研究LDPC码编、译码原理的基础上，提出了基于生成矩阵的非正规LDPC码的构造方法，进一步学