CN107301432B - 基于时频分析的自适应辐射源调制识别方法 - Google Patents

Abstract

本发明提供的是一种基于时频分析的自适应辐射源调制识别方法。一，对接收的辐射源信号使用时频分布进行时频分析，将辐射源信号从时域信号转换为时间‑频率二维图像；二，使用图像处理技术降低计算复杂度和特征维数，经过归一化、二值化、图像细化图像预处理操作，增强信号特征信息在图像中的比重；三，联合二阶四阶矩估计方法，利用自适应主成分分析算法，对预处理后的图像进行图像形状特征提取；四，使用LIBSVM分类器识别辐射源信号的调制方式。本发明即能有效避免低信噪比信号特征缺失，又能避免高信噪比信号特征冗余的情况，同时不影响调制识别率。

Description

基于时频分析的自适应辐射源调制识别方法

技术领域

本发明涉及的是一种辐射源信号脉内调制识别方法。

背景技术

辐射源信号脉内调制识别是辐射源的重要特征之一，是对辐射源进行接收、分选和定位等非合作通信方式时的重要参数。电磁威胁环境的信号密度高达每秒120万个脉冲以上，雷达工作频率覆盖范围达到了0.1～20GHz，并正在向0.05～140GHz扩展，雷达信号波形在时频等多个域中同时变化。辐射源信号的调制识别面临着电磁环境复杂化、调制样式多样化、接收信噪比低等问题。

近年来关于雷达辐射源信号识别的大量研究成果表明，脉内调制识别将有望在新型复杂体制雷达辐射源信号识别技术与装备方面获得重要突破。现有的雷达辐射源信号脉内调制识别方法主要有小波及小波包变换法、小波脊线法、相像系数法、经验模式分解法、相对无模糊相位重构法、熵特征法、围线积分双谱法、复杂度特征法、包络法、分形维数法、模糊函数主脊切面法、瞬时频率派生特征法和模糊函数主脊切面特征法等。现有的方法都是对所有不同信噪比情况进行统一处理，为使不同调制信号之间差异明显将采样信号做某种变换，再使用其相应的特征提取方法对变换后的信号进行特征提取，进而完成对信号的调制方式的识别。这些方法在选用信号脉内特征时，都按照同一特征提取方式提取在不同噪声情况的信号脉内特征，没有探讨高信噪比和低信噪比时对采样信号做相应变换，噪声能量对信号调制特征提取的影响。

发明内容

本发明的目的在于提供一种能够解决提取信号调制特征时高信噪比特征冗余、低信噪比特征缺失问题的基于时频分析的自适应辐射源调制识别方法。

本发明的目的是这样实现的：

首先，对接收的辐射源信号使用时频分布进行时频分析，将辐射源信号从时域信号转换为时间-频率二维图像；

然后使用图像处理技术降低计算复杂度和特征维数，经过归一化、二值化、图像细化图像预处理操作，增强信号特征信息在图像中的比重；

随后，联合二阶四阶矩估计方法，利用自适应主成分分析算法，对预处理后的图像进行图像形状特征提取；

最后，使用LIBSVM分类器识别辐射源信号的调制方式。

本发明的优点在于提出了一种运用主成分分析算法对雷达辐射源信号时频分析图像进行特征提取，通过自适应主成分分析算法提取不同信噪比下关键特征，在信噪比为0dB时，基于自适应PCA识别算法识别率达到98％以上，较之Hu矩和伪Zernike矩有12dB左右的提升，同时自适应PCA处于PCA-15和PCA-20之间，表明自适应PCA即能有效避免低信噪比信号特征缺失，又能避免高信噪比信号特征冗余的情况，同时不影响调制识别率。

附图说明

图1是本发明信号调制识别的流程图。

图2是本发明图像中任一像素P的八邻域图。

图3是本发明二阶四阶矩估计信噪比均值示意图。

图4是本发明自适应PCA及不同PCA特征提取识别结果示意图。

图5是本发明自适应PCA及不变矩特征提取识别结果示意图。

具体实施方式

下面举例对本发明做更详细的描述。

步骤1：选用Choi-Williams分布(CWD)对辐射源信号s(t)进行时频分析。对信号瞬时相关函数R(t,τ)的时间变量t做傅立叶逆变换，信号s(t)的模糊函数χ_s(τ,ξ)的定义为式(1)，

式中，τ为时移，ξ为频移。

选取

为时频分布的核函数，式中，τ为时移，ξ为频移，σ为常数(σ越大，自项的分辨率越高；σ越小，对交叉项的抑制越大，实际应用中视信号特点取自项分辨率和交叉项抑制的折中)。将模糊函数χ_s(τ,ξ)和核函数φ(τ,ξ)代入Cohen类时频分布统一表示形式：

得Choi-Williams分布，即

式中，t是时间，ω是频率，τ是时移，μ是积分σ变量。

引入指数核函数的CWD分布可有效抑制模数平面中的由互模数函数形成的交叉项，辐射源信号经过CWD分布时频变换后得到时频图像。

步骤2：对步骤1得到的时频图像进行图像预处理。

各个信号的时频图像灰度值的动态范围是不一样的，像素点灰度值大的数据对分类识别有着大的影响，为了减少数据间的不平衡性，需要对时频图像灰度值进行归一化。设

和σ²分别是所有像素点灰度值的均值和方差。

归一化后像素点的灰度值为：

采用最大方差比的阈值确定方法，提高图像二值化的准确度。

使用Hilditch经典细化算法对二值图像进行细化。目标图像任一像素P相邻8个像素的八邻域如图2所示，若目标像素同时满足以下条件即可删除像素：

i.P₀+P₂+P₄+P₆≤3；

ii.N_C＝1；

iii.P_K(0≤k≤7)中至少有一个目标像素为1；

iv.P₂＝1或N_C2＝1，N_C2为假定P₂＝0时P的联结数；

v.P₄＝1或N_C4＝1，N_C4为假定P₄＝0时P的联结数；

该细化算法每次扫描删除图像上目标的轮廓像素(满足条件i至v)，直到图像上不存在可删除的轮廓像素为止，最终得到辐射源信号的时频骨架图像。

步骤3：运用二阶四阶矩估计法对接收到的辐射源信号x(n)＝s(n)+w(n)进行信噪比估计。其中s(n)＝Ae^jφ(n)为辐射源信号，w(n)是均值为0，方差为σ²的复高斯白噪声，雷达信号和高斯白噪声之间相互独立。

用接收信号序列的时间平均来近似二阶量和四阶量，其中辐射源信号x(n)的二阶矩和四阶矩近似表达式如下：

联立由上两式，解得A²、σ²的估计值为：

定义辐射源信号信噪比估计值为

联立由上两式解得：

步骤4：辐射源信号经过步骤1得到的时频图像，经过步骤2的预处理后可得到辐射源信号的时频骨架图像，将该图像作为原始样本构成一个数据矩阵为：

其协方差矩阵为R＝XX^T，可对该协方差矩阵作特征值分解：

R_M×M＝U∧U^T (12)

其中，T表示转置，∧为协方差矩阵的特征值对角阵，U为相应的特征矩阵，对时频图像作如下变换：

P_M×N＝U^TX＝[p₁,p₂,…,p_M]^T (13)

其中，P为时频图像二值矩阵的主成分，p₁是第一主成分，p_j为第j主成分，选取前k个主成分，构成时频图像的骨架特征矩阵。

步骤5：在低信噪比条件下，由于噪声的影响图片上的骨架特征也会掺入杂质，而对于高信噪比过多的提取成分分量，会弱化了区分性能强的成分分量。依据步骤3解得辐射源信号信噪比估计值

实现在不同信噪比下选取不同的主成分分量的自适应主成分分析算法。

步骤6：使用台湾大学林智仁(Lin Chih-Jen)教授等开发设计的一个简单、易于使用和快速有效的SVM模式识别与回归的软件包LIBSVM，采用交叉验证的方法确定支持向量机本身的参数调节过程。选取支持向量机分类器的核函数，并根据训练样本的形态特征确定核函数的参数，以训练支持向量机分类器。其中，支持向量机分类器的核函数为高斯核函数，即公式(14)：

k(x₁,x₂)＝exp(-γ*|x₁-x₂|²) (14)

具体而言，选择上述的高斯核函数作为核函数，核函数参数的选择采用交叉验证的方法。其中，C(惩罚因子)从2^-5～2⁵(成倍增加)，γ(核函数参数)从2^-5～2⁵(成倍增加)，对于某个(C,γ)，选择5折交叉验证，即将训练样本分成5份，挑选1份作为测试样本，剩下4份作为训练样本，重复5次，计算平均正确率，以此衡量当前(C,γ)的性能，选择性能相对最好的(C,γ)作为最终的参数。

最后，将步骤5在不同信噪比条件下，自适应得到的时频图像特征向量输入训练好的支持向量机，得到输入该系统的辐射源信号脉内调制方式。

实施例

1、实验参数设置如下：对LFM、BPSK、COSTAS、FRANK四种调制方式的雷达辐射源信号进行仿真分析，其中载频为20MHz，采样频率为100MHz，脉冲宽度为10us，BPSK为7位Barker码。对于每一种待分类信号，在信噪比为-2～20dB之间，每隔2dB进行仿真试验，每种信号在当前信噪比下产生200个样本，其中每个信号随机选取140个样本进行SVM的训练，剩余的60个样本用于SVM分类器识别率的测试。试验结果如图3至5所示。

其中图3表示辐射源信噪比估计值的确定，根据步骤3估计AWGN信道下，以BPSK调制信号为例，待估计信噪比范围设为[-10，20]dB，得到二阶四阶矩估计法的信噪比估计均值。

其中图4表示主成分分析选取不同数量特征以LIBSVM为分类器的识别结果，基于PCA-10(选取10个PCA特征)、PCA-15(选取15个PCA特征)、PCA-20(选取20个PCA特征)及自适应PCA这四种不同特征提取情况下的辐射源信号调制识别率变化曲线。

其中图5表示自适应主成分分析及不变矩特征以LIBSVM为分类器的识别结果，基于自适应PCA、Hu矩、伪Zernike矩这三种不同特征提取情况下的辐射源信号调制识别率变化曲线。

2、实验及结果分析

仿真结果表明，使用相同结构的LIBSVM分类器情况下，基于PCA提取的图像特征、基于Hu矩提取的图像特征和基于伪Zernike矩提取的图像特征在信噪比大于0dB的条件下识别效果都较好，基本达到了90％以上。但是另外两种特征提取方法识别率明显低于基于PCA的特征提取方法，同时，使用主成分分量特征越多，识别率越高，基于信噪比的自适应PCA特征提取的识别率在各个信噪比情况下都有较高的识别率，达到了97％以上。由于Hu矩和伪Zernike矩都是提取图像中的不变矩特征，将不变矩特征应用于基于时频分析的信号调制识别中，得到的识别效果明显弱于基于PCA图像特征的调制识别。PCA-15、PCA-20相比，在高信噪比情况下选用PCA-20明显出现了特征冗余，PCA-10、PCA-15、PCA-20相比，在低信噪比情况下选用PCA-10、PCA-15明显出现了特征缺失，基于信噪比的自适应PCA特征提取实现了在对应信噪比选取合适数量的PCA特征，达到了较高的识别率。

本发明提供了一种基于不同信噪比自适应提取PCA特征的调制识别新方法。具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式。本实施例中未明确的各组成部分均可用现有技术加以实现。

Claims (1)

1.一种基于时频分析的自适应辐射源调制识别方法，包括如下步骤：

步骤一，对接收的辐射源信号使用时频分布进行时频分析，将辐射源信号从时域信号转换为时间-频率二维图像；

步骤二，使用图像处理技术降低计算复杂度和特征维数，经过归一化、二值化、图像细化图像预处理操作，增强信号特征信息在图像中的比重；

步骤三，联合二阶四阶矩估计方法，利用自适应主成分分析算法，对预处理后的图像进行图像形状特征提取；

步骤四，使用LIBSVM分类器识别辐射源信号的调制方式；

其特征是：

步骤一具体包括：

选用Choi-Williams分布对辐射源信号s(t)进行时频分析，对信号瞬时相关函数R(t,τ)的时间变量t做傅立叶逆变换，信号s(t)的模糊函数χ_s(τ,ξ)的定义为式(1)，

式中，τ为时移，ξ为频移；

选取

为时频分布的核函数，式中，τ为时移，ξ为频移，σ为常数；将模糊函数χ_s(τ,ξ)和核函数φ(τ,ξ)代入Cohen类时频分布统一表示形式：

得Choi-Williams分布，即

式中，t是时间，ω是频率，τ是时移，μ是积分σ变量；

辐射源信号经过CWD分布时频变换后得到时频图像；

步骤二具体包括：

对时频图像灰度值进行归一化，设

和σ²分别是所有像素点灰度值的均值和方差，

归一化后像素点的灰度值为：

采用最大方差比的阈值确定方法，提高图像二值化的准确度，

使用Hilditch经典细化算法对二值图像进行细化，若目标像素同时满足以下条件即删除像素：

i.P₀+P₂+P₄+P₆≤3；

ii.N_C＝1；

iii.P_K,0≤k≤7中至少有一个目标像素为1；

iv.P₂＝1或N_C2＝1，N_C2为假定P₂＝0时P的联结数；

v.P₄＝1或N_C4＝1，N_C4为假定P₄＝0时P的联结数；

该细化算法每次扫描删除图像上目标的轮廓像素，直到图像上不存在可删除的轮廓像素为止，最终得到辐射源信号的时频骨架图像；

步骤三具体包括：

运用二阶四阶矩估计法对接收到的辐射源信号x(n)＝s(n)+w(n)进行信噪比估计，其中s(n)＝Ae^jφ(n)为辐射源信号，w(n)是均值为0，方差为σ²的复高斯白噪声，雷达信号和高斯白噪声之间相互独立；

用接收信号序列的时间平均来近似二阶量和四阶量，其中辐射源信号x(n)的二阶矩和四阶矩近似表达式如下：

联立由上两式，解得A²、σ²的估计值为：

定义辐射源信号信噪比估计值为

联立由上两式解得：

辐射源信号经过步骤1得到的时频图像，经过步骤2的预处理后得到辐射源信号的时频骨架图像，将该图像作为原始样本构成一个数据矩阵为：

其协方差矩阵为R＝XX^T，对该协方差矩阵作特征值分解：

R_M×M＝U∧U^T (12)

其中，T表示转置，∧为协方差矩阵的特征值对角阵，U为相应的特征矩阵，对时频图像作如下变换：

P_M×N＝U^TX＝[p₁,p₂,…,p_M]^T (13)

其中，P为时频图像二值矩阵的主成分，p₁是第一主成分，p_j为第j主成分，选取前k个主成分，构成时频图像的骨架特征矩阵；

依据步骤3解得辐射源信号信噪比估计值

实现在不同信噪比下选取不同的主成分分量的自适应主成分分析算法；

步骤四具体包括：

选取支持向量机分类器的核函数，并根据训练样本的形态特征确定核函数的参数，以训练支持向量机分类器，其中，支持向量机分类器的核函数为高斯核函数，即公式(14)：

k(x₁,x₂)＝exp(-γ*|x₁-x₂|²) (14)

具体而言，选择上述的高斯核函数作为核函数，核函数参数的选择采用交叉验证的方法，其中，惩罚因子C从2^-5～2⁵成倍增加，核函数参数γ从2^-5～2⁵成倍增加，对于某个(C,γ)，选择5折交叉验证，即将训练样本分成5份，挑选1份作为测试样本，剩下4份作为训练样本，重复5次，计算平均正确率，以此衡量当前(C,γ)的性能，选择性能相对最好的(C,γ)作为最终的参数；

最后，将在不同信噪比条件下，自适应得到的时频图像特征向量输入训练好的支持向量机，得到输入的辐射源信号脉内调制方式。

Images