一种基于分数阶傅里叶变换域模态分解的高分辨ISAR成像方法

本发明属于isar成像，特别是涉及一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法。

背景技术：

1、逆合成孔径雷达(isar)作为一种新体制雷达，是一种重要的目标特征提取工具，能够全天候、全方位的对目标区域进行探测搜索，同时可以对飞机、导弹、舰船、卫星等运动目标进行成像以及特征提取和属性区分，在战略防御、战术武器、反卫星等军事领域以及未来的空中、空间交通管制等民用领域中都有重要的应用价值。

2、在isar成像中，最为普遍的成像方法是基于傅里叶变换的距离-多普勒成像方法。该方法要求运动目标各散射点回波的多普勒频率是恒定的或经补偿后是恒定的，才能得到清晰的成像。但在实际情况下，目标的运动状态常伴有机动性，各散射点回波的多普勒频率通常是时变的，导致基于传统傅里叶变换的距离-多普勒成像方法得到的图像非常模糊，无法用于目标识别。

3、一种有效的解决办法就是引入时频分析，即对回波信号进行时频分析，得到每一时刻散射点的瞬时多普勒频率，此即距离-瞬时多普勒成像方法。常见的时频分析方法有短时傅里叶变换(stft)和维纳-魏格尔分布(wvd)。前者可以理解为信号时域加窗后的傅里叶变换，其时频分辨率受制于窗函数，时频分辨率不能同时提升。当提高时间分辨率时，其频率分辨率就会降低；而当提高频率分辨率时，其时间分辨率就会降低。后者具有较高的时频分辨率，但其存在交叉干扰项，限制了实际应用。在wvd基础上，一种常用的抑制交叉项的方法便是伪魏格纳分布(pwvd)，但其是以牺牲时频分辨率来换取交叉项的抑制。鉴于此，本发明将提出一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法。

4、为了便于分析，首先介绍传统距离-多普勒isar成像方法的基本过程。传统距离-多普勒isar成像方法是利用傅里叶变换(ft)获取散射点回波的多普勒信息，进而得到距离-多普勒图像。具体地说，对已经经过距离跟踪和多普勒跟踪后数据g＝(gmn)m×n(其中m表示距离单元，n表示脉冲单元)不同距离下l的对脉冲单元进行傅里叶变换进而可以得到一个二维的图像：r＝(rmω)m×n，其中对于任意1≤m≤m，有(rmω)1×n＝fft((gmn)1×n)，fft表示快速傅里叶变换，ω代表多普勒频移。相比较传统的方法，在基于时频分析(tfa)的isar成像中，用tfa替换ft处理数据g＝(gmn)m×n，具体是对不同距离m下做关于脉冲单元n的时频分析，并得到一个三维的距离-时间-多普勒立方图：q＝(qmnω)m×n×n，其中对于任意1≤m≤m，对数据(gmn)1×n进行时频分析，这里可以得到一组时间-频率二维数据(qmnω)1×n×n。

5、进一步地，还需引入离散分数阶傅里叶变换的定义。任意能量有限时间序列的分数阶傅里叶变换定义为：

6、

7、式中，表示分数阶傅里叶变换算子，核函数kα(u,n)的表示式为

8、

9、式中，k∈z，α表示分数阶傅里叶变换的角度，变量u通常称为分数阶频率，其所在的坐标轴通常称为分数阶傅里叶变换域(简称分数域)。相应地，分数阶傅里叶变换的逆变换的公式为：

10、

11、式中，上标符号*表示共轭运算。特别地，当α＝π/2时，分数阶傅里叶变换便退化为传统傅里叶变换。

12、而符号函数表示如下，

13、

14、此外对于能量有限时间序列它的wvd可以表示如下，

15、

16、这里

17、基于此，下面将详细介绍本发明所提的一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法。

技术实现思路

1、本发明目的是为了解决现有技术中的问题，提供了一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法。

2、本发明是通过以下技术方案实现的，本发明提出一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法，所述方法包括如下步骤：

3、步骤s1，获取经过距离跟踪和多普勒跟踪的isar回波数据g＝(gmn)m×n；初始化循环变量s＝0；

4、步骤s2，更新循环距离变量s＝s+1，对一个距离数值s下的数据(gsn)1×n进行分数阶傅里叶变换域模态分解，得到一组在分数域带限的模态分量{ck(n)}＝{c1(n),c2(n)···ck(n)}，k是模态分量数量，且对于任意1≤n≤n满足

5、步骤s3，对{ck(n)}求取它其中各个元素的wvd，对任意的1≤k≤k，有ck(n)的wvd为wvd(ck(n))＝wck(n,ω)；

6、步骤s4，将{wck(n,ω)}＝{wc1(n,ω),wc2(n,ω)···wck(n,ω)}中元素进行合成，得到距离单元为s下的时频分析结果(qsnω)1×n×n，其中对于任意的1≤n≤n和1≤ω≤n，

7、步骤s5，判断s≥m，若是，则转到下一步骤，否则转到步骤s2；

8、步骤s6，输出距离-时间-多普勒立方图q＝(qmnω)m×n×n,在时间上的一处进行采样并输出二维图像。

9、进一步地，所述分数阶傅里叶变换域模态分解具体实现过程如下:

10、步骤s2.1，对待分解信号f(n)＝(gsn)1×n进行关于变量n的离散分数阶傅里叶变换，并将其分数阶傅里叶变换记为fα(u)；初始化参数模态数量k，离散分数阶傅里叶变换的角度α和惩罚项权重η，模态分量的离散分数傅里叶变换中心分数阶频率拉格朗日乘子的分数阶傅里叶变换及计数参数l＝k，更新步长ρ，循环次数截止条件υ，误差截止条件ε和τ，最终输出本征模态分量{ck(n)}＝{c1(n),c2(n)···ck(n)}，其中1≤n≤n；

11、步骤s2.2，更新外循环计数参数l＝l-1，周期计数参量γ,k＝0，并依据信息熵的计算更新α的值，其中并将α赋给

12、步骤s2.3，更新内循环周期数γ＝γ+1；

13、步骤s2.4，更新k＝k+1；

14、步骤s2.5，更新

15、步骤s2.6，更新

16、步骤s2.7，比较k和l+1的值的大小，当k小于l+1时，说明没有完成所有模态的分解此时转到步骤s2.4，否则进行下一步；

17、步骤s2.8，更新

18、步骤s2.9，重复步骤s2.4至步骤s2.8直到满足截止条件：

19、

20、步骤s2.10，更新k＝0；

21、步骤s2.11，根据式求取其能量聚集的角度

22、步骤s2.12，比较与τ，如果时，则将γ＝0，并回到步骤s2.3，否则进行下一步；

23、步骤s2.13，将的离散分数阶傅里叶逆变换输出为ck-l(n)，此时如果l＞1将信号更新为输入到步骤s2.2中，否则将的离散分数阶傅里叶逆变换输出为ck(n)，完成分解。

24、本发明的有益效果为：

25、本发明提出了一种基于分数阶傅里叶变换域模态分解的高分辨isar成像方法。该方法首先对运动目标散射点回波中每个距离单元的信号进行分数阶傅里叶变换域模态分解，得到信号各个模态分量；然后对每个模态分量进行魏格纳-维尔分布，并将各个模态分量的魏格纳-维尔分布累加，最终得到一个时间-距离-瞬时多普勒立方图，那么抽取时间上的一帧即可得到目标的距离-瞬时多普勒图像。与传统isar成像方法相比，本发明方法构造的时频分析不存在交叉项，且具有较高的时频分辨率，有利于进一步提升isar成像性能。