CN103760545A - 合成孔径雷达中子带串扰对成像性能影响的分析方法 - Google Patents

合成孔径雷达中子带串扰对成像性能影响的分析方法 Download PDF

Info

Publication number
CN103760545A
CN103760545A CN201410019949.6A CN201410019949A CN103760545A CN 103760545 A CN103760545 A CN 103760545A CN 201410019949 A CN201410019949 A CN 201410019949A CN 103760545 A CN103760545 A CN 103760545A
Authority
CN
China
Prior art keywords
subband
delta
decoy
data
synthetic
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.)
Pending
Application number
CN201410019949.6A
Other languages
English (en)
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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201410019949.6A priority Critical patent/CN103760545A/zh
Publication of CN103760545A publication Critical patent/CN103760545A/zh
Priority to EP14169232.7A priority patent/EP2762918B1/en
Pending legal-status Critical Current

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/40Means for monitoring or calibrating
    • 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/023Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
    • 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/003Bistatic radar systems; Multistatic radar systems
    • 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

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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明是一种合成孔径雷达中子带串扰对成像性能影响的分析方法,该方法包括步骤如下:步骤S1:构建带有串扰的子带信号合成后的距离向脉冲响应模型,获得子带串扰对成像性能影响的信息;步骤S2:对带有串扰的子带信号进行子带合成仿真分析处理,获得距离向脉冲响应模型的仿真结果,用以验证距离向脉冲响应模型的正确性,并构建子带交叠度与假目标幅度的关系曲线;步骤S3:利用多子带单发单收实测数据模拟多子带多发多收实测数据,并对模拟的多子带多发多收实测数据成像,用以验证子带串扰对模拟的多子带多发多收实测数据成像性能影响的信息。

Description

合成孔径雷达中子带串扰对成像性能影响的分析方法
技术领域
本发明属于多发多收多子带合成孔径雷达技术领域,涉及多发多收合成孔径雷达的信号串扰问题,尤其涉及一种采用步进频线性调频信号的多发多收合成孔径雷达中的子带信号串扰的分析方法。
背景技术
近年来实现二维高分辨率宽测绘带成为星载合成孔径雷达追求的首要目标。方位向分辨率与测绘带宽的矛盾可以由方位向单发多收(SIMO)技术(多通道技术)解决。而距离向高分辨率的实现需要发射大带宽信号,大带宽信号的发射和接收对硬件的要求很高。当硬件水平不能达到要求时,可发射多个步进频的线性调频信号,接收和解调后在地面合成宽带信号。为了不抵消掉方位向多通道技术实现的等效脉冲重复频率的提高,采用多发多收(MIMO)的天线结构,其中每个孔径同时发射不同载频的子带信号。由于各个子带同时发射和接收,如果子带频谱不能完全分离开,将会引起子带间的串扰,可见子带串扰是由子带信号频谱的交叠造成的。由于截断线性调频信号频谱的吉布斯效应,信号带宽外仍存在较小的频谱分量,故即使发射信号带宽等于载频间隔(理论上无交叠),仍存在子带串扰。
子带串扰问题是多发多收多子带合成孔径雷达面临的主要问题。串扰对成像性能的影响分析至关重要,它能为发射信号频谱交叠度的选择提供依据,也能为后期的目标检测提供先验知识。
发明内容
(一)要解决的技术问题
有鉴于此,本发明的主要目的在于从数学推导、仿真实验和实测数据实验三方面得出子带串扰对多发多收多子带合成孔径雷达成像性能的影响特征,为此提供一种合成孔径雷达中子带串扰对成像性能影响的分析方法。
(二)技术方案
为达到上述目的,本发明的多发多收多子带合成孔径雷达中子带串扰对成像性能影响的分析方法包括步骤如下:
步骤S1:构建带有串扰的子带信号合成后的距离向脉冲响应模型,获得子带串扰对成像性能影响的信息;
步骤S2:对带有串扰的子带信号进行子带合成仿真分析处理,获得距离向脉冲响应模型的仿真结果,用以验证距离向脉冲响应模型的正确性,并构建子带交叠度与假目标幅度的关系曲线;
步骤S3:利用多子带单发单收实测数据模拟多子带多发多收实测数据,并对模拟的多子带多发多收实测数据成像,用以验证子带串扰对模拟的多子带多发多收实测数据成像性能影响的信息。
(三)有益效果
本发明公开了多发多收多子带合成孔径雷达中子带串扰对成像性能的影响分析结果,构建了带有串扰的子带回波子带合成后的脉冲响应模型;对该模型进行仿真实验,得到交叠度与子带串扰造成的假目标幅度的关系曲线;利用多子带单发单收实测数据模拟多子带多发多收模式,并对其成像,验证了子带串扰对实测数据成像性能的影响。
本发明构建了带有子带串扰的距离向脉冲响应函数,得知子带交叠度越大,串扰的影响越大,脉冲响应函数中的假目标幅度越大,同时计算出了假目标的具体位置;仿真实验考虑了信号频谱的吉布斯效应,给出不同带宽下交叠度与假目标幅度的关系曲线;实测数据实验给出了31.5MHz子带带宽下,子带交叠5%和子带无交叠(仍存在吉布斯效应造成的子带串扰)时的成像结果,进一步验证了子带串扰对实测数据成像性能的影响。本发明为多子带并发模式的子带交叠度选择提供了依据,同时发明中假目标位置的确定也为后期目标检测提供了先验知识。
附图说明
图1a为本发明完整构思的流程图;
图1b为本发明中多发多收多子带模式下的两发四收几何模型和收发时序图;
图2为本发明中有子带交叠时的回波分离示意图;
图3为本发明中无子带交叠时的回波分离示意图;
图4a为本发明中不同交叠度下,子带合成后的脉冲压缩结果原图;
图4b为本发明中不同交叠度下,子带合成后的脉冲压缩结果中真实目标的放大图;
图5为本发明中不同合成带宽下,子带交叠度与假目标幅度的关系曲线;
图6为本发明中各个处理步骤后的等效相位中心分布示意图;
图7为本发明中多发多收多子带合成孔径雷达的成像处理框图;
图8为本发明中子带交叠5%时的实测数据成像结果;
图9为本发明中子带交叠0%(即只考虑吉布斯效应引起的串扰)时的实测数据成像结果。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图1a为本发明对多发多收多子带合成孔径雷达中子带串扰的分析方法的流程图;该方法包括步骤如下:
步骤S1:构建带有串扰的子带信号合成后的距离向脉冲响应模型,获得子带串扰对成像性能影响的信息;
步骤S2:对带有串扰的子带信号进行子带合成仿真分析处理,获得距离向脉冲响应模型的仿真结果,用以验证距离向脉冲响应模型的正确性,并构建子带交叠度与假目标幅度的关系曲线;
步骤S3:利用多子带单发单收实测数据模拟多子带多发多收实测数据,并对模拟的多子带多发多收实测数据成像,用以验证子带串扰对模拟的多子带多发多收实测数据成像性能影响的信息。
步骤S1中,构建距离向脉冲响应模型的步骤如下:
步骤S11:根据串扰的频谱分量,得出带有串扰的子带信号频谱;
步骤S12:对带有串扰的子带信号频谱在距离向匹配滤波,获得距离压缩后的子带信号频谱;
步骤S13:对子带信号频谱在距离频域补零,获得距离向升采样的子带信号频谱;
步骤S14:对子带信号频谱做逆傅里叶变换,获得时域子带信号,对时域子带信号乘以对应的频移函数,获得频移后的时域子带信号;
步骤S15:将频移后的时域子带信号相干叠加,获得带有串扰的子带信号合成后的距离向脉冲响应模型。
图1b为本发明中多发多收多子带模式下的两发四收几何模型和收发时序图。图1b中打为轨道高度,vs为平台运动速度。发射时,子孔径1发射子带1,子孔径2发射子带2;接收时,每个子孔径同时接收两个子带回波。则每次发射接收可形成四个等效相位中心,子带1对应的等效相位中心位置如图1b中①所示,子带2对应的等效相位中心位置如图1b中②所示。
有子带交叠时的回波分离示意图如图2所示,以得到子带1的基带回波信号为例,将接收到的回波信号用子带1的载频f1解调,则子带1的中心频率变为O,子带2的中心频率变为f2-f1,其中f2为子带2的载频;再经过一个模拟低通滤波器(其带宽大于信号带宽),如图中梯形所示;再经过AD采样并压缩后传到地面;对读取的基带回波信号在数字域截频,数字截频滤波器如图中矩形框所示。截频滤波器的使用既可以将串扰的影响控制到最小,又可以使无串扰的合成频谱幅度连续。由于三子带或更多子带问的串扰分析与两子带情况类似,本发明只给出两子带问的串扰分析方法。
考虑两子带的情况,理论上,带有串扰后的子带1基带信号频谱为:
S 1 ′ ′ ( f ) = S 1 ( f ) + ΔS 2 ( f )
= exp { - j 2 πτ a ( f 1 + f ) } · exp { - jπ f 2 k r } · rect ( f + ΔB / 4 Δf + ΔB / 2 ) - - - ( 1 )
+ exp { - j 2 πτ a ( f 2 + f - Δf ) } · exp { - jπ ( f - Δf ) 2 k r } · rect ( f - Δf 2 + ΔB 4 ΔB / 2 )
公式(1)中f为距离向频率,f1为子带1的载频,f2为子带2的载频,j为虚数单位,π为圆周率,τa为当前方位位置对应的延时,τ为距离向时间变量,Δf为频率步进,ΔB=Bs-Δf为子带交叠部分的信号带宽,Bs为单子带带宽,kr为发射脉冲的调频率。第一项S1(f)为无串扰时合成的理想时域波形,第二项ΔS2(f)为子带2对子带1的串扰,串扰部分的带宽为ΔB/2。
同样的,带有串扰的子带2的基带信号频谱为:
Figure BDA0000457502630000051
= exp { - j 2 πτ a ( f 2 + f ) } · exp { - jπ f 2 k r } · rect ( f - ΔB / 4 Δf + ΔB / 2 ) - - - ( 2 )
+ exp { - j 2 πτ a ( f 1 + f - Δf ) } · exp { - jπ ( f - Δf ) 2 k r } · rect ( f - Δf 2 + ΔB 4 ΔB / 2 ) .
公式(2)中ΔS1(f)为子带1对子带2的串扰。
匹配滤波、升采样和频移后,这两个子带信号的频谱为:
S 1 ″ ( f + Δf 2 ) = exp { - j 2 πτ a ( f c + f ) }
× [ rect ( f + ΔB / 4 + Δf / 2 Δf + ΔB / 2 ) + exp { j 2 πf Δf k r } · rect ( f + ΔB 4 ΔB / 2 ) ] - - - ( 3 )
S 2 ″ ( f - Δf 2 ) = exp { - j 2 πτ a ( f c + f ) }
× [ rect ( f - ΔB / 4 - Δf / 2 Δf + ΔB / 2 ) + exp { j 2 πf Δf k r } · rect ( f - ΔB 4 ΔB / 2 ) ] . - - - ( 4 )
其中,fc为合成信号的载频。
公式(3)和(4)相加得到叠加后的频谱为:
Figure BDA0000457502630000061
× [ rect ( f Δf + B s ) + exp { j 2 πf Δf k r } · rect ( f + ΔB 4 ΔB / 2 ) + exp { - j 2 πf Δf k r } · rect ( f - ΔB 4 ΔB / 2 ) ] . - - - ( 5 )
傅里叶逆变换后得到合成后的时域波形,步骤S1中所述距离向脉冲响应模型
Figure BDA0000457502630000067
表示如下:
S w ′ ′ ( τ ) = B w · exp { - j 2 πτ a f c } · sin c { B w ( τ - τ a ) } + ΔB 2 exp { - j 2 πτ a f c } - - - ( 6 )
× exp { jπ ( ΔB 2 τ a - ΔB 2 Δf k r ) - j ΔB 4 τ } sin c { ΔB 2 ( τ - ( τ a - Δf k r ) ) } + exp { jπ ( ΔB 2 Δf k r - ΔB 2 τ a ) + j ΔB 4 τ } sin c { ΔB 2 ( τ - ( τ a + Δf k r ) ) }
其中在距离向脉冲响应模型
Figure BDA0000457502630000066
的加号前为第一项、加号后为第二项,Bw为合成后信号的总带宽;第一项为无串扰时合成的理想时域波形,则合成后的距离分辨率提高为c/(2Bw),c为光速;第二项为串扰造成的附加项,则串扰造成了假目标,假目标的相对幅度与交叠度正比,假目标的位置在离真实目标±c·Δf/(2kr)处。
即使同时发射的各个子带频谱没有交叠,由于信号频谱的吉布斯效应,同时接收的子带信号频谱也存在串扰,如图3所示为本发明中无子带交叠时的回波分离示意图。但这时串扰的频谱幅度小(-6dB以下),串扰的影响小。理论推导无法考虑信号频谱的吉布斯效应(数学表达式复杂),将在仿真实验和实测数据实验中进一步精确地验证子带串扰对成像性能的影响。
步骤S2中所述子带合成仿真分析处理的步骤包括:
步骤S21:根据信号频谱的吉布斯效应,构建带有串扰的两子带基带脉冲信号;
步骤S22:对两子带基带脉冲信号进行子带合成,检测子带合成信号中的假目标,获得假目标的幅度;
步骤S23:改变两子带基带脉冲信号的子带带宽和频谱交叠度,获得不同带宽下假目标幅度;
步骤S24:绘出假目标幅度与频谱交叠度的关系曲线,得到的关系曲线示出:子带带宽越大,吉布斯效应引起的假目标幅度越小;子带频谱的交叠度越大,假目标幅度越大,由此得到的假目标幅度与脉冲响应模型中的假目标幅度越接近。
对于合成带宽为60MHz的情况,图4a为本发明中不同交叠度下,子带合成后的距离向脉冲压缩结果原图,图4b为本发明中不同交叠度下,子带合成后的距离向脉冲压缩结果中真实目标的放大图;交叠度为10%时,单子带带宽33MHz,频率步进30MHz,kr=1.5×1012。假目标的幅度为20·1g(ΔB/Bw)=20·1g(1.5/60)=-32dB。交叠度为0%时,单子带带宽30MHz,频率步进30MHz,吉布斯效应造成的假目标幅度为-45dB。假目标的位置在离真实目标±cΔf/2kr=±3km处。
对于不同的合成带宽,6dB以下的频谱分量占总的信号频谱能量的比例不同,吉布斯效应导致的子带串扰的影响不同。下面给出几种合成信号带宽下,交叠度(重叠率)与假目标幅度的关系,如图5所示。由图5可见,合成带宽为120M时,吉布斯效应造成的假目标幅度为-47dB,合成带宽为240MHz时,吉布斯效应造成的假目标为-49.5dB,合成带宽为480MHz时,吉布斯效应造成的假目标幅度为-53dB。故合成信号带宽越大,吉布斯效应造成的串扰对成像性能的影响越小。另外,对于相同的合成带宽,交叠度越大,子带串扰造成的假目标幅度越大。例如对于两子带合成480MHz带宽,如果要求假目标幅度不超过-50dB,那么子带频谱交叠度不能超过0.8%即20MHz。
步骤S3中所述利用多子带单发单收实测数据模拟多子带多发多收实测数据的步骤如下:
步骤S31:对一个子带数据进行不同的方位延时,对方位延时后的子带数据降采样,获得这个子带的一发多收数据;
步骤S32:对其它的子带数据进行不同的方位延时,对方位延时后的子带数据降采样,获得其它子带的一发多收数据;
步骤S33:结合所有子带的一发多收数据,并加入串扰频谱分量,获得模拟的多子带多发多收实测数据;
步骤S34:对模拟的多子带多发多收实测数据进行成像,获得带有串扰的实测数据图像。
步骤S3中对模拟的多子带多发多收实测数据成像的步骤包括:
步骤S3a:对单个子带的一发多收数据进行多通道重建,获得单子带的方位向均匀分布数据;
步骤S3b:对齐各子带方位向均匀分布数据的等效相位中心;
步骤S3c:对等效相位中心对齐后的子带数据进行子带合成,得到子带合成后的数据;
步骤S3d:对子带合成数据进行常规成像,获得子带串扰对模拟的多子带多发多收实测数据成像性能的影响。
由于目前国内无多发多收多子带的实测合成孔径雷达数据,为了验证子带串扰对实测数据成像性能的影响,用顺发的两子带实测数据模拟两发两收的多发多收多子带实测数据,模拟方法为:
对接收到的各子带信号在进行相应的方位延时后,进行两倍降采样。子孔径n接收的子带m对应的方位延时函数Hm,n(fη)表示如下:
Hm,n(fη)=exp{-j2πfη·tm,n}  (7)其中fη为方位多普勒频率,n表示子孔径个数的序号,m表示子带个数的序号,tm,n表示子带m子孔径n对应的方位延时。对于两子孔径两子带的情况,tm,n表示如下:
t 1,1 = 0 t 1,2 = L sub 2 v s = 3 4 · 1 2 PRF t 2,1 = t 1,2 - 1 4 PRF t 2,2 = t 1,2 + t 2,1 - - - ( 8 )
其中Lsub和PRF分别为模拟的多发多收多子带合成孔径雷达的天线子孔径长度和雷达脉冲重复频率,再根据式(1)和式(2)人为地加入串扰频谱分量,就模拟出了多发多收多子带回波。模拟处理前和处理后的等效相位中心分布如图6第一行和第二行所示,其中①表示子带1的等效相位中心,②表示子带2的等效相位中心。
接下来以两发四收两子带为例简要介绍多发多收多子带合成孔径雷达的成像流程,如图7所示:
步骤S3a对单个子带的一发两收数据进行多通道重建,获得单子带的方位向均匀分布数据。多通道重建后的等效相位中心分布如图6第三行所示;
步骤S3b对齐各子带方位向均匀分布数据的等效相位中心,使得用于合成的子带相位中心重合。这一步通过对子带2乘函数Hc(fη)实现,Hc(fη)的表达式为:
Hc(fη)=exp{-j2πfηΔη}  (9)其中
Figure BDA0000457502630000091
对齐后的等效相位中心分布如图6第四行所示;
步骤S3c对等效相位中心对齐后的子带数据进行子带合成,得到子带合成后的数据。每个方位时刻的子带合成方法与S1中距离向脉冲响应模型的构建方法一样。子带合成后的等效相位中心分布如图6第五行所示。
步骤S4d对子带合成数据进行常规成像,获得高分辨率合成孔径雷达图像。从获得的雷达图像中可以观测到子带串扰对成像质量的影响。
现有两子带的顺发实测数据,该顺发实测数据的子带带宽为31.5MHz,频率步进为30MHz,子带交叠度为5%。按照上述的处理方法得到模拟的多发多收多子带回波数据,再根据图7的成像流程成像,得到的图像如图8所示子带交叠5%的实测数据成像结果。为了模拟吉布斯效应造成的串扰对成像性能的影响,首先对单子带信号频谱利用一个吉布斯形状的滤波器滤波,得到30MHz带宽的子带信号频谱。再按照上述处理方法得到模拟的多发多收多子带回波数据,再根据图7的成像流程成像,得到的图像如图9所示子带交叠0%(即只考虑吉布斯效应引起的串扰)时的实测数据成像结果。从图8和图9可以看出,在交叠度为5%时,假目标在图像中非常明显;在交叠度为0%时,吉布斯效应引起的假目标在图像中不是很明显,但在暗背景中,仍可以区分。对于更大的合成带宽,吉布斯效应引起的假目标在图像中会更不明显。可见,多发多收多子带模式下,子带频谱无交叠时的假目标幅度最小,但仍然不能完全忽略。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内。

Claims (6)

1.合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,该方法包括步骤如下:
步骤S1:构建带有串扰的子带信号合成后的距离向脉冲响应模型,获得子带串扰对成像性能影响的信息;
步骤S2:对带有串扰的子带信号进行子带合成仿真分析处理,获得距离向脉冲响应模型的仿真结果,用以验证距离向脉冲响应模型的正确性,并构建子带交叠度与假目标幅度的关系曲线;
步骤S3:利用多子带单发单收实测数据模拟多子带多发多收实测数据,并对模拟的多子带多发多收实测数据成像,用以验证子带串扰对模拟的多子带多发多收实测数据成像性能影响的信息。
2.根据权利要求1所述合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,构建距离向脉冲响应模型的步骤如下:
步骤S11:根据串扰的频谱分量,得出带有串扰的子带信号频谱;
步骤S12:对带有串扰的子带信号频谱在距离向匹配滤波,获得距离压缩后的子带信号频谱;
步骤S13:对子带信号频谱在距离频域补零,获得距离向升采样的子带信号频谱;
步骤S14:对子带信号频谱做逆傅里叶变换,获得时域子带信号,对时域子带信号乘以对应的频移函数,获得频移后的时域子带信号;
步骤S15:将频移后的时域子带信号相干叠加,获得带有串扰的子带信号合成后的距离向脉冲响应模型。
3.根据权利要求1所述合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,所述距离向脉冲响应模型
Figure FDA0000457502620000011
表示如下:
S w ′ ( τ ) = B w · exp { - j 2 πτ a f c } · sin c { B w ( τ - τ a ) }
+ ΔB 2 exp { - j 2 πτ a f c } · exp { jπ ( ΔB 2 τ a - ΔB 2 Δf k r ) - j ΔB 4 τ } · sin c { ΔB 2 ( τ - ( τ a - Δf k r ) ) } + exp { jπ ( ΔB 2 Δf k r - ΔB 2 τ a ) + j ΔB 4 τ } · sin c { ΔB 2 ( τ - ( τ a + Δf k r ) ) } 其中在距离向脉冲响应模型
Figure FDA0000457502620000021
的加号前为第一项、加号后为第二项,Bw为合成后信号的总带宽,j为虚数单位,τa为当前方位位置对应的距离延时,fc为合成信号的载频,τ为距离向时间变量,ΔB为子带交叠部分的信号带宽,Δf为频率步进,kr为发射脉冲的调频率;第一项为无串扰时合成的理想时域波形,则合成后的距离分辨率提高为c/(2Bw),c为光速;第二项为串扰造成的附加项,则串扰造成了假目标,假目标的相对幅度与交叠度正比,假目标的位置在离真实目标±c·Δf/(2kr)处。
4.根据权利要求1所述合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,所述子带合成仿真分析处理的步骤包括:
步骤S21:根据信号频谱的吉布斯效应,构建带有串扰的两子带基带脉冲信号;
步骤S22:对两子带基带脉冲信号进行子带合成,检测子带合成信号中的假目标,获得假目标的幅度;
步骤S23:改变两子带基带脉冲信号的子带带宽和频谱交叠度,获得不同带宽、不同交叠度下的假目标幅度;
步骤S24:绘出假目标幅度与频谱交叠度的关系曲线,得到的关系曲线示出:子带带宽越大,吉布斯效应引起的假目标幅度越小;子带频谱的交叠度越大,假目标幅度越大,由此得到的假目标幅度与脉冲响应模型中的假目标幅度越接近。
5.根据权利要求1所述合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,所述利用多子带单发单收实测数据模拟多子带多发多收实测数据的步骤如下:
步骤S31:对一个子带数据进行不同的方位延时,对方位延时后的子带数据降采样,获得这个子带的一发多收数据;
步骤S32:对其它的子带数据进行不同的方位延时,对方位延时后的子带数据降采样,获得其它子带的一发多收数据;
步骤S33:结合所有子带的一发多收数据,并加入串扰频谱分量,获得模拟的多子带多发多收实测数据;
步骤S34:对模拟的多子带多发多收实测数据进行成像,获得带有串扰的实测数据图像。
6.根据权利要求1所述合成孔径雷达中子带串扰对成像性能影响的分析方法,其特征在于,对模拟的多子带多发多收实测数据成像的步骤包括:
步骤S3a:对单个子带的一发多收数据进行多通道重建,获得单子带的方位向均匀分布数据;
步骤S3b:对齐各子带方位向均匀分布数据的等效相位中心;
步骤S3c:对等效相位中心对齐后的子带数据进行子带合成,得到子带合成后的数据;
步骤S3d:对子带合成数据进行常规成像,获得子带串扰对模拟的多子带多发多收实测数据成像性能影响的信息。
CN201410019949.6A 2014-01-16 2014-01-16 合成孔径雷达中子带串扰对成像性能影响的分析方法 Pending CN103760545A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201410019949.6A CN103760545A (zh) 2014-01-16 2014-01-16 合成孔径雷达中子带串扰对成像性能影响的分析方法
EP14169232.7A EP2762918B1 (en) 2014-01-16 2014-05-21 Method for analyzing effect of sub-band interference on imaging performance in synthetic aperture radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410019949.6A CN103760545A (zh) 2014-01-16 2014-01-16 合成孔径雷达中子带串扰对成像性能影响的分析方法

Publications (1)

Publication Number Publication Date
CN103760545A true CN103760545A (zh) 2014-04-30

Family

ID=50527812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410019949.6A Pending CN103760545A (zh) 2014-01-16 2014-01-16 合成孔径雷达中子带串扰对成像性能影响的分析方法

Country Status (2)

Country Link
EP (1) EP2762918B1 (zh)
CN (1) CN103760545A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338731A (zh) * 2016-08-22 2017-01-18 西安电子科技大学 一种多子带并发的mimo‑sar雷达成像方法
CN108646242A (zh) * 2018-05-10 2018-10-12 北京航空航天大学 一种针对复杂目标的多子波段雷达数据融合成像方法
CN109274628A (zh) * 2017-07-17 2019-01-25 普天信息技术有限公司 一种多子带系统中的上行业务发送方法及装置
CN110988875A (zh) * 2019-11-14 2020-04-10 西安空间无线电技术研究所 一种超高分辨率星载sar系统多子带频谱拼接合成方法
CN111419211A (zh) * 2019-01-09 2020-07-17 成都心吉康科技有限公司 抑制ecg信号中伪吉布斯效应的方法及装置
CN112098954A (zh) * 2020-08-14 2020-12-18 西安电子科技大学 基于非线性调频信号逼近期望功率谱的发射波形优化方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6400494B2 (ja) * 2015-01-23 2018-10-03 株式会社東芝 レーダ装置及びそのレーダ信号処理方法
CN106483511B (zh) * 2015-08-28 2018-12-11 中国电力科学研究院 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法
CN111638493B (zh) * 2020-05-14 2023-04-11 湖北航天技术研究院总体设计所 一种sar成像导引头干扰对抗仿真系统
CN112272064B (zh) * 2020-09-29 2021-07-06 电子科技大学 一种合作式mimo雷达的检测概率和互信息计算方法
CN112346021B (zh) * 2020-10-30 2023-06-30 西安电子科技大学 Fda-mimo雷达基于频率步进量控制的抗干扰方法
CN112946645B (zh) * 2021-01-29 2022-10-14 北京理工大学 一种无人机载超宽带sar自聚焦方法
CN113176568B (zh) * 2021-03-31 2023-11-10 中国科学院空天信息创新研究院 一种回波信号分离方法、装置及存储介质
CN113203992A (zh) * 2021-04-29 2021-08-03 电子科技大学 一种多基sar的抗欺骗性干扰方法
CN114252855B (zh) * 2022-03-01 2022-05-06 南京雷电信息技术有限公司 一种同时sar/gmt回波模拟方法及其系统
CN115128608B (zh) * 2022-08-31 2022-11-18 南京隼眼电子科技有限公司 车载多普勒分集mimo合成孔径雷达时域成像方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070152874A1 (en) * 2005-12-30 2007-07-05 Woodington Walter G Reducing undesirable coupling of signal(s) between two or more signal paths in a radar system
CN101452073A (zh) * 2007-11-30 2009-06-10 清华大学 基于多发多收频分雷达的宽带信号合成方法
EP2228667A1 (en) * 2009-03-11 2010-09-15 BAE Systems PLC Sensor for determining velocity
CN102393513A (zh) * 2011-08-31 2012-03-28 中国电子科技集团公司第三十八研究所 一种基于自然分布场景和稀少定标器的极化定标技术
CN102792180A (zh) * 2010-03-08 2012-11-21 罗伯特·博世有限公司 用于减小数字信号对hf雷达信号的串扰的雷达传感器及其控制方法
US20130236115A1 (en) * 2012-03-09 2013-09-12 U.S. Government As Represented By The Secretary Of The Army Method and System for Recovery of Missing Spectral Information in Wideband Signal

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102013216461A1 (de) * 2013-08-20 2014-06-26 Deutsches Zentrum für Luft- und Raumfahrt e.V. Synthetik-Apertur-Radarverfahren

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070152874A1 (en) * 2005-12-30 2007-07-05 Woodington Walter G Reducing undesirable coupling of signal(s) between two or more signal paths in a radar system
CN101452073A (zh) * 2007-11-30 2009-06-10 清华大学 基于多发多收频分雷达的宽带信号合成方法
EP2228667A1 (en) * 2009-03-11 2010-09-15 BAE Systems PLC Sensor for determining velocity
CN102792180A (zh) * 2010-03-08 2012-11-21 罗伯特·博世有限公司 用于减小数字信号对hf雷达信号的串扰的雷达传感器及其控制方法
CN102393513A (zh) * 2011-08-31 2012-03-28 中国电子科技集团公司第三十八研究所 一种基于自然分布场景和稀少定标器的极化定标技术
US20130236115A1 (en) * 2012-03-09 2013-09-12 U.S. Government As Represented By The Secretary Of The Army Method and System for Recovery of Missing Spectral Information in Wideband Signal

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XIULIAN LUO等: "Influences of Channel Errors and Interference on the OFDM-MIMO SAR", 《2013 IEEE RADAR CONFERENCE》 *
邓云凯等: "一种基于频域子带合成的多发多收高分辨率SAR成像算法", 《电子与信息学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338731A (zh) * 2016-08-22 2017-01-18 西安电子科技大学 一种多子带并发的mimo‑sar雷达成像方法
CN106338731B (zh) * 2016-08-22 2019-08-02 西安电子科技大学 一种多子带并发的mimo-sar雷达成像方法
CN109274628A (zh) * 2017-07-17 2019-01-25 普天信息技术有限公司 一种多子带系统中的上行业务发送方法及装置
CN109274628B (zh) * 2017-07-17 2021-08-03 普天信息技术有限公司 一种多子带系统中的上行业务发送方法及装置
CN108646242A (zh) * 2018-05-10 2018-10-12 北京航空航天大学 一种针对复杂目标的多子波段雷达数据融合成像方法
CN108646242B (zh) * 2018-05-10 2022-02-08 北京航空航天大学 一种针对复杂目标的多子波段雷达数据融合成像方法
CN111419211A (zh) * 2019-01-09 2020-07-17 成都心吉康科技有限公司 抑制ecg信号中伪吉布斯效应的方法及装置
CN111419211B (zh) * 2019-01-09 2022-09-30 成都心吉康科技有限公司 抑制ecg信号中伪吉布斯效应的方法及装置
CN110988875A (zh) * 2019-11-14 2020-04-10 西安空间无线电技术研究所 一种超高分辨率星载sar系统多子带频谱拼接合成方法
CN112098954A (zh) * 2020-08-14 2020-12-18 西安电子科技大学 基于非线性调频信号逼近期望功率谱的发射波形优化方法
CN112098954B (zh) * 2020-08-14 2024-04-30 西安电子科技大学 基于非线性调频信号逼近期望功率谱的发射波形优化方法

Also Published As

Publication number Publication date
EP2762918A3 (en) 2014-10-22
EP2762918B1 (en) 2016-01-06
EP2762918A2 (en) 2014-08-06

Similar Documents

Publication Publication Date Title
CN103760545A (zh) 合成孔径雷达中子带串扰对成像性能影响的分析方法
CN103728618B (zh) 一种高分辨率、宽测绘带的星载sar体制实现方法
CN104166141B (zh) 基于子带合成的多输入多输出合成孔径雷达系统设计方法
CN104714231B (zh) 一种基于完全互补序列和相位补偿的mimo sar成像方法
CN106405541A (zh) 全相参连续波多普勒雷达及其测距测速方法
CN103605112B (zh) 多发多收干涉合成孔径雷达时频二维信号波形设计方法
CN103472449B (zh) 一种基于mimo图像域的bp宽带合成方法
CN103308892B (zh) 一种多扫频雷达信号的产生及压缩方法
CN107678028A (zh) 低信噪比条件下的微波凝视关联成像方法
CN107544068A (zh) 一种基于频域bp的图像域宽带合成方法
CN103412295A (zh) 基于回波精确模型的高速机动弱目标检测方法
CN103969644A (zh) 一种多通道调频连续波sar成像方法
CN106932778A (zh) 方位多通道fmcw sar滑动聚束成像方法
CN103454630A (zh) 一种基于多阵元发射技术的超宽带三维成像方法
CN109239675A (zh) 一种基于同带宽不同调频斜率的分布式星载雷达lfm波形确定方法
CN106842198A (zh) 一种基于凸优化理论的子带融合成像算法
Yu et al. Ground moving target motion parameter estimation using Radon modified Lv's distribution
CN106680791A (zh) 基于宽带扫频数据的雷达回波仿真方法
CN107102328B (zh) 基于fpga的实时成像信号处理方法及fpga
CN105510878A (zh) 一种分频段处理的多水下自主航行器被动相互定位方法
CN110286360B (zh) 一种基于固定距离延迟的星载sar回波模拟及成像方法
CN103926586B (zh) 一种使用发射子阵的mimo阵列测深方法
Karakasiliotis et al. Two-dimensional ISAR model and image reconstruction with stepped frequency-modulated signal
CN106093941A (zh) 一种基于fmcw雷达系统的冰层截面成像方法
CN102998659B (zh) 基于脉间调制的多普勒频谱赋形方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140430