CN111257878A - 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法 - Google Patents

基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法 Download PDF

Info

Publication number
CN111257878A
CN111257878A CN202010222643.6A CN202010222643A CN111257878A CN 111257878 A CN111257878 A CN 111257878A CN 202010222643 A CN202010222643 A CN 202010222643A CN 111257878 A CN111257878 A CN 111257878A
Authority
CN
China
Prior art keywords
pulse
signal
frequency
time
phi
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.)
Granted
Application number
CN202010222643.6A
Other languages
English (en)
Other versions
CN111257878B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202010222643.6A priority Critical patent/CN111257878B/zh
Publication of CN111257878A publication Critical patent/CN111257878A/zh
Application granted granted Critical
Publication of CN111257878B publication Critical patent/CN111257878B/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
    • G01S13/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • G01S13/9011SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
    • 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
    • G01S13/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9094Theoretical aspects
    • 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
    • G01S7/414Discriminating targets with respect to background clutter
    • 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
    • G01S7/418Theoretical aspects

Landscapes

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

Abstract

本发明提出了一种基于俯仰维频率脉内扫描的高分宽幅SAR的波形设计方法。旨在通过俯仰维频率脉内扫描技术来完成单通道高分宽幅SAR的波形设计。实施步骤为:设置输入参数;构建发射天线阵列;合成发射波形;提取发射波形中的俯仰维频率脉内扫描信号的时‑空方向图;计算时‑空方向图的瞬时波束指向角;计算发射天线阵列的发射天线数量;计算发射信号的信号带宽;计算发射信号的中心频率和相邻发射天线之间的时间延迟;获取俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果。本发明通过设计发射波形,来控制合成方向图的瞬时波束指向角实现脉内波束连续扫描,利用频域滤波等效实现空域滤波完成距离模糊抑制,工程实现较为简单。

Description

基于俯仰维频率脉内扫描高分宽幅SAR的波形设计方法
技术领域
本发明属于雷达信号处理技术领域,涉及一种高分宽幅SAR的波形设计方法,具体涉及一种基于俯仰维频率脉内扫描的高分宽幅SAR的波形设计方法。
背景技术
合成孔径雷达(SAR)由于具有全天时,全天候可对地面场景进行成像的特点,在军事领域和民事领域都有着广泛的应用。在一些应用场景下,如全球动态观测,地物目标分类和识别等,要求SAR具有高分宽幅的成像能力。但传统的SAR成像中,高分辨率和宽测绘带由于对脉冲重复频率PRF有相反的要求,是一对矛盾。高分辨率要求合成孔径时间长,所形成的多普勒带宽足够大,为了避免方位模糊,PRF与多普勒带宽需要满足奈奎斯特采样定理,即高的方位分辨率需要高PRF。宽测绘带成像需要满足测绘带内的所有回波可以在一个脉冲重复周期内返回来避免距离模糊。所以宽测绘带成像需要低的PRF。PRF的选择难以同时满足距离无模糊和方位无模糊的要求。传统成像模式只能在一方面满足要求,如聚束模式可以实现高分成像,但牺牲了成像幅宽,扫描模式可以实现宽幅成像,但牺牲了方位分辨率。
传统的成像模式都是单通道模式,只能在高分和宽幅之间做一个折中的选择来满足不同的应用场景。为了同时兼顾方位高分辨和距离幅宽,需要在多普勒域或距离域进行模糊抑制。国内外提出了多种基于多通道体制和数字波束形成DBF相结合的高分宽幅成像系统。利用DBF进行多普勒模糊抑制或距离模糊抑制。基于DBF的模糊抑制技术取决于独立空间的样本数量,因此接收通道数越多,模糊抑制效果越好。随着对成像性能要求的提高,多通道SAR所需的通道数越来越多,设备量大。针对这个问题,Federica Bordoni等人于2019年在IEEE Geoscience and Remote Sensing Letters上发表了名称为“Multifrequency Subpulse SAR:Exploiting Chirp Bandwidth for an IncreasedCoverage”的文章,公开了一种多频子脉冲技术用于高分宽幅成像,该方法利用多子脉冲占据不同的频段来抑制距离模糊,可以不依赖于DBF技术实现高分宽幅成像。但上述方法对波束的控制是从子脉冲到子脉冲的离散波束控制,需要在子脉冲之间做快速的相位变换,在工程实现时存在困难和限制。
发明内容
本发明的目的在于针对上述已有技术的不足,提出基于俯仰维频率脉内扫描的高分宽幅SAR信号的波形设计方法,旨在实现高分宽幅成像时,可以不依赖于DBF技术实现距离模糊抑制,并且是脉内波束连续扫描,工程实现较为简单。
为实现上述目的,本发明采取的技术方案包括如下步骤:
(1)设置输入参数:
根据SAR成像的指标要求,设置SAR成像的测绘带大小为W、距离分辨率为ρr、方位分辨率为ρa,测绘带近端入射角
Figure BDA0002426622140000021
SAR平台的高度和速度分别为H和V;
(2)构建发射天线阵列Z:
构建包括周期性排布的N个发射天线的发射天线阵列Z=[z1,z2,…,zn,…,zN],每个发射天线zn上连接有时间延迟线TTDn,相邻发射天线之间的距离和时间延迟分别为d和τ,其中,N≥2;
(3)合成发射波形s(t,φ):
(3a)设发射天线z1为参考发射天线,并获取第n个发射天线zn的延迟时间(n-1)·τ,以及发射机产生的发射信号g(t):
Figure BDA0002426622140000022
其中,t表示快时间,K表示发射信号g(t)的调频率,K=B/TP,B和TP分别表示发射信号g(t)的带宽和脉宽,fc表示发射信号g(t)的中心频率;
(3b)通过(n-1)·τ和g(t),计算发射机产生的发射信号g(t)经过TTDn时间延迟后的时间延迟信号sn(t),得到时间延迟信号集合S:
S=[s1(t),s2(t),…,sn(t),…,sN(t)]
sn(t)=g(t-(n-1)·τ);
(3c)对S中的所有时间延迟信号在远场进行合成,得到目标的俯仰角为φ的发射波形s(t,φ);
(4)提取发射波形s(t,φ)中的俯仰维频率脉内扫描信号的时-空方向图pe(t,φ):
计算发射波形s(t,φ)的模sabs(t,φ),并将sabs(t,φ)和发射信号g(t)的商作为时-空方向图pe(t,φ);
(5)计算时-空方向图pe(t,φ)的瞬时波束指向角φpeak(t):
令时-空方向图pe(t,φ)的分母为0,求解瞬时波束指向角φpeak(t),即
Figure BDA0002426622140000031
其中,
Figure BDA0002426622140000032
round(·)表示取整操作,0≤t≤TP
(6)计算发射天线阵列Z的发射天线数量N:
(6a)根据SAR成像的测绘带W和测绘带近端的入射角
Figure BDA0002426622140000033
计算测绘带远端的入射角
Figure BDA0002426622140000034
Figure BDA0002426622140000035
(6b)根据SAR平台速度V和SAR成像方位分辨率ρa,计算发射信号g(t)的脉冲重复频率Fr
Figure BDA0002426622140000036
(6c)根据测绘带远端的入射角
Figure BDA0002426622140000037
和发射信号g(t)的脉冲重复频率Fr,计算俯仰维波束宽度θ3dB
Figure BDA0002426622140000038
(6d)根据俯仰维波束宽度θ3dB,计算发射天线阵列Z的发射天线数量N:
Figure BDA0002426622140000041
其中,ceil(·)表示向上取整;
(7)计算发射信号g(t)的信号带宽B:
(7a)根据距离分辨率ρr,确定测绘带的远端俯仰角φfar分配的信号带宽B(φfar):
B(φfar)=2c/ρr
(7b)根据测绘带的近端入射角
Figure BDA0002426622140000042
和远端入射角
Figure BDA0002426622140000043
确定俯仰维频率脉内扫描的波束扫描范围θscan
Figure BDA0002426622140000044
(7c)根据波束扫描范围θscan,计算发射信号g(t)的信号带宽B:
Figure BDA0002426622140000045
(8)计算发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ:
(8a)设定波束扫描范围θscan的起始角度和终止角度分别为θ1和θ2,求解时-空方向图pe(t,φ)的瞬时波束指向角分别为θ1和θ2时的m1和m2
Figure BDA0002426622140000046
(8b)根据给定的期望的发射信号g(t)的中心频率为fc0,并通过假定m1=m2,求解期望的时间延迟τ0
Figure BDA0002426622140000047
(8c)根据期望的时间延迟τ0,得出m1和m2
Figure BDA0002426622140000051
Figure BDA0002426622140000052
将m1和m2带入方程组
Figure BDA0002426622140000053
求解出发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ;
(9)获取俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果:
根据计算得到的发射信号g(t)的信号带宽B和中心频率fc,发射天线阵列Z的发射天线数量N和相邻发射天线之间的时间延迟τ,带入俯仰维频率脉内扫描信号的时-空方向图pe(t,φ),得到俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果。
本发明与现有技术相比,具有如下优点:
本发明通过构造俯仰维频率脉内扫描阵列,在每个发射天线之后连接一个TTD,通过设计俯仰维频率脉内扫描阵列的阵元数,阵元间延迟,发射信号的带宽,载频和脉宽参数来控制俯仰维频率扫描阵列的合成方向图的瞬时波束指向角,实现脉内波束连续扫描。测绘带内距离模糊分量被隔离到不同的距离频段,可通过频域滤波等效实现空域滤波。相比于多通道SAR实现高分宽幅成像,所需的接收通道较少,设备量小,不存在通道间的幅相误差;相比于多频子脉冲实现单通道波束扫描,本发明是通过波形设计来实现脉内波束连续扫描,不需要做子脉冲分割处理和子脉冲间的波束转向控制,工程实现简单。
附图说明
图1是本发明的实现流程图;
图2是本发明的俯仰维频率脉内扫阵列的几何模型示意图;
图3是本发明的时-空方向图的仿真结果图;
图4是本发明的测绘带内的距离分辨率的仿真结果图;
图5是本发明测绘带内的点目标及其距离模糊分量分布的仿真结果图;
图6是本发明测绘带内的距离模糊比的仿真结果图。
具体实施方式
下面结合附图和具体实施例,对本发明作进一步的描述。
参照图1,本发明的具体实现步骤如下:
步骤1)设置输入参数:
根据SAR成像的指标要求,设置SAR成像的测绘带大小为W、距离分辨率为ρr、方位分辨率为ρa,测绘带近端入射角
Figure BDA0002426622140000061
SAR平台的高度和速度分别为H和V。本实施例中设定成像的测绘带大小W=100km,距离分辨率ρr=0.75m,方位分辨率ρa=0.75m,测绘带近端入射角
Figure BDA0002426622140000062
SAR平台的高度和速度分别为H=600km和V=7500m/s。
步骤2)构建发射天线阵列Z:
构建包括周期性排布的N个发射天线的发射天线阵列Z=[z1,z2,…,zn,…,zN],每个发射天线zn上连接有时间延迟线TTDn,相邻发射天线之间的距离和时间延迟分别为d和τ,其中,N≥2。
步骤3)合成发射波形s(t,φ):
(3a)设发射天线z1为参考发射天线,并获取第n个发射天线zn的延迟时间(n-1)·τ,以及发射机产生的发射信号g(t):
Figure BDA0002426622140000063
其中,t表示快时间,K表示发射信号g(t)的调频率,K=B/TP,B和TP分别表示发射信号g(t)的带宽和脉宽,fc表示发射信号g(t)的中心频率。参考发射天线选取的不同,会对合成发射波形s(t,φ)有影响。为方便起见,本实施例选择第一个发射天线z1为参考发射天线。
(3b)通过(n-1)·τ和g(t),计算发射机产生的发射信号g(t)经过TTDn时间延迟后的时间延迟信号sn(t),得到时间延迟信号集合S:
S=[s1(t),s2(t),…,sn(t),…,sN(t)]
sn(t)=g(t-(n-1)·τ);
(3c)对S中的所有时间延迟信号在远场进行合成,得到目标的俯仰角为φ的发射波形s(t,φ):
Figure BDA0002426622140000071
其中,c表示光速。在远场对所有时间延迟信号进行合成时,则对于目标来说,波前是平面波,各阵元间的波程差和阵元间距d相关。
步骤4)提取发射波形s(t,φ)中的俯仰维频率脉内扫描信号的时-空方向图pe(t,φ):
计算发射波形s(t,φ)的模sabs(t,φ),并将sabs(t,φ)和发射信号g(t)的商作为时-空方向图pe(t,φ)。对于在远场合成的信号s(t,φ)中包含发射信号和俯仰维频率脉内扫描阵列的时-空方向图,可以通过将发射信号提取分离,从而得到时-空方向图的表达式pe(t,φ):
sabs(t,φ)=|s(t,φ)|,
Figure BDA0002426622140000072
其中,
Figure BDA0002426622140000073
步骤5)计算时-空方向图pe(t,φ)的瞬时波束指向角φpeak(t):
令时-空方向图pe(t,φ)的分母为0,求解瞬时波束指向角φpeak(t),即
Figure BDA0002426622140000074
上式可以改写成
(τ·c-d·sinφ)/λ(t,φ)=m0+m1(t,φ)
其中,
Figure BDA0002426622140000081
m1(t,φ)=τ·K·(t+(N-1)·d·sinφ/2c)-d·sinφ/λ(t,φ)。
当俯仰维角度覆盖范围远小于(-90°,90°),时延τ满足τ<<1/B,所以|m1(t,φ)|<1。因此m的取值为最接近m0的整数,可以取
Figure BDA0002426622140000082
其中,round(·)表示取整操作,0≤t≤TP
从时-空方向图pe(t,φ)的表达式可以看出,时-空方向图接近于Sinc函数形式,则时-空方向图的瞬时波束指向角的位置即为pe(t,φ)的峰值位置,即令分母为0求解得到瞬时波束指向角的表达式为
Figure BDA0002426622140000083
其中,f(t)=fc+K·(t-TP/2)为发射信号g(t)的瞬时频率。
步骤6)计算发射天线阵列Z的发射天线数量N:
(6a)根据SAR成像的测绘带W和测绘带近端的入射角
Figure BDA0002426622140000084
计算测绘带远端的入射角
Figure BDA0002426622140000085
Figure BDA0002426622140000086
根据图2所示的俯仰维频率脉内扫阵列的几何模型示意图,可以看出,测绘带远端的入射角
Figure BDA0002426622140000087
可以通过测绘带W和测绘带近端的入射角
Figure BDA0002426622140000088
并利用三角关系求解。本实施例中,测绘带的远端入射角
Figure BDA0002426622140000089
(6b)根据SAR平台速度V和SAR成像方位分辨率ρa,计算发射信号g(t)的脉冲重复频率Fr
Figure BDA0002426622140000091
在实现高分宽幅成像时,为了达到高的方位分辨率ρa,需要大的多普勒带宽Bd,而方位分辨率ρa和多普勒带宽Bd之间的关系是ρa=V/Bd。为了避免多普勒模糊,则脉冲重复频率Fr和多普勒带宽之间应满足奈奎斯特采样定律,考虑10%的余量,因此可以得到脉冲重复频率Fr的计算公式,这个重频也需要同时满足可以避免星下点回波。本实施例中,脉冲重复频率最后选取为Fr=10.76KHz。
(6c)根据测绘带远端的入射角
Figure BDA0002426622140000092
和发射信号g(t)的脉冲重复频率Fr,计算俯仰维波束宽度θ3dB
Figure BDA0002426622140000093
本实施例利用俯仰维频率脉内扫描实现高分宽幅成像的波形设计是利用高的脉冲重复频率来满足高方位分辨率时,无多普勒模糊,则在整个测绘带内会存在距离模糊。为了保证目标点和其对应的距离模糊分量能被分开,则其对应的主瓣需要被完全分开。因为成像测绘带远端及其模糊分量之间的俯仰角变化范围相比于测绘带近端及其模糊分量之间的俯仰角变化范围要小,所以可以利用测绘带远端及其一次模糊分量之间的俯仰角变化范围作为俯仰维频率脉内扫描阵列的俯仰维波束主瓣宽度。本实施例中得到的俯仰维波束3dB宽度为θ3dB=0.85°。
(6d)根据俯仰维波束宽度θ3dB,计算发射天线阵列Z的发射天线数量N:
Figure BDA0002426622140000094
其中,ceil(·)表示向上取整。在实际应用中,由于阵元个数较多,所需的时间延迟线会较多,可以使用子阵技术,子阵内的各阵元间没有时间延迟,各子阵间有时间延迟。本实施例中,发射天线数量N=25,子阵内有4个阵元。
步骤7)计算发射信号g(t)的信号带宽B:
(7a)根据距离分辨率ρr,确定测绘带的远端俯仰角φfar分配的信号带宽B(φfar):
B(φfar)=2c/ρr
本实施例中的俯仰维频率脉内扫描阵列,每个发射天线发射的是LFM信号,相邻发射天线之间有时间延迟,时域延迟等效于频域加权。因此,不同的脉内时间对应不同的瞬时频率,不同的瞬时频率会被辐射到不同的俯仰角上。在测绘带内,每个点目标仅能被一部分脉内时间照射到,对应所分配到的信号带宽也只占据总的信号带宽的一部分。本实施例中,俯仰维波束从测绘带的远端扫描到近端,远端分配到的信号带宽最少,因此要保持整个测绘带内的距离分辨率,只需要测绘带远端的带宽能满足距离分辨率即可。本实施例中,按照距离分辨率,得到测绘带远端能分配到的信号带宽B(φfar)=200MHz。
(7b)根据测绘带的近端入射角
Figure BDA0002426622140000101
和远端入射角
Figure BDA0002426622140000102
确定俯仰维频率脉内扫描的波束扫描范围θscan
Figure BDA0002426622140000103
为了保证测绘带边缘部分也能分配到足够的信号带宽来满足距离分辨率,应该使俯仰维频率脉内扫描范围大于测绘带内的入射角变化范围,对波束扫描范围放宽一个主瓣宽度。本实施例中计算的波束扫描范围θscan=9.03°
(7c)根据波束扫描范围θscan,计算发射信号g(t)的信号带宽B:
Figure BDA0002426622140000104
采用俯仰维脉内频率扫描,所以测绘带内的某个点目标只会分配到信号的一部分带宽,总的信号带宽和波束扫描范围和波束主瓣宽度有关系。本实施例中,信号带宽B=1200MHz,此值略大于计算得到的信号带宽,可以让测绘带内的距离分辨率更好。
步骤8)计算发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ:
(8a)设定波束扫描范围θscan的起始角度和终止角度分别为θ1和θ2,求解时-空方向图pe(t,φ)的瞬时波束指向角分别为θ1和θ2时的m1和m2
Figure BDA0002426622140000111
脉冲起始时刻t=0,脉内扫描开始,对应的瞬时波束指向角为θ1,脉冲终止时刻t=TP,脉内扫描结束,对应的瞬时波束指向角为θ2,将此条件带入瞬时波束指向角,便可以得到上述方程组。由方程组可以看出,需要通过优化τ、fc、TP和B来实现对θ1和θ2的控制。但通过上述方程组,只能求解两个参数,因此需要在τ、fc、TP和B这四个参数中先确定两个参数,余下的再由上述方程组求解。上述方程组是分式方程,如果求解参数TP,可能存在无解或增根,因此通常将TP设为已知值。为了保持测绘带内的距离分辨率,对信号带宽B有限制。所以在本发明中,通过优化τ和fc来实现对θ1和θ2的控制。本实施例中设定脉宽TP=20us。
(8b)根据给定的期望的发射信号g(t)的中心频率为fc0,并通过假定m1=m2,求解期望的时间延迟τ0
Figure BDA0002426622140000112
求解上述方程时需要确定m1和m2的值,而在实际中m1和m2是未知的。因此,我们需要根据已经确定的TP和B,并给定期望的工作频率fc0,来求解期望的时间延迟τ0。期望的工作频率fc0可以选定为发射信号的中心频率。本实施例中,fc0=10GHz。
(8c)根据期望的时间延迟τ0,得出m1和m2
Figure BDA0002426622140000121
Figure BDA0002426622140000122
将m1和m2带入方程组
Figure BDA0002426622140000123
求解出发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ。本实施例中,求解出中心频率fc=10.706GHz,时间延迟τ=0.28ns。
步骤9)获取俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果:
根据计算得到的发射信号g(t)的信号带宽B和中心频率fc,发射天线阵列Z的发射天线数量N和相邻发射天线之间的时间延迟τ,带入俯仰维频率脉内扫描信号的时-空方向图pe(t,φ),得到俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果。
下面结合仿真实验,对本发明的技术效果作进一步描述。
1.仿真条件和内容:
假设SAR的平台高度为H=600km,速度为V=7500m/s,地面测绘带宽度为W=100km,距离分辨率ρr=0.75m,方位分辨率ρa=0.75m,测绘带近端的入射角为
Figure BDA0002426622140000124
表1列出了在上述参数设计下,根据本发明的俯仰维频率脉内扫描设计方法所获得的阵列参数和波形参数。软件环境:MATLAB仿真软件。
载频 10.706GHz 带宽 1200MHz
脉宽 20us 阵元间延迟 0.28ns
方位天线长度 1.5m 俯仰维阵元个数 25
重频 10.76KHz 俯仰维阵元间距 0.06m
速度 7500km/s 平台高度 600km
表1阵列参数和波形参数设计
软件环境:MATLAB仿真软件。
仿真1,对本发明的时空方向图进行仿真,其结果如图3所示;
仿真2,对本发明的测绘带内的距离分辨率进行仿真,其结果如图4所示;
仿真3,对本发明的测绘带内的点目标及其距离模糊分量分布进行仿真,其结果如图5所示;
仿真4,对本发明的测绘带内的距离模糊比进行仿真,其结果如图6所示。
2.仿真结果分析:
参照图3,图中的横坐标表示脉内时间,纵坐标表示时-空方向图的幅值。从图中可以看出,脉内任意时间的俯仰维频率脉内扫描的时-空方向图近似为Sinc函数。
参照图4,图中的横坐标表示测绘带内的入射角变化,纵坐标表示斜距分辨率。由图中可以看出整个测绘带内的距离分辨率都优于设定的0.75m的指标,在测绘带的近端距离分辨率要好于测绘带远端的距离分辨率,这是因为脉内扫描从测绘带远端开始,终止在测绘带近端,测绘带远端分配到的信号带宽较小。
参照图5,其中图5(a)表示测绘带近端的点目标及其距离模糊分量分布,横轴表示脉内时间,纵轴表示时-空方向图幅值的平方,实线表示测绘带近端的点目标,虚线表示测绘带近端的距离模糊分量。由图中可以看出,在脉冲末端扫描到测绘带近端,并且测绘带近端的点目标和其模糊分量可以分离。图5(b)表示测绘带中部的点目标及其距离模糊分量分布,横轴表示脉内时间,纵轴表示时-空方向图幅值的平方,实线表示测绘带中部的点目标,虚线表示测绘带中部的距离模糊分量。图5(c)表示测绘带远端的点目标及其距离模糊分量分布,横轴表示脉内时间,纵轴表示时-空方向图幅值的平方,实线表示测绘带远端的点目标,虚线表示测绘带远端的距离模糊分量。图5表明,可以通过波形设计来控制瞬时波束指向角指向测绘带内的不同俯仰角处,实现俯仰维频率脉内扫描。
参照图6,图中的横坐标表示测绘带内的入射角变化,纵坐标表示距离模糊比的大小。从图中可以看出在测绘带内,随着入射角的增大,距离模糊比变差,但即便是最坏的情况,在测绘带远端的距离模糊比为24dB,也能满足常规的SAR成像的距离模糊隔离要求。
综上,本发明所提出的方法将高分宽幅成像中出现的距离模糊可以通过距离频域滤波等效实现空域滤波,而不需要依赖DBF技术,通过脉内频率连续扫描,而不需要分割子脉冲和子脉冲间波束控制,工程实现简单。

Claims (3)

1.一种基于俯仰维频率脉内扫描的高分宽幅SAR的波形设计方法,其特征在于,包括如下步骤:
(1)设置输入参数:
根据SAR成像的指标要求,设置SAR成像的测绘带大小为W、距离分辨率为ρr、方位分辨率为ρa,测绘带近端入射角
Figure FDA0002426622130000011
SAR平台的高度和速度分别为H和V;
(2)构建发射天线阵列Z:
构建包括周期性排布的N个发射天线的发射天线阵列Z=[z1,z2,…,zn,…,zN],每个发射天线zn上连接有时间延迟线TTDn,相邻发射天线之间的距离和时间延迟分别为d和τ,其中,N≥2;
(3)合成发射波形s(t,φ):
(3a)设发射天线z1为参考发射天线,并获取第n个发射天线zn的延迟时间(n-1)·τ,以及发射机产生的发射信号g(t):
Figure FDA0002426622130000012
其中,t表示快时间,K表示发射信号g(t)的调频率,K=B/TP,B和TP分别表示发射信号g(t)的带宽和脉宽,fc表示发射信号g(t)的中心频率;
(3b)通过(n-1)·τ和g(t),计算发射机产生的发射信号g(t)经过TTDn时间延迟后的时间延迟信号sn(t),得到时间延迟信号集合S:
S=[s1(t),s2(t),…,sn(t),…,sN(t)]
sn(t)=g(t-(n-1)·τ);
(3c)对S中的所有时间延迟信号在远场进行合成,得到目标的俯仰角为φ的发射波形s(t,φ);
(4)提取发射波形s(t,φ)中的俯仰维频率脉内扫描信号的时-空方向图pe(t,φ):
计算发射波形s(t,φ)的模sabs(t,φ),并将sabs(t,φ)和发射信号g(t)的商作为时-空方向图pe(t,φ);
(5)计算时-空方向图pe(t,φ)的瞬时波束指向角φpeak(t):
令时-空方向图pe(t,φ)的分母为0,求解瞬时波束指向角φpeak(t),即
Figure FDA0002426622130000021
其中,
Figure FDA0002426622130000022
round(·)表示取整操作,0≤t≤TP
(6)计算发射天线阵列Z的发射天线数量N:
(6a)根据SAR成像的测绘带W和测绘带近端的入射角
Figure FDA0002426622130000023
计算测绘带远端的入射角
Figure FDA0002426622130000024
Figure FDA0002426622130000025
(6b)根据SAR平台速度V和SAR成像方位分辨率ρa,计算发射信号g(t)的脉冲重复频率Fr
Figure FDA0002426622130000026
(6c)根据测绘带远端的入射角
Figure FDA0002426622130000027
和发射信号g(t)的脉冲重复频率Fr,计算俯仰维波束宽度θ3dB
Figure FDA0002426622130000031
(6d)根据俯仰维波束宽度θ3dB,计算发射天线阵列Z的发射天线数量N:
Figure FDA0002426622130000032
其中,ceil(·)表示向上取整;
(7)计算发射信号g(t)的信号带宽B:
(7a)根据距离分辨率ρr,确定测绘带的远端俯仰角φfar分配的信号带宽B(φfar):
B(φfar)=2c/ρr
(7b)根据测绘带的近端入射角
Figure FDA0002426622130000033
和远端入射角
Figure FDA0002426622130000034
确定俯仰维频率脉内扫描的波束扫描范围θscan
Figure FDA0002426622130000035
(7c)根据波束扫描范围θscan,计算发射信号g(t)的信号带宽B:
Figure FDA0002426622130000036
(8)计算发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ:
(8a)设定波束扫描范围θscan的起始角度和终止角度分别为θ1和θ2,求解时-空方向图pe(t,φ)的瞬时波束指向角分别为θ1和θ2时的m1和m2
Figure FDA0002426622130000037
(8b)根据给定的期望的发射信号g(t)的中心频率为fc0,并通过假定m1=m2,求解期望的时间延迟τ0
Figure FDA0002426622130000041
(8c)根据期望的时间延迟τ0,得出m1和m2
Figure FDA0002426622130000042
Figure FDA0002426622130000043
将m1和m2带入方程组
Figure FDA0002426622130000044
求解出发射信号g(t)的中心频率fc和相邻发射天线之间的时间延迟τ;
(9)获取俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果:
根据计算得到的发射信号g(t)的信号带宽B和中心频率fc,发射天线阵列Z的发射天线数量N和相邻发射天线之间的时间延迟τ,带入俯仰维频率脉内扫描信号的时-空方向图pe(t,φ),得到俯仰维频率脉内扫描的高分宽幅SAR的波形设计结果。
2.根据权利要求1所述的基于俯仰维频率脉内扫描的高分宽幅SAR的波形设计方法,其特征在于,步骤(3c)中所述的发射波形s(t,φ)的表达式为:
Figure FDA0002426622130000045
其中,c表示光速。
3.根据权利要求1所述的基于俯仰维频率脉内扫描的高分宽幅SAR的波形设计方法,其特征在于,步骤(4)中所述的发射波形s(t,φ)的模sabs(t,φ),以及时-空方向图pe(t,φ),其表达式分别为:
sabs(t,φ)=|s(t,φ)|
Figure FDA0002426622130000051
其中,
Figure FDA0002426622130000052
CN202010222643.6A 2020-03-26 2020-03-26 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法 Active CN111257878B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010222643.6A CN111257878B (zh) 2020-03-26 2020-03-26 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010222643.6A CN111257878B (zh) 2020-03-26 2020-03-26 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法

Publications (2)

Publication Number Publication Date
CN111257878A true CN111257878A (zh) 2020-06-09
CN111257878B CN111257878B (zh) 2022-10-21

Family

ID=70944505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010222643.6A Active CN111257878B (zh) 2020-03-26 2020-03-26 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法

Country Status (1)

Country Link
CN (1) CN111257878B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208798A (zh) * 2019-05-27 2019-09-06 西安空间无线电技术研究所 一种高分宽幅星载马赛克sar成像处理方法及系统
CN112083416A (zh) * 2020-09-18 2020-12-15 电子科技大学 一种运动平台扫描雷达超分辨成像视场选择方法
CN112904340A (zh) * 2021-03-02 2021-06-04 中国电子科技集团公司第三十八研究所 一种合成孔径雷达宽幅扫描实现方法和装置
CN113030970A (zh) * 2021-03-17 2021-06-25 西安电子科技大学 基于方位频扫的二维扫描高分宽幅sar的波形设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2743727A2 (en) * 2014-01-16 2014-06-18 Institute of Electronics, Chinese Academy of Sciences Method for implementing high-resolution wide-swath spaceborne SAR system
CN103885037A (zh) * 2014-04-11 2014-06-25 中国人民解放军国防科学技术大学 基于空时编码的mimo-sar信号发射及接收方法
CN107942327A (zh) * 2017-10-20 2018-04-20 西安电子科技大学 基于脉冲相位编码的单通道hrws‑sar成像方法
EP3432027A1 (en) * 2017-07-20 2019-01-23 Airbus Defence and Space GmbH High resolution wide swath synthetic aperture radar system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2743727A2 (en) * 2014-01-16 2014-06-18 Institute of Electronics, Chinese Academy of Sciences Method for implementing high-resolution wide-swath spaceborne SAR system
CN103885037A (zh) * 2014-04-11 2014-06-25 中国人民解放军国防科学技术大学 基于空时编码的mimo-sar信号发射及接收方法
EP3432027A1 (en) * 2017-07-20 2019-01-23 Airbus Defence and Space GmbH High resolution wide swath synthetic aperture radar system
CN107942327A (zh) * 2017-10-20 2018-04-20 西安电子科技大学 基于脉冲相位编码的单通道hrws‑sar成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU NAN ET AL.: "Signal Modeling and Analysis for Elevation Frequency Scanning HRWS SAR", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
叶俊: "一种高动态平台SAR成像自适应波形设计", 《中国科技信息》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208798A (zh) * 2019-05-27 2019-09-06 西安空间无线电技术研究所 一种高分宽幅星载马赛克sar成像处理方法及系统
CN110208798B (zh) * 2019-05-27 2021-06-11 西安空间无线电技术研究所 一种高分宽幅星载马赛克sar成像处理方法及系统
CN112083416A (zh) * 2020-09-18 2020-12-15 电子科技大学 一种运动平台扫描雷达超分辨成像视场选择方法
CN112904340A (zh) * 2021-03-02 2021-06-04 中国电子科技集团公司第三十八研究所 一种合成孔径雷达宽幅扫描实现方法和装置
CN112904340B (zh) * 2021-03-02 2023-06-20 中国电子科技集团公司第三十八研究所 一种合成孔径雷达宽幅扫描实现方法和装置
CN113030970A (zh) * 2021-03-17 2021-06-25 西安电子科技大学 基于方位频扫的二维扫描高分宽幅sar的波形设计方法

Also Published As

Publication number Publication date
CN111257878B (zh) 2022-10-21

Similar Documents

Publication Publication Date Title
CN111257878B (zh) 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法
JP6437924B2 (ja) 高分解能ストリップマップsarイメージング
US6633253B2 (en) Dual synthetic aperture radar system
CN107741586B (zh) 基于DBF-TOPS加权的星载Ka InSAR信号处理方法
CN109581352B (zh) 一种基于毫米波雷达的超分辨测角系统
CN107942295B (zh) 一种前视阵列sar系统的稀疏天线
CN110412570B (zh) 基于空间脉冲相位编码的hrws-sar成像方法
CN114002673A (zh) 星载无源sar的非合作信号感知系统与多维参数估计方法
Li et al. Radar and inverse scattering
Villano et al. Advanced spaceborne SAR systems with planar antenna
He et al. Ground moving target detection with nonuniform subpulse coding in SAR system
US6650274B1 (en) Radar imaging system and method
CN113030970A (zh) 基于方位频扫的二维扫描高分宽幅sar的波形设计方法
CN114924269B (zh) 一种基于星载f-scan sar的距离向模糊度分析方法
CN115436940A (zh) 一种稀疏滑动聚束sar成像模式实现方法及装置
Zhang et al. Multi-beam Doppler beam sharpening approach for airborne forward-looking radar imaging
Baumgartner et al. Simultaneous low PRF GMTI and HRWS SAR imaging without changing the system operation mode
Bordoni et al. Beam-switch wide swath mode for interferometrically compatible products of the post-sentinel HRWS SAR system
Nguyen et al. Cross-platform Range Ambiguity Rejection Algorithm for Multistatic SAR Constellations in Elevation
CN112180368A (zh) 一种数据处理方法、设备、系统及存储介质
Nan et al. Intrapulse Continuous Azimuth Frequency Scanning-Based Spotlight SAR
Prats et al. A TOPSAR processing algorithm based on extended chirp scaling: Evaluation with TerraSAR-X data
US20030132874A1 (en) Radar system and method
Chen et al. Investigation on an ultra–wide-swath, multiple-elevation-beam SAR based on sweep-PRI mode
Vu Ultrawideband-Ultrawidebeam Synthetic Aperture Radar–Signal Processing and Applications

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