主权项 |
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" />对目标平动分量进行高精度补偿,完成包络对齐。 |