发明名称 |
一种无人机固定半径盘旋领航方法 |
摘要 |
本发明公开了一种无人机固定半径盘旋领航方法,根据既定的盘旋半径、期望的盘旋方式以及由传感器获得的飞机位置、高度及地速信息,连续实时计算飞机期望轨迹的侧向偏离及侧偏移速度等导航参数,并同由姿态运动传感器获得的飞机姿态运动信息一同输入至侧向控制回路得到舵偏度指令,最终实现引导飞机沿着期望航线飞行。应用本发明提供的领航方法能够实现引导飞机进行固定半径盘旋飞行;避免了折线飞行过程中频繁的直飞和转弯动作切换,飞机机体所承受的过载较小,在风干扰情况下,飞机仍然按照期望轨迹飞行。 |
申请公布号 |
CN101261131A |
申请公布日期 |
2008.09.10 |
申请号 |
CN200810103302.6 |
申请日期 |
2008.04.02 |
申请人 |
北京航空航天大学 |
发明人 |
王养柱;崔中兴;方晓星;田波;王瑛 |
分类号 |
G01C21/20(2006.01);G01C21/00(2006.01) |
主分类号 |
G01C21/20(2006.01) |
代理机构 |
北京永创新实专利事务所 |
代理人 |
周长琪 |
主权项 |
1、一种无人机固定半径盘旋领航方法,其特征在于包括如下步骤:步骤一:计算飞机当前位置子午面内地球主曲率半径R<sub>m</sub>以及当前位置子午面垂直平面内的地球主曲率半径R<sub>n</sub>,然后计算飞机盘旋时的圆心位置与当前飞机位置的地理经度偏差量ΔL及地理纬度偏差量ΔB:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>=</mo><msub><mi>R</mi><mi>a</mi></msub><mo>×</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><mo>×</mo><mi>f</mi><mo>+</mo><mn>3</mn><mo>×</mo><mi>f</mi><mo>×</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>B</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>n</mi></msub><mo>=</mo><msub><mi>R</mi><mi>a</mi></msub><mo>×</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>f</mi><mo>×</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>B</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths><maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>ΔL</mi><mo>=</mo><mfrac><mrow><mi>Rturn</mi><mo>·</mo><msub><mi>V</mi><mi>dn</mi></msub></mrow><mrow><msqrt><msubsup><mi>V</mi><mi>dn</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>V</mi><mi>de</mi><mn>2</mn></msubsup></msqrt><mo>·</mo><msub><mi>R</mi><mi>m</mi></msub><mo>·</mo><mi>cos</mi><mi>B</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi>ΔB</mi><mo>=</mo><mfrac><mrow><mi>Rturn</mi><mo>·</mo><msub><mi>V</mi><mi>de</mi></msub></mrow><msqrt><msubsup><mi>V</mi><mi>dn</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>V</mi><mi>de</mi><mn>2</mn></msubsup><mo>·</mo><msub><mi>R</mi><mi>n</mi></msub></msqrt></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>上式中的参数依据WGS_84坐标系,地球长半轴R<sub>a</sub>=6378137.0m,地球椭偏度f=0.003352811;Rturn为既定的盘旋半径,V<sub>dn</sub>为飞机的北向地速,V<sub>de</sub>为飞机的东向地速,B为飞机的地理纬度;步骤二:根据收到的盘旋指令的不同解算相应的盘旋圆心地理经度L<sub>C</sub>和地理纬度B<sub>C</sub>;根据步骤一中的ΔB和ΔL,计算盘旋圆心地理经度L<sub>C</sub>和盘旋圆心地理纬度B<sub>C</sub>为:<img file="A20081010330200023.GIF" wi="597" he="138" /><img file="A20081010330200024.GIF" wi="606" he="139" />其中B为飞机的地理纬度,L为飞机的地理经度;步骤三:根据盘旋圆心地理纬度B<sub>C</sub>计算盘旋圆心处的地心半径B<sub>p</sub>、地心纬度B<sub>ECC</sub>及飞机的地心纬度B<sub>EC</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>R</mi><mi>p</mi></msub><mo>=</mo><msub><mi>R</mi><mi>a</mi></msub><mo>·</mo><mo>{</mo><mn>1</mn><mo>+</mo><mfrac><msup><mi>f</mi><mn>2</mn></msup><mn>2</mn></mfrac><mo>·</mo><mo>[</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>f</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>·</mo><mi>tg</mi><msub><mi>B</mi><mi>C</mi></msub></mrow></mfrac><mo>-</mo><mn>1</mn><mo>]</mo><mo>}</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>B</mi><mi>EC</mi></msub><mo>=</mo><mi>arctg</mi><mo>[</mo><msup><mrow><mo>(</mo><mfrac><mrow><mi>H</mi><mo>+</mo><msub><mi>R</mi><mi>b</mi></msub></mrow><mrow><mi>H</mi><mo>+</mo><msub><mi>R</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>·</mo><mi>tgB</mi><mo>]</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>B</mi><mi>ECC</mi></msub><mo>=</mo><mi>arctg</mi><mo>[</mo><msup><mrow><mo>(</mo><mfrac><mrow><mi>H</mi><mo>+</mo><msub><mi>R</mi><mi>b</mi></msub></mrow><mrow><mi>H</mi><mo>+</mo><msub><mi>R</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>·</mo><msub><mi>tgB</mi><mi>C</mi></msub><mo>]</mo></mrow></math>]]></maths>其中,地球短半轴R<sub>b</sub>=6356752.3m,长半轴R<sub>a</sub>=6378137.0m,H和B分别为飞机高度和纬度;步骤四:计算飞机到盘旋圆心的距离D和飞机的侧偏距D<sub>Z</sub>:<maths num="0006"><![CDATA[<math><mrow><mi>D</mi><mo>=</mo><msqrt><msup><mi>Δx</mi><mn>2</mn></msup><mo>+</mo><mi>Δ</mi><msup><mi>y</mi><mn>2</mn></msup></msqrt><mo>,</mo><mi></mi></mrow></math>]]></maths>其中<maths num="0007"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>Δx</mi><mo>=</mo><msub><mi>R</mi><mi>p</mi></msub><mo>·</mo><mi>cos</mi><msub><mi>B</mi><mi>ECC</mi></msub><mo>·</mo><mi>ΔL</mi></mtd></mtr><mtr><mtd><mi>Δy</mi><mo>=</mo><msub><mi>R</mi><mi>p</mi></msub><mo>·</mo><mrow><mo>(</mo><msub><mi>B</mi><mi>EC</mi></msub><mo>-</mo><msub><mi>B</mi><mi>ECC</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>D<sub>Z</sub>=D-Rturn;步骤五:计算从盘旋圆心到飞机当前位置的向量与东向的夹角α、侧偏移速度D<sub>Zd</sub>:α=arctg(Δy/Δx)<img file="A20081010330200035.GIF" wi="898" he="138" />步骤六:将步骤四及步骤五中得到的飞机的侧偏距D<sub>Z</sub>及侧偏移速度D<sub>Zd</sub>输出至侧向控制回路得到舵偏度指令,最终实现引导飞机沿着期望航线飞行。 |
地址 |
100083北京市海淀区学院路37号 |