[发明专利]一种基于时频图像特征的脉内调制识别方法

权利要求书

1.一种基于时频图像特征的脉内调制识别方法，其特征在于，包括如下步骤：

(1)对雷达信号s(t)进行Wigner-Vill分布时频变换，其中t为时间，得到反映信号能量随时间和频率的时频分布图；

对相关函数作滑窗处理，得到时变的局部相关函数公式，

当窗函数取时间冲击函数δ(t)，不加限制，而在时域取瞬时值，

对时变局部相关函数作Fourier变换，即可得到WVD，

设有a个单分量信号s_ε,ε＝1,2,…,a，得到多分量信号的WVD，

雷达信号经过Wigner-Vill分布时频变换后得到时频图像；

(2)对步骤(1)得到的时频图像进行图像预处理；

对时频图像灰度值进行归一化；设和σ²分别是所有像素点灰度值的均值和方差；

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

采用自适应维纳滤波器初步去除时频图像的噪声点，对图像进行增强；采用形态学方法对灰度图像进行处理；形态学处理是应用具有一定形态的结构元素对灰度图像进行腐蚀和膨胀的操作，其中使用结构元素α_b对输入元素F的灰度膨胀记为

选择半径为3的菱型结构元素对时频图像进行闭运算；

通过检测信号时频分布的起止频率，将没有信号分布的图像区域剪切掉；

采用最近邻域插值法归一化时频图像的纵横比；

将原图定义为一个二维函数F(x,y)，(x,y)点的函数值对应像素点的灰度值，f(x′,y′)为归一化纵横比后的目标图；

Fw、Fh设为原图的宽度和高度；fw、fh设为目标图的宽度和高度；原图坐标(x,y)和目标图坐标(x′,y′)之间公式为：

(3)对预处理后的图像进行插值打散原轨迹；

(4)将步骤(3)中被打散成单个点的轨迹类比成目标跟踪中的目标位置，并引入联合概率数据关联来识别出不同雷达信号的轨迹；

引入确认矩阵；确认矩阵定义为：

Ω＝(ω_jb),j＝1,2,…,m_k；b＝0,1,…,B

其中ω_jb表示观测j是否落入目标b的确认门之内；b＝0表示没有目标，此时对应的Ω的列的元素都是1；设在杂波环境中已有B个目标，则他们的状态方程和测量方程分别表示为：

JPDA的状态更新方程为：

其中，表示k时刻目标b的状态矢量，表示目标b的状态预测矢量，V′(k)为联合新息：

其中，β_jb表示目标b与观测j的关联概率，V′_j(k)表示滤波新息；

关联概率表示为：

其中，表示在可行事件θ^j中，观测j是否源于目标b；如果源于目标则等于1，否则为0；在k时刻联合事件θ的条件概率表示为

其中，σ_b为目标检测指示器，如果目标在θ中与观测关联，则σ_b的值等于1，反之为零；Φ是虚假观测事件的个数，表示目标b的检测概率，V表示航迹有效门体积；

(5)运用中心矩提取时频图像的形状特征；

设大小为M*N的二值目标图像表示为：f(x,y)∈{1,0}，则图像的(p+q)阶原点矩m_pq定义为：

图像的(p+q)阶中心矩μ_pq定义为：

其中p＝1且q＝0时原点矩为m₁₀，p＝0且q＝0时原点矩为m₀₀，则表示水平方向上的质心；p＝0且q＝1时原点矩为m₀₁，则表示垂直方向的质心；

其中p＝0且q＝2时μ₀₂表示图像在垂直方向上的伸展度；p＝2且q＝0时μ₂₀表示图像在水平方向上的伸展度；p＝1且q＝1时μ₁₁表示图像的倾斜度；p＝0且q＝3时μ₀₃表示图像在垂直方向上的重心偏移度；p＝1且q＝2时μ₁₂表示图像垂直伸展的均衡程度；由于时频图的差异性主要体现在垂直方向的频域上，故舍弃μ₀₂，而构造出特征向量；

(6)选取支持向量机分类器的核函数，并根据训练样本的形态特征确定核函数的参数，以训练支持向量机分类器；

其中，支持向量机分类器的核函数为高斯核函数，高斯核函数为：h(x₁,x₂)＝exp(-γ*|x₁-x₂|²)选择任意两类设计1个二分类SVM，总共设计n(n-1)/2个二分类SVM；然后将这个二分类SVM构成一个有向无循环图，有向无循环图有n(n-1)/2个非叶子节点和n个叶子节点，每个非叶子节点表示一个二分类SVM分类器，并与下一层的两个节点相连，每个叶子节点表示一个输出；

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