首页 > 专利 > 南京信息职业技术学院 > 一种应用混沌振子的逆合成孔径雷达距离像增强方法专利详情

一种应用混沌振子的逆合成孔径雷达距离像增强方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-04-08
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2019-08-20
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-04-22
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-04-08
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910276584.8 申请日 2019-04-08
公开/公告号 CN110058231B 公开/公告日 2022-04-22
授权日 2022-04-22 预估到期日 2039-04-08
申请年 2019年 公开/公告年 2022年
缴费截止日
分类号 G01S13/90G01S7/41 主分类号 G01S13/90
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 2
权利要求数量 3 非专利引证数量 1
引用专利数量 0 被引证专利数量 0
非专利引证 1、CN 105425238 A,2016.03.23CN 108375416 A,2018.08.07CN 108663605 A,2018.10.16CN 109031299 A,2018.12.18罗文茂.基于双扩展型Duffing振子的脉冲信号检测方法《.系统仿真学报》.2019,全文.;
引用专利 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 南京信息职业技术学院 当前专利权人 南京信息职业技术学院
发明人 罗文茂、姜敏敏、陈雪娇、顾艳华 第一发明人 罗文茂
地址 江苏省南京市栖霞区仙林大学城文澜路99号 邮编 210023
申请人数量 1 发明人数量 4
申请人所在省 江苏省 申请人所在市 江苏省南京市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
南京天翼专利代理有限责任公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
奚铭
摘要
一种基于混沌振子的逆合成孔径雷达距离像增强方法,属于雷达信号处理领域。对于逆合成孔径雷达ISAR来说,在低信噪比情况下,其距离像可能会淹没在噪声中无法观察,也无法进行后续的目标运动参数估计和成像。本发明方法将叠加噪声的距离像信号经过处理后,输入耦合混沌振子系统中,达到抑制噪声、增强距离像信号的目的,具有推广应用的价值。
  • 摘要附图
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
  • 说明书附图:图1
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
  • 说明书附图:图2
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
  • 说明书附图:图3
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
  • 说明书附图:图4
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
  • 说明书附图:图5
    一种应用混沌振子的逆合成孔径雷达距离像增强方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-04-22 授权
2 2019-08-20 实质审查的生效 IPC(主分类): G01S 13/90 专利申请号: 201910276584.8 申请日: 2019.04.08
3 2019-07-26 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种应用混沌振子的逆合成孔径雷达距离像增强方法,其特征是包括如下步骤:
1)当逆合成孔径雷达ISAR发射宽带线性调频LFM脉冲串信号时,将ISAR接收到的目标回波信号混频,去除载频;
2)对去除载频后的混频信号在快时间维进行傅里叶变换,得到基带频域信号,设成像积累时间内接收了M个回波信号,每个回波信号采样点数为N,则形成M*N维的基带频域信号矩阵;
3)根据目标的速度范围,选择基带频域信号矩阵的前L行,对该L行的数据进行keystone变换,则前L个回波包络走动被校正,对keystone变换后得到的L*N维数据矩阵,按行在快时间维逆傅里叶变换得到对齐的目标一维距离像;
4)将步骤3)得到的对齐的L*N维距离像数据矩阵排列为一维向量,具体的排列方法为:
将该矩阵的每一列取出,首尾相接,形成一个一维向量;
5)将步骤4)得到的一维向量求模后,输入耦合Duffing振子降噪,耦合Duffing振子表示为如下二阶常微分方程组:
式中,x1是第一振子的状态变量,x2是第二振子的状态变量,ξ是振子的阻尼系数,k、q分别为振子间的耦合系数,fcos(t)是周期驱动力,f是周期驱动力的幅度,s(t)是待去噪的输入信号,振子状态变量的初始值可以任取,所述常微分方程组可由定步长四阶龙格‑库塔法求解;
该耦合Duffing振子系统的降噪输出是振子状态变量的差值x1‑x2,所述降噪输出即为输入信号s(t)的去噪结果;
6)将步骤5)得到的降噪后的一维向量重新恢复为L*N维的矩阵,其恢复方法是步骤4)的反向过程;
7)将恢复后的L*N维矩阵的所有行相加,得到一个一维向量,并将该一维向量作为一个新的(M‑L+1)*N维观测图像的第一行;
8)将基带频域信号矩阵的第2至第L+1行取出,重复步骤3)~7)的操作,得到(M‑L+1)*N维观测图像的第二行,以此类推,得到整个(M‑L+1)*N维观测图像,该观测图像就是增强后的ISAR距离像。

2.根据权利要求1所述的一种应用混沌振子的逆合成孔径雷达距离像增强方法,其特征是步骤3)中,L的选取不小于30,同时为确保目标在L个回波时间内的运动速度近似为恒定,L的取值小于100。

3.根据权利要求1所述的一种应用混沌振子的逆合成孔径雷达距离像增强方法,其特征是步骤5)中,振子参数分别取:ξ=0.7、k=10、q=10、f=1.7。
说明书

技术领域

[0001] 本发明属于雷达技术领域,涉及雷达信号处理,为一种应用混沌振子的逆合成孔径雷达距离像增强方法。

背景技术

[0002] 为了降低被雷达探测的概率,许多军事飞行器正致力于缩小雷达散射面积,所以逆合成孔径雷达(ISAR)将不可避免需要在低信噪比下成像,这将大大增加ISAR成像的难度。
[0003] 常用的ISAR成像方法有距离多普勒(RD)算法、基于运动参数估计的成像算法两大类。在高信噪比条件下,ISAR成像的方法已经很成熟了。在低信噪比条件下,也提出了一些ISAR运动补偿算法,比如:高阶相位项参数估计实现运动补偿、修正的Keystone变换方法完成包络对齐、图像熵法实现相位自聚焦、联合自聚焦等。但是,现有的低信噪比运动补偿方法的抗噪性能是有限的,信噪比进一步降低,现有的方法都将失效。本发明的目的就是要进一步提升ISAR的抗噪性能。
[0004] 混沌作为非线性系统的一种独特现象,虽其动力学本质仍有待揭示,但其理论已经开始在很多工程问题中广泛应用。在信号处理领域,混沌理论已经在保密通信、混沌神经网络、混沌优化、图像处理、信号检测等方面取得了成功应用。在雷达领域,目前有混沌脉冲雷达、目标微多普勒信号检测、混沌杂波、雷达图像后处理等方面的应用。本发明则是将混沌理论应用于雷达成像环节中。
[0005] 对于发射线性调频(LFM)脉冲串信号的ISAR,在目标机动飞行且速度不是太高的条件下,本发明提出了一种基于混沌振子降噪的ISAR距离像增强算法。

发明内容

[0006] 本发明要解决的技术问题为:对于低信噪比条件下,ISAR的距离像淹没于噪声下,如何抑制噪声,增强ISAR距离像。现有文献并无该种方法,因为点目标的回波脉冲压缩后是一个冲激信号,对于这种一个采样点宽度的冲激信号,没有任何方法能增强。本发明利用特定方法绕过了这一难点,实现了ISAR距离像的增强。
[0007] 本发明提供的一种应用混沌振子的ISAR距离像增强方法,包括如下步骤:
[0008] 1)当逆合成孔径雷达ISAR发射宽带线性调频LFM脉冲串信号时,将ISAR接收到的目标回波信号混频,去除载频;
[0009] 2)去除载频后的混频信号在快时间维进行傅里叶变换,得到基带频域信号,设成像积累时间内接收了M个回波信号,每个回波信号采样点数为N,则形成M*N维的基带频域信号矩阵;
[0010] 3)根据目标的速度范围,选择基带频域信号矩阵的前L行,对该L行的数据进行keystone变换,则前L个回波包络走动被校正,对keystone变换后得到的L*N维数据矩阵,按行在快时间维逆傅里叶变换得到对齐的目标一维距离像;
[0011] 4)将步骤3)得到的对齐的L*N维距离像数据矩阵排列为一维向量,具体的排列方法为:将该矩阵的每一列取出,首尾相接,形成一个一维向量;
[0012] 5)将步骤4)得到的一维向量求模后,输入耦合Duffing振子降噪,耦合Duffing振子表示为如下二阶常微分方程组:
[0013]
[0014] 式中,x1是第一振子的状态变量,x2是第二振子的状态变量,ξ是振子的阻尼系数,k、q分别为振子间的耦合系数,fcos(t)是周期驱动力,f是周期驱动力的幅度,s(t)是待去噪的输入信号,振子状态变量的初始值可以任取,所述常微分方程组可由定步长四阶龙格‑库塔法求解;
[0015] 该耦合Duffing振子系统的降噪输出是振子状态变量的差值x1‑x2,所述降噪输出即为输入信号s(t)的去噪结果;
[0016] 6)将步骤5)得到的降噪后的一维向量重新恢复为L*N维的矩阵,其恢复方法是步骤4)的反向过程;
[0017] 7)将恢复后的L*N维矩阵的所有行相加,得到一个一维向量,并将该一维向量作为一个新的(M‑L+1)*N维观测图像的第一行;
[0018] 8)将基带频域信号矩阵的第2至第L+1行取出,重复步骤3)~7)的操作,得到(M‑L+1)*N维观测图像的第二行,以此类推,得到整个(M‑L+1)*N维观测图像,该观测图像就是增强后的ISAR距离像。
[0019] 作为优选方式,步骤3)中,L的选取不小于30,同时为确保目标在L个回波时间内的运动速度近似为恒定,L的取值小于100。步骤5)中,振子参数分别取:ξ=0.7、k=10、q=10、f=1.7。
[0020] 本发明的有益效果在于:本发明方法将叠加噪声的距离像信号经过处理后,输入耦合混沌振子系统中,达到抑制噪声、增强距离像信号的目的。对比常规ISAR的距离‑多普勒成像算法(R‑D成像算法)的脉冲压缩距离像,本发明方法能够实现ISAR距离像6dB的信噪比提升,具有推广应用的价值。

实施方案

[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] 以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,将该方法的修改、变化或延伸到其它应用场景,都应涵盖在本发明的包含范围之内。

附图说明

[0021] 图1为本发明的方法流程图。
[0022] 图2为本发明的方法验证时的仿真飞机模型。
[0023] 图3为仿真飞机模型在没有噪声时经过脉冲压缩后的距离像。
[0024] 图4为仿真飞机模型在有噪声时经过脉冲压缩后的距离像。
[0025] 图5为仿真飞机模型在有噪声时经过本发明方法增强后的距离像。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号