CN105717508A - 一种基于发射波形方位向调制的机载雷达前视成像方法 - Google Patents

一种基于发射波形方位向调制的机载雷达前视成像方法 Download PDF

Info

Publication number
CN105717508A
CN105717508A CN201610065516.3A CN201610065516A CN105717508A CN 105717508 A CN105717508 A CN 105717508A CN 201610065516 A CN201610065516 A CN 201610065516A CN 105717508 A CN105717508 A CN 105717508A
Authority
CN
China
Prior art keywords
airborne radar
signal
modulation
echo
waveform
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
CN201610065516.3A
Other languages
English (en)
Other versions
CN105717508B (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 CN201610065516.3A priority Critical patent/CN105717508B/zh
Publication of CN105717508A publication Critical patent/CN105717508A/zh
Application granted granted Critical
Publication of CN105717508B publication Critical patent/CN105717508B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/411Identification of targets based on measurements of radar reflectivity

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

本发明公开了一种基于发射波形方位向调制的机载雷达前视成像方法,其主要思路为:获取方位向调制波形的调制序列波形,并对机载雷达在距离向发射的线性调频信号在方位向进行缓慢调制后,再接收经过方位向缓慢调制后的发射信号波形的回波信号,并对所述回波信号依次进行距离向脉冲压缩和匹配滤波,得到距离脉压后的数据,进而计算得到距离徙动校正后的数据,然后再计算得到机载雷达的多普勒中心估计值后,对所述线性调频信号的回波数据进行多普勒中心补偿,进而计算得到第m个距离单元的脉冲表达式后进行方位向匹配滤波处理,获得到匹配滤波后的机载雷达成像,将所述匹配滤波后的机载雷达成像转化为大地坐标系下,即可得到前视场景下的机载雷达成像。

Description

一种基于发射波形方位向调制的机载雷达前视成像方法
技术领域
本发明属于机载雷达成像技术领域,特别涉及一种基于发射波形方位向调制的机载雷达前视成像方法,适用于机载雷达的对地侦察攻击、自主导航、海洋搜索以及导弹末端匹配制导,也可适用于非前视情况下的地形测绘和广域监视。
背景技术
机载雷达通过发射大时宽带宽积信号提高机载雷达成像的距离向分辨率,并利用机载雷达的载机平台与地面非合作目标的相对运动带来的多普勒频域提高机载雷达成像的方位向分辨率,具有全天时、全天候、高分辨的特点。但是,当机载雷达工作在前视条件下时,机载雷达的载机航迹方向与机载雷达的天线波束指向重合,再加上机载雷达回波信号的多普勒变化率非常小,致使机载雷达成像的方位分辨率急剧下降。
目前,常用的合成孔径雷达(SAR)和多普勒波束锐化(DBS)技术等方法失效,分别会形成机载雷达成像的前视探测盲区,使得无法对机载雷达的前视场景进行成像。对于机载雷达的载机着陆时对前方跑道的成像,或者导弹下落时对打击区域的成像,通常采用实波束成像方法或单脉冲前视成像方法完成,但实波束成像方法或单脉冲前视成像方法的方位分辨率要么依赖于波束主瓣宽度,要么对回波信号的信噪比要求较高,尤其复杂环境下时更难满足实时性要求。
发明内容
针对上述现有技术的不足,本发明的目的在于提出一种基于发射波形方位向调制的机载雷达前视成像方法,该种基于发射波形方位向调制的机载雷达前视成像方法将发射的机载雷达信号波形在方位向进行调制,用于增加发射信号的方位向带宽,并且在接收机载雷达回波信号时进行方位向的匹配滤波,最终获得前视场景下的机载雷达高分辨成像。
为实现上述技术目的,本发明采用如下技术方案予以实现。
一种基于发射波形方位向调制的机载雷达前视成像方法,包括以下步骤:
步骤1,获取方位向调制波形,得到所述方位向调制波形的调制序列波形;
步骤2,机载雷达在距离向发射线性调频信号,同时根据所述方位向调制波形的调制序列波形对所述线性调频信号在方位向进行缓慢调制,得到经过方位向缓慢调制后的发射信号波形其中,表示快时间,tn表示第n时刻;
步骤3,机载雷达接收第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号,并对所述第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号依次进行距离向的脉冲压缩和匹配滤波处理,得到第q个相干积累时间间隔对应的距离脉压回波信号;其中,q∈{1,2,…,Q},Q表示机载雷达天线在一个周期内扫描的相干积累时间间隔总个数;
步骤4,对第q个相干积累时间间隔对应的距离脉压回波信号在距离频域内乘以距离走动校正因子H(f,t),得到第q个相干积累时间间隔对应的距离徙动校正回波信号;其中,f表示距离频率,t表示慢时间;
步骤5,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行相关处理并求其平均值,获得第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd
步骤6,根据第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行多普勒中心补偿,得到第q个相干积累时间间隔对应的多普勒中心补偿回波信号;
步骤7,重复执行步骤3-步骤6,直到得到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,此时得到第1个相干积累时间间隔对应的多普勒中心补偿回波信号到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,并将其按照Q个相干积累时间间隔各自经过方位向缓慢调制后的发射信号波形的回波信号接收顺序进行存储,得到机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号,其中每一个相干积累时间间隔对应的回波信号包含N个脉冲,并将机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号划分为M个距离单元,然后计算得到第m个距离单元的脉冲表达式S(m);m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;
步骤8,根据所述方位向调制波形的调制序列波形,对M个距离单元各自对应的脉冲表达式分别进行方位向匹配滤波处理,最终得到匹配滤波后的机载雷达成像
步骤9,将所述匹配滤波后的机载雷达成像转化为大地坐标系下,即可得到前视场景下的机载雷达成像。
本发明的有益效果为:
第一,本发明创造性地提出了一种方位调制的SBAMIR成像框架,突破了传统的机载雷达成像对多普勒相位历史的依赖,不仅能够在多普勒相位历史变化很小的情况下实现前视场景下机载雷达的高分辨成像,还能够实现非前视场景下机载雷达的高分辨成像,使得本发明的方法的应用范围更加广泛;
第二,本发明方法在具体实现时不需要对现有的机载雷达进行硬件改进,只需通过发射波形的改变便能够在实际工程中进行实现应用。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是本发明的一种基于发射波形方位向调制的机载雷达前视成像方法流程示意图;
图2是本发明SBAMIR的成像示意图;其中,v表示机载雷达的载机平台速度;
图3是20dB条件下SBAMIR框架下的未进行多普勒中心补偿的前视成像结果;
图4是20dB条件下SBAMIR框架下的经过多普勒中心补偿后的前视成像结果;
图5是本发明方法在0dB条件下的前视成像结果图。
具体实施方式
参照图1,为本发明的一种基于发射波形方位向调制的机载雷达前视成像方法流程示意图,该种基于发射波形方位向调制的机载雷达前视成像方法,包括以下具体步骤:
步骤1,获取方位向调制波形,得到所述方位向调制波形的调制序列波形。
具体地,获取方位向调制波形,该组方位向调制波形可以是线性调频信号、非线性调频信号,也可以是编码信号。由于线性调频信号、非线性调频信号和编码信号在原理上具有一致性,所以此处设定该组方位向调制波形为编码信号,并得到编码信号的调制序列波形,记为w(tn),其表达式为:
w ( t n ) = Σ n = 0 N - 1 c n u ( t n - n τ )
其中,u(tn)表示第n时刻的二项码子脉冲包络,且u(tn)=rect(tn/τ),rect(·)表示矩形窗函数,τ表示二项码子脉冲的时间宽度,cn表示第n时刻的子脉冲幅度值,并且第n时刻的子脉冲幅度值为+1或-1;N表示编码信号的二项码调制序列长度,tn表示第n时刻。
步骤2,机载雷达在距离向发射线性调频信号,同时根据所述方位向调制波形的调制序列波形对所述线性调频信号在方位向进行缓慢调制,得到经过方位向缓慢调制后的发射信号波形其中,表示快时间,tn表示第n时刻。
具体地,机载雷达在距离向发射线性调频信号,同时根据所述方位向调制波形的调制序列波形w(tn)对所述线性调频信号在方位向进行缓慢调制,使得方位向的线性调频信号幅度按照步骤1得到的所述调制序列波形w(tn)变化。为了保证所述线性调频信号在距离向和方位向的同步,所述线性调频信号的方位调制信号的慢时间采样间隔,必须与所述线性调频信号的脉冲重复周期一致,进而得到经过方位向缓慢调制后的发射信号波形其表达式为:
s ( t ^ , t n ) = Σ n = 0 N - 1 c n · u ( t n - n τ ) · r e c t ( t ^ T p ) · exp [ j 2 π ( f c t ^ + γ 2 t ^ 2 ) ]
其中,γ表示发射的线性调频信号的调频率,表示快时间,tn表示第n时刻,fc表示发射的线性调频信号的载频,Tp表示每一个脉冲的持续时间,cn表示第n时刻的子脉冲幅度值,并且第n时刻的子脉冲幅度值为+1或-1;u(tn)表示第n时刻的二项码子脉冲包络,rect(·)表示矩形窗函数,n∈{1,2,…,N},N表示编码信号的二项码调制序列长度,n跟时刻有关。
设定机载雷达发射线性调频信号的脉冲重复频率为tr,为了保证机载雷达在距离向发射线性调频信号与方位向发射线性调频信号同步一致,此处设定二项码子脉冲的时间宽度τ与机载雷达发射线性调频信号的脉冲重复频率为tr相等,因此得到经过方位向缓慢调制后的发射信号波形另一个表达式
s ^ ( t ^ , t n ) = Σ n = 0 N - 1 c n · u ( t n - nt r ) · r e c t ( t ^ T p ) · exp [ j 2 π ( f c t ^ + γ 2 t ^ 2 ) ]
步骤3,机载雷达接收第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号,并对所述第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号依次进行距离向的脉冲压缩和匹配滤波处理,得到第q个相干积累时间间隔对应的距离脉压回波信号;其中,q∈{1,2,…,Q},Q表示机载雷达天线在一个周期内扫描的相干积累时间间隔总个数。
步骤4,对第q个相干积累时间间隔对应的距离脉压回波信号在距离频域内乘以距离走动校正因子H(f,t),得到第q个相干积累时间间隔对应的距离徙动校正回波信号;其中,f表示距离频率,t表示慢时间。
具体地,所述距离走动校正因子H(f,t),其表达式为:
H ( f , t ) = exp ( - j 4 π f R ( t ) c )
其中,R(t)≈R0-v×t×cosθ0,v表示机载雷达的载机平台速度,θ0表示机载雷达的波束视线与机载雷达的载机运动方向之间的夹角,R0表示机载雷达与散射点目标之间的初始斜距,f表示距离频率,t表示慢时间,c表示光速。
步骤5,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行相关处理并求其平均值,获得第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd
具体地,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行相关处理并求其平均值,并结合惯导信息获得第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd,其表达式为:
f d = a n g l e { E [ s ( m , k ) s * ( m , k - l ) ] } 2 π l f r + f r × f l o o r [ f d c _ I N S / f r ]
其中,s(m,k)表示第m个距离单元、第k个脉冲的回波信号离散形式,m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;fr表示机载雷达的脉冲重复频率,fdc_INS表示由惯导信息估计的多普勒中心频率,()*表示共轭,angle(·)表示求角度,floor[·]表示四舍五入操作,k∈{1,…,N},N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;为了方便,通常情况下取l=1。
步骤6,根据第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行多普勒中心补偿,得到第q个相干积累时间间隔对应的多普勒中心补偿回波信号。
具体地,根据第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd,对所述第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号进行多普勒中心补偿,即将第q个相干积累时间间隔对应的距离徙动校正回波信号乘以多普勒中心补偿函数,得到第q个相干积累时间间隔对应的多普勒中心补偿回波信号,所述多普勒中心补偿函数为exp(-j2πfdk/fr),其中,k∈{1,…,N},N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;fr表示机载雷达的脉冲重复频率。
步骤7,重复执行步骤3-步骤6,直到得到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,此时得到第1个相干积累时间间隔对应的多普勒中心补偿回波信号到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,并将其按照Q个相干积累时间间隔各自经过方位向缓慢调制后的发射信号波形的回波信号接收顺序进行存储,得到机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号,其中每一个相干积累时间间隔对应的回波信号包含N个脉冲,并将机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号划分为M个距离单元,然后计算得到第m个距离单元的脉冲表达式S(m);m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数。
具体地,第m个距离单元的脉冲表达式S(m)为:
S(m)=[s(m,q1),s(m,q2),…,s(m,qk),…,s(m,Nscan)]T
其中,s(m,qk)表示第m个距离单元中第q个相干积累时间间隔内的第k个脉冲,m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;qk∈{q1,q2,…,Nscan},qk=q×(N-1)+k,Nscan=Q×N,Q表示机载雷达天线在一个周期内扫描的相干积累时间间隔总个数;N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;Nscan表示机载雷达在一个天线扫描周期内扫描的脉冲总个数。
步骤8,根据所述方位向调制波形的调制序列波形,对M个距离单元各自对应的脉冲表达式分别进行方位向匹配滤波处理,最终得到匹配滤波后的机载雷达成像
具体地,所述匹配滤波后的成像的表达式为:
s ( t ^ , t n ) = As a ( t n - θ ( t n ) ω ) sin c [ B ( t ^ - 2 R ( t ) c ) ]
其中,sinc[·]表示辛格函数,R(tn)表示tn时刻机载雷达与目标的瞬时斜距,θ(tn)表示tn时刻机载雷达与目标的瞬时角度,tn表示第n时刻,ω表示机载雷达天线扫描的角速度,B表示机载雷达发射的雷达信号带宽,sa(t)表示二项码序列的相关函数,A表示机载雷达所在场景的散射点系数。
步骤9,将所述匹配滤波后的机载雷达成像转化为大地坐标系下,即可得到前视场景下的机载雷达成像。
本发明的优点可通过以下仿真实验作进一步说明。
1)仿真实验环境:
实验环境:MATLABR2009b,Intel(R)Pentium(R)2CPU2.7GHz,Window7旗舰版。
2)仿真实验内容及结果分析:
在仿真实验环境下,分别应用本发明的SBAMIR成像方法进行成像,前视场景是一个包含了107散射点的复杂场景,该复杂场景形成了机载雷达的飞机轮廓;并且该复杂场景中有许多散射点是非常靠近的,其中最近的相邻散射点之间只有0.05°;图3是20dB条件下SBAMIR框架下的未进行多普勒中心补偿的前视成像结果;图4是20dB条件下SBAMIR框架下的经过多普勒中心补偿后的前视成像结果;图5是本发明方法在0dB条件下的前视成像结果图。其中,成像角度范围分别是机载雷达的载机正前方正负5°的范围。
从图3中可以看出,由于受多普勒频移的影响,直接的方位匹配滤波的前视成像结果具有一定的模糊性。
从图4中可以看出,经过多普勒相位补偿之后,本发明方法的前视成像结果从机载雷达的载机轮廓的细节信息能够非常明显的显示出来。
从图5中可以看出,即使是在0dB情况下,本发明的SBAMIR方法依旧可以对场景中的飞机轮廓进行很好的成像。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (9)

1.一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,包括以下步骤:
步骤1,获取方位向调制波形,得到所述方位向调制波形的调制序列波形;
步骤2,机载雷达在距离向发射线性调频信号,同时根据所述方位向调制波形的调制序列波形对所述线性调频信号在方位向进行缓慢调制,得到经过方位向缓慢调制后的发射信号波形其中,表示快时间,tn表示第n时刻;
步骤3,机载雷达接收第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号,并对所述第q个相干积累时间间隔内经过方位向缓慢调制后的发射信号波形的回波信号依次进行距离向的脉冲压缩和匹配滤波处理,得到第q个相干积累时间间隔对应的距离脉压回波信号;其中,q∈{1,2,…,Q},Q表示机载雷达天线在一个周期内扫描的相干积累时间间隔总个数;
步骤4,对第q个相干积累时间间隔对应的距离脉压回波信号在距离频域内乘以距离走动校正因子H(f,t),得到第q个相干积累时间间隔对应的距离徙动校正回波信号;其中,f表示距离频率,t表示慢时间;
步骤5,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行相关处理并求其平均值,获得第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd
步骤6,根据第q个相干积累时间间隔对应的机载雷达多普勒中心估计值fd,对第q个相干积累时间间隔对应的距离徙动校正回波信号进行多普勒中心补偿,得到第q个相干积累时间间隔对应的多普勒中心补偿回波信号;
步骤7,重复执行步骤3-步骤6,直到得到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,此时得到第1个相干积累时间间隔对应的多普勒中心补偿回波信号到第Q个相干积累时间间隔对应的多普勒中心补偿回波信号,并将其按照Q个相干积累时间间隔各自经过方位向缓慢调制后的发射信号波形的回波信号接收顺序进行存储,得到机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号,其中每一个相干积累时间间隔对应的回波信号包含N个脉冲,并将机载雷达在一个扫描周期Q个相干积累时间间隔内的回波信号划分为M个距离单元,然后计算得到第m个距离单元的脉冲表达式S(m);m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;
步骤8,根据所述方位向调制波形的调制序列波形,对M个距离单元各自对应的脉冲表达式分别进行方位向匹配滤波处理,最终得到匹配滤波后的机载雷达成像
步骤9,将所述匹配滤波后的机载雷达成像转化为大地坐标系下,即可得到前视场景下的机载雷达成像。
2.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤1中,所述方位向调制波形还包括:所述方位向调制波形为线性调频信号、非线性调频信号或编码信号。
3.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤1中,所述方位向调制波形的调制序列波形,记为w(tn),其表达式为:
w ( t n ) = Σ n = 0 N - 1 c n u ( t n - n τ )
其中,u(tn)表示第n时刻的二项码子脉冲包络,且u(tn)=rect(tn/τ),rect(·)表示矩形窗函数,τ表示二项码子脉冲的时间宽度,cn表示第n时刻的子脉冲幅度值,N表示编码信号的二项码调制序列长度,tn表示第n时刻。
4.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤2中,所述经过方位向缓慢调制后的发射信号波形其表达式为:
s ( t ^ , t n ) = Σ n = 0 N - 1 c n · u ( t n - n τ ) · r e c t ( t ^ T p ) · exp [ j 2 π ( f c t ^ + γ 2 t ^ 2 ) ]
其中,γ表示发射的线性调频信号的调频率,表示快时间,tn表示第n时刻,fc表示发射的线性调频信号的载频,Tp表示每一个脉冲的持续时间,cn表示第n时刻的子脉冲幅度值,u(tn)表示第n时刻的二项码子脉冲包络,rect(·)表示矩形窗函数,n∈{1,2,…,N},N表示编码信号的二项码调制序列长度。
5.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤4中,所述距离走动校正因子H(f,t),其表达式为:
H ( f , t ) = exp ( - j 4 π f R ( t ) c )
其中,R(t)≈R0-v×t×cosθ0,v表示机载雷达的载机平台速度,θ0表示机载雷达的波束视线与机载雷达的载机运动方向之间的夹角,R0表示机载雷达与散射点目标之间的初始斜距,f表示距离频率,t表示慢时间,c表示光速。
6.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤5中,所述机载雷达的多普勒中心估计值fd,其表达式为:
f d = a n g l e { E [ s ( m , k ) s * ( m , k - l ) ] } 2 π l f r + f r × f l o o r [ f d c _ I N S / f r ]
其中,s(m,k)表示第m个距离单元、第k个脉冲的回波信号离散形式,m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;fr表示机载雷达的脉冲重复频率,fdc_INS表示由惯导信息估计的多普勒中心频率,()*表示共轭,angle(·)表示求角度,floor[·]表示四舍五入操作,k∈{1,…,N},N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;l表示设定的正整数。
7.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤6中,所述对一个周期内机载雷达接收到的所述线性调频信号的回波数据进行多普勒中心补偿,具体为:
将距离徙动校正后的数据乘以多普勒中心补偿函数,得到经过多普勒中心补偿后的回波数据,所述多普勒中心补偿函数为exp(-j2πfdk/fr),其中,k∈{1,…,N},N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;fr表示机载雷达的脉冲重复频率。
8.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤7中,所述第m个距离单元的脉冲表达式S(m)为:
S(m)=[s(m,q1),s(m,q2),…,s(m,qk),…,s(m,Nscan)]T
其中,s(m,qk)表示第m个距离单元中第q个相干积累时间间隔内的第k个脉冲,m∈{1,2,…,M},M表示机载雷达在一个扫描周期内接收到的回波信号包含的距离单元总个数;qk∈{q1,q2,…,Nscan},qk=q×(N-1)+k,Nscan=Q×N,Q表示机载雷达天线在一个周期内扫描的相干积累时间间隔总个数;N表示机载雷达在一个相干积累时间间隔内的脉冲总个数,N也表示编码信号的二项码调制序列长度;Nscan表示机载雷达在一个天线扫描周期内扫描的脉冲总个数。
9.如权利要求1所述的一种基于发射波形方位向调制的机载雷达前视成像方法,其特征在于,在步骤8中,所述匹配滤波后的成像的表达式为:
s ( t ^ , t n ) = As a ( t n - θ ( t n ) ω ) sin c [ B ( t ^ - 2 R ( t ) c ) ]
其中,sinc[·]表示辛格函数,R(tn)表示tn时刻机载雷达与目标的瞬时斜距,θ(tn)表示tn时刻机载雷达与目标的瞬时角度,tn表示第n时刻,ω表示机载雷达天线扫描的角速度,B表示机载雷达发射的雷达信号带宽,sa(t)表示二项码序列的相关函数,A表示机载雷达所在场景的散射点系数。
CN201610065516.3A 2016-01-29 2016-01-29 一种基于发射波形方位向调制的机载雷达前视成像方法 Active CN105717508B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610065516.3A CN105717508B (zh) 2016-01-29 2016-01-29 一种基于发射波形方位向调制的机载雷达前视成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610065516.3A CN105717508B (zh) 2016-01-29 2016-01-29 一种基于发射波形方位向调制的机载雷达前视成像方法

Publications (2)

Publication Number Publication Date
CN105717508A true CN105717508A (zh) 2016-06-29
CN105717508B CN105717508B (zh) 2019-02-19

Family

ID=56155372

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610065516.3A Active CN105717508B (zh) 2016-01-29 2016-01-29 一种基于发射波形方位向调制的机载雷达前视成像方法

Country Status (1)

Country Link
CN (1) CN105717508B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106707278A (zh) * 2017-01-18 2017-05-24 北京无线电测量研究所 一种基于稀疏表示的多普勒波束锐化成像方法及装置
CN108572362A (zh) * 2018-01-30 2018-09-25 哈尔滨工业大学 一种tws雷达空时联合关联跟踪方法及装置
CN110568410A (zh) * 2019-10-09 2019-12-13 上海无线电设备研究所 一种空间频率色散的微波雷达超分辨方法
CN111381217A (zh) * 2020-04-01 2020-07-07 上海无线电设备研究所 基于低精度惯性导航系统的弹载sar运动补偿方法
WO2021175168A1 (zh) * 2020-03-02 2021-09-10 华为技术有限公司 信号处理方法及装置
CN113820673A (zh) * 2021-10-19 2021-12-21 中安锐达(北京)电子科技有限公司 一种基于伺服转速的机扫雷达方位波束中心的修正方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008134224A (ja) * 2006-08-31 2008-06-12 Sierra Nevada Corp 三次元レーダー像のレンダリングのための装置及び方法
CN102331577A (zh) * 2011-09-21 2012-01-25 北京理工大学 一种适用于地球同步轨道sar的改进ncs成像算法
CN103308913A (zh) * 2013-07-12 2013-09-18 西安电子科技大学 一种高速飞行器载双天线前视sar解模糊算法
CN103885058A (zh) * 2014-02-26 2014-06-25 西安电子科技大学 一种利用稀疏统计特性的机载雷达前视超分辨成像方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008134224A (ja) * 2006-08-31 2008-06-12 Sierra Nevada Corp 三次元レーダー像のレンダリングのための装置及び方法
CN102331577A (zh) * 2011-09-21 2012-01-25 北京理工大学 一种适用于地球同步轨道sar的改进ncs成像算法
CN103308913A (zh) * 2013-07-12 2013-09-18 西安电子科技大学 一种高速飞行器载双天线前视sar解模糊算法
CN103885058A (zh) * 2014-02-26 2014-06-25 西安电子科技大学 一种利用稀疏统计特性的机载雷达前视超分辨成像方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HONG MENG CHEN等: "Sparse super-resolution imaging for airborne single channel forward-looking radar in expanded beam space via Ip regularisation", 《ELECTRNICS LETTERS》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106707278A (zh) * 2017-01-18 2017-05-24 北京无线电测量研究所 一种基于稀疏表示的多普勒波束锐化成像方法及装置
CN106707278B (zh) * 2017-01-18 2019-12-06 北京无线电测量研究所 一种基于稀疏表示的多普勒波束锐化成像方法及装置
CN108572362A (zh) * 2018-01-30 2018-09-25 哈尔滨工业大学 一种tws雷达空时联合关联跟踪方法及装置
CN108572362B (zh) * 2018-01-30 2020-07-31 哈尔滨工业大学 一种tws雷达空时联合关联跟踪方法及装置
CN110568410A (zh) * 2019-10-09 2019-12-13 上海无线电设备研究所 一种空间频率色散的微波雷达超分辨方法
CN110568410B (zh) * 2019-10-09 2021-08-31 上海无线电设备研究所 一种空间频率色散的微波雷达超分辨方法
WO2021175168A1 (zh) * 2020-03-02 2021-09-10 华为技术有限公司 信号处理方法及装置
CN111381217A (zh) * 2020-04-01 2020-07-07 上海无线电设备研究所 基于低精度惯性导航系统的弹载sar运动补偿方法
CN113820673A (zh) * 2021-10-19 2021-12-21 中安锐达(北京)电子科技有限公司 一种基于伺服转速的机扫雷达方位波束中心的修正方法

Also Published As

Publication number Publication date
CN105717508B (zh) 2019-02-19

Similar Documents

Publication Publication Date Title
CN105717508A (zh) 一种基于发射波形方位向调制的机载雷达前视成像方法
CN106970386B (zh) 一种雷达多普勒波束锐化的优化方法
US6911933B1 (en) Dynamic logic algorithm used for detecting slow-moving or concealed targets in synthetic aperture radar (SAR) images
CN106872974B (zh) 基于高超声速平台双通道雷达的高精度运动目标成像方法
NO164623B (no) Fremgangsmaate og system for ved luftbaaren radar aa danne et avstand/asimut-bilde av et maal-skip.
CN102707269B (zh) 一种机载雷达距离走动校正方法
CN106597437B (zh) 高频机载调频连续波sar运动补偿及成像处理方法
CN103983974B (zh) 双站调频连续波合成孔径雷达成像方法
CN102778681B (zh) 固定发射站双基地前视合成孔径雷达成像方法
CN110488230B (zh) 一种双多基前视成像方法和装置
CN104330779B (zh) 机载合成孔径雷达运动误差补偿方法
CN105607055A (zh) 一种基于天线方向图的机载雷达单脉冲前视成像方法
CN104865573A (zh) 一种双基前视高机动平台sar扩展场景成像优化方法
CN107765226A (zh) 一种sar卫星雷达回波模拟方法、系统和介质
CN113589285B (zh) 一种飞行器sar实时成像方法
CN105182340A (zh) 一种双基地前视sar运动补偿方法
CN105487074A (zh) 一种双基地合成孔径雷达数值距离多普勒成像方法
CN109001700B (zh) 一种实现目标轮廓重构的雷达前视成像方法
Wang et al. Detection and tracking of UAVs using interferometric radar
CN105022060A (zh) 针对快速空天目标的步进isar成像方法
CN102176013B (zh) 一种弹载非理想凝视聚束式sar混合度提取方法
CN106526544B (zh) 基于高超声速平台的mimosar杂波抑制方法
CN103235308B (zh) 雷达正前视扫描相干成像方法
CN104360333A (zh) 一种同时校正一阶和二阶距离徙动的相参积累检测方法
CN105223573A (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
GR01 Patent grant
GR01 Patent grant