基于能力的国防资源分配方法
源自:系统工程与电子技术
作者:张玉婷 杨镜宇
摘要
针对体系建设发展过程中国防资源有限、分配困难问题, 考虑国防总经费预算受限、年度资源平均分配、体系能力满足能力需求等方面的约束, 首先采用数学规划法, 构建基于能力的国防资源静态分配模型, 将资源分配到不同能力领域的体系建设发展中。其次, 给出了基于遗传算法的模型求解算法。最后, 通过示例验证所构建的资源分配模型能够合理地进行国防资源分配, 且计算结果高效、可行,为体系建设发展中的资源分配问题提供了一定的理论支撑。
关键词
资源分配 ; 数学规划法 ; 体系能力 ; 遗传算法
0 引言
资源分配问题源于经济学, 旨在将有限的资源, 如人力、物力等, 通过一定的规则分配给众多对象, 进行资源优化配置, 以提升经济效益。作战体系结构复杂、规模庞大, 在有限的国防预算约束下, 对其进行资源分配及统筹规划是体系建设发展中的一项重要内容, 对于提升国防资源利用率、提高军事效益、形成强大的军队战斗力具有重要的现实意义。
现实生活中, 各行各业都存在资源分配问题, 许多学者也进行了一定程度的应用研究, 并取得了一些成果。文献[1]提出了一种与资源管理相关的混合整数非线性数学模型, 为救灾行动等不确定场景下的资源分配选址问题提供了帮助。文献[2]针对土地资源分配问题, 根据当前政策和未来趋势, 构建“粮食-水-能源”关系框架, 对光伏产量、可行性、利润、生物质产量进行了全面评估。文献[3]基于“灰色-自回归移动平均模型”算法, 构建改进的集装箱维修费用分配模型。文献[4]在贪婪调度算法的基础上, 研究下一代光学数据中心网络资源调度问题, 引入并行技术, 加快了资源调度速率。文献[5]为高通量多波束通信卫星系统构建资源分配模型, 并通过遗传算法对模型进行求解。文献[6]基于改进萤火虫优化算法, 对云计算中的动态资源分配问题进行了研究, 提出一种基于负载平衡优化的改进Firefly算法, 提高了资源利用效率和生产力。文献[7]为武器装备经费分配问题构建数学规划模型, 并通过差分进化方法求解模型。文献[8]使用机器学习算法预测了医院每日急诊就诊资源分配策略。文献[9]对战争不确定因素下的国防和民用支出资源分配问题进行了研究。文献[10]通过引入Agent构建优化框架, 为解决军事规划中的成本效益优化问题提供了思路。文献[11]围绕战略管理环节评估, 提出一种资源配置需求的形式化描述和探索性方法。文献[12]利用强化学习的方法, 设计了一种基于异步深度强化学习算法的资源分配方法, 实现了资源的自动规划。文献[13]开发了一种通过统计矩传播不确定性的方法, 优化了建模和分析中的系统资源分配方案。文献[14]给出了一种基于Web的实时资源分配工具, 辅助指挥员完成复杂的资源分配和运输任务。
可见, 当前资源分配问题的研究思路和方法众多, 且在军事方面已有学者研究了火力资源分配、装备经费分配、卫星资源分配等[15-20]相关问题, 但对于联合作战体系建设发展层面的资源分配问题尚缺乏具体的模型和算法。
本文从体系建设发展的顶层设计角度出发, 借鉴“基于能力的规划(capability-based planning, CBP)”[21-22]思想, 在国防预算的约束下, 通过分析各体系能力需求, 构建体系资源分配的数学规划模型, 并采用遗传算法对模型进行求解, 得到体系建设发展的帕累托最优资源分配方案及体系总能力最大值。
1 基于数学规划法的国防资源静态分配模型
1.1 问题描述
体系建设发展中基于能力的国防资源分配问题可描述为: 在一个五年规划内, 在总预算B有限的约束下, 以体系能力C作为输入, 如何对N种不同类型的体系/系统进行合理的资源配置, 达到预算最低、体系总能力最高的目标。基于能力的国防资源分配问题如图 1所示。
图1 基于能力的国防资源分配问题
1.2 符号和决策变量说明
1.3 基于数学规划法的资源分配模型
1.3.1 条件声明
(1) 假定能力生成过程为理想化的累积生成过程, 所有种类的能力构成体系总能力。
(2) 假定能力值为标量值, 具备可加性。
(3) 假定国防经费在一个五年规划内按年度平均分配。
(4) 假定在一个五年规划内, 每个二级体系仅有一次建设发展机会, 即选择五年内的某一年进行发展。
1.3.2 模型构建
针对体系建设发展中的国防资源分配问题, 采用数学规划法构建模型。
(1) 目标函数: 体系建设发展的体系总能力达到最大值
(1)
式(1)为目标函数, 其中表示在一个五年规划内所有体系对应的所有能力, 即体系总能力, 目标是寻求该值的最大值。
(2) 约束1: 所有体系的实际建设发展费用不超过总预算
(2)
式中: 表示所有体系中的所有二级体系的实际建设发展费用需求之和,表示一个五年规划内所有体系中的所有二级体系的实际建设发展费用需求之和, 该值不超过一个五年规划的总预算。
(3) 约束2: 年度经费平均分配, 且在允许范围内波动
(3)
式中: 表示分配给每个年度用于体系建设发展的经费总数, 该值在[B(1-δ)/5, B(1+δ)/5]范围内波动。
(4) 约束3: 一级体系对应的能力加权和满足最低需求
(4)
式中: 表示第u年, 一级体系Si所对应的第k种体系能力加权和, 该值不小于最低需求值Luk, Si。
(5) 约束4: 五年规划内的第u年, 所有体系对应的k种能力加权和满足最低需求
(5)
式中: 表示第u年, 所有一级体系所包含的二级体系对应的第k种能力加权和, 该值不小于最低需求值Laku。
(6) 约束5: 一个五年规划内, 每个二级体系仅发展1次
(6)
式中: 表示五年内每个二级体系是否得到发展, 其值为1表示仅发展一次。
2 基于遗传算法的模型求解算法
2.1 遗传算法
遗传算法是Holland[23]受自然进化理论启发而提出的一种搜索算法, 是启发式算法的一种, 通过模仿自然选择和繁殖过程, 进行染色体选择、交叉、变异等运算, 并不断迭代寻找问题的全局最优解。相比传统搜索和优化算法, 遗传算法有许多优点, 如能够高效处理参数规模大、数学表达式复杂等问题。遗传算法的流程如图 2所示。
图2 遗传算法流程图
国内外学者运用遗传算法解决了许多现实问题。文献[24]提出一种基于遗传算法的增强型资源分配算法, 实现了云数据中心资源的高效分配。文献[25]设计了一种具有调度规则和双亲进化的改进遗传算法, 有效逼近了网络效用和收敛速率。文献[26]基于遗传算法开发了一种可靠的柔性供应链网络设计模型, 能够在一定程度上增强战略能力。文献[27]针对颌面部弹性体颜色配方问题, 通过遗传算法开发了一种临床应用程序, 辅助临床医生制定颜色并降低成本。文献[28]对P2P网络借贷问题进行了研究, 运用遗传算法优化了传统网贷预警模型, 提升了预警效果。文献[29]采用遗传算法, 对数据预处理过程的特征选择和特征提取步骤进行了改进。
本文提出的国防资源分配模型中的数学表达式复杂, 用传统算法求解较困难, 容易导致组合爆炸, 因此本文采用遗传算法进行编程, 进行模型求解。为了使算法与模型较好地融合, 对遗传算法进行相关设计, 如确定初始种群、确定适应度函数、确定终止条件, 进行编码和解码、选择、交叉、变异等。
2.2 编码和解码
通过将编码过程转化为计算机可识别的语言, 本文将上述数学模型转化为遗传空间的染色体。本文的决策变量为xuSij, 其表示第Sij个二级体系, 第u年是否被投资, 为0-1变量, 若使用二进制0-1编码, 则染色体的编码所需要的决策变量个数为
(7)
考虑到约束5, 每个二级体系仅能发展1次, 故可以使用实数编码, 将决策变量中的参数u转换为1-5的实数, 即针对一个二级体系, 分配一个1-5的实数决定其在哪一年进行发展, 于是决策变量个数可缩小至之前的20%。
例如, 采用0-1编码: (xSij1, xSij2, xSij3, xSij4, xSij5), 则S11在第2年投资为(0, 1, 0, 0, 0), S12在第3年投资为(0, 0, 1, 0, 0), S13在第1年投资为(1, 0, 0, 0, 0);采用实数编码: 则S11在第2年投资为2, S12在第3年投资为3, S13在第1年投资为1。
2.3 求解步骤
步骤1 确定初始种群
初始种群是种群迭代的基础, 其数量不能过少或过多, 否则容易陷入局部最优或增加计算难度。此处设定初始种群规模为200。
步骤2 确定适应度函数
通过适应度函数能够判断染色体的优劣, 本文将目标函数映射为适应度函数进行评估, 即体系总能力最大为目标函数, 目标函数值越大, 适应度越大; 反之亦然。
步骤3 选择
选择算子指从群体中选择优秀个体、淘汰劣质个体, 使得有用的遗传信息保留下来, 提高收敛性。其过程需要根据由适应函数得到的适应度值的大小来实现, 即不断选择优秀基因进行遗传, 从而完成种群迭代。本文采用较常用的轮盘赌法进行选择操作, 即个体被选择的概率与个体的适应度有关, 适应度越高的个体被选择的概率越大, 适应度越低的个体被选择的概率越小。被选择的个体进一步进行遗传操作。
步骤4 交叉
交叉算子是遗传算法的核心, 指由父代个体形成新个体的操作, 其过程需按规定选择交叉点, 以及交换交叉点的基因片段。本文采用的交叉策略为两点交叉和均匀交叉两种方式结合的交叉。
两点交叉指在父代中随机选择两个位置进行交叉, 若父代1=1-2-3-4-5-4-3-2-1-2, 父代2=5-5-5-5-4-5-4-3-2-1。假设选择的交叉点位置为2和4, 则子代1=1-5-5-5-5-4-3-2-1-2, 子代2=5-2-3-4-4-5-4-3-2-1;均匀交叉是个体中各个位置的基因都以相同的概率参与交叉, 即在每个基因位置上都进行判定, 如果随机数的大小满足概率要求, 则该位置进行交叉, 否则判定下一个交叉位置。
步骤5 变异
变异算子指替换父代染色体基因片段, 产生新的子代染色体序列。本文采用的变异方法为单点变异, 即在满足变异概率的前提下, 随机选择一个变异位置, 将其基因随机转换为1-5的值。变异操作后应计算父代、子代染色体的适应度, 若子代较优, 则变异成功, 父代由子代替换; 若父代较优, 则变异失败, 继续迭代, 直至满足终止条件。
步骤6 确定终止条件
进化次数限制、计算资源限制、最优解限制、适应度限制等均可设置为终止条件。本文设置迭代次数为10 000, 通过不断迭代达到预设的迭代次数并得到最优解时, 结束迭代, 输出迭代结果。
3 示例分析
以某个五年规划中的联合作战体系建设发展资源分配问题为例进行示例分析。
3.1 参数设置
表1 每种能力最低需求
一级体系 | 能力1 | 能力2 | 能力3 | 能力4 | 能力5 |
S1 | 3/4/4/6/7 | 0/2/4/6/11 | 1/5/9/16/16 | 2/5/8/10/15 | 3/3/9/14/17 |
S2 | 2/3/7/10/16 | 2/2/2/6/11 | 4/8/12/17/18 | 0/1/3/3/10 | 3/3/8/9/16 |
S3 | 1/4/4/8/14 | 2/7/10/13/20 | 2/6/9/12/19 | 0/4/4/7/10 | 1/4/9/15/19 |
S4 | 3/5/6/10/17 | 1/3/4/6/14 | 3/7/10/11/11 | 2/6/11/16/22 | 2/3/5/10/14 |
表2 每年每种能力加权和最低需求
年数 | 能力1 | 能力2 | 能力3 | 能力4 | 能力5 |
1 | 1.064 | 6.651 | 24.244 | 3.320 8 | 13.219 |
2 | 21.359 | 20.463 | 24.857 | 24.034 | 20.035 |
3 | 27.890 | 22.866 | 41.073 | 41.858 | 22.267 |
4 | 33.356 | 35.244 | 46.107 | 68.330 | 49.581 |
5 | 46.276 | 39.166 | 61.315 | 92.205 | 67.561 |
3.2 结果分析
设置初始种群规模为200, 迭代次数为10 000, 利用算法迭代求解模型, 输出一个五年规划内每个二级体系Sij得到发展的年份, 如表 3所示。
表3 二级体系发展年份
j | i | |||||||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | |
1 | 4 | 2 | 3 | 5 | 1 | 3 | 4 | 4 | 3 | 2 | - | - | - | - |
2 | 1 | 5 | 2 | 5 | 3 | 3 | 1 | 3 | 4 | 2 | 1 | 3 | - | - |
3 | 4 | 2 | 1 | 1 | 2 | 1 | 3 | 1 | 2 | 4 | 1 | 1 | 5 | 1 |
4 | - | 5 | 5 | 2 | 1 | 2 | 5 | 4 | 2 | 3 | 2 | 3 | 2 | 5 |
可见, 每个一级体系对应的二级体系在一个五年规划内分别发展一次, 表格中数字代表得到发展的年份。输出模型求解迭代图如图 3所示。
图3 模型求解迭代图
由图 3可知, 当种群迭代至大约5 500次后, 目标函数值保持不变, 得到帕累托最优解, 即体系总能力最大值。输出一个五年规划内每个二级体系每年的实际发展费用, 如表 4所示。
表4 二级体系发展费用
二级体系 | 第1年 | 第2年 | 第3年 | 第4年 | 第5年 |
S11 | - | - | - | 0.242 7 | - |
S12 | - | 0.263 2 | - | - | - |
S13 | - | - | 0.252 7 | - | - |
S14 | - | - | - | - | 0.239 1 |
S15 | 0.284 2 | - | - | - | - |
S16 | - | - | 0.256 6 | - | - |
S17 | - | - | - | 0.219 3 | - |
S18 | - | - | - | 0.245 4 | - |
S19 | - | - | 0.245 9 | - | - |
S1, 10 | - | 0.272 4 | - | - | - |
S21 | 0.235 1 | - | - | - | - |
S22 | - | - | - | - | 0.189 8 |
S23 | - | 0.214 6 | - | - | - |
S24 | - | - | - | - | 0.194 4 |
S25 | - | - | 0.208 9 | - | - |
S26 | - | - | 0.201 9 | - | - |
S27 | 0.202 2 | - | - | - | - |
S4, 12 | - | - | 0.184 | - | - |
S4, 13 | - | 0.197 9 | - | - | - |
S4, 14 | - | - | - | - | 0.168 9 |
可见, 设置的49个二级体系在其对应的年份内均得到发展。输出一个五年规划内每个一级体系每年的发展费用, 如表 5所示。
表5 一级体系每年的发展费用
费用 | 第1年 | 第2年 | 第3年 | 第4年 | 第5年 |
S1 | 0.284 2 | 0.535 6 | 0.755 2 | 0.704 7 | 0.239 1 |
S2 | 0.652 3 | 0.443 8 | 0.815 2 | 0.185 3 | 0.384 2 |
S3 | 1.401 2 | 0.571 6 | 0.180 9 | 0.336 8 | 0.158 9 |
S4 | 0.215 2 | 0.997 5 | 0.373 8 | 0.175 2 | 0.333 9 |
综上所述, 该示例结果表明, 在一个五年规划内, 在达到总能力最大目标的同时, 体系建设发展的总费用为9.744 6亿元, 未超出总预算10亿的约束范围, 因此通过该示例能够验证所提方法的正确性与可行性。
4 结论
资源分配是联合作战体系建设发展中的重要问题, 合理配置资源、实现资源有效利用, 对于提高军事能力意义重大。本文采用数学规划的方法, 基于体系能力对联合作战体系资源分配问题构建模型, 并运用遗传算法重复迭代、求解模型, 得到模型的帕累托最优解, 通过示例验证了所提方法能够完成五年规划内的国防资源统筹分配, 进一步支撑联合作战体系建设发展的顶层规划。
声明:公众号转载的文章及图片出于非商业性的教育和科研目的供大家参考和探讨,并不意味着支持其观点或证实其内容的真实性。版权归原作者所有,如转载稿涉及版权等问题,请立即联系我们删除。