CN1395400A - 无需先验知识的常用数字信号调制方式的自动识别方法 - Google Patents

无需先验知识的常用数字信号调制方式的自动识别方法 Download PDF

Info

Publication number
CN1395400A
CN1395400A CN 02123627 CN02123627A CN1395400A CN 1395400 A CN1395400 A CN 1395400A CN 02123627 CN02123627 CN 02123627 CN 02123627 A CN02123627 A CN 02123627A CN 1395400 A CN1395400 A CN 1395400A
Authority
CN
China
Prior art keywords
signal
phase
parameter
thresholding
frequency
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.)
Granted
Application number
CN 02123627
Other languages
English (en)
Other versions
CN1141824C (zh
Inventor
曹志刚
杨志俊
范海波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CNB021236275A priority Critical patent/CN1141824C/zh
Publication of CN1395400A publication Critical patent/CN1395400A/zh
Application granted granted Critical
Publication of CN1141824C publication Critical patent/CN1141824C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明属于通信技术领域,涉及一种无需先验知识的常用数字信号调制方式的自动识别方法。本发明通过A1参数首先分离出噪声,再识别出CW信号,其后对信号进行预处理,分离出模拟信号或未知的非恒包络信号,随后用瞬时频率归一化参数中的Cp将信号分为调频信号和调相信号,调频信号用瞬时频率归一化参数中的峰值数Nf进行分类,调相信号通过平方谱和四次方谱进行分类。通过对BPSK、QPSK、OQPSK、pi/4QPSK、MSK、FSK、FM、CW的仿真表明:其总体识别率达95%以上。本发明具备算法简单、识别速度快捷、识别率高并能达到实用的特点。

Description

无需先验知识的常用数字信号调制方式的自动识别方法
技术领域
本发明属于通信技术领域,特别涉及常用数字信号调制方式的自动识别方法。
背景技术
信号的自动调制识别技术有很高的实用价值。在民用方面,有关当局为了实施有效的无线电频谱管理,常常需要监视民用信号的传输,以便对它们保持控制或者发现和监视未注册登记的发射机;此时,信号认证、干扰识别等都涉及调制的自动识别。在军事和国家安全方面,调制自动识别技术应用更为广泛。任何国家的通信情报系统为了获取通信情报,首先要澄清信号的调制方式,之后才能实施正确解调以及随后的信息处理和分析;在电子战中,为实施电子对抗、电子反对抗、威胁探测、报警、目标捕获和寻的等,都需要通过调制识别技术查清相关通信或电子信号的参数和性质。
现有的通信系统都是根据特定用途和特定应用环境来设计的。不同的通信系统间存在着工作频段、调制方式的差异,导致不同系统间难以实现互通,形成了通信装备品种杂、系列多、互通差、协同难的局面,很难适应未来信息时代方便快捷进行信息交换的需要。针对这样一个现状,1992年5月,Joe.Mitola在美国国家远程系统会议上首次提出了“软件无线电(Software Radio)”的概念。软件无线电系统的目标之一就是通信自动化。通过信号调制方式的自动识别及其它一些相应手段,最终可以实现不同通信系统间的互连互通。因此,信号的自动调制识别是一项富于挑战性的课题。
埃及开罗军事技术学院电子电气工程系E.E.Azzouz和英国格拉斯哥市斯特拉斯克莱德大学(University of Strathclyde)电子电气工程系A.K.Nandi联合编著的《Automatic Modulation Recognition of Communication Signals》一书详细论述了对模拟信号和数字信号的识别方法;另外该书还评述了近年来发表的许多自动调制识别方法。其中的数字调制识别方法工作流程如图1所示,包括以下步骤:
1)假设处理数据是数字调制信号
2)计算γmax,如果其值小于门限 那么判定是移频键控(FSK)信号,则计算σaf,如果大于该门限,则将其它信号归为一类,计算σap
3)计算σaf,如果其值小于门限
Figure A0212362700052
那么判定为二频移频键控(FSK2)信号,反之为四频移频键控(FSK4)信号;
4)计算σap,如果其值大于门限
Figure A0212362700053
那么判定为四相移相键控(PSK4)信号,反之则计算σdp
5)计算σdp,如果其值大于门限
Figure A0212362700061
判定为二相移相键控(PSK2)信号,反之,为ASK信号;
6)对ASK信号计算σaa,如果其值大于门限
Figure A0212362700062
判定为四进制幅度键控(ASK4)信号,反之,为二进制幅度键控(ASK2)信号。
该方法各步骤相关参数说明如下:
1.该方法中零中心归一化瞬时幅度的功率谱密度的最大值γmax为:γmax=max{DFT[acn(i)]}2/Ns式中,Ns为取样点数,acn(i)为零中心归一化瞬时幅度,acn(i)=an(i)-1,an(i)=a(i)/ma m a = 1 N s Σ i = 1 Ns a ( i ) 为瞬时幅度a(i)的均值。
参数γmax主要用来区分是恒包络信号还是非恒包络的信号,因为对调频(FM)、载波(CW)、FSK信号其瞬时幅度为常数(恒定不变),所以它的零中心归一化瞬时幅度acn(i)≈0,对应其谱密度也就为零。而对调幅(AM)、残余边带(DSB)、移相键控(PSK)等信号,由于其瞬时幅度不为恒定值,所以谱密度不为零。假设判决门限为 判决规则如下: γ max ≤ t γ max 时,判为恒包络信号; γ max > t γ max 时,判为非恒包络信号
2.该方法中的非弱信号段零中心归一化瞬时频率的绝对标准差σaf为: σ af = 1 c Σ a n ( i ) > a t f N 2 ( i ) - ( 1 c Σ a n ( i ) > a t | f N ( i ) | ) 2 式中,at是判断弱信号段的一个幅度门限,C是在全部取样数据Ns中属于非弱信号的个数,fN(i)=fm(i)/fc,fm(i)=f(i)-mf m f = 1 N s Σ i = 1 Ns f ( i ) , fc为数字信号的中频,f(i)为信号的瞬时频率。
该参数σaf用来区分是2FSK信号还是4FSK信号。因为2FSK信号的瞬时频率只有两个值,所以它的零中心归一化瞬时频率的绝对值是常数,则其标准差σaf≈0,而4FSK信号的瞬时频率有四个值,所以它的零中心归一化瞬时频率的绝对值不为常数,合理设定判决门限 就可以区分这两种信号。
3.该方法中的非弱信号段零中心瞬时相位非线性分量的绝对标准差σap为: σ ap = 1 c Σ a n ( i ) > a t φ NL 2 ( i ) - ( 1 c Σ a n ( i ) > a t | φ NL ( i ) | ) 2 式中,φNL(i)是经零中心化处理后瞬时相位的非线性分量,在载波完全同步时,φNL(i)=φ(i)-φ0,其中 φ 0 = 1 N s Σ i = 1 Ns φ ( i ) . σap用来区分有相位变化的信号中DSB、BPSK信号还是其它信号。因为DSB和BPSK的φNL只有-π/2和π/2两个取值,不含绝对相位信息。通过选取合适的门限 就可以区分这两类信号。
4.该方法中的非弱信号段零中心瞬时相位非线性分量的直接标准差σdp为: σ dp = 1 c Σ a n ( i ) > a t φ NL 2 ( i ) - ( 1 c Σ a n ( i ) > a t φ NL ( i ) ) 2 式中,σdp用来区分调制信号有无相位变化,可以门限设 区分PSK、FSK、FM、DSB、SSB与AM、ASK、CW调制信号。
5.该方法中的中心归一化瞬时幅度的绝对标准差σaa为: σ aa = 1 N s Σ i = 1 Ns a cn 2 ( i ) - ( 1 N s Σ i = 1 Ns | a cn ( i ) | ) 2 式中,σaa主要用来区分是2ASK信号还是4ASK信号。因为2ASK信号的幅度绝对值是一常数,不含幅度信息,所以有σaa≈0。而4ASK信号的幅度绝对值不是常数,仍含有幅度信息,所以σaa≠0。假设其判决门限为
Figure A0212362700077
该方法存在以下不足:
1.该方法的算法依赖于先验知识,比如说信号准确的载波频率、载波相位、波特率、码元同步及信号采样率与信号码元成整数倍关系等;但实际情况这些参数很难得到,或者只能估计得到准确性较差值,这些参数的误差对信号特征及其判决影响很大。
2.侧重于理论研究,并没有考虑实际调制信号的调制参数的影响,因此其建立的数学模型本身就有较大的误差,可行性不佳。
3.该方法及其它很多文献可以识别信号调制种类有限,要求接收的必须是其讨论范围中的信号,当未知信号是该范围中的调制方式时,可以识别,而当未知信号不是讨论范围中的调制方式时,会误判。
4.并没有对噪声进行有效的分离,这点很重要,因为在实际的通信信号识别中,噪声的分离是第一步的,如果不对其进行分离,无法实现识别的自动化。
5.没有对使用最多的QPSK、OQPSK及pi/4QPSK信号进行分类,而只是简单地将其分为一类。
发明内容
本发明的目的是为克服已有技术的不足,提出一种无需先验知识的常用数字信号调制方式的自动识别方法,使其具有算法简单、识别速度快捷、识别率高并能达到实用的特点。
本发明提出的一种无需先验知识的常用数字信号调制方式的自动识别方法,其特征在于,包括以下步骤:
1)接收待处理数据;
2)分离噪声:计算加性高斯白噪声的A1、A2参数,判断A1、A2的值是否大于门限TA1,如果大于TA1,那么该数据是通信信号,进入下一步处理;如果小于门限TA1,那么判定该数据是噪声,结束流程;
3)计算载波因子的Rcw参数,判断Rcw值是否大于门限TRcw,如果大于该门限,判定该数据为CW信号,结束流程;否则,仍作为未知信号进入下一步处理;
4)估计该未知信号的带宽及载频,对信号进行带通滤波,并将该信号变频到合适的频段;
5)对该未知信号计算R参数;
6)如果R参数的值大于门限TR,将该未知信号归类为模拟信号或未知的非恒包络信号,结束流程;如果小于门限TR,则仍作为未知信号进入下一步处理;
7)对该未知信号计算瞬时频率归一化参数中的Cp,如果小于门限TCp,将信号归类为调频信号,包括调频及移频键控信号;反之,归类为调相信号;
8)对调频信号计算信号瞬时频率归一化参数中的峰值数Nf,如果等于2,则信号为FSK信号,否则为FM信号,结束流程;
9)对调相信号进行平方处理,计算功率谱;
10)计算平方处理后功率谱的载波因子参数Rcw1,如果大于门限T2S,判定为BPSK信号,结束流程;反之,进入下一步处理;
11)计算信号四次方功率谱的N4参数,如果N4=2,则为pi/4 QPSK信号,结束流程;如果N4=0,则归类为未知的恒包络信号,结束流程;如果N4=1,则为QPSK和OQPSK信号,进入下一步处理;
12)计算信号平方功率谱的速率指示参数F2r,如果F2r=1,判定为QPSK信号,否则为OQPSK信号,结束处理流程。
本发明有以下几个方面优点:
1.本发明的方法不依赖于先验知识,包括信号滚降系数或调制指数、信噪比、准确载波频率、载波相位、码元同步等;全面考虑了最常用调制参数对信号特征及调制识别的影响,建立的数学模型与实际非常相近,识别率高;
2.提出一组鲁棒性强的特征参数用于识别,其中首次提出了瞬时频率归一化分布参数和AWGN因子参数,其特征明显,分类效果好;
3.充分考虑了信号识别的完备性,即:当未知信号是识别集合中的元素时,能正确识别,当未知信号不是识别集合中的元素时,也能给出相应的分类指示;
4.有效地对噪声进行识别分离。
5.处理的数据特征提取和算法所用的判决规则简单,关于调制方式的判决可以在很短的时间内完成,完全可以在线分析。
本发明达到的良好效果:
在现代通信系统中最常用的是数字调制技术而非模拟调制技术,因此,对数字调制方式的识别是现代调制识别的主流。本发明研究了基于谱分析的通信信号调制方式自动识别,提出了一组具有鲁棒性强的特征参数,本发明的特点是无需先验知识,能对加性高斯白噪声(AWGN)信道中的常用通信调制方式:BPSK、QPSK、OQPSK、pi/4 QPSK、最小移频键控(MSK)、FSK、FM、CW等信号进行自动识别,并对其它种类的调制信号进行分类。仿真结果表明:其总体识别率达95%以上。
附图说明
图1为已有的一种数字调制自动识别方法流程图。
图2为本发明的数字调制自动识别方法流程图。
具体实施方式
本发明提出的一种无需先验知识的常用数字信号调制方式的自动识别方法,结合附图及实施例详细说明如下:
本发明的调制方式的自动识别方法总体流程如图2所示,包括以下步骤:
1.接收待处理数据;
2.分离噪声:计算加性高斯白噪声(AWGN)的A1、A2参数,判断A1、A2的值是否大于门限TA1,如果大于TA1,那么该数据是通信信号,进入下一步处理;如果小于门限TA1,那么判定该数据是噪声,结束流程;
3.计算载波因子的Rcw参数,判断Rcw值是否大于门限TRcw,如果大于该门限,判定该数据为CW信号,结束流程;否则,仍作为未知信号进入下一步处理;
4.估计该未知信号的带宽及载频,对信号进行带通滤波,并将该信号变频到合适的频段;
5.对该未知信号计算R参数;
6.如果R参数的值大于门限TR,将该未知信号归类为模拟信号或未知的非恒包络信号,结束流程;如果小于门限TR,则仍作为未知信号进入下一步处理;
7.对该未知信号计算信号瞬时频率归一化参数中的Cp,如果小于门限TCp,将信
号归类为调频信号,包括调频(FM)及移频键控(FSK)信号;反之,如果大于门限TCp, 将该信号归类为调相信号,包括BPSK、QPSK、OQPSK、pi/4 QPSK、8相移相键控(8PSK)等;
8.对调频信号计算信号瞬时频率归一化参数中的峰值数Nf,如果等于2,则信号为FSK信号,否则为FM信号,结束流程;
9.对调相信号进行平方处理,计算功率谱;
10.计算平方处理后功率谱的CW因子参数Rcw1,如果大于门限T2S,判定为BPSK信号,结束流程;反之,进入下一步处理;
11.计算信号四次方功率谱的N4参数,如果N4=2,则为pi/4 QPSK信号,结束流程;如果N4=0,则归类为未知的恒包络信号,结束流程;如果N4=1,则为QPSK和OQPSK信号,进入下一步处理;
12.计算信号平方功率谱的速率指示参数F2r,如果F2r=1,判定为QPSK信号,否则为OQPSK信号,结束处理流程。
上述各流程步骤的详细实现方法说明如下:
1、加性高斯白噪声(AWGN)因子参数A1、A2的计算方法:
A1=(E2-E1)/E1、A2=(E3-E2)/E1。其中E1是未知信号功率谱的均值, E 1 = 1 N Σ n = 1 n F ( n ) , F(n)是信号的功率谱;E2是高于E1的功率谱分量值部分的均值;E3是高于E2的功率谱分量值部分的均值,计算经过大量的仿真统计表明:当处理信号的样点数大于2048点时,为AWGN时,其值均在0.7~1.5范围内;当为信号时,A1、A2之值大于2。
2、载波因子Rcw的计算方法对未知信号功率谱F(n)(1≤n≤N),搜索出最大值Fmax及次大值Fsec定义最大值和次大值的比为载波因子Rcw,Rcw=Fmax/Fsec。;显然,当Rcw的值大于门限TRcw时,可判定其为载波信号CW。
3、带宽、载频估计及数字变频实现方法的详细流程如下:
(1)用计算AWGN因子时得到的信号功率谱F(n)的均值E1及小于E1部分的功率谱再次求均值E4;
(2)设定门限Tsn=(E1+E4)*0.5,作为带内信号和噪声的判决门限;
(3)对功率谱F(n)从序号1开始按顺序往上搜索,如果F(n)中连续有四次大于Tsn,则记该位置为pos1;
(4)对功率谱F(n)从序号N/2开始往下搜索F(n),如果F(n)中连续有四次大于Tsn,则记该位置为pos2;
(5)设定起始数字频率F1=pos1-(pos2-pos1)*0.15;
(6)设定截止数字频率F2=pos2+(pos2-pos1)*0.15;
(7)设定信号数字带宽BW=F2-F1,数字载频F0=(F1+F2)*0.5;
(8)在频域对信号进行理想带通滤波,即:对小于F1和大于F2部分的信号FFT得到的复序列S(n)置0;
(9)对信号进行频域数字变频,即:对S(n)进行相应的移位,使其满足,F0=0.15*N;其中N为信号处理长度;
(10)对S(n)进行IFFT变换,得到滤波、变频后的时域信号s(n),后续步骤的所有处理均是对该信号s(n)进行的。
4、R参数的计算方法
R参数反应信号包络的变化程度,R参数的推导及其特性在Y.T.CHAN and L.G.GADBOIS“Identification of The Modulation Type of Signal”Signal Processing 16(1989),pp149~154,North-Holland.文献中有详细的说明,该文献给出R参数的如下特征:AM、单边带调幅(SSB)、DSB等包络变化剧烈的信号,R参数的值也较大;FM、FSK、PSK、CW等信号包络变化较小,R参数的值也较小;通过R参数,可以很好地将这两类信号分类。
5、瞬时频率归一化参数中Cp及Nf的计算方法
未知信号因其准确的载频未知,因此其瞬时相位、瞬时频率存在误差而不可信,而作为统计量的经过数字滤波的信号瞬时频率分量出现概率能较好地反映信号的性能。本发明首次提出该参数,其实现方法如下:
(1)信号s(n)=x(n)+i*y(n)的瞬时相位的公式如下: φ [ n ] = arctan ( y ( n ) x ( n ) ) ( 1 ≤ n ≤ N )
(2)对瞬时相位进行相位去卷叠。相位去卷叠算法是给模2π相位序列φ(n)加上如下校正相位序列Ck(n)
Figure A0212362700112
其中Ck(0)=0,得到去卷叠后的相位序列θ(n)=φ(n)+Ck(n);
(3)去卷叠后的相位序列消除线性相位,得到信号的瞬时频率f1(n): f 1 ( n ) = f 11 ( n ) - 1 N Σ 1 N f 11 ( n ) , 其中,f1(n)为相位序列的差分序列f11(n)=θ(n)-θ(n-1);
(4)对f1(n)进行绝对值大于π/2的统计,得到的统计次数Cp小于门限TCp时,判为没有相位跳变的信号,当Cp大于门限时,判为移相信号。
(5)对没有相位跳变的信号的f1(n)中绝对值大于π/2的值进行插值处理,消除因噪声突变引起的相位突变,插值处理后得到f2(n),插值处理的方法是:
(6)对f2(n)进行归一化处理得到信号的归一化频率f(n): f ( n ) = f 2 ( n ) Max i = 1 N ( abs ( f 2 ( i ) ) )
(7)对f(n)进行分布统计,得到峰值个数Nf。
6、四次方功率谱的N4参数的计算方法信号四次方处理具体步骤如下:
(1)对复信号s(n)=x(n)+i*y(n);实部和虚部分别平方,得到新的序列x2(n)和y2(n),令s2(n)=x2(n)+i*y2(n),对其FFT变换,得到S2(n)=X2(n)+i*Y2(n)
(2)对复数序列S2(n)中序号小于F1的部分清0,以消除信号滚降对以后处理的影响。
(3)对复数序列S2(n),进行频域数字下变频,即:按式S2D(n)=S2(n-F1)对该复数序列移位,得到新的序列S2D(n)(0<n<N+1)
(4)对S2D(n)进行IFFT变换,得到复时域信号s2(n)=x2(n)+i*yd(n)信号。
(5)对sd(n)的实部和虚部分别平方,得到s4(n)=x4(n)+i*y4(n)
(6)对s4(n)进行FFT变换,得到频域信号S4(n)=X4(n)+i*Y4(n)
(7)对S4(n)求功率谱P4(n)=(X4(n))2+(Y4(n))2
(8)对F4(n)序列中序号小于F1的部分清0;
(9)搜索F4(n)序列,得到极大值MAX,次大值SEC,和第三大值THD
(10)如果MAX/THD足够大,大于门限TN4,证明有一条离散谱分量,N4=1;
(11)如果SEC/THD也大于门限TN4,证明有两条离散谱分量,N4=2;
(12)否则,没有足够大的离散谱分量,N4=0;
7、速率指示参数Rate的计算方法速率指示参数Rate的具体算法如下:
(1)对处理数据s(n)=x(n)+i*y(n),求包络 p 2 ( n ) = x 2 ( n ) + y 2 ( n )
(2)对包络数据求对数p1(n)=ln(p2(n))
(3)对p1(n)序列FFT运算,得到复数序列P1(n)
P1(n)=P1r(n)+i*P1i(n)
(4)对P1(n)求模 P ( n ) = P 1 r 2 ( n ) + P 1 i 2 ( n )
(5)从0到BW点,对P(n)搜索,求出极大值max及其的位置pos
(6)对极大值附近区域求平均 comp = 0.1 * ( Σ pos - 9 pos - 5 P ( i ) + Σ pos + 5 pos + 9 P ( i ) )
(7)再与极大值进行比较max/comp,如果其值大于门限TRate,说明有码元速率指示F2r=1,否则F2r=0。
本发明提出的一种无需先验知识的常用数字信号调制方式的自动识别方法实施例的步骤如上所述,不再复述。
本实施例以实际通信系统中最常用的BPSK、QPSK、OQPSK、pi/4 QPSK、FSK、FM、CW等信号及AWGN为例,进行调制方式的自动识别计算机仿真,使用随机序列作为数字调制信号的信源。调制信号码元速率、载波频率、取样频率进行了归一化,分别选定为1、1.5、10。用均值为0正态分布的随机数发生器产生加性高斯白噪声AWGN序列。FSK信号调制指数从0.5到2,步进为0.25;PSK信号的滚降系数从0.3到1.0,步进为0.1;SNR变化范围从6dB到24Db,步进为2dB。每个识别样本用4096个样点。每个识别样本仿真400次。仿真平台操作系统是WIN2000,采用VC进行编程,其中,嵌套了intel公司的信号处理软件包。仿真中各判决门限的值是基于大量仿真统计得到的,其具体设定如下:
加性高斯白噪声的A1、A2参数的判决门限TA1=2;
载波因子的Rcw门限TRcw=100
R参数值的门限TR=0.62
信号瞬时频率归一化参数中的Cp的门限TCp=5
平方处理后功率谱的载波因子参数Rcw1的门限T2S=5
信号平方功率谱的速率指示参数的门限TRate=20
PSK信号中离散谱分量个数Nf检测门限TN4=3

Claims (4)

1、一种无需先验知识的常用数字信号调制方式的自动识别方法,其特征在于,包括以下步骤:
1)接收待处理数据;
2)分离噪声:计算加性高斯白噪声的A1、A2参数,判断A1、A2的值是否大于门限TA1,如果大于TA1,那么该数据是通信信号,进入下一步处理;如果小于门限TA1,那么判定该数据是噪声,结束流程;
3)计算载波因子的Rcw参数,判断Rcw值是否大于门限TRcw,如果大于该门限,判定该数据为CW信号,结束流程;否则,仍作为未知信号进入下一步处理;
4)估计该未知信号的带宽及载频,对信号进行带通滤波,并将该信号变频到合适的频段;
5)对该未知信号计算R参数;
6)如果R参数的值大于门限TR,将该未知信号归类为模拟信号或未知的非恒包络信号,结束流程;如果小于门限TR,则仍作为未知信号进入下一步处理;
7)对该未知信号计算信号瞬时频率归一化参数中的Cp,如果小于门限TCp,将信号归类为调频信号,包括调频及移频键控信号;反之,如果大于门限TCp,将该信号归类为调相信号;
8)对调频信号计算信号瞬时频率归一化参数中的峰值数Nf,如果等于2,则信号为FSK信号,否则为FM信号,结束流程;
9)对调相信号进行平方处理,计算功率谱;
10)计算平方处理后功率谱的载波因子参数Rcw1,如果大于门限T2S,判定为BPSK信号,结束流程;反之,进入下一步处理;
11)计算信号四次方功率谱的N4参数,如果N4=2,则为pi/4QPSK信号,结束流程;如果N4=0,则归类为未知的恒包络信号,结束流程;如果N4=1,则为QPSK和OQPSK信号,进入下一步处理;
12)计算信号平方功率谱的速率指示参数F2r=1,判定为QPSK信号,否则为OQPSK信号,结束整个处理流程。
2、如权利要求1所述的自动识别方法,其特征在于,所说的加性高斯白噪声因子参数A1、A2实现方法为:
A1=(E2-E1)/E1、A2=(E3-E2)/E1,其中E1是未知信号功率谱的均值, E 1 = 1 N Σ n = 1 n F ( n ) , F(n)是信号的功率谱;E2是高于E1的功率谱分量值部分的均值;E3是高于E2的功率谱分量值部分的均值。
3、如权利要求1所述的自动识别方法,其特征在于,所说的瞬时频率归一化参数中Cp及峰值数Nf的具体实现方法为:
(1)信号s(n)=x(n)+i*y(n)的瞬时相位的公式如下: φ [ n ] = arctan ( y ( n ) x ( n ) ) ; ( 1 ≤ n ≤ N )
(2)对瞬时相位进行相位去卷叠,相位去卷叠算法是给模2π相位序列φ(n)加上如下校正相位序列Ck(n):其中Ck(0)=0,得到去卷叠后的相位序列θ(n)=φ(n)+Ck(n);
(3)去卷叠后的相位序列消除线性相位,得到信号的瞬时频率f1(n): f 1 ( n ) = f 11 ( n ) - 1 N Σ 1 N f 11 ( n ) , 其中,f1(n)为相位序列的差分序列f11(n)=θ(n)-θ(n-1);
(4)对f1(n)进行绝对值大于π/2的统计,得到的统计次数Cp小于门限TCp时,判为没有相位跳变的信号,当Cp大于门限TCp时,判为移相信号;
(5)对没有相位跳变的信号的f1(n)中绝对值大于π/2的值进行插值处理,消除因噪声突变引起的相位突变,插值处理后得到f2(n),插值处理的方法是:
Figure A0212362700036
(6)对f2(n)进行归一化处理得到信号的归一化频率f(n): f ( n ) = f 2 ( n ) Max i = 1 N ( abs ( f 2 ( i ) ) )
(7)对f(n)进行分布统计,得到其峰值数Nf。
4、如权利要求1所述的自动识别方法,其特征在于,所说的速率指示参数Rate的具体算法如下:
(1)对处理数据s(n)=x(n)+i*y(n),求包络数据 p 2 ( n ) = x 2 ( n ) + y 2 ( n )
(2)对包络数据求对数p1(n)=ln(p2(n))
(3)对p1(n)序列进行FFT运算,得到复数序列P1(n)
P1(n)=P1r(n)+i*P1i(n)
(4)对P1(n)求模 P ( n ) = P 1 r 2 ( n ) + P 1 i 2 ( n )
(5)从0到BW点,对P(n)搜索,求出极大值max及其的位置pos
(6)对极大值附近区域求平均 comp = 0.1 * ( Σ pos - 9 pos - 5 P ( i ) + Σ pos + 5 pos + 9 P ( i ) )
(7)再与极大值进行比较max/comp,如果其值大于门限TRate,则有码元速率指示
F2r=1,否则F2r=0。
CNB021236275A 2002-07-05 2002-07-05 无需先验知识的常用数字信号调制方式的自动识别方法 Expired - Fee Related CN1141824C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB021236275A CN1141824C (zh) 2002-07-05 2002-07-05 无需先验知识的常用数字信号调制方式的自动识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB021236275A CN1141824C (zh) 2002-07-05 2002-07-05 无需先验知识的常用数字信号调制方式的自动识别方法

Publications (2)

Publication Number Publication Date
CN1395400A true CN1395400A (zh) 2003-02-05
CN1141824C CN1141824C (zh) 2004-03-10

Family

ID=4745169

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB021236275A Expired - Fee Related CN1141824C (zh) 2002-07-05 2002-07-05 无需先验知识的常用数字信号调制方式的自动识别方法

Country Status (1)

Country Link
CN (1) CN1141824C (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1707610B (zh) * 2004-06-04 2012-02-15 本田研究所欧洲有限公司 对两个谐波信号的共同起源的确定
CN102664724A (zh) * 2011-09-20 2012-09-12 北京华力创通科技股份有限公司 数字通信系统中的数字信号的处理方法和装置
CN102707157A (zh) * 2012-05-18 2012-10-03 天津理工大学 一种基于功率谱的单频脉冲信号参数估值方法
CN102749615A (zh) * 2012-07-11 2012-10-24 天津理工大学 一种信号识别的方法
CN103199945A (zh) * 2013-03-29 2013-07-10 南京信息职业技术学院 一种低信噪比条件下认知无线电信号调制方式的识别方法
CN104022837A (zh) * 2014-05-28 2014-09-03 广州海格通信集团股份有限公司 一种基于特征提取的通信信号调制方式进行识别的方法
CN104270331A (zh) * 2014-03-14 2015-01-07 南京信息职业技术学院 基于hough变换的bpsk信号盲识别结果的有效性评估方法
CN104507106A (zh) * 2014-12-31 2015-04-08 中国电子科技集团公司第五十研究所 8psk信号和pi/4-dqpsk信号的识别方法
CN104639268A (zh) * 2014-12-31 2015-05-20 中国电子科技集团公司第五十研究所 Qpsk信号和oqpsk信号的识别方法
CN105790848A (zh) * 2014-12-23 2016-07-20 深圳市中兴微电子技术有限公司 一种载波相位估计方法和装置
CN106330805A (zh) * 2016-08-29 2017-01-11 重庆会凌电子新技术有限公司 一种信号调制方式自动识别方法及系统
CN106357575A (zh) * 2016-10-17 2017-01-25 中国电子科技集团公司第五十四研究所 一种多参数联合估计的干扰类型识别方法
CN108205079A (zh) * 2016-12-16 2018-06-26 北京普源精电科技有限公司 利用频谱仪检测2fsk信号的方法及装置
CN108282428A (zh) * 2018-01-23 2018-07-13 大连理工大学 一种无需先验知识的常用通信信号调制方式的自动识别方法
CN108965181A (zh) * 2018-09-25 2018-12-07 电子科技大学 一种用于数据链信号侦收系统的信号调制方式识别方法
CN110247868A (zh) * 2019-06-26 2019-09-17 北京邮电大学 一种卫星信号盲解调方法、装置及电子设备
CN111444805A (zh) * 2020-03-19 2020-07-24 哈尔滨工程大学 一种基于改进的多尺度小波熵数字信号调制识别方法
CN117201249A (zh) * 2023-11-01 2023-12-08 中孚安全技术有限公司 一种信号调制方式识别方法、系统及装置

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1707610B (zh) * 2004-06-04 2012-02-15 本田研究所欧洲有限公司 对两个谐波信号的共同起源的确定
CN102664724A (zh) * 2011-09-20 2012-09-12 北京华力创通科技股份有限公司 数字通信系统中的数字信号的处理方法和装置
CN102664724B (zh) * 2011-09-20 2015-03-25 北京华力创通科技股份有限公司 数字通信系统中的数字信号的处理方法和装置
CN102707157A (zh) * 2012-05-18 2012-10-03 天津理工大学 一种基于功率谱的单频脉冲信号参数估值方法
CN102707157B (zh) * 2012-05-18 2015-10-28 天津理工大学 一种基于功率谱的单频脉冲信号参数估值方法
CN102749615B (zh) * 2012-07-11 2015-08-05 天津理工大学 一种信号识别的方法
CN102749615A (zh) * 2012-07-11 2012-10-24 天津理工大学 一种信号识别的方法
CN103199945A (zh) * 2013-03-29 2013-07-10 南京信息职业技术学院 一种低信噪比条件下认知无线电信号调制方式的识别方法
CN103199945B (zh) * 2013-03-29 2016-06-01 南京信息职业技术学院 一种低信噪比条件下认知无线电信号调制方式的识别方法
CN104270331A (zh) * 2014-03-14 2015-01-07 南京信息职业技术学院 基于hough变换的bpsk信号盲识别结果的有效性评估方法
CN104022837B (zh) * 2014-05-28 2016-08-17 广州海格通信集团股份有限公司 一种基于特征提取的信号调制方式识别方法
CN104022837A (zh) * 2014-05-28 2014-09-03 广州海格通信集团股份有限公司 一种基于特征提取的通信信号调制方式进行识别的方法
CN105790848B (zh) * 2014-12-23 2018-06-08 深圳市中兴微电子技术有限公司 一种载波相位估计方法和装置
CN105790848A (zh) * 2014-12-23 2016-07-20 深圳市中兴微电子技术有限公司 一种载波相位估计方法和装置
CN104507106A (zh) * 2014-12-31 2015-04-08 中国电子科技集团公司第五十研究所 8psk信号和pi/4-dqpsk信号的识别方法
CN104639268A (zh) * 2014-12-31 2015-05-20 中国电子科技集团公司第五十研究所 Qpsk信号和oqpsk信号的识别方法
CN104507106B (zh) * 2014-12-31 2018-07-27 中国电子科技集团公司第五十研究所 8psk信号和pi/4-dqpsk信号的识别方法
CN104639268B (zh) * 2014-12-31 2017-04-12 中国电子科技集团公司第五十研究所 Qpsk信号和oqpsk信号的识别方法
CN106330805B (zh) * 2016-08-29 2019-08-09 重庆会凌电子新技术有限公司 一种信号调制方式自动识别方法及系统
CN106330805A (zh) * 2016-08-29 2017-01-11 重庆会凌电子新技术有限公司 一种信号调制方式自动识别方法及系统
CN106357575A (zh) * 2016-10-17 2017-01-25 中国电子科技集团公司第五十四研究所 一种多参数联合估计的干扰类型识别方法
CN108205079A (zh) * 2016-12-16 2018-06-26 北京普源精电科技有限公司 利用频谱仪检测2fsk信号的方法及装置
CN108282428A (zh) * 2018-01-23 2018-07-13 大连理工大学 一种无需先验知识的常用通信信号调制方式的自动识别方法
CN108965181A (zh) * 2018-09-25 2018-12-07 电子科技大学 一种用于数据链信号侦收系统的信号调制方式识别方法
CN110247868A (zh) * 2019-06-26 2019-09-17 北京邮电大学 一种卫星信号盲解调方法、装置及电子设备
CN111444805A (zh) * 2020-03-19 2020-07-24 哈尔滨工程大学 一种基于改进的多尺度小波熵数字信号调制识别方法
CN111444805B (zh) * 2020-03-19 2023-03-17 哈尔滨工程大学 一种基于改进的多尺度小波熵数字信号调制识别方法
CN117201249A (zh) * 2023-11-01 2023-12-08 中孚安全技术有限公司 一种信号调制方式识别方法、系统及装置
CN117201249B (zh) * 2023-11-01 2024-02-20 中孚安全技术有限公司 一种信号调制方式识别方法、系统及装置

Also Published As

Publication number Publication date
CN1141824C (zh) 2004-03-10

Similar Documents

Publication Publication Date Title
CN1141824C (zh) 无需先验知识的常用数字信号调制方式的自动识别方法
CN101834819B (zh) 基于平行判决的模拟数字混合调制方式识别装置和数字调制方式识别装置
CN107948107B (zh) 一种基于联合特征的数字调制信号分类方法
CN108718288B (zh) 基于卷积神经网络的数字信号调制模式识别方法
CN104301056B (zh) 一种基于信号特征分析的频谱监测方法
CN102710572B (zh) 一种通信信号的特征提取与调制识别方法
CN111800359B (zh) 一种通信信号调制方式的识别方法、装置、设备及介质
CN104468001B (zh) 基于无线电信号频谱特征模板的信号识别方法及系统
CN111814777B (zh) 一种基于特征量分级的调制样式识别方法
CN1750523A (zh) 连续相位调制自适应跳频系统的信道质量估计方法
CN106330805B (zh) 一种信号调制方式自动识别方法及系统
CN108282428A (zh) 一种无需先验知识的常用通信信号调制方式的自动识别方法
Ahmadi Using fuzzy clustering and TTSAS algorithm for modulation classification based on constellation diagram
CN104363194A (zh) 基于波形变换的psk调制识别方法
CN108206797A (zh) 一种基于软件无线电的通信信号调制方式自适应识别方法
JP2001127810A (ja) 変調方式認識装置およびその方法
CN108270703A (zh) 一种通信信号数字调制类型识别方法
CN111371510A (zh) 一种复杂电磁环境下的通信测试方法
CN106453174B (zh) 基于信号调制识别的认知无线网络资源分配方法
CN111245756B (zh) 基于级联svm和全数字接收机的复合信号调制识别方法
CN103546418A (zh) 基于峰均比的调制方式的识别方法及装置
Ye et al. Digital modulation classification using temporal waveform features for cognitive radios
Li et al. A novel modulation classification for PSK and QAM signals in wireless communication
Kumar SoC implementation of a modulation classification module for cognitive radios
Chuangzhan et al. Modulation recognition method of communication signals based on correlation characteristics

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee