调制识别一直都是非合作通信的重要一环,实际中大部分都采用统计识别的模式来完成调制识别[1]。这种方法需要先对信号的特征进行提取,再通过分类器分类[2]。因此目前对调制方式识别的研究方向大致分2类:一类是将研究的重点放在新特征的提取上,包括高阶累积量[3]、循环谱[4]、小波特征[5]、星座图特征[6]等;另一类是采用机器学习作为分类器,比如随机森林[7]、支持向量机[8]、卷积神经网络[9]、AdaBoost[10]等,这些分类器通过大量的数据训练与测试,与传统决策树分类器相比,在使用相同特征参数的条件下可以获得更高的正确识别率。近年来,随着特征的种类和提取方法越来越完善,一些基于经典特征提取算法的新特征被相继提出并应用于调制识别中。文献[11-12]中作者用循环谱特征图像分析方法,通过汉明距离测度枚举训练和测试数据所得到的特征参数之间的差异,成功实现了多种调制信号的分类工作。文献[13]将信号的各种熵特征与特征选择算法相结合,运用5种不同的集成学习分类器完成了对多种数字调制信号的分类工作。文献[14]提出一种压缩时间高阶循环累积量。文献[15]为了解决循环累积量特征在高采样率下使用受限问题,将压缩感知框架引入调制识别体系中,并且将文献[14]中的新特征作为对照目标进行研究。然而,传统的瞬时特征和决策树分类器由于原理简单、容易硬件实现,目前仍然被广泛使用并被继续深入研究。本文针对2ASK、4ASK、8ASK、2PSK、4PSK、8PSK这6种数字调制方式,在文献[16]提出的瞬时特征基础上,利用MPSK调制的非线性相位取值的差异性提取特征,完成了4PSK与8PSK的类内识别;同时从瞬时频率的角度提取特征完成了4ASK和8ASK的调制识别。
1 调制识别流程基于统计模式理论的调制识别方法主要包含3个核心技术单元,即信号预处理、特征参数提取和分类器判决。本文的重点在于特征参数提取,因此文中将着重介绍新提出的区分4PSK和8PSK的瞬时特征——分段相位均值最大值,并从理论和仿真的角度完成对此特征的分析。
1.1 调制信号瞬时特征设
$x\left( t \right) = a\left( t \right)\cos \left( {2{\text{π}}{f_{\rm{c}}}\left( t \right)t + {\varphi _{{\rm{NL}}}}\left( t \right)} \right)$ |
式中:
$z\left( t \right) = x\left( t \right) + H\left[ {x\left( t \right)} \right]$ |
式中
则信号的瞬时幅度为
${A_{{\rm{instan}}}} = \sqrt {{x^2}\left( t \right) + {H^2}\left[ {x\left( t \right)} \right]} $ |
瞬时相位为
$\varphi \left( t \right) = \arctan \frac{{H\left[ {x\left( t \right)} \right]}}{{x\left( t \right)}}$ | (1) |
鉴于由式(1)得到的信号相位限制在
对于待识别的6种调制方式,一般通过提取信号的瞬时幅度、瞬时频率、瞬时相位等5个瞬时特征来完成调制识别。对于本文所使用的传统特征参数,因为在其他论文中也多次提及,所以这里只列出所使用的特征,并不进行详细介绍。
1)归一化中心化瞬时幅度功率谱密度最大值:
${R_{\max }} = {{\max {{\left| {{\rm{FFT}}\left( {{a_{{\rm{cn}}}}\left( i \right)} \right)} \right|}^2}} / {{N_{\rm{s}}}}}$ |
2)非弱信号段中心化瞬时相位非线性分量绝对值的标准偏差
${\sigma _{{\rm{ap}}}} = \sqrt {\frac{1}{C}\left( {\sum\limits_{{a_n}\left( i \right) > {a_t}} {\phi _{{\rm{NL}}}^2\left( i \right)} } \right) - {{\left( {\frac{1}{C}\sum\limits_{{a_n}\left( i \right) > {a_t}} {\left| {{\phi _{{\rm{NL}}}}\left( i \right)} \right|} } \right)}^2}} $ |
3)非弱信号段中心化瞬时相位非线性分量的标准偏差
${\sigma _{{\rm{dp}}}} = \sqrt {\frac{1}{C}\left( {\sum\limits_{{a_n}\left( i \right) > {a_t}} {\phi _{{\rm{NL}}}^2\left( i \right)} } \right) - {{\left( {\frac{1}{C}\sum\limits_{{a_n}\left( i \right) > {a_t}} {{\phi _{{\rm{NL}}}}\left( i \right)} } \right)}^2}} $ |
该特征表征信号瞬时相位的变化情况。
4)中心化瞬时频率的均方值:
$\mu _{20}^f = {\rm{E}}{\left[ {{f_{\rm{N}}}\left( i \right)} \right]^2}$ |
式中
该特征可以完成4ASK、8ASK这2种ASK信号的类内识别。对于ASK信号而言,只有符号0对应的调制信号是没有载频的,其余符号对应的调制信号都存在载频,只是幅度信息不一样。这样在求瞬时频率时,符号0对应的瞬时频率为0,其余的符号所对应的就是正弦波2个相邻采样点之间的相位差值。而从统计角度分析,2ASK取每个符号的概率为0.5,4ASK为0.25,8ASK为0.125。当每个符号等概率出现时,2ASK、4ASK、8ASK这3种调制方式里符号0出现的概率依次是0.5、0.25、0.125。对于瞬时频率来说,符号0出现的概率越大,相应的瞬时频率的离散程度就越大。所以2ASK的该特征值最大,4ASK其次,8ASK相对来说最小。此特征在不同信噪比下的仿真结果如图1所示,从图中可以看出,在信噪比大于10 dB时,4ASK与8ASK的特征值差异明显,可以用此特征区分。
![]() |
Download:
|
图 1 |
对于4PSK和8PSK这2种利用多个不同相位来传递基带序列的调制方式,目前已有的调制识别方法是利用高阶累积量[3],或者使用瞬时相位平方的标准差去区分。一般高阶累积量需要数量较多的码元,才能将相位调制所具有的特征展现出来。文献[18]提出用特征瞬时相位非线性分量平方值的标准偏差
![]() |
Download:
|
图 2 |
![]() |
Download:
|
图 3 |
由图2可知,此特征并不稳定,在实际应用中很难找到比较合适的门限来区分这2种调制方式。而从图3可以看出,在对此特征值多次取平均后,可以在信噪比高于5 dB时完成4PSK与8PSK的区分。针对该特征波动范围大,不适合用决策树分类这些缺点,本文提出了一种新的特征提取方法,即在已知符号速率的前提下,针对瞬时相位特征进行提取,并以此完成PSK的类内识别。
对于2PSK、4PSK、8PSK这3种相位调制,在一个符号内,3种相位调制的相位分量保持不变,在符号变化的位置,其相位值也发生了突变。在不考虑相位突变位置所产生的尖峰毛刺前提下,2PSK瞬时相位的最大值为
但是,在有噪声的情况下,每一个符号内的相位信息受噪声的干扰较大,图4所示是4PSK调制在信噪比为5 dB时的非线性相位。从图中可以看出,噪声的存在破坏了信号原本的特征,无法准确地确定相位的最大值。
![]() |
Download:
|
图 4 信噪比为5 dB下4PSK相位 |
因此,为了降低噪声的干扰,在假设已经完成符号速率估计的情况下,对于计算得到的信号非线性相位序列进行分段,每段包含
在无噪声的条件下,仿真得到的3种MPSK调制的该特征取值如图5所示。图中2PSK、4PSK、8PSK的特征值分别为5、6、7,与上述的理论值存在差异。
![]() |
Download:
|
图 5 SPAM随仿真次数变化曲线 |
仿真得到2PSK相位如图6所示。从图中可以看出2PSK的2个相位取值是
![]() |
Download:
|
图 6 2PSK相位 |
当采样速率是载频的8倍时,4PSK和8PSK仿真得到的瞬时相位如图7、8所示。从图中可以看出,在8倍采样率的条件下,4PSK的部分相位取值为1.568 5、3.142 8、4.697 4、6.28,与理论值相差
![]() |
Download:
|
图 7 4PSK相位 |
![]() |
Download:
|
图 8 8PSK相位 |
当采样速率分别为载频的16、32、64、128倍时,4PSK和8PSK相位取值如图9、10所示,部分相位取值已在图中标出。由图可知,对于4PSK和8PSK调制,采样速率越高,相位的取值越接近理论值。但因为在实际应用中并不会用过高的采样率进行采样,所以本文后续仿真都采用8倍采样进行,并考虑采样速率导致的相位偏差对特征值提取的影响。
![]() |
Download:
|
图 9 4PSK相位 |
![]() |
Download:
|
图 10 8PSK相位 |
对于4PSK和8PSK,由于8倍采样条件下会有一定误差,因此可以看到2种相位调制的相位取值都会达到
![]() |
Download:
|
图 11 SPAM随信噪比变化曲线 |
分类器采用传统的决策树分类器,此分类器因其结构简单、易理解、易扩展等优点而被广泛使用。针对本文的种调制信号,分类流程如下。
1)如果
2)对于MASK信号:如果
3)对于MPSK信号:如果
根据上述的特征和识别流程,在Matlab仿真下得到的6种调制方式的识别率如图12所示。仿真条件为基带符号速率
![]() |
Download:
|
图 12 决策树分类识别率 |
为了更好地评价本文所提出的新特征在对4PSK和8PSK这2种信号的分类性能,与文献[17]中所构造的具有相同功能的高阶累积量特征A1进行对照仿真实验。在不同信噪比下,此2种特征对以上MPSK信号的调制方式正确识别率如表1所示。
![]() |
表 1 2种特征提取算法识别率对照 |
根据表1中结果数据可知,当信噪比在−2 dB以下时,2种特征均已失效,出现了严重的误识别现象。但是当信噪比为2 dB时,本文提出的SPAM特征明显好于经典的高阶累积量特征A1,且随着信噪比增加,SPAM特征对于2种信号的识别率均大于97%,而高阶累积量特征A1由于受到码元数量的限制,在较高信噪比条件下的分类性能稍逊于SPAM特征。
4 结论本文针对待识别的6种调制方式,重点研究了用来区分4PSK和8PSK所提取出的新特征SPAM的计算流程,同时分析了特征SPAM的取值与理论值之间的差异以及产生原因。与现有特征
然而,码元速率的估计误差可能会影响此特征的分类性能,所以这将是下一步需要深入研究的问题。
[1] |
李蒙蒙. 基于LabVIEW的数字信号调制识别[D]. 北京邮电大学, 2018: 1-5.
(![]() |
[2] |
吴效威. 基于邻域粗集神经网络的通信信号调制识别研究[D]. 太原理工大学, 2011: 1-3.
(![]() |
[3] |
GHODESWAR S, POONACHA P G. An SNR estimation based adaptive hierarchical modulation classification method to recognize M-ary QAM and M-ary PSK signals[C]//Proceedings of the 2015 3rd International Conference on Signal Processing. Chennai, India, 2015.
(![]() |
[4] |
杜盼盼. 基于高阶累积量与循环谱分析的数字调制信号识别[D]. 兰州: 兰州理工大学, 2018: 2-6.
(![]() |
[5] |
HASSANPOUR S, PEZESHK A M, BEHNIA F. A robust algorithm based on wavelet transform for recognition of binary digital modulations[C]//Proceedings of the 2015 38th International Conference on Telecommunications and Signal Processing. Prague, Czech Republic, 2015.
(![]() |
[6] |
LI Qiang, SHEN Dong, WANG Fei. MQAM modulation recognition based on AP clustering method[J]. MATEC web of conferences, 2016, 44: 01002. DOI:10.1051/matecconf/20164401002 (![]() |
[7] |
LAU K, SALIBIAN-BARRERA M, LAMPE L. Modulation recognition in the 868 MHz band using classification trees and random forests[J]. AEU - international journal of electronics and communications, 2016, 70(9): 1321-1328. DOI:10.1016/j.aeue.2016.07.001 (![]() |
[8] |
HASSANPOUR S, PEZESHK A M, BEHNIA F. Automatic digital modulation recognition based on novel features and support vector machine[C]//Proceedings of the 2016 12th International Conference on Signal-image Technology & Internet-based Systems. Naples, Italy, 2016.
(![]() |
[9] |
KARRA K, KUZDEBA S, PETERSEN J. Modulation recognition using hierarchical deep neural networks[C]//Proceedings of 2017 IEEE International Symposium on Dynamic Spectrum Access Networks. Piscataway, USA, 2017.
(![]() |
[10] |
CHIKHA W B, CHAOUI S, ATTIA R, et al. Performance of AdaBoost classifier in recognition of superposed modulations for MIMO TWRC with physical-layer network coding[C]//Proceedings of the 2017 25th International Conference on Software, Telecommunications and Computer Networks. Split, Croatia, 2017: 1–5.
(![]() |
[11] |
YAN Xiao, FENG Guoyu, WU H C, et al. Innovative robust modulation classification using graph-based cyclic-spectrum analysis[J]. IEEE communications letters, 2017, 21(1): 16-19. DOI:10.1109/LCOMM.2016.2618868 (![]() |
[12] |
YAN Xiao, LIU Guannan, WU H C, et al. New automatic modulation classifier using cyclic-spectrum graphs with optimal training features[J]. IEEE communications letters, 2018, 22(6): 1204-1207. DOI:10.1109/LCOMM.2018.2819991 (![]() |
[13] |
ZHANG Zhen, LI Yibing, JIN Shanshan, et al. Modulation signal recognition based on information entropy and ensemble learning[J]. Entropy, 2018, 20(3): 198. DOI:10.3390/e20030198 (![]() |
[14] |
LIM C W, WAKIN M B. Compressive temporal higher order cyclostationary statistics[J]. IEEE transactions on signal processing, 2015, 63(11): 2942-2956. DOI:10.1109/TSP.2015.2415760 (![]() |
[15] |
XIE Lijin, WAN Qun. Cyclic feature-based modulation recognition using compressive sensing[J]. IEEE wireless communications letters, 2017, 6(3): 402-405. DOI:10.1109/LWC.2017.2697853 (![]() |
[16] |
AZZOUZ E E, NANDI A K. Automatic identification of digital modulation types[J]. Signal processing, 1995, 47(1): 55-69. DOI:10.1016/0165-1684(95)00099-2 (![]() |
[17] |
戈彤彤. 基于多维度特征的数字通信信号调制识别技术研究[D]. 哈尔滨: 哈尔滨工程大学, 2019: 11-29.
(![]() |
[18] |
DAHAP B I, LIAO Hongshu. Advanced algorithm for automatic modulation recognition for analogue & digital signals[C]//Proceedings of 2015 International Conference on Computing, Control, Networking, Electronics and Embedded Systems Engineering. Khartoum, Sudan, 2015.
(![]() |