[0026] 下面具体说明本发明的实施。
[0027] 步骤1.当ISAR发射宽带线性调频LFM脉冲串信号时,将ISAR接收到的目标回波信号混频,去除载频;
[0028] 雷达回波经过混频后的基带信号可以写成:
[0029]
[0030] 式中,i表示第i个目标散射点;m表示第m个脉冲;Ai为目标散射点的散射强度;p(·)函数表示信号复包络;为快时间;tm为慢时间;Ri(tm)为目标离雷达的距离函数;c为光速;fc为载波频率。
[0031] 步骤2.对混频后的基带信号在快时间维进行傅里叶变换,得到基带频域信号。假设成像积累时间内接收了M个回波,每个回波采样点数为N,则形成M*N维的基带频域信号矩阵。此时目标的距离像已经可以得到,但是在低信噪比下,距离像可能会淹没在噪声中无法观察,也无法进行后续的目标运动参数估计和成像;
[0032] 混频后的信号对快时间 做傅里叶变换得到的基带频域信号可以表示为:
[0033]
[0034] 式中,f为频率;P(f)为信号复包络的频谱函数。
[0035] 仅考虑二次项的目标距离为:
[0036]
[0037] 式中,Ri0为第i个散射点的初始距离;vi为第i个散射点的初始速度;ai为第i个散射点的加速度。
[0038] 将上式代入基带频域信号,可得:
[0039]
[0040] 上式的第三个相位项是由加速度引起的高次项,在少数几次回波间隔时间内可以认为其与f无关,其对于信号的影响体现在初始相位上。第二个相位项包含了多普勒频率,但是多普勒频率随f而变化,可以采用keystone变换来消除这种影响。
[0041] 步骤3.根据目标的速度范围,选择基带频域信号矩阵的前L行,对该L行的数据进行keystone变换,则前L个回波包络走动被校正。对keystone变换后的L*N维数据矩阵,按行(快时间维)逆傅里叶变换得到对齐的目标一维距离像。L同时表示了回波数,根据目标的速度范围选取L属于本领域的常识,不再详述其原理,本发明中优选L不小于30,但是也不能取太大,尤其是目标在成像时间内速度变化时,因为取太大时,目标在L个回波时间内的运动速度不能近似为恒定,建议L小于100。
[0042] 对于Keystone变换,可以指定新的时间变量:
[0043]
[0044] 将上式代入基带频域信号,并忽略二次项,可以得到时间尺度变换后的信号为:
[0045]
[0046] 从上式可以看出,keystone变换利用线性坐标变换消除了多普勒频率和f的耦合,校正了目标的距离走动。由于keystone变换是时间坐标的线性变换,高斯白噪声通过keystone变换后统计特性不变。
[0047] 步骤4.将步骤3得到的对齐的L*N维距离像数据矩阵排列为一维向量,具体的排列方法为:将该矩阵的每一列取出首尾相接,形成一个一维向量;
[0048] 矩阵重排的方法如下所示。
[0049]
[0050] 步骤5.将步骤4得到的一维向量求模后,输入耦合Duffing振子降噪。耦合Duffing振子具体表示为如下二阶常微分方程组:
[0051]
[0052] 式中,x1是第一振子的状态变量,x2是第二振子的状态变量,ξ是振子的阻尼系数,k、q分别为振子间的耦合系数,fcos(t)是周期驱动力,f是周期驱动力的幅度,s(t)是待去噪的输入信号。振子参数分别取:ξ=0.7、k=10、q=10、f=1.7。振子状态变量的初始值可以任取。该常微分方程组可由定步长四阶龙格‑库塔法求解。该系统的降噪输出是振子状态变量的差值,如:x1‑x2,该输出即为输入信号s(t)的去噪结果;
[0053] 步骤6.将步骤5得到的降噪后的一维向量重新恢复为L*N维的矩阵,其恢复方法是步骤4的反向过程;
[0054] 步骤7.将恢复后的L*N维矩阵的所有行相加,得到一个一维向量,并将该一维向量作为一个新的(M‑L+1)*N维观测图像的第一行;
[0055] 步骤8.将基带频域信号矩阵的第2至第L+1行取出,重复步骤3~7的操作,得到(M‑L+1)*N维观测图像的第二行。以此类推,可以得到整个(M‑L+1)*N维观测图像,该观测图像就是增强后的ISAR距离像。
[0056] 图1为本发明的方法流程图。
[0057] 本发明通过仿真算例来显示本发明的实施,图2为仿真算例的飞机目标几何结构。该仿真为一个距离雷达50公里的飞机目标,飞机的运动参数为:径向初始速度‑300m/s(朝
2
向雷达运动)、加速度150m/s 、绕几何中心的转动速度为4°/s。雷达参数为:雷达发射LFM脉冲串,脉冲宽度为5μs,带宽500MHz,载频为10GHz,快时间采样点数2000点,脉冲重复周期
500μs,成像时间为1s,共2000个回波,回波信噪比设置为‑23dB。
[0058] 图3为没有噪声的情况下脉冲压缩后的距离像,该图是从2000个回波中等间隔抽取出200个生成的,从中可以看出目标运动引起的距离像走动和弯曲。
[0059] 图4为混入噪声后的距离像,该图也是从2000个回波中等间隔抽取出200个生成的,可以看出距离像很模糊。
[0060] 图5为混沌振子增强后的距离像,该图也是从2000个回波中等间隔抽取出200个生成的,可以看出距离像变清晰了。在该仿真中,步骤3中的L取的是50。评估图4和图5的图像熵,可以得出图5的信噪比比图4提升了6dB左右。
[0061] 以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,将该方法的修改、变化或延伸到其它应用场景,都应涵盖在本发明的包含范围之内。