基于毫米波雷达的近场ISAR干涉转台成像方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112764031A(43)申请公布日2021.05.07(21)申请号202011572206.3(22)申请日2020.12.28(71)申请人中国海洋大学地址266100山东省青岛市崂山区松岭路238号(72)发明人王超张彦敏王运华(74)专利代理机构青岛海昊知识产权事务所有限公司37201代理人张中南邱岳(51)Int.Cl.G01S13/90(2006.01)权利要求书3页说明书7页附图7页(54)发明名称基于毫米波雷达的近场ISAR干涉转台成像方法(57)摘要基于毫米波雷达的ISAR干涉转台成像算法，利用毫米波雷达的一发双收模式对近场转台目标干涉成像，其中发射天线A发射调频连续波信号，接收天线B和C接收信号，分别对接收天线B和接收天线C得到的回波数据进行脉压处理得到两幅目标回波的一维距离像，对接收天线B的回波数据计算得到的一维距离像进行距离向插值、相位补偿最后相干叠加得到BP成像结果得到目标的二维图像，对得到两幅一维距离像进行方位向压缩得到两幅距离多普勒图像，在将两幅图像进行干涉处理得到干涉相位，反演目标的第三维高度信息。本发明实现了近场条件下反演目标的三维尺寸，扩展了算法的应用场景，能够有效的提高对目标的识别能力。CN112764031ACN112764031A权利要求书1/3页1.基于毫米波雷达的近场ISAR干涉转台成像方法，其特征在于包括以下步骤：步骤1、利用正对目标放置的毫米波雷达发射天线A和接收天线B、C发射并接收目标的回波信号；所述毫米波雷达发射调频连续波信号，雷达发射信号的总帧数为M，每帧包括N个调频连续波信号，脉冲信号重复时间为T，第n个脉冲时刻时间为tn＝nT；步骤2、对天线B、C接收到的目标回波信号做脉冲压缩处理，得到目标的一维距离像；雷达得到的目标回波信号是由多个散射点组成的，对于目标回波信号上的每一个散射点P，以RAP(tn)、RBP(tn)和RCP(tn)分别表示散射点P到发射天线A、接收天线B和接收天线C的在tn时刻的瞬时距离；定义RB(tn)与RC(tn)分别如下：RB(tn)＝(RAP(tn)+RBP(tn))/2RC(tn)＝(RAP(tn)+RCP(tn))/2接收天线B和接收天线C接收到的回波信号SB(tk,tn)、SC(tk,tn)分别为：其中，tk为距离向快时间，tn为方位向慢时间、即步骤1中的第n个脉冲时刻时间，σp为散射点P的后向散射系数，c为真空中电磁波的传播速度，Tp发射信号的脉冲宽度，fc为调频连续波信号的载波频率，K为调频连续波信号的调频率；对两接收天线接收到的回波信号进行距离向脉冲压缩，得到一维距离像SBrp(frp,tn)、SCrp(frp,tn)分别为：其中，frp为距离向脉冲压缩后得到的距离频率，其与距离向中的距离相互对应，σpf为脉冲压缩后的后向散射系数；从而将由单个散射点的一维距离向扩展得到整个目标的一维距离像；步骤3、计算步骤2中得到的一维距离像的距离向分辨率和方位向分辨率；根据预先设置的调频连续波信号以及毫米波雷达的系统参数计算距离分辨率dr：2CN112764031A权利要求书2/3页其中，Fs为采样频率，As为采样点数；利用毫米波雷达系统的转台参数计算方位向分辨率da：其中ω为转台的旋转角速度，λ为发射信号电磁波的波长；步骤4、对接收天线B处理得到的目标一维距离像数据进行沿距离向插值，然后进行相位补偿，相干叠加得到BP成像结果：根据步骤2得到的一维距离像数据截取包含目标的成像区域，并根据此成像区域划分成像网格(x,y)，其中x为划分后的成像网格的方位向位置，y为划分后的成像区域的距离向位置，且划分后的成像网格区域包含目标区域，划分后的成像网格的采样间隔均小于步骤3中计算得到的距离向分辨率以及方位向分辨率；对截取的、包含目标的一维距离像数据沿距离向使用线性插值法进行插值升采样；记插值后的一维距离像数据起始时刻对应的方位向角度为0°，以此计算转台目标在tn时刻的旋转累加角度以转台旋转中心为参考系坐标原点，计算划分后的成像区域所有网格点(x,y)到转台旋转中心的距离Rxy：以转台旋转中心为参考系坐标原点，计算划分后的成像区域所有网格点(x,y)在tn时刻的位置坐标：得到tn时刻，划分后的成像区域所有网格点(x,y)到接收天线B的瞬时距离其中，为接收天线B到转台旋转中心的距离，为固定参数；对截取的、包含包含目标的一维距离像数据进行相位补偿，相干叠加得到BP成像结果；步骤5、将从接收天线B和接收天线C分别处理得到的一维距离像数据，进行方位向处理，分别得到接收天线B和接收天线C的距离多普勒成像结果，再将成像结果的相位分别表示为：3CN112764031A权利要求书3/3页其中t