CN107370566A - 一种误码条件下的删余Turbo码盲识别方法 - Google Patents
一种误码条件下的删余Turbo码盲识别方法 Download PDFInfo
- Publication number
- CN107370566A CN107370566A CN201710633678.7A CN201710633678A CN107370566A CN 107370566 A CN107370566 A CN 107370566A CN 201710633678 A CN201710633678 A CN 201710633678A CN 107370566 A CN107370566 A CN 107370566A
- Authority
- CN
- China
- Prior art keywords
- mrow
- code
- mtd
- matrix
- msub
- 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.)
- Granted
Links
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/0059—Convolutional codes
-
- 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/0001—Systems modifying transmission characteristics according to link quality, e.g. power backoff
- H04L1/0036—Systems modifying transmission characteristics according to link quality, e.g. power backoff arrangements specific to the receiver
- H04L1/0038—Blind format detection
-
- 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/0061—Error detection codes
-
- 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/0071—Use of interleaving
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Error Detection And Correction (AREA)
Abstract
本发明属于编码参数识别技术领域，公开了一种误码条件下的删余Turbo码盲识别方法，在伽罗华域GF2上对矩阵进行下三角变换，相比于传统矩阵求秩的方法，本发明可以在很大程度上消除误码对码长及交织长度识别的影响。本发明在对分组起点与寄存器长度进行识别时，对哈达玛变换求解出的校验关系进行置信度评价，大大提高了分组起点和寄存器长度的识别率。本发明提出的识别交织映射关系的方法，采用先剔除误码后识别的方法具有很强的抗误码性能，并且针对交织关系的求解方法具有复杂度低，易于实现的特点。
Description
技术领域
本发明属于编码参数识别技术领域，尤其涉及一种误码条件下的删余Turbo码盲识别方法。
背景技术
1993年，Berrou和Glavieux等人提出了Turbo码(PCCC)，其本质是将卷积码和交织器结合在一起，典型的turbo码编码器通常由两个分量编码器、交织器以及删余复用模块组成。其中分量编码器通常是递归系统卷积码(Recursive System Convolutional，RSC)。两个分量编码器分别为RSC1和RSC2，信息序列X0分为两路，第一路直接经过分量编码器RSC1进行卷积码编码形成校验序列X1。第二路首先经过交织器，然后再通过RSC2编码器进行卷积码编码形成X2。为了提高码率，删余Turbo码还具有删余结构，删余矩阵将两个校验序列删余后再与X0经过复接形成最终的Turbo码序列。在非合作通信中，如何在仅获得编码后的比特流的情况下成功的译码得到截获的信息非常困难，需要首先对编码参数进行识别，只有在正确估计出对方的编码参数时才能进行译码，进而得到信源发送的信息。目前针对Turbo码的识别主要局限在无误码和非删余模式下的Turbo码。张永光等曾在其专利《一种删余Turbo码的编码参数盲识别方法》中提出了一种删余Turbo码的识别方法，但其主要采用求秩的方法求取码长和起始点，该方法并不具备抗误码性能，当接收的码流中出现误码时将使该方法失效。
综上所述，现有技术存在的问题是：目前删余Turbo码的识别容错性差，删余Turbo码的交织关系难以重建。
发明内容
针对现有技术存在的问题，本发明提供了一种误码条件下的删余Turbo码盲识别方法。
本发明是这样实现的，一种误码条件下的删余Turbo码盲识别方法，所述误码条件下的删余Turbo码盲识别方法包括以下步骤：
第一步，通过截获的删余Turbo码比特流排列成分析矩阵，利用矩阵分析的方法识别出码长及交织长度，分离出删余Turbo码的信息位比特流和RSC1路校验位比特流构造2/3码率的删余卷积码；
第二步，通过对删余卷积码进行分析，识别出校验向量和码组起点，进一步识别出删余Turbo码的分量编码器的寄存器长度以及生成矩阵；
第三步，通过抽取Turbo码信息位比特流的和RSC2路校验位比特流识别出交织起始点；
第四步，在对交织映射关系进行识别时，求解RSC2路校验位比特流校验比特与信息位之间的校验关系，由于交织器的存在，RSC2路校验位比特流校验比特并不能和信息位比特流信息位形成稳定且唯一的校验关系；
第五步，求解出一组校验向量，对比该组校验向量与得到的校验序列，可得到交织映射位置之间的约束关系，通过深度优先搜索，恢复出交织映射关系。
进一步，所述误码条件下的删余Turbo码盲识别方法具体包括以下步骤：
步骤一，对删余Turbo码数据流构造分析矩阵，排列成p行q列的分析矩阵Cp×q，根据校验位与信息位的约束关系，采用矩阵变换的方法求得码长n和交织长度L；
步骤二，提取删余Turbo码数据流的信息位以及RSC1分量编码器产生的校验位，构造删余卷积码，根据构造出的删余卷积码识别出校验序列h以及校验序列的置信度θ，根据置信度分析出删余Turbo码的码组起点shift、校验序列h、生成多项式序列G以及分量编码器的寄存器长度m；
步骤三，提取Turbo码的信息位与RSC2分量编码器产生的校验位，构造先交织再编码的删余卷积码，根据约束关系遍历起始点求得Turbo码的交织起始点head，并从head开始截取Turbo码，去除无效比特；
步骤四，构造码字矩阵，根据线性约束关系剔除矩阵中含误码的行，得到无误码的分析矩阵；
步骤五，根据得到的分析矩阵，求解RSC2路校验比特与信息位之间的校验关系，求解出一组校验向量，对比校验向量与得到的校验序列h，得到交织映射位置之间的约束关系，通过深度优先搜索，恢复出交织映射关系。
进一步，所述码长求取方法包括：
1)构造分析矩阵Cp×q，取q＝10，p＝q+50；
2)对分析矩阵Cp×q进行下三角变换，得到变换矩阵Mq，取Mq的第q+1至p行构成矩阵Hq+1→p；
3)记其中Bk表示Hq+1→p第k列的列重，k＝1～q；
4)记集合其中card{·}表示集合的势；
5)若集合不为空，将对应的列值q及对应的秩亏d保存到集合colum_set中；
6)取q＝q+1，若q＜＝100返回1)，否则执行7)；
7)对colum_set所有相邻值求最大公约数并记录最大公约数出现的次数，取出现次数最多的最大公约数即为码长n；
8)对colum_set的秩亏求二阶导数，提取秩变化率较大的位置，求取秩亏较大的q值的最大公约数，取出现次数最多的最大公约数，为交织块长Lblock，交织长度L＝Lblock/2。
进一步，所述码组起点shift、校验序列h、生成多项式序列G以及分量编码器的寄存器长度m的识别方法包括：
1)初始化m＝2，shift＝1，
2)从shift处开始截取删余Turbo码比特流，得到Turbo_cut，提取Turbo_cut的信息比特以及RSC1分量编码器的校验比特，得到2/3码率的删余卷积码比特流Conv1；
3)对Conv1进行哈达玛变换，得到h和置信度θ。令shift＝shift+1，若shift＞4，令m＝m+1，若m＞8，执行4)，否则返回2)；
4)选取置信度θ最高的h，为所识别校验序列，shift和m为正确的码组起点和寄存器长度；
5)识别出的2/3码率卷积码的校验序列h表示为：
h＝[hm hm-1 hm-2 h0]＝[hm，1 hm，2 hm，3 hm-1，1 hm-1，2 hm-1，3…h0，1 h0，2 hm，3]；
其相应的校验多项式矩阵表示为：HD＝[h1(D)h2(D)h3(D)]，其中：
6)设未删余Turbo码的信息支路与RSC1分量编码器支路的合成序的卷积码的基本生成矩阵为：
g∞＝[g0 g1…gm]＝[g0，1 g0，2 g0，3 g1，1 g1，2 g1，3…gm，1 gm，2 gm，3]；
其相应的生成多项式矩阵表示为：GD＝[g1(D)g2(D)]，其中则其等效的2/4码率的卷积码的生成多项式矩阵表示为：
其中其中k＝0，1；
7)遍历删余模式Punc，根据删余模式中删余的位置，删除Gp_tmp中相应位置的列，形成删余生成多项式矩阵Gp；
8)根据Gp·HD T＝0解方程求出未知系数gi，j，完成生成多项式G的识别，若返回7)直至G不为空，完成生成多项式的识别。
进一步，所述交织起始点识别方法具体包括：
1)提取Turbo码的信息位与RSC2分量编码器产生的校验位；
2)构造先交织再编码的删余卷积码，依次截取Turbo码的前head＝1～L*3/2个比特，并构造矩阵Cp*q(q＝L*3/4，p＝q+20)；
3)通过下三角变换求解出每个分析矩阵的秩亏；
4)遍历处理接收到的所有比特块，秩亏的和最大的起始点head即为正确的交织起始点。
进一步，所述步骤四具体包括：
1)构造一个q行L×4列的码字矩阵，通过列变换，列化简为0；
2)筛选出列重低于10的列，根据1的位置，找到对应的含误码的行，并剔除该行。
进一步，所述步骤五具体包括：
1)从剔除误码后的Turbo码中提取信息位比特，构造一个L×4行L×2列的信息位矩阵；
2)提取RSC2路校验位构造一个L×4行，L列的校验位矩阵；
3)循环求解校验位矩阵中的相邻width列与信息位矩阵的校验关系，其中width＝length(h)/3，得到多个校验向量，校验向量中1的位置与步骤二中校验向量中1的位置具有对应关系，通过对应关系求出每个交织位置的可能的交织序号或者交织序号的集合；
4)通过深度优先搜索，恢复出交织映射关系。
本发明的目的在于提供一种使用所述误码条件下的删余Turbo码盲识别方法。
本发明的优点及积极效果为：本发明在伽罗华域GF2上对矩阵进行下三角变换，相比于传统矩阵求秩的方法，本发明可以在很大程度上消除误码对交织长度识别的影响。本发明在对RSC编码器结构进行识别时，对哈达玛变换求解出的校验关系进行置信度评价，根据置信度对校验向量与分组起点进行筛选，相比传统直接根据校验向量结构进行判断的方法在抗误码性能上具有很大提升。本发明提供的方法识别率明显优于传统方法。针对删余Turbo码交织映射关系难以重建，本发明提供一种识别交织映射关系的方法，该方法采用先剔除误码后识别的方式，具有很强的容错能力，并且针对交织关系的恢复具有复杂度低，易于实现的特点。
附图说明
图1是本发明中的方法与秩准则法识别率对比曲线
图2是本发明中的方法与传统方法在识别编码器结构时的识别率对比曲线
图3是本发明中的方法在识别删余Turbo码时的总识别率曲线
图4是本发明实施例提供的误码条件下的删余Turbo码盲识别方法流程图。
图5是本发明实施例提供的误码条件下的删余Turbo码盲识别方法的实现流程图。
图6是本发明实施例提供的Turbo码码长以及交织长度识别的子流程图。
图7是本发明实施例提供的Turbo码的码组起点、校验序列、生成多项式序列以及分量编码器的寄存器长度识别的子流程图。
图8是本发明实施例提供的恢复Turbo码交织映射关系子流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明在伽罗华域GF2上对矩阵进行下三角变换，相比于传统矩阵求秩的方法，可以在很大程度上消除误码对交织长度识别的影响，交织长度的识别结果如图1所示。本发明在对RSC编码器结构进行识别时，对哈达玛变换求解出的校验关系进行置信度评价，根据置信度对校验向量与分组起点进行筛选，相比传统直接根据校验向量结构进行判断的方法在抗误码性能上具有很大提升。如图2所示，本发明中的方法识别率明显优于传统方法。针对删余Turbo码交织映射关系难以重建，本发明提出了一种新型的识别交织映射关系的方法，该方法采用先剔除误码后识别的方式，具有很强的容错能力，并且针对交织关系的恢复具有复杂度低，易于实现的特点。图3给出了删余Turbo码误码率与识别率的关系图。
下面结合附图对本发明的应用原理作详细的描述。
如图4所示，本发明实施例提供的误码条件下的删余Turbo码盲识别方法包括以下步骤：
S101：通过截获的删余Turbo码比特流排列成分析矩阵，利用矩阵分析的方法识别出码长及交织长度，分离出删余Turbo码的信息位比特流和RSC1路校验位比特流构造2/3码率的删余卷积码；
S102：通过对删余卷积码进行分析，识别出校验向量和码组起点，可识别出删余Turbo码的分量编码器的寄存器长度以及生成矩阵；
S103：通过抽取Turbo码信息位比特流的和RSC2路校验位比特流识别出交织起始点；
S104：在对交织映射关系进行识别时，求解RSC2路校验位比特流校验比特与信息位之间的校验关系，由于交织器的存在，RSC2路校验位比特流校验比特并不能和信息位比特流信息位形成稳定且唯一的校验关系；
S105：求解出一组校验向量，对比该组校验向量与步上面得到的校验序列，可得到交织映射位置之间的约束关系，通过深度优先搜索，可恢复出交织映射关系。
下面结合附图对本发明的应用原理作进一步的描述。
本发明提出一种误码条件下的删余Turbo码的盲识别方法，在未知任何先验知识的情况下，准确识别出码长n，寄存器长度，生成多项式矩阵，交织长度，交织映射关系等编码参数，主要包括以下步骤：
步骤1，对删余Turbo码数据流构造分析矩阵，将其排列成p行q列的分析矩阵Cp×q，根据校验位与信息位的约束关系，采用矩阵变换的方法求得码长n和交织长度L。
步骤2，提取删余Turbo码数据流的信息位以及RSC1分量编码器产生的校验位，构造删余卷积码，根据构造出的删余卷积码识别出校验序列h以及校验序列的置信度θ，根据置信度可进一步分析出删余Turbo码的码组起点shift、校验序列h、生成多项式序列G以及分量编码器的寄存器长度m。
步骤3，提取Turbo码的信息位与RSC2分量编码器产生的校验位，构造先交织再编码的删余卷积码，根据约束关系遍历起始点即可求得Turbo码的交织起始点head，并从head开始截取Turbo码，去除无效比特。
步骤4，构造码字矩阵，根据线性约束关系剔除矩阵中含误码的行，得到无误码的分析矩阵。
步骤5，根据步骤4得到的分析矩阵，求解RSC2路校验比特与信息位之间的校验关系，由于交织器的存在，RSC2路的校验比特并不能和信息位形成稳定且唯一的校验关系，最终可求解出一组校验向量，对比该组校验向量与步骤2得到的校验序列h，即可得到交织映射位置之间的约束关系，通过深度优先搜索，即可恢复出交织映射关系。
参考图5，本发明实施例的具体实现步骤如下：
步骤1：识别码长，交织长度。参考图6，本步骤的具体实现流程如下：
1.1)初始化Turbo码参数，码长n＝0，码字起始点shift＝1，交织起始点head＝1，删余模式交织映射关系最大交织长度Lmax
1.2)构造分析矩阵Cp×q，取q＝10，p＝q+50，对分析矩阵Cp×q进行下三角变换，得到变换矩阵Mq，取Mq的第q+1至p行构成矩阵Hq+1→p，记其中Bk表示Hq+1→p第k列的列重k＝1～q，记集合其中card{·}表示集合的势。若集合不为空，将对应的列值q与秩亏d保存到集合colum_set中。取q＝10～100，得到集合colum_set；
1.3)对colum_set所有相邻q值求最大公约数并记录最大公约数出现的次数，取出现次数最多的最大公约数即为码长n。
1.4)对colum_set中的秩亏求二阶导数，并提取秩变化率较大的位置，求取秩亏较大的q值的最大公约数，取出现次数最多的最大公约数即为交织块长Lblock，此时交织长度L＝Lblock/2。
步骤2，识别码组起点shift，校验序列h，生成多项式G，寄存器长度m。参考图7，本步骤的具体实现如下：
2.1)遍历m＝2～mmax，shift＝1～nmax从shift处开始截取删余Turbo码比特流，得到Turbo_cut。提取Turbo_cut的信息比特以及RSC1分量编码器的校验比特，得到2/3码率的删余卷积码比特流Conv1。
2.2)对Conv1进行哈达玛变换，得到h和h的置信度θ，选取遍历过程中置信度最高的h即为校验序列，此时的shift和m即为正确的码组起点和寄存器长度。识别出的2/3码率卷积码的校验序列可表示为：h_estimate＝[hm hm-1 hm-2 h0]＝[hm，1 hm，2 hm，3 hm-1，1 hm-1，2hm-1，3…h0，1 h0，2 hm，3]，其相应的校验多项式矩阵可表示为：HD＝[h1(D)h2(D)h3(D)]，其中：
2.3)设未删余Turbo码的信息支路与RSC1分量编码器支路的合成序的卷积码的基本生成矩阵为：
g∞＝[g0 g1…gm]＝[g0，1 g0，2 g0，3 g1，1 g1，2 g1，3…gm，1 gm，2 gm，3]
其相应的生成多项式矩阵可表示为：GD＝[g1(D)g2(D)]，其中则其等效的2/4码率的卷积码的生成多项式矩阵可表示为：
其中其中k＝0，1。遍历删余模式Punc，根据删余模式中删余的位置，删除Gp_tmp中相应位置的列，形成删余生成多项式矩阵Gp。根据Gp·HD T＝0解方程求出未知系数gi，j，即完成生成多项式G的识别，此时识别出G以及删余模式Punc。
步骤3，识别交织起点。具体流程如下：
3.1)提取Turbo码的信息位与RSC2分量编码器产生的校验位，构造先交织再编码的删余卷积码，依次截取Turbo码的前head＝1～L*3/2个比特，并构造矩阵Cp*q(q＝L*3/4，p＝q+20)，通过下三角变换求解出每个分析矩阵的秩亏。
3.2)遍历处理接收到的所有比特块，秩亏的和最大的起始点head即为正确的交织起始点。
步骤4，剔除误码
4.1)构造一个q行L*4列(其中L为交织长度，q可取列数的4倍)的码字矩阵。
4.2)通过列变换，化简矩阵，找到列重w＜10的列，然后根据该列中1出现的行，该行即为含误码的行，剔除该行得到矩阵Turbo_matrix。
步骤5，恢复交织映射关系
5.1)首先提取Turbo_matrix中的奇数列，与前L*4行，构造一个L*4行L*2列的信息位矩阵；
5.2)提取矩阵Turbo_matrix中RSC2路校验位构造一个L*4行，L列的校验位矩阵；
5.3)循环求解校验位矩阵中的相邻width列与信息位矩阵的校验关系，其中width＝length(h)/3，从而得到多个校验向量，这些校验向量中1的位置与步骤2中校验向量中1的位置具有对应关系，通过对应关系可以求出每个交织位置的可能的交织序号或者交织序号的集合；
5.4)通过深度优先搜索，即可恢复出交织映射关系。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种误码条件下的删余Turbo码盲识别方法，其特征在于，所述误码条件下的删余Turbo码盲识别方法包括以下步骤：
第一步，通过截获的删余Turbo码比特流排列成分析矩阵，利用矩阵分析的方法识别出码长及交织长度，分离出删余Turbo码的信息位比特流和RSC1路校验位比特流构造2/3码率的删余卷积码；
第二步，通过对删余卷积码进行分析，识别出校验向量和码组起点，进一步识别出删余Turbo码的分量编码器的寄存器长度以及生成矩阵；
第三步，通过抽取Turbo码信息位比特流的和RSC2路校验位比特流识别出交织起始点；
第四步，在对交织映射关系进行识别时，求解RSC2路校验位比特流校验比特与信息位之间的校验关系，由于交织器的存在，RSC2路校验位比特流校验比特并不能和信息位比特流信息位形成稳定且唯一的校验关系；
第五步，求解出一组校验向量，对比该组校验向量与得到的校验序列，可得到交织映射位置之间的约束关系，通过深度优先搜索，恢复出交织映射关系。
2.如权利要求1所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述误码条件下的删余Turbo码盲识别方法具体包括以下步骤：
步骤一，对删余Turbo码数据流构造分析矩阵，排列成p行q列的分析矩阵Cp×q，根据校验位与信息位的约束关系，采用矩阵变换的方法求得码长n和交织长度L；
步骤二，提取删余Turbo码数据流的信息位以及RSC1分量编码器产生的校验位，构造删余卷积码，根据构造出的删余卷积码识别出校验序列h以及校验序列的置信度θ，根据置信度分析出删余Turbo码的码组起点shift、校验序列h、生成多项式序列G以及分量编码器的寄存器长度m；
步骤三，提取Turbo码的信息位与RSC2分量编码器产生的校验位，构造先交织再编码的删余卷积码，根据约束关系遍历起始点求得Turbo码的交织起始点head，并从head开始截取Turbo码，去除无效比特；
步骤四，构造码字矩阵，根据线性约束关系剔除矩阵中含误码的行，得到无误码的分析矩阵；
步骤五，根据得到的分析矩阵，求解RSC2路校验比特与信息位之间的校验关系，求解出一组校验向量，对比校验向量与得到的校验序列h，得到交织映射位置之间的约束关系，通过深度优先搜索，恢复出交织映射关系。
3.如权利要求2所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述码长求取方法包括：
1)构造分析矩阵Cp×q，取q＝10，p＝q+50；
2)对分析矩阵Cp×q进行下三角变换，得到变换矩阵Mq，取Mq的第q+1至p行构成矩阵Hq+1→p；
3)记其中Bk表示Hq+1→p第k列的列重，k＝1～q；
4)记集合其中card{·}表示集合的势；
5)若集合不为空，将对应的列值q及对应的秩亏d保存到集合colum_set中；
6)取q＝q+1，若q＜＝100返回1)，否则执行7)；
7)对colum_set所有相邻值求最大公约数并记录最大公约数出现的次数，取出现次数最多的最大公约数即为码长n；
8)对colum_set中的秩亏求二阶导数，并提取秩变化率较大的位置，求取秩亏较大的q值的最大公约数，最后取出现次数最多的最大公约数为交织块长Lblock，则交织长度L＝Lblock/2。
4.如权利要求2所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述码组起点shift、校验序列h、生成多项式序列G以及分量编码器的寄存器长度m的识别方法包括：
1)初始化m＝2，shift＝1，
2)从shift处开始截取删余Turbo码比特流，得到Turbo_cut，提取Turbo_cut的信息比特以及RSC1分量编码器的校验比特，得到2/3码率的删余卷积码比特流Conv1；
3)对Conv1进行哈达玛变换，得到h和置信度θ；令shift＝shift+1，若shift＞4，令m＝m+1，若m＞8，执行4)，否则返回2)；
4)选取置信度θ最高的h，为所识别校验序列，其中shift和m为正确的码组起点和寄存器长度；
5)识别出的2/3码率卷积码的校验序列h表示为：
h＝[hmhm-1hm-2h0]＝[hm,1hm,2hm,3hm-1,1hm-1,2hm-1,3…h0,1h0,2hm,3]；
其相应的校验多项式矩阵表示为：HD＝[h1(D)h2(D)h3(D)]，其中：
<mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>h</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mi>D</mi>
<mi>j</mi>
</msup>
<mo>;</mo>
</mrow>
6)设未删余Turbo码的信息支路与RSC1分量编码器支路的合成序的卷积码的基本生成矩阵为：
g∞＝[g0 g1 … gm]＝[g0,1 g0,2 g0,3 g1,1 g1,2 g1,3 … gm,1 gm,2 gm,3]；
其相应的生成多项式矩阵表示为：GD＝[g1(D) g2(D)]，其中则其等效的2/4码率的卷积码的生成多项式矩阵表示为：
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mo>_</mo>
<mi>t</mi>
<mi>m</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mo>&CenterDot;</mo>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>&CenterDot;</mo>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>&CenterDot;</mo>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中其中k＝0,1；
7)遍历删余模式Punc，根据删余模式中删余的位置，删除Gp_tmp中相应位置的列，形成删余生成多项式矩阵Gp；
8)根据Gp·HD T＝0解方程求出未知系数gi,j，完成生成多项式G的识别，若返回7)直至G不为空，完成生成多项式的识别。
5.如权利要求2所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述交织起始点识别方法具体包括：
1)提取Turbo码的信息位与RSC2分量编码器产生的校验位；
2)构造先交织再编码的删余卷积码，依次截取Turbo码的前head＝1～L*3/2个比特，并构造矩阵Cp*q(q＝L*3/4,p＝q+20)；
3)通过下三角变换求解出每个分析矩阵的秩亏；
4)遍历处理接收到的所有比特块，秩亏的和最大的起始点head即为正确的交织起始点。
6.如权利要求2所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述步骤四具体包括：
1)构造一个q行L×4列的码字矩阵，通过列变换，列化简为0；
2)筛选出列重低于10的列，根据1的位置，找到对应的含误码的行，并剔除该行。
7.如权利要求2所述的误码条件下的删余Turbo码盲识别方法，其特征在于，所述步骤五具体包括：
1)从剔除误码后的Turbo码中提取信息位比特，构造一个L×4行L×2列的信息位矩阵；
2)提取RSC2路校验位构造一个L×4行，L列的校验位矩阵；
3)循环求解校验位矩阵中的相邻width列与信息位矩阵的校验关系，其中width＝length(h)/3，得到多个校验向量，校验向量中1的位置与步骤二中校验向量中1的位置具有对应关系，通过对应关系求出每个交织位置的可能的交织序号或者交织序号的集合；
4)通过深度优先搜索，恢复出交织映射关系。
8.一种使用权利要求1～7任意一项所述误码条件下的删余Turbo码盲识别方法的编码器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710633678.7A CN107370566B (zh) | 2017-07-28 | 2017-07-28 | 一种误码条件下的删余Turbo码盲识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710633678.7A CN107370566B (zh) | 2017-07-28 | 2017-07-28 | 一种误码条件下的删余Turbo码盲识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107370566A true CN107370566A (zh) | 2017-11-21 |
| CN107370566B CN107370566B (zh) | 2020-07-14 |
Family
ID=60308243
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710633678.7A Active CN107370566B (zh) | 2017-07-28 | 2017-07-28 | 一种误码条件下的删余Turbo码盲识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107370566B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110535478A (zh) * | 2019-09-27 | 2019-12-03 | 电子科技大学 | 一种DVB-RCS2协议中双输入类Turbo码闭集识别方法 |
| CN110690907A (zh) * | 2019-09-27 | 2020-01-14 | 电子科技大学 | 一种已知支路信息turbo码删除模式估计方法 |
| CN112165338A (zh) * | 2020-09-30 | 2021-01-01 | 电子科技大学 | 一种卷积码随机交织序列交织关系的估计方法 |
| CN112165337A (zh) * | 2020-09-30 | 2021-01-01 | 电子科技大学 | 基于线性约束的卷积码随机交织序列交织关系估计方法 |
| CN113078914A (zh) * | 2021-03-26 | 2021-07-06 | 重庆邮电大学 | 对单个非完整Turbo码分量编码器参数盲识别方法 |
| CN116566404A (zh) * | 2023-07-11 | 2023-08-08 | 北京谷数科技股份有限公司 | 删余Turbo码交织映射关系的确定方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103312457A (zh) * | 2013-05-09 | 2013-09-18 | 西安电子科技大学 | 卷积码编码参数全盲识别方法 |
| CN104467875A (zh) * | 2014-12-09 | 2015-03-25 | 山东大学 | 一种rs码与删余卷积码级联码的参数盲识别方法 |
| CN104683072A (zh) * | 2015-03-26 | 2015-06-03 | 山东大学 | 一种删余turbo码分量编码器的参数盲识别方法 |
| CN105391455A (zh) * | 2015-10-31 | 2016-03-09 | 电子科技大学 | 一种归零Turbo码起点及深度盲识别方法 |
| US20160127032A1 (en) * | 2014-11-05 | 2016-05-05 | Electronics And Telecommunications Research Institute | Low-rate data transmission in lte based satellite ratio interface |
-
2017
- 2017-07-28 CN CN201710633678.7A patent/CN107370566B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103312457A (zh) * | 2013-05-09 | 2013-09-18 | 西安电子科技大学 | 卷积码编码参数全盲识别方法 |
| US20160127032A1 (en) * | 2014-11-05 | 2016-05-05 | Electronics And Telecommunications Research Institute | Low-rate data transmission in lte based satellite ratio interface |
| CN104467875A (zh) * | 2014-12-09 | 2015-03-25 | 山东大学 | 一种rs码与删余卷积码级联码的参数盲识别方法 |
| CN104683072A (zh) * | 2015-03-26 | 2015-06-03 | 山东大学 | 一种删余turbo码分量编码器的参数盲识别方法 |
| CN105391455A (zh) * | 2015-10-31 | 2016-03-09 | 电子科技大学 | 一种归零Turbo码起点及深度盲识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张永光: ""一种Turbo码编码参数的盲识别方法"", 《西安电子科技大学学报(自然科学版)》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110535478A (zh) * | 2019-09-27 | 2019-12-03 | 电子科技大学 | 一种DVB-RCS2协议中双输入类Turbo码闭集识别方法 |
| CN110690907A (zh) * | 2019-09-27 | 2020-01-14 | 电子科技大学 | 一种已知支路信息turbo码删除模式估计方法 |
| CN110535478B (zh) * | 2019-09-27 | 2023-02-07 | 电子科技大学 | 一种DVB-RCS2协议中双输入类Turbo码闭集识别方法 |
| CN110690907B (zh) * | 2019-09-27 | 2023-04-25 | 电子科技大学 | 一种已知支路信息turbo码删除模式估计方法 |
| CN112165338A (zh) * | 2020-09-30 | 2021-01-01 | 电子科技大学 | 一种卷积码随机交织序列交织关系的估计方法 |
| CN112165337A (zh) * | 2020-09-30 | 2021-01-01 | 电子科技大学 | 基于线性约束的卷积码随机交织序列交织关系估计方法 |
| CN112165337B (zh) * | 2020-09-30 | 2024-01-26 | 电子科技大学 | 基于线性约束的卷积码随机交织序列交织关系估计方法 |
| CN113078914A (zh) * | 2021-03-26 | 2021-07-06 | 重庆邮电大学 | 对单个非完整Turbo码分量编码器参数盲识别方法 |
| CN113078914B (zh) * | 2021-03-26 | 2023-08-08 | 重庆邮电大学 | 对单个非完整Turbo码分量编码器参数盲识别方法 |
| CN116566404A (zh) * | 2023-07-11 | 2023-08-08 | 北京谷数科技股份有限公司 | 删余Turbo码交织映射关系的确定方法及装置 |
| CN116566404B (zh) * | 2023-07-11 | 2023-09-19 | 北京谷数科技股份有限公司 | 删余Turbo码交织映射关系的确定方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107370566B (zh) | 2020-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107370566A (zh) | 一种误码条件下的删余Turbo码盲识别方法 | |
| CN105720992A (zh) | 一种极化码的简化译码方法 | |
| CN109936377A (zh) | 一种分段crc辅助的极化码编码和译码方法 | |
| CN107147400A (zh) | 一种ira码开集盲识别方法 | |
| CN106059712A (zh) | 一种高误码的任意码率卷积码编码参数盲识别方法 | |
| CN104242957B (zh) | 译码处理方法及译码器 | |
| CN102315911B (zh) | 一种低密度奇偶校验码编码方法及装置 | |
| CN105978577A (zh) | 一种基于比特翻转的串行列表译码算法 | |
| CN102244554B (zh) | 一种删余Turbo码编码参数的盲识别方法 | |
| CN107453807B (zh) | 一种大气光通信信道模型的极化方法、装置及电子设备 | |
| CN103236900B (zh) | 一种串行级联Turbo码交织器参数盲估计方法 | |
| CN101373978B (zh) | 一种Turbo码译码方法以及装置 | |
| CN109981224A (zh) | 一种深空通信信道编译码系统及其方法 | |
| CN111726125B (zh) | 一种基于模板匹配的含错级联码识别方法 | |
| CN102835062B (zh) | 确定码距的方法及计算装置 | |
| CN117097441B (zh) | 基于数据分析的载波通信系统传输效率优化方法 | |
| CN114595427A (zh) | 基于非自回归模型的含缺失值车辆轨迹序列填补修复方法 | |
| Anderson et al. | An optimal circular Viterbi decoder for the bounded distance criterion | |
| CN103856218B (zh) | 译码处理方法及译码器 | |
| CN116232340B (zh) | 一种卫星通信信号的稀疏校验矩阵参数估计方法及装置 | |
| CN102244553B (zh) | 一种非归零Turbo码编码参数的盲识别方法 | |
| CN109361401A (zh) | 一种用于随钻测量传输系统的极化信道编码方法 | |
| CN102244521A (zh) | 一种归零Turbo码编码参数的盲识别方法 | |
| CN108833057B (zh) | 一种基于3gpp 5g-nr极化码的分布式crc处理方法 | |
| CN105610550B (zh) | 一种用于电力线载波通信的Viterbi译码方法 |
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 |