基于图的度特征的BPSK及QPSK信号调制识别方法技术

本发明专利技术针对BPSK和QPSK两种调制信号的识别问题，提出了一种基于图的度特征的BPSK及QPSK信号调制识别方法。首先获取待识别信号的平方谱，而后加矩形窗得到截断平方谱。之后，将截断平方谱进行图域转换，将该图各顶点度之和作为识别统计量，并设定适当门限，将识别统计量和门限相比较，以实现对BPSK和QPSK信号的调制识别。仿真结果表明，本发明专利技术在低信噪比条件下，对BPSK和QPSK两种调制信号能够进行更有效地识别，且由于图转换的输入信号点数较少，算法复杂度较低。法复杂度较低。法复杂度较低。

【技术实现步骤摘要】
基于图的度特征的BPSK及QPSK信号调制识别方法


[0001]本专利技术属于信号识别与处理领域，具体涉及一种基于图的度特征的BPSK及QPSK信号调制识别方法。

技术介绍

[0002]雷达、通信信号的调制方式识别广泛应用于军事和民用领域。在军事领域如电子战中，信号调制方式识别是获取敌方雷达性能参数的重要前提；而在民用领域，在无线通信频道资源有限的前提下，接收端实现对信号的调制方式识别可以更大程度地节省信道。BPSK和QPSK信号是雷达信号中最常用的两种相位编码调制信号。
[0003]近年来，一种新的基于图的信号处理算法在信号处理等方面得到较广泛的应用。特别地，Kun Yan等人于2017年提出一种将时间序列信号转换为图拓扑结构的处理算法(Yan K，Wu H C，Xiao H，et al.Novel Robust Band
‑
Limited Signal Detection Approach Using Graphs[J].IEEE Communications Letters，2017.)，该算法的主要思想是将噪声变成完全图，带噪信号转换为非完全图，以检测图的完全连接性来实现对信号的有无的检测。如果按照这个框架，将其应用于信号调制识别中，势必需要将BPSK或QPSK信号中的某一个时域或者其他变换域形式变成噪声的形式，而另一个变成非噪声的形式。其中一个可行的思路是：将BPSK信号及QPSK信号进行平方运算后，去掉若干较大值，可以得到各自的平方修正谱，并分别转换成完全图与非完全图，可将调制识别问题转化为完全图检测问题。但由于较大的谱线中包含了大部分信号分量的信息，这样处理，将损失掉这部分重要信息，导致在信噪比较低时，性能表现欠佳。此外，去掉多少根这样的大谱线本身也是一个难以折中的技术问题。

技术实现思路

[0004]本专利技术针对现有技术中的不足，提供一种基于图的度特征的BPSK及QPSK信号调制识别方法。该方法将截断平方谱进行图转换，选择特定的特征量及门限，完成BPSK及QPSK两种调制信号的识别，算法的计算复杂度低，且在低信噪比时识别正确率较高。
[0005]为实现上述目的，本专利技术采用以下技术方案：
[0006]基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于，包括如下步骤：
[0007]步骤1：计算待识别信号的平方谱；
[0008]步骤2：将平方谱的峰值谱线作为中心谱线，加矩形窗得到截断平方谱；
[0009]步骤3：对截断平方谱进行图转换，并提取该图所有顶点度之和作为识别特征量；
[0010]步骤4：设定相应的门限；
[0011]步骤5：将识别特征量和设定的门限相比较，进而识别出BPSK信号和QPSK信号。
[0012]为优化上述技术方案，采取的具体措施还包括：
[0013]进一步地，在步骤1中，设待识别信号为x(n)，将其作平方运算并做DFT变换，取频谱模的平方，得到平方谱，记为Y(k)＝(|DFT[x2(n)]|)2，k＝0，1，...，N
‑
1，其中N为信号的样
本点数。
[0014]进一步地，在步骤2中，对平方谱Y(k)，搜索出平方谱的峰值谱线位置为k＝k
max
，其中k表示谱线位置，k
max
表示峰值谱线位置，设定矩形窗长度为2d，以k＝k
max
处峰值谱线为中心加窗得到：
[0015][0016]将Y
′
(k)中为0的点数删去后得到截断平方谱为B(m)＝Y
′
(k)，0≤m≤2d。
[0017]进一步地，当信号的样本点数N为1024时，d取50。
[0018]进一步地，步骤3具体包括如下步骤：
[0019]步骤3.1：将截断平方谱B(m)作归一化处理，得到其归一化频谱：
[0020][0021]设定量化级数q，对归一化后的频谱进行均匀量化，即当i/q＜B
′
(m)＜i+1/q，0≤i≤q
‑
1时，量化后的频谱为U(m)＝i+1；
[0022]将U(m)转换到图域，构成图G(V，E)，其中图的顶点集V表示量化级{1，2，...，q}的映射，V＝{v1，v2，...v
q
}；图的边集合E＝{e
α，β
|ν
α
∈V，ν
β
∈V}，e
α，β
表示图的两个顶点之间的边；构成图G(V，E)的具体做法是：对每一个量化后样本U(m)，逐个遍历其与U(m+1)的电平关系，当存在v
α
到v
β
的电平跳变时，则两个顶点相连，e
α，β
＝1；反之，则两个顶点无连接，e
α，β
＝0；
[0023]步骤3.2：计算G(V，E)的度矩阵D，并提取其对角线元素构成图的度向量d＝(d1，d2，...，d
j
，...，d
q
)，其中d
j
为第j个顶点上连接的边数之和，计算度的和sum(d)作为识别特征量。
[0024]进一步地，在步骤4中，设定BPSK和QPSK信号的识别门限λ
evt
，当信号的样本点数N为1024，图顶点个数取10时，λ
evt
取10。
[0025]进一步地，在步骤5中，若识别特征量小于设定的门限，则信号为BPSK调制信号；否则为QPSK调制信号。
[0026]本专利技术的有益效果是：本专利技术对观测信号的平方谱适当加窗后进行图转换，并根据图的度矩阵和确定识别特征量，完成对BPSK信号和QPSK信号的识别，其中直接利用平方谱，不需要进行去大值的修正过程，可以较好地保留信号分量的信息。相比于传统的基于修正谱识别算法而言，本专利技术直接将两类信号平方的功率谱进行图域变换，有效地保存了信号的主要信息，在低信噪比条件下可对BPSK和QPSK信号进行有效识别，且受参数变化影响较小，具有一定的鲁棒性，算法效率较高。此外，本专利技术与现有的基于完全图检测的识别处理算法不同，从随机图理论的角度分析信号的差别、定义图域特征，进一步扩展了现有图域处理方法的应用领域，丰富该类方法的处理手段。
附图说明
[0027]图1是本专利技术的识别方法流程图。
[0028]图2表示不同信噪比下BPSK与QPSK信号截断平方谱生成图的度之和的均值。
具体实施方式
[0029]现在结合附图对本专利技术作进一步详细的说明。
[0030]图1是基于图的度特征的BPSK及QPSK信号调制识别方法的流程图。首先获取待识别信号的平方谱，而后将平方谱加窗得到截断平方谱。之后，将其转换到图域，以该图度的和作为识别统计量，设定适当门限，将识别统计量和门限相比较，实现对BPSK和QPSK信号的调制识别。仿真结果表明，在无信号先验信息的条件下，本专利技术可对BPSK和QPSK两种调制信本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于，包括如下步骤：步骤1：计算待识别信号的平方谱；步骤2：将平方谱的峰值谱线作为中心谱线，加矩形窗得到截断平方谱；步骤3：对截断平方谱进行图转换，并提取该图所有顶点度之和作为识别特征量；步骤4：设定相应的门限；步骤5：将识别特征量和设定的门限相比较，进而识别出BPSK信号和QPSK信号。2.如权利要求1所述的基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于：在步骤1中，设待识别信号为x(n)，将其作平方运算并做DFT变换，取频谱模的平方，得到平方谱，记为Y(k)＝(|DFT[x2(n)]|)2，k＝0，1，...，N
‑
1，其中N为信号的样本点数。3.如权利要求1所述的基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于：在步骤2中，对平方谱Y(k)，搜索出平方谱的峰值谱线位置为k＝k
max
，其中k表示谱线位置，k
max
表示峰值谱线位置，设定矩形窗长度为2d，以k＝k
max
处峰值谱线为中心加窗得到：将Y
′
(k)中为0的点数删去后得到截断平方谱为B(m)＝Y(k)，0≤m≤2d。4.如权利要求3所述的基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于：当信号的样本点数N为1024时，d取50。5.如权利要求3所述的基于图的度特征的BPSK及QPSK信号调制识别方法，其特征在于：步骤3具体包括如下步骤：步骤3.1：将截断平方谱B(m)作归一化处理，得到其归一化频谱：设定量化级数q，对归一化后的频谱进行均匀量化，即当i/q＜B
′

【专利技术属性】
技术研发人员：杨莉，胡国兵，赵敦博，姜志鹏，
申请(专利权)人：金陵科技学院，
类型：发明
国别省市：