CN104639268B - Qpsk信号和oqpsk信号的识别方法 - Google Patents
Qpsk信号和oqpsk信号的识别方法 Download PDFInfo
- Publication number
- CN104639268B CN104639268B CN201410853901.5A CN201410853901A CN104639268B CN 104639268 B CN104639268 B CN 104639268B CN 201410853901 A CN201410853901 A CN 201410853901A CN 104639268 B CN104639268 B CN 104639268B
- Authority
- CN
- China
- Prior art keywords
- signal
- phase difference
- oqpsk
- instantaneous
- qpsk
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000010363 phase shift Effects 0.000 title abstract 4
- 238000005070 sampling Methods 0.000 claims description 25
- 238000011084 recovery Methods 0.000 claims description 9
- 239000000969 carrier Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 7
- 238000010606 normalization Methods 0.000 abstract 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000005034 decoration Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
一种QPSK信号和OQPSK信号的识别方法包括：从下变频后的I/Q信号提取出瞬时相位；基于所述瞬时相位计算相邻码元的相位差；基于所述相邻码元的相位差计算中心归一化瞬时相位差；基于所述中心归一化瞬时相位差计算所述中心归一化瞬时相位差的绝对值标准偏差；根据绝对值标准偏差进行信号识别，绝对值标准偏差大于零为QPSK信号，绝对值标准偏差等于零为OQPSK信号。本发明的技术方案可以简便有效的识别QPSK信号和OQPSK信号。
Description
技术领域
本发明涉及无线电监测领域，特别涉及一种QPSK信号和OQPSK信号的识别方法。
背景技术
在无线电监测领域，经常需要对非法无线电服务、干扰源等信号进行识别，只有识别出信号，才能为频谱管理提供依据。目前已有的调制识别方法主要是基于提取特征进行识别的方法，提取特征主要有以下几种：包含幅度、频率和相位的瞬时特征，统计量特征及变换域特征等。其中统计量特征(如高阶累积量)和变换域(如小波变换)特征工程计算复杂，实时性差。目前比较有效的识别特征方法是英国学者A.K.Nandi和Azzouz等人提出的数字调制信号识别算法(DMRAs)，采用的就是瞬时包络、相位和频率等关键特征参数作为调制识别的特征，该算法计算速度快，分类效果明显，工程上较易实现，但在某些信号之间没有相应的特征提取算法，如QPSK和OQPSK之间的识别，QPSK和OQPSK的相位特性完全相同，频谱特性也非常相似，DMRAs没有给出相应的特征提取参数。
因此如何简便有效的识别QPSK信号和OQPSK信号成为目前亟待解决的问题之一。
发明内容
本发明的技术方案解决的技术问题是如何简便有效的识别QPSK信号和OQPSK信号。
为了实现上述发明目的，本发明提供了一种QPSK信号和OQPSK信号的识别方法，包括：
从下变频后的I/Q信号提取出瞬时相位其中I(i)、Q(i)为正交下变频后的I路信号和Q路信号，0＜i＜Ns，Ns为采样点数；
基于所述瞬时相位计算相邻码元的相位差0＜i＜Nr，Nr为符号数，Nr＝Ns/(fs/rs)，rs为符号率，fs为采样率；
基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)：
基于所述中心归一化瞬时相位差ΔφNL(i)计算所述中心归一化瞬时相位差ΔφNL(i)的绝对值标准偏差σaΔp：其中，an(i)是经过定时恢复后信号的瞬时归一化幅度，ma是经过定时恢复后的信号的瞬时幅度的平均值，ta是判断弱信号段的一个幅度的判决门限电平，C是在集合{φNL(i)}中an(i)＞ta的an(i)信号的个数，0＜i＜Nr；
根据所述绝对值标准偏差σaΔp进行信号识别，所述绝对值标准偏差σaΔp＞0为QPSK信号，所述绝对值标准偏差σaΔp＝0为OQPSK信号。
可选的，所述瞬时相位反正切采用CORDIC算法。
可选的，所述基于所述瞬时相位计算相邻码元的相位差包括对所述I/Q信号进行位同步，以等同于符号率rs的速率在最佳时刻定时恢复出码元后，计算出相应的瞬时相位，并得到所述相邻码元的相位差
可选的，所述最佳时刻根据采用的位同步算法确定。
可选的，所述基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)为当载波同步时，计算所述中心归一化瞬时相位差ΔφNL(i)。
可选的，所述I/Q信号是将待识别的信号经过A/D采样变成数字信号，对所述数字信号进行正交数字下频所得。
可选的，所述将待识别的信号经过A/D采样变成数字信号包括将所述待识别的信号模拟下变频到A/D采样的范围，并满足过采样定律或带通采样定律。
本发明采取以上技术方案，与现有技术相比，具有以下优点：
利用QPSK信号和OQPSKI信号相位跳变的不同，简化QPSK信号和OQPSK信号的识别步骤，并且本发明的技术方案操作简单，易于操作实现，实现速度快，结果精确。
附图说明
图1是本发明实施例提供的QPSK信号和OQPSK信号的识别方法的流程图；
图2是本发明实施例提供的QPSK信号的相位跳转图；
图3是本发明实施例提供的OQPSK信号的相位跳转图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。
在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
图1是本发明实施例提供的QPSK信号和OQPSK信号的识别方法的流程图，下面结合图1(图1中仅以文字描述，未显示数学标识)详细说明。
步骤S1，从下变频后的I/Q信号提取出瞬时相位
步骤S2，基于所述瞬时相位计算相邻码元的相位差
步骤S3，基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)；
步骤S4，基于所述中心归一化瞬时相位差ΔφNL(i)计算ΔφNL(i)的绝对值标准偏差σaΔp，根据σaΔp进行信号识别，σaΔp＞0为QPSK信号，σaΔp＝0为OQPSK信号。
所述步骤S1中，所述I/Q信号是将待识别的信号经过A/D采样变成数字信号，对所述数字信号进行正交数字下频所得。所述将待识别的信号经过A/D采样变成数字信号包括将所述待识别的信号模拟下变频到A/D采样的范围，并满足过采样定律或带通采样定律。
所述步骤S1中，所述瞬时相位其中I(i)、Q(i)为正交下变频后的I路信号和Q路信号，0＜i＜Ns，Ns为采样点数；所述瞬时相位反正切采用CORDIC算法。
所述步骤S2中，所述相邻码元的相位差 Nr为符号数，Nr＝Ns/(fs/rs)，rs为符号率，fs为采样率；所述基于所述瞬时相位计算相邻码元的相位差包括对所述I/Q信号进行位同步，以等同于符号率rs(信号源的码元速率，简称码元速率)的速率在最佳时刻定时恢复出码元后(所述最佳时刻根据采用的位同步算法确定，例如，采用Gardner算法时，在经过成形滤波后的码元的最大值时刻恢复出码元，具体参考Gardner算法或其他位同步算法)，计算出相应的瞬时相位，并得到所述相邻码元的相位差
所述步骤S3中，基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)：当载波同步时，计算中心归一化瞬时相位差其中，Nr为符号数，Nr＝Ns/(fs/rs)，rs为符号率，fs为采样率。
所述步骤S4中，基于所述中心归一化瞬时相位差ΔφNL(i)计算ΔφNL(i)的绝对值标准偏差σaΔp：其中， an(i)是经过定时恢复后信号的瞬时归一化幅度，ma是经过定时恢复后的信号的瞬时幅度的平均值，ta是判断弱信号段的一个幅度的判决门限电平，C是在集合{φNL(i)}中an(i)＞ta的an(i)信号的个数，0＜i＜Nr；根据σaΔp进行信号识别，σaΔp＞0为QPSK信号，σaΔp＝0为OQPSK信号。
下面举例详细说明：
目前的信号识别方法为数字调制信号识别算法(也称DMRAs算法)，DMRAs调制识别的原理如下：被识别的信号经过A/D采样(如果信号频率太大，硬件中需要有射频前端电路，首先把被识别的信号模拟下变频到A/D可采样的范围，能满足过采样定律或带通采样定律)变成数字信号，数字信号经过正交数字下变频后，形成正交I/Q信号，然后利用I/Q信号计算幅度、频率和相位瞬时特征参数，再以相位瞬时特征参数作为输入，利用DMRAs算法计算信号识别特征参数，最后利用信号识别特征参数进行信号识别。
利用I/Q信号计算幅度、频率及相位等瞬时特征参数的公式如下：
其中a(i)为离散瞬时幅度，为瞬时相位，f(i)为瞬时频率。
信号识别特征参数是利用瞬时特征参数计算的，经典的DMRAs信号识别特征参数主要有以下几个：
(1)零中心归一化瞬时幅度之谱密度的最大值γmax；
(2)零中心归一化非弱信号段瞬时幅度的标准偏差σa；
(3)零中心归一化非弱信号段瞬时幅度绝对值的标准偏差σaa；
(4)零中心非弱信号段瞬时相位非线性分量的标准偏差σdp；
(5)零中心非弱信号段瞬时相位非线性分量的绝对值的标准偏差σap；
(6)零中心归一化非弱信号段瞬时频率绝对值的标准偏差σaf；
(7)参量φa的非线性分量的绝对值标准偏差σaap；
(8)参量fa的非线性分量的绝对值标准偏差σaaf；
(9)谱对称性P；
(10)零中心归一化瞬时幅度的紧致性(四阶矩)
(11)零中心归一化瞬时频率的紧致性(四阶矩)
由于QPSK信号和OQPSK信号的频谱图和星座图相同，用上述11个公式无法进行信号识别。
本发明的技术方案是利用QPSK信号和OQPSK信号相位跳变的不同(如图2和图3所示，对σap进行改进，完成QPSK信号和OQPSK信号之间的信号识别。
DMRAs的σap公式如下：
其中Ns为采样点数an(i)是经过定时恢复后信号的瞬时归一化幅度，ma是经过定时恢复后的信号的瞬时幅度的平均值；ta是判断弱信号段的一个幅度的判决门限电平，C是在集合{φNL(i)}中an(i)＞ta的an(i)信号的个数；0＜i＜Nr，Nr为符号数；φNL(i)是经中心化处理后的瞬时相位的非线性分量(也称中心归一化瞬时相位)。在载波同步时，φNL(i)计算如下：
上两个公式中，为瞬时相位，ta是判断弱信号段的一个幅度判决门限电平，低于此电平门限时，信号的相位对噪声很敏感，会造成较大的相位失真。C是在数据集合{φNL(i)}中an(i)＞ta的信号的个数，即属于非弱信号值的个数。
σap在信号识别中用来区分4PSK和2PSK信号的。对于2PSK信号，因为其只有两个相位值，分别取0和π，其中心归一化瞬时相位φNL(i)的绝对值为常数，即2PSK信号不包含绝对相位信息(σap＜t(σap))。而对于4PSK，因为其瞬时相位有四个值，φNL(i)取绝对值后仍带有相位信息，即本身包含有绝对和直接信号的相位信息(σap＞t(σap))。因此，该特征值可被用于区分具有绝对相位信息的信号(4PSK)和不含绝对相位信息的信号(2PSK)。但是σap不能直接区分QPSK信号和OQPSK信号。
QPSK相位跳转方式有{0°,90°,180°,-90°}共4种情况，当产生180度的跳转时，会使信号的旁瓣增大，影响频谱性能。为了解决这个问题，改进后的OQPSK调制方式将I/Q两路信号产生一个码元的时间间隔，有效避免了两路码元信号同时改变，从而使相位跳转只有{90°,-90°}，解决了180°的大幅跳转问题，图2和图3给出了QPSK和OQPSK的相位跳转图。
由于QPSK相邻码元相位有3个跳变值，最大跳变为±π，而OQPSK的相邻码元相位有2个跳变值，最大跳变为±π/2。根据QPSK和OQPSK的这个特征，可以参考DMRAs算法中经典特征参数σap，改进后用于QPSK和OQPSK之间的识别，具体方法如下：
从下变频后的I/Q信号提取出瞬时相位
用于计算σap，反正切可以用CORDIC算法实现。为了计算相邻码元相位的跳变值(也就是相位差)，还要对I/Q信号进行位同步，以等同于符号率rs的固定速率在最佳时刻定时恢复出码元后，计算出相应的从而得到相邻码元的相位差即为：
将φNL(i)中的瞬时相位替换为相邻码元相位的跳变值可以得到中心归一化瞬时相位差ΔφNL(i)：
Nr为符号数，Nr＝Ns/(fs/rs)，rs为符号率，fs为采样率，Ns为采样点数。将σap中的φNL(i)用ΔφNL(i)替换，可以得到：
其中an(i)是经过定时恢复后信号的瞬时归一化幅度，ma是经过定时恢复后的信号的瞬时幅度的平均值。
对于OQPSK信号，因为只有±π/2，其中心归一化瞬时相位差ΔφNL(i)的绝对值为常数，即OQPSK信号不包含绝对相位差信息，σaΔp＝0。而对于QPSK，有四个值，ΔφNL(i)取绝对值后仍带有相位信息，即本身包含有绝对和直接信号的相位信息，σaΔp＞0。因此，利用该特征值可用于区分QPSK和OQPSK信号。
本发明采取以上技术方案，与现有技术相比，具有以下优点：
利用QPSK信号和OQPSK信号相位跳变的不同，简化QPSK信号和OQPSK信号的识别步骤，并且本发明的技术方案操作简单，易于操作实现，实现速度快，结果精确。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种QPSK信号和OQPSK信号的识别方法，其特征在于，包括：
从下变频后的I/Q信号提取出瞬时相位
其中I(i)、Q(i)为正交下变频后的I路信号和Q路信号，0＜i＜Ns，Ns为采样点数；
基于所述瞬时相位计算相邻码元的相位差 0＜i＜Nr，Nr为符号数，Nr＝Ns/(fs/rs)，rs为符号率，fs为采样率；
基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)：0＜i＜Nr；
基于所述中心归一化瞬时相位差ΔφNL(i)计算所述中心归一化瞬时相位差ΔφNL(i)的绝对值标准偏差σaΔp：
其中，
an(i)是经过定时恢复后信号的瞬时归一化幅度，ma是经过定时恢复后的信号的瞬时幅度的平均值，ta是判断弱信号段的一个幅度的判决门限电平，C是在数据集合{φNL(i)}中an(i)＞ta的an(i)信号的个数，0＜i＜Nr，a(i)为离散瞬时幅度；
根据所述绝对值标准偏差σaΔp进行信号识别，所述绝对值标准偏差σaΔp＞0为QPSK信号，所述绝对值标准偏差σaΔp＝0为OQPSK信号。
2.如权利要求1所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述瞬时相位反正切采用CORDIC算法。
3.如权利要求1所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述基于所述瞬时相位计算相邻码元的相位差包括对所述I/Q信号进行位同步，以等同于符号率rs的速率在最佳时刻定时恢复出码元后，计算出相应的瞬时相位，并得到所述相邻码元的相位差
4.如权利要求3所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述最佳时刻根据采用的位同步算法确定。
5.如权利要求1所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述基于所述相邻码元的相位差计算中心归一化瞬时相位差ΔφNL(i)为当载波同步时，计算所述中心归一化瞬时相位差ΔφNL(i)。
6.如权利要求1所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述I/Q信号是将待识别的信号经过A/D采样变成数字信号，对所述数字信号进行正交数字下变频所得。
7.如权利要求6所述的QPSK信号和OQPSK信号的识别方法，其特征在于，所述将待识别的信号经过A/D采样变成数字信号包括将所述待识别的信号模拟下变频到A/D采样的范围，并满足过采样定律或带通采样定律。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410853901.5A CN104639268B (zh) | 2014-12-31 | 2014-12-31 | Qpsk信号和oqpsk信号的识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410853901.5A CN104639268B (zh) | 2014-12-31 | 2014-12-31 | Qpsk信号和oqpsk信号的识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104639268A CN104639268A (zh) | 2015-05-20 |
| CN104639268B true CN104639268B (zh) | 2017-04-12 |
Family
ID=53217634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410853901.5A Active CN104639268B (zh) | 2014-12-31 | 2014-12-31 | Qpsk信号和oqpsk信号的识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104639268B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1395400A (zh) * | 2002-07-05 | 2003-02-05 | 清华大学 | 无需先验知识的常用数字信号调制方式的自动识别方法 |
| CN104022837A (zh) * | 2014-05-28 | 2014-09-03 | 广州海格通信集团股份有限公司 | 一种基于特征提取的通信信号调制方式进行识别的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7751503B2 (en) * | 2006-04-27 | 2010-07-06 | Jennic Limited | Method for acquiring timing and carrier synchronization of offset-QPSK modulated signals |
-
2014
- 2014-12-31 CN CN201410853901.5A patent/CN104639268B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1395400A (zh) * | 2002-07-05 | 2003-02-05 | 清华大学 | 无需先验知识的常用数字信号调制方式的自动识别方法 |
| CN104022837A (zh) * | 2014-05-28 | 2014-09-03 | 广州海格通信集团股份有限公司 | 一种基于特征提取的通信信号调制方式进行识别的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104639268A (zh) | 2015-05-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jajoo et al. | Blind signal modulation recognition through clustering analysis of constellation signature | |
| Jajoo et al. | Blind signal PSK/QAM recognition using clustering analysis of constellation signature in flat fading channel | |
| CN103780462B (zh) | 基于高阶累积量和谱特征的卫星通信信号调制识别方法 | |
| Headley et al. | Asynchronous classification of digital amplitude-phase modulated signals in flat-fading channels | |
| US8275077B1 (en) | Coherent demodulation of ais-GMSK signals in co-channel | |
| CN111711589B (zh) | 一种基于椭圆球面波信号的连续相位调制解调方法 | |
| KR101722505B1 (ko) | 입력 신호의 변조 형태 인식 방법 및 장치 | |
| Spooner et al. | Modulation recognition using second-and higher-order cyclostationarity | |
| KR101828790B1 (ko) | 주파수 편이 변조 신호의 수신 방법 및 장치 | |
| CN111935046B (zh) | 一种低复杂度的频移键控信号符号率估计方法 | |
| Majhi et al. | Novel blind modulation classification of circular and linearly modulated signals using cyclic cumulants | |
| EP3188428A2 (en) | Ultra low power wideband asynchronous binary phase shift keying demodulation circuit using first order sideband filters aligned at zero degree phase | |
| Flohberger et al. | Modulation classifier for signals used in satellite communications | |
| Bari et al. | Simple features for separating CPFSK from QAM and PSK modulations | |
| Haq et al. | Recognition of digital modulated signals based on statistical parameters | |
| Ghauri | KNN based classification of digital modulated signals | |
| CN104507106B (zh) | 8psk信号和pi/4-dqpsk信号的识别方法 | |
| CN104639268B (zh) | Qpsk信号和oqpsk信号的识别方法 | |
| CN108270703A (zh) | 一种通信信号数字调制类型识别方法 | |
| Bari et al. | Robust recognition of linear and nonlinear digital modulations of RRC pulse trains | |
| Yadav et al. | Modulation scheme detection of blind signal using constellation graphical representation | |
| Jajoo et al. | Blind signal digital modulation classification through k-medoids clustering | |
| US20180191540A1 (en) | Method and device for receiving a signal the phase or frequency of which is modulated by a sequence of two-state symbols | |
| Hatoum et al. | GENERALIZED WAVELET-BASED SYMBOL RATE ESTIMATION FOR LINEAR SINGLECARRIER MODULATION IN BLIND ENVIRONMENT | |
| KR102097366B1 (ko) | 구간합 궤환을 이용하는 반송파 주파수 오프셋 추정 방법 및 그 수신 장치 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |