发明名称 基于模糊虚拟力的无人机航路规划方法
摘要 解决基于虚拟力航路规划方法出现的局部极小问题,并实现规划过程中实时的自适应规划参数设置,本发明采用基于模糊虚拟力(FVF)的无人机航路规划方法,其采用贝叶斯网络和模糊逻辑推理相结合进行实时的自适应规划参数设置,并提出威胁合并法解决虚拟力法的局部极小问题。本发明的基于模糊虚拟力的无人机航路规划方法包含以下步骤:1)设置无人机航路规划的初始条件,包括规划起始点、目标点、威胁分布及属性;2)设置无人机航路规划的迭代步长;3)设置规划参数,从而确定虚拟斥力系数与虚拟引力系数之间的关系;4)进行航路规划;5)判断是否进入局部极小,若是,则进行威胁合并,若否,则继续进行航路规划直至目标点。
申请公布号 CN102269593B 申请公布日期 2014.03.12
申请号 CN201010187959.2 申请日期 2010.06.01
申请人 北京航空航天大学 发明人 董卓宁;陈宗基;周锐;李卫琪
分类号 G01C21/24(2006.01)I 主分类号 G01C21/24(2006.01)I
代理机构 代理人
主权项 1.一种基于模糊虚拟力的无人机航路规划方法,其特征在于包含以下步骤:1)设置无人机航路规划的初始条件,包括规划起始点、目标点、威胁分布及属性;2)设置无人机航路规划的迭代步长;3)设置规划参数k,从而确定虚拟斥力系数与虚拟引力系数之间的关系,其中,k=G<sub>R</sub>/G<sub>A</sub>,G<sub>A</sub>表示虚拟引力系数,G<sub>R</sub>表示虚拟斥力系数;4)进行航路规划,航路规划坐标变化量Δx和Δy为:<![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;x</mi><mo>=</mo><mi>&delta;</mi><mo>&CenterDot;</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>F</mi><mi>Ax</mi></msub><mo>+</mo><mi>&Sigma;</mi><msub><mi>F</mi><mi>Rx</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;y</mi><mo>=</mo><mi>&delta;</mi><mo>&CenterDot;</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>F</mi><mi>Ay</mi></msub><mo>+</mo><mi>&Sigma;</mi><msub><mi>F</mi><mi>Ry</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>其中,<![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>F</mi><mi>Ax</mi></msub><mo>=</mo><msub><mi>G</mi><mi>A</mi></msub><mo>&CenterDot;</mo><mi>cos</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>A</mi></msub><mo>)</mo></mrow><mo>/</mo><msup><msub><mi>R</mi><mi>A</mi></msub><mn>2</mn></msup></mtd></mtr><mtr><mtd><msub><mi>F</mi><mi>Ay</mi></msub><mo>=</mo><msub><mi>G</mi><mi>A</mi></msub><mo>&CenterDot;</mo><mi>sin</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>A</mi></msub><mo>)</mo></mrow><mo>/</mo><msup><msub><mi>R</mi><mi>A</mi></msub><mn>2</mn></msup></mtd></mtr></mtable></mfenced></math>]]></maths>F<sub>Ax</sub>和F<sub>Ay</sub>分别表示虚拟引力在x轴和y轴的投影,R<sub>A</sub>是当前点和目标点间的距离,θ<sub>A</sub>是当前点与目标点连线和x轴的夹角;<![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>F</mi><mi>Rx</mi></msub><mo>=</mo><mo>-</mo><msub><mi>G</mi><mi>R</mi></msub><mo>&CenterDot;</mo><mi>cos</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>R</mi><mi>R</mi></msub><mo>/</mo><msub><mi>r</mi><mn>0</mn></msub><mo>)</mo></mrow></msup></mtd></mtr><mtr><mtd><msub><mi>F</mi><mi>Ry</mi></msub><mo>=</mo><mo>-</mo><msub><mi>G</mi><mi>R</mi></msub><mo>&CenterDot;</mo><mi>sin</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>R</mi><mi>R</mi></msub><mo>/</mo><msub><mi>r</mi><mn>0</mn></msub><mo>)</mo></mrow></msup></mtd></mtr></mtable></mfenced></math>]]></maths>F<sub>Rx</sub>和F<sub>Ry</sub>分别表示虚拟斥力在x轴和y轴的投影,R<sub>R</sub>是当前点和威胁间的距离,r<sub>0</sub>是常量可设置为威胁的半径,θ<sub>R</sub>是当前点与威胁连线和x轴的夹角;δ为航路规划的迭代步长,且α=[(F<sub>Ax</sub>+∑F<sub>Rx</sub>)<sup>2</sup>+(F<sub>Ay</sub>+∑F<sub>Ry</sub>)<sup>2</sup>]<sup>-1/2</sup>;5)判断是否进入局部极小,若是,则按照步骤5.1)和5.2)进行威胁合并,若否,则继续进行航路规划直至目标点;5.1)根据任意两个威胁所对应的规划参数k<sub>1</sub>、k<sub>2</sub>模糊推理得到两个威胁间CritiDis为临界间距;定义<img file="FSB0000112956800000014.GIF" wi="710" he="144" />为威胁间距,其中TTDis为两个威胁中心之间的距离,GapDis为两个威胁影响范围间的最短距离,r和R分别为两个威胁的影响半径;5.2)将ThrtDis进行归一化处理,使之处于[0,1]之间,比较任意两个威胁对应的ThrtDis和依据重心法解模糊后的CritiDis之间关系,得到合并标志CmbVal为:<![CDATA[<math><mrow><mi>CmbVal</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>ThrtDis</mi><mo>&le;</mo><mi>CritiDis</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>ThrtDis</mi><mo>></mo><mi>CritiDis</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>合并标志CmbVal为1的两个威胁构成一个威胁组。
地址 100191 北京市海淀区学院路37号