发明名称 FFS空间随机三阵元分布探测方法
摘要 本发明一种FFS空间随机三阵元分布探测方法,属于无线电无源探测技术领域。每个子母弹的子弹丸单独装载一个阵元天线,每个阵元独自工作,通过自身装载的GPS定位系统记录自身地球坐标,并不断记录目标源发射来电磁波相位,并将位置坐标和目标相位信息向周围阵元发射,同时接受来自其他阵元的数据信号。由于指定频率目标源到达不同的阵元存在一个波程差,这个波程差导致了各个弹丸接收阵元间的相位差。每个阵元根据通信数据中的相位差和位置坐标计算出目标信号的入射俯仰角θ及信号的入射方位角<img file="DBM2016033001160000011.GIF" wi="46" he="54" />。本方法提高了探测精度,可适用于引信探测中的阵元位置随机分布情况下对指定电磁辐射源相对位置的探测。<pb pnum="1" />
申请公布号 CN105659860B 申请公布日期 2014.02.12
申请号 CN201010052006.5 申请日期 2010.12.13
申请人 北京理工大学 发明人 潘曦;徐立新;崔占忠;张帆
分类号 G01S5/04(2006.01)I 主分类号 G01S5/04(2006.01)I
代理机构 代理人
主权项 FFS空间随机三阵元分布探测方法,其特征在于:具体实现过程为:在探测过程中,三个炮弹阵元相对独立工作,沿逆时针方向分别设三个炮弹阵元为S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>;以炮弹阵元S<sub>1</sub>作为研究目标,令三个炮弹阵元的外接圆圆心为坐标原点O,将由O点指向S<sub>1</sub>的向量设为X轴,外接圆所在的平面为XOY平面,垂直于该平面的法线方向为Z轴,并根据右手定则确定Y轴,最终确定系统空间直角坐标系;S<sub>1</sub>通过自身坐标定位确定自身在空间中的坐标,同时接收来自目标源的相位信息;并将坐标信息和相位信息编码后以广播的形式向外发射;S<sub>1</sub>还接收来自另外两个炮弹阵元S<sub>2</sub>,S<sub>3</sub>的信息,其中包含另外两个炮弹阵元的坐标GPS位置信息,以及它们接收到的来自信号源的相位信息;设以上三个炮弹阵元在系统坐标系中的坐标分别是:S<sub>1</sub>(x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>),S<sub>2</sub>(x<sub>2</sub>,y<sub>2</sub>,z<sub>2</sub>),S<sub>3</sub>(x<sub>3</sub>,y<sub>3</sub>,z<sub>3</sub>),得出炮弹阵元间的距离为:<maths id="cmaths0001" num="0001"><math><![CDATA[<mrow><msub><mi>S</mi><mn>1</mn></msub><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><mi>a</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mn>2</mn></msub><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FBM2016033001180000031.GIF" wi="932" he="124" /></maths><maths id="cmaths0002" num="0002"><math><![CDATA[<mrow><msub><mi>S</mi><mn>3</mn></msub><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><mi>b</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mn>3</mn></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FBM2016033001180000032.GIF" wi="926" he="132" /></maths><maths id="cmaths0003" num="0003"><math><![CDATA[<mrow><msub><mi>S</mi><mn>3</mn></msub><msub><mi>S</mi><mn>1</mn></msub><mo>=</mo><mi>c</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>-</mo><msub><mi>y</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FBM2016033001180000033.GIF" wi="920" he="119" /></maths>三角形ΔS<sub>1</sub>S<sub>2</sub>S<sub>3</sub>的半周长为:<img file="FBM2016033001180000034.GIF" wi="352" he="118" />ΔS<sub>1</sub>S<sub>2</sub>S<sub>3</sub>面积为<img file="FBM2016033001180000035.GIF" wi="640" he="116" />外接圆的半径为:<img file="FBM2016033001180000036.GIF" wi="176" he="93" />今原点O到三个顶点S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>的向量之间的夹角∠S<sub>1</sub>OS<sub>2</sub>,∠S<sub>2</sub>OS<sub>3</sub>,∠S<sub>3</sub>OS<sub>1</sub>分别为θ<sub>1</sub>,θ<sub>2</sub>,θ<sub>3</sub>,则:<maths id="cmaths0004" num="0004"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><mn>2</mn><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mi>a</mi><mrow><mn>2</mn><mi>r</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>&theta;</mi><mn>2</mn></msub><mo>=</mo><mn>2</mn><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mi>b</mi><mrow><mn>2</mn><mi>r</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>&theta;</mi><mn>3</mn></msub><mo>=</mo><mn>2</mn><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mi>c</mi><mrow><mn>2</mn><mi>r</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FBM2016033001180000037.GIF" wi="1270" he="165" /></maths>三个炮弹阵元S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>分别距离原点O的波程差为:<img file="FBM2016033001180000038.GIF" wi="1648" he="134" />故相邻两炮弹阵元间的波程差为:<img file="FBM2016033001180000039.GIF" wi="954" he="113" /><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mn>1</mn><mo>,</mo><mn>2</mn></mrow></msub><mo>=</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>&times;</mo><msub><msup><mi>x</mi><mo>*</mo></msup><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>s</mi></msub><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>&Delta;</mi><mn>12</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FBM2016033001180000043.GIF" wi="803" he="113" /></maths><maths id="cmaths0006" num="0006"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mn>2</mn><mo>,</mo><mn>3</mn></mrow></msub><mo>=</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>&times;</mo><msub><msup><mi>x</mi><mo>*</mo></msup><mn>3</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>s</mi></msub><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>&Delta;</mi><mn>23</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FBM2016033001180000044.GIF" wi="813" he="116" /></maths><maths id="cmaths0007" num="0007"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mn>3</mn><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>&times;</mo><msub><msup><mi>x</mi><mo>*</mo></msup><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>s</mi></msub><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>&Delta;</mi><mn>31</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FBM2016033001180000045.GIF" wi="800" he="113" /></maths>其中,E代表数学期望,“*”代表复共轭运算,P<sub>s</sub>代表信号功率,λ为信号波长,x<sub>1</sub>(t)、x<sub>2</sub>(t)、x<sub>3</sub>(t)分别为S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>的接收到的相位信息;并且有:<img file="FBM2016033001180000046.GIF" wi="1029" he="142" /><img file="FBM2016033001180000047.GIF" wi="1042" he="135" />其中,求解互相关函数R<sub>1,2</sub>,R<sub>2,3</sub>的具体方法为:接收机的三个炮弹阵元接收到目标发射出的高频电磁辐射波后,分别经过变频转换为低频同相信号后各分为I路和Q路两路信号输出,I,Q两路信号相互正交,用一定的采样频率进行A/D采样,后分别将采样的数据存储在I<sub>1</sub>,Q<sub>1</sub>,I<sub>2</sub>,Q<sub>2</sub>,I<sub>3</sub>,Q<sub>3</sub>数组中;则每个炮弹阵元信号表示为Ch1=I<sub>1</sub>+j·Q<sub>1</sub>Ch2=I<sub>2</sub>+j·Q<sub>2</sub>Ch3=I<sub>3</sub>+j·Q<sub>3</sub>有<img file="FBM2016033001180000048.GIF" wi="1310" he="372" /><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>I</mi><mn>2</mn></msub><mo>=</mo><mi>A</mi><mi>r</mi><mi>g</mi><mrow><mo>(</mo><msub><mi>R</mi><mrow><mn>2</mn><mo>,</mo><mn>3</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>r</mi><mi>g</mi><mo>&lsqb;</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>&times;</mo><msub><msup><mi>x</mi><mo>*</mo></msup><mn>3</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo><mi>A</mi><mi>r</mi><mi>g</mi><mo>&lsqb;</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>Ch</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><msup><mi>Ch</mi><mo>*</mo></msup><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>A</mi><mi>r</mi><mi>g</mi><mo>{</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>I</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>j</mi><mo>&CenterDot;</mo><msub><mi>Q</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&lsqb;</mo><msub><mi>I</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>j</mi><mo>&CenterDot;</mo><msub><mi>Q</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FBM2016033001180000051.GIF" wi="1355" he="378" /></maths>公式中N为采样点数;在上述的无源测向过程中,互相关函数R<sub>1,2</sub>,R<sub>2,3</sub>由接收机的接收数据直接求出,r,θ<sub>1</sub>,θ<sub>2</sub>,θ<sub>3</sub>由三个炮弹阵元S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>的GPS接收坐标信息计算得出;联立以上两个方程,可推出存在未知变量θ、<img file="FBM2016033001180000052.GIF" wi="38" he="45" />满足非线性方程组;在求解未知参量俯仰角θ,入射方位角<img file="FBM2016033001180000053.GIF" wi="37" he="43" />的过程中,由于非线性方程组的复杂程度,无法求解出满足条件的解的精确形式;因此,在探测系统给定求解精度ε的约束条件下,用数值迭代的方法,求出数值解θ<sup>*</sup>,<img file="FBM2016033001180000054.GIF" wi="82" he="66" />使其满足<img file="FBM2016033001180000055.GIF" wi="554" he="86" />pb pnum="3" />
地址 100081 北京市海淀区中关村南大街5号