无线电引信信号识别和重构系统及方法与流程

本发明属于雷达通信领域，具体涉及一种无线电引信信号识别和重构系统及方法。

背景技术：

现代战争中，为对敌方雷达引信实施干扰，需要准确获取引信信号的相应参数，这样才能做到有的放矢，在电子对抗中占据主动。然而在复杂战场环境中对雷达引信信号进行探测时，极有可能接收到各种不同调制制式的雷达引信信号，这就需要采取相应的方法识别出调制制式，方便利用各种对应方法提取雷达引信信号的参数。

在无线电引信信号领域，目前的识别方法主要存在以下问题：

(1)大多数的识别方法都是针对某一种特定的信号类型，在先验条件较少条件下的盲识别方法比较匮乏。

(2)大多数的识别方法是单一地基于信号时频分布而进行识别，这种方法在信噪比较为恶劣的条件下，难以保证识别准确率，尤其是对于脉冲串及伪随机脉位调制信号等的识别，情况不太理想。

(3)当前的识别方法均是在软件平台上的仿真，而没有在硬件的实验平台上的具体实践工作，这就导致这些结果具有一定的局限性，同实际的应用还有一定距离。

技术实现要素：

本发明的目的在于提供一种无线电引信信号识别和重构系统及方法，解决当前方法识别制式类型少、只能对某种特定信号识别、缺乏硬件平台以及参数提取准确率较低等问题。

实现本发明目的的技术方案为：一种无线电引信信号识别和重构系统，包括搭建在labview平台的信号产生与发射模块、信号接收模块、信号制式识别模块、信号参数提取和重构模块；

信号产生与发射模块用于各种制式信号的生成和射频信号的发射；

信号接收模块用于射频信号的接收，并将其转化为数字信号供信号制式识别模块识别；

信号制式识别模块用于对接收到的信号进行识别，判断其调制制式；

信号参数提取和重构模块用于根据判断出的信号制式进行参数提取，再根据所得参数进行重构。

一种无线电引信信号识别和重构方法，包括以下步骤：

步骤1，设置需发射信号的制式以及各项参数，系统生成该信号的时域波形图，将中频调制后的射频信号发射出去；

步骤2，接收射频信号，并将其转化为数字信号；

步骤3，对数字信号进行识别，判断其调制制式；

步骤4，根据判断出的信号制式，选择与之对应的方法进行参数提取，再根据所得参数构建出一个新的信号。

与现有技术相比，本发明的显著优点为：

本发明搭建了一个完善的制式识别、参数提取及重构系统，具有环境模拟度强，识别速度快，参数提取及恢复准确率高等特性。

附图说明

图1是pxie5673射频信号发生器示意图。

图2是pxie5663射频信号矢量分析仪示意图。

图3是不同体制下无线电引信信号的识别流程图。

图4是线性调频最优阶数分数阶傅里叶变换示意图。

具体实施方式

一种无线电引信信号识别和重构系统，包括搭建在labview平台的信号产生与发射模块、信号接收模块、信号制式识别模块、信号参数提取和重构模块；

信号产生与发射模块用于各种制式信号的生成和射频信号的发射；

信号接收模块用于射频信号的接收，并将其转化为数字信号供信号制式识别模块识别；

信号制式识别模块用于对接收到的信号进行识别，判断其调制制式；

信号参数提取和重构模块用于根据判断出的信号制式进行参数提取，再根据所得参数进行重构。

一种基于上述系统的无线电引信信号识别和重构方法，包括以下步骤：

步骤1，设置需发射信号的制式以及各项参数，系统生成该信号的时域波形图，将中频调制后的射频信号发射出去；

步骤2，接收射频信号，并将其转化为数字信号；

步骤3，对数字信号进行识别，判断其调制制式；

步骤4，根据判断出的信号制式，选择与之对应的方法进行参数提取，再根据所得参数构建出一个新的信号。

进一步的，如图1所示，步骤1具体为：

信号产生与发射模块将labview软件界面上所产生的数字信号接入i-q信号发生器；

将输出的中频信号与本振源提供的载频信号一同进入上变频器进行变频；

将变频器所产生的射频信号经天线发出，即为该模块产生的射频信号。

进一步的，如图2所示，步骤2具体为：

信号接收模块将步骤1产生的射频信号经天线接收下来，与本振源提供的载频信号一同进入下变频器，下变频后得到一模拟中频信号；

将该模拟中频信号输入中频数字化器得到数字信号。

进一步的，如图3所示，步骤3通过占空比的求取，判断其类型为脉冲串或连续波；对于连续波信号，由信号在频谱上的3db带宽b3db，区分单频信号和频率受调制信号；针对单频信号，根据瞬时相位特性分为多普勒连续波和伪码调相信号；针对调频信号，根据线性调频因子值分为线性调频和正弦调频信号；针对线性调频信号，分为锯齿波线性调频信号和三角波线性调频信号；具体为：

(1)脉冲串信号和连续波信号识别：

假设侦收到的引信信号为s[n]，通过下式计算其包络波形：

as[n]＝|s[n]|

当信号无噪时，很容易看出信号是否为脉冲串；然而当信号参杂噪声时，其包络会在一定范围内波动，通过设定门限阈值pd来判断该信号是否为脉冲串信号：

当信噪比较低时，as1[n]中会出现不规则的跳变点，从而造成识别错误；为避免这种情况发生，在波形发生跳变后连续观察多组跳变后的波形值，若连续多点始终保持跳变后的值，则认为该跳变位置为脉冲沿，否则认为是噪声引起的跳变；具体过程可用下式表示：

对于as2[n]，计算信号占空比τ，当占空比大于阈值τd时，判为连续波信号，否则判为脉冲串信号；

(2)单频信号和频率受调制信号识别：

对侦收到的引信信号进行傅立叶变换，并求出信号3db带宽b3db，设定3db带宽阈值bd，若b3db＞bd，侦收到的信号为调频类信号；若b3db≤bd，则侦收到的信号为单频信号；

(3)连续波多普勒信号和伪码调相信号识别：

a)对侦收到的引信信号进行hilbert变换，求取信号的瞬时相位

b)对瞬时相位进行一次差分得

c)为减小噪声对的影响，做一个滑动窗口，每次滑动一个数据单位，在这个窗口下依次对求平均，命名平均值为j＝1,2…p，p为滑动次数；

d)检测的范围，即选取相位阈值若则侦收到的信号为伪码调相信号；若则侦收的信号为连续波多普勒信号；

(4)线性调频信号和正弦调频信号识别：

针对线性调频和正弦调频，由分数阶傅里叶变换对线性调频信号的聚敛特性进行识别；判断方法如下：

首先，截取信号周期长度的二分之一；

然后，在[0,2]内以δp为步进，取n个点；分别对信号进行n次傅里叶变换，得到对应个数的分数阶域谱frftp(f)并找出每个谱上的最大值mp＝max(frftp(f))，在所有阶数所对应的mp中找出最大值max(mp)，定义σf＝max(mp)，得到与该最大值对应的阶数pmax＝arg(max(mp))；

设置线性调频因子阈值σd来判断该信号是否为线性调频信号，若σf＞σd，侦收到的信号为线性调频信号；若σf≤σd，则侦收到的信号为正弦调频信号；

(5)锯齿波线性调频和三角波线性调频信号识别：

首先，通过对信号做短时傅立叶变换得到一个时频能量分布的三维图形；

然后，对其提取时频脊线得到一个二维上的瞬时频率图；由图中最大值所处横坐标t和采样长度l(td＝0.7*l)之间关系，判断该信号是否为锯齿波调频信号；具体过程为：设一个横坐标阈值td，若t＞td，侦收到的信号为锯齿波调频信号；若t≤td，则侦收到的信号为三角波调频信号。

进一步的，步骤4中参数提取的具体过程为：

(1)锯齿波线性调频信号：

锯齿连续波线性调频信号表示为：

其中，f0为载频频率，tm为调制周期，k为调频斜率；

a.调制周期tm

采用自相关函数提取其调制周期；

当-tm≤τ＜0时，τ为位移量，自相关函数表示为：

当0≤τ＜tm时，自相关函数表示为：

于是，单周期基带线性调频信号自相关函数的幅值为：

由上式可知，当τ＝0和τ＝tm时，可得到|ru(τ)|的最大值，因此，通过检测连续波线性调频信号自相关峰值出现的位置，即可得到信号的调制周期；

b.调制斜率k

结合图4，当分数阶傅里叶变换选取的阶数与线性调频信号的调频斜率相对应时，线性调频信号可在对应的分数阶频域达到最大的能量积累；

通过不断改变frft的阶数p＝-2α/π，α为旋转角度，计算线性调频信号对应的分数阶域谱frftp(f)，对比求取所有分数阶谱的最大幅值mp＝max(frftp(f))并记录该最大幅值对应的阶数pmax，最后最优阶数pmax和调频斜率k之间的对应关系估计信号的调频斜率,fs为采样率，n为信号点数；

c.载波频率f0

求得斜率k之后，在原始信号上乘以exp(-jπkt²)因子，即得到单载频信号；将去斜率的信号进行傅立叶变换，求出信号的载频：

f0＝arg{max{fft[s(t)·exp(-jπkt²)]}}

(2)正弦调频连续波信号：

正弦调频信号表示为：

ssfm＝a·exp{j2π[f0(ntm)+m·cos(2π·fsfm·ntm+ψsfm)]}

其中，tm为调制周期，fo为载频，m为信号的调频指数，fsfm为调制频率，ψsfm为调频初始相位；

a.调制周期tm提取过程同锯齿波线性调频信号；

b.调频初始相位及调制指数

设定变量ψ∈[0,2π]和fs为采样率，并建立一个参考信号：

式中，a为幅值，ts为采样周期；

与未知信号相乘后，做傅里叶变换得函数p：

此时可以对这两个参数均设定相应的取值范围，利用二次搜索估计参数，然而同时对这两个参数进行精确搜索会造成极大的运算量，因此可以采用先进行一次二维搜索，再进行一次一维搜索的方法来提高效率：

首先使m的搜索步进较大，ψ的搜索步进较小，进行二维搜索，估计出再以为确定值，使m的搜索步进较小，进行一维搜索；

c.载波频率f0

通过傅立叶变换得到载频：

(3)单频连续波信号

单载频连续波信号表示为：

sdo＝exp{j[2πf0(nts)]+ψ0}

其中，n为信号点数，ts为采样周期，fo为载频，ψ0为初始相位；

a.载频估计

(4)单频伪码调相信号

单载伪码调相信号表示为：

其中，ts为采样周期，fo为载频，ψo为初始相位，q是伪随机序列的长度，tc为伪码码元宽度，cq是伪随机序列，取值为1或-1；

a.采用自相关函数提取采样周期ts；

b.载波频率f0

c.伪随机序列

可知未经调相的载波信号取其共轭与原信号相乘得到去载频信号：

对进行阈值判断，确定伪随机序列；

(5)脉冲串信号：

脉冲串信号表示为：

其中，n为信号点数，tr为脉冲重复周期；

a.脉冲重复周期tr提取过程同上；

b.脉冲序列

先将信号做求模处理，再按照从大到小的顺序提取信号的前几个最大值，求其平均值作为判断阈值，大于等于该阈值记为1，小于该阈值记为0；统计出所有值为1处的个数，其与信号周期长度的比值即为脉冲的占空比；

由以下关系来确定高频脉冲起始位置ds和终止位置de：

ds:x[n-1]≤x[n]x[n]≥x[n+1]

de:x[n-1]≥x[n]x[n]≤x[n+1]

(6)三角波线性调频的参数提取过程与锯齿波相同，只是在求取调制斜率这一步骤上，只取信号长度的一半，后一半的斜率则取其相反数即可。