一种基于多维特征提取的通信调制信号识别方法与流程

本发明涉及通信领域，具体涉及一种基于多维特征提取的通信调制信号识别方法。

背景技术：

随着通信技术日益进步，通信环境渐渐变得复杂，通信信号的调制模式也变得更加多元化。为了信息可以更加有效、可靠地在信道中传输，信号会采用各种复杂严密的调制方式，包含不同的参数、不同的调制模式。由于这些内涵丰富的调制信号的产生，无线电也不再具有单一功能、单一频段的特征了。在没有充足的先决条件的情况下将无线通信信号的调制方式通过特定的方法识别出来的过程，就是信号调制的分类识别。接收机收到的是掺杂在一起的混合信号。通信信号所处的外界环境中存在着噪声的干扰，导致调制信号的模式识别也变得日益困难，此时如何克服外界的干扰信号，如何从接收的信号中顺利地提取出我们所需要的有用信号是一个重要的问题。通过学者们的不断研究，调制信号的分类识别方法趋于稳定和成熟，包括基于最大似然理论的方法、基于参数提取的方法等。

现有的基于最大似然理论的常见识别方法如下：基于最大似然理论将cpm分解和em结合，且信道参数是可估计的^[1]；基于似然函数在加性高斯白噪声下识别不同阶数的调相信号的方法^[2]；利用对数似然方法在加性高斯白噪声的环境下对正交调幅信号识别的方法^[3]；i-q数据在理想情况下对调制信号的识别具有统计方面的意义^[4]；利用了检测理论来识别调制信号模式，得出其似然比值并且以阈值的检测来代替原本的识别问题^[5]。

常见的基于特征参数提取的调制识别方法有：利用零中心归一化瞬时频率绝对值方差和瞬时频率发散度等方法来实现模拟调制信号和数字调制信号同时识别的方法^[6]，方法中还结合中心化非线性相位标准差、零中心归一化瞬时幅度绝对值方差和归瞬时幅度等参数共同识别调制信号的调制模式；利用非线性变换来识别调制信号的方法^[7]，首先对截获的未知调制信号进行预处理，接着仅用该信号的差分序列方差这一个参数作为信号特征来区分模拟和数字信号；利用haar小波变换得出信号符号率，从而得到信号参数特征进行模式识别^[8]，该方法有利于滤除通信媒介中的无用信息，并且可以保证一定的准确率；利用高阶累积量混合分形盒维数识别出了二进制调幅信号和多进制调频信号^[9]，避免了其他方法无法一步识别出此几种类型的调制信号，且此方法易于理解，结果准确。

引证文献：

[1]吴斌,袁亚博,汪勃.基于记忆因子的连续相位调制信号最大似然调制识别[j].电子与信息学报,2016,38(10):2546-2552.

[2]huancy,polydorosa.likelihoodmethodsformpskmodulationclassification[j].ieeetransactionsoncommunications,1995,43(234):1493-1504.

[3]yangy,liuch,soongtw.alog-likelihoodfunction-basedalgorithmforqamsignalclassification[j].signalprocessing,1998,70(1):61-71.

[4]weiw,mendeljm.maximum-likelihoodclassificationfordigitalamplitude-phasemodulations[j].ieeetransactionsoncommunications,2000,48(2):189-193.

[5]李侃,孙进平,邸杰.一种psk/qam调制方式识别的似然比方法[j].电子测量技术,2007,30(2):6-7.

[6]azzouzee,nandiak.automaticmodulationrecognitionofcommunicationsignals[m].kluweracademicpublishers,1996.

[7]姚亚峰,陈建文,黄载禄.模拟与数字调制方式的非线性变换识别方法[j].电讯技术,2005,45(1):23-26.

[8]chanyt,plewsjw,hokc.symbolrateestimationbythewavelettransform[c]circuitsandsystems,1997.iscas'97.,proceedingsof1997ieeeinternationalsymposiumon.ieee,1997:177-180vol.1.

[9]党月芳,徐启建,张杰,等.高阶累积量和分形理论在信号调制识别中的应用研究[j].信号处理,2013,29(6):761-765.

技术实现要素：

本发明的目的是根据上述现有技术的不足之处，提供一种基于多维特征提取的通信调制信号识别方法，该方法结合小波变换、高阶累积量、瞬时参数特征这三种算法的多维特征提取识别方法，可有效识别信号的调试方式。

为了实现上述目的，本发明提供了一种基于多维特征提取的通信调制信号识别方法，其包括以下步骤：

(s1)通过haar小波变换得到接收信号的符号率信息以及码元的突变信息；

(s2)计算接收信号的高阶累积量；

(s3)计算接收信号的瞬时频率特征参数；

(s4)根据接收信号的符号率、码元突变信息、高阶累积量以及瞬时频率特征参数，采用决策树的方式判断接收信号的调制方式。

本发明的进一步改进在于，对接收信号进行haar小波变换包括以下步骤：

(s11)对接收信号进行下变频预处理，使得高频调制信号频谱搬移至基带；

(s12)对预处理后的信号进行haar小波变换，根据变换结果得到符号率信息；

(s13)对预处理后的信号进行傅里叶变换，根据变换结果得到码元的突变信息。

本发明的进一步改进在于，计算接收信号的高阶累积量的过程中，分别计算接收信号的各阶累积量。

本发明的进一步改进在于，计算接收信号的瞬时频率特征参数的过程中，定义中心归一化瞬时频率丰富度μf⁴²，并将其作为瞬时频率特征参数，瞬时频率丰富度μf⁴²的计算公式为：

式中f(i)为接收信号。

本发明的进一步改进在于，判断判断接收信号的调制方式的过程中，采用如下方式：

(1)根据符号率信息判断接收信号传输的是数字信号还是模拟信号；或/和根据码元的突变信息判断收信号传输的是数字信号还是模拟信号；

(2)对于接收信号的高阶累积量，根据接收信号的各阶累积量判断接收信号的调制方式；

(3)根据接收信号的瞬时频率特征参数，判断接收信号的调频信号的进制。

本发明的优点是：本项发明从提高算法计算效率出发，在基于特征提取的识别方法下，提出结合小波变换、高阶累积量、瞬时参数特征这三种算法的多维特征提取识别方法。此方法易于理解、精确度高、计算量低，易于实现。

附图说明

图1为dsb小波变换频谱图；

图2为2psk小波变换频谱图；

图3为4psk小波变换频谱图

图4为2fsk小波变换频谱图

图5为4fsk小波变换频谱图

图6为8fsk小波变换频谱图

图7为高阶累积量参数f1；

图8为高阶累积量参数f2；

图9为调频信号的μf⁴²趋势图。

具体实施方式

以下结合附图通过实施例对本发明的特征及其它相关特征作进一步详细说明，以便于同行业技术人员的理解：

实施例：如图1所示，本发明的实施例包括一种基于多维特征提取的通信调制信号识别方法，该方法包括以下步骤：

(s1)通过haar小波变换得到接收信号的符号率信息以及码元的突变信息；

(s2)计算接收信号的高阶累积量；

(s3)计算接收信号的瞬时频率特征参数；

(s4)根据接收信号的符号率、码元突变信息、高阶累积量以及瞬时频率特征参数，采用决策树的方式判断接收信号的调制方式。

步骤(s1)主要采用小波变换的手段。小波变换是一个分析数据非常灵活的工具，用小波变换对信号进行分析能够利用其窗函数宽度随意调整的特点对信号进行不同精度的观测，提高使用效率。

连续小波变换的公式如下所示：

上式中a是尺度因子，b是位移因子。g(t)∈l²(r)是小波基函数，该函数满足：

上式中g(t)是小波基函数，g(ω)为g(t)的傅里叶变换，g(b,a)(t)为g(t)通过线性的伸缩平移变换而得到的函数集，具体如下所示：

常见的一维小波有haar小波、mexicanhat小波和morlet小波。由于haar小波对于信号时域、频域的局部突变信息均有较为精准的检测能力，故在信号处理的研究和应用中主要采用haar小波来提取信息的局部特征。本实施例中选用haar小波来提取信号的时域跳变信息，以此对待识别信号码元速率的存在与否进行判断。同时小波变换还可得到待测信号的符号率有关信息，从而判断出该信号是模拟信号还是数字信号。然而对于已经被调制过的高频信号来说，单纯地对信号进行小波变换并不适用于所有的频率，故采用将归一化后的接收信号的频谱搬移到基带的方法，使小波变换更具频率兼容性。步骤(s1)的算法具体步骤如下：

(s11)对接收信号进行下变频预处理，使高频调制信号频谱搬移至基带；

(s12)对预处理后的信号进行haar小波变换，根据变换结果得到符号率信息；

(s13)对预处理后的信号进行傅里叶变换，使码元的突变更直观的体现出来，并根据傅里叶变换的变换结果得到码元的突变信息。如图1至6所示，图1为一种模拟调制信号的小波变换频谱图，图2至6均为数字调制信号的小波变换频谱图。根据图1至6的对比可知，模拟调制信号与数字调制信号的小波变换频谱图有显著的差异。

如图7、8所示，在上述步骤(s2)中，计算出的高阶累积量可作为判断依据。高斯白噪声在高阶累积量的阶数多于二阶时的累积量数值为0，因此高阶累积量可抑制其他方法中无法避免的高斯噪声，在近年来被广大学者所熟知，并广泛应用于调制识别中。因调制信号的高阶累积量数值只取决于信号本身的调制方式和信号的功率，故可对同一个信号的不同阶数累积量进行简单运算，得出一个各信号统一的比值，以区分不同调制信号。基于高阶累积量的识别方法具体步骤如下：

(1)设x(k)为零均值独立同分布的复随机序列，其p阶混合矩为mpq＝e{x(k)^p-q[x^*(k)]^q}，则其高阶累积量如下所示：

c20＝cum[x(k)x(k+l1)](4)

c21＝cum[x(k)x^*(k+l1)](5)

c40＝cum[x(k)x(k+l1)x(k+l2)x(k+l3)](6)

c41＝cum[x(k)x(k+l1)x(k+l2)x^*(k+l3)](7)

c42＝cum[x(k)x(k+l1)x^*(k+l2)x^*(k+l3)](8)

当l1＝l2＝l3＝...＝0时，上述五个式子可简化为：

c20＝m20(9)

c21＝m21(10)

c41＝m41-3m21m20(12)

因待测信号和识别方法选取的原因，本发明所需的高阶累积量只到四阶，所以更高阶的累积量理论知识就不再过多赘述。将各信号通过上述理论计算出各阶累积量，以便于下一步对信号特征的处理。

(2)常见的数字调制信号高阶累积量结果如下表所示，其中e表示调制信号自身的功率。

表1数字调制信号的高阶累积量值

通过观察数据，本文构造了fr＝[f1,f2]的特征矩阵来识别五种数字信号的调制模式，其中f1和f2的具体含义为：

f1＝c40/c42(14)

f2＝c41/c42(15)

(3)根据特征向量fr的值来分析不同的数字调制信号，即可看出非常明显的差异，结果如下：

若fr＝[1,1]，则可以判定待测调制信号是2psk；若fr＝[1,0]，则可以判定待测调制信号是4psk；若fr＝[0,0]，则可以判定待测调制信号是mfsk。由理论结果可知高阶累积量的矩阵构造方法的可以将数字调制信号分为2psk、4psk、2fsk/4fsk/8fsk三类。

如图9所示，在上述步骤(s3)中，得出的瞬时频率特征参数可作为数字调频信号的进制数的判断依据。想要识别不同进制的调频信号，就要用和信号频率相关的特征参数。不难想到，同为调频信号的2fsk、4fsk和8fsk，唯一的区别就在于码元取值的数量和调制信号中不同频率载波的数量。故本发明在此问题上采用瞬时频率丰富度的方法。定义中心归一化瞬时频率丰富度μf⁴²：

式中

由方差的计算公式可以看出，参数μf⁴²的值是随着mfsk信号进制数的增加而随之增大的，刚好符合理论分析得出的结论，故直接对信号进行此计算即可。

在上述的步骤(s4)中，判断接收信号的调制方式的过程中，采用如下方式：

(1)根据符号率信息判断接收信号传输的是数字信号还是模拟信号；或/和根据码元的突变信息判断收信号传输的是数字信号还是模拟信号；对于模拟调制信号，码元不跳变，也没有稳定的符号率；对于数字调制信号符号率应该是稳定的；

(2)对于接收信号的高阶累积量，根据接收信号的各阶累积量判断接收信号的调制方式；具体的判断过程参见步骤(s2)的详细描述；该过程适用于判断数字调制信号的调制类型；但是本步骤无法分辨调频调制的进制，例如2fsk/4fsk/8fsk；

(3)根据接收信号的瞬时频率特征参数，判断接收信号的调频信号的进制。公式(17)可得出一个数值，根据该数值可得出调频信号的进制；本步骤是步骤(2)的补充。

本发明深入研究了通信调制信号特征提取与分类技术，并且分析了各种调制信号的产生机理。通过查阅大量资料和文献，结合相关的理论背景和实际需要，提出了一种基于多维特征提取的通信调制信号识别方法。针对本文的仿真结果，结论如下：

下变频预处理将中高频调制信号频谱搬移至基带，可解决中频信号的小波变换中数字调制信号的傅里叶变换无法直接与模拟调制信号区分开来的问题，即数字调制信号无法呈现谱峰等距分布的特征，频谱图几乎与模拟调制信号一致。

直接对基带信号进行一次小波变换即可区分数字调制信号和模拟调制信号，不用对中高频信号进行二次小波变换或者对中高频信号进行成型整顿后再进行一次小波变换。下变频处理后的小波变换相比于其他基于小波变换的算法，在效率和精确度上都有很大的提高。

单纯地对信号进行普通的四阶、六阶甚至更高阶的累积量运算也难以直接对信号进行分类。本文通过对各信号理论高阶累积量数值的对比分析，构造出的基于高阶累积量的特征向量，就可轻松将复杂的问题简单化。该方法只需对待测信号进行最高阶为四阶的累积量运算，再将所得的累计量数值进行一些特定的简单运算，就可将待测信号做出分类，在降低了理解难度的同时提高了运算效率。

在对比了大量文献中提出的方法后得出，瞬时频率丰富度是对调频信号进行类内区分的最简便的方法。该方法依据阶数越高的调频信号包含的载频种类越多，在进行瞬时频率分析时，其丰富度数值就越高这一理论，对不同阶数的调频信号进行处理和对比，通过设定合适的阈值，来区分各阶信号。此方法具有精准的识别效果。

以上的本发明实施方式，并不构成对本发明保护范围的限定。任何在本发明的精神和原则之内所作的修改、同替换和改进等，均应包含在本发明的保护范围之内。