本发明属于通信侦察接收机领域，特别涉及一种基于侦察接收机的调制模式识别方法。
背景技术：
随着通信技术的高速发展,实际通信环境中信号日趋密集。不同功率、带宽和调制模式的通信信号同时存在,使得实际通信环境愈加复杂。在信号侦察中，现代侦察接收机必须具备在复杂情况下动态变化的信号环境中实时对信号进行检测和分类识别的能力。而在实际远距离侦察中，接收端接收到的信号已经十分微弱，再加上环境和仪器噪声的影响，使得信号的分析识别更加困难。
目前主要的识别算法包括频谱检测法、星座图法和基于高阶累计量的识别方法等。频谱检测法由于实现简单、复杂度低，是调制模式识别算法中一种广泛使用的识别方法，但当输入信号存在干扰时，频谱检测识别率很难满足实际需求，而且实际识别中bpsk和qpsk信号频谱相似，无法通过频谱检测区分。星座图可以反映不同调制模式信号的结构关系，通过信号星座图分布特征进行调制模式识别是最稳健的方法之一。但当输入信号存在输入信号幅度、载波频率动态变化星座图对信号进行识别会出现大量虚警和误识。基于高阶累计量的方法由于其抗多径干扰的能力得到广泛的关注和研究，但是实际识别过程中ask和bpsk以及mfsk之间高阶累计量相同,无法区分，且高阶累计量算法复杂度高，不易于硬件实现。
总的来说，目前识别算法主要针对解调后的基带信号进行研究，对于实际环境中接收到的输入信号幅度、载波频率高度动态变化的信号，识别效果并不好。在处理能量较低的信号时，由于仪器本身和环境噪声对信号影响很大，以上这些方法均不能对目标信号做出准确判断，无法满足实际工程要求。并且目前的识别算法实现复杂度高、资源占用过大、不适于硬件实现。
技术实现要素：
本发明的目的是针对现有技术的不足，对于实际环境中存在的bpsk、qpsk、fsk信号在载波频率和信号能量高动态变化时难以识别的问题，提出了一种基于侦察接收机的调制模式识别方法。
本发明解决其技术问题所采用的技术方案包括以下步骤：
步骤1、将输入信号通过自动增益控制处理(agc),将不同能量的信号幅值统一；
步骤2、通过fft对agc处理后的信号做频率粗估计，估计频率与输入信号下变频；
步骤3、对下变频之后信号做位同步处理，消除收发端时钟偏移；
步骤4、通过载波同步跟踪剩余频偏，解调出基带信号；
步骤5、利用频率跟踪曲线区分出psk和fsk信号，利用基带信号区分bpsk和qpsk信号。
步骤1具体实现如下：
1-1.对输入信号幅值估计，之后与参考电平r比较；
1-2.将比较后的值与收敛因子α相乘后作对数和指数处理，通过对数函数和指数函数的非线性特性来控制反馈函数；
1-3.将输入信号x(n)与反馈函数相乘得到增益控制后的输出信号，递推公式如下:
y(n)＝x(n)exp(a(n))(1)
a(n+1)＝a(n)+α[r-y(n)](2)
其中a(n)为控制信号，r为参考电平，α为收敛因子。此时信号已经被统一到同一幅值。
步骤2具体实现如下：
2-1.将输入信号实时分组，每组1024个数据点，存储到寄存器中；
2-2.对每组数据做1024点fft，取频谱的最大值和第二大值对应横坐标的平均值作为估计频率；
2-3.当第五组数据输入时，释放之前保存的第一组数据，将第五组数据保存到第一个寄存器中，如此四组数据为一次循环，循环32次做128组fft；
2-4.将每次fft的估计频率加权平均，得到最后的载波频率估计值；
2-5.将agc处理后数据与fft估计频率做下变频。此过程中可以将信号写成两路形式，表达式如下
s(t)＝i(t)cos(2πfct)+q(t)sin(2πfct)+n(t)(3)
其中i(t)、q(t)分别为输入s(t)的i、q两路基带数据，fc表示其载波频率，n(t)为服从n(0,σ2)的高斯白噪声。
下变频后，暂时不考虑噪声，信号表达式如下：
i(k)＝ikcos(2πkδf/fs)-qksin(2πkδf/fs)(4)
q(k)＝iksin(2πkδf/fs)+qkcos(2πkδf/fs)(5)
其中，i(k)为同相支路，q(k)为正交支路；ik、qk为支路系数，k为第k个采样点。
步骤3具体实现如下：
3-1.对输入信号做线性内插，公式如下
y(k)＝x(mk)+μk[x(mk+1)-x(mk)](6)
其中μk为环路估计误差；
3-2.通过同步误差检测器检测接收信号的定时误差值，公式如下
3-3.使用预滤波进一步减小进入环路滤波的误差抖动，传输函数如下
其中a为滤波因子，且a＜1。
3-4.对信号相位使用二阶环路滤波滤除高频分量，同时对相位误差进行调节以使系统达到稳定。
3-5.以位定时nco溢出时刻作为采样允许时钟产生的标志，此时将环路估计误差反馈到内插滤波器。
3-6.重复步骤3-1～3-5，直到位同步环路稳定收敛。
步骤4具体实现如下：
4-1.位同步后信号与载波同步中本地nco反馈频率下变频后低通滤波。此时信号i路q路表达式如下
其中
4-2.通过鉴相方程计算出信号相位，计算过程如下
其中sign(*)为符号函数。
将式(9)(10)代入式(11)有
4-3.对信号相位使用二阶环路滤波滤除高频分量，同时对相位误差进行调节以使系统达到稳定。
4-4.更新nco中本地频率。
4-5.重复步骤4-1～4-4，直到载波同步环路收敛稳定，当环路稳定时
步骤5具体实现如下：
5-1.将接收到的频偏跟踪数据实时分组，判断该组数据的稳定性，当检测到跟踪载频出现较大跳变时，则将该组信号判定为fsk信号。当跟踪频率稳定在一个固定值时，进入后续判断；
5-2.将接收到的解调数据实时分组，每组128个点，统计i，q两路中符号相同的数据个数，将识别阈值设为100，如果超过设定的阈值即判定为bpsk信号，反之低于阈值则判定为qpsk信号。
本发明有益效果如下：
1.本发明针对先验信息完全未知的信号，在全盲的情况下对信号进行调制模式识别。
2.本发明能够消除信号能量和载波频率对调制模式识别的影响，提高了低能量低信噪比的信号的识别率。
3.在信号输入幅度和频率动态变化的的情况下，本发明能够对bpsk、qpsk以及fsk信号在接收机同步系统中进行准确地判断，识别率高于传统识别方法。
4.本发明只需要对两种特征参量进行识别判定，减少了计算量和系统资源占用，易于硬件实现。
附图说明
图1本发明中侦察接收机系统流程图
图2本发明中增益自动控制模块算法框图；
图3本发明中多窗谱加权fft流程图；
图4本发明中位同步结构图；
图5本发明中载波同步结构图；
图6二阶环路滤波器结构图；
图7fsk信号载波同步频偏跟踪曲线；
图8psk信号载波同步频偏跟踪曲线；
图9bpsk信号解调后基带信号；
图10qpsk信号解调后基带信号；
具体实施方式
下面结合附图对本发明的具体实施例做进一步的说明。
如图1所示，一种基于侦察接收机的调制模式识别方法，具体实现过程如下：
步骤1、将输入信号通过自动增益控制处理(agc),将不同能量的信号幅值统一，如图2所示；
步骤2、通过fft对agc处理后的信号做频率粗估计，估计频率与输入信号下变频；
步骤3、对下变频以后信号做位同步处理，消除收发端时钟偏移，具体过程如图4所示；
步骤4、通过载波同步跟踪剩余频偏，解调出基带信号，如图5所示；
步骤5、利用频偏跟踪曲线区分psk和fsk信号，利用基带信号区分bpsk和qpsk信号。
步骤1具体实现如下：
1-1.输入信号幅值估计后与参考电平比较；
1-2.将比较后值与收敛因子相乘后作对数和指数处理，通过对数函数和指数函数的非线性特性来控制反馈函数；
1-3.将输入信号与反馈函数相乘得到增益控制后的输出信号，递推公式如下
y(n)＝x(n)exp(a(n))(1)
a(n+1)＝a(n)+α[r-y(n)](2)
其中a(n)为控制信号，r为参考电平，α为收敛因子。此时信号已经被统一到同一幅值。
步骤2具体实现如下：
2-1.将输入信号实时分组，每组1024个数据点，存储到寄存器中；
2-2.对每组数据做1024点fft，取频谱的最大值和第二大值对应横坐标的平均值作为估计频率；
2-3.当第五组数据输入时，释放之前保存的第一组数据，将第五组数据保存到第一个寄存器中，如此四组数据为一次循环，循环32次做128组fft；
2-4.将每次fft的估计频率加权平均，得到最后的载波频率估计值；
2-5.将agc后信号与fft估计频率做下变频。将信号写成两路形式，表达式如下
s(t)＝i(t)cos(2πfct)+q(t)sin(2πfct)+n(t)(3)
其中i(t)、q(t)分别为信号s1(t)的i、q两路基带数据，fc表示其载波频率，n(t)为服从n(0,σ2)的高斯白噪声。
下变频后，暂时不考虑噪声，信号表达式如下
i(k)＝ikcos(2πkδf/fs)-qksin(2πkδf/fs)(4)
q(k)＝iksin(2πkδf/fs)+qkcos(2πkδf/fs)(5)
其中，i(k)为同相支路，q(k)为正交支路；ik、qk为支路系数，k为第k个采样点。
步骤3具体实现如下：
3-1.对输入信号做线性内插，公式如下
y(k)＝x(mk)+μk[x(mk+1)-x(mk)](6)
其中μk为环路估计误差；
3-2.通过同步误差检测器检测接收信号的定时误差值，公式如下
3-3.使用预滤波进一步减小进入环路滤波的误差抖动，传输函数如下
其中a为滤波因子，且a＜1。
3-4.对信号相位使用二阶环路滤波滤除高频分量，同时对相位误差进行调节使系统达到稳定，如图6所示。其中c1，c2为环路系数，由以下公式计算
其中
3-5.以位定时nco溢出时刻作为采样允许时钟产生的标志，此时将环路估计误差反馈到内插滤波器。
3-6.重复步骤3-1～3-5，直到位同步环路稳定收敛。
步骤4具体实现如下：
4-1.位同步后信号与载波同步nco反馈频率下变频后低通滤波。此时信号i路q路表达式如下
其中
4-2.通过鉴相方程计算出信号相位，计算过程如下
其中sign(*)为符号函数，i′(k)、q′(k)分别为信号下变频后的同相和正交支路。
将式(10)(11)代入式(12)有
4-3.对信号相位使用二阶环路滤波滤除高频分量，同时对相位误差进行调节使系统达到稳定，如图6所示。其中c1，c2为环路系数，由以下公式计算
其中
4-4.更新nco中本地频率，返回步骤4-1直到载波同步环路收敛稳定，当环路稳定时
步骤5具体实现如下：
5-1.将接收到的频偏跟踪数据实时分组，判断该组数据的稳定性，当检测到跟踪载频出现较大跳变时，如图7所示，则将该组信号判定为fsk信号。当跟踪频率稳定在一个固定值时，如图8所示，则判定为psk信号，进入后续判断；
5-2.将接收到的psk信号解调数据实时分组，每组128个点，统计i，q两路中符号相同的数据个数，将识别阈值设为100，如果超过设定的阈值即判定为bpsk信号，如图9所示；反之低于阈值则判定为qpsk信号，如图10所示。
实施例：
采用中国科技电子集团提供的空口信号输入侦察接收机系统，经过安捷伦矢量信号发生器m9381a射频调制，信号频率由1720m下变频到25m～75mhz，信号能量为-10dbm～-60dbm，输入300组不同调制方式的信号，每组128个数据点。
先将输入信号通过agc处理，统一输入信号的幅度。agc处理后，对信号做fft频率粗估计，以纠正大频偏。接下来将信号经过位同步，消除收发端的延迟和时钟偏移。最后将信号通过载波跟踪环，跟踪剩余频差并解调出基带信号。在调制模式识别中，先通过频率跟踪曲线对fsk和psk进行区分，之后通过解调出的基带信号对bpsk和qpsk信号进行区分。通过统计实验结果，本发明对于300组-45dbm、-55dbm和-60dbm的输入信号，识别率分别为99.66％、98.69％和95.10％，而采用高阶累积量进行识别，识别率分别为90.38％、88.65％和83.33％。
实施例结果表明，本专利的基于侦察接收机的调制模式识别算法在在降低了系统复杂度和运算量的同时，保证了识别率，且更加适合硬件实现。
最后需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附的权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。