CN108712233B - 一种基于两边类低密度奇偶校验码的物理层网络编码方法 - Google Patents

一种基于两边类低密度奇偶校验码的物理层网络编码方法 Download PDF

Info

Publication number
CN108712233B
CN108712233B CN201810474633.4A CN201810474633A CN108712233B CN 108712233 B CN108712233 B CN 108712233B CN 201810474633 A CN201810474633 A CN 201810474633A CN 108712233 B CN108712233 B CN 108712233B
Authority
CN
China
Prior art keywords
code
information
tet
point
variable
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810474633.4A
Other languages
English (en)
Other versions
CN108712233A (zh
Inventor
游莹
陈平平
林捷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujian Normal University
Original Assignee
Fujian Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujian Normal University filed Critical Fujian Normal University
Priority to CN201810474633.4A priority Critical patent/CN108712233B/zh
Publication of CN108712233A publication Critical patent/CN108712233A/zh
Application granted granted Critical
Publication of CN108712233B publication Critical patent/CN108712233B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0056—Systems characterized by the type of code used
    • H04L1/0057—Block codes
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105—Decoding
    • H03M13/1111—Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1125—Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using different domains for check node and bit node processing, wherein the different domains include probabilities, likelihood ratios, likelihood differences, log-likelihood ratios or log-likelihood difference pairs
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/27—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes using interleaving techniques
    • H03M13/2703—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes using interleaving techniques the interleaver involving at least two directions
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0041—Arrangements at the transmitter end
    • H04L1/0043—Realisations of complexity reduction techniques, e.g. use of look-up tables
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0052—Realisations of complexity reduction techniques, e.g. pipelining or use of look-up tables
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0076—Distributed coding, e.g. network coding, involving channel coding

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明涉及一种基于两边类低密度奇偶校验码的物理层网络编码方法，包括以下步骤：步骤S1：根据网络层编码模型的特点，将TET‑LPDC码和PNC码结合，编码出TET‑LPDC‑PNC码；步骤S2：结合TET‑LDPC码的度分布，对TET‑LDPC‑PNC码进行译码。本发明将两边类型低密度奇偶校验码与物理层网络编码进行联合设计，根据TET‑LDPC码“两边”的结构特点，从“结点”的角度对TET‑LDPC码进行联合编译码设计，在译码是采用一种串行似然信息译码的方式，该方法使得编译码复杂度都大大减少。

Description

一种基于两边类低密度奇偶校验码的物理层网络编码方法
技术领域
本发明涉及通信的物理层网络编码和信道编译码技术领域，特别是一种基于两边类低密度奇偶校验码的物理层网络编码方法。
背景技术
1948年Shannon在其论文《A mathematical theory of communication》中提出了信道编码理论，从此开创了信道编码理论成为科学家们不懈研究的主题。2000 年，香港中文大学的Ahlswede R、蔡宁等在论文《Network information flow》中提出了网络编码的概念后，引起人们的广泛关注。由于网络编码能够增加网络吞吐量的特性，越来越多的人开始研究网络编码，并取得了许多新成果。
在很长一段时间里，信道编码和网络编码都是分开设计和研究。2006年，技术人员把信道编码和网络编码结合在一起，进一步提出了物理层网络编码 (physical-layernetwork coding，简称PNC)技术。物理层网络编码技术已被应用在一些信道编码上。但是物理层网络编码却还没有与多边类型低密度奇偶校验码(即MET-LDPC码)结合设计。多边类型低密度奇偶校验码作为一种更加一般化的LDPC码，不仅从结构上包容了许多种类的LDPC码，而且在具有很好的纠错性能和较低的编码实现复杂度。其中两边类型低密度奇偶校验码 (Two-Edge Type Low Density Parity Check Codes，简称两边类型LDPC码或 TET－LDPC码)在码字设计阶段就不仅引入了删余变量点，还有效结合了度为 2的边类型，取得更好的编码增益。
发明内容
有鉴于此，本发明的目的是提出一种基于两边类低密度奇偶校验码的物理层网络编码方法，针对TET-LDPC码所具有的特殊结构特点和性能优势，将 TET-LDPC码和物理层网络编码进行联合设计，构造了一种基于两边类型低密度奇偶校验码的物理层网络编码，简称TET-LDPC-PNC码。
本发明采用以下方案实现：一种基于两边类低密度奇偶校验码的物理层网络编码方法，包括以下步骤：
步骤S1：根据网络层编码模型的特点，将TET-LPDC码和PNC码结合，编码出TET-LPDC-PNC码；
步骤S2：结合TET-LDPC码的度分布，对TET-LDPC-PNC码进行译码。
进一步地，所述步骤S1中，网络层编码模型为双向中继的三节点系统，包括第一源节点A、第二源节点B和中继点R，其中第一源节点A和第二源节点 B在中继点R交换彼此的信息；
其中，中继点R接收到的信息通过以下方式获得：
设定XA和XB分别是第一源节点A和第二源节点B待传出的信息比特 (XA∈{0,1}，XB∈{0,1})，XA和XB经过各自的编码器，即TET-LDPC码的编码器，得到编码后的码字CA和CB(CA∈{0,1},CB∈{0,1})，CA和CB经过BPSK调制后得到DA和DB(DA∈{-1,1},DB∈{-1,1})，调制后的信号经过高斯白噪声信道AWGN 后发送给中继点R，则中继点接收到的信息为：
YR＝DA+DB+η
其中η为高斯白噪声信道AWGN的噪声。
进一步地，所述步骤S2中，译码的过程具体包括以下步骤：
步骤S21：定义变量，具体为：
假设信源长度为K，码字总长为N，信噪比为SNR，码率为Rc，迭代次数为iterativeNum，
则采用BPSK调制时，信道的单边功率谱N0为：
Figure BDA0001664103440000031
步骤S22：计算信道信息，具体为：
假设经过中继端点R后得到信道信息是pLch0和pLch2，
对于未删余的变量点，它的信道信息为：
pLch0[i]＝(-4*(1.0-YR))/N0
pLch2[i]＝(-4*(1.0+YR))/N0
0≤i≤N-1
对于已删余的变量点，它的信道信息为：
pLch0＝log(0.5)
pLch2＝log(0.5)
假设边度类型为r(0,1)x(0,2)的变量点的信道信息为pLch_D2_0和pLch_D2_2；假设边度类型为r(1,0)x(q,0)的变量点的信道信息为pLch_D1_0和pLch_D1_2；把经过中继端点R后得到信道信息pLch0和pLch2分别赋给它们；
步骤S23：译码迭代开始，设置迭代次数iterativeNum，设置计数值count＝0；
步骤S24：对边度类型r(0,1)x(0,2)的变量点进行译码，具体为：
对于第i个变量点，与它相连的分别为第i-1个和第i个校验点，用pLef0[i] 和pLef2[i]表示第i-1个校验点传给它的信息，则：
Figure BDA0001664103440000041
Figure BDA0001664103440000042
用pLeb0[i]和pLeb2[i]表示第i个校验点传给它的信息，则：
Figure BDA0001664103440000043
Figure BDA0001664103440000044
步骤S25：对边度类型为x(1,2)的校验点进行译码，具体为：
检验点包括2个输入和1个输出，假设检验点的2个输入为A和B，其中， A的似然信息包含分别为0、2码字时的A0、A2；B的似然信息包含分别为0、 2码字时B0、B2；假设检验点的一个输出为C，通过下列公式计算C的似然信息C0和C2：
Figure BDA0001664103440000045
Figure BDA0001664103440000046
其中
Figure BDA0001664103440000047
两个符号表示在校验点进行的计算；
接着计算D2类的边经过校验点后传到D1类边的似然信息pLex0和pLex2，如下式所示：
Figure BDA0001664103440000051
Figure BDA0001664103440000052
步骤S26：对边度类型为r(1,0)x(q,0)的变量点进行译码：
假设步骤S25中的似然信息pLex0和pLex2的信息经过交织器解交织后的信息为pLex_deinter0和pLex_deinter2，经过第i条D1类边传入交织器前的信息为pLox_deinter0和pLox_deinter2，则根据变量点的边度特点进行如下计算：
Figure BDA0001664103440000053
Figure BDA0001664103440000054
步骤S27：进入交织器：把从D1类边传出来的似然信息pLox_deinter0和 pLox_deinter2送入交织器后，得到传给校验点的似然信息pLox0和pLox2；
步骤S28：判断是否继续迭代：如果计数变量count<iterativeNum，则 count＝count+1，返回步骤S24，继续计算，否则进入步骤S29；
步骤S29：对最后得到的信源点的似然信息行判决，具体为：
假设最后得到的信源点的似然信息是pLLR0和pLLR2，根据以下公式计算信源点的似然信息：
Figure BDA0001664103440000061
Figure BDA0001664103440000062
对pLLR0和pLLR2进行判决，当第i个信源变量点的pLLR0[i]和pLLR2[i] 的值，满足pLLR0[i]≤0和pLLR2[i]≤0时，则判定第i个信源变量点为码字1，其他情况下判定为码字0。
上述步骤清晰地说明了TET-LDPC-PNC码的编译码过程，本发明提供的方法针对TET-LDPC码所具有的特殊结构特点和性能优势，将TET-LDPC码和物理层网络编码进行联合设计，构造了一种基于两边类型低密度奇偶校验码的物理层网络编码，简称TET-LDPC-PNC码。在构造过程中，本发明首先在两个终端节点采用相同信道编码，即TET-LDPC码，设计出性能优异的TET-LDPC-PNC 码；接着在中继端采用一个译码器对接收到的叠加信息进行译码，译码时采用串行的似然信息译码算法，就可以译码得到该信道编码的有效码字。该方法采用相同编码器，可以大大的简化编码复杂度，同时在译码时采用串行的似然信息译码算法，也大大简化译码复杂度，用加运算替代原物理层网络编码中的乘运算，大大减少了译码的计算量。
与现有技术相比，本发明有以下有益效果：
(1)本发明是针对TET-LDPC码的“两边”的结构特点提出的，将其和网络编码相结合，构成性能优异的TET-LDPC-PNC码。
(2)本发明设计的TET-LDPC码具有很好的性能优势，其编译码复杂度都大大减少。
(3)本发明从一个全新的视角来设计物理层网络编译码设计方案，即从“多边”和“结点”角度出发对变量点进行编译码设计。
附图说明
图1是本发明实施例的TET-LDPC码的编码框图。
图2是本发明实施例的TET-LDPC-PNC码的编码框图。
图3是本发明实施例的TET-LDPC-PNC的译码框图。
图4是本发明实施例的边度类型为r(0,1)x(0,2)的变量点的边链接情况。
图5是本发明实施例的边度类型为x(1,2)的校验点的边链接图。
图6是本发明实施例中当q＝2，K分别为512、1000、4000时，TET-LDPC-PNC 码的BER性能曲线。
图7是本发明实施例中当q＝3，K分别为512、1000、4000时，TET-LDPC-PNC 码的BER性能曲线。
图8是本发明实施例中当q＝5，K分别为512、1000、4000时，TET-LDPC-PNC 码的BER性能曲线。
图9是本发明实施例中当K＝512时，一组码率灵活调节的TET-LDPC-PNC码的 BER性能曲线。
具体实施方式
下面结合实施方式与实施例对本发明做进一步说明。
本实施例提供了一种基于两边类低密度奇偶校验码的物理层网络编码方法，包括以下步骤：
步骤S1：根据网络层编码模型的特点，将TET-LPDC码和PNC码结合，编码出TET-LPDC-PNC码；
步骤S2：结合TET-LDPC码的度分布，对TET-LDPC-PNC码进行译码。
在本实施例中，TET-LDPC码与传统的LDPC码不同的是，TET-LDPC码的度分布表示方法是从“结点”的角度来描述，可以通过两个多项式表示，即：
Figure BDA0001664103440000081
μ(x):＝x(1,2) (公式2)
其中第一个多项式对应于变量结点，第二个多项式对应于校验结点。具体参数可以用参数表描述。
本实施例采用的TET-LPDC码的参数表如下：
表1 TET-LDPC码参数
Figure BDA0001664103440000082
其中：
矢量
Figure BDA0001664103440000083
表示边的类型，矢量d的取值个数nT表示该码字中具有的边的种类数；di表示第i类型的边，1≤i≤nT。由上述参数表可知，矢量d 的取值个数nT＝2，因此这是一个TET-LDPC码。其中第一类型d1＝q的边用变量 D1表示，第二类型d2＝2的边用变量D2表示。
矢量
Figure 1
表示该变量结点经过了的信道，在上述参数表中nτ的取值为 2，用b＝(1,0)表示经过该信道的变量点被删余了。
假如用
Figure BDA0001664103440000091
表示变量所对应的接收分布ν(r,x)是变量结点的分布情况，νb,d表示度类型为(b,d)的变量结点个数与码长的比值，μ(x)是校验结点的分布情况，μd表示度类型为d的校验结点个数与码长的比值。
在本实施例中，所述步骤S1中，网络层编码模型为双向中继的三节点系统，包括第一源节点A、第二源节点B和中继点R，其中第一源节点A和第二源节点B在中继点R交换彼此的信息；
其中，根据图1所示TET-LDPC码的编码框图，对TET-LDPC码进行编码设计，同时，结合网络层编码模型的特点，根据图2的编码框图实现 TET-LPDC-PNC码的编码，则中继点R接收到的信息通过以下方式获得：
设定XA和XB分别是第一源节点A和第二源节点B待传出的信息比特 (XA∈{0,1}，XB∈{0,1})，XA和XB经过各自的编码器，即TET-LDPC码的编码器 (如图1所示)，得到编码后的码字CA和CB(CA∈{0,1},CB∈{0,1})，CA和CB经过BPSK调制后得到DA和DB(DA∈{-1,1},DB∈{-1,1})，调制后的信号经过高斯白噪声信道AWGN后发送给中继点R，则中继点接收到的信息为：
YR＝DA+DB+η (公式3)
其中η为高斯白噪声信道AWGN的噪声。
在本实施例中，结合TET-LDPC码的度分布，按照图3给出总的译码框图对TET-LDPC-PNC进行译码，其中在变量点译码和校验点译码时采用的译码方法不同。
因此，所述步骤S2中，译码的过程具体包括以下步骤：
步骤S21：定义变量，具体为：
假设信源长度为K，码字总长为N，信噪比为SNR，码率为Rc，迭代次数为iterativeNum，
则采用BPSK调制时，信道的单边功率谱N0为：
Figure BDA0001664103440000101
步骤S22：计算信道信息，具体为：
假设经过中继端点R后得到信道信息是pLch0和pLch2，
对于未删余的变量点，它的信道信息为：
Figure BDA0001664103440000102
对于已删余的变量点，它的信道信息为：
Figure BDA0001664103440000103
假设边度类型为r(0,1)x(0,2)的变量点的信道信息为pLch_D2_0和pLch_D2_2；假设边度类型为r(1,0)x(q,0)的变量点的信道信息为pLch_D1_0和pLch_D1_2；把经过中继端点R后得到信道信息pLch0和pLch2分别赋给它们；
步骤S23：译码迭代开始，设置迭代次数iterativeNum，设置计数值count＝0；
步骤S24：对边度类型r(0,1)x(0,2)的变量点进行译码，图4显示了边度类型为的变量点的边链接情况，根据图4可以计算出与之相连的校验点通过D2类边传给它的似然信息，具体为：
对于第i个变量点，与它相连的分别为第i-1个和第i个校验点，用pLef0[i] 和pLef2[i]表示第i-1个校验点传给它的信息，则：
Figure BDA0001664103440000111
用pLeb0[i]和pLeb2[i]表示第i个校验点传给它的信息，则：
Figure BDA0001664103440000112
步骤S25：对边度类型为x(1,2)的校验点进行译码，具体为：
根据表1和图5所示，检验点包括2个输入和1个输出，假设检验点的2 个输入为A和B，其中，A的似然信息包含分别为0、2码字时的A0、A2；B 的似然信息包含分别为0、2码字时B0、B2；假设检验点的一个输出为C，通过下列公式计算C的似然信息C0和C2：
Figure BDA0001664103440000113
Figure BDA0001664103440000114
其中
Figure BDA0001664103440000115
两个符号表示在校验点进行的计算；
接着，根据图5边度类型为x(1,2)的校验点的边链接图，计算D2类的边经过校验点后传到D1类边的似然信息pLex0和pLex2，如下式所示：
Figure BDA0001664103440000121
步骤S26：对边度类型为r(1,0)x(q,0)的变量点进行译码：
假设步骤S25中的似然信息pLex0和pLex2的信息经过交织器解交织后的信息为pLex_deinter0和pLex_deinter2，经过第i条D1类边传入交织器前的信息为pLox_deinter0和pLox_deinter2，则根据变量点的边度特点进行如下计算：
Figure BDA0001664103440000122
步骤S27：进入交织器：把从D1类边传出来的似然信息pLox_deinter0和 pLox_deinter2送入交织器后，得到传给校验点的似然信息pLox0和pLox2；
步骤S28：判断是否继续迭代：如果计数变量count<iterativeNum，则 count＝count+1，返回步骤S24，继续计算，否则进入步骤S29；
步骤S29：对最后得到的信源点的似然信息行判决，具体为：
假设最后得到的信源点的似然信息是pLLR0和pLLR2，根据以下公式计算信源点的似然信息：
Figure BDA0001664103440000131
对pLLR0和pLLR2进行判决，当第i个信源变量点的pLLR0[i]和pLLR2[i] 的值，满足pLLR0[i]≤0和pLLR2[i]≤0时，则判定第i个信源变量点为码字1，其他情况下判定为码字0。
根据上述方法，本实施例中设置了具体的TET-LDPC码，通过对多边参数表中的q值进行调整，可以得到不同码率的两边类型低密度奇偶校验码。当所取得q值不一样时，如果从结点角度描述的度分布来描述上述的TET-LDPC码，则多项式分别为：
q＝2时：
Figure BDA0001664103440000132
q＝3时：
Figure BDA0001664103440000133
q＝5时：
Figure BDA0001664103440000134
根据上述的多项式，构造信源长度K分别为512、1000、4000的TET-LDPC 码，再在AWGN信道和BPSK调制技术下，按照本实施例提供的TET-LDPC-PNC 码的编译码方法，搭建TET-LDPC-PNC码的仿真平台。
译码时，设置迭代次数为iterativeNum＝5。按照以下步骤进行实施：
step 1：设置信源长度为K＝512、1000、4000。设置信噪比为SNR＝(0、1、2、3、4……)，计算码率为Rc，根据公式4计算每个SNR下的N0；
step 2：计算信道信息，根据公式5计算出经过中继端点R后得到信道信息是pLch0和pLch2，并对两种边类相连的变量点进行赋值，其中与D2类边相连的变量点的信道信息为pLch_D2_0和pLch_D2_2；与D1类边相连的变量点的道信息为pLch_D1_0和pLch_D1_2；
step 3：设置计数值count＝0；
step4：对边度类型r(0,1)x(0,2)，即与D2类边相连的变量点进行译码，根据公式7和公式8分别计算出pLef0，pLef2，pLeb0，pLeb2；
step5：对边度类型为x(1,2)，即与D1和D2类边相连的校验点进行译码，根据公式10计算出pLex0和pLex2；
step6：对边度类型为r(1,0)x(q,0)，即与D1类边相连的变量点进行译码；将pLex0 和pLex2解交织后得到pLex_deinter0和pLex_deinter2，根据公式11计算经过第i条D1类边传入交织器前的信息为pLox_deinter0和pLox_deinter2；
step7：把从D1类边传出来的似然信息pLox_deinter0和pLox_deinter2送入交织器后，就得到了传给校验点的似然信息pLox0和pLox2；
step8：判断是否继续进行迭代：计数变量count<iterativeNum，则count＝count+1，返回step4，继续计算；直到count＝iterativeNum时，进入步骤 step9；
step9：判决：根据公式12计算最后的到信源点的似然信息是pLLR0和 pLLR2，并进行判决，最后统计BER的值。
经过上述的实施步骤，可以得到一组TET-LDPC-PNC码的性能曲线。
如图6所述，当q＝2，K分别为512、1000、4000时，TET-LDPC-PNC码的BER性能曲线。其中，x轴表示信噪比，y轴表示误比特率。带△的曲线表示q＝2时，信源长度为512的TET-LDPC-PNC码的误比特率仿真曲线。带*的曲线表示q＝2时，信源长度为1000的TET-LDPC-PNC码的误比特率仿真曲线。带○的曲线表示q＝2时，信源长度为4000的TET-LDPC-PNC码的误比特率仿真曲线。
如图7所述，当q＝3，K分别为512、1000、4000时，TET-LDPC-PNC码的BER性能曲线。其中，x轴表示信噪比，y轴表示误比特率。带△的曲线表示q＝3时，信源长度为512的TET-LDPC-PNC码的误比特率仿真曲线。带*的曲线表示q＝3时，信源长度为1000的TET-LDPC-PNC码的误比特率仿真曲线。带○的曲线表示q＝3时，信源长度为4000的TET-LDPC-PNC码的误比特率仿真曲线。
如图8所述，当q＝5，K分别为512、1000、4000时，TET-LDPC-PNC码的BER性能曲线。其中，x轴表示信噪比，y轴表示误比特率。带△的曲线表示q＝5时，信源长度为512的TET-LDPC-PNC码的误比特率仿真曲线。带*的曲线表示q＝5时，信源长度为1000的TET-LDPC-PNC码的误比特率仿真曲线。带○的曲线表示q＝5时，信源长度为4000的TET-LDPC-PNC码的误比特率仿真曲线。
综上所述，可以看到在BER＝10-6的情况下，以上参数设计的TET-LDPC-PNC 码都未出现错误平台。
另外，在AWGN信道下，当对两边类型的参数进行改变时，可以得到一组码率灵活调节的TET-LDPC-PNC码。当K＝512时，一组码率灵活调节的 TET-LDPC-PNC码的BER性能曲线如图9所示，图中，x轴表示信噪比，y轴表示误比特率。带△的曲线表示Rc＝0.5时，信源长度为512的TET-LDPC-PNC 码的误比特率仿真曲线。带*的曲线表示Rc＝0.3时，信源长度为512的TET-LDPC-PNC码的误比特率仿真曲线。带○的曲线表示Rc＝0.2时，信源长度为512的TET-LDPC-PNC码的误比特率仿真曲线。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件和软件结合来完成。以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于两边类低密度奇偶校验码的物理层网络编码方法，其特征在于：包括以下步骤：
步骤S1：根据网络层编码模型的特点，将TET-LDPC码和物理层网络编码结合，编码出TET-LDPC-PNC码；
所述网络层编码模型为双向中继的三节点系统，包括第一源节点A、第二源节点B和中继点R，其中第一源节点A和第二源节点B在中继点R交换彼此的信息；其中，中继点R接收到的信息通过以下方式获得：设定XA和XB分别是第一源节点A和第二源节点B待传出的信息比特(XA∈{0,1}，XB∈{0,1})，XA和XB经过各自的编码器，即TET-LDPC码的编码器，得到编码后的码字CA和CB(CA∈{0,1},CB∈{0,1})，CA和CB经过BPSK调制后得到DA和DB(DA∈{-1,1},DB∈{-1,1})，调制后的信号经过高斯白噪声信道AWGN后发送给中继点R，则中继点接收到的信息为：
YR＝DA+DB+η
其中η为高斯白噪声信道AWGN的噪声；
步骤S2：结合TET-LDPC码的度分布，对TET-LDPC-PNC码进行译码，所述译码的过程具体包括以下步骤：
步骤S21：定义变量，具体为：
假设信源长度为K，码字总长为N，信噪比为SNR，码率为Rc，迭代次数为iterativeNum，
则采用BPSK调制时，信道的单边功率谱N0为：
Figure FDA0002654377860000011
步骤S22：计算信道信息，具体为：
假设经过中继端点R后得到信道信息是pLch0和pLch2，
对于未删余的变量点，它的信道信息为：
pLch0[i]＝(-4*(1.0-YR))/N0
pLch2[i]＝(-4*(1.0+YR))/N0
0≤i≤N-1
对于已删余的变量点，它的信道信息为：
pLch0＝log(0.5)
pLch2＝log(0.5)
假设边度类型为r(0,1)x(0,2)的变量点的信道信息为pLch_D2_0和pLch_D2_2；假设边度类型为r(1,0)x(q,0)的变量点的信道信息为pLch_D1_0和pLch_D1_2；把经过中继端点R后得到信道信息pLch0和pLch2分别赋给它们；
步骤S23：译码迭代开始，设置迭代次数iterativeNum，设置计数值count＝0；
步骤S24：对边度类型r(0,1)x(0,2)的变量点进行译码，具体为：
对于第i个变量点，与它相连的分别为第i-1个和第i个校验点，用pLef0[i]和pLef2[i]表示第i-1个校验点传给它的信息，则：
Figure FDA0002654377860000021
Figure FDA0002654377860000022
用pLeb0[i]和pLeb2[i]表示第i个校验点传给它的信息，则：
Figure FDA0002654377860000031
Figure FDA0002654377860000032
步骤S25：对边度类型为x(1,2)的校验点进行译码，具体为：
检验点包括2个输入和1个输出，假设检验点的2个输入为A和B，其中，A的似然信息包含分别为0、2码字时的A0、A2；B的似然信息包含分别为0、2码字时B0、B2；假设检验点的一个输出为C，通过下列公式计算C的似然信息C0和C2：
Figure FDA0002654377860000033
Figure FDA0002654377860000034
其中
Figure FDA0002654377860000035
两个符号表示在校验点进行的计算；
接着计算D2类的边经过校验点后传到D1类边的似然信息pLex0和pLex2，如下式所示：
Figure FDA0002654377860000036
步骤S26：对边度类型为r(1,0)x(q,0)的变量点进行译码：
假设步骤S25中的似然信息pLex0和pLex2的信息经过交织器解交织后的信息为pLex_deinter0和pLex_deinter2，经过第i条D1类边传入交织器前的信息为pLox_deinter0和pLox_deinter2，则根据变量点的边度特点进行如下计算：
Figure FDA0002654377860000041
步骤S27：进入交织器：把从D1类边传出来的似然信息pLox_deinter0和pLox_deinter2送入交织器后，得到传给校验点的似然信息pLox0和pLox2；
步骤S28：判断是否继续迭代：如果计数变量count<iterativeNum，则count＝count+1，返回步骤S24，继续计算，否则进入步骤S29；
步骤S29：对最后得到的信源点的似然信息行判决，具体为：
假设最后得到的信源点的似然信息是pLLR0和pLLR2，根据以下公式计算信源点的似然信息：
Figure FDA0002654377860000042
Figure FDA0002654377860000043
对pLLR0和pLLR2进行判决，当第i个信源变量点的pLLR0[i]和pLLR2[i]的值，满足pLLR0[i]≤0和pLLR2[i]≤0时，则判定第i个信源变量点为码字1，其他情况下判定为码字0。
CN201810474633.4A 2018-05-17 2018-05-17 一种基于两边类低密度奇偶校验码的物理层网络编码方法 Active CN108712233B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810474633.4A CN108712233B (zh) 2018-05-17 2018-05-17 一种基于两边类低密度奇偶校验码的物理层网络编码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810474633.4A CN108712233B (zh) 2018-05-17 2018-05-17 一种基于两边类低密度奇偶校验码的物理层网络编码方法

Publications (2)

Publication Number Publication Date
CN108712233A CN108712233A (zh) 2018-10-26
CN108712233B true CN108712233B (zh) 2020-11-20

Family

ID=63869147

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810474633.4A Active CN108712233B (zh) 2018-05-17 2018-05-17 一种基于两边类低密度奇偶校验码的物理层网络编码方法

Country Status (1)

Country Link
CN (1) CN108712233B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109639394B (zh) * 2018-11-16 2021-06-01 福建师范大学福清分校 多边类型低密度奇偶校验码的分边类中继译码方法
CN109560820B (zh) * 2018-11-23 2021-07-13 福州大学 一种联合信道编码和物理层网络编码的译码简化方法
CN110611552B (zh) * 2019-10-14 2022-01-14 福建师范大学福清分校 衰落信道下两边类型物理层网络编码的均衡随机删余方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102349257B (zh) * 2009-01-14 2015-02-25 汤姆森特许公司 设计用于多边型低密度奇偶校验编码调制的多路分用器的方法和装置
CN102118231B (zh) * 2011-04-13 2013-04-17 厦门大学 多边类型低密度奇偶校验码的码率自适应编码方法
DE102015116157A1 (de) * 2015-09-24 2017-03-30 Intel IP Corporation Mobilkommunikationssystem unter Verwendung von Subcodiertechniken
CN107733440B (zh) * 2016-08-12 2022-12-02 中兴通讯股份有限公司 多边类型结构化ldpc处理方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Generalized Sum-Product Algorithm for Joint Channel Decoding and Physical-Layer Network Coding in Two-Way Relay Systems";Dirk Wübben;《2010 IEEE Global Telecommunications Conference GLOBECOM 2010》;20101231;第1-5页 *
"多边类型LDPC码的算法研究及其编码器实现";谢伟;《万方学位论文数据库》;20111229;第1-22页 *

Also Published As

Publication number Publication date
CN108712233A (zh) 2018-10-26

Similar Documents

Publication Publication Date Title
Lee New rate-compatible punctured convolutional codes for Viterbi decoding
Alamdar-Yazdi et al. A simplified successive-cancellation decoder for polar codes
CN101039119B (zh) 编码与解码的方法及系统
Barbulescu Iterative decoding of turbo codes and other concatenated codes
CN107248866B (zh) 一种降低极化码译码时延的方法
Garcia-Frias Decoding of low-density parity-check codes over finite-state binary Markov channels
Jeong et al. SC-Fano decoding of polar codes
CN108712233B (zh) 一种基于两边类低密度奇偶校验码的物理层网络编码方法
CN104025459A (zh) 译码处理方法及译码器
CZ407397A3 (cs) Paralelní zřetězené konvoluční kódy s koncovými bity a jejich dekodéry
EP1119109A1 (en) Bandwidth-efficient concatenated trellis-coded modulation decoder and decoding method thereof
CN112564716B (zh) 一种基于剪枝迭代的pc-scma系统联合译码方法
WO2008048723A2 (en) Method and system for improving decoding efficieny in wireless receivers
KR100445899B1 (ko) 터보 코드 인코더 및 그의 부호율 감소 방법
CN108134612A (zh) 纠正同步与替代错误的级联码的迭代译码方法
CN117040545A (zh) 一种基于ldpc码的bicm-id迭代接收方法
CN103457612A (zh) 针对里德所罗门－卷积级联码的迭代软判决译码方法
US7170946B2 (en) System and method for reviving catastrophic codes
Liang et al. Rateless polar-spinal coding scheme with enhanced information unequal error protection
Hanif et al. An efficient serially concatenated polar code with unequal error protection property
RU2301492C2 (ru) Способ передачи голосовых данных в цифровой системе радиосвязи и устройство для его осуществления
CN101150377A (zh) 用于低密度奇偶校验编码的32apsk系统的比特映射方案
US7225392B2 (en) Error correction trellis coding with periodically inserted known symbols
Hebbes et al. Genetic algorithms for turbo codes
Zheng et al. A Universal List Decoding Algorithm with Application to Decoding of Polar Codes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant