发明名称 基于近场SBR算法及双站测试数据的RCS外推方法
摘要 本发明涉及目标电磁散射特性测量技术。所要解决的技术问题是提供基于近场SBR算法及双站测试数据的RCS外推方法,可以解决目标大尺寸部件间多次耦合散射的近场修正问题,提高目标RCS外推数据的精度。其特征在于:紧密结合电磁散射测量与建模两种手段,根据目标几何外形特征,选取目标尺寸较长的两个方向为入射方向,通过分别测量两个入射方向下目标双站散射数据,并利用提取的两组线状散射中心分别拟合,再基于SBR算法原理,计算两组散射中心间的多次散射贡献,可外推大角度范围内不同距离下的RCS结果,且可有效修正常规外推方法难以修正的多次散射误差。<pb pnum="1" />
申请公布号 CN106342320B 申请公布日期 2014.08.20
申请号 CN201110011754.3 申请日期 2011.07.07
申请人 上海无线电设备研究所 发明人 梁子长;岳慧;陈奇平
分类号 G06F19/00(2011.01)I;G01S17/00(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 上海航天局专利中心 31107 代理人 张绪成
主权项 基于近场SBR算法及双站测试数据的RCS外推方法,其特征在于:根据目标几何外形特征,选取目标尺寸较长的两个方向为入射方向,测量获取两个入射角下的目标双站RCS测试数据,再通过最小二乘估计法分别获取两入射角下目标线状散射中心幅相分布数据,在此基础上,采用近场SBR算法计算两组线状散射中心间的多次散射贡献,以提高目标RCS外推数据的精度;包括以下步骤:步骤一:测量入射方向的选择选取目标尺寸较长的两个方向为入射方向,两选定方向间的夹角记为θ<sub>12</sub>;同时,利用激光全站仪测定上述选定方向上目标的横向长度,分别记为L<sub>1</sub>、L<sub>2</sub>;步骤二:一个入射方向下目标双站RCS测量测量时将发射天线固定,将接收天线安装在直线扫描架上不同于发射天线高度的位置,使其进行水平扫描,并调整观测距离,这里要求接收天线相对的散射角范围大于20°;再将目标装订,并根据步骤一选定的入射方向设定目标方位,并固定目标位置;然后,利用扫描架带动接收天线等间隔运动,并依次测量获取目标的双站散射数据V<sub>t1</sub>(y),对应的散射角范围记为‑θ<sub>s1</sub>≤θ≤θ<sub>s2</sub>;步骤三:另一个入射方向下目标双站RCS测量根据步骤一选定的另一入射方向,将目标旋转一定角度‑θ<sub>12</sub>,并固定目标位置;<math><![CDATA[<mrow><msub><mi>l</mi><mrow><mi>t</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>E</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>E</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>t</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>V</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>l</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FWW0000000068000000021.GIF" wi="1188" he="124" /></maths><maths id="cmaths0002" num="0002"><math><![CDATA[<mrow><msub><mi>l</mi><mrow><mi>t</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>t</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>V</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>V</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>l</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FWW0000000068000000022.GIF" wi="1062" he="125" /></maths>式中,l<sub>0</sub>(y)为定标体的双站散射复函数;步骤五:目标线状散射中心参数提取假定在一入射方向下,目标由一组沿直线等间隔分布的散射中心表示,其幅度相位记为a<sub>1</sub>(ξ<sub>i</sub>),ξ<sub>i</sub>=L<sub>1</sub>(i‑1)/(N<sub>1</sub>‑1),i=1,2,...,N<sub>1</sub>;并将其拟合结果与测试数据l<sub>t1</sub>(y)的偏差表示为:<maths id="cmaths0003" num="0003"><math><![CDATA[<mrow><mi>e</mi><mi>r</mi><mi>r</mi><mo>=</mo><mi>&Sigma;</mi><mo>|</mo><msub><mi>l</mi><mrow><mi>t</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mfrac><mrow><mi>exp</mi><mo>&lsqb;</mo><msub><mi>ikR</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><msup><mo>|</mo><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FWW0000000068000000023.GIF" wi="1335" he="154" /></maths>式中,R<sub>1</sub>为各散射中心至观测点间的距离;k为入射电磁波波数;再利用最小二乘估计法获取式(5)取最小值时的散射点幅度相位参数,记为a<sub>1m</sub>(ξ<sub>i</sub>);同理,可获取另一入射方向下线状散射中心分布参数a<sub>2m</sub>(ξ<sub>i</sub>);<math><![CDATA[<mrow><msub><mi>l</mi><mrow><mi>f</mi><mi>i</mi><mi>r</mi><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mfrac><mrow><mi>exp</mi><mo>&lsqb;</mo><msub><mi>ikR</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><msub><mi>R</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mn>0</mn><mo>&le;</mo><mi>&theta;</mi><mo>&le;</mo><msub><mi>&theta;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>/</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mi>a</mi><mo>)</mo></mrow></mrow>]]></math><img file="FWW0000000068000000031.GIF" wi="1403" he="137" /></maths><maths id="cmaths0005" num="0005"><math><![CDATA[<mrow><msub><mi>l</mi><mrow><mi>f</mi><mi>i</mi><mi>r</mi><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><msub><mi>a</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mfrac><mrow><mi>exp</mi><mo>&lsqb;</mo><msub><mi>ikR</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mrow><msub><mi>R</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><msub><mi>&theta;</mi><mn>12</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>/</mo><mn>2</mn><mo>&le;</mo><mi>&theta;</mi><mo>&le;</mo><msub><mi>&theta;</mi><mi>12</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FWW0000000068000000032.GIF" wi="1402" he="135" /></maths>d)计算线状散射中心组间二次散射贡献:<img file="FWW0000000068000000033.GIF" wi="1602" he="168" />e)合成外推计算的后向RCS:σ=4πR<sup>2</sup>|l<sub>first</sub>(θ)+l<sub>second</sub>(θ)|<sup>2</sup>      (8)式中,R为观测点到目标中心或转台中心的距离。pb pnum="3" />
地址 200090 上海市杨浦区黎平路203号