发明名称 一种基于DEM数据的宽带雷达相干杂波仿真方法
摘要 本发明属于雷达仿真技术领域,具体涉及一种基于DEM数据的宽带雷达相干杂波仿真方法。其包括如下步骤:载入仿真参数;读入第n个航迹点数据;确定雷达平台位置和波束照射方向;结合天线方向图,确定雷达照射区域;以dθ和<img file="DDA00002255411100011.GIF" wi="76" he="56" />分别为俯仰角和方位角的角度间隔沿方位向和俯仰向进行地面散射单元划分,计算散射单元面积;根据雷达方程,计算单元的回波电压功率和幅度;根据得到的杂波幅度,结合LFM信号特性,计算散射单元的回波电压信号;计算散射单元回波在雷达回波数据的位置;对所有散射单元的回波电压信号进行相干合成;判断n<N是否成立,如果成立,那么仿真结束;否则重复步骤二。该方法在考虑地面起伏变化情况下,不进行复杂的电磁散射计算。
申请公布号 CN103729485A 申请公布日期 2014.04.16
申请号 CN201210390929.0 申请日期 2012.10.15
申请人 中国航天科工集团第二研究院二〇七所 发明人 耿旭朴;崔亮;董纯柱
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 1.一种基于DEM数据的宽带雷达相干杂波仿真方法,其特征在于:该方法包括如下步骤:步骤一:载入仿真参数,包括雷达平台运动航迹、速度、天线波束指向、地面网格DEM数据,以及雷达的波长、脉宽、带宽、作用距离和天线方向图;步骤二:读入第n个航迹点的数据;步骤三:确定雷达平台位置和波束照射方向;步骤四:结合天线方向图,确定雷达的照射区域;步骤五:以dθ和<img file="FDA00002255410800011.GIF" wi="60" he="47" />分别为俯仰角和方位角的角度间隔沿方位向和俯仰向进行地面散射单元划分,散射单元的面积为<img file="FDA00002255410800012.GIF" wi="1322" he="106" />其中,θ和<img file="FDA00002255410800013.GIF" wi="30" he="34" />分别表示球坐标系下的俯仰角和方位角,r表示雷达到散射单元的距离;步骤六:根据雷达方程,计算单元的回波电压功率和幅度:<maths num="0001"><![CDATA[<math><mrow><msub><mi>P</mi><mi>r</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>r</mi></msub><msup><mi>G</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msup><mi>&lambda;</mi><mn>2</mn></msup><mi>&sigma;</mi></mrow><mrow><msup><mrow><mo>(</mo><mn>4</mn><mi>&pi;</mi><mo>)</mo></mrow><mn>3</mn></msup><msup><mi>r</mi><mn>4</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><msqrt><msub><mi>P</mi><mi>r</mi></msub></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:P<sub>t</sub>表示雷达发射功率,单位W;G(θ)表示在地面单元方向的天线增益;λ表示雷达工作波长,单位m;σ表示地面单元的散射截面积;r表示地面单元到雷达的距离,单位m;步骤七:根据步骤六得到的杂波幅度A<sub>i</sub>,结合LFM信号特性,按下式计算散射单元的回波电压信号:<img file="FDA00002255410800016.GIF" wi="1587" he="150" />其中,<img file="FDA00002255410800017.GIF" wi="42" he="38" />表示初始相位,k<sub>r</sub>表示LFM信号的调制率,f<sub>d</sub>表示多普勒频率;步骤八:计算散射单元回波在雷达回波数据的位置m<sub>i</sub>=(int){(r-R<sub>start</sub>)/bin+0.5}    (6)其中,(int)表示取整运算,R<sub>start</sub>表示雷达的起始采样距离,bin表示距离门的宽度;步骤九:对所有散射单元的回波电压信号,进行相干合成,第n个航迹点的回波信号为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>S</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤十:判断n<N(N表示全部航迹点数)是否成立,如果成立,那么仿真结束;否则重复步骤二。
地址 100854 北京市海淀区永定路50号