一种飞机机翼ISAR回波振动补偿与整体聚焦成像方法

一种飞机机翼isar回波振动补偿与整体聚焦成像方法
技术领域
1.本发明涉及isar成像技术领域，具体涉及一种飞机机翼isar回波振动补偿与整体聚焦成像方法。


背景技术：

2.逆合成孔径雷达(isar)具有全天候、全天时获取远距离运动目标高分辨率图像的能力，被广泛应用于军事和民用领域。isar成像通过发射宽带信号以及脉冲压缩技术获得距离向高分辨率，同时利用目标与雷达的相对姿态变化来实现方位高分辨成像，实现非合作目标的二维成像。isar图像中包含目标形状、尺寸、姿态等信息，为非合作目标的识别提供重要支撑。然而在气流影响下机翼存在振动的飞机目标无法视为刚体目标进行isar成像处理，现有的刚体目标成像算法难以获得聚焦良好的飞机目标成像结果，无法为后续目标结构分析、分类和识别等提供支撑。而且，随着分辨率的提升，机翼振动对飞机目标聚焦成像的影响更加明显。
3.在isar成像过程中，主体散射点与振动散射点的多普勒谱存在差异，振动散射点由于振动产生的多普勒表现为随时间的正弦变化曲线，而主体散射点多普勒不随时间变化，表现为直线。有学者通过振动部件和主体部件在时频域的不同表现形式，将振动部件回波与主体回波分离，对主体回波应用刚体目标isar成像方法进行成像处理，对振动散射点回波应用逆radon变换等方法估计振动参数并做成像处理。然而，当距离压缩后的雷达回波信噪比较低时，难以在时频域中提取到散射点的时频变化曲线，而且基于逆radon变换的非刚体目标成像方法对于主体与微动部件是在不同表达域内进行成像的，不利于目标整体结构的分析和后续对目标的识别，使得算法应用具有一定的局限性。


技术实现要素：

4.为了克服现有非刚体目标isar成像方法的不足，本发明利用粒子群优化算法良好的全局优化性能，提出一种飞机机翼isar回波振动分量补偿与整体聚焦成像方法，能精确估计机翼振动散射点的振动参数，对飞机目标isar回波中机翼振动分量进行补偿，从而可以应用刚体目标isar成像方法实现机身和机翼的同时聚焦成像，为目标结构分析和识别提供支撑。
5.本发明的技术思路为：
6.首先对输入的回波数据进行一阶keystone校正，并对二次相位进行调频率补偿，提高飞机主体散射点的成像聚焦度；然后，根据距离单元回波信号，计算相应的方位成像结果的能量和熵值，筛选出含有振动散射点的回波距离单元；接着，对筛选出的含有振动散射点的回波距离单元应用pso算法搜索每个振动散射点的振动参数，直至该距离单元回波剩余能量小于预先设定的门限，并构造相位补偿因子，对距离单元回波进行补偿，分离含正弦调制项的散射点回波数据，实现强散射点的分离与正弦分量的去除；不断迭代搜索过程，直至满足迭代终止条件，输出去除振动分量的距离单元回波；最后，不断迭代循环搜索筛选出
来的距离单元回波，实现整体回波的振动分量去除，实现飞机目标的整体方位聚焦成像。
7.为了达到上述目的，本发明采用以下技术方案予以实现。
8.一种飞机机翼isar回波振动补偿与整体聚焦成像方法，包括以下步骤：
9.步骤1，对回波数据依次进行一阶keystone校正、调频率估计与补偿，得到补偿后的回波数据；
10.步骤2，对所述补偿后的回波数据进行方位压缩并计算所有距离单元回波的能量和熵值，去除能量和熵值较低的距离单元后，筛选出含有振动散射点的距离单元回波；
11.步骤3，从所述含有振动散射点的距离单元回波中选取任一距离单元回波数据，建立以振动参数为优化变量，以振动分量补偿后的距离单元回波方位成像结果的熵值为目标函数的无约束优化求解模型；
12.利用粒子群优化算法估计振动参数，并在方位频域分离该振动参数对应的振动散射点回波信号，通过估计的振动参数构造相位补偿函数，采用所述相位补偿函数对所述振动散射点信号进行补偿，获得去除振动分量后的距离单元回波信号；
13.步骤4，对所述含有振动散射点的距离单元回波中的每个距离单元回波数据迭代应用步骤3，获得飞机目标去除机翼振动分量的整体回波信号；
14.步骤5，对所述飞机目标去除机翼振动分量的整体回波信号进行方位压缩处理，输出飞机目标整体聚焦的isar图像。
15.本发明技术方案的特点和进一步的改进为：
16.(1)步骤1具体包含以下子步骤：
17.子步骤1.1，假设雷达发射线性调频信号，机翼上含有k1个散射点，机身含有k2个散射点，其距离向采样点数为m，方位向脉冲数为n，经过平动补偿和去斜处理后的回波数据s(f
r
,t)为；
[0018][0019]
其中，t为方位慢时间，rect(
·
)为矩形窗函数，σ为后向散射系数，j为虚数单位，f
c
为信号载频，f
r
为距离频率，b为信号带宽，c为电磁波传播速度；r
p
(t)为第p个机翼散射点与雷达的瞬时距离，r
q
(t)为第q个机身散射点与雷达的瞬时距离，表达式分别如下：
[0020][0021]
r
q
(t)＝x
q
cos(ωt+α)cosβ+y
q
sin(ωt+α)cosβ
[0022]
≈(x
q
cosα+y
q
sinα)cosβ+(y
q
cosα
‑
x
q
sinα)cosβωt
[0023]
其中，(x
p
,y
p
)为第p个机翼散射点相对于飞机中心的坐标，α和β分别表示飞机与雷达连线在雷达测量坐标系中的方位角和俯仰角，ω为目标转动角速度，θ0为机翼最大振动角度，f
v
为机翼振动频率，为机翼振动初相，(x
q
,y
q
)为第q个机身散射点相对于飞机中心
的坐标；
[0024]
子步骤1.2，在相干积累时间内，由于转动分量的影响，散射点距离包络跨多个距离单元，使用一阶keystone变换对整体回波进行越距离单元徙动校正：
[0025]
f
c
τ＝(f
c
+f
r
)t
[0026]
其中，τ为虚拟慢时间，经过上述变量代换可完成徙动校正，使得同一主体散射点处于同一距离单元内；
[0027]
子步骤1.3，使用方位调频率估计方法，补偿回波沿方位向的二次相位变化，得到补偿后的回波数据为s：
[0028]
s＝s1(r,τ)+s2(r,τ)
[0029]
其中，s1(r,τ)为调频率补偿之后的机身回波，s2(r,τ)为调频率补偿之后的机翼回波；
[0030]
其中，调频率补偿之后的机身回波为：
[0031][0032]
式中，sinc(
·
)为辛格函数，r为回波矩阵的距离单元的位置，a和b分别表示与散射点有关的距离单元位置和方位单元位置，表达式如下所示：
[0033]
a＝(xcosα+ysinα)cosβ
[0034]
b＝(ycosα
‑
xsinα)cosβω
[0035]
式中，x和y分别表示目标坐标系下散射点的横向坐标和纵向坐标；
[0036]
调频率补偿之后的机翼回波为：
[0037][0038]
其中，f为经过补偿后的正弦相位的变化频率，a为与散射点坐标相关的幅度值，f的表达式如下所示：
[0039][0040]
(2)步骤2包含以下子步骤：
[0041]
子步骤2.1，对所述补偿后的回波数据进行方位傅里叶变换计算，得到方位成像结果q；并按照下述公式计算每个距离单元回波的能量p
i
和熵值e
i
：
[0042][0043][0044]
其中，|
·
|表示取模运算，i＝1,2,
…
,m表示回波距离单元标号，q
i
为1
×
n的向量，表示第i个距离单元回波，q
i
(n)表示q
i
中第n个采样点的值，p
i
和e
i
分别表示第i个距离单元
回波的能量和熵值；
[0045]
子步骤2.2，将所有距离单元回波的能量从低到高排序，设置能量阈值，剔除能量低于所述能量阈值的距离单元回波；
[0046]
子步骤2.3，设置熵值门限，再将熵值低于所述熵值门限的距离单元回波剔除，则剩余的距离单元就是筛选出来的含有振动散射点的距离单元回波。
[0047]
(3)子步骤2.2中，所述能量阈值为：
[0048]
计算所有距离单元回波的能量的平均能量e
a
，所述平均能量e
a
的50％即为能量阈值。
[0049]
(4)步骤3包含以下子步骤：
[0050]
子步骤3.1，将机翼振动散射点的振动参数建模为利用振动参数构造相位补偿因子pha：
[0051][0052]
子步骤3.2，使用相位补偿因子对某一距离单元回波s
i
进行补偿，得到振动分量补偿后的距离单元回波s
′
i
：
[0053]
s
′
i
＝s
i
·
pha
[0054]
子步骤3.3，计算振动分量补偿后的距离单元回波s
′
i
方位成像结果的熵值；
[0055]
以所述振动分量补偿后的距离单元回波s
′
i
方位成像结果的熵值为优化的目标函数，以振动参数为待优化的变量，建立振动参数估计的无约束优化求解模型如下：
[0056][0057]
其中，为振动参数估计值，e
i
(θ)为计算第i个距离单元回波在振动参数θ构造的补偿因子补偿后方位成像结果的熵值；
[0058]
子步骤3.4，采用粒子群优化算法求解所述无约束优化求解模型，估计一个距离单元中最强振动散射点的振动参数；
[0059]
子步骤3.5，利用估计的振动参数构造方位图像域滤波器，将该散射点回波从距离单元回波中分离，得到振动散射点回波信号；
[0060]
子步骤3.6，利用估计的振动参数构造相位补偿函数，采用所述相位补偿函数对所述振动散射点回波信号进行补偿，并累加至输出的去除振动分量后的距离单元回波信号；
[0061]
子步骤3.7，若距离单元回波剩余能量小于门限阈值δ或者迭代次数达到最大迭代次数时，则停止迭代，获得去除振动分量后的距离单元回波信号；否则，跳转至子步骤3.4继续执行。
[0062]
(5)子步骤3.4具体包含以下子步骤：
[0063]
子步骤3.4.1，基于粒子群优化算法构建所述振动参数估计的无约束优化求解模型，在三变量的优化求解问题中，建立包含d个粒子的种群，其中，第m个粒子的位置和速度分别为x
m
＝{x
m1
,x
m2
,x
m3
}和v
m
＝{v
m1
,v
m2
,v
m3
}，每个粒子位置分别表示待估计的机翼振动参数；所有粒子的位置构成解空间x＝{x1,x2,
…
,x
d
}，粒子本身历史最优位置为pbest
m
＝{pbest
m1
,pbest
m2
,pbest
m3
}，种群历史最优位置为gbest＝{gbest1,gbest2,
gbest3}；粒子的速度和位置更新公式如下所示：
[0064][0065][0066]
其中，和分别表示第l次迭代时第m个粒子的速度和位置，w为惯性系数，c1和c2表示加速常数，rand1和rand2表示均匀分布于区间[0,1]的随机数，表示第l次迭代时第m个粒子的个体最优位置，gbest
l
表示第l次迭代时种群的最优位置；
[0067]
子步骤3.4.2，参数初始化：设置最大迭代次数g
max
和当前迭代次数l＝1，确定粒子群规模d，粒子速度边界v
max
；对d个粒子在解空间内随机初始化，设定粒子速度为
‑
v
max
和v
max
区间随机分布，避免初始速度越界；最后，随机初始化每个粒子的个体最优位置pbest
i
以及种群最优位置gbest；
[0068]
子步骤3.4.3，针对每个粒子，计算基于该粒子对应的振动参数相位补偿后的方位成像结果的熵值
[0069]
子步骤3.4.4，更新粒子的pbest
m
和种群的gbest；若粒子当前熵值小于粒子最小熵值，则记录当前熵值并更新为当前位置；比较所有粒子熵值的最小值与全局最优熵值的大小，若则记录其熵值并更新gbest
l+1
为当前最优位置；
[0070]
子步骤3.4.5，更新粒子的速度和位置：在搜索过程中若出现粒子速度越界情况则需进行越界处理；即若v
m
>v
max
，则令v
m
＝v
max
；若v
m
<
‑
v
max
，则令v
m
＝
‑
v
max
；
[0071]
子步骤3.4.6，判断迭代终止条件：若当前迭代次数l>g
max
，则算法终止，输出振动参数估计的最优值否则，跳转至子步骤3.4.4继续执行。
[0072]
(6)步骤4包含以下子步骤：
[0073]
子步骤4.1，用集合t表示所有的筛选出的含有振动散射点的距离单元回波，并将筛选出的含有振动散射点的距离单元回波集合以标号升序排列，并使用新的标号表示为t＝{t1,t2,
…
,t
l
}，并初始化迭代次数k＝1；
[0074]
子步骤4.2，判断距离单元回波数据t
k
是否属于集合t，若属于集合t，则执行子步骤4.3；否则，结束迭代，跳至步骤5；
[0075]
子步骤4.3，利用步骤3获得去除振动分量的距离单元回波信号t
k
′
，替换原始回波对应位置的回波信号，并更新迭代次数，令k增加1，跳至子步骤4.2。
[0076]
(7)步骤5具体为：
[0077]
对所述飞机目标去除机翼振动分量的整体回波信号进行方位傅里叶变换，得到飞机目标整体聚焦的isar图像。
[0078]
与现有技术相比，本发明的有益效果为：
[0079]
(1)本发明方法为：首先对输入的回波数据进行越距离单元徙动校正和调频率补偿，提高飞机主体散射点的成像聚焦度；接着，对飞机目标回波进行方位成像处理，计算所
有距离单元的能量和熵值，进而筛选出含机翼振动散射点的距离单元；然后，在筛选出的距离单元回波中迭代应用粒子群优化方法依次搜索每个振动散射点的振动参数，直至该距离单元回波剩余能量小于预先设定的门限，构造每个散射点的振动相位补偿因子，实现机翼振动散射点正弦调制相位的去除；之后，将上述步骤应用于筛选出的每个距离单元回波，得到机翼回波振动分量补偿后的回波数据；最后，通过方位压缩处理即可得到飞机目标的整体聚焦成像结果。本发明能在飞机含有机翼振动情况下，克服机翼散射点振动引起的图像方位散焦现象，实现飞机整体的聚焦成像，成像质量优于传统刚体目标isar成像方法的效果。
[0080]
(2)本发明中粒子群优化算法模型简单，全局搜索能力强，计算效率高，且在同一个距离单元内，将主体散射点看作振动散射点的特例，实现振动散射点的参数估计和振动分量补偿。
[0081]
(3)振动散射中心参数估计之前利用能量和熵值门限，筛选含有振动散射点的距离单元，计算效率更高。
[0082]
(4)传统方法将微动散射点作为干扰项去除或者将主体与微动散射点分别成像，难以进行目标整体结构的分析。本方法利用粒子群优化方法解决了飞机目标机翼散射点振动运动导致刚体目标成像方法难以实现聚焦成像的问题，通过对飞机目标机翼振动参数的估计和补偿，将飞机目标非刚体运动形式转化为刚体运动，从而实现飞机目标主体和机翼整体聚焦成像。
附图说明
[0083]
下面结合附图和具体实施例对本发明做进一步详细说明。
[0084]
图1为飞机机翼isar回波振动补偿与整体聚焦成像方法的流程图；
[0085]
图2(a)为本发明中基于pso的距离单元振动散射点振动参数估计流程图；图2(b)为本发明中距离单元回波振动分量补偿流程图；图2(c)为本发明中飞机目标整体聚焦成像流程图；
[0086]
图3(a)为在机翼振动情况下刚体目标isar成像方法获得的飞机成像结果图；图3(b)为应用本发明获得的飞机目标聚焦成像结果。
具体实施方式
[0087]
下面将结合实施例对本发明的实施方案进行详细描述，但是本领域的技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限制本发明的范围。
[0088]
参考图1，一种飞机机翼isar回波振动补偿与整体聚焦成像方法，包括以下步骤：
[0089]
步骤1，对回波数据依次进行一阶keystone校正、调频率估计与补偿，得到补偿后的回波数据。
[0090]
具体的，步骤1包含以下子步骤：
[0091]
子步骤1.1，假设雷达发射线性调频信号，机翼上含有k1个散射点，机身含有k2个散射点，其距离向采样点数为m，方位向脉冲数为n，经过平动补偿和去斜处理后的回波数据s(f
r
,t)为；
[0092][0093]
其中，t为方位慢时间，rect(
·
)为矩形窗函数，σ为后向散射系数，j为虚数单位，f
c
为信号载频，f
r
为距离频率，b为信号带宽，c为电磁波传播速度。r
p
(t)为第p个机翼散射点与雷达的瞬时距离，r
q
(t)为第q个机身散射点与雷达的瞬时距离，表达式分别如下：
[0094][0095]
r
q
(t)＝x
q
cos(ωt+α)cosβ+y
q
sin(ωt+α)cosβ
[0096]
≈(x
q
cosα+y
q
sinα)cosβ+(y
q
cosα
‑
x
q
sinα)cosβωt
[0097]
其中，(x
p
,y
p
)为第p个机翼散射点相对于飞机中心的坐标，α和β分别表示飞机与雷达连线在雷达测量坐标系中的方位角和俯仰角，ω为目标转动角速度，θ0为机翼最大振动角度，f
v
为机翼振动频率，为机翼振动初相，(x
q
,y
q
)为第q个机身散射点相对于飞机中心的坐标。
[0098]
子步骤1.2，在相干积累时间内，由于转动分量的影响，散射点距离包络跨多个距离单元，使用一阶keystone变换对整体回波进行越距离单元徙动校正。
[0099]
f
c
τ＝(f
c
+f
r
)t
[0100]
其中，τ为虚拟慢时间，经过上述变量代换可完成徙动校正，使得同一主体散射点处于同一距离单元内。
[0101]
子步骤1.3，使用方位调频率估计方法，补偿回波沿方位向的二次相位变化，使得机身散射点能够聚焦成像。调频率补偿之后的整体回波数据为s。
[0102]
s＝s1+s2[0103]
其中，s1为调频率补偿之后的机身回波，s2为调频率补偿之后的机翼回波。为体现回波与距离方位之间的关系，机身回波数据s1也可表示为s1(r,τ)。
[0104][0105]
其中，sinc(
·
)为辛格函数，r为回波矩阵的距离单元的位置。a和b分别表示与散射点有关的距离单元位置和方位单元位置，表达式如下所示：
[0106]
a＝(xcosα+ysinα)cosβ
[0107]
b＝(ycosα
‑
xsinα)cosβω
[0108]
其中，x和y分别表示目标坐标系下散射点的横向坐标和纵向坐标。
[0109]
机身回波数据s2也可表示为s2(r,τ)。
[0110][0111]
其中，f
v
为机翼振动频率，f为经过补偿后的正弦相位变化频率，a为与散射点坐标相关的幅度值，f的表达式如下所示：
[0112][0113]
步骤2，对所述补偿后的回波数据进行方位压缩并计算所有距离单元回波的能量和熵值，去除能量和熵值较低的距离单元后，筛选出含有振动散射点的距离单元回波。
[0114]
利用补偿后的回波数据的稀疏性和机翼振动散射点的散焦现象，筛选出含有振动散射点的距离单元回波信号，减少算法搜索范围，提高算法的效率，具体包含以下子步骤：
[0115]
子步骤2.1，对补偿后的回波数据s经过方位傅里叶变换(fft)计算，得到方位成像结果q。并按照下述公式计算每个距离单元回波的能量p
i
和熵值e
i
。
[0116][0117][0118]
其中，|
·
|表示取模运算，i＝1,2,
…
,m表示回波距离单元标号，q
i
为1
×
n的向量，表示第i个距离单元回波，q
i
(n)表示q
i
中第n个采样点的值，p
i
和e
i
分别表示第i个距离单元回波的能量和熵值。
[0119]
子步骤2.2，将所有距离单元回波的能量从低到高排序，剔除若干个能量最小的距离单元，一般认为能量低于距离单元平均能量e
a
的50％时认为其不含有散射点。
[0120]
子步骤2.3，将剩余距离单元内挑出若干个熵值最小的距离单元，由于主体经过运动补偿聚焦度良好，熵值较低，而机翼部分由于存在振动分量，成像结果存在严重散焦，因此通过设置距离单元熵值门限，将熵值小于该门限的距离单元剔除，则剩余的距离单元就是筛选出来的含有振动散射点的距离单元。
[0121]
步骤3，从所述含有振动散射点的距离单元回波中选取任一距离单元回波数据，建立以振动参数为优化变量，以振动分量补偿后的距离单元回波方位成像结果的熵值为目标函数的无约束优化求解模型；利用粒子群优化算法估计振动参数，并在方位频域分离该振动参数对应的振动散射点回波信号，通过估计的振动参数构造相位补偿函数，采用所述相位补偿函数对所述振动散射点信号进行补偿，获得去除振动分量后的距离单元回波信号。
[0122]
具体的，参考图2(b)，建立筛选结果中的某一距离单元的熵值函数，将振动散射点多普中的正弦参数对应的isar图像散焦问题转化为正弦参数估计的无约束优化问题，并利用粒子群优化算法(pso)算法搜索振动参数并补偿振动分量，具体包含以下子步骤：
[0123]
子步骤3.1，将机翼振动散射点的振动参数建模为利用振动参数构造相位补偿因子pha。
[0124][0125]
子步骤3.2，使用相位补偿因子对某一距离单元回波s
i
进行补偿，得到振动分量补偿后的距离单元回波s
′
i
。
[0126]
s
′
i
＝s
i
·
pha
[0127]
子步骤3.3，计算振动分量补偿后的距离单元回波s
′
i
方位成像结果的熵值。若正弦参数得到准确估计，则散射点的正弦分量去除，其方位成像聚焦度最高，熵值最小。因此，可以振动补偿后方位成像结果的熵值为优化的目标函数，以振动参数为待优化的变量，建立振动参数估计的无约束优化求解模型如下：
[0128][0129]
其中，为振动参数估计值，e
i
(θ)为计算第i个距离单元回波在参数θ构造的补偿因子补偿后方位成像结果的熵值。
[0130]
子步骤3.4，参考图2(a)，采用pso算法求解上述振动参数估计的无约束优化问题，估计一个距离单元中最强振动散射点的振动参数。具体如下：
[0131]
子步骤3.4.1，构建基于pso的振动参数优化求解模型，在三变量的优化求解问题中，建立包含d个粒子的种群，其中，第m个粒子的位置和速度分别为x
m
＝{x
m1
,x
m2
,x
m3
}和v
m
＝{v
m1
,v
m2
,v
m3
}，每个粒子位置分别表示待估计的机翼振动参数。所有粒子的位置构成解空间x＝{x1,x2,
…
,x
d
}，粒子本身历史最优位置为pbest
m
＝{pbest
m1
,pbest
m2
,pbest
m3
}，种群历史最优位置为gbest＝{gbest1,gbest2,gbest3}。粒子的速度和位置更新公式如下所示：
[0132][0133][0134]
其中，和表示第l次迭代时第m个粒子的速度和位置，w为惯性系数，c1和c2表示加速常数，rand1和rand2表示均匀分布于区间[0,1]的随机数。表示第l次迭代时第m个粒子的个体最优位置，gbest
l
表示第l次迭代时种群的最优位置。
[0135]
子步骤3.4.2，参数初始化：设置最大迭代次数g
max
和当前迭代次数l＝1，确定粒子群规模d，粒子速度边界v
max
。对d个粒子在解空间内随机初始化，设定粒子速度为
‑
v
max
和v
max
区间随机分布，避免初始速度越界。最后，随机初始化每个粒子的个体最优位置pbest
i
以及种群最优位置gbest。
[0136]
子步骤3.4.3，针对每个粒子，计算基于该粒子对应的振动参数相位补偿后的方位成像结果的熵值
[0137]
子步骤3.4.4，更新粒子的pbest
m
和种群的gbest。若粒子当前熵值小于粒子最小熵值，则记录当前熵值并更新为当前位置；比较所有粒子熵值的最小值与全局最
优熵值的大小，若则记录其熵值并更新gbest
l+1
为当前最优位置。
[0138]
子步骤3.4.5，更新粒子的速度和位置：在搜索过程中若出现粒子速度越界情况则需进行越界处理。即若v
m
>v
max
，则令v
m
＝v
max
；若v
m
<
‑
v
max
，则令v
m
＝
‑
v
max
。
[0139]
子步骤3.4.6，判断迭代终止条件：若当前迭代次数l>g
max
，则算法终止，输出振动参数估计的最优值否则，跳转至子步骤3.4.4继续执行。
[0140]
子步骤3.5，利用估计的振动参数构造方位图像域滤波器，将该散射点回波从距离单元回波中分离，得到振动散射点回波信号。
[0141]
子步骤3.6，利用估计的振动参数构造相位补偿函数，对强散射点的回波进行振动分量补偿，并累加至输出的距离单元回波信号。
[0142]
子步骤3.7，若距离单元回波剩余能量小于门限阈值δ或者迭代次数达到最大迭代次数时，则停止迭代，输出振动补偿后的距离单元回波信号，否则，跳转至子步骤3.4继续执行。
[0143]
步骤4，对所述含有振动散射点的距离单元回波中的每个距离单元回波数据迭代应用步骤3，获得飞机目标去除机翼振动分量的整体回波信号。
[0144]
具体的，参考图2(c)，步骤4包含以下子步骤：
[0145]
子步骤4.1，用集合t表示所有的筛选出的距离单元回波，并将筛选出的距离单元回波集合以标号升序排列，并使用新的标号表示为t＝{t1,t2,
…
,t
l
}，并初始化迭代次数k＝1。
[0146]
子步骤4.2，判断距离单元回波数据t
k
是否属于集合t，若属于集合t，则执行子步骤4.3，否则，结束迭代，调至步骤5。
[0147]
子步骤4.3，利用步骤3获得去除振动分量的距离单元回波数据t
k
′
，替换原始回波对应位置的回波数据，并更新迭代次数k＝k+1，调至子步骤4.2。
[0148]
步骤5，对所述飞机目标去除机翼振动分量的整体回波信号进行方位压缩处理，输出飞机目标整体聚焦的isar图像。
[0149]
具体的，对所述飞机目标去除机翼振动分量的整体回波信号进行方位傅里叶变换(fft)计算，得到飞机目标整体聚焦的isar图像。
[0150]
进一步的，以下通过仿真实验对本发明实施例方法上述有益效果进行验证：
[0151]
实验一：
[0152]
1、仿真条件
[0153]
雷达发射信号处于ka频段，带宽为4ghz，距离采样点数为512，积累脉冲数为1024，飞机目标最大振动幅度为1度，振动频率为1hz，振动初相为0度；
[0154]
2、仿真实验内容及结果分析
[0155]
本发明实施例中，pso算法的粒子种群规模为100，最大迭代次数是500次。图3(a)展示了含机翼振动的飞机目标使用刚体isar成像方法的成像结果，可以看出，isar图像中机翼散射点处方位向成像结果出现散焦，聚焦度差。图3(b)为应用本发明处理后的成像结果，消除了机翼方位向散焦的现象，实现了主体与机翼的整体聚焦成像。
[0156]
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过
程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
[0157]
虽然，本说明书中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。