发明名称 基于自适应误差剔除的逆合成孔径雷达包络对齐方法
摘要 本发明公开了一种基于自适应误差剔除的逆合成孔径雷达包络对齐方法。主要解决现有技术运算量大和对齐精度低之问题,其实现过程为:由相邻相关法得相邻走动量估计值和绝对走动量估计值;对绝对走动量估计值的最小二乘拟合值求一阶差分得相邻走动量拟合值,由拟合值与估计值之差得到相邻走动量拟合误差并计算其标准差,定义大于3倍标准差的拟合误差为突跳误差,并将对应的相邻走动量估计值用拟合值替换,对其累加得新绝对走动量估计值,重复进行拟合、替换和更新至无突跳误差检出;根据新绝对走动量估计值的最小二乘拟合值完成包络对齐。本发明具有运算量小和包络对齐精度高之优点,用于在低信噪比条件下估计并补偿逆合成孔径雷达回波的平动分量。
申请公布号 CN102262222B 申请公布日期 2013.11.06
申请号 CN201110155651.4 申请日期 2011.06.10
申请人 西安电子科技大学 发明人 廖桂生;杨志伟;刘志凌;李延;曾操;徐青;束宇翔;何嘉懿
分类号 G01S7/02(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/02(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于自适应误差剔除的逆合成孔径雷达包络对齐方法,包括如下步骤: (1)用相邻相关包络对齐算法依次计算相邻雷达回波间走动量的相邻回波走动量估计值δ: 1a)根据相邻相关包络对齐算法,计算δ的第k个数据值δ<sub>k(k+1)</sub>: <img file="FDA00003191191400011.GIF" wi="783" he="118" />式中s<sub>k</sub>(t)和s<sub>k+1</sub>(t)分别为第k次回波信号和第k+1次回波信号,k=1,...,N-1,N为观测时间内的回波脉冲数; 1b)根据数据值δ<sub>k(k+1)</sub>,k=1,...,N-1,得到相邻回波走动量估计值δ表达式为: δ=[δ<sub>12</sub>,δ<sub>23</sub>...,δ<sub>(N-1)N</sub>]<sup>T</sup>, 其中上标T表示转置操作; (2)累加δ中的数据得到各次回波的绝对走动量估计值为D,并初始化误差剔除后新的绝对走动量估计值D'=D; (3)采用最小二乘多项式拟合绝对走动量估计值D,得到相应拟合值<img file="FDA00003191191400012.GIF" wi="117" he="93" />3a)设t<sub>m</sub>=[t<sub>1</sub>,...,t<sub>N</sub>]<sup>T</sup>为慢时间向量,其中t<sub>k</sub>=k·PRI为雷达发射脉冲对应的慢时间,k=1,...,N,N为观测时间内的回波脉冲数,PRI为脉冲重复周期,上标T表示转置操作;设<img file="FDA00003191191400013.GIF" wi="452" he="79" />为相应的慢时间矩阵,其中M为拟合多项式阶数,<img file="FDA00003191191400014.GIF" wi="342" he="86" />l=1,...,M,e=[1,...,1]<sup>T</sup>为单位矢量;设<img file="FDA00003191191400015.GIF" wi="504" he="94" />为拟合多项式系数向量的最小二乘估计值,其中<img file="FDA00003191191400016.GIF" wi="50" he="89" />对应<img file="FDA00003191191400017.GIF" wi="56" he="84" />的系数,l=1,...,M,<img file="FDA00003191191400018.GIF" wi="56" he="93" />对应单位矢量e的系数;3b)根据最小二乘理论,利用如下公式计算拟合多项式系数向量的最小二乘估计值<img file="FDA00003191191400019.GIF" wi="109" he="92" /><img file="FDA000031911914000110.GIF" wi="388" he="88" />进而计算绝对走动量估计值D的拟合值<img file="FDA00003191191400021.GIF" wi="90" he="93" />为:<img file="FDA00003191191400022.GIF" wi="303" he="93" />(4)对拟合值<img file="FDA00003191191400023.GIF" wi="86" he="87" />求一阶差分,得到对应的雷达相邻回波走动量拟合值<img file="FDA00003191191400024.GIF" wi="68" he="74" />计算相邻回波走动量估计值δ与拟合值<img file="FDA00003191191400025.GIF" wi="32" he="73" />之差可得拟合误差Δδ:(5)设δ<sub>k(k+1)</sub>、<img file="FDA00003191191400026.GIF" wi="127" he="93" />和Δδ<sub>k(k+1)</sub>分别为所述δ、<img file="FDA00003191191400027.GIF" wi="42" he="77" />和Δδ中的第k个数据值,其中k=1,...,N-1,N为观测时间内的回波脉冲数;根据标准差公式:<img file="FDA00003191191400028.GIF" wi="742" he="170" />计算拟合误差Δδ中数据的标准差<img file="1.GIF" wi="60" he="31" />,若<img file="2.GIF" wi="153" he="39" />,)则认为δ<sub>k(k+1)</sub>存在突跳误差,用<img file="FDA00003191191400029.GIF" wi="127" he="96" />替换δ<sub>k(k+1)</sub>以实现误差剔除,并对剔除所有突跳误差后的相邻回波走动量估计值δ'中数据进行累加,得到新的绝对走动量估计值D',式中Δδ<sub>i(i+1)</sub>为所述拟合误差Δδ中的第i个数据值,其中i=1,...,N-1,N为观测时间内的回波脉冲数,<img file="FDA000031911914000210.GIF" wi="434" he="149" />为拟合误差Δδ中数据的均值;令D=D'并重复执行步骤(2)至步骤(5),直至无突跳误差检出; (6)对误差剔除后新的绝对走动量估计值D'进行最小二乘多项式拟合,得到相应拟合值<img file="FDA000031911914000211.GIF" wi="118" he="88" />6a)根据最小二乘理论,利用如下公式计算与该D'对应的拟合多项式系数向量的最小二乘估计值<img file="FDA000031911914000212.GIF" wi="104" he="94" /><img file="FDA000031911914000213.GIF" wi="419" he="92" />其中上标T表示转置操作; 6b)计算误差剔除后新的绝对走动量估计值D'的拟合值<img file="FDA000031911914000214.GIF" wi="88" he="87" />为:<img file="FDA000031911914000215.GIF" wi="277" he="93" />其中C为慢时间矩阵; (7)根据拟合值<img file="FDA000031911914000216.GIF" wi="88" he="87" />对目标平动分量进行高精度补偿,完成包络对齐。
地址 710071 陕西省西安市太白南路2号