CN110261852B - 一种频率分集逆合成孔径雷达成像系统及方法 - Google Patents

一种频率分集逆合成孔径雷达成像系统及方法 Download PDF

Info

Publication number
CN110261852B
CN110261852B CN201910221288.8A CN201910221288A CN110261852B CN 110261852 B CN110261852 B CN 110261852B CN 201910221288 A CN201910221288 A CN 201910221288A CN 110261852 B CN110261852 B CN 110261852B
Authority
CN
China
Prior art keywords
frequency
target
imaging
radar
baseband signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910221288.8A
Other languages
English (en)
Other versions
CN110261852A (zh
Inventor
廖可非
杜毅
欧阳缮
李晶晶
刘扬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN201910221288.8A priority Critical patent/CN110261852B/zh
Publication of CN110261852A publication Critical patent/CN110261852A/zh
Application granted granted Critical
Publication of CN110261852B publication Critical patent/CN110261852B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种频率分集逆合成孔径雷达成像系统和方法,包括发射部分和接收部分,发射部分包括频率选择单元、基带信号产生器、上变频器和发射天线,频率选择单元的输出分为两路,一路连接基带信号生成器,一路送至成像算法单元,频率选择单元在每个观测时刻从频率偏置集合中按照线性或随机等规律选择单个频率偏置作为基带信号产生器的频率,经过上变频器进行载波调制,得到射频发射信号由发射天线发射出;接收部分包括接收天线、下变频器、数字采样单元和成像算法单元,通过收到的回波信号经过下变频器和数字采样单元,恢复为携带目标信息的基带信号,送到成像算法单元进行处理。本发明在每个观测时间点仅发射一个单频信号,降低了系统和成像方法的复杂度。

Description

一种频率分集逆合成孔径雷达成像系统及方法
技术领域
本发明涉及雷达信号处理技术领域,尤其涉及一种频率分集逆合成孔径雷达成像系统及方法。
背景技术
随着宽带雷达技术的不断发展和应用,雷达成像技术在国防和民用领域取得了重大的应用和发展,它不仅大大提高了雷达获取信息的能力,而且为非合作目标识别提供了一种重要方式,是雷达发展史上的一个里程碑。ISAR成像技术可以全天候的对空中目标进行远距离的二维成像,是现代战场战略防御系统中一种极其重要的目标识别手段。在获得目标ISAR像之后,可以根据目标的大致尺寸与形状能够初步确定其类别,不需要样本库训练,这对提高部队快速反应能力,提高部队的指挥自动化水平、攻防能力、国土防空反导、反舰以及战略警能力等都具有十分重要的作用。现有的ISAR成像系统需要发射宽带信号以获得有较高分辨率的目标距离像,这对发射/接收机有较高的要求,也大大的增加了系统的复杂度和造价。
但是,现有的宽带ISAR成像方法需要经过脉冲压缩、插值、重采样过程,在成像算法上较复杂。
发明内容
针对现有技术的不足,本发明所解决的技术问题是如何降低ISAR成像系获得有较高分辨率的目标距离像的复杂度。
为解决上述技术问题,本发明采用的技术方案是一种频率分集逆合成孔径雷达成像系统,包括发射部分和接收部分,所述发射部分包括频率选择单元、基带信号产生器、上变频器和发射天线,所述频率选择单元的输出分为两路,一路连接基带信号生成器,一路送至成像算法单元,所述频率选择单元在每个观测时刻从频率偏置集合中按照线性或随机等规律选择单个频率偏置作为基带信号产生器的频率,经过上变频器进行载波调制,得到射频发射信号由发射天线发射出去;所述接收部分包括接收天线、下变频器、数字采样单元和成像算法单元,通过收到的回波信号经过下变频器和数字采样单元,恢复为携带目标信息的基带信号,送到成像算法单元进行处理。
本发明解决的另一个问题是频率分集逆合成孔径雷达成像方法。
为此,本发明提供的技术方案是一种频率分集逆合成孔径雷达成像方法,包括以下步骤:
(1)雷达以一定的时间间隔观测目标,其观测时间为ti i=1,2,...,N;
(2)在初始观测时刻ti=1,频率选择单元从频率偏置集合{0Δf,1Δf,…(L-1)Δf}按照线性或随机等规则选择单个频率:fi作为单频基带信号,其中Δf为固定频率增量;
(3)基带信号产生器产生频率为fi的单频信号Bi(t),Bi(t)=cos(2πfit),t为时间;
(4)基带信号Bi(t)通过上变频器,得到射频发射信号Si(t),由发射天线发射出去,发射信号可表示为:
Si(t)=cos(2π(f0+fi)t)
(5)经目标反射后,假设目标由M个散射点组成,经过接收天线接收、下变频和数字采样,目标回波信号可表示为:
Figure BDA0002003686520000031
其中m表示第m个散射点,σm表示第m个散射点的散射系数,其中Rm(ti)表示ti时刻时,第m个散射点与雷达的距离,c表示光速;
(6)将成像区域按X轴和Y轴划分网格(像素点),对成像区域中每一个网格节点(x,y)进行相应的时延补偿,对其中某一个网格点的时延补偿可表示为:
Figure BDA0002003686520000032
其中R(x,y)(ti)表示ti时刻时,网格点(x,y)与雷达的距离,得到成像场景中所有网格点上时延补偿后的回波矩阵
Figure BDA0002003686520000034
(7)令ti=ti+1,重复上述步骤(2)-(6),直到i=N,得到所有回波矩阵
Figure BDA0002003686520000035
(8)将N个时刻获得的回波矩阵
Figure BDA0002003686520000036
相加,获得目标二维像:
Figure BDA0002003686520000033
与现有技术相比,本发明的有益效果:
(1)在使用单个雷达收发单元的情况下,可用目标与雷达相对运动产生逆合成孔径,并使用多个单频信号完成频率分集宽带信号的合成,获得目标二维成像,降低了系统复杂度,在匀速直线目标成像场景中尤为有效。
(2)采用一个雷达收发单元分时工作,每个观测时间点只需发射一个单频信号,相比于传统的宽带信号,降低了对系统设备的要求。
(3)目标成像方法的发射信号形式简单,无需脉冲压缩,即可得到距离向分辨率,且一次补偿相位即可同时得到目标的距离向、方位向图像。相比于宽带信号ISAR后向投影算法,不需要脉冲压缩、插值、重采样,降低了在成像算法上的复杂度。
本发明相比于现有的宽带ISAR成像系统和成像方法,在每个观测时间点仅发射一个单频信号,降低了系统和成像方法的复杂度。
附图说明
图1为雷达系统与目标几何结构关系图;
图2为本发明系统结构示意图;
图3为成像模型图;
图4为本发明方法流程示意图;
图5为多目标仿真场景图;
图6为多目标仿真成像结果图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明,但不是对本发明的限定。
雷达系统与目标几何结构关系,如图1所示,由1个雷达收发单元(空心正方形)和目标所构成,雷达在不同时刻ti向目标发射单频信号。
图2示出了一种频率分集逆合成孔径雷达成像系统,包括发射部分和接收部分,所述发射部分包括频率选择单元、基带信号产生器、上变频器和发射天线,所述频率选择单元的输出分为两路,一路连接基带信号生成器,一路送至成像算法单元,所述频率选择单元在每个观测时刻从频率偏置集合中按照线性或随机等规律选择单个频率偏置作为基带信号产生器的频率,经过上变频器进行载波调制,得到射频发射信号由发射天线发射出去;所述接收部分包括接收天线、下变频器、数字采样单元和成像算法单元,通过收到的回波信号经过下变频器和数字采样单元,恢复为携带目标信息的基带信号,送到成像算法单元进行处理,得到的成像模型如图3所示。
图4示出了本发明方法流程,一种频率分集逆合成孔径雷达成像方法,包括以下步骤:
(1)雷达以一定的时间间隔观测目标,其观测时间为ti i=1,2,...,N;
(2)在初始观测时刻ti=1,频率选择单元从频率偏置集合{0Δf,1Δf,…(L-1)Δf}按照线性或随机等规则选择单个频率:fi作为单频基带信号,其中Δf为固定频率增量;
(3)基带信号产生器产生频率为fi的单频信号Bi(t),Bi(t)=cos(2πfit),t为时间;
(4)基带信号Bi(t)通过上变频器,得到射频发射信号Si(t),由发射天线发射出去,发射信号可表示为:
Si(t)=cos(2π(f0+fi)t)
(5)经目标反射后,假设目标由M个散射点组成,经过接收天线接收、下变频和数字采样,目标回波信号可表示为:
Figure BDA0002003686520000051
其中m表示第m个散射点,σm表示第m个散射点的散射系数,其中Rm(ti)表示ti时刻时,第m个散射点与雷达的距离,c表示光速;
(6)将成像区域按X轴和Y轴划分网格(像素点),对成像区域中每一个网格节点(x,y)进行相应的时延补偿,对其中某一个网格点的时延补偿可表示为:
Figure BDA0002003686520000061
其中R(x,y)(ti)表示ti时刻时,网格点(x,y)与雷达的距离,得到成像场景中所有网格点上时延补偿后的回波矩阵
Figure BDA0002003686520000063
(7)令ti=ti+1,重复上述步骤(2)-(6),直到i=N,得到所有回波矩阵
Figure BDA0002003686520000064
(8)将N个时刻获得的回波矩阵
Figure BDA0002003686520000065
相加,获得目标二维像:
Figure BDA0002003686520000062
实验场景:
频率分集逆合成孔径雷达成像系统采用如图1所示的目标与雷达结构图,目标以100m/s的速度向Y轴的正方向运动,雷达一共观测241次数。雷达信号载频f0=10GHz,频率增量Δf=200KHz,频率增量个数N=241。
成像区域范围:X轴方向为-10m~10m,Y轴方向为-10m~10m,在区域等间距划分一个100×100的网格,相邻网络节点间隔为0.202m,网格中(50,50)的位置,即为目标中心(即坐标轴原点)。现将5个点目标分别设置在(0,0),(0,-4.04),(0,4.04),(-4.04,0),(4.04,0),假设目标散射系数为1,目标如图5所示。
实验结果:
图6为采用本文方法的仿真雷达探测与成像结果,比较图5和图6可以看出,仿真结果不仅重建了各目标的相对位置,也准确的重建了目标的相对散射系数。
图6中的距离向表示雷达与目标的连线方向、方位向为与雷达目标连线垂直的方向,信号强度为回波矩阵多次相加,各网格节点上的数值大小。
与现有技术相比,本发明有益效果通过实验结果进行进一步说明。
实验场景:
频率分集逆合成孔径雷达成像系统采用如图1所示的目标与雷达结构图,目标以100m/s的速度向Y轴的正方向运动,雷达一共观测241次数。雷达信号载频f0=10GHz,频率增量Δf200KHz,频率增量个数N=241。
成像区域范围:X轴方向为-10m~10m,Y轴方向为-10m~10m,在区域等间距划分一个100×100的网格,相邻网络节点间隔为0.202m,网格中(50,50)的位置,即为目标中心(即坐标轴原点)。现将5个点目标分别设置在(0,0),(0,-4.04),(0,4.04),(-4.04,0),(4.04,0),假设目标散射系数为1,目标如图5所示。
实验结果:
图6为采用本文方法的仿真雷达探测与成像结果,比较图5和图6可以看出,仿真结果不仅重建了各目标的相对位置,也准确的重建了目标的相对散射系数。从成像实验结果可以看出,本文所提的频率分集逆合成孔径雷达系统在所提的成像算法下,获得了较好的目标二维像,相比于宽带雷达成像系统,降低了系统和成像方法的复杂度。
以上结合附图对本发明的实施方式做出了详细说明,但本发明不局限于所描述的实施方式。对于本领域技术人员而言,在不脱离本发明的原理和精神的情况下,对这些实施方式进行各种变化、修改、替换和变型仍落入本发明的保护范围内。

Claims (1)

1.一种频率分集逆合成孔径雷达成像方法,其特征在于,包括以下步骤:
(1)雷达以一定的时间间隔观测目标,其观测时间为tii=1,2,...,N;
(2)在初始观测时刻ti=1,频率选择单元从频率偏置集合{0Δf,1Δf,…(L-1)Δf}按照线性或随机规则选择单个频率:fi作为单频基带信号,其中Δf为固定频率增量;
(3)基带信号产生器产生频率为fi的单频信号Bi(t),Bi(t)=cos(2πfit),t为时间;
(4)基带信号Bi(t)通过上变频器,得到射频发射信号Si(t),由发射天线发射出去,发射信号可表示为:
Si(t)=cos(2π(f0+fi)t)
(5)经目标反射后,假设目标由M个散射点组成,经过接收天线接收、下变频和数字采样,目标回波信号可表示为:
Figure FDA0003613086060000011
其中m表示第m个散射点,σm表示第m个散射点的散射系数,其中Rm(ti)表示ti时刻时,第m个散射点与雷达的距离,c表示光速;
(6)将成像区域按X轴和Y轴划分网格像素点,对成像区域中每一个网格像素点(x,y)进行相应的时延补偿,对其中某一个网格像素点的时延补偿可表示为:
Figure FDA0003613086060000021
其中R(x,y)(ti)表示ti时刻时,网格点(x,y)与雷达的距离,得到成像场景中所有网格像素点上时延补偿后的回波矩阵
Figure FDA0003613086060000022
(7)令ti=ti+1,重复上述步骤(2)-(6),直到i=N,得到所有回波矩阵
Figure FDA0003613086060000023
(8)将N个时刻获得的回波矩阵
Figure FDA0003613086060000024
相加,获得目标二维像:
Figure FDA0003613086060000025
CN201910221288.8A 2019-03-22 2019-03-22 一种频率分集逆合成孔径雷达成像系统及方法 Active CN110261852B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910221288.8A CN110261852B (zh) 2019-03-22 2019-03-22 一种频率分集逆合成孔径雷达成像系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910221288.8A CN110261852B (zh) 2019-03-22 2019-03-22 一种频率分集逆合成孔径雷达成像系统及方法

Publications (2)

Publication Number Publication Date
CN110261852A CN110261852A (zh) 2019-09-20
CN110261852B true CN110261852B (zh) 2022-06-17

Family

ID=67913484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910221288.8A Active CN110261852B (zh) 2019-03-22 2019-03-22 一种频率分集逆合成孔径雷达成像系统及方法

Country Status (1)

Country Link
CN (1) CN110261852B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111273269B (zh) * 2020-02-18 2022-06-03 桂林电子科技大学 基于ipso-bp的频率分集阵列的雷达目标定位方法
CN112924969B (zh) * 2021-02-01 2023-07-21 桂林电子科技大学 频率扫描阵列逆合成孔径雷达目标成像方法
CN113093136B (zh) * 2021-03-31 2022-06-10 桂林电子科技大学 一种频率分集阵列雷达去目标位置模糊成像方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100785473B1 (ko) * 2006-10-18 2007-12-13 삼성전자주식회사 멀티-밴드 ofdm 스킴에서 수신 신호의 주파수 옵셋을보상하는 수신기 및 방법
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001074829A (ja) * 1999-09-06 2001-03-23 Mitsubishi Electric Corp Fmパルスドップラレーダ用送受信装置
CN101846741B (zh) * 2010-05-07 2012-11-21 北京航空航天大学 一种逆合成孔径雷达成像数据段选择方法
JP6304257B2 (ja) * 2013-09-12 2018-04-04 パナソニック株式会社 レーダ装置、車両及び移動体速度検出方法
US10502824B2 (en) * 2015-11-09 2019-12-10 Infineon Technologies Ag Frequency modulation scheme for FMCW radar
CN105589070B (zh) * 2015-12-11 2017-12-05 桂林电子科技大学 基于频率分集阵列的雷达目标成像方法
CN106646394B (zh) * 2017-02-28 2018-03-02 湖北工业大学 一种fda‑mimo雷达中基于循环优化的发射波束矩阵设计方法
CN107132536B (zh) * 2017-04-10 2019-10-11 中国科学院国家空间科学中心 一种消除目标微动对雷达成像干扰的方法
CN107356921B (zh) * 2017-08-11 2020-09-11 桂林电子科技大学 一种频率分集阵列雷达基于一次频偏目标定位的方法
CN108732545A (zh) * 2018-03-06 2018-11-02 电子信息系统复杂电磁环境效应国家重点实验室 基于多通道混频匹配滤波的频控阵信号接收方法及结构
CN108717179A (zh) * 2018-03-06 2018-10-30 电子信息系统复杂电磁环境效应国家重点实验室 基于多波形自适应处理的频率分集雷达及方法
CN108776337B (zh) * 2018-04-24 2021-11-05 桂林电子科技大学 Mimo-fda探地雷达近目标二维成像方法
CN108983234B (zh) * 2018-06-28 2022-03-18 山西大学 太赫兹isar成像系统及基于后向投影的图像重建方法
CN109031296B (zh) * 2018-07-26 2020-07-07 中国人民解放军国防科技大学 基于随机间歇采样线性调频信号的宽带雷达成像方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100785473B1 (ko) * 2006-10-18 2007-12-13 삼성전자주식회사 멀티-밴드 ofdm 스킴에서 수신 신호의 주파수 옵셋을보상하는 수신기 및 방법
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法

Also Published As

Publication number Publication date
CN110261852A (zh) 2019-09-20

Similar Documents

Publication Publication Date Title
CN111537966B (zh) 一种适用于毫米波车载雷达领域的阵列天线误差校正方法
CN110261852B (zh) 一种频率分集逆合成孔径雷达成像系统及方法
WO2021062914A1 (zh) 一种基于mimo体制的提高雷达角度分辨率的布局及方法
CN105445701B (zh) Ddma‑mimo雷达目标的单脉冲角度估计方法
CN110007303B (zh) 频率分集阵列合成孔径三维成像雷达系统及其成像方法
CN108226871B (zh) 基于毫米波避撞雷达的分集相控阵设计方法
CN113314832B (zh) 一种毫米波车载mimo雷达天线阵列装置及设计方法
CN107064880A (zh) 分布式多基雷达收发波束同步和波束控制的高精度方法
CN105301589B (zh) 高分辨宽测绘带sar地面运动目标成像方法
CN109597076B (zh) 用于地基合成孔径雷达的数据处理方法及装置
CN113030946B (zh) 二次雷达探测方法、装置、设备、系统、介质及程序产品
CN108562894B (zh) 雷达波束指向与发射功率的分配方法
CN109655819B (zh) 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法
CN115128608B (zh) 车载多普勒分集mimo合成孔径雷达时域成像方法
JP6716352B2 (ja) レーダシステム及びそのレーダ信号処理方法
CN114185038B (zh) 一种基于毫米波雷达网实现超高角度分辨力的方法及系统
CN106019280A (zh) 基于距离多普勒校正的fmcw sar成像方法及装置
CN109633587B (zh) 一种组网雷达信号带宽的自适应调整方法
CN107783111B (zh) 一种基于最大熵准则的雷达前视超分辨成像方法
CN114185047B (zh) 一种基于最优极坐标变换的双基sar动目标重聚焦方法
CN116148891A (zh) 基于顺轨干涉的星地双基海面二维流场测量方法
CN112946650B (zh) 一站固定双站低频超宽带sar运动目标检测与成像方法
CN114994676A (zh) 一种一站固定式双站低频超宽带sar运动目标成像方法
CN113960597A (zh) 一种多目标isar成像时间点自动选取方法
CN114966684A (zh) 群多基合成孔径雷达快速分解投影成像方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant