主权项 |
1.一种弹载SAR子孔径前斜视高阶非线性调频变标成像方法,其弹载SAR工作在条带模式,弹载平台高度为H,以速度v沿X轴匀速直线飞行,θ<sub>0</sub>为波束射线指向的斜视角,R<sub>0</sub>为波束中心线扫过目标时的斜距,t<sub>m</sub>为方位慢时间,以弹载平台位于O点的时刻为方位慢时间的起点,R<sub>B</sub>为场景中某一点目标与弹载SAR的最近距离;点目标到弹载SAR的瞬时斜距为<maths num="0001"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mi>cos</mi><msub><mi>θ</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mi>sin</mi><msub><mi>θ</mi><mn>0</mn></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>vt</mi><mi>m</mi></msub><mo>-</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths>其中X<sub>n</sub>=vt<sub>n</sub>,t<sub>n</sub>为天线波束中心穿越目标点的时刻;其特征在于,包括以下步骤:步骤1,弹载SAR发射的是线性调频信号,则其基频回波信号为S<sub>0</sub>(t<sub>r</sub>,t<sub>m</sub>),其中t<sub>r</sub>为距离快时间,t<sub>m</sub>为方位慢时间;步骤2,对基频回波信号S<sub>0</sub>(t<sub>r</sub>,t<sub>m</sub>)进行距离向傅里叶变换处理,得到距离频域方位时域信号S<sub>1</sub>(f<sub>r</sub>,t<sub>m</sub>),其中f<sub>r</sub>为距离频率;步骤3,对距离频域方位时域信号S<sub>1</sub>(f<sub>r</sub>,t<sub>m</sub>)进行距离向脉冲压缩处理,得到距离脉压信号S<sub>2</sub>(f<sub>r</sub>,t<sub>m</sub>);步骤4,对距离脉压信号S<sub>2</sub>(f<sub>r</sub>,t<sub>m</sub>)进行线性距离走动校正,得到线性距离走动校正后的信号S<sub>3</sub>(f<sub>r</sub>,t<sub>m</sub>);步骤5,采用驻定相位原理,对线性距离走动校正后的信号S<sub>3</sub>(f<sub>r</sub>,t<sub>m</sub>)作方位向傅里叶变换,得到二维频域信号S<sub>4</sub>(f<sub>r</sub>,f<sub>a</sub>),其中f<sub>a</sub>为方位频率;步骤6,对二维频域信号S<sub>4</sub>(f<sub>r</sub>,f<sub>a</sub>)进行距离弯曲校正和二次距离压缩,得到解耦后的二维频域信号S<sub>5</sub>(f<sub>r</sub>,f<sub>a</sub>),然后对解耦后的二维频域信号S<sub>5</sub>(f<sub>r</sub>,f<sub>a</sub>)进行距离向逆傅里叶变换,得到距离时域方位频域信号S<sub>6</sub>(t<sub>r</sub>,f<sub>a</sub>);步骤7,对距离时域方位频域信号S<sub>6</sub>(t<sub>r</sub>,f<sub>a</sub>)在f<sub>a</sub>=0处进行方位向五阶泰勒级数展开,并对方位向四阶、五阶相位进行补偿,得到补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>);步骤8,对补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>)作方位向逆傅里叶变换,得到补偿完方位向四阶、五阶相位的时域信号S<sub>10</sub>(t<sub>r</sub>,t<sub>m</sub>),<img file="FDA0000483718310000021.GIF" wi="1481" he="202" />其中函数<img file="FDA0000483718310000022.GIF" wi="422" he="140" />B<sub>r</sub>为发射信号的频带,w<sub>a</sub>(·)为方位窗函数,t<sub>r</sub>为距离快时间,t<sub>m</sub>为方位慢时间,c为光速,exp(·)为指数函数,t<sub>n</sub>为天线波束中心穿越目标点的时刻,R<sub>0</sub>为波束中心线扫过目标时的斜距,<img file="FDA0000483718310000025.GIF" wi="57" he="48" />为常数项,f<sub>DR</sub>为多普勒调频率,f<sub>DT</sub>为方位向三次项系数;对f<sub>DR</sub>和f<sub>DT</sub>分别作泰勒级数展开得<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>f</mi><mi>DR</mi></msub><mo>≈</mo><msub><mi>f</mi><mi>DRref</mi></msub><mo>+</mo><msub><mi>k</mi><mi>SAC</mi></msub><msub><mi>t</mi><mi>n</mi></msub><mo>+</mo><msub><mi>k</mi><mi>TAC</mi></msub><msubsup><mi>t</mi><mi>n</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>DT</mi></msub><mo>≈</mo><msub><mi>f</mi><mi>DTref</mi></msub><mo>+</mo><msub><mi>k</mi><mi>SAT</mi></msub><msub><mi>t</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>其中f<sub>DRref</sub>为多普勒调频率的参考值,k<sub>SAC</sub>t<sub>n</sub>为多普勒调频率的一阶空变项,<img file="FDA0000483718310000026.GIF" wi="119" he="66" />为多普勒调频率的二阶空变项,f<sub>DTref</sub>为方位三次相位的参考值,k<sub>SAT</sub>t<sub>n</sub>为方位三次相位的一阶空变项;提取为多普勒调频率的一阶空变项k<sub>SAC</sub>t<sub>n</sub>、多普勒调频率的二阶空变项<img file="FDA0000483718310000027.GIF" wi="144" he="67" />为方位三次相位的一阶空变项k<sub>SAT</sub>t<sub>n</sub>;步骤9,对补偿完方位向四阶、五阶相位的距离时域方位频域信号S<sub>9</sub>(t<sub>r</sub>,f<sub>a</sub>)引入用于消除多普勒调频率的一阶空变项k<sub>SAC</sub>t<sub>n</sub>、多普勒调频率的二阶空变项<img file="FDA0000483718310000028.GIF" wi="119" he="67" />和方位三次相位的一阶空变项k<sub>SAT</sub>t<sub>n</sub>的方位高阶非线性变标扰动因子H<sub>ncs</sub>以消除调频率和相位的误差,再对引入方位高阶非线性变标扰动因子H<sub>ncs</sub>后的回波信号进行方位向逆傅里叶变换,得到消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>);步骤10,消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)中含有因引入方位高阶非线性调频变标扰动因子H<sub>ncs</sub>产生的高次相位,对消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)补偿剩余高次相位,其方位剩余高次相位补偿因子H<sub>amz</sub>为<maths num="0003"><![CDATA[<math><mrow><msub><mi>H</mi><mi>amz</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mi>jπ</mi><mfrac><msub><mi>K</mi><mi>SAC</mi></msub><mn>3</mn></mfrac><msubsup><mi>t</mi><mi>m</mi><mn>3</mn></msubsup><mo>-</mo><mi>jπ</mi><mfrac><msub><mi>K</mi><mi>TAC</mi></msub><mn>6</mn></mfrac><msubsup><mi>t</mi><mi>m</mi><mn>4</mn></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>将消除方位空变后的时域信号S<sub>11</sub>(t<sub>r</sub>,t<sub>m</sub>)与方位剩余高次相位补偿因子H<sub>amz</sub>相乘,得到补偿完方位剩余高次相位的时域信号S<sub>12</sub>(t<sub>r</sub>,t<sub>m</sub>);步骤11,对补偿完方位剩余高次相位的时域信号S<sub>12</sub>(t<sub>r</sub>,t<sub>m</sub>)进行方位向的去斜处理,然后作方位向傅里叶变换,得到方位频域聚焦信号S<sub>13</sub>(t<sub>r</sub>,f<sub>a</sub>),从而将图像聚焦在方位频域上。 |