一种二进制本原 BCH 码盲识别方法.pdf

(19)中华人民共和国国家知识产权局*CN103401567A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103401567103401567A(43)申请公布日2013.11.20(21)申请号201310343717.1(22)申请日2013.08.08(71)申请人山东大学地址250100山东省济南市历城区山大南路27号(72)发明人马丕明王丹杨勇(74)专利代理机构济南金迪知识产权代理有限公司37219代理人许德山(51)Int.Cl.H03M13/15(2006.01)权利要求书4页权利要求书4页说明书8页说明书8页(54)发明名称一种二进制本原BCH码盲识别方法(57)摘要一种二进制本原BCH码盲识别方法，属于信道编码盲识别技术领域，采用线性矩阵分析法或比特相关性算法识别马场，再识别起点，线性矩阵分析法可在误码率较高时正确识别码长小于等于15的BCH码，比特相关性算法可快速识别出码率不太高的码字码长。码根特性算法通过统计码根的分布判断码长是否正确，可以识别出高码率的码字码长。识别起点时采用公因式算法。公因式算法的依据是无误码情况下所有码字的公因式是生成多项式或生成多项式的倍数。识别生成多项式时采用码根算法，最后遍历本原多项式译码，根据误码率得到本原多项式。本方法运算速度快，正确率高，容错性好，适用于所有二进制本原BCH码。CN103401567ACN1034567ACN103401567A权利要求书1/4页1.一种二进制本原BCH码盲识别方法，包括码长识别、起点识别、生成多项式识别和本原多项式识别，其具体步骤如下：〈1〉码长识别1)初始化参数，可能的起点：shiftg＝-1，取值范围为[-1∞]，shiftg的类型gchoice＝-1，取值范围为{-1,1}，识别的码长：ng＝-1，取值范围为[-1∞]，从本原多项式文件中读入本原多项式表，接收到的码字序列为R，长度为Rlength；2)用线性矩阵分析法识别码长是否小于等于15，定义q为3～50内的循环变量，将接收到的数据按行写入50×q的矩阵M中，如果M的秩不等于q，则保留q，计算保留数即q在3～50的循环中保留的满足上述条件的q的最大公因式作为识别的码长ng，如果没有保留数，则没有识别出码长；3)判断ng是否满足下列两个条件：ng+1是2的整数幂；ng＞0，如果ng不满足这两个条件，则进入步骤4），如果ng满足这两个条件，则识别到的码长就是ng，码长识别过程结束；4)若ng+1不是2的整数幂并且ng≤0，码长识别如下：（a）定义阶数m＝5，偏移量shift＝0，数据序列指针R0＝R；（b）设码长n＝2m-1,R＝R0+shift，使用比特相关性算法判断R中数据的码长是否是n：i.定义矩阵XN×n，矩阵函数为XN×n，将接受到的码字序列R的前N*n比特按行写入矩阵XN×n，然后将矩阵XN×n转置得到Xtn×N，Xtn×N的列数集合为N＝{0,1,…,n-1}，Xtn×N中信息集的列数集合为I1×n，Xtn×N冗余列数集合为J1×(N-n)，计算得到的码字个数num＝0，码重分布向量为W1×(N+1)，将其初始化为全0；ii.寻找集合I1×n，I∪J＝N,其中符号∪、∩、分别表示集合相加、取交集、空集合，I、J是两个集合，集合元素是矩阵Xtn×N的列数，Xt(I)表示I1×n中的元素对应的列构成的方阵，Xt(J)表示J1×(N-n)中元素对应的列构成的矩阵，如果Xt(I)能够通过行化简得到单位阵，则I1×n为正确的信息集，记录下Xt(I)的行化简步骤，对Xt(J)实施相同的行化简步骤，执行步骤iii；如果找不到I1×n使得Xt(I)能够通过行化简得到单位阵，则认为n为接收到的码字序列的码长，ng＝n，gchoice＝1，shiftg＝shift；iii.将{0,1,…,n-1}随机平分成两个向量Z1(1×k1)，Z2(1×k2)，k1+k2＝n，k1、k2分别是向量Z1,Z2的列数，Z1中元素对应的行组成的矩阵为Λ1，Z2中元素对应的行组成的矩阵为Λ2，定义sigma＝log(n)/log(2)，随机选取sigma个J1×(N-n)中的元素构成向量L，选取Λ1中的任意2行，共有种取法，对于遍取的所有的Λ1中的任意2行，将其在二元域上相加，得到向量Λ1|L，Λ1|(J\L)，其中Λ1|L表示L中的元素对应的Λ1的列构成的矩阵，J\L表示J-L,，L∈J,同样计算Λ2|L，Λ2|(J\L)；iv.比较所有的Λ1|L和Λ2|L，如果Λ1|L＝Λ2|L，则Λ1+Λ2的码重为weight＝weight((Λ1+Λ2)|(J\L))+4，W(weight)＝W(weight)+1，num＝num+1;其中：(Λ1+Λ2)是将Λ1、Λ2两个向量相加，(Λ1+