五边类型低密度奇偶校验码的外部信息转移图的构造方法与流程

本发明属于数字通信
技术领域：
，具体涉及一种五边类型低密度奇偶校验码的外部信息转移图的构造方法。
背景技术：
：1948年香农《amathematicaltheoryofcommunication》的这篇论文的提出，标志着信息论的产生。人们开始不断深入研究信道编码理论，接着研究出不少性能靠近香农限的好码。其中值得引起关注的是tomrichardson等在2004年论文中提出的一种五边类型低密度奇偶校验码(t.j.richardsonandr.l.urbanke.five-edgetypeldpccodes)。五边类型低密度奇偶校验码(five-edge_type_ldpc_codes，即fet-ldpc码)具有独特的因子图结构，其编译码复杂度低、结构灵活且性能优异。首先，fet-ldpc码在其因子图中引入五种边类，使其结构更加灵活。其次，fet-ldpc码还在其因子图的结构中引入了度为2的变量节点，度为1的变量节点和删余变量节点。这些结构的设计，使fet-ldpc码的编译码更加灵活且复杂度低，码字的译码性能也有所提高。fet-ldpc码的译码由以度为1的变量节点为主的单校验(即spc)检测器和以度为2的变量节点为主的累积码(即acc)检测器相互迭代完成。在信道编码的研究过程中，码字的迭代译码结构引起人们的广泛关注。为了对迭代译码结构的收敛行为进行更好的评估，可以用码字的译码门限值来反映码字性能的好坏。近年来，密度进化(densityevolution，即de)算法被用于分析码字的门限值。de算法已被成功地用在ldpc码、fet-ldpc码等中。但是在实际应用中，de算法非常高的计算复杂度限制了它的发展。为了解决这一问题，stephantenbrink在其论文《“convergenceofiterativedecoding》中提出了一种新的方法，即外部信息转移(extrinsicinformationtransfer，即exit)图。exit图描述了迭代译码过程中输入信息与输出信息之间的关系。解码轨迹显示了exit转移图中各解码器之间的外部信息交换。exit图作为一种新的分析工具，已经被成功的用在了并行级联码(pc码)、重复累积码(ra码)、低密度奇偶校验码(ldpc码)和原模图ldpc码中。但是作为性能更加优异和结构更加灵活的fet-ldpc码，目前的研究中还缺乏对fet-ldpc码外部信息转移图的设计。如果仅仅把传统的exit方法直接运用在fet-ldpc码中，对其码字的性能分析并不十分准确。fet-ldpc码中度为2的变量节点，度为1的变量节点和删余变量节点使其特有的exit图的设计变得更加灵活。技术实现要素：本发明的目的之一在于克服以上缺点，提供一种五边类型低密度奇偶校验码的外部信息转移图的构造方法，既能对fet-ldpc码的性能分析更为准确，同时能降低门限值估计的误差以及降低计算的复杂度。为了解决上述技术问题，本发明提供了一种五边类型低密度奇偶校验码的外部信息转移图的构造方法，包括以下步骤：步骤0、计算高斯白噪声信道的方差计算公式为：其中，r为所述五边类型低密度奇偶校验码的码率，snr为信噪比；步骤1、计算变量节点的信道似然信息的方差计算公式为：若五边类型低密度奇偶校验码中的变量节点通过信道时，若五边类型低密度奇偶校验码中的变量节点未通过信道时，步骤2、计算所述五边类型低密度奇偶校验码对应的spc检测器的fexit曲线，即fexit图的上界曲线upperbound(iav，iev)；步骤3、计算所述五边类型低密度奇偶校验码对应的acc检测器的fexit曲线，即fexit图的下界曲线lowerbound(iec，iac)；步骤4、根据所述fexit图的上界曲线以及所述fexit图的下界曲线，构造所述五边类型低密度奇偶校验码的fexit的曲线带。进一步地，所述的五边类型低密度奇偶校验码的外部信息转移图的构造方法，还包括步骤5、根据所述fexit的上界曲线以及所述fexit图的下界曲线，估算fet-ldpc码的门限。进一步地，所述“计算所述五边类型低密度奇偶校验码对应的spc检测器的fexit曲线，即fexit图的上界曲线upperbound(iav，iev)”，包括以下步骤：步骤2.1、将t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav初始化为0；步骤2.2、计算t3边上由信息位变量节点v2产生的流向校验节点c1(v2→c1)的似然信息与v2之间的平均互信息iev3，计算公式为：其中，dv2表示与信息位变量节点v2相连的t2类边的度，j函数为单调递增函数；步骤2.3、初始化spc迭代计数器it1为0；步骤2.4、计算t3边上从信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的先验互信息iac3，计算公式为：iac3＝iev3；步骤2.5、计算t4边上从单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的先验互信息iac4，计算公式为：其中，iev4表示t4边上由单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的平均互信息，j函数为单调递增函数；步骤2.6、计算t3边上由校验节点c1产生的流向信息位变量节点v2(c1→v2)的似然信息与v2之间的平均互信息iec3，计算公式为：其中，dc3表示与校验节点c1相连的t3类边的度，dc4表示与校验节点c1相连的t4类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.7、计算t3边上从校验节点c1传给信息位变量节点v2(c1→v2)的似然信息与v2之间的先验互信息iav3，计算公式为：iav3＝iec3；步骤2.8、更新t3边上由信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的平均互信息为iev3，计算公式为：其中，dv2表示与信息位变量节点v2相连的t2类边的度；dv3表示与信息位变量节点v2相连的t3类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.9、根据公式it1＝it1+1计算，再判断若it1<3，则跳至步骤2.4，否则跳至步骤2.10。步骤2.10、计算t2边上由信息位变量节点v2传给校验节点c2和c3(v2→c2、c3)的似然信息与v2之间的平均互信息iev2，计算公式为：其中，dv2表示与信息位变量节点v2相连的t2类边的度，dv3表示与信息位变量节点v2相连的t3类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.11、计算t1边上由信息位变量节点v1传给校验节点c2和c3(v1→c2、c3)的似然信息与v1之间的平均互信息iev1，计算公式为：其中，dv1表示与信息位变量节点相连的t1类边的度，j-1函数是j函数的反函数；步骤2.12、计算变量节点v1中t1类边的边数占信息位变量节点总边数的比例r1，计算公式为：计算变量节点v2中t2类边的边数占信息位变量节点总边数的比例r2，计算公式为：其中，dv1表示与信息位变量节点相连的t1类边的度，dv2表示与信息位变量节点v2相连的t2类边的度，v1表示与t1类边相连的信息位变量节点v1的个数占总的变量节点个数的比例，v2表示与t2类边相连的信息位变量节点v2的个数占总的变量节点个数的比例；计算由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev，计算公式为：iev＝r1×iev1+r2×iev2；步骤2.13、根据公式iav＝iav+0.0001进行计算，再判断如果iav≤1，则跳转至步骤2.2，否则跳转至步骤2.14；步骤2.14、根据上述步骤中求得的所有iav和iev取值的对应关系，得到坐标轴中fexit图的上界曲线upperbound(iav，iev)，其中，横轴坐标为t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav，纵轴坐标为由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev。进一步地，所述“计算所述五边类型低密度奇偶校验码对应的acc检测器的fexit曲线，即fexit图的下界曲线lowerbound(iec，iac)”，包括以下步骤：步骤3.1、将t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac初始化为0；步骤3.2、计算t5边上由校验节点c2和c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的平均互信息iec5，计算公式为：其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.3、初始化acc迭代计数器it2为0；步骤3.4、计算t5边上从校验节点c2、c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的先验互信息iav5，计算公式为：iav5＝iec5；步骤3.5、计算t5边上从累积码v4传给校验节点c2、c3(v4→c2、c3)的似然信息与v4之间的平均互信息iev5，计算公式为：其中j函数为单调递增函数；步骤3.6、计算t5边上从累积码v4传给校验节点c2和c3(v4→c2、c3)的似然信息与v4之间的先验互信息iac5，计算公式为：iac5＝iev5；步骤3.7、更新t5边上由校验节点c2和c3传给累积码v4的似然信息与v4之间的平均互信息为iec5，计算公式为：其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，dc5表示与校验节点c2和c3相连的t5类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.8、根据公式it2＝it2+1计算，并判断若it2<3，则跳至步骤3.4，否则跳至步骤3.9；步骤3.9、更新由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec，计算方式为：其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.10、根据公式iac＝iac+0.0001计算，再判断如果iac≤1，则跳转至步骤3.2，否则跳转至步骤3.11；步骤3.11、根据上述步骤求得所有的iac和iec的对应关系，利用线性插值法求出当iec为已知iec时与之对应的iac的具体值iac，计算公式为：iac＝interp1(iec,iac,iec,'linear'),其中0≤iec≤1；步骤3.12、根据iec和iac的取值，得到坐标轴中fexit的下界曲线lowerbound(iec，iac)，其中，横轴坐标为由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec，纵轴坐标为t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac。进一步地，所述“根据所述fexit的上界曲线以及所述fexit图的下界曲线，估算fet-ldpc码的门限”，包括以下步骤：步骤5.1、初始化高斯白噪声方差σn为σ，初始化最小差值max_dif为100，初始化iav＝iec＝0，初始化门限取值范围的最大值σmax；步骤5.2、根据所述步骤2的方法求得iev；根据所述步骤3的方法求得iac；步骤5.3、计算iev与iac的差值dif，计算公式为：dif＝iev-iac；步骤5.4、判断若dif小于等于0，则跳至步骤5.5，否则进一步判断如果dif小于等于max_dif，则将max_idf设置为dif，并将门限值threshold设置为σn，否则如果dif大于max_dif，跳至步骤5.5；步骤5.5、根据公式σn＝σn+0.0001进行计算，并判断若σn≤σmax，则跳至步骤5.2；否则跳至步骤5.6；步骤5.6、获取当前的threshold值作为fet-ldpc码的门限值。综上所述，本发明技术方案的有益效果有:1.根据五边类型低密度奇偶校验码的结构以及译码迭代过程的特点，分别对五边类型低密度奇偶校验码的spc检测器和acc检测器中输入信息和输出信息的关系进行计算，构建fexit曲线图，提供更为准确性能分析依据。2.通过fexit曲线图估计五边类型低密度奇偶校验码的门限值误差很小，同时计算复杂度也比现有的de算法更低，更加适用于实际的应用场景中。附图说明图1是本发明的fet-ldpc码fexit互信息流图。图2是本发明的多边类型低密度奇偶校验码的外部信息转移图实现流程图。图3是本发明计算fexit图的上界曲线upperbound的步骤流程图。图4是本发明计算fexit图的下界曲线lowerbound的步骤流程图。图5是本发明fet-ldpc码门限值的计算步骤流程图。图6是本发明实施例在snr＝0.5db时，码率r＝0.5时，fet-ldpc码对应的fexit图的上界曲线。图7是本发明实施例在不同snr下，码率r＝0.5时，fet-ldpc码对应的fexit图的上界曲线。图8是本发明实施例在snr＝0.5db时，码率r＝0.5时，fet-ldpc码对应的fexit图的下界曲线。图9是本发明实施例在不同snr下，码率r＝0.5时，fet-ldpc码对应的fexit图的下界曲线。图10是本发明实施例在snr＝0.5db时，码率r＝0.5时，由fexit图的上界曲线和下界曲线构成的fet-ldpc码的fexit图。图11是本发明门实施例在限值σn＝0.9700时，得到的fet-ldpc码的fexit图。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。五边类型低密度奇偶校验码(fet-ldpc码)的结构特点是在其因子图中引入了五种类型的边，分别用t1、t2、t3、t4和t5表示。根据变量节点与不同种类的边的连接情况，变量节点被划分成四类，分别用v1、v2、v3和v4表示，其中v1和v2是信息位变量节点，v3是单校验码，v4是累积码。同时，根据校验节点与不同种类的边的连接情况，校验节点被划分成三类，分别用c1、c2和c3表示，其中c1通过t3和t4这两种边把信息位变量节点v2和单校验码v3连接起来，c2和c3通过t1、t2、t5这三种类的边把信息位变量节点v1、v2和累积码v4连接起来。因此，fet-ldpc码的因子图可以通过下述两个多项式来描述：其中公式1对应于fet-ldpc码的变量节点，公式2对应于fet-ldpc码的校验节点。公式中的相关变量含义为：rm表示变量节点是否经过删余信道，当m＝0时，表示该变量节点经过删余信道；当m＝1时，表示该变量节点经过awgn信道。vi表示变量节点vi的数量占总的变量节点数量的比例，1≤i≤4。表示变量节点与五种边的连接情况，其中x＝(x1,x2,x3,x4,x5)，dv＝(dv1,dv2,dv3,dv4,dv5)，dvt表示变量节点与第t种类边相连的边度，1≤t≤5，且有dv4＝1，dv5＝2。例如表示变量节点v2与第2种类边相连的度为dv2和第3种类边相连的度为dv3。cj表示校验节点cj的数量占总的变量节点数量的比例，1≤j≤3。表示校验节点与五种边的连接情况，其中x＝(x1,x2,x3,x4,x5)，dc＝(dc1,dc2,dc3,dc4,dc5)，dct表示该校验节点与第t种边相连的边度，1≤t≤5，其中dc4＝1，dc5＝2。根据fet-ldpc码的上述结构特点，可以知道fet-ldpc码的译码迭代过程如下：通过t3类边和t4类边，信息位变量节点v2、校验节点c1和单校验码v3之间形成一个spc迭代，我们可以把这一部分称为spc检测器。通过t1类边、t2类边和t5类边，信息位变量节点v1和v2、校验节点c2、c3和累积码v4之间形成一个acc迭代，我们可以把这一部分称为acc检测器。spc检测器和acc检测器之间又构成一个大迭代。迭代完成后，对信息位v1和v2进行判决，译码结束。根据fet-ldpc码的迭代过程，本发明提出的一种fet-ldpc码的多边类型低密度奇偶校验码的外部信息转移(five-edge-type-ldpcextrinsicinformationtransfer，即fexit)图的构造方法，如图1是本发明的fet-ldpc码fexit互信息流图，由图中可以看出fet-ldpc码的fexit信息流的转移方向，fexit中的信息流主要通过spc检测器和acc检测器之间迭代传输得到。在spc检测器中，第一、互信息通过t1、t2边从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)；第二、互信息通过t3边由信息位变量节点v2传给校验节点c1(v2→c1)；第三、似然信息通过t4边从校验节点c1传给单校验码v3(c1→v3)；第四、互信息从单校验码v3传给校验节点c1(v3→c1)；第五、互信息通过t3边从校验节点c1流向信息位变量节点v2(c1→v2)。这就完成一个spc检测器迭代。在acc检测器中，第一、互信息通过t1、t2边从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)；第二、互信息通过t5边由校验节点c2和c3传给累积码v4(c2、c3→v4)；第三、互信息通过t5边从累积码v4传给校验节点c2、c3(v4→c2、c3)；最后互信息由校验节点c2和c3传给信息位变量节点v1和v2。这就完成一个acc检测器迭代。如图2，是本发明的多边类型低密度奇偶校验码的外部信息转移图实现流程图，包括以下步骤：步骤0、计算高斯白噪声信道的方差计算公式为：其中，r为所述五边类型低密度奇偶校验码的码率，snr为信噪比；步骤1、计算变量节点的信道似然信息的方差计算公式为：若五边类型低密度奇偶校验码中的变量节点通过信道时，若五边类型低密度奇偶校验码中的变量节点未通过信道时，在信道编码中的专业术语中，若五边类型低密度奇偶校验码中的变量节点不经过信道传输，也可以称为该五边类型低密度奇偶校验码中的变量节点被删余。步骤2、计算所述五边类型低密度奇偶校验码对应的spc检测器的fexit曲线，即fexit图的上界曲线upperbound(iav，iev)；在这种情况下，本方法不仅要关注fet-ldpc码中的信息位变量节点，还要关注单校验码v3给信息位变量节点v2所提供的先验信息，即t3类边和t4类边上的信息。如图3，是本发明计算fexit图的上界曲线upperbound的步骤流程图，包括以下步骤:步骤2.1、将t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav初始化为0；步骤2.2、计算t3边上由信息位变量节点v2产生的流向校验节点c1(v2→c1)的似然信息与v2之间的平均互信息iev3，计算公式为：其中，dv2表示与信息位变量节点v2相连的t2类边的度，j函数为单调递增函数；步骤2.3、初始化spc迭代计数器it1为0；步骤2.4、计算t3边上从信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的先验互信息iac3，计算公式为：iac3＝iev3；步骤2.5、计算t4边上从单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的先验互信息iac4，计算公式为：其中，iev4表示t4边上由单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的平均互信息，j函数为单调递增函数；步骤2.6、计算t3边上由校验节点c1产生的流向信息位变量节点v2(c1→v2)的似然信息与v2之间的平均互信息iec3，计算公式为：其中，dc3表示与校验节点c1相连的t3类边的度，dc4表示与校验节点c1相连的t4类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.7、计算t3边上从校验节点c1传给信息位变量节点v2(c1→v2)的似然信息与v2之间的先验互信息iav3，计算公式为：iav3＝iec3；步骤2.8、更新t3边上由信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的平均互信息为iev3，计算公式为：其中,dv2表示与信息位变量节点v2相连的t2类边的度，dv3表示与信息位变量节点v2相连的t3类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.9、根据公式it1＝it1+1计算，再判断若it1<3，则跳至步骤2.4，否则跳至步骤2.10。步骤2.10、计算t2边上由信息位变量节点v2传给校验节点c2和c3(v2→c2、c3)的似然信息与v2之间的平均互信息iev2，计算公式为：其中,dv2表示与信息位变量节点v2相连的t2类边的度，dv3表示与信息位变量节点v2相连的t3类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤2.11、计算t1边上由信息位变量节点v1传给校验节点c2和c3(v1→c2、c3)的似然信息与v1之间的平均互信息iev1，计算公式为：其中，dv1表示与信息位变量节点相连的t1类边的度，j-1函数是j函数的反函数；步骤2.12、计算变量节点v1中t1类边的边数占信息位变量节点总边数的比例r1，计算公式为：计算变量节点v2中t2类边的边数占信息位变量节点总边数的比例r2，计算公式为：其中，dv1表示与信息位变量节点相连的t1类边的度，dv2表示与信息位变量节点v2相连的t2类边的度，v1表示与t1类边相连的信息位变量节点v1的个数占总的变量节点个数的比例，v2表示与t2类边相连的信息位变量节点v2的个数占总的变量节点个数的比例；r1+r2＝1；计算由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev，计算公式为：iev＝r1×iev1+r2×iev2；步骤2.13、根据公式iav＝iav+0.0001进行计算，再判断如果iav≤1，则跳转至步骤2.2，否则跳转至步骤2.14；步骤2.14、根据上述步骤中求得的所有iav和iev取值的对应关系，得到坐标轴中fexit图的上界曲线upperbound(iav，iev)，其中，横轴坐标为t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav，纵轴坐标为由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev。步骤3、计算所述五边类型低密度奇偶校验码对应的acc检测器的fexit曲线，即fexit图的下界曲线lowerbound(iec，iac)；这时需要关注校验节点c2和c3和累积码v4之间的输入输出关系，即t1、t2、t5边上信息之间的输入输出关系。如图4，是本发明计算fexit图的下界曲线lowerbound的步骤流程图，包括以下步骤：步骤3.1、将t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac初始化为0；步骤3.2、计算t5边上由校验节点c2和c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的平均互信息iec5，计算公式为：其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.3、初始化acc迭代计数器it2为0；步骤3.4、计算t5边上从校验节点c2、c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的先验互信息iav5，计算公式为：iav5＝iec5；步骤3.5、计算t5边上从累积码v4传给校验节点c2、c3(v4→c2、c3)的似然信息与v4之间的平均互信息iev5，计算公式为：其中j函数为单调递增函数；步骤3.6、计算t5边上从累积码v4传给校验节点c2和c3(v4→c2、c3)的似然信息与v4之间的先验互信息iac5，计算公式为：iac5＝iev5；步骤3.7、更新t5边上由校验节点c2和c3传给累积码v4的似然信息与v4之间的平均互信息为iec5，计算公式为：其中，其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，dc5表示与校验节点c2和c3相连的t5类边的度，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.8、根据公式it2＝it2+1计算，并判断若it2<3，则跳至步骤3.4，否则跳至步骤3.9；步骤3.9、更新由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec，计算方式为：其中，dc1表示与校验节点c2和c3相连的t1和t2类边的度之和，j函数为单调递增函数，j-1函数是j函数的反函数；步骤3.10、根据公式iac＝iac+0.0001计算，再判断如果iac≤1，则跳转至步骤3.2，否则跳转至步骤3.11；步骤3.11、根据上述步骤求得所有的iac和iec的对应关系，利用线性插值法求出当iec为已知iec时与之对应的iac的具体值iac，计算公式为：iac＝interp1(iec,iac,iec,'linear'),其中0≤iec≤1；步骤3.12、根据iec和iac的取值，得到坐标轴中fexit的下界曲线lowerbound(iec，iac)，其中，横轴坐标为由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec，纵轴坐标为t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac。步骤4、根据所述fexit图的上界曲线以及所述fexit图的下界曲线，构造所述五边类型低密度奇偶校验码(fet-ldpc)的fexit的曲线带，从而准确的分析fet-ldpc码的性能。在具体的实施例中，本发明的五边类型低密度奇偶校验码(fet-ldpc)的fexit的曲线构造方法，还可以包括以下步骤：步骤5、根据所述fexit的上界曲线以及所述fexit图的下界曲线，估算fet-ldpc码的门限，如图5，是本发明fet-ldpc码门限值的计算步骤流程图，包括以下步骤：步骤5.1、初始化高斯白噪声方差σn为σ，初始化最小差值max_dif为100，初始化iav＝iec＝0，初始化门限取值范围的最大值σmax；步骤5.2、根据所述步骤2的方法求得iev；根据所述步骤3的方法求得iac；步骤5.3、计算iev与iac的差值dif，计算公式为：dif＝iev-iac；步骤5.4、判断若dif小于等于0，则跳至步骤5.5，否则进一步判断如果dif小于等于max_dif，则将max_idf设置为dif，并将门限值threshold设置为σn，否则如果dif大于max_dif，跳至步骤5.5；步骤5.5、根据公式σn＝σn+0.0001进行计算，并判断若σn≤σmax，则跳至步骤5.2；否则跳至步骤5.6；步骤5.6、获取当前的threshold值作为fet-ldpc码的门限值。上述步骤清晰地说明了本申请fexit图的实现步骤，下面，以一具体的实施例来说明本发明的fexit图的实现方法，并对fet-ldpc码的门限值进行估计。该实施例以一种fet-ldpc码为例，该码字的度分布可以表示为：同时，该码字的参数表可以表示如下：根据上述的多项式构造fet-ldpc码，并在awgn信道和bpsk调制技术下，按照本发明专利提出的fet-ldpc码的输入输出函数进行计算，具体步骤如下：步骤0、设置该fet-ldpc码的码率为r＝0.5，信噪比为snr，则高斯白噪声信道的方差为：步骤1、计算变量节点v1、v2、v3、v4的信道似然信息的方差，根据该码字的参数表可知，v1、v3、v4是通过信道的变量节点，v2是被删余的变量节点(即未通过信道的变量节点)，因此计算公式为：步骤2、计算fet-ldpc码对应的spc检测器的fexit曲线，即fexit图的上界曲线upperbound。步骤2.1、初始化t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav。iav＝0。步骤2.2、根据fet-ldpc码的度分布可知，与信息位变量节点v2相连的t2类边的度dv2＝3。因为信息位变量节点v2被删余，所以其计算t3边上由信息位变量节点v2产生的流向校验节点c1(v2→c1)的似然信息与v2之间的平均互信息iev3：步骤2.3、spc迭代开始，定义it1表示迭代次数，设置it1＝0。步骤2.4、计算t3边上从信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的先验互信息iac3：iac3＝iev3。步骤2.5、计算t4边上从单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的先验互信息iac4：其中，iev4表示t4边上由单校验码v3传给校验节点c1(v3→c1)的似然信息与v3之间的平均互信息。步骤2.6、由fet-ldpc码的度分布可知，与校验节点c1相连的t3类边的度dc3＝3，与校验节点c1相连的t4类边的度dc4＝1。计算t3边上由校验节点c1产生的流向信息位变量节点v2(c1→v2)的似然信息与v2之间的平均互信息为iec3：步骤2.7、计算t3边上从校验节点c1传给信息位变量节点v2(c1→v2)的似然信息与v2之间的先验互信息iav3：iav3＝iec3。步骤2.8、更新t3边上由信息位变量节点v2传给校验节点c1(v2→c1)的似然信息与v2之间的平均互信息为iev3。步骤2.9、令it1＝it1+1，若it1<3，则跳至步骤2.4；否则跳至步骤2.10。步骤2.10、计算t2边上由信息位变量节点v2传给校验节点c2和c3(v2→c2、c3)的似然信息与v2之间的平均互信息iev2：步骤2.11、与信息位变量节点v1相连的t1类边的度dv1＝3。计算t1边上由信息位变量节点v1传给校验节点c2和c3(v1→c2、c3)的似然信息与v1之间的平均互信息iev1：步骤2.12、根据fet-ldpc码的参数表，计算t1类边的边数占信息位变量节点总边数的比例r1：计算t2类边的边数占信息位变量节点总边数的比例r2：计算由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev：iev＝0.6×iev1+0.4×iev2。步骤2.13、iav＝iav+0.0001，如果iav≤1，则跳转至步骤2.2。否则跳转至步骤2.14步骤2.14、根据上述步骤中求得的所有iav和iev取值的对应关系，就可以画出fexit图中的上界曲线upperbound(iav，iev)，图6给出了在snr＝0.5db，码率r＝0.5时，fet-ldpc码对应的spc检测器的fexit曲线，即fexit图的上界曲线upperbound。其中，横轴坐标为t1、t2边上从校验节点c2和c3传给信息位变量节点v1和v2(c2、c3→v1、v2)的似然信息与v1和v2的先验互信息iav，纵轴坐标为由信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1、v2之间的平均互信息iev。如图7，本发明实施例在不同snr下，码率r＝0.5时，该fet-ldpc码对应的fexit图的上界曲线。步骤3、计算fet-ldpc码对应的acc检测器的fexit曲线，即fexit图的下界曲线lowerbound。步骤3.1、初始化t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac为0。步骤3.2、由fet-ldpc码的度分布可知，与校验节点c2和c3相连的t1和t2类边的度之和dc1＝3。计算t5边上由校验节点c2和c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的平均互信息iec5：步骤3.3、acc迭代开始，定义it2表示迭代次数，设置it2＝0。步骤3.4、计算t5边上从校验节点c2、c3传给累积码v4(c2、c3→v4)的似然信息与v4之间的先验互信息iav5：iav5＝iec5步骤3.5、计算t5边上从累积码v4传给校验节点c2、c3(v4→c2、c3)的似然信息与v4之间的平均互信息iev5：其中步骤3.6、计算t5边上从累积码v4传给校验节点c2和c3(v4→c2、c3)的似然信息与v4之间的先验互信息iac5：iac5＝iev5步骤3.7、更新t5边上由校验节点c2和c3传给累积码v4的似然信息与v4之间的平均互信息为iec5：步骤3.8、令it2＝it2+1，若it2<3，则跳至步骤3.4；否则跳至步骤3.9。步骤3.9、更新由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec：步骤3.10、iac＝iac+0.0001，如果iac≤1，则跳转至步骤3.2。否则跳转至步骤3.11。步骤3.11、根据上述步骤求得所有的iac和iec的对应关系，利用线性插值法求出，当iec为已知iec时，与之对应的iac的具体值iac。其中0≤iec≤1。iac＝interp1(iec,iac,iec,'linear')；步骤3.12、根据iec和iac的取值画出fexit的下界曲线lowerbound(iec，iac)，如图8所示，给出了在snr＝0.5db，码率r＝0.5时，fet-ldpc码对应的acc检测器的fexit曲线，即fexit图的下界曲线lowerbound。此时，横轴坐标为由校验节点c2和c3传给信息位变量节点v1和v2的似然信息与v1、v2之间的平均互信息为iec，纵轴坐标为t1、t2边上从信息位变量节点v1和v2传给校验节点c2和c3(v1、v2→c2、c3)的似然信息与v1和v2之间的先验互信息iac。如图9，是本发明实施例在不同snr下，码率r＝0.5时，fet-ldpc码对应的fexit图的下界曲线。步骤4、根据步骤2得到的spc检测器的fexit曲线(即fexit图的上界曲线upperbound(iav，iev))和步骤3得到的acc检测器的fexit曲线(即fexit图的下界曲线lowerbound(iec，iac))，构造fet-ldpc码的fexit的曲线带，如图10所示，是本发明实施例在snr＝0.5db，码率r＝0.5时，由fexit图的上界曲线和下界曲线构成的fet-ldpc码的fexit图，其中两条曲线中间的区域为译码通道。当两条曲线仅在互信息等于1时相交，就说明fet-ldpc码的译码器可以成功收敛。对σn的值进行改变，可以改变中间区域的面积大小。当上界曲线和下界曲线十分接近时，此时的σn也就是fet-ldpc码的门限值。步骤5、根据步骤2和步骤3得到的fexit的上界曲线upperbound(iav，iev)和下界曲线lowerbound(iec，iac)，估算fet-ldpc码的门限。步骤5.1、设置高斯白噪声方差σn＝0.9400，最小差值max_dif＝100，初始值iav＝iec＝0。步骤5.2、根据步骤2和步骤3求出iev和iac。步骤5.3、计算差值dif：dif＝iev-iac。步骤5.4、若dif≤0，则跳至步骤5.5；否则如果dif≤max_dif，则max_dif＝dif，且门限值threshold＝σn；否则如果dif>max_dif，跳至步骤5.5。步骤5.5、σn＝σn+0.0001，若σn≤0.9999，则跳至步骤5.2；否则跳至步骤5.6。步骤5.6、停止迭代计算，将当前的threshold值0.9700作为该fet-ldpc码的门限值。如图11所示，是本发明门实施例在限值σn＝0.9700时，得到的fet-ldpc码的fexit图，此时fexit图的上界曲线upperbound和下界曲线lowerbound最接近。下表给出用fexit图、传统exit图和de算法对该fet-ldpc码的门限值进行估计的结果。fexitde传统exit门限值(db)0.97000.96000.3812从表中可以看到，由于传统exit图没有考虑到多边等特殊结构对fet-ldpc码的影响，因此直接用传统exit图对fet-ldpc码进行估计，其估计的到门限值不准确。传统exit图对门限的估计结果与计算复杂度较高的de算法的估计结果相差很远，而fexit图与de算法估计的结果的误差只有0.01db。因此可以看出发明设计的fexit图的估计出的结果较传统exit图更加准确，同时估计门限的步骤更加简单清晰，更适合用来对fet-ldpc码进行分析。上述具体实施方式只是对本发明的技术方案进行详细解释，本发明并不只仅仅局限于上述实施例，凡是依据本发明原理的任何改进或替换，均应在本发明的保护范围之内。当前第1页12