卷积码编码参数全盲识别方法.pdf

(19)中华人民共和国国家知识产权局*CN103312457A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103312457103312457A(43)申请公布日2013.09.18(21)申请号201310169040.4(22)申请日2013.05.09(71)申请人西安电子科技大学地址710071陕西省西安市太白南路2号(72)发明人阔永红宋超陈健杨龙任超(74)专利代理机构陕西电子工业专利中心61205代理人王品华朱红星(51)Int.Cl.H04L1/00(2006.01)权权利要求书2页利要求书2页说明书5页说明书5页附图5页附图5页(54)发明名称卷积码编码参数全盲识别方法(57)摘要本发明公开了一种卷积码编码参数全盲识别方法，主要解决现有技术的未全盲识别，计算复杂度高的问题。其技术方案是：（1）根据卷积码编码器与信息码字之间的相互校验关系求解出码长n、输入信息位长度k及寄存器长度m；（2）当码率k/n是1/n时，依次取卷积码比特流中码字的相邻两位比特构造分析矩阵C，对其使用快速Walsh-Hadamard变换求得生成矩阵G；（3）当码率k/n是(n-1)/n时，取卷积码比特流中码字的第1位到第(k+1)位比特构造分析矩阵C，对其使用快速Walsh-Hadamard变换求得基本校验矩阵H，完成卷积码编码参数的全盲识别。本发明较好地解决了未知任何先验知识及含误码情形下的卷积码的编码参数全盲识别问题，简化运算复杂度，可用于多点广播通信以及智能通信。CN103312457ACN1032457ACN103312457A权利要求书1/2页1.一种卷积码编码参数全盲识别方法，包括如下步骤：1）对进入接收机的射频信号进行预处理，得到卷积码比特流Bs；2）通过卷积码比特流Bs构造出行数为L，列数为l的码率分析矩阵Rl，列数l从2遍历到43，行数L＝l+50，求解出每个码率分析矩阵Rl的秩r，若r＜l，则将当前码率分析矩阵Rl的列数l并入到列集合Q中，将当前码率分析矩阵Rl的秩r并入到秩集合R中；3）求解出列集合Q中所有列数l的最大公约数，即为卷积码的码长n，求解出秩集合R中任意两个相邻秩r的差，即为卷积码的信息位长度k，根据这些求得卷积码的寄存器长度，式中要求列数l与秩r选至同一个码率分析矩阵Rl；4）根据卷积码的码长n、信息位长度k、寄存器长度m和已知卷积码比特流Bs，对卷积码的基本校验矩阵H和生成矩阵G进行识别：4a）当卷积码的码率k/n等于1/n时,按如下步骤识别生成矩阵G：4a1）在伽罗华域GF(2)上，假设接收的卷积码比特流Bs为：c1,1c1,2...c1,nc2,1c2,2...c2,n......ci,1ci,2...ci,n......其中ci,n表示第i个码字的第n位比特，i∈[1,+∞)；4a2）取卷积码比特流Bs中每个码字的第一位和第二位比特构造列数为p＝2·(m+1)，行数为q＝p2的分析矩阵C，形式如下：,其中i∈[1,+∞)；4a3）将分析矩阵C的每个行向量以行为单位转化为十进制数并放入集合I中，集合Ipp中数“j”出现的次数，即为频率向量va的第2-j个元素的值，j=0，...,2-1，从而构造出频率向量va；4a4）对频率向量va进行快速Walsh-Hadamard变换，得到解向量B；4a5）将解向量B中最大元素对应的除第1列外的列值减1，并转化为二进制形式，作为地址向量vb，将地址向量vb中的元素逆序排列得逆地址向量vb'，假设逆地址向量(1.1)vb'=[v1v2v3...v2m+1v2m+2]，则卷积码码字第一位比特的生成矩阵g=[v1v3...v2m+1]和第二(1.2)位比特的生成矩阵g=[v2v4...v2m+2];4a6）取卷积码比特流Bs中每个码字的第二位和第三位的比特，重复步骤4a2）-4a5）的过程，直到识别出第n-1位比特的生成矩阵g(1,n-1)和第n位比特的生成矩阵g(1,n)，完成对卷积码生成矩阵G=[g(1,1);g(1,2);...g(1,n-1);g(1,n)]的识别。4b）当卷积码的码率k/n等于(n-1)/n时，按如下步骤识别基本校验矩阵H：4b1）在伽罗华域GF(2)上，假设接收的卷积码比特流Bs为：c1,1c1,2...c1,nc2,1c2,2...c2,n......ci,1ci,2...ci,n......其中ci,n表示第i个码字的第n位比特，i∈[1,+∞)；4b2）取卷积码比特流Bs中每个码字的第一位到第(k+1)位比特构造列数为p＝n·(m+1)，行数为q＝p2的分析矩阵C，形式如下：2CN103312457A权利要求书2/2页其中i∈[1,+∞)；4b3）按步骤4a3）对分析矩阵C作相同处理构造出频率向量va；4b4）