发明名称 风场影响下微型无人飞行器的多分辨率路径规划方法
摘要 本发明公开了一种风场影响下微型无人飞行器的多分辨率路径规划方法,该方法有六大步骤:步骤一:获取飞行器路径规划所需要的环境区域;步骤二:选取小波变换的基函数;步骤三:构建风场影响下的飞行环境代价函数;步骤四:利用Haar小波网格分解法重构多分辨率飞行环境代价函数;步骤五:利用多分辨率环境代价函数值生成环境加权图;步骤六:寻找最优路径。本发明方法在确定环境区域不同分辨率等级过程中考虑了机载计算机的计算能力、传感器的探测范围、障碍物的分布、风场的分布及强度等因素,提高了路径规划的效率,且该方法简单可靠实用,适合于风场影响下的低空微型无人飞行器在线路径规划。
申请公布号 CN102622653A 申请公布日期 2012.08.01
申请号 CN201210046577.7 申请日期 2012.02.27
申请人 北京航空航天大学 发明人 龚光红;高栋栋;王江云;韩亮
分类号 G06Q10/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京永创新实专利事务所 11121 代理人 周长琪
主权项 1.一种风场影响下微型无人飞行器的多分辨率路径规划方法,其特征在于,本方法包括如下步骤:步骤1:获取飞行器路径规划所需要的环境区域;具体本步骤包括如下步骤:步骤1.1:将飞行器所处的三维环境空间区域,从高度方向上分解为不同层次的二维欧氏空间,所得到的二维欧氏空间就是所需要的环境区域;步骤1.2:将得到的二维欧氏空间根据地形高程和飞行器可飞行高度,划分为障碍区和非障碍区,并确定飞行器的初始位置和目标位置;对原始数字地形高程图中第i点位置X<sub>i</sub>的划分如下:<img file="FDA0000138455160000011.GIF" wi="1189" he="107" />其中,<img file="FDA0000138455160000012.GIF" wi="39" he="40" />表示障碍物空间,<img file="FDA0000138455160000013.GIF" wi="45" he="40" />表示非障碍物空间,q<sub>i</sub>表示位置X<sub>i</sub>的地形高程,H表示飞行器的可飞行高度;步骤1.3:将环境区域归一化为单位方形区域,并进行网格划分;步骤2:选取Haar小波变换基函数,构建二维空间的Haar小波基函数,并用二维空间的小波变换基函数对二元函数f(x,y)近似逼近,得到:<img file="FDA0000138455160000014.GIF" wi="708" he="68" />(2)<img file="FDA0000138455160000015.GIF" wi="775" he="70" />其中,(x,y)表示二维空间的位置变量,m,n为二维空间两个方向的位置标号,<img file="FDA0000138455160000016.GIF" wi="37" he="40" />表示整数,Jmin表示最低分辨率等级;<img file="FDA0000138455160000017.GIF" wi="106" he="68" />为近似逼近中二维空间Haar尺度函数<img file="FDA0000138455160000018.GIF" wi="243" he="67" />的系数,<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><msub><mi>J</mi><mi>min</mi></msub></msubsup><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></msubsup><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></msubsup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msubsup><mi>&Phi;</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><msub><mi>J</mi><mi>min</mi></msub></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>dx dy</mi><mo>;</mo></mrow></math>]]></maths><img file="FDA00001384551600000110.GIF" wi="99" he="68" />为近似逼近中二维空间Haar母小波函数<img file="FDA00001384551600000111.GIF" wi="215" he="68" />的系数,<img file="FDA00001384551600000112.GIF" wi="840" he="69" />步骤3:考虑飞行器与障碍物之间的距离、风场的强度、以及风场的梯度三种因素的影响,构建飞行环境代价函数f(X):f(X)=α·rm<sub>1</sub>(X)+β·rm<sub>2</sub>(X)+γ·rm<sub>3</sub>(X)         (3)其中,X表示当前飞行器的位置坐标,X=(x,y),α,β和γ是调节三种因素在代价函数中相对重要性的参数,取值范围为0~1之间;rm<sub>1</sub>(X)表示飞行器与障碍物之间的距离因素对飞行环境的影响,rm<sub>2</sub>(X)表示风场的强度对飞行环境的影响,rm<sub>3</sub>(X)表示风场的梯度对飞行环境的影响;步骤4:利用Haar小波网格分解法重构多分辨率飞行环境代价函数,包括如下步骤:步骤4.1:根据原始数字地形高程图的分辨率和机载计算机的计算资源确定环境区域表示的最高分辨率等级J<sub>max</sub>,根据规划问题的具体需求确定环境区域表示的最低分辨率等级J<sub>min</sub>;步骤4.2:将任意一个分辨率等级的网格单元、同一层次分辨率区域范围、以及同一层次分辨率区域中心坐标进行公式化表示;在分辨率等级为j的情况下,任一网格单元<img file="FDA0000138455160000021.GIF" wi="97" he="68" />表示为:<img file="FDA0000138455160000022.GIF" wi="1201" he="68" />其中,分辨率等级j满足条件:J<sub>min</sub>≤j≤J<sub>max</sub>;<img file="FDA0000138455160000023.GIF" wi="165" he="64" />分别表示在分辨率等级为j情况下,m方向、n方向的空间间隔,<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>h</mi><mi>m</mi><mi>j</mi></msubsup><mo>=</mo><mo>[</mo><mi>m</mi><mo>/</mo><msup><mn>2</mn><mi>j</mi></msup><mo>,</mo><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msup><mn>2</mn><mi>j</mi></msup><mo>]</mo><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>h</mi><mi>n</mi><mi>j</mi></msubsup><mo>=</mo><mo>[</mo><mi>n</mi><mo>/</mo><msup><mn>2</mn><mi>j</mi></msup><mo>,</mo><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><msup><mn>2</mn><mi>j</mi></msup><mo>]</mo><mo>;</mo></mrow></math>]]></maths>在分辨率等级为j的情况下,同一层次分辨率区域中心坐标的集合R(X,r<sup>j</sup>)为:<img file="FDA0000138455160000026.GIF" wi="1271" he="62" />其中,r<sup>j</sup>表示在同一层次分辨率等级为j的情况下的区域范围,r<sup>j</sup>>r<sup>j+1</sup>,Y为某个位置坐标,<img file="FDA0000138455160000027.GIF" wi="54" he="36" />表示环境区域;步骤4.3:将步骤3得到的飞行环境代价函数,采用步骤2中的二维空间的Haar小波变换基函数近似逼近,得到多分辨率飞行环境代价函数,得到:<img file="FDA0000138455160000028.GIF" wi="699" he="136" /><img file="FDA0000138455160000029.GIF" wi="1263" he="220" />其中,N(J<sub>min</sub>)表示在分辨率等级为J<sub>min</sub>的情况下的位置标号的取值集合,根据式(7)得到:<img file="FDA00001384551600000210.GIF" wi="1091" he="62" />G(j)表示在分辨率等级为j的情况下,考虑了障碍物和风场区域后,m方向的位置标号取值序列,根据式(8)得到:<img file="FDA00001384551600000211.GIF" wi="908" he="68" />(8)<img file="FDA00001384551600000212.GIF" wi="681" he="67" />H(j)表示在分辨率等级为j情况下,考虑了障碍物和风场区域后,n方向的位置标号取值序列,根据式(9)得到:<img file="FDA00001384551600000213.GIF" wi="892" he="67" />(9)<img file="FDA00001384551600000214.GIF" wi="672" he="68" />其中,式(8)和(9)中的z<sub>x</sub>,z<sub>y</sub>分别表示风场中心位置的横、纵坐标;步骤5:确定各网格单元的飞行环境代价函数值,将环境区域根据障碍物与非障碍物的表示值,以及飞行环境代价函数值生成环境加权图,将多分辨率路径规划问题转化为图的搜索问题;步骤6:基于A*或Dijkstra搜索方法寻找最优路径,所得到的最优路径就是为飞行器所规划的路径。
地址 100191 北京市海淀区学院路37号