一种基于粒子群算法的机载武器系统精度分配方法

本发明涉及武器分析与分配领域，尤其是一种机载武器系统精度分配方法。

背景技术：

1、精度是衡量机载武器系统性能的重要指标，是影响装备作战效能的关键因素。机载武器系统精度分析和分配是机载武器系统总体方案论证、工程研制和设计定型工作的重要组成部分，在武器系统分析中占有重要的地位，贯穿于武器装备研制的全过程中。传统的精度分析和分配主要依靠飞行试验数据，受经济性、研制周期的制约，仅仅依靠试验的方法无法满足现代高科技条件下武器攻击精度分析的要求。同时，传统的精度分配方法主要是按照等作用原则、等精度原则等进行分配，考虑的因素过少，只适用于结构简单的系统。机载武器系统结构复杂，而且具有较强的非线性，因此传统方法难以求解机载武器系统精度分配问题。

2、粒子群算法是近年来提出的一种启发式搜索算法，通过种群粒子中个体的交互作用来寻找复杂问题空间中的优化解。该算法随机产生一个初始种群并赋予每个微粒一个随机速度，在飞行过程中，粒子的飞行速度和轨迹通过自己及同伴的飞行经验来动态调整，整个群体有飞向更好搜索区域的能力。与传统精度分配方法相比，粒子群算法具有自适应、计算简便、收敛速度快等优势，求解复杂问题时也能够具有较高的效率和准确性，能够高效、快速地求解出精度分配结果。

技术实现思路

1、为了克服现有技术的不足，本发明提供一种基于粒子群算法的机载武器系统精度分配方法。

2、本发明解决其技术问题所采用的技术方案的具体步骤如下：

3、步骤1.数学模型建立：

4、机载武器系统的输入参数为u＝(u1,u2,...,un)t，其中n为输入参数的维数，输出参数为w＝(xe,ze)t，xe、ze分别为命中点在北天东地理坐标系的x轴、z轴的坐标；建立机载武器系统数学模型为w＝f(u)，并确定输入误差σu＝(σ1,σ2,...,σn)t，以均方差表示输入误差的大小；

5、步骤2.利用蒙特卡洛法进行精度分析：

6、步骤3.采用单因素方差分析方法进行误差源敏感性分析：

7、计算各输入误差的f检验值，设置显著水平并查表求出临界值fα，若f≥fα则认为该输入误差源对系统精度的影响显著，否则认为影响不显著；

8、对所有输入误差源，均采用单因素方差分析方法进行敏感性分析，将对系统精度影响显著的输入误差源数量记为n1，对应的均方差记为

9、步骤4.粒子群初始化：

10、将σs作为决策变量，进行粒子群初始化；设置粒子位置和速度的取值范围以及粒子群规模m，随机产生m个n1维向量x1,x2,…,xm作为粒子的初始位置，并生成粒子对应的初始速度向量v1,v2,…,vm，通过精度分析计算出机载武器系统的cep，以cep最小为优化目标，建立目标函数fitness＝min{cep}，以每个粒子在历史迭代中的目标函数最小值为当前的个体最优值，以所有粒子在历史迭代中的目标函数最小值为整个种群的全局最优值；

11、步骤5.粒子速度和位置更新：

12、通过速度更新公式和位置更新公式对粒子的速度和位置进行更新；

13、步骤6.个体最优值和全局最优值更新：

14、根据每个粒子的位置重新计算粒子的目标函数值，设第l个粒子在第t次迭代时的目标函数值为fitnessl(t)，个体最优值为fitnesspbestl，个体最优值为粒子在历史迭代过程中的目标函数最小值，全局最优值为整个种群的所有粒子在历史迭代中的目标函数最小值，整个种群的全局最优值为fitnessgbest，若满足fitnessl(t)<fitnesspbestl，则令fitnesspbestl＝fitnessl(t)，pl＝xl；若满足fitnessl(t)<fitnessgbest，则令fitnessgbest＝fitnessl(t)，pg＝xl；若两个条件均不满足，则执行步骤7；

15、步骤7.判断终止条件：

16、设置精度指标要求cepmax，判断两个终止条件是否满足，两个终止条件为粒子的目标函数值fitnessl(t)≤cepmax和迭代次数t≥tmax，若满足至少一个终止条件，则仿真结束，令σs＝pg，其他误差源取值不变，输出所有误差均方差作为最终的精度分配结果；若两个终止条件均不满足，令t加1，返回步骤5继续进行更新。

17、所述步骤2中，利用蒙特卡洛法进行精度分析的步骤为：

18、根据机载武器系统的输入参数理论值u0＝(u10,u20,...,un0)t及输入误差σu，随机产生p组正态分布随机数ui＝(u1i,u2i,...,uni)t，i＝1,2,...,p，从而得到输入参数的抽样值，其中p为蒙特卡洛统计次数；将输入参数理论值u0＝(u10,u20,...,un0)t和抽样值分别代入机载武器系统数学模型中，计算得到输出参数的理论值w0＝(xe0,ze0)t和实际值wi＝(xei,zei)t，i＝1,2,...,p，从而求出命中点坐标的误差为：

19、△xi＝xei-xe0

20、△zi＝zei-ze0

21、通过下列公式计算输出误差：

22、

23、

24、

25、

26、

27、其中，分别为命中点纵向误差均值和横向误差均值，σx、σz分别为命中点纵向误差均方差和横向误差均方差，cep(circular error probable)为圆概率偏差。

28、所述步骤3中，将第d项输入误差σd,d＝1,2,...,n平均划分为r个不同的误差水平，在每个误差水平下进行s次重复试验，每次试验通过精度分析，计算出机载武器系统的cep，将在第j个误差水平下，第k次重复试验时计算得到的cep记为ajk，j＝1,2,…,r；k＝1,2,…,s，由下式对样本进行统计处理：

29、

30、其中，sst为总离差平方和，ssa为组间离差平方和，sse为组内离差平方和，dft为总体自由度，dfa为组间自由度，dfe为组内自由度，msa为组间均方，mse为组内均方，为总平均值，为组内平均值；

31、经过统计处理后，根据样本数据计算得到关于误差σd的f检验值为：

32、

33、设置显著水平为α，通过查f分布表求出临界值若则认为该输入误差源对机载武器系统精度的影响显著；若则认为影响不显著。

34、所述步骤5中，通过速度更新公式和位置更新公式对粒子的速度和位置进行更新：

35、vl(t+1)＝wvl(t)+λ(pl(t)-xl(t))+(1-λ)(pg(t)-xl(t))

36、xl(t+1)＝xl(t)+vl(t+1)

37、其中，xl、vl分别为第l个粒子的位置和速度向量，l＝1,2,...,m，t为当前迭代次数，w为惯性权重，pl为粒子在历史中的最好位置，pg为整个种群的最好位置，λ为调整因子；

38、所述步骤5中，惯性权重和调整因子的计算公式如下：

39、

40、

41、其中，wmax、wmin分别为惯性权重的最大值和最小值，tmax为最大迭代次数。

42、本发明的有益效果在于采用粒子群算法对机载武器系统进行精度分配，自适应地调整误差源指标，能够得到满足精度要求的分配结果，提高机载武器系统的精度，而且还可以节约成本，最大限度地提升机载武器系统的作战效能；同时，本发明对传统粒子群算法进行了改进，在速度更新公式中引入线性惯性权重和调整因子，能够动态调整局部寻优能力和全局寻优能力，使得粒子群的搜索能力得到有效增强，可以更加快速地给出求解结果，节省了计算时间。