一种卫星链路线性分组码的盲识别方法.pdf

本发明涉及一种卫星链路线性分组码的盲识别方法，该方法包括下列顺序的步骤：基于线性分组码码重分布概率和循环码的最大公因式次数分布概率估计码组长度和起始点；基于部分Walsh-Hadamard变换法及其改进算法估计校验矩阵和生成矩阵。本发明采用Walsh-Hadamard变换的方法识别分组码的校验矩阵和生成矩阵，Walsh-Hadamard变换值表示其对应的地址向量作为方程组解向量时，方程组中成立方程的个数和不成立方程个数之差，因此Walsh-Hadamard变换值的最大值对应的地址向量即为方程组的解向量，该方法能够适应在误码率较高的情况。

权利要求书一种卫星链路线性分组码的盲识别方法，该方法包括下列顺序的步骤：（1）基于线性分组码码重分布概率和循环码的最大公因式次数分布概率估计码组长度和起始点；（2）基于部分Walsh‑Hadamard变换法及其改进算法估计校验矩阵和生成矩阵。根据权利要求1所述的卫星链路线性分组码的盲识别方法，其特征在于：所述的码重分布概率采用统计方法估计，变换码组长度为和起始点为当两个码重分布之间的相关系数ρ最小时，表示两个码重分布概率之间的差异越大，则此时的和是真实的码组长度和起始点。根据权利要求1所述的卫星链路线性分组码的盲识别方法，其特征在于：所述的最大公因式次数分布概率采用统计方法估计，变换码组长度为和起始点为当两个最大公因式次数分布概率之间的相关系数ρ最小时，表示两个最大公因式次数分布概率之间的差异越大，则此时的和是真实的码组长度和起始点。根据权利要求1所述的卫星链路线性分组码的盲识别方法，其特征在于：所述的Walsh‑Hadamard变换采用蝶形运算，设参加Walsh‑Hadamard变换的位数为nw，其运算量为次加/减运算，状态个数为对(n,k)线性分组码，Walsh‑Hadamard变换及其改进算法的状态个数和运算量如表1所示：表1三种Walsh‑Hadamard变换方法比较根据权利要求2所述的卫星链路线性分组码的盲识别方法，其特征在于：所述基于码重分布概率估计码组长度和起始点的步骤为：a）初始化估计的码组长度和起始点b）对接收到的序列进行分组；c）计算各码组的重量；d）统计各重量的码组数量，并除以码组总数得到码重分布概率；e）利用实际序列的码重分布概率和随机序列的理论码重分布概率之间的相关系数估计码组长度和起始点。根据权利要求3所述的卫星链路线性分组码的盲识别方法，其特征在于：所述基于最大公因式次数分布概率估计码组长度和起始点采用以下步骤：a)依据估计的码组长度和起始点对接收到的序列进行分组；b)计算各码组与其循环右移形成的码字的最大公因式次数，循环右移的位数取3位；c)统计各最大公因式次数的码组数量，并除以码组总数得到分布概率；d)计算最大公因式次数分布概率相关系数e）寻找最小对应的根据权利要求4所述的卫星链路线性分组码的盲识别方法，其特征在于：当码长较短，即n＜24时，采用部分Walsh‑Hadamard变换法识别线性分组码的校验矩阵，其步骤为：a)读取码元数据并对接收到的数据进行分组；b)取出系数矩阵的前k列得到C0，并初始化i=1；c)取出系数矩阵的第k+i列得到C′,进行状态统计；d)利用部分Walsh‑Hadamard变换，计算不同的解向量H′i，并判断循环次数i是否小于n‑k，若判断结果为否，则执行下一步骤，否则i=i+1,并返回步骤c)；e)将H′i和单位矩阵组成校验矩阵。根据权利要求4所述的卫星链路线性分组码的盲识别方法，其特征在于：当码长较长，即n≥24时，采用改进的部分Walsh‑Hadamard变换法识别线性分组码的校验矩阵，其步骤为：a)读取码元数据并对接收到的数据进行分组；b)取出系数矩阵的前k列得到C0，并初始化i=1；c)取出系数矩阵的第k+i列得到C′,进行状态统计；d)C0分解成前后两部分即C0＝[C1 C2],让C2参加Walsh‑Hadamard变换，得到解向量，并判断循环次数i是否小于n‑k，若判断结果为否，则执行下一步骤，否则i=i+1,并返回步骤c)；e)将二进制向量和解向量按照顺序组合起来得到一组解向量H′i；f)将H′i和单位矩阵组成校验矩阵。

说明书一种卫星链路线性分组码的盲识别方法
技术领域
本发明涉及智能移动通信、多点广播通信和非协作通信领域领域，尤其是一种卫星链路线性分组码的盲识别方法。
背景技术
在数字通信系统中，将信息从信源传送给信宿的过程中，通信信道中的噪声和干扰会不可避免地对传输信息产生不同程度的干扰，信道编码技术主动地在所传输的信息中增加一些冗余，使其具有自动检错或纠错能力来克服干扰。信道编码按照编码方法主要可分为线性分组码、卷积码、Turbo码和LDPC码等，分组码是最早得到研究、也是成果最丰富、应用最广泛的一类码，其中线性分组码是分组码中最重要的一类码，它有着明显的数学结构，是其他编码方式的基础，它也广泛应用在卫星通信系统中。
在未来的智能移动通信、多点广播通信等领域中，自适应调制编码是复杂信道下实现信息高效可靠传输的重要手段之一，根据信道变化动态改变信道编码方式，可以获得最优的通信效率和服务质量。在这种通信环境中，一般难以通过协议实现多方通信单元的同步联络，由于调制和编码的格式是随着信道的不同质量情况而变化，因此需要接收方能够仅通过部分信号数据实现信道编码参数的快速盲识别，达到智能通信的目的。非协作通信领域传输过程由于受信道、时延等因素的影响，可能造成相关控制信息不能及时准确地到达接收方，接收方为了实时获取传输信息，就需要研究信道编码的盲识别技术。目前对信道编码的识别研究大多集中在卷积码的识别上，而对分组码识别的研究较少，只能识别低码率的线性分组码，并且适应的误码率较低。
发明内容
本发明的目的在于提供一种在误码情况下，能够实现不同码长线性分组码和循环码识别的卫星链路线性分组码的盲识别方法。
为实现上述目的，本发明采用了以下技术方案：一种卫星链路线性分组码的盲识别方法，该方法包括下列顺序的步骤：
（1）基于线性分组码码重分布概率和循环码的最大公因式次数分布概率估计码组长度和起始点；
（2）基于部分Walsh‑Hadamard变换法及其改进算法估计校验矩阵和生成矩阵。
由上述技术方案可知，本发明将分组码分为循环码和非循环码两种情况，对于循环码来说，可采用码字与其循环右移后码字的最大公因式次数分布概率相关法识别码组长度和起始点，由于最大公因式次数分布概率与循环码的码率没有关系，因此这种该方法对高码率的循环码仍然有效；采用Walsh‑Hadamard变换的方法识别分组码的校验矩阵和生成矩阵，Walsh‑Hadamard变换值表示其对应的地址向量作为方程组解向量时，方程组中成立方程的个数和不成立方程个数之差，因此Walsh‑Hadamard变换值的最大值对应的地址向量即为方程组的解向量，该方法能够适应在误码率较高的情况。
附图说明
图1是基于码重分布概率相关法进行线性分组码码组长度和起始点估计流程图；
图2是基于最大公因式次数分布概率相关法进行循环码码组长度和起始点估计流程图；
图3、4均为线性分组码校验矩阵识别的流程图；
图5是不同k值时(15,11)线性分组码部分Walsh‑Hadamard变换的结果图；
图6是不同i值时(15,11)线性分组码部分Walsh‑Hadamard变换的结果图；
图7是解向量矩阵；
图8是(63,36)线性分组码的改进部分Walsh‑Hadamard变换的结果图。
具体实施方式
一种卫星链路线性分组码的盲识别方法，该方法包括下列顺序的步骤：（1）基于线性分组码码重分布概率和循环码的最大公因式次数分布概率估计码组长度和起始点；（2）基于部分Walsh‑Hadamard变换法及其改进算法估计校验矩阵和生成矩阵。
如图1所示，根据线性分组码的码重分布概率与随机序列的码重分布概率有很大差距的特性，利用二者之间的差距估计码组长度和起始点。所述的码重分布概率采用统计方法估计，变换码组长度为和起始点为当两个码重分布之间的相关系数ρ最小时，表示两个码重分布概率之间的差异越大，则此时的和是真实的码组长度和起始点。
所述基于码重分布概率估计码组长度和起始点的步骤为：
a）初始化估计的码组长度和起始点
码组长度的变化范围可以设定为1~nmax，nmax为可能最大的码组长度。nmax选取的过小会使真实的码组长度不在变化范围内，选取的过大会增加计算量，实际应用时应综合考虑。实际中，分组码的码长一般不会超过255，起始点的变换范围可设为
b）对接收到的序列进行分组；
c）计算各码组的重量；
若Ai是分组码(n,k)中汉明重量为i的码字数目，则集合{A0,A1,…,An}称为该分组码的重量分布，简称码重分布。
d）统计各重量的码组数量，并除以码组总数得到码重分布概率；
若码的重量分布为{Ai}＝{A0,A1,…,An}，则其码重分布概率为