主权项 |
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>′</mo></msubsup><mo>-</mo><msub><mover><mi>x</mi><mo>‾</mo></mover><mi>j</mi></msub></mrow><msub><mi>σ</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>‾</mo></mover><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>Σ</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>σ</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>Σ</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>‾</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>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>μ</mi><mi>ik</mi></msub><mo>)</mo></mrow><mi>h</mi></msup><msub><mi>d</mi><mi>ik</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>μ</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>μ</mi><mi>ik</mi></msub><mo>)</mo></mrow><mrow><mi>c</mi><mo>×</mo><mi>m</mi></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>μ</mi><mn>11</mn></msub></mtd><mtd><msub><mi>μ</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>μ</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>μ</mi><mn>21</mn></msub></mtd><mtd><msub><mi>μ</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>μ</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>μ</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>μ</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>μ</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:最终可以得到每个无人机的归属类号,将同一类号的无人机归到一起,确定为一个编队,得到所有的编队划分结果。 |