CN104052702B - 一种复杂噪声下数字调制信号的识别方法 - Google Patents
一种复杂噪声下数字调制信号的识别方法 Download PDFInfo
- Publication number
- CN104052702B CN104052702B CN201410279349.3A CN201410279349A CN104052702B CN 104052702 B CN104052702 B CN 104052702B CN 201410279349 A CN201410279349 A CN 201410279349A CN 104052702 B CN104052702 B CN 104052702B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- signal
- sigma
- component
- 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
Landscapes
- Error Detection And Correction (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
一种复杂噪声下数字调制信号的识别方法,所述方法包括以下步骤:对接收到的信号的观测向量x做分数低阶快速独立成分分析,将接收到的信号分离为发送信号和Alpha稳定分布噪声;对分离出的发送信号s(n)做基于插值的总体平均LMD,将发送信号分解成多个分量;提取第一个分类特征:分量的分段瞬时频率标准差σf和设定相应的门限δ1;提取第二个分类特征:分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3;利用门限δ1将信号集{MSK、2ASK、QPSK、16QAM}分为{MSK}和{2ASK、QPSK、16QAM}两类,利用门限δ2和δ3将信号集{2ASK、QPSK、16QAM}中的信号识别出来。本发明可以对Alpha稳定分布噪声下的信号在低信噪比环境下具有较好的识别性能和良好的稳健性。
Description
技术领域
本发明属于通信技术领域,具体涉及一种复杂噪声下数字调制信号的方法。可用于Alpha稳定分布噪声下MSK信号、2ASK信号、QPSK信号和16QAM信号的识别。
背景技术
数字调制信号的识别是在接收方未知发送方信号调制方式的前提下,从给定的信号备择集中确定发送信号的调制方式,从而为后续的解调工作提供所需信息。信号的调制识别虽然最早起源于军事侦察,但随着识别技术、软件无线电等技术的发展,调制识别技术在民用通信领域也得到了广泛的应用。为了便于理论的分析,传统的信号处理中通常假设噪声为高斯模型。然而,实验研究表明,大部分无线电波信道中的人为电磁噪声、自然噪声以及两者的联合噪声显著偏离高斯分布,这些噪声具有显著的尖峰脉冲特性和比较厚的概率密度函数拖尾。研究发现可以采用Alpha稳定分布模型来描述这些脉冲性很强的复杂噪声,近年来很多学者将复杂噪声建模为Alpha稳定分布来进行研究。因此,研究Alpha稳定分布噪声环境下的数字信号的调制识别技术具有实际的工程意义和很好的应用价值。
目前,一些学者已对Alpha稳定分布噪声环境下的数字调制信号的识别进行了研究,但研究还是较少。Fanggang Wang等人利用K-S检测算法对4QAM,16QAM和64QAM信号在Alpha稳定分布噪声下进行识别。但该方法在低信噪比条件下,识别性能不理想(FanggangWang,Xiaodong Wang.Fast and Robust Modulation Classifi cation via Kolmogorov-Smirnov Test[J].IEEE TRANSACTIONS ON COMMUNICATIONS,2010,58(8):2324-2332.)。杨伟超等人以信号相位的分形盒维数作为识别特征,对Alpha噪声下信号进行识别。但该方法仅能在一定混合信噪比范围内适用且识别性能较差(杨伟超,赵春晖,成宝芝.Alpha稳定分布噪声下的通信信号识别[J].应用科学学报,2010,28(2):111-114.)。赵春晖和刘明骞等人利用含有Alpha噪声的信号的分数低阶循环谱的相干系数构造特征实现信号的系别。但该方法中循环谱指数b的设定因缺少智能方法而导致普适性较差(赵春晖,杨伟超,杜宇.采用分数低阶循环谱相干系数的调制识别[J].应用科学学报,2011,29(6):565-570.刘明骞,李兵兵,曹超凤,李钊.认知无线电中非高斯噪声下数字调制信号识别方法[J].通信学报,2014,35(1):82-88.)。赵春晖等人利用广义二阶循环谱中构造特征进行识别。但该方法计算复杂度高(赵春晖,杨伟超,马爽.基于广义二阶循环统计量的通信信号调制识别算法研究[J].通信学报,2011,32(1):144-150.)。刘明骞等人提取广义分数阶傅里叶变换的零中心归一化瞬时幅度谱密度的最大值和分数低阶Wigner-Ville分布幅度的最大值作为特征进行识别。但该方法在低信噪比条件下,识别性能不理想(刘明骞,李兵兵,曹超凤.非高斯噪声下数字调制信号识别方法[J].电子与信息学报,2013,35(1):85-91.)。赵春晖等人利用广义四次方谱进行了调制识别的研究,但该方法在低信噪比条件下QPSK的识别性能很差(赵春晖,杨伟超.Alpha稳定分布噪声下MPSK信号调制识别算法研究[J].沈阳大学学报,2013,25(1):10-14.)。
发明内容
本发明的目的是克服上述已有技术的不足,提出了一种有效识别复杂噪声下数字调制信号的方法,以提高在复杂噪声-Alpha稳定分布噪声环境下数字调制信号的识别率。本发明选取MSK、2ASK、QPSK和16QAM这4种数字调制信号作为待识别信号集。
为了实现上述目的,本发明采用的技术方案如下:
一种复杂噪声下数字调制信号的识别方法,所述方法包括以下步骤:
S1对接收到的信号的观测向量x做分数低阶快速独立成分分析,将接收到的信号分离为发送信号和Alpha稳定分布噪声;
S2对分离出的发送信号s(n)做基于插值的总体平均LMD,将发送信号分解成多个分量;
S3提取第一个分类特征:分量的分段瞬时频率标准差σf和设定相应的门限δ1;
S4提取第二个分类特征:分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3;
S5利用门限δ1将信号集{MSK、2ASK、QPSK、16QAM}分为{MSK}和{2ASK、QPSK、16QAM}两类,利用门限δ2和δ3将信号集{2ASK、QPSK、16QAM}中的信号识别出来。
在上述技术方案的基础上,所述接收信号的观测向量x的分数低阶快速独立成分分析按以下进行:
1)将观测向量x减去均值进行中心化,利用主要分量分析对中心化后的观测向量进行分数低阶预白化处理,得到白化矩阵v,其中,PCA中用到的分数低阶相关矩阵定义为
其中,xi(n)和xj(n)分别为第i路和第j路观测信号,且1≤i,j≤2,*代表取共轭,对分数低阶相关矩阵进行特征值分解可以得到特征值矩阵D与特征向量矩阵V,由特征值矩阵和特征向量矩阵可得预白化矩阵M=D-1/2VT,利用预白化矩阵将观测数据x向信号子空间投影,得到白化矩阵v=Mx;
2)对权向量w0进行随机初始化,同时初始化序列号k使k=1;
3)进行权值向量迭代
其中,和为一阶统计量,g(·)是对比函数,η=Ε[wTxg(wTx)],且在权值的迭代过程中参数λ是可变的,且λ<1。
4)利用范数归一化权值向量wk,wk=wk/||wk||;
5)若充分接近于1,算法得到收敛,wk就是最终的解混矩阵,执行6),否则重复步骤3)与步骤4);
6)Y=wkx,其中x为观测信号,Y为分离后的二维信号。
在上述技术方案的基础上,所述分离出的发送信号s(n)的基于插值的总体平均LMD按以下进行:
1)给发送信号s(n)加上一组白噪声v(n),得到一个总体S(n):
S(n)=s(n)+v(n)
2)确定总体S(n)所有的局部极值点ni,计算相邻两个极值点ni和ni+1的平均值mi,即
利用三次样条插值确定局部均值函数m11(n);
3)采用局部极值点ni计算包络估计值ai
利用三次样条插值确定包络估计函数a11(n);
4)从总体S(n)中分离分离出局部均值函数m11(n),得到
h11(n)=S(n)-m11(n)
5)用h11(n)除以包络估计函数a11(n),以对h11(n)进行解调,得到
其中,
h11(n)=s(n)-m11(n)
h12(n)=s11(n)-m12(n)
h1n(n)=s1(n-1)(n)-m1n(n)
式中:
迭代终止的条件为:
1-Δ≤a1n(n)≤1+Δ
6)把迭代过程中产生的所有包络估计函数相乘便可以得到包络信号,其包括,
7)将包络信号a1(n)和纯调频信号s1n(n)相乘便可以得到原始信号的第1个分量
PF1(n)=a1(n)s1n(n)
它包含了原始信号中最高的频率成分,是一个单分量的调幅-调频信号,其瞬时幅值就是包络信号a1(n),其瞬时频率f1(n)则可由纯调频信号s1n(n)求出,即:
8)将第1个PF分量PF1从发送信号s(n)中分离出来,得到一个新的信号u1(n),将u1(n)作为原始数据重复步骤1)至步骤7),并循环k次,直到uk(n)为一个单调函数;
u1(n)=s(n)-PF1(n)
u2(n)=u1(n)-PF2(n)
uk(n)=uk-1(n)-PFk(n)
至此,将信号s(n)分解为k个PF分量和一个单调函数uk(n)之和,即
9)给发送信号s(n)加入不同的白噪声ni(n),重复步骤1)-8):
Si(n)=s(n)+ni(n)
分解后得到的各个总体的PF分量组为PFip
10)取各个组PF的均值作为最终的PF组:
其中,N表示表示添加噪声的次数,最后发送信号的LMD分解可表示为:
在上述技术方案的基础上,所述分量的分段瞬时频率标准差σf和设定相应的门限δ1按以下进行:
分量的分段频率标准差σf定义为:对于信号的分量中的纯调频信号s1n(n),计算其瞬时频率f1(n),将得到的瞬时频率序列分成M段,每段n个点,然后计算每段的瞬时频率的均值,最后求M段频率的标准差,即
其中,为第i段瞬时频率的均值,f1i(j)为分量的第i段数据的第j个点的瞬时频率值;为分量的瞬时频率的均值;
门限设定如下:
其中min(σf-MSK)表示MSK信号的特征值σf的最小值,max(σf-2ASK,σf-QPSK,σf-16QAM)为2ASK、QPSK和16QAM信号的特征值σf的最大值。
在上述技术方案的基础上,所述分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3按以下进行:
的分段瞬时幅值标准差σa定义为:对于信号的PF1分量,将其a1(n)分成M段,每段n个点,然后计算每段的瞬时幅值的均值,最后求M段幅值的标准差,即
其中,为第i段瞬时幅值的均值,a1i(j)为分量的第i段数据的第j个点的瞬时幅值;为PF1分量的瞬时幅值的均值;
门限设定如下:
其中min(σf-2ASK)表示2ASK信号的特征值σa的最小值,max(σa-16QAM)和min(σa-QPSK)分别为16QAM信号的特征值σa的最大值和最小值,max(σa-QPSK)为QPSK信号的特征值σa的最大值。
本发明有益效果在于:
1、本发明可以识别多种不同调制类型的Alpha稳定分布噪声下的信号;
2、本发明在低信噪比环境下具有较好的识别性能;
3、本发明对Alpha稳定分布噪声的特征指数α有良好的稳健性;
4、在相同的仿真实验环境和相同的码元速率、载波频率、采样频率、采样点数和信噪比等信号参数设置条件下,本发明比现有的方法具有更高的识别率。
附图说明
图1为本发明的流程示意图;
图2为本发明未考虑滚降滤波时在不同混合信噪比下信号的正确识别率;
图3中为本发明未考虑滚降滤波时在不同特征指数α下信号的正确识别率;
图4中为本发明考虑滚降滤波时在不同混合信噪比下信号的正确识别率;
图5中为本发明在相同的仿真实验环境和信号参数设置下,本发明与传统的基于分数低阶循环谱相干系数方法和基于FRFT和FLOWVD方法的性能对比图。
具体实施方式
本发明的具体实现步骤如下:
如图1所示,本发明为一种复杂噪声下数字调制信号的识别方法,所述方法包括以下步骤:
S1对接收到的信号的观测向量x做分数低阶快速独立成分分析,将接收到的信号分离为发送信号和Alpha稳定分布噪声;
需要说明的是,所述接收信号的观测向量x的分数低阶快速独立成分分析按以下进行:
1)将观测向量x减去均值进行中心化,利用主要分量分析对中心化后的观测向量进行分数低阶预白化处理,得到白化矩阵v,其中,PCA中用到的分数低阶相关矩阵定义为
其中,xi(n)和xj(n)分别为第i路和第j路观测信号,且1≤i,j≤2,*代表取共轭,对分数低阶相关矩阵进行特征值分解可以得到特征值矩阵D与特征向量矩阵V,由特征值矩阵和特征向量矩阵可得预白化矩阵M=D-1/2VT,利用预白化矩阵将观测数据x向信号子空间投影,得到白化矩阵v=Mx;
2)对权向量w0进行随机初始化,同时初始化序列号k使k=1;
3)进行权值向量迭代
其中,和为一阶统计量,g(·)是对比函数,η=Ε[wTxg(wTx)],且在权值的迭代过程中参数λ是可变的,且λ<1。
4)利用范数归一化权值向量wk,wk=wk/||wk||;
5)若充分接近于1,算法得到收敛,wk就是最终的解混矩阵,执行6),否则重复步骤3)与步骤4);
6)Y=wkx,其中x为观测信号,Y为分离后的二维信号。
S2对分离出的发送信号s(n)做基于插值的总体平均LMD,将发送信号分解成多个分量;
需要说明的是,所述分离出的发送信号s(n)的基于插值的总体平均LMD按以下进行:
1)给发送信号s(n)加上一组白噪声v(n),得到一个总体S(n):
S(n)=s(n)+v(n)
2)确定总体S(n)所有的局部极值点ni,计算相邻两个极值点ni和ni+1的平均值mi,即
利用三次样条插值确定局部均值函数m11(n);
3)采用局部极值点ni计算包络估计值ai
利用三次样条插值确定包络估计函数a11(n);
4)从总体S(n)中分离分离出局部均值函数m11(n),得到
h11(n)=S(n)-m11(n)
5)用h11(n)除以包络估计函数a11(n),以对h11(n)进行解调,得到
其中,
h11(n)=s(n)-m11(n)
h12(n)=s11(n)-m12(n)
h1n(n)=s1(n-1)(n)-m1n(n)
式中:
迭代终止的条件为:(建议此处直接采用当满足1-Δ≤a1n(n)≤1+Δ时的状态,因此,建议改变迭代终止的条件)
1-Δ≤a1n(n)≤1+Δ
6)把迭代过程中产生的所有包络估计函数相乘便可以得到包络信号,其包括,
7)将包络信号a1(n)和纯调频信号s1n(n)相乘便可以得到原始信号的第1个分量
PF1(n)=a1(n)s1n(n)
它包含了原始信号中最高的频率成分,是一个单分量的调幅-调频信号,其瞬时幅值就是包络信号a1(n),其瞬时频率f1(n)则可由纯调频信号s1n(n)求出,即:
8)将第1个PF分量PF1从发送信号s(n)中分离出来,得到一个新的信号u1(n),将u1(n)作为原始数据重复步骤1)至步骤7),并循环k次,直到uk(n)为一个单调函数;
u1(n)=s(n)-PF1(n)
u2(n)=u1(n)-PF2(n)
uk(n)=uk-1(n)-PFk(n)
至此,将信号s(n)分解为k个PF分量和一个单调函数uk(n)之和,即
9)给发送信号s(n)加入不同的白噪声ni(n),重复步骤1)-8):
Si(n)=s(n)+ni(n)
分解后得到的各个总体的PF分量组为PFip
10)取各个组PF的均值作为最终的PF组:
其中,N表示表示添加噪声的次数,最后发送信号的LMD分解可表示为:
S3提取第一个分类特征:分量的分段瞬时频率标准差σf和设定相应的门限δ1;
需要说明的是,所述分量的分段瞬时频率标准差σf和设定相应的门限δ1按以下进行:
分量的分段频率标准差σf定义为:对于信号的分量中的纯调频信号s1n(n),计算其瞬时频率f1(n),将得到的瞬时频率序列分成M段,每段n个点,然后计算每段的瞬时频率的均值,最后求M段频率的标准差,即
其中,为第i段瞬时频率的均值,f1i(j)为分量的第i段数据的第j个点的瞬时频率值;为分量的瞬时频率的均值;
门限设定如下:
其中min(σf-MSK)表示MSK信号的特征值σf的最小值,max(σf-2ASK,σf-QPSK,σf-16QAM)为2ASK、QPSK和16QAM信号的特征值σf的最大值。
S4提取第二个分类特征:分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3;
需要说明的是,所述分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3按以下进行:
的分段瞬时幅值标准差σa定义为:对于信号的PF1分量,将其a1(n)分成M段,每段n个点,然后计算每段的瞬时幅值的均值,最后求M段幅值的标准差,即
其中,为第i段瞬时幅值的均值,a1i(j)为分量的第i段数据的第j个点的瞬时幅值;为PF1分量的瞬时幅值的均值;
门限设定如下:
其中min(σf-2ASK)表示2ASK信号的特征值σa的最小值,max(σa-16QAM)和min(σa-QPSK)分别为16QAM信号的特征值σa的最大值和最小值,max(σa-QPSK)为QPSK信号的特征值σa的最大值。
S5利用门限δ1将信号集{MSK、2ASK、QPSK、16QAM}分为{MSK}和{2ASK、QPSK、16QAM}两类,利用门限δ2和δ3将信号集{2ASK、QPSK、16QAM}中的信号识别出来。
本发明的复杂噪声为Alpha稳定分布噪声。为了从不同侧面评估方法的性能,下面的仿真实验采用信号的类型为MSK信号、2ASK信号、QPSK信号和16QAM信号,并进行500次Monte Carlo实验。识别的评估标准为正确识别率。
为了测试混合信噪比对未考虑升余弦滚降滤波的信号识别性能的影响,BPSK信号、QPSK信号、8PSK信号和MSK信号的参数设置如下:载波频率为fc=30kHz,码元速率为rb=40000Baud,采样频率fs=120KHz,数据长度为2000。Alpha稳定噪声的特征指数为1.5。如图2所示,当固定其他参数时,各个信号的正确识别率随着混合信噪比的增大而相应的增大。当混合信噪比大于-8dB时,MSK和2ASK信号的正确识别率接近100%,QPSK和16QAM信号的正确识别率大于95%。由此说明本发明发放在低信噪比Alpha稳定分布噪声环境下具有较好的识别性能。
为了测试不同的特征指数α对识别性能的影响,BPSK信号、QPSK信号、8PSK信号和MSK信号的参数设置如下:未考虑升余弦滚降滤波的影响,载波频率为fc=3.0kHz,码元速率为rb=1200Baud,采样频率fs=19.2KHz,数据长度为3000。Alpha噪声的特征指数α的变化范围是1到2,间隔为0.1。如图3所示,MSK、2ASK、QPSK和16QAM信号的正确识别率均随着噪声的特征指数α的逐渐增大而有所提升。当α≥1.2时,四种数字调制信号的正确识别率均大于90%,当α≥1.5时,信号的识别性能受特征指数α的影响较小。并且在高斯情况下,即α=2时,本发明方法也具有非常好的识别性能。
为了测试混合信噪比对升余弦滚降系数为0.35的信号识别性能的影响,BPSK信号、QPSK信号、8PSK信号和MSK信号的参数设置如下:载波频率为fc=30kHz,码元速率为rb=40000Baud,采样频率fs=120KHz,数据长度为2000。Alpha稳定噪声的特征指数为1.5。如图4所示,在考虑滚降滤波的情况下,4种信号的正确识别率均随着混合信噪比的增大而相应的提高。当混合信噪比大于0dB时,MSK信号和2ASK信号的正确识别率接近100%,QPSK信号的正确识别率大于96%,16QAM信号的正确识别率大于92%。由此说明本发明方法具有较强的鲁棒性。
为了进一步说明本发明的优越性,在相同的仿真实验环境和信号参数设置下,本发明方法与刘明骞等人的基于分数低阶循环谱的识别方法和刘明骞等人的基于分数阶傅里叶变换(FRFT)和分数低阶Wigner-Ville(FLOWVD)方法进行对比试验。如图5所示,本发明方法的识别性能均优于传统的基于分数低阶循环谱方法和基于FRFT和FLOWVD的识别方法。
对于本领域的技术人员来说,可根据以上描述的技术方案以及构思,做出其它各种相应的改变以及变形,而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。
Claims (1)
1.一种复杂噪声下数字调制信号的识别方法,其特征在于:所述方法包括以下步骤:
S1对接收到的信号的观测向量x做分数低阶快速独立成分分析,将接收到的信号分离为发送信号和Alpha稳定分布噪声;所述接收信号的观测向量x的分数低阶快速独立成分分析按以下进行:
1)将观测向量x减去均值进行中心化,利用主要分量分析对中心化后的观测向量进行分数低阶预白化处理,得到白化矩阵v,其中,独立成分分析(PCA)中用到的分数低阶相关矩阵定义为
<mrow>
<mo><</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>></mo>
<mi>p</mi>
</msub>
<mo>=</mo>
<mi>E</mi>
<mo>{</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mrow>
<mi>p</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<msubsup>
<mi>x</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
<mo>}</mo>
<mo>,</mo>
<mn>0</mn>
<mo><</mo>
<mi>p</mi>
<mo><</mo>
<mi>&alpha;</mi>
<mo>&le;</mo>
<mn>2</mn>
</mrow>
其中,xi(n)和xj(n)分别为第i路和第j路观测信号,且1≤i,j≤2,*代表取共轭,α为alpha稳定分布噪声的特征指数,对分数低阶相关矩阵进行特征值分解可以得到特征值矩阵D与特征向量矩阵V,由特征值矩阵和特征向量矩阵可得预白化矩阵M=D-1/2VT,利用预白化矩阵将观测数据x向信号子空间投影,得到白化矩阵v=Mx;
2)对权向量w0进行随机初始化,同时初始化序列号k使k=1;
3)进行权值向量迭代
<mrow>
<msub>
<mi>w</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>v</mi>
<mi>g</mi>
<mo>(</mo>
<mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<mi>v</mi>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<msub>
<mi>&eta;w</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<msup>
<mi>g</mi>
<mo>&prime;</mo>
</msup>
<mo>(</mo>
<mrow>
<msubsup>
<mi>w</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<mi>v</mi>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,和为一阶统计量,g(·)是对比函数,且在权值的迭代过程中参数λ是可变的,且λ<1;
4)利用范数归一化权值向量wk,wk=wk/||wk||;
5)若充分接近于1,算法得到收敛,wk就是最终的解混矩阵,执行6),否则重复步骤3)与步骤4);
6)Y=wkx,其中x为观测信号,Y为分离后的二维信号;
S2对分离出的发送信号s(n)做基于插值的总体平均局部均值分解(LMD);
将发送信号s(n)分解为k个PF分量和一个单调函数uk(n)之和,给发送信号s(n)加入不同的白噪声,分解后得到的各个总体的PF分量组为PFip,
取各个分量组PFip的均值作为最终的PF组:
<mrow>
<mover>
<mrow>
<msub>
<mi>PF</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>PF</mi>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mover>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mo>~</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
其中,N表示表示添加噪声的次数,最后发送信号的LMD分解可表示为:
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<mrow>
<msub>
<mi>PF</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
S3提取第一个分类特征:分量的分段瞬时频率标准差σf和设定相应的门限δ1;所述分量的分段瞬时频率标准差σf和设定相应的门限δ1按以下进行:
分量的分段频率标准差σf定义为:对于信号的分量中的纯调频信号,计算其瞬时频率f1(n),将得到的瞬时频率序列分成M段,每段n个点,然后计算每段的瞬时频率的均值,最后求M段频率的标准差,即
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>~</mo>
</mover>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>E</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,为第i段瞬时频率的均值,f1i(j)为分量的第i段数据的第j个点的瞬时频率值;为分量的瞬时频率的均值;
门限设定如下:
<mrow>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>f</mi>
<mo>-</mo>
<mi>M</mi>
<mi>S</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>f</mi>
<mo>-</mo>
<mn>2</mn>
<mi>A</mi>
<mi>S</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>f</mi>
<mo>-</mo>
<mi>Q</mi>
<mi>P</mi>
<mi>S</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>f</mi>
<mo>-</mo>
<mn>16</mn>
<mi>Q</mi>
<mi>A</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
其中min(σf-MSK)表示MSK信号的特征值σf的最小值,max(σf-2ASK,σf-QPSK,σf-16QAM)为2ASK、QPSK和16QAM信号的特征值σf的最大值;
S4提取第二个分类特征:分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3;所述分量的分段瞬时幅值标准差σa和设定相应的门限δ2和δ3按以下进行:
的分段瞬时幅值标准差σa定义为:对于信号的分量,将其包络信号分成M段,每段n个点,然后计算每段的瞬时幅值的均值,最后求M段幅值的标准差,即
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>~</mo>
</mover>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>E</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,为第i段瞬时幅值的均值,a1i(j)为分量的第i段数据的第j个点的瞬时幅值;为PF1分量的瞬时幅值的均值;
门限设定如下:
<mrow>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mn>2</mn>
<mi>A</mi>
<mi>S</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mn>16</mn>
<mi>Q</mi>
<mi>A</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&delta;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mn>16</mn>
<mi>Q</mi>
<mi>A</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mi>Q</mi>
<mi>P</mi>
<mi>S</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
其中min(σa-2ASK)表示2ASK信号的特征值σa的最小值,max(σa-16QAM)和min(σa-16QAM)分别为16QAM信号的特征值σa的最大值和最小值,max(σa-QPSK)为QPSK信号的特征值σa的最大值;
S5利用门限δ1将信号集{MSK、2ASK、QPSK、16QAM}分为{MSK}和{2ASK、QPSK、16QAM}两类,利用门限δ2和δ3将信号集{2ASK、QPSK、16QAM}中的信号识别出来。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410279349.3A CN104052702B (zh) | 2014-06-20 | 2014-06-20 | 一种复杂噪声下数字调制信号的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410279349.3A CN104052702B (zh) | 2014-06-20 | 2014-06-20 | 一种复杂噪声下数字调制信号的识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104052702A CN104052702A (zh) | 2014-09-17 |
CN104052702B true CN104052702B (zh) | 2017-12-08 |
Family
ID=51505078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410279349.3A Active CN104052702B (zh) | 2014-06-20 | 2014-06-20 | 一种复杂噪声下数字调制信号的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104052702B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017178075A1 (en) * | 2016-04-15 | 2017-10-19 | Huawei Technologies Co., Ltd. | Method for generating pulse shapes with adjustable length, orthogonality and localization properties |
CN106067004B (zh) * | 2016-05-30 | 2019-07-12 | 西安电子科技大学 | 一种脉冲噪声下数字调制信号的识别方法 |
CN106597243B (zh) * | 2017-02-14 | 2018-12-07 | 吴笃贵 | 一种基于局部放电全息数据的概率特征参数提取方法 |
CN107158684B (zh) * | 2017-04-11 | 2019-02-15 | 石磊 | 一种用于大学生智能体育体能综合训练控制系统 |
CN107271978B (zh) * | 2017-05-18 | 2020-08-28 | 西安电子科技大学 | 一种多个异构卫星下基于Rao检测的目标探测方法 |
CN107332632B (zh) * | 2017-06-24 | 2020-10-09 | 西安电子科技大学 | MSK调制方式下Alpha稳定分布噪声参数的估计方法 |
CN107458596A (zh) * | 2017-08-16 | 2017-12-12 | 重庆华凤衣道文化创意有限公司 | 一种龙形太阳能与空气能复合飞行器 |
CN108354620A (zh) * | 2018-02-28 | 2018-08-03 | 罗兴梅 | 一种用于评估心理神经活动的生理指标检测系统 |
CN108492834A (zh) * | 2018-02-28 | 2018-09-04 | 湖南城市学院 | 一种声乐教学演唱技能检测系统 |
CN108712808A (zh) * | 2018-04-04 | 2018-10-26 | 湖南城市学院 | 一种智能环境艺术装饰灯光控制系统 |
CN108627621A (zh) * | 2018-04-25 | 2018-10-09 | 遵义市产品质量检验检测院 | 一种智能食品质量快速检测系统 |
CN108519745A (zh) * | 2018-05-03 | 2018-09-11 | 广东水利电力职业技术学院(广东省水利电力技工学校) | 一种运动控制器系统 |
CN109165545A (zh) * | 2018-07-05 | 2019-01-08 | 上海电机学院 | 一种新的基于分形维数的信号特征提取方法 |
CN110046676B (zh) * | 2019-04-26 | 2021-08-03 | 厦门大学 | 一种基于声暴露级的海洋工程信号自动提取方法 |
CN112631147B (zh) * | 2020-12-08 | 2023-05-02 | 国网四川省电力公司经济技术研究院 | 一种面向脉冲噪声环境的智能电网频率估计方法及系统 |
CN112769724B (zh) * | 2021-04-07 | 2021-07-02 | 广州慧睿思通科技股份有限公司 | 对讲机的信号类型的识别方法、装置、设备及存储介质 |
CN113518049B (zh) * | 2021-04-13 | 2024-04-26 | 江苏师范大学 | 一种基于分数低阶极坐标和深度学习的调制识别方法 |
CN115225438B (zh) * | 2022-07-07 | 2023-05-26 | 金陵科技学院 | 基于分段线性压缩量化的bpsk及qpsk信号调制识别方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326975A (zh) * | 2013-07-15 | 2013-09-25 | 西安电子科技大学 | 一种Alpha稳定分布噪声下数字调制信号识别方法 |
CN103457890A (zh) * | 2013-09-03 | 2013-12-18 | 西安电子科技大学 | 一种有效识别非高斯噪声下数字调制信号的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007027839A2 (en) * | 2005-09-01 | 2007-03-08 | University Of Florida Research Foundation, Inc. | Device and methods for enhanced matched filtering based on correntropy |
-
2014
- 2014-06-20 CN CN201410279349.3A patent/CN104052702B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326975A (zh) * | 2013-07-15 | 2013-09-25 | 西安电子科技大学 | 一种Alpha稳定分布噪声下数字调制信号识别方法 |
CN103457890A (zh) * | 2013-09-03 | 2013-12-18 | 西安电子科技大学 | 一种有效识别非高斯噪声下数字调制信号的方法 |
Non-Patent Citations (2)
Title |
---|
认知无线电中非高斯噪声下数字调制信号识别方法;刘明骞等;《通信学报》;20140131;全文 * |
非高斯噪声下数字调制信号识别方法;刘明骞等;《电子与信息学报》;20130131;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104052702A (zh) | 2014-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104052702B (zh) | 一种复杂噪声下数字调制信号的识别方法 | |
US10985955B2 (en) | Method for automatically identifying modulation mode for digital communication signal | |
CN106130942B (zh) | 一种基于循环谱的无线通信信号调制识别及参数估计方法 | |
CN104038454B (zh) | 一种Alpha稳定分布噪声下PSK信号的载波频率估计方法 | |
CN103780462B (zh) | 基于高阶累积量和谱特征的卫星通信信号调制识别方法 | |
CN107948107B (zh) | 一种基于联合特征的数字调制信号分类方法 | |
CN103326975B (zh) | 一种Alpha稳定分布噪声下数字调制信号识别方法 | |
CN104301056B (zh) | 一种基于信号特征分析的频谱监测方法 | |
CN103457890A (zh) | 一种有效识别非高斯噪声下数字调制信号的方法 | |
CN105119862A (zh) | 一种信号调制类型识别方法和系统 | |
CN105678273A (zh) | 射频指纹识别技术瞬态信号的起始点检测算法 | |
CN102882820A (zh) | 认知无线电中非高斯噪声下数字调制信号识别方法 | |
CN103997475B (zh) | 一种识别Alpha稳定分布噪声下数字调制信号的方法 | |
CN111585922A (zh) | 一种基于卷积神经网络的调制方式识别方法 | |
CN106452673B (zh) | 多传感器接收中基于特征向量的最大似然调制识别方法 | |
CN101764786A (zh) | 基于聚类算法的mqam信号识别方法 | |
Ahmadi | Using fuzzy clustering and TTSAS algorithm for modulation classification based on constellation diagram | |
CN113225282A (zh) | 一种基于bp神经网络的通信信号调制识别方法 | |
CN105516036B (zh) | 一种cpfsk调制识别方法 | |
CN101895494A (zh) | 基于随机共振预处理的数字调制方式自动识别方法 | |
CN111027614B (zh) | 一种噪声增强射频指纹识别方法及装置 | |
CN101764784B (zh) | 多径信道下基于图像处理的正交幅度调制类内识别方法 | |
CN108111213A (zh) | 一种用于多天线的频谱感知方法 | |
Hu et al. | Modulation classification in alpha stable noise | |
CN105763499B (zh) | 一种基于信号瞬时频率滑动平均的cpfsk调制识别方法 |
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 |