CN112924969B - 频率扫描阵列逆合成孔径雷达目标成像方法 - Google Patents
频率扫描阵列逆合成孔径雷达目标成像方法 Download PDFInfo
- Publication number
- CN112924969B CN112924969B CN202110137176.1A CN202110137176A CN112924969B CN 112924969 B CN112924969 B CN 112924969B CN 202110137176 A CN202110137176 A CN 202110137176A CN 112924969 B CN112924969 B CN 112924969B
- Authority
- CN
- China
- Prior art keywords
- target
- frequency
- radar
- scanning array
- array
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种频率扫描阵列逆合成孔径雷达目标成像方法,将频率分集的思想应用在使用频率扫描阵列的ISAR成像雷达上,通过发射多个单频信号并结合频率合成技术,实现用于频率扫描阵列的ISAR目标成像方法,既可以形成波束指向目标,又可以克服无法发射宽带信号的限制,实现对运动目标的二维成像。
Description
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种频率扫描阵列逆合成孔径雷达(Inverse synthetic aperture radar,ISAR)目标成像方法。
背景技术
频率扫描阵列雷达是一种电扫描成像雷达,通过改变工作频率而使天线单元间的相位发生变化,能够实现阵列天线波束扫描及指向目标。因为其天线波束变换灵活和硬件成本低等优势,频率扫描阵列雷达在目标识别方面呈现着重要的作用。将频率扫描阵列与逆合成孔径雷达系统相结合,可形成一种高增益的成像体制,然而频率扫描阵列由于各辐射单元间馈线产生的相位差不能匹配发射或接收信号频率的变化,因此该成像体制不能发射或接收宽带信号,这使得在使用频率扫描阵列时,不能进行需要发射宽带信号的ISAR成像。
针对上述问题,本发明将频率分集的思想应用在使用频率扫描阵列的ISAR成像雷达上,通过发射多个单频信号并结合频率合成技术,实现用于频率扫描阵列的ISAR目标成像方法,既可以形成波束指向目标,又可以克服无法发射宽带信号的限制,实现对运动目标的二维成像。
发明内容
本发明所要解决的技术问题是在使用频率扫描阵列进行ISAR成像时,如何克服无法发射宽带信号的限制,实现ISAR成像。为此,本发明提供一种基于频率扫描阵列的逆合成孔径雷达目标成像方法,该方法通过发射不同频率的单频信号,可同时实现波束指向目标及宽带信号合成,实现对运动目标的二维ISAR成像。
频率扫描阵列逆合成孔径雷达目标成像方法,包括以下步骤:
步骤1,初始化频率扫描阵列参数;
频率扫描阵列各辐射单元等间距排列,辐射单元个数为M,相邻两个辐射单元间距为d,各辐射单元由馈线L串联连接,天线扫描角度与发射信号的频率关系为:
式中θ为扫描角度,c为光速,f为频率扫描阵列发射信号频率,d为相邻两个辐射单元间距,l为相邻辐射单元的馈线长度,表示介质波长,εr表示相对介电常数;
步骤2,确定扫描频率fn;
在tn时刻,目标运动到某一位置,此时目标与法线的夹角为θn,利用几何关系可得:
式中v是目标的运动速度,α是飞行角度与水平方向的夹角,R0是目标与雷达的垂直距离;
假设目标做平稳匀速直线运动,即α=0°,令k=0,上两式联立可得目标运动位置角度与频率之间的关系:
步骤3,对目标进行观测;
雷达在不同观测时刻,向目标场景发射不同频率的单频信号,假设当前观测时间为tn,所需发射的频率值为fn,此时发射天线发射的信号可表示为:
x(tn)=s(tn)exp{j2πfntn} n=1,2,...,N
其中s(tn)为信号的复包络;
步骤4,获取回波信号;
选择阵列中第一个辐射单元为参考单元,对成像区域中某一目标散射点p进行观测,该散射点沿阵列射线的法向夹角为θp(tn),与雷达的距离为Rp(tn),则第m个辐射单元的回波信号时延τm为:
扫描信号经过目标反射,由接收装置接收并进行数字采样、混频后,得到tn时刻的回波信号:
其中σp表示该散射点的散射系数,m=1,2,...,M,M为辐射单元个数,由于发射信号是单频信号,s(tn)的变化可忽略不计;
步骤5,更新观测时间为tn+1,重复上述步骤(2)~(4),直至得到N个观测时刻的回波矩阵:
Xe(tn)=[xe(t1) xe(t2) … xe(tN)]
步骤6,使用后向投影成像算法,对成像区域进行网格划分,并对每个网格节点(x,y)进行对应的相位补偿,该点相对于第m个辐射单元的相位补偿项为:
其中Rx,y是网格节点(x,y)到雷达的距离,θx,y是该节点相对于雷达法线方向的角度,则tn时刻进行相位补偿后的回波矩阵可表示为:
步骤7,将N个观测时刻进行相位补偿的回波矩阵相参累加,获得目标二维像:
与现有技术相比,本发明的有益效果:
(1)建立了频率扫描阵列ISAR成像体制,使得频率的改变同时实现波束指向和宽带信号合成的功能,可实现对目标进行大转角ISAR二维成像。
(2)结合频率分集的思想,发射多个单频信号合成宽带信号,克服ISAR成像发射信号的大带宽限制,相比于相控阵ISAR成像,使用频率扫描阵列ISAR成像可降低阵列成本和复杂度。
附图说明
图1为本发明成像方法流程示意图。
图2为实施例频率扫描阵列结构模型图。
图3为实施例频率扫描阵列ISAR体制模型图。
图4为实施例成像场景目标散射点分布示意图。
图5为基于本发明成像方法的目标ISAR成像二维仿真结果图。
具体实施方式
下面结合实施例和附图对本发明内容进行进一步说明,但不是对本发明的限定:
实施例
参照图1,频率扫描阵列逆合成孔径雷达目标成像方法,包括如下步骤:
步骤1,初始化频率扫描阵列参数;
如图2所示,频率扫描阵列各辐射单元等间距排列,辐射单元个数为M,相邻两个辐射单元间距为d,各辐射单元由馈线L串联连接,天线扫描角度与发射信号的频率关系为:
式中θ为扫描角度,c为光速,f为发射频率,l为相邻辐射单元的馈线长度,表示介质波长,εr表示相对介电常数;
步骤2,确定扫描频率fn;
如图3所示,在tn时刻,目标运动到某一位置,此时目标与法线的夹角为θn,利用几何关系可得:
式中v是目标的运动速度,α是飞行角度与水平方向的夹角,R0是目标与雷达的垂直距离;
假设目标做平稳匀速直线运动,即α=0°,令k=0,上两式联立可得目标运动位置角度与频率之间的关系:
步骤3,对目标进行观测;
雷达在不同观测时刻,向目标场景发射不同频率的单频信号;假设当前观测时间为tn,所需发射的频率值为fn,此时发射天线发射的信号可表示为:
x(tn)=s(tn)exp{j2πfntn}n=1,2,...,N,
其中s(tn)为信号的复包络;
步骤4,获取回波矩阵;
选择阵列中第一个辐射单元为参考单元,对成像区域中某一目标散射点p进行观测,该散射点沿阵列射线的法向夹角为θp(tn),与雷达的距离为Rp(tn),则第m个辐射单元的回波信号时延τm为:
扫描信号经过目标反射,由接收装置接收并进行数字采样、混频后,得到tn时刻的回波信号:
其中σp表示该散射点的散射系数,m=1,2,...,M,M为辐射单元个数,由于发射信号是单频信号,s(tn)的变换可忽略不计;
步骤5,更新观测时间为tn+1,重复上述步骤(2)~(4),直至得到N个观测时刻的回波矩阵:
Xe(tn)=[xe(t1) xe(t2) … xe(tN)]
步骤6,使用后向投影成像算法,对成像区域进行网格划分,并对每个网格节点(x,y)进行对应的相位补偿,该点相对于第m个辐射单元的相位补偿项为:
其中Rx,y是网络节点(x,y)到雷达的距离,θx,y是该节点相对于雷达法线方向的角度,则tn时刻进行相位补偿后的回波矩阵可表示为:
步骤7,将N个观测时刻进行相位补偿的回波矩阵相参累加,获得目标二维像:
本发明有益效果可通过实验仿真进一步说明:
(1)实验场景设置
频率扫描阵列结构图如图2所示,其中辐射单元数M=15,阵元间距d=λ/2=0.1m,馈线l=0.1m,介电常数εr=2.3,雷达信号基准载频为f0=1.5GHz,合成带宽为1GHz。
频率扫描阵列ISAR成像雷达系统采用如图3所示的雷达收发体制,目标以3m/s的速度匀速向X轴正方向移动,距离雷达初始距离为R0=2×103m,目标相对于雷达转过的角度为30°,雷达共观测401次。
成像区域范围:X轴方向为0~100m,Y轴方向为0~100m,在区域等间距划分一个100×100的网格,相邻网络节点间隔为1m,网格点(50,50)为目标中心点。假设目标做刚体运动,将目标等效为5个散射点,分别设置在(50,50),(50,20),(50,80),(20,50),(80,50),目标在场景中的分布示意图如图4所示。
(2)仿真结果
图5为采用本发明所提出的频率扫描阵列ISAR成像体制雷达目标成像的仿真结果,与图4进行比较可以得出,使用频率扫描阵列ISAR体制对运动目标进行二维成像,仿真结果可以重建目标散射点的相对位置。仿真成像结果表明,本发明所提的频率扫描阵列ISAR目标成像方法,通过发射多个单频信号并结合频率合成技术,可以实现对运动目标的二维成像,解决了频率扫描阵列不能发射宽带信号的不足,相比于宽带相控阵ISAR成像,降低了系统复杂度与成本。
Claims (1)
1.频率扫描阵列逆合成孔径雷达目标成像方法,其特征在于,包括以下步骤:
步骤1,初始化频率扫描阵列参数;
频率扫描阵列各辐射单元等间距排列,辐射单元个数为M,相邻两个辐射单元间距为d,各辐射单元由馈线L串联连接,天线扫描角度与发射信号的频率关系为:
式中θ为扫描角度,c为光速,f为频率扫描阵列发射信号频率,d为相邻两个辐射单元间距,l为相邻辐射单元的馈线长度,表示介质波长,εr表示相对介电常数;
步骤2,确定扫描频率fn;
在tn时刻,目标运动到某一位置,此时目标与法线的夹角为θn,利用几何关系可得:
式中v是目标的运动速度,α是飞行角度与水平方向的夹角,R0是目标与雷达的垂直距离;
假设目标做平稳匀速直线运动,即α=0°,令k=0,上两式联立可得目标运动位置角度与频率之间的关系:
步骤3,对目标进行观测;
雷达在不同观测时刻,向目标场景发射不同频率的单频信号,假设当前观测时间为tn,所需发射的频率值为fn,此时发射天线发射的信号表示为:
x(tn)=s(tn)exp{j2πfntn} n=1,2,...,N
其中s(tn)为信号的复包络;
步骤4,获取回波信号;
选择阵列中第一个辐射单元为参考单元,对成像区域中某一目标散射点p进行观测,该散射点沿阵列射线的法向夹角为θp(tn),与雷达的距离为Rp(tn),则第m个辐射单元的回波信号时延τm为:
扫描信号经过目标反射,由接收装置接收并进行数字采样、混频后,得到tn时刻的回波信号:
其中σp表示该散射点的散射系数,m=1,2,...,M,M为辐射单元个数,由于发射信号是单频信号,s(tn)的变化忽略不计;
步骤5,更新观测时间为tn+1,重复上述步骤(2)~(4),直至得到N个观测时刻的回波矩阵:
Xe(tn)=[xe(t1) xe(t2)…xe(tN)]
步骤6,使用后向投影成像算法,对成像区域进行网格划分,并对每个网格节点(x,y)进行对应的相位补偿,该点相对于第m个辐射单元的相位补偿项为:
其中Rx,y是网格节点(x,y)到雷达的距离,θx,y是该节点相对于雷达法线方向的角度,则tn时刻进行相位补偿后的回波矩阵表示为:
步骤7,将N个观测时刻进行相位补偿的回波矩阵相参累加,获得目标二维像:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110137176.1A CN112924969B (zh) | 2021-02-01 | 2021-02-01 | 频率扫描阵列逆合成孔径雷达目标成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110137176.1A CN112924969B (zh) | 2021-02-01 | 2021-02-01 | 频率扫描阵列逆合成孔径雷达目标成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112924969A CN112924969A (zh) | 2021-06-08 |
CN112924969B true CN112924969B (zh) | 2023-07-21 |
Family
ID=76169339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110137176.1A Active CN112924969B (zh) | 2021-02-01 | 2021-02-01 | 频率扫描阵列逆合成孔径雷达目标成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112924969B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017198162A1 (zh) * | 2016-04-29 | 2017-11-23 | 深圳市太赫兹科技创新研究院有限公司 | 基于合成孔径雷达成像的三维图像重建方法及装置 |
CN110007303A (zh) * | 2019-04-22 | 2019-07-12 | 桂林电子科技大学 | 频率分集阵列合成孔径三维成像雷达系统及其成像方法 |
CN110261852A (zh) * | 2019-03-22 | 2019-09-20 | 桂林电子科技大学 | 一种频率分集逆合成孔径雷达成像系统及方法 |
-
2021
- 2021-02-01 CN CN202110137176.1A patent/CN112924969B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017198162A1 (zh) * | 2016-04-29 | 2017-11-23 | 深圳市太赫兹科技创新研究院有限公司 | 基于合成孔径雷达成像的三维图像重建方法及装置 |
CN110261852A (zh) * | 2019-03-22 | 2019-09-20 | 桂林电子科技大学 | 一种频率分集逆合成孔径雷达成像系统及方法 |
CN110007303A (zh) * | 2019-04-22 | 2019-07-12 | 桂林电子科技大学 | 频率分集阵列合成孔径三维成像雷达系统及其成像方法 |
Non-Patent Citations (3)
Title |
---|
An Anti-jamming Method Against Frequency Diverse Array for ISAR by Spatial Location Feature Recognition;Zhidong Liu等;《2020 IEEE 11th Sensor Array and Multichannel Signal Processing Workshop (SAM)》;全文 * |
Two-Dimensional Imaging Using Frequency Diversity with Inverse Synthetic Aperture;Yi Du等;《2019 6th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR)》;全文 * |
分布式MIMO-ISAR子图像融合方法;谷文堃等;《雷达学报》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112924969A (zh) | 2021-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Charvat et al. | An ultrawideband (UWB) switched-antenna-array radar imaging system | |
CN108983234B (zh) | 太赫兹isar成像系统及基于后向投影的图像重建方法 | |
JP5080795B2 (ja) | イメージング装置 | |
Wang | Large-area remote sensing in high-altitude high-speed platform using MIMO SAR | |
CN112526512B (zh) | 大功率大口径宽带毫米波空馈相控阵雷达系统及成像方法 | |
Klare et al. | First experimental results with the imaging MIMO radar MIRA-CLE X | |
CN109597076B (zh) | 用于地基合成孔径雷达的数据处理方法及装置 | |
WO2023015067A2 (en) | Millimeter-wave massive mmo fmcw radar with binary-phase-coded ofdm | |
CN113835087B (zh) | 一种脉间跳频积累的随机辐射雷达高分辨成像方法 | |
CN112924969B (zh) | 频率扫描阵列逆合成孔径雷达目标成像方法 | |
CN111146560B (zh) | 一种复合馈源抛物柱面天线及探测卫星 | |
Cetinkaya et al. | Millimeter-wave MIMO array based on semi-circular topology | |
Lim et al. | Azimuth beam pattern synthesis for airborne SAR system optimization | |
CN113064165A (zh) | 一种扫描雷达俯仰-方位二维超分辨方法 | |
JP3649565B2 (ja) | 合成開口レーダ装置 | |
Schuetz et al. | A promising outlook for imaging radar: Imaging flash radar realized using photonic spatial beam processing | |
CN110471059A (zh) | 一种小型化毫米波成像雷达装置 | |
JP2019120659A (ja) | レーダ装置及びレーダ方式 | |
US20210239789A1 (en) | Conical scan weather radar | |
RU2682169C1 (ru) | Способ повышения разрешающей способности рлс по азимуту и дальности и уменьшения времени сканирования наземных объектов при посадке самолета и приемное устройство, реализующее этот способ | |
Wang et al. | A new airborne Ka-band DBF-SAR system and flight test results | |
CN114740468B (zh) | 一种宽波束凝视体制稀疏阵同时多目标瞬时成像雷达系统 | |
Loomis | Digital beamforming-a retrospective | |
CN114660595B (zh) | 一种发射涡旋电磁波的频率分集圆阵雷达成像方法 | |
Knott et al. | RADAR systems-technology and challenges |
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 |