一种bch码做外码的级联码识别交织长度的盲识别方法
【专利摘要】一种BCH码做外码的级联码识别交织长度的盲识别方法，属信道编码盲识别【技术领域】，通过计算机进行数据读入、分析及计算处理，首先建立一个矩阵，通过线性矩阵分析找出对应的最小秩时所存在的最大列数z，再以z或z的整倍数做为列数建立新的矩阵，求出不同列数的秩，若满足所有的以秩为纵轴，列数为横轴的点斜率均相等，即可求出交织长度。此方法可以解决BCH作为外码时交织长度的识别问题，同时也提升了识别速度。
【专利说明】—种BCH码做外码的级联码识别交织长度的盲识别方法
【技术领域】
[0001]本发明涉及数字通信系统中的一种BCH码做外码的级联码识别交织长度的盲识别方法，属于信道编码盲识别【技术领域】。
【背景技术】
[0002]级联编码是一种通过短分量码构造强纠错能力的长码的有效技术。使用非二进制码作为外码、二进制码作为内码的级联码已广泛应用于通信和数据存储系统中，大多数应用中，使用RS码或BCH码作为外码，卷积码作为内码。
[0003]由于级联码较好的编码增益和较低的实现复杂度在各种通信场合中得到了广泛的应用。张永光、楼才义著的《信道编码及其识别分析》中给出定理:对经(n，k)线性分组码+分组交织编码后的数据序列所构成的P X q矩阵，若q为L或L的整数倍，则单位化后左上角单位阵的维数相等，且此时矩阵的秩不等于列数q。此定理对于RS码做外码的级联码交织长度的盲识别适用，然而有时却不适用于BCH码做外码的级联码交织长度的识别。且因为其对每一个不同的列数均要建立矩阵并进行化简，导致识别速度降低。因此，针对上述问题本文提出了一种针对BCH码做外码的级联码交织长度盲识别的算法。
【发明内容】
[0004]为了克服现有技术存在的缺陷和不足，本发明提供了一种BCH码做外码的级联码识别交织长度的盲识别方法。
[0005]本发明的技术方案如下:
[0006]一种BCH码做外码的级联码识别交织长度的盲识别方法，通过计算机进行数据读入、分析及计算处理，该方法步骤如下:
[0007](I)由计算机从待识别数据读入部分待识别的数据，设置一个表示列数循环时其最大值的变量q_max，其值为100 ；
[0008](2)设列数q的初值为10，此值应小于可能出现的最小的交织长度，其中交织长度等于BCH码长乘以交织宽度，BCH码阶数范围为3?8，交织宽度最小为2，因此最小交织长度为(23-1) X2 = 14 ；
[0009](3)建立一个pXq大小的矩阵，其中P是行数，q是列数，且p>q ;
[0010](4)对所建立的矩阵从左到右按列化简，若对角线上元素为1，则将此行依次与其下方每一行进行模二加运算；如果对角线元素为0，则寻找该列对角线下方的非零元素所在行，将非零元素所在行与当前行互换，再执行上述化简；如果对角线下方元素全为0，则不再化简；
[0011](5)计算化简后矩阵的秩，若秩不等于列数时则记录此时矩阵的列数和秩，否则不记录;
[0012](6)将列数q加I,当q>q_max时,转入步骤(7),否则转入步骤(3)；
[0013](7)比较保存的所有秩，求出最小的秩r_min ;在此基础上，统计所有满足该最小秩的列数，求出该最小秩所对应的最大的列数Z ;
[0014](8)设循环变量i的初值为2 ；
[0015](9)使用待识别数据重新建立一个pXq大小的矩阵，其中P为行数，q为列数(P和q跟上面的是一样的变量，只是赋值不同)，且q的值为⑶中i的值与(7)中所求出的最大列数z的值的乘积，同时使p>q ;
[0016](10)对(9)中所建立的矩阵进行步骤⑷中的化简运算；
[0017](11)计算化简后矩阵的秩，若秩不等于上面所求的最小秩r_min，则记录此时矩阵的列数和秩，将列数和秩两个数作为一组，依次记为1^、r2q2、r3q3L ;然后i加I,若i小于10则转入步骤(9)，否则转入步骤(12)，若秩等于上面所求出的最小秩r_min，则将最大列数z更新为此时的列数q，且给i赋值为2，转入步骤(8)；.r,-r min
[0018](12)设另一个循环变量flag为2 ;设比值变量rat1，使其
【权利要求】
1.一种BCH码做外码的级联码识别交织长度的盲识别方法，通过计算机进行数据读入、分析及计算处理，该方法步骤如下: (1)由计算机从待识别数据读入部分待识别的数据，设置一个表示列数循环时其最大值的变量q_max，其值为100 ； (2)设列数q的初值为10，此值应小于可能出现的最小的交织长度，其中交织长度等于BCH码长乘以交织宽度，BCH码阶数范围为3?8，交织宽度最小为2，因此最小交织长度为(23-1) X2 = 14 ； (3)建立一个pXq大小的矩阵，其中P是行数，q是列数，且p>q ; (4)对所建立的矩阵从左到右按列化简，若对角线上元素为1，则将此行依次与其下方每一行进行模二加运算；如果对角线元素为0，则寻找该列对角线下方的非零元素所在行，将非零元素所在行与当前行互换，再执行上述化简；如果对角线下方元素全为0，则不再化简； (5)计算化简后矩阵的秩，若秩不等于列数时则记录此时矩阵的列数和秩，否则不记录; (6)将列数q加I,当q>q_max时,转入步骤(7),否则转入步骤(3)； (7)比较保存的所有秩，求出最小的秩r_min;在此基础上，统计所有满足该最小秩的列数，求出该最小秩所对应的最大的列数z ； (8)设循环变量i的初值为2； (9)使用待识别数据重新建立一个PXq大小的矩阵，其中P为行数，q为列数(P和q跟上面的是一样的变量，只是赋值不同)，且q的值为⑶中i的值与(7)中所求出的最大列数z的值的乘积，同时使p>q ； (10)对(9)中所建立的矩阵进行步骤⑷中的化简运算； (11)计算化简后矩阵的秩，若秩不等于上面所求的最小秩r_min，则记录此时矩阵的列数和秩，将列数和秩两个数作为一组，依次记为--、--2、ιγι3? ;然后i加I,若i小于10则转入步骤(9)，否则转入步骤(12)，若秩等于上面所求出的最小秩r_min，则将最大列数z更新为此时的列数q，且给i赋值为2，转入步骤(8)；
/1 - r mi η (12)设另一个循环变量flag为2;设比值变量rat1，使其'_> = J^=^，其中F1、
仏-2Qi为步骤(11)中的A、q1 ；
厂" —厂/7 I (13)计冑nS 一 ；由于 flag 设为 2,故 rflag、rflairl、qflag、qflairl 分别对应为步骤(11)中的 r2、q2、q1 ；
厂" —厂// I (14)若ng _π 8一等于rat1,则flag加I,如果flag小于9则转入步骤(13)，若
iIflag flag-V? g _n ^一不等于rat1, q_max加50,转入步骤(2),若flag等于9则交织长度应等于此时最大列数Z。
【文档编号】H04L1/00GK104184557SQ201410459604
【公开日】2014年12月3日 申请日期:2014年9月11日 优先权日:2014年9月11日
【发明者】马丕明, 张丽媛, 杨勇 申请人:山东大学