发明名称 一种多无人机编队划分方法
摘要 本发明提供了一种多无人机编队划分方法,依次获得空战各个目标态势信息数据,并进行标准化处理;初始化隶属矩阵,使其满足约束条件;计算聚类中心;计算价值函数J,若满足计算终止的跳出条件,则算法停止,否则计算新的隶属矩阵,并返回上一步,重新计算聚类中心;满足跳出条件后,得到最终隶属矩阵;根据最终的隶属矩阵内各数据值对多目标进行聚类划分。本发明有效地降低多目标编队执行任务分配时的解算维数,减少运算量,提高解算速度。
申请公布号 CN103324830A 申请公布日期 2013.09.25
申请号 CN201310194166.7 申请日期 2013.05.23
申请人 西北工业大学 发明人 张堃;张才坤;周德云;冯琦;马云红
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.一种多无人机编队划分方法,其特征在于包括下述步骤:步骤1:获得各个无人机态势信息数据;并采用下式标准化<maths num="0001"><![CDATA[<math><mrow><msub><mi>x</mi><mi>kj</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>x</mi><mi>kj</mi><mo>&prime;</mo></msubsup><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub></mrow><msub><mi>&sigma;</mi><mi>j</mi></msub></mfrac></mrow></math>]]></maths>式中,x'<sub>kj</sub>表示第k个无人机的第j个属性值的原始数据,x<sub>kj</sub>表示归一化后的第k个无人机的第j个属性值,1≤k≤m,1≤j≤n,m表示无人机数目,n表示无人机属性,分别取无人机的经度、纬度、高度,<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>x</mi><mi>kj</mi></msub><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>&sigma;</mi><mi>j</mi></msub><mo>=</mo><msqrt><mfrac><mn>1</mn><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>kj</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow></math>]]></maths>步骤2:用值在0,1间的随机数初始化隶属矩阵U=(μ<sub>ik</sub>)<sub>c×m</sub>,使其满足约束条件<img file="FDA00003233492300014.GIF" wi="231" he="149" />μ<sub>ik</sub>∈(0,1),μ<sub>ik</sub>表示第k个无人机属于第i个类的隶属度值,1≤i≤c,c表示聚类数目,<maths num="0004"><![CDATA[<math><mrow><mi>c</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mo>[</mo><mi>m</mi><mo>/</mo><mn>4</mn><mo>]</mo><mo>,</mo></mtd><mtd><mi>m</mi><mo>%</mo><mn>4</mn><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mo>[</mo><mi>m</mi><mo>/</mo><mn>4</mn><mo>]</mo><mo>+</mo><mn>1</mn><mo>,</mo></mtd><mtd><mi>else</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>[]表示取整,%表示取余;步骤3:计算c个聚类中心<img file="FDA00003233492300016.GIF" wi="376" he="278" />式中,1≤i≤c,h为模糊权重系数,取2-6之间,V<sub>i</sub>表示各个属性的第i个聚类中心,V<sub>i</sub>=(v<sub>i1</sub>,v<sub>i2</sub>,…,v<sub>in</sub>),X<sub>k</sub>表示第k架无人机各个属性值,X<sub>k</sub>=(x<sub>k1</sub>,x<sub>k2</sub>,…,x<sub>kn</sub>);步骤4:计算价值函数<maths num="0005"><![CDATA[<math><mrow><mi>J</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>ik</mi></msub><mo>)</mo></mrow><mi>h</mi></msup><msub><mi>d</mi><mi>ik</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>ik</mi></msub><mo>)</mo></mrow><mi>h</mi></msup><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>k</mi></msub><mo>-</mo><msub><mi>V</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>如果J小于阀值ε,其中阈值ε=0.000001,或J相对上次价值函数值的改变量小于ε,则算法停止,跳到步骤6;式中,d<sub>jk</sub>=||X<sub>k</sub>-V<sub>i</sub>||<sup>2</sup>=(X<sub>k</sub>-V<sub>i</sub>)<sup>T</sup>(X<sub>k</sub>-V<sub>i</sub>);步骤5:重新计算隶属矩阵U=(μ<sub>ik</sub>)<sub>c×m</sub>,<img file="FDA00003233492300018.GIF" wi="433" he="165" />返回步骤3;步骤6:得到最终的隶属矩阵<maths num="0006"><![CDATA[<math><mrow><mi>U</mi><mo>=</mo><msub><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>ik</mi></msub><mo>)</mo></mrow><mrow><mi>c</mi><mo>&times;</mo><mi>m</mi></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&mu;</mi><mn>11</mn></msub></mtd><mtd><msub><mi>&mu;</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&mu;</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mn>21</mn></msub></mtd><mtd><msub><mi>&mu;</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&mu;</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mi></mi></mtd><mtd><mi></mi></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mi></mi></mtd><mtd><mi></mi></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&mu;</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>&mu;</mi><mi>cm</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>步骤7:计算第k个无人机类属的标号<img file="FDA00003233492300021.GIF" wi="359" he="102" />式中,c<sub>k</sub>表示,1≤i≤c,<img file="FDA00003233492300022.GIF" wi="232" he="92" />表示1≤i≤c时,μ<sub>ik</sub>最大时的i的取值步骤8:最终可以得到每个无人机的归属类号,将同一类号的无人机归到一起,确定为一个编队,得到所有的编队划分结果。
地址 710072 陕西省西安市友谊西路127号