基于最小多面体模型的ldpc码线性规划译码方法
【专利摘要】本发明公开了一种基于最小多面体模型的LDPC码线性规划译码方法，主要解决现有LDPC码线性规划译码中译码速度慢和信息传递类译码中存在错误平层的问题。其实现方案是：首先通过分解校验节点的方法将LDPC码的最大似然译码松弛为基于最小多面体的线性规划LP模型，然后利用基于最小多面体的LP模型中矩阵的稀疏性和正交性，建立增广拉格朗日函数并采用交替方向乘子法ADMM算法进行迭代求解得到译码的码字。本发明与现有基于ADMM算法的LP译码方法相比，在不降低LP误码性能的前提下，提高了译码速度，与置信传播BP译码方法比较，在高信噪比下没有出现错误平台，可用于通信技术领域，以提高通信系统译码模块的效率。
【专利说明】
基于最小多面体模型的LDPC码线性规划译码方法
技术领域
[0001] 本发明属于通信领域，特别涉及一种对低密度奇偶校验LDPC码的译码方法，可用 于磁存储、光纤通信和卫星数字视频等领域。
【背景技术】
[0002] 低密度奇偶校验码LDPC是目前一种能够逼近香农信道容量限的最佳编码方案之 一，受到国内外学者研究的关注，并被广泛应用在各种通信领域中。通常采用置信传播BP算 法对LDPC码进行译码。由于BP算法会受到码字结构图中存在的许多有害特性的影响，包括 伪码字、瞬子、陷阱集合等，因而在高信噪比区域，往往受到错误平层的影响，误码率随着信 噪比的增加几乎不再下降。受困于此，在误码性能要求较高的系统中，例如磁存储和光纤通 信中，LDPC码的性能仍不足以满足系统的需求。
[0003] 为了解决上述问题，随后Feldman等首次将最大似然ML译码模型松弛为线性规划 LP译码模型，成功应用于二元线性分组码的译码，奠定了LP译码强大的数学理论支撑。同 时，Feldman还证明了LP译码具有ML特性，码字独立特性和全零假设等良好特性；同时可通 过伪码字图、最小分式距离等分析译码性能;对于校验矩阵中存在短环的LDPC码，LP译码算 法可通过增加冗余校验节点，消除短环影响，提高译码性能。LP译码优点很多，但是译码复 杂度很高，求解困难，严重阻碍了它在实际场景中的应用。
[0004] 为了解决这一难题，Taghavi和Siegel通过向LP模型中增加有效约束，研究出一种 自适应的线性规划ALP译码方法。Xiaojie Zhang等基于这一算法，融合更好的有效割生成 和搜索算法，提出了迭代形式的自适应的线性规划ALP译码算法，不仅降低了复杂度，还提 升了译码性能。另外，Kai Yang等通过精心的度分解设计，研究出一种全新的线性规划译码 模型，相比于Feldman的基本多面体模型，复杂度大大降低，被称为最小多面体MP。
[0005] 以上模型大都通过标准的线性规划工具，如CVX、CPLEX等，通过调用单纯形法或内 点法等方法求解线性规划模型。Barman等将传统的线性规划译码模型与交替方向乘子法 ADMM相结合所提出的迭代式投影译码算法，是目前最好的LP译码方法之一，但其不足是译 码速度较慢。
【发明内容】
[0006] 本发明的目的在于针对上述现有技术的不足，提出一种基于最小多面体模型的 LDPC码线性规划译码方法，以在不降低LP译码性能的前提下，进一步提高LDPC码译码速率， 满足现代无线通信系统的需求。
[0007] 本发明的基本思路是:将LDPC码的最大似然译码通过校验节点分解的方法松弛为 基于最小多面体的LP模型;利用基于最小多面体的LP模型具有稀疏性和正交性特性，采用 分布式并行快速算法ADMM对基于最小多面体的LP模型进行求解，以提高LDPC码译码速率。 其技术方案包括如下：
[0008] (1)将最大似然ML译码模型转化为线性规划LP译码：
[0009] 根据线性规划的定义，利用对数似然将最大似然ML译码模型转化为如下线性规划 LP译码：
[0010] 目标函数:minyTx n
[0011] 约束条件十2 = 0，y/ = l，2,…，m '1、
[0012] xG{0,l}n
[0013] 其中，y表示对数似然比向量，Xi表示发送的第i个码元，i = l ,2, . . . ,n，n表示码 元的总个数，j = l，2, ...，m表示第j个校验节点，m表示校验节点的总个数，x表示译码的码 n 字，（h#)mXn表示m Xn的校验矩阵第j行第i列的数，[/VK ? 2 = 0表示校验方程，"?"表示模 H 2运算。
[0014] (2)分解校验节点，使每个子校验节点的度为3,对每个子校验节点利用奇偶校验 方程构造最小多面体，得到如下最小多面体C:
[0015] C={(xi，X2，X3)} <2>
[0016] 约束条件:xi+X2+X3<2，
[0017] -Xl-X2+X3^;0,
[0018] Xl-X2-X3^;0,
[0019] -Xl+X2-X3^;0,
[0020] XiG[0,l],i = l,2,3
[0021] 其中，X1表示最小多面体的第1个码元变量，X2表示最小多面体的第2个码元变量， X3表示最小多面体的第3个码元变量。
[0022] (3)建立最小多面体的LP译码模型并建立增广拉格朗日函数
[0023] (3a)根据步骤(2)构造的最小多面体，将模型〈1>松驰为如下基于最小多面体的LP 译码模型： 目标函数：tttin
[0024] <3> 约束条件:Ad<b, 0<4 <1
[0025] 其中，q表示扩展后对数似然比向量，T表示矩阵的转置，d表示扩展后的码字，A表 示系数矩阵，b表示系数向量；
[0026] (3b)对基于最小多面体的LP译码模型进行变形，即对式〈3>的不等式约束条件增 加辅助变量w将其转化成等式约束：
[0027] 目标函数:min qTd
[0028] 约束条件:Ad+w = b, <4>
[0029] O^d^l,
[0030] w^O
[0031] (3c)对式〈4>建立增广拉格朗日函数：
<5：>
[0033] 其中，U(d，w，A)表示拉格朗日函数，A表示拉格朗日对偶变量，y表示惩罚参数， ||f + w 表示Ad+w-b的2-范数平方。
[0034] (4)利用ADMM算法对式〈5>中扩展后的码字d，辅助变量w，拉格朗日对偶变量A进行 循环迭代求解，直到满足迭代终止条件，得到最优的扩展码字d'并从中提取出译码的码字 x*〇
[0035] 本方法与现有技术相比具有以下优点：
[0036] 1 ?降低了误码率。
[0037]常用的BP译码由于受到码字结构图中伪码字、陷阱集合等有害特性的影响，在高 信噪比下会出现误码平台，误码率不再降低;本发明与BP译码相比，利用凸优化理论，使码 字具有最大似然ML特性、独立特性和全零特性，减少了码字固有特性的影响，在高信噪比下 不会出现误码平台，保持着较好的瀑区性能，极大地降低了误码率。
[0038] 2.提高了译码速度。
[0039] 传统的基于ADMM算法的线性规划LP译码是线性规划LP译码，比常用的BP译码误码 率性能好，但是译码过程中需要调用投影算法，大大降低了译码速度;本方法相比调用投影 算法的ADMM算法译码，减少了投影算法的操作，且充分利用了矩阵的稀疏性和正交性，因而 在不降低误码性能的前提下，极大地提高了译码速度。
【附图说明】
[0040]图1是本发明的实现流程图；
[0041 ]图2是本发明中进行最小多面体分解原理图；
[0042] 图3是用本发明与现有译码方法对不同规则LDPC码译码的误比特率比较图；
[0043] 图4是用本发明与现有译码方法对不同规则LDPC码译码的平均译码时间比较图。
【具体实施方式】
[0044]以下结合附图对本发明的实施例及效果作进一步详细描述。
[0045] 本实施例是对规则的LDPC码进行信道译码。
[0046] 参照图1，本实例的实现步骤如下：
[0047] 步骤1:根据线性规划的标准型，利用对数似然比向量将最大似然ML译码模型转化 为线性规划LP译码。
[0048] (la)假设发送端发送的二进制低密度奇偶校验码LDPC码字为x={X1，…， Xl，…， xn}，该码字对应的奇偶检验矩阵为H= (hji)mXn，经过噪声为n = {m,…，m,…，nn}的加性高 斯白噪声AWGN信道后，接收的码字为r= {ri，-_，:^，…，rn}，其中xi表示发送的第i个码元，ri 表示接收的第i个码元，i = l，2, . . .，n，n表示码元的总个数，j = l，2, . . .，m表示第j个校验 节点，m表示校验节点的总个数，表示mXn的校验矩阵第j行第i列的数；
[0049] (lb)计算对数似然比向量：y =[ yi, ? ? ?，y i, ? ? ?，yn]T:
[0050] 其中第i个对数似然y i为：
[0052]在加性高斯白噪声AWGN信道中，噪声是均值为0,方差为f的高斯随机变量，服从 正态分布，所以有 <：6:> <7、>
[0055]其中，e表示指数，No表示高斯白噪声功率谱密度；[0056] 根据式〈6>、式〈7>得到：
[0058]最终得到对数似然比向量如下：
[0060] 其中，T表示转置；
[0061] (1幻利用二进制低密度奇偶校验码〇^(：的码字1￡{〇，1}"满足奇偶校验方程 //x ? 2 = 0的条件，得到码字满足每一行奇偶校验的约束条件为： n
[0062] IXV 十 2 = 0
[0063] j = l，2, . . .，m n
[0064] 其中，xG {〇，1}n表示n维向量x中的元素等于〇或1，? 2 = 0表示校验方程，" /-I 由"表示模2运算，E表示求和；
[0065] 得到线性规划LP模型：
[0066]目标函数:minyTx n
[0067]约束条件：Z心不' 十2 = o, vy =-l，2，..,".?n ^8/
[0068] xE{〇，l}n
[0069]步骤2:分解校验节点，使每个子校验节点的度为3,对每个度为3的子校验节点利 用奇偶校验方程构造最小多面体。
[0070] (2a)对于每个度为3的子校验节点，利用奇偶校验方程将一个码字集合E表示为：
[0071] E={(xi,X2JX3)}
[0072] 约束条件：xi+X2+X3<2，
[0073] -xi-X2+X3^0,
[0074] xi-X2~X3^0?
[0075] -Xl+X2-X3^;0,
[0076] xiG {0，1}，i = 1，2,3
[0077] 其中，X1表示最小多面体的第1个码元变量，X2表示最小多面体的第2个码元变量， X3表示最小多面体的第3个码元变量。
[0078] (2b)对于度大于3的校验节点进行分解使每个子校验节点的度为3,本步骤的具体 实现如下：
[0079] (2bl)假设第j个校验节点的度为\多3,则第j个约束条件为
[0080] + …+ + …+ ) ? 2 = 0 ; <9 >
[0081 ]根据hji的取值为0或1，将式〈9>化简为：
[0082] (A:,) +……+ A、）.?2 = 0
[0083] 其中，\表示与第j个校验节点相连的第Ap个码元，AP=1，2,...人，九表示序列 号；
[0084] (2b2)对于与第j个校验节点相连的\个码元，增加辅助码元以分解校验节点：
[0085] 令分解之前码元的个数l(()) = k，对于第v次分解，令增加的辅助码元个数
[0086] 当l(v)为奇数时，满足如下条件：
[0087] (〇) +4;-1} +4r>)?2 = 0， g = l，2，_"，/(r) -1 ,
[0088] 4X)， W.)，
[0089] xxr@2=〇?
[0090] 当1(v)为偶数时，满足如下条件：
[0091 ] ? 2 = 0 , g = l,2,...,/('°, ，v)
[0092] H.rf 十 2 = 0;
[0093] 则每次分解后子校验节点的度都为3,其中与子校验节点相连的码元集合为
表不上界，l〇g2表不以2为底的对数，g表不序列号；
[0094]根据因子图中节点和边的关系，对度为\多3的校验节点，增加\-3个辅助变量，通 过分解得到\_2个度为3的子校验节点；
[0095]参照图2,校验节点的度为6,与校验节点相连的原始码元分别为if1，xf，xf， aP，Af"，，分解校验节点使每个子校验节点的度为3,则需要增加3个辅助码元，分别 为xP，411，4"，将相邻的两个原始码元xf和xf与增加的辅助码元相结合，构成一个 码元集合{(AfUi'Af)}使子校验节点的度为3;将相邻的两个原始码元，和#与增加的 辅助码元相结合，构成一个码元集合使子校验节点的度为3;将相邻的两 个原始码元;#和与增加的辅助码元xp相结合，构成一个码元集合使子 校验节点的度为3;将增加的3个辅助码元xP，X!11，#相结合，构成一个码元集合 {(g'x丨"，x 3(1))}使子校验节点的度为3;对于度为6的校验节点，通过分解得到4个度为3的子 校验节点；
[0096] (2c)将变量的取值范围Xi G {〇，1}松弛为线性约束Xi G [0，1 ]，得到如下最小多面 体C:
[0097] C={(xi，X2，X3)} <10>
[0098] 约束条件：xi+X2+X3<2,
[0099] -xi-X2+X3^；0,
[0100] Xl-X2-X3^;0,
[0101] -Xl+X2-X3^;0,
[0102] XiG[0,l],i = l,2,3,
[0103] 其中，[0,1]表示0到1。
[0104] 步骤3:建立最小多面体的LP译码模型并建立增广拉格朗日函数。
[0105] (3a)根据步骤2构造的最小多面体，建立最小多面体的LP译码模型：
[0106] (3ai)定义= I^(a；-3)为辅助变量的总个数，r; = 2^(1, -2)为分解出的最小 多面体的总个数，将原始变量x和辅助变量合并扩展为d = [X，i/lxrJ7'，将对数似然比向量 扩展为:，则式〈8>中的目标函数转化为min qTd，其中，表示1行ra列的向量， 九ra表示1行r。列的向量全为0;
[0107] (3a2)对于第Y。个最小多面体，假设与扩展后的码字d相连的码元变量为<， ，则定义其相对应的矩阵为Qi，根据式〈1〇>利用线性方程组的矩阵形式，令不等式 右侧的值用向量g表示，即g=.t2,0s0s0f，不等式左侧的系数用矩阵F表示，即 ~ 1 1 1' F= | j j，则第y。个最小多面体的矩阵形式为FQ; ;对于r。个最小多面体，令系 -1 1 -1 数矩阵4= [(FQf s(FQys f.，(FQrt ff,系数向量6 = 1? kxl。则式〈8>中的约束条件 转化为Ad彡b，0彡d彡1，其中，^4表示第y。个最小多面体中与扩展后的码字d相连的第一个 码元，表示第y。个最小多面体中与扩展后的码字d相连的第二个码元，表示第y。个最 小多面体中与扩展后的码字d相连的第三个码元，矩阵Q&中只有(U/:! )，（2,J;2)，（3, 对应位置的元素为1，其他元素为零，"表示笛卡尔积，:kx!表示长度为r。的元素全为1， FQ\表示第y。个最小多面体的系数矩阵，由于系数矩阵F中有12个非零元素且任意两列是 相互正交的，且系数矩阵与系数矩阵F相比，只是增加了（n+ r a-3)个全零列向量，所以 FQ:中只有12个非零元素且任意两列也是相互正交的；由于系数矩阵A是由r c个最小多面 体直接级联得到的，无需改变FQk.中任意两列的正交关系，所以系数矩阵A具有正交性；由 于系数矩阵A有4 r。X (n+r a)个元素，其中只有12 r。个非零元素，所以系数矩阵A具有稀疏 性；
[0108] (3a3)根据(3al)和(3a2)，得到最小多面体的LP译码模型： 目标函数：.min xx
[0109] .... ：<ll> 约束条件：Ad<b, 0<d<l
[0110] (3b)对最小多面体LP译码模型进行变形，即对式〈11>的不等式约束条件增加辅助 变量w将其转化成等式约束：
[0111] 目标函数:min qTd
[0112] 约束条件:Ad+w = b, <12>
[0113] O^d^l,
[0114] w^0
[0115] (3c)对式〈12>建立增广拉格朗日函数：
[0116] = qTd + f{Ad + w -h) + -^llj
[0117]其中，U(d，w，A)表示拉格朗日函数，A表示拉格朗日对偶变量，y表示惩罚参数， + wA【表示Ad+w-b的2 -范数平方。
[0118] 步骤4:利用ADMM算法对扩展后的码字d，辅助变量w，拉格朗日对偶变量A进行循环 迭代求解，直到满足迭代终止条件，得到最优的扩展码字d'并从中提取出译码的码字/。
[0119] (4a)利用ADMM算法的如下迭代更新公式，求解第k+1次迭代后的扩展后的码字dk +1，辅助变量#+1，拉格朗日对偶变量Ak+1:
[0120] </'' ' = arg min L,,(d. w1'. /') <13x 0<^<l
[0121] '1 = argmin L.Xd1' 1, w, /A) \M>
[0122] Ak+1 = Ak+y(Adk+1+wk+1-b) <15>
[0123] (4al)对扩展后的码字d进行更新，即固定辅助变量wk和拉格朗日对偶变量Ak，对式 〈13>中扩展后的码字d求导，并令导数等于零，得到更新后的码字d k+1:
<16>
[0125]其中，表示在超立方体[0,ir~上的投影操作，在求解式〈16>时，利用系数 矩阵A的稀疏性和正交性大大降低了计算复杂度，提高了译码速度；
[0126] (4a2)对辅助变量w进行更新，即固定更新后的码字dk+1和拉格朗日对偶变量Ak，对 式〈14>中的辅助变量w求导，并令导数等于零，得到更新后的辅助变量wk+1: <17>
[0128] 其中，nw>Q表示在w>0上的投影操作；
[0129] (4a3)根据式〈15>利用(4al)更新后的码字dk+1和(4a2)更新后的辅助变量w k+1得到 更新后的拉格朗日对偶变量Ak+1;
[0130] (4b)定义第k+1次迭代后原始残差Rk+1 = Adk+1+wk+1-b，对偶残差Sk+1 = wk+1-wk，在迭 代求解过程中，当原始残差2-范数的平方和对偶残差2-范数的平方同时小于 阈值1(T5时停止迭代，得到最优的扩展码字d'并从中提取出译码的码字，。
[0131] 本发明的效果通过以下仿真结果进一步说明：
[0132] 仿真方法:本发明、现有基于ADMM算法的线性规划LP译码方法、置信传播BP译码方 法。
[0133] 仿真1:用本发明和现有的两种方法分别对不同的规则LDPC码进行译码，比较其误 比特率BER，其结果如图3所示；
[0134] 由图3可见，采用本发明和现有基于ADMM算法的LP译码方法对（160，80)规则LDPC 码进行译码，得到的误比特率BER曲线基本吻合；同样，采用这两种译码方式对(512,256)规 则LDPC码进行译码，得到的误比特率BER曲线也基本重合;说明本发明在误码性能方面可以 到达和现有基于AD丽算法的LP译码方法一样的效果；同样，对（160，80 )、（ 512，256)规则 LDPC码分别采用BP译码方法，都会出现了错误平层，而本发明仍保持较好的瀑区性能，没有 出现错误平层，说明本发明保持了 LP译码低错误平层的优点；
[0135] 仿真2:用本发明和现有的两种方法分别对不同的规则LDPC码进行译码，比较其平 均译码时间，其结果如图4所示；
[0136] 由图4可见，采用本发明和现有基于ADMM算法的LP译码方法对（160，80)规则LDPC 码进行译码，得到本发明的平均译码时间短；同样，采用这两种译码方式对(512,256)规则 LDPC码进行译码，得到本发明的平均译码时间短;说明本发明的平均译码速度比现有基于 ADMM算法的LP译码方法的平均译码速度快；同样，对（160，80)、（512，256)规则LDPC码分别 采用BP译码方法，得到本发明的平均译码时间比BP译码方法的平均译码时间短，说明本发 明是快速的有效译码方法。
【主权项】
1.基于最小多面体模型的LDPC码线性规划译码方法包括： (1) 将最大似然ML译码模型转化为线性规划LP译码： 根据线性规划的定义，利用对数似然将最大似然ML译码模型转化为如下线性规划LP译 码：其中，γ表示对数似然比向量，xi表示发送的第i个码元，i = l，2,. . .，n，n表示码元的总 个数，j = l，2, ...，m表示第j个校验节点，m表示校验节点的总个数，X表示译码的码字， (1^)ηΧη表示mXn的校验矩阵第j行第i列的数，$>7Λ?2 = 0表示校验方程，" ? "表示模2 卜1 运算。 (2) 分解校验节点，使每个子校验节点的度为3,对每个子校验节点利用奇偶校验方程 构造最小多面体，得到如下最小多面体C: C= {(xi,X2,X3)} <2> 约束条件:Xl+X2+X3<2 , -Χ1-Χ2+Χ3^;0 , Χ1-Χ2-Χ3^?Ο, -Χ1+Χ2-Χ3^;0, Χ?￡[0,1],? = 1,2,3 其中，XI表示最小多面体的第1个码元变量，Χ2表示最小多面体的第2个码元变量，Χ3表 示最小多面体的第3个码元变量。 (3) 建立最小多面体的LP译码模型并建立增广拉格朗日函数 (3a)根据步骤(2)构造的最小多面体，将模型〈1>松驰为如下基于最小多面体的LP译 码模型：<：3.> 其中，q表示扩展后对数似然比向量，T表示矩阵的转置，d表示扩展后的码字，A表示系 数矩阵，b表示系数向量； (3b)对基于最小多面体的LP译码模型进行变形，即对式〈3>的不等式约束条件增加辅 助变量w将其转化成等式约束：(3c)对式〈4>厘11削悝份朗Η圈数： <：4>其中，Lw(d，w，A)表示拉格朗日函数，λ表示拉格朗日对偶变量，μ表示惩罚参数， + W .-.6忙表示Ad+w-b的2-范数平方。 (4)利用ADMM算法对式〈5>中扩展后的码字d，辅助变量w，拉格朗日对偶变量λ进行循环 迭代求解，直到满足迭代终止条件，得到最优的扩展码字d'并从中提取出译码的码字/。2. 根据权利要求1所述的方法，其中步骤(1)中利用对数似然比将最大似然ML译码模型 转化为LP译码，按如下步骤进行： (la) 假设发送端发送的二进制低密度奇偶校验码LDPC码字为χ= {χι，…，Xi，…，χη}，该 码字对应的奇偶检验矩阵为H= (hji)mXn，经过噪声为n= {m，…，m，…，ηη}的加性高斯白噪 声AWGN信道后，接收的码字为r={ri,…，ri，···，r n}，ri表示接收的第i个码元； (lb) 计算对数似然比向量γ = [ γι, · · ·，yi, · · ·，γη]τ:其中，No为高斯白噪声功率谱密度； (lc) 利用二进制低密度奇偶校验码LDPC的码字xe {〇，l}n满足奇偶校验方程 /&領2= 0的条件，得到码字满足每一行奇偶校验的约束条件为： J .丄.5〃*r 从码字集合中寻找一个码字使目标函数最小的码字，得到线性规划模型： 目标函数:η?ηγτχ 约束条件 ?-1xe {0，ι}ηο3. 根据权利要求1所述的方法，其中步骤(2)对每个子校验节点利用奇偶校验方程构造 最小多面体，按如下步骤进行： (2a)对于每个度为3的子校验节点，利用奇偶校验方程将一个码字集合Ε表示为： E= {(X1,X2,X3)} 约束条件:Xl+X2+X3<2 , -Χ1-Χ2+Χ3^;0 , Χ1-Χ2-Χ3^;0, -Χ1+Χ2-Χ3^;0, Xie{〇，l}，i = l，2,3 (2b)对于度大于3的校验节点进行分解使每个子校验节点的度为3,假设第j个校验节 <7> 点的度为多3，则第j个约束条件为 根据h#的取值为0或1，将式〈7>化简为： 其中，l表示与第j个校細·ρ好、相进tw弟兀，λρ=i，z，…，λ」，八表示序列号； 根据因子图中节点和边的关系，对度为\多3的校验节点，增加\-3个辅助变量，通过分 解得到\_2个度为3的子校验节点； (2c)将变量的取值范围Xi e {〇，1}松弛为线性约束Xi e [〇，1 ]，得到如下最小多面体c: C= {(X1,X2,X3)} 约束条件:Xl+X2+X3<2 , -Χ1-Χ2+Χ3^;0 , Χ1-Χ2-Χ3^?Ο, -Χ1+Χ2-Χ3^;0, xie[0,l],i = l,2,3〇4. 根据权利要求1所述的方法，其中步骤(3a)中建立最小多面体的LP译码模型，按如下 步骤进行： (3al)定义1；=[；；从-3)为辅助变量的总个数，Γ^Σ;(4-2)为分解出的最小多面 体的总个数，将原始变量X和辅助变量合并扩展为d = [X, ?lxr/，将对数似然比向量扩展 为F [7,九F，则式〈1>中的目标函数转化为min qTd; (3a2)对于Γ。个最小多面体，利用线性方程组的矩阵形式，用矩阵A表示不等式左侧的 系数，用向量b表示不等式右侧的值，则式〈1>中的约束条件转化为Ad<b，0<d彡1; (3a3)根据(3al)和(3a2)，得到最小多面体的LP译码模型： 目标函数:min qTd 约束条件:Ad彡b，0彡d彡1。5. 根据权利要求1所述的方法，其中步骤(4)中利用ADMM算法求解码字集合，按如下步 骤进行： (4a)利用Α·Μ算法的如下迭代更新公式，求解第k+Ι次迭代后的扩展后的码字dk+1，辅 助变量wk+1，拉格朗日对偶变量Ak+1:<7> <8/ Ak+1 = Ak+y(Adk+1+wk+1-b) <9> (4al)对扩展后的码字d进行更新，即固定辅助变量wk和拉格朗日对偶变量，对式〈7> 中扩展后的码字d求导，并令导数等于零，得到更新后的码字dk+1:<10> 其中，示在超立方体[〇,丨]"11 '上的投影操作，T表示矩阵的转置； (4a2)对辅助变量w进行更新，即固定更新后的码字dk+1和拉格朗日对偶变量，对式〈8> 中的辅助变量w求导，并令导数等于零，得到更新后的辅助变量wk+1:<11> - μ 其中，IIwX)表示在w>0上的投影操作； (4a3)根据式〈9>利用(4al)更新后的码字dk+1和(4a2)更新后的辅助变量wk+1得到更新 后的拉格朗日对偶变量11{+1; (4b)定义第k+Ι次迭代后原始残差Rk+1 = Adk+1+wk+1-b，对偶残差Sk+1 = wk+1-wk，在迭代求 解过程中，当原始残差2-范数的平方||#13和对偶残差2-范数的平方同时小于阈值 10-5时停止迭代，得到最优的扩展码字d'并从中提取出译码的码字，。
【文档编号】H03M13/11GK105959015SQ201610255059
【公开日】2016年9月21日
【申请日】2016年4月22日
【发明人】王勇超, 白晶, 杜倩
【申请人】西安电子科技大学