发明名称 一种提高地面运动目标检测性能的多基线设计方法
摘要 本发明针对多通道单一SAR平台,分析基线长度对于动目标检测的影响,提出一种多基线的天线设计方法,并推导出相应的多通道ATI通用算法。本发明的基线设计方法增加系统选择基线的灵活性,可以根据不同需求,提高动目标检测性能,为后续处理提供丰富的干涉相位信息;与之匹配的ATI通用算法中适用于所有天线模式为一发多收的多通道GMTI系统,且涉及到的相位滤波器形式简单且有效,具有更好的通用性和更强的实用性。
申请公布号 CN102955158A 申请公布日期 2013.03.06
申请号 CN201210181199.3 申请日期 2012.06.04
申请人 北京航空航天大学 发明人 李景文;王冠琼;于泽;杨威;曾虹程
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 代理人
主权项 1.一种提高地面运动目标检测性能的多基线设计方法包括:一、设计天线构型,实现多基线要求;二、根据所设计天线构型,设计与之匹配的多通道ATI通用算法。具体内容如下:一、设计天线构型,实现多基线要求步骤一:(图1)设计天线为一发多收模式,全孔径发射,n个子孔径接收。称作通道i,i=0,±1,±2...。通道间的基线长度设置的原则是有长有短,并非均匀分配,且存在互质的基线。该设计有利于动目标检测性能的提高(延长盲速周期,扩大非模糊速度),以及多通道ATI通用算法的推导(将会在下一个内容中有所体现)。步骤二:设置通道0作为参考通道,其作用在推导多通道ATI通用算法时得以体现。步骤三:计算任意两个通道i和通道k相对于通道0的基线长度(即相位中心间距)为d<sub>i0</sub>,d<sub>k0</sub>。步骤四:计算通道i相对与通道k的基线长度d<sub>ik</sub>=d<sub>i0</sub>-d<sub>k0</sub>。步骤五:选择两个互质的基线d<sub>ik</sub>,d<sub>jr</sub>组成基线对。可以根据具体需求:(a)检测慢速目标或对测速精度要求较高时,采用长基线对的检测结果;(b)检测速度较快目标或对杂波抑制要求较高时,采用短基线对的检测结果;(c)当长短基线对的数据都可以采纳时,可根据实际情况或检测结果加权各基线的数据,得到较为理想的检测结果。二、根据第一部分天线构型,设计与之匹配的多通道ATI通用算法,具体内容为:步骤一:读入SAR系统的相关参数,包括:雷达高度H,雷达速度V<sub>a</sub>,雷达中心频率f<sub>0</sub>,光速c,雷达斜视角<img file="FDA00001722731500021.GIF" wi="54" he="37" />脉冲时宽T<sub>r</sub>,波束中心距离R<sub>c</sub>,距离向采样率F<sub>r</sub>,距离向调频率K<sub>r</sub>,天线方位向长度D,基线单位长度B。步骤二:根据天线配置方式、雷达与地面目标几何关系,推导SAR回波模型,具体为:(a)天线的配置方式为:全孔径发射信号,各子孔径同时接收信号。(b)设置通道0的相位中心设在天线中心(为便于推导所设,并非一定要位于天线中心)。天线延方位向排列,各通道相位中心用O<sub>i</sub>表示,获取任一通道i和通道0的基线长度d<sub>i0</sub>。(c)计算雷达平台到场景中心的最短距离为R<sub>b</sub>,获取雷达方位向速度V<sub>a</sub>。零时刻,天线的中间O<sub>0</sub>位于原点;目标位于P(x<sub>0</sub>,y<sub>0</sub>,h),方位向速度和对地(对斜面)距离向速度分别为v<sub>x</sub>,v<sub>y</sub>(v<sub>r</sub>)。(d)t<sub>a</sub>=t时刻,天线运动到O′<sub>i</sub>点,目标运动到P<sub>t</sub>点,计算O′<sub>i</sub>与P<sub>t</sub>距离R<sub>i</sub>(t):<maths num="0001"><![CDATA[<math><mrow><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>b</mi></msub><mo>+</mo><msub><mi>v</mi><mi>r</mi></msub><mi>t</mi><mo>+</mo><mfrac><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>d</mi></mrow><mrow><mi>i</mi><mn>0</mn></mrow></msub><mrow><mo>+</mo><mi>t</mi></mrow><mrow><mo>(</mo><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>V</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>v</mi><mi>r</mi></msub><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>特别地,发射脉冲与目标的距离为R<sub>0</sub>(t):<maths num="0002"><![CDATA[<math><mrow><msub><mi>R</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>b</mi></msub><mo>+</mo><msub><mi>v</mi><mi>r</mi></msub><mi>t</mi><mo>+</mo><mfrac><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mi>t</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>V</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>v</mi><mi>r</mi></msub><mi>t</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>(e)获取任一通道i接收到的SAR回波信号形式:<maths num="0003"><![CDATA[<math><mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mn>0</mn></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j&pi;</mi><msub><mi>K</mi><mi>r</mi></msub><msup><mi>t</mi><mn>2</mn></msup><mo>}</mo><mo>&times;</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三:对步骤一中SAR回波信号进行距离压缩,即乘以exp{jπK<sub>r</sub>t<sup>2</sup>},得到压缩后的信号形式:<maths num="0004"><![CDATA[<math><mrow><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mi>R</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>将式(1)(2)带入式(4)即为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mi>t</mi><mo>+</mo><mfrac><mrow><msup><mrow><msubsup><mi>V</mi><mi>a</mi><mn>2</mn></msubsup><mi>t</mi></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><msub><mrow><mn>2</mn><mi>V</mi></mrow><mi>a</mi></msub><msub><mi>x</mi><mn>0</mn></msub><mi>t</mi><mo>+</mo><msup><msub><mrow><mn>2</mn><mi>d</mi></mrow><mrow><mi>i</mi><mn>0</mn></mrow></msub><mn>2</mn></msup><mo>+</mo><msub><mrow><mn>2</mn><mi>d</mi></mrow><mrow><mi>i</mi><mn>0</mn></mrow></msub><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><msub><mrow><mn>2</mn><mi>d</mi></mrow><mrow><mi>i</mi><mn>0</mn></mrow></msub><msub><mi>V</mi><mi>a</mi></msub><mi>t</mi></mrow><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤四:通道i接收到的回波相对于参考通道0延时为τ<sub>i0</sub>=d<sub>i0</sub>/V<sub>a</sub>,做延时处理得到:<maths num="0006"><![CDATA[<math><mrow><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mi>t</mi><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><msub><mi>&tau;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><msup><mrow><msubsup><mi>V</mi><mi>a</mi><mn>2</mn></msubsup><mi>t</mi></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><msub><mrow><mn>2</mn><mi>V</mi></mrow><mi>a</mi></msub><msub><mi>x</mi><mn>0</mn></msub><mi>t</mi><mo>+</mo><msup><msub><mrow><mn>2</mn><mi>d</mi></mrow><mrow><mi>i</mi><mn>0</mn></mrow></msub><mn>2</mn></msup></mrow><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤五:根据步骤四中信号形式,推导通道i与方位时间无关的相位补偿函数:<maths num="0007"><![CDATA[<math><mrow><msub><mi>C</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mo>&times;</mo><mfrac><msup><msub><mi>d</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mn>2</mn></msup><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤六:根据步骤四中信号形式,推导每个通道与方位时间的二次项有关的相位补偿函数,作用相当于方位压缩:<maths num="0008"><![CDATA[<math><mrow><mi>A</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mo>&times;</mo><mfrac><msup><mrow><msup><msub><mi>V</mi><mi>a</mi></msub><mn>2</mn></msup><mi>t</mi></mrow><mn>2</mn></msup><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤七:将步骤五、六的相位补偿函数与步骤四中的信号做相乘处理后,信号形式为:<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>i</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mi>t</mi><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><msub><mi>&tau;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><msub><mrow><mn>2</mn><mi>V</mi></mrow><mi>a</mi></msub><msub><mi>x</mi><mn>0</mn></msub><mi>t</mi></mrow><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤八:按照步骤二到步骤七推导其他通道k的处理过程,处理后得到信号形式:<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>k</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><msub><mi>&tau;</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>R</mi></mrow><mi>b</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mi>t</mi><mo>+</mo><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><msub><mi>&tau;</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><msub><mrow><mn>2</mn><mi>V</mi></mrow><mi>a</mi></msub><msub><mi>x</mi><mn>0</mn></msub><mi>t</mi></mrow><msub><mi>R</mi><mi>b</mi></msub></mfrac><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤九:将各个通道的信号分别作方位向FFT后,选择两个通道进行干涉处理,使得到的干涉相位中涉及到的若干基线长度可以组成互质基线对:<maths num="0011"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>S</mi><mi>ik</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>S</mi><mi>k</mi></msub><msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>*</mo></msup><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>}</mo><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><msub><mi>&pi;v</mi><mi>r</mi></msub><mn></mn></mrow><mi>&lambda;</mi></mfrac><mfrac><msub><mi>d</mi><mi>ik</mi></msub><msub><mi>V</mi><mi>a</mi></msub></mfrac><mo>}</mo></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>jr</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>S</mi><mi>r</mi></msub><msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>*</mo></msup><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mrow><mn>2</mn><mi>v</mi></mrow><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><msub><mrow><mn>4</mn><mi>&pi;v</mi></mrow><mi>r</mi></msub><mi>&lambda;</mi></mfrac><mfrac><msub><mi>d</mi><mi>jr</mi></msub><msub><mi>V</mi><mi>a</mi></msub></mfrac><mo>}</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中关于互质基线对d<sub>ik</sub>、d<sub>jr</sub>的干涉相位:<maths num="0012"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&Delta;&phi;</mi><mi>ik</mi></msub><mo>=</mo><mfrac><msub><mrow><mn>4</mn><mi>&pi;v</mi></mrow><mi>r</mi></msub><mi>&lambda;</mi></mfrac><mfrac><msub><mi>d</mi><mi>ik</mi></msub><msub><mi>V</mi><mi>a</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>&Delta;&phi;</mi><mi>jr</mi></msub><mo>=</mo><mfrac><msub><mrow><mn>4</mn><mi>&pi;v</mi></mrow><mi>r</mi></msub><mi>&lambda;</mi></mfrac><mfrac><msub><mi>d</mi><mi>jr</mi></msub><msub><mi>V</mi><mi>a</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤十:由干涉相位计算得到运动目标速度:<maths num="0013"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>v</mi><mi>rik</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&lambda;V</mi><mi>a</mi></msub><msub><mi>&Delta;&phi;</mi><mi>ik</mi></msub></mrow><msub><mrow><mn>4</mn><mi>&pi;d</mi></mrow><mi>ik</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>rjr</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&lambda;V</mi><mi>a</mi></msub><msub><mi>&Delta;&phi;</mi><mi>jr</mi></msub></mrow><msub><mrow><mn>4</mn><mi>&pi;d</mi></mrow><mi>jr</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤十一:根据不同基线对所对应的动目标检测结果,采用:(a)长基线对对应的动目标检测结果用于检测慢速目标,选择该基线对内共有的检测速度作为慢速目标速度,以排除模糊速度(虚假目标);(b)短基线对对应的动目标检测结果用于检测速度较快目标,选择该基线对内共有的检测速度作为速度较快目标速度,以排除模糊速度(虚假目标)。
地址 100191 北京市海淀区学院路37号