一种基于KL散度的通信信号图域特征迭代提取方法技术

本发明专利技术公开了一种基于KL散度的通信信号图域特征迭代提取方法，利用通信信号的循环谱，在保证算法鲁棒性的前提下实现了特征序列的自动构建；具体讲，本发明专利技术首先通过图域映射理论，将通信信号的循环谱转换为一系列邻接矩阵，并提取邻接矩阵中所有元素构建特征序列备选集；然后对每一种调制类型，计算其在特征序列备选集中每一个索引处相对于其他调制类型的KL散度并相加，求得属于该调制类型的KL散度，根据每种调制类型的KL散度大小确定提取特征的顺序；按照顺序依次选择KL散度最大的索引作为相应调制类型的特征，每次提取一个特征都将其从特征序列备选集中删除，直到所有调制类型的特征序列构建完成。

【技术实现步骤摘要】
一种基于KL散度的通信信号图域特征迭代提取方法
本专利技术属于信号处理
，更为具体地讲，涉及一种基于KL散度的通信信号图域特征迭代提取方法。
技术介绍
自动调制分类(AMC)可以在很少或没有先验知识的情况下识别接收信号的调制类型，广泛应用于军事和民用通信。典型的自动调制识别方法通常分为两类：基于最大似然的方法(ML)和基于特征提取的方法(FB)。基于最大似然的方法是一种基于假设检验的理论，通过接收信号的似然函数，将似然比与一个门限值做比较做出判决，这种方法可以得到贝叶斯意义上的最优解，但是也存在诸多弊端；基于特征识别的方法包括特征提取(FE)和模式识别(PR)两个阶段，特征提取阶段从接收到的未知信号中提取若干参考特征，然后在模式识别阶段根据提取的特征判断信号的调制类型，这种方法虽然不是最优的但实施效率相比前者较高。但两种方法都需要系统提供较高的运算能力，难以用于一些实时性要求较高而系统资源受限的特殊应用场合；现有识别方法在处理实际无线通信信号时性能严重恶化，实际工程应用中鲁棒性差。基于图域的自动调制分类(AMCG)第一次将AMC变换到图形域，并且已经实现了比现有PR和基于LB的决策理论算法更优的性能。该方法利用调制信号的循环谱，依据循环频率将循环谱映射到图域，构建为加权有向环，手动记录邻接矩阵次对角线上的非零项，这些非零项被构建为有效特征参数。然而在AMCG中的整个图域特征构建是通过人工进行的，计算十分繁琐，工作量大，如果不适当地选择特征序列，容易造成较大的误差，通常会影响识别效果。这就需要一种科学地选择特征的方法，用于AMCG特征的自动构建。KL散度用来表示两个随机分布的差异，通过计算某种调制类型相对于其他调制类型的KL散度并根据其大小进行排序，可以实现图域信号特征的自动构建。采取迭代提取特征的方法，可以保证所有被提取的特征都是唯一的，而且当有新调制类型加入时，只需构建新调制类型的特征序列而不影响已有调制信号的特征序列。在避免人工参与特征构建的同时保证AMCG的鲁棒性和延展性。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种基于KL散度的通信信号图域特征迭代提取方法，在保证算法鲁棒性的前提下实现无人工参与的特征序列自动构建，并且可以延展到任意调制类型。为实现上述专利技术目的，本专利技术一种基于KL散度的通信信号图域特征提取方法，其特征在于，包括以下步骤：(1)、通信信号图域映射(1.1)、设置调制类型候选集Mdef，Mdef＝{M1,M2,...,MK}，其中，Mk表示第k种调制类型，k＝1,2,...,K，K表示调制类型总数；(1.2)、利用FAM算法计算通信信号x(n)的循环谱其中，α为循环频率，α的取值为[α1,α2,…,αp]，p为循环频率的取值个数，f为x(n)的频率，再对循环谱进行归一化和量化处理得到谱(1.3)、在谱中，根据的对称性，在每种调制类型下取α和f均为正的四分之一谱映射到图域，得到一个图集其中，表示第k种调制类型下循环频率为ατ时的一个图，τ＜p；将图集中的每一个图转换为一个邻接矩阵，建立起相应的邻接矩阵集由图转换成的邻接矩阵；同理，建立其余调制类型下的邻接矩阵集；(2)、构建特征序列备选集在所有邻接矩阵集中，取所有邻接矩阵的元素加入到特征序列备选集Idef中；Idef＝{β1,β2,…,βi,…,βI}其中，βi表示第i个元素，i＝1,2,...,I，I表示元素的最大数量；(3)、利用迭代法构建特征序列(3.1)、计算每种调制类型下的KL散度设置重复计算次数M；在第k种调制类型下第m次计算时，m＝1,2,…,M，根据Idef中索引i找到对应邻接矩阵中的一个元素，其值记为重复计算M次，得到第k种调制类型下的M个值同理，分别得到其余调制类型下的M个值；在第i个索引处，统计M次计算时第k种调制类型下的的概率分布：其中，其中，xk,i是每次计算时的随机变量；再对xk,i取绝对值|xk,i|，得：同理，计算出第i个索引处M次重复计算时，其余调制类型下的概率分布将当前第k种调制类型下的概率分布记为然后计算第i个索引处第k种调制类型相对与其他调制类型的联合KL散度同理，按照上述方法可以求得每一个索引i＝1,2...,I处第k种调制类型相对于其他调制类型的联合KL散度将所有i＝1,2...,I处的KL散度相加，得到第k种调制类型下的KL散度Ψk；同理，按照上述方法可以求得其余调制类型下的KL散度Ψ1、Ψ2、...、ΨK；(3.2)、确定特征提取顺序将所有调制类型下的KL散度进行升序排列，并作为特征提取顺序；(3.3)、迭代法构建特征序列从KL散度最小的调制类型开始依次提取特征，每种调制类型提取N个特征；1)、在KL散度最小的调制类型中，比较每一个索引i处的KL散度的大小，再从中提取KL散度最大的一个索引对应的元素作为该调制类型的一个特征；2)、将提取的特征从特征序列备选集Idef中删除；3)、重复步骤1)～3)直到该调制类型的第N个特征提取完毕，再用这提取的N个特征构建出该调制类型下的图域特征序列；4)、当上一种调制类型下的图域特征序列构建完成后，进行下一种调制类型下的图域特征序列构建，直到所有调制类型下的图域特征序列构建完成；(4)、当新调制类型加入时构建图域特征序列当新的调制类型加入调制类型候选集后，直接计算该种调制类型在剩余特征序列备选集中所有索引处的KL散度，按照步骤(3.3)所述方法构建该种调制类型下的图域特征序列。本专利技术的专利技术目的是这样实现的：本专利技术一种基于KL散度的通信信号图域特征迭代提取方法，利用通信信号的循环谱，在保证算法鲁棒性的前提下实现了特征序列的自动构建；具体讲，本专利技术首先通过图域映射理论，将通信信号的循环谱转换为一系列邻接矩阵，并提取邻接矩阵中所有元素构建特征序列备选集；然后对每一种调制类型，计算其在特征序列备选集中每一个索引处相对于其他调制类型的KL散度并相加，求得属于该调制类型的KL散度，根据每种调制类型的KL散度大小确定提取特征的顺序；按照顺序依次选择KL散度最大的索引作为相应调制类型的特征，每次提取一个特征都将其从特征序列备选集中删除，直到所有调制类型的特征序列构建完成。同时，本专利技术一种基于KL散度的通信信号图域特征迭代提取方法还具有以下有益效果：(1)、避免人工参与构建特征序列，不会因为人为地选择到不适当的特征而影响识别效果；且相比人工记录的方法，自动构建的特征序列包含较少的特征，不会导致任何的计算复杂性，但是达到更优的性能(2)、构建的特征序列在基于AMCG算法的整个分类过程中可以保持不变，并且可以忽略因训练信号传输符号随机性而导致的多种调制类型的调制特征序列的不一致性；(3)、提取的每一个特征都是唯一的，不与其他调制方式的特征重复；当有新的调制类型加入时，只需要提取新调制类型的特征，不会影响构建原有调制类型的特征序列。附图说明图1是本专利技术基于KL散度的通信信号图域特征迭代提取方法流程图；图2是BPSK调制类型下的图域映射示意图。具体实施方式下面结合附图对本专利技术的具体实施方式进行描述，以便本领域的技术人员更好地理解本专利技术。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本专利技术的主要内容时，这些描述在这里将被忽略本文档来自技高网...

【技术保护点】
1.一种基于KL散度的通信信号图域特征迭代提取方法，其特征在于，包括以下步骤：(1)、通信信号图域映射(1.1)、设置调制类型候选集Mdef，Mdef＝{M1,M2,...,MK}，其中，Mk表示第k种调制类型，k＝1,2,...,K，K表示调制类型总数；(1.2)、利用FAM算法计算通信信号x(n)的循环谱

【技术特征摘要】
1.一种基于KL散度的通信信号图域特征迭代提取方法，其特征在于，包括以下步骤：(1)、通信信号图域映射(1.1)、设置调制类型候选集Mdef，Mdef＝{M1,M2,...,MK}，其中，Mk表示第k种调制类型，k＝1,2,...,K，K表示调制类型总数；(1.2)、利用FAM算法计算通信信号x(n)的循环谱其中，α为循环频率，α的取值为[α1,α2,…,αp]，p为循环频率的取值个数，f为x(n)的频率，再对循环谱进行归一化和量化处理得到谱(1.3)、在谱中，根据的对称性，在每种调制类型下取α和f均为正的四分之一谱映射到图域，得到一个图集其中，表示第k种调制类型下循环频率为ατ时的一个图，τ＜p；将图集中的每一个图转换为一个邻接矩阵，建立起相应的邻接矩阵集由图转换成的邻接矩阵；同理，建立其余调制类型下的邻接矩阵集；(2)、构建特征序列备选集在所有邻接矩阵集中，取所有邻接矩阵的元素加入到特征序列备选集Idef中；Idef＝{β1,β2,…,βi,…,βI}其中，βi表示第i个元素，i＝1,2,...,I，I表示元素的最大数量；(3)、利用迭代法构建特征序列(3.1)、计算每种调制类型下的KL散度设置重复计算次数M；在第k种调制类型下第m次计算时，m＝1,2,…,M，根据Idef中索引i找到到对应邻接矩阵中的一个元素，其值记为重复计算M次，得到第k种调制类型下的M个值同理，分别得到其余调制类型下的M个值；在第i个索引处，统计M次计算时第k种调制类型下的的概率分布：其中，其中，xk,i是每次迭代时的随机变量；再对xk,i取绝对值|xk,i|，得：...

【专利技术属性】
技术研发人员：阎啸，王茜，张国玉，吴孝纯，刘冠男，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51