基于Bayes混合先验分布的兵棋推演规则建模方法*
薛 辉 1，2，刘铁林1，苏小波1
（1.陆军工程大学石家庄校区，石家庄 050003；2.空军石家庄飞行学院，石家庄 050071）
摘 要： 针对兵棋推演胜负裁决规则建立难的问题，军事专家利用经验数据总结战争胜负规律，数学家、工程师利用军事运筹分析方法计算战争胜负规律，为了综合利用两种方法，利用Bayes混合先验分布将历史战争规律总结法与军事运筹分析法合理有效融合，对胜负裁决规则进行建模研究；通过在贝塔分布中引入继承因子，将兰彻斯特方程生成的军事运筹分析数据作为先验信息，与历史战争的经验规律数据作为后验信息进行结合。引入的混合贝塔先验分布方法，有效避免了两种方法互相淹没信息的问题；建立的胜负概率裁决模型，有效地提高了兵棋推演中胜负裁决概率的准确度及可信度，使其更加符合战争规律。
关键词： Bayes方法，混合Beta分布，兵棋推演，胜负概率，准确度
0 引言
战争作为复杂巨系统，影响其胜负的因素很多，武器装备先进程度、人员能力素质水平、士气、指挥员指挥艺术、作战样式、国家的政治、经济等等，尤其是信息化战争的出现，使得战争更加扑朔迷离。兵棋推演作为战争的一种模拟工具，近年来备受军事领域推崇，兵棋是通过对历史的更深理解，尝试推断未来，是历史和科学的结合，是智慧的较量，是纸做的时间机器，它是分析战争的重要手段，是一种研究战争的利器［1］。历史指运用数理统计的方法对实战、演习等数据及经验的总结，科学即运用科学的方法和理论探索战争的基本规律。胜负裁决规则作为兵棋推演的核心内容决定着模拟效果的好坏，决定着兵棋推演的结果是否符合战争规律。
当前兵棋系统的裁决规则主要通过对战争历史数据广泛采集和统计分析，通过经验数据寻找战争规律，该方法最易于被军事专家及指挥人员接受。而数学家、工程师最为推崇的是作战运筹分析，通过对武器装备战技指标的分析，运用数学方法建立作战模型，通过数学分析得出相应的结果进而揭示战争的一般规律，该方法在“二战”时的工业化背景中取得了明显成效。通过总结不难发现，当影响战争胜负的因素比较单一且易于量化时，适合运筹分析，当影响战争胜负的因素比较多元且仅能定性分析时，适合经验规律总结。日益先进的武器装备作为内因和经济全球化、精细化分工作为外因，决定着未来战争样式必然复杂多变，影响战争胜负因素必然多元化、复杂化且量化评估难。
2.2 信度与效度分析 Cronbach’s 系数为0.873，显示该测试可靠；KMO值为0.866，Bartlett球形检验卡方值为2 968.74(υ=45，P<0.001)，说明本次问卷的内在信度和结构效度均符合一般要求。
本文的研究思路是运用运筹分析的方法建立兵棋推演胜负裁决的基础作战模型，通过Bayes方法及混合Beta分布，使其与历史战争经验总结的规律进行合理有效的融合，对兵棋推演中胜负裁决的概率进行调整，进而有效提高胜负裁决概率的准确度。
贝叶斯（Bayes）方法可充分有效利用先验信息，通过对先验信息的收集、挖掘和加工，使之形成先验分布而参加到统计推断中来，以提高统计推断的准确度。贝塔分布（Beta Distribution）作为伯努力分布和二项分布的共轭先验分布的密度函数，在机器学习和数理统计学中都有重要应用，同样胜负型事件的先验分布也可为贝塔分布。
1 兵棋推演胜负裁决规则建模方法
如果仅由人员、武器装备影响战争胜负，可以根据人员训练水平、武器装备战技指标确定战斗力指数后，按照军事运筹分析方法，建立数学模型分析战争胜负规律。
在模拟机械化步兵部队在对阵地防御之敌机械化步兵部队的进攻作战行动中，设定以下辅助规则：1）进攻一方的剩余战斗力低于初始战斗力的2/3时即认为进攻失败；2）防守一方的剩余战斗力低于初始战斗力的2/5时即认为防守失败。该规则的制定虽然符合大多数战斗规律，但事实上战争中还有很多意外的情况会发生，敌对双方的胜负不一定是按照上述约定进行，其胜负规律具有一定的随机性，下面假设红方为进攻方，蓝方为防守方，将达到此临界值之前的胜负概率设置如下。
1.1 依据战术基本原则建立胜负裁决规则模型
原假设H0：X与Y是相同的总体，检验统计量为
原料经处理后，将根霉发酵剂加入已经被加水冷却至30 ℃的原料中，先酿成酒酿，再加水(1∶4)，保温28～33 ℃。拌麸皮、谷糠，接种醋酸菌种制成醋醅，发酵成熟后，经过淋醋灭菌制成成品。
农民捕捉市场信息的能力不足，农产品增产不增收、市场决定生产、生产服务市场，这是社会主义市场经济的基本规律。农业生产经营活动只有通过市场机制来组织和调配生产要素，才能得以持续健康地发展。农民种什么、养什么，不是政府确定的，也不是农民随心所欲确定的，而是由国内外市场供求状况决定的。近年来，一些地方在结构调整中不研究市场规律，而是搞“政绩工程”、“面子工程”，造成各种“大战”，如红麻大战、柑桔大战、苹果大战、蚕茧大战等此起彼伏，造成农业生产能力的大起大落和农民投资的巨大损失[5]。
求解以上一阶常系数线性微分方程组可得：
混合Beta验前分布的使用是对经典统计方法和传统Bayes方法的合理折衷。如果存在多种相关信息可以利用，如战争的样式、力量规模、武器装备性能相近等，可以构造如下的混合验前分布，不同的ρi表示历史战争经验规律总结和军事运筹分析在不同方面的相似性和继承性。
为便于计算，当假定一方剩余战斗力值进行分析时，知道 x（t）或 y（t）可利用以下公式计算另一方的剩余战斗力值
依据战术基本原则和战争规律，在双方剩余兵力不同阶段建立相应的胜负裁决概率模型：
1）当 x（t）/x0>2/3 且 y（t）/y0>2/5 时
从军事战术基本原则分析，相比于对方战斗力损耗越慢取胜的把握越大，即剩余战斗力比值大的一方胜率要大于剩余战斗力值的一方，为使双方胜负概率充满整个概率空间，通过各方剩余战斗力值与初始值的比值 y（t）/y0与 x（t）/x0进行归一化处理后，作为该方在相应t时刻的胜率。具体如下所示：
1）投用了装置1.0 MPa蒸汽专线，在4月12日停工第二天，配合动力厂投用进装置蒸汽专线，确保装置吹扫蒸汽供应充足。
2）当 x（t）/x0≤2/3 或者 y（t）/y0≤2/5 时
他是喜欢易非的，打心眼里喜欢，可是，喜欢又有多爱呢？爱，是要很多很多的喜欢的，可他和易非，还没有上升到那个程度。
按照战术基本原则，任何一方先达到临界值，即判定进攻或防守失败，暂时不考虑双方同时达到临界值的可能性。
1.2 模拟双方士气建立胜负裁决规则模型
综合考虑双方人员、装备水平及训练情况后，用α，β表示双方士气，其取值范围为［0，1］，当取0值时即认为“不战而逃”，取1值时认为“誓死战斗到底”，则将双方的胜负概率进行如下调整：
按照以上模型，兵棋推演中胜负裁决规则比较接近战争规律，这与孙子兵法中的“车杂而乘之，卒善而养之，是谓胜敌而益强”所描述的战争规律是一致的。即战争规律是双方的胜负概率随着时间的推进在不断地发生变化，且胜率大者愈大，胜率小者愈小。
许钧是我国翻译界领航的翻译家、翻译实践者和翻译理论家。他做了许多翻译实践，“在翻译的同时，还关注相关的理论问题，提出问题，也自己解决问题，从翻译实践和翻译理论两个方面做出自己的贡献”[1]他从翻译研究的角度出发，结合自己的翻译实践进行深人的思考，为推动翻译学的学科建设做出了巨大贡献，开拓了翻译研究新的视野。
2 基于Bayes混合先验的胜负裁决规则建模方法研究
在不考虑气候环境、后勤供应保障、指挥艺术等等且武器装备战斗力占主要决定因素的战斗时，上节所建立的模型可以最大程度反映战争规律，但如果影响因素多元且定量化分析难时，需要借鉴历史战争经验总结的规律，并将两种方法合理结合，提高战争模拟的可信度。
用x0表示红方初始战斗力值，x（t）表示对抗中红方t时刻的战斗力值；y0表示蓝方初始战斗力值，y（t）表示对抗中蓝方t时刻的战斗力值，a表示蓝方单个作战单元对红方作战单元在单位时间内的损耗系数，b表示红方单个作战单元对蓝方作战单元在单位时间内的损耗系数，我们是通过兰彻斯特方程［2］求解相应时刻红蓝对抗双方的剩余战斗力，如下所示：
当前国内兵棋系统采取的胜负裁决规则模型都是在经典频率学派的理论指导下得出的，即历史经验、演习数据等的统计分析，它是对所有影响战争胜负因素聚合后的分析，是整体涌现的结果。本文将其作为后验样本信息，将上节建立的模型作为先验信息，即假设上节的模型抓住了影响战争胜负的主要因素，但要通过后验样本信息对该信念进行不断调整，如果后验样本信息支持则增强信念，如果后验样本信息反对则降低信念。
用规范的语言描述，即贝叶斯学派注重于参数空间，描述的是观察者知识状态在新的观测发生后如何更新，即观察者持有某个前置信念（prior belief），通过观测获得统计证据（evidence），通过满足一定条件的逻辑一致推断得出的关于该陈述的合理性，从而得出后置信念（posterior belief）来最好地表征观测后的知识状态（state of knowledge）。严谨的数学定义如下所示：
ρ作为继承因子，反映了历史战争经验规律总结和军事运筹分析方面的相似程度，可以用统计数据或军事指挥员给出；为了体现军事指挥艺术与科学理论方法的完美结合，继承因子ρ可以由军事指挥员给出。美军的兵棋推演系统将自由式（依靠军事指挥员经验和指挥艺术）兵棋与严格式（依靠数学模型及算法）兵棋巧妙地结合进行研讨式兵棋的开发，与本文的混合Beta先验分布的思路不谋而合。由专家给出ρ的方法有时很难做到准确，文献［10］利用两样本的卡方拟合优度检验，给出了确定ρ的方法。（1-ρ）称为更新因子，反映了兵棋系统在继承历史战争经验时引入的差异性、创新性。
2.流行情况。鲤鱼、鲫鱼、鲢鱼、鳙鱼、金鱼、鳊鱼等多种鱼类从苗种到成鱼都会感染此病，尤以草鱼、青鱼最容易感染。对草鱼、青鱼的苗种危害最为严重。水温5℃以上就开始发病，鱼密度越大、水温越高、水质越差，越容易流行暴发此病。
贝塔分布中的参数可以理解为伪计数，伯努力分布的似然函数可以表示为一次事件发生的概率，它与贝塔分布有同样的形式，因此，可以用贝塔分布作为其先验分布。
政策性农业保险费用是由中央和地方财政统一承担，中央占40%，省级占20%，市县占20%，农户占20%。如果财政部门补贴资金由于当年的财政发生问题而没有发放，保险公司对农户的理赔将无法顺利进行，这会对农户正常利益造成影响。甘肃省属于经济欠发达地区，经济建设发展速度过慢。由于近几年各种农业保险险种的增加和范围的扩大，很容易导致财政预算过大、无法及时发放的问题，不但会影响保险公司进行理赔，也会对地区创新保险险种造成影响。如果政府在当年财政预算中对于农业保险的补助相应减少，也会影响保险的发展。
式中，（0<θ<1，α>0，β>0），其中 α，β 为先验分布的超参数。
2.1 不考虑差异性的胜负裁决规则建模
用SH表示红方在t时刻的胜率的期望，战争的胜负与成败型产品类型，其先验分布为贝塔分布Beta（α，β），则
“17栋老林家小子。老魏他们有意思，一到高考，就跟矿里学生家长鼓动，一定要报矿治、化工专业，说陵矿子弟，都有责任擎起父辈的旗帜，都要子承父业，女娃一样，要替父从军。”
兵棋推演系统中规定其胜算概率越高，则其方差就越小。k为调整系数取值0.01~0.05，k值越小其散布越小。
有了先验分布，结合历史战争的统计数据（n，x），n为总数，x为红方胜利数，f=n-x为红方失败数，可以得到后验分布 Beta（θ|a+x，b+f）密度函数为
这样就调整更新了兵棋推演中胜负裁决概率θ的分布，使其更加准确，更加合理。事实上这么做就是将历史战争规律的总结融入到军事运筹分析中，是军事专家与工程师、数学专家所各自推崇战争模拟方法的融合，但是两者之间各自所占比重是必须考虑的问题，主要有两点原因：一是样本量的不同必然导致先验信息与后验信息的互相淹没；二是需要科学确定两者对整体规律的贡献权重。
2.2 基于混合Beta分布胜负裁决规则建模
2.1 中建立的数学模型调整更新了胜负裁决概率，但是如果先验信息过多就会淹没现场信息，且利用上式的Beta分布进行Bayes胜负裁决概率调整，只考虑了兵棋对历史战争经验的继承性，而忽略了兵棋能够推演模拟未来新战争样式的独特方面，即任何兵棋系统都会有其差异性、创新性，可以采用均匀分布来描述兵棋推演的差异性、创新性，因此，采用混合Beta分布来描述兵棋推演胜负裁决规则的继承和引入的差异性、创新性。文献［4-9］解决成败型产品可靠性评估时均采用了混合Beta先验分布。
“十三五”全国油气资源评价全面展开。完成了致密油气、页岩气、煤层气等关键评价参数研究。经济可采性、生态环境允许程度评价方法研究取得重要进展。初步完成了海域深水(300米以深)和银额盆地油气资源潜力评价。
借助德国卡尔蔡司偏光显微镜观察C/C试样的金相结构.在MM-W1型立式万能摩擦磨损试验机上进行摩擦磨损试验(见图1).C/C试样的对磨销为40Gr钢,硬度HRC45,规格为Φ5.5 mm×10 mm,表面粗糙度Ra0.8µm.
定义1［3］在获得样本x后，参数θ的后验分布（posterior distribution）就是在给定X=x条件下θ的条件分布，记为π（θ|x）。对于有密度的情形，它的密度函数为
设有历史样本信息（m，x）来自总体X，m为试验数，x为成功数，f=m-x为失败数；样本信息（n，y）来自总体Y，n为试验数，y为成功数，f'=n-y为失败数。
经典兰彻斯特方程假定，对抗双方均处于敌方火力范围内，由于战斗紧迫可以不考虑支援部队，双方的战斗力完全取决于双方的作战单元数和单个作战单元的战斗力值。
K依分布收敛到自由度为1的卡方分布。给定检验水平α，有
式中，χ12（α）为自由度为1的卡方分布的α分位点。式（16）中要求，x，y，f，f'都应大于 5，若不能满足这一条件，Yates［11］给出了对K的一个修正：
以下使用的K均为式（18）中的定义，式（18）中的K近似服从自由度为1的卡方分布。上述检验问题中，一个略小于χ12（α）的K和一个远小于χ12（α）的K，意义有所不同，后者支持原假设的理由更强烈一些，令
称为该检验的拟合优度，Q（K）越大，支持原假设的证据就越强。文献［10］指出从Q（K）的概念来看，它可以作为两总体的相似程度的度量，那么，ρ和Q（K）是具有相互联系和影响的，记此种关系为
f取决于现场样本信息和历史样本信息的情况，可以通过在各种可能的概率下产生现场样本和历史样本相同容量的随机数来模拟计算，大多数情况下可以令
通过变换可得到如下兰彻斯特平方律方程
根据贝叶斯公式，结合混合先验分布得到的后验概率模型为：
使用兵棋推演的目的是检验及调整作战计划，提高打赢概率。从“预己从严、料敌从宽”的角度讲，制定的作战计划应该确保有最低的胜算，因此，军事指挥员更关心的是给定可信度（置信度）γ后，最低的胜率值SHL，以便定下决心是否执行该作战方案。在给定可信度γ后，SHL可由下式解出：
3 仿真算例分析
3.1 参数假设
假设按照军事运筹分析的方法，通过第1节的方法，运用兰彻斯特方程计算得出红方胜率的期望值SH=0.8，即红方胜率为0.8，蓝方胜率为0.2。
多年的实践表明，流域和区域“一湖两河”水行政执法联合巡查的持续开展，为水行政执法工作提供了合作平台，有效整合了流域与区域执法力量，丰富了执法手段和方式，发挥了流域与区域两方面的积极性，有利于及时、集中查处省际边界违法行为和复杂疑难案件，保障了边界地区与重点河湖和谐稳定，取得了流域与区域优势互补、事半功倍的效果；同时促进了各级队伍的日常巡查活动，各级水行政执法队伍密切配合，营造了声势，扩大了影响，树立了水行政执法队伍的威信。
通过对历史相似战争案例广泛采集和统计分析后，红蓝双方的胜负情况分别为：HS=4，LS=1。在给定置信水平γ=0.8的情况下，采用经典的频率统计方法对战争规律进行总结，取
红方的胜率为θ=0.51。
3.2 计算分析
按照军事运筹分析的方法，无论实战规律是否与数学模型推演结果一致，兵棋推演系统仍然没有必要修改红方胜率SH=0.8。同样，按照对战争经验总结，兵棋推演系统仍然没有必要修改红方胜率SH=0.51。
按照第2节建立的模型，基于实战规律作为小样本，在有先验信息的情况下，按照基于混合Beta分布的胜负裁决概率模型，假设调整系数k=0.02，利用式（13），可计算得到贝塔分布超参数的拟合值α=31.2，β=7.8。
两种样本信息之间的继承因子可以由卡方拟合优度检验方法求得，也可由专家给出。假如由军事指挥员根据战争经验给出继承因子ρ=0.7，则在给定的置信水平γ=0.8的情况下，根据式（23）和式（24）得出红方的胜率为0.751 2。
通过以上对比，折衷后的胜负概率值较前面两者更为稳健，通过卡方拟合优度检验法确定的继承因子，充分体现了历史战争经验规律总结法和军事运筹分析法两者之间的继承性以及差异性，使得兵棋推演包含两者共同优点，有效地提高了模拟程度。
由式子可以看出，在l1与l2一定及在同一个人的情况下臀部运动轨迹与θ1、θ2的有关.所以要首先研究人坐起过程中θ1、θ2的变化过程.
4 结论
利用混合贝塔分布的贝叶斯方法将历史战争经验与军事运筹分析方法结合，对兵棋推演的胜负裁决规则建模方法暂没有相关文献进行研究，本文尚属首次提出，通过算例分析，证明该方法确实有效，得出的胜负裁决概率合理可靠，能够有效促进兵棋推演系统更好地模拟战争规律，对于提高战斗力、增强打赢本领具有一定的军事意义。下一步可以通过Bayes网络建模，使用实时更新的数据库、知识库不断训练参数，动态调整胜负裁决规则。
表1 3种方法的比较（置信度γ=0.7）
参考文献：
［1］DUNIGAN J F.The complete wargames handbook［M］.New York：William Morrow﹠Company，1997.
［2］陈向勇，井元伟，李春吉，等.一种基于Lanchester方程的交战取胜最优策略 ［J］. 控制与决策，2011，26（6）：945-948.
［3］韦来生，张伟平.贝叶斯分析［M］.合肥：中国科学技术大学出版社，2013.
［4］张士峰，樊树江，张金槐.成败型产品可靠性的Bayes评估［J］.兵工学报，2001，22（2）：238-240.
［5］董海平，蔡瑞娇，王玮.基于混合Beta分布的火箭弹射座椅可靠性评估［J］.航空学报，2009，30（2）：232-235.
［6］张文杰，杨华波，张士峰.基于Bayes混合验前分布的成败型产品可靠性评估 ［J］. 兵工学报，2016，37（3）：505-511.
［7］明志茂，陶俊勇，陈循，等.基于混合Beta分布的成败型产品Bayes可靠性鉴定试验方案研究［J］.兵工学报，2008，29（2）：204-207.
［8］杨华波，蔡洪，张士峰.基于混合验前分布的制导精度评定方法［J］.航空学报，2009，30（5）：855-860.
［9］田福平，汶博，郑鹏鹏.基于贝叶斯网络的作战目标评估［J］.火力与指挥控制，2017，42（2）：79-82.
［10］王玮，周海云，尹国举.使用混合Beta分布的Bayes方法［J］.系统工程理论与实践，2005，25（9）：142-144.
［11］FISHER R A，YATES F.Statistical tables for biological agriculturaland medicalresearch ［M］.Edinburgh：OLiver﹠Boyd，1957.
Research on Modeling Methods for Wargames Deduction Rules Based on Bayesian Hybrid Prior Distribution
XUE Hui1，2，LIU Tie-lin1，SU Xiao-bo1
（1.Shijiazhuang Campus of Army Engineering University，Shijiazhuang 050003，China；2.Shijiazhuang Flying College of the PLA Air Force，Shijiazhuang 050071，China）
Abstract： To solve the problem of Wargames which is difficult to establish outcome adjudication rule，military experts use empirical data to summarize the laws of war，mathematicians and engineers use the military operation analysis to calculate the laws of war.The bayesian mixed prior distribution is used to combine the historical war rules and the military operational analysis method reasonable and effectively for studying modeling methods of outcome adjudication rule.The paper merges the different means by introducing inheritance factor in the beta distribution，which takes military operation analysis data that established by Lanchester equation as prior information，the other takes historical war experience as posterior information.The introduction of the mixed beta prior distribution methods can avoid the problem of that prior data inundate the posterior data effectively.This paper established the model of outcome adjudication rule，which can not only effectively improve the accuracy but also enhanced the credibility.what’s more，which makes wargames deduction more in line with the laws of war.
Key words： bayesian method，mixed beta distribution，wargames deduction，outcome probability，accuracy
中图分类号 ：TJ01；E919
文献标识码： A
doi： 10.3969/j.issn.1002-0640.2019.06.021
引用格式： 薛辉，刘铁林，苏小波.基于Bayes混合先验分布的兵棋推演规则建模方法［J］.火力与指挥控制，2019，44（6）：106-110.
文章编号 ：1002-0640（2019）06-0106-05
收稿日期： 2018-05-05 修回日期：2018-06-07
*基金项目： 国家社会科学基金军事学基金资助项目（15GJ003-185）
作者简介 ：
薛 辉（1987- ），男，河北鹿泉人，博士研究生。研究方向：装备保障指挥理论与应用。
刘铁林（1971- ），男，河北满城人，教授，博士生导师。研究方向：军事装备学。
Citation format： XUE H，LIU T L，SU X B.Research on modeling methods for wargames deduction rules based on bayesian hybrid prior distribution ［J］.Fire Control﹠Command Control，2019，44（6）：106-110.
标签：Bayes方法论文; 混合Beta分布论文; 兵棋推演论文; 胜负概率论文; 准确度论文; 陆军工程大学石家庄校区论文; 空军石家庄飞行学院论文;