基于马尔可夫链的装备技术成熟度评估方法论文

基于马尔可夫链的装备技术成熟度评估方法

李亮，宋敬华，郭齐胜

（陆军装甲兵学院，北京 100072）

摘要：针对传统技术成熟度评估中存在的时滞问题，基于有限状态齐次马尔可夫链模型，对装备技术成熟度、集成成熟度和系统成熟度的初始状态、演化过程和稳态分布进行研究，降低了时滞带来的决策误差，通过示例进行了演示验证，示例显示建立的方法科学有效、通用性强，具有实际意义和推广价值。

关键词：马尔可夫链，时滞，技术成熟度，集成成熟度，系统成熟度，评估方法

0 引言

评估在日常生活中扮演着重要的角色。人们经常希望通过比较不同事物的优劣以作出正确的决策。对于构成装备及装备体系的关键技术也是如此，成熟度评估是评价不同关键技术好坏的重要手段之一。

现代武器装备技术水平进行综合评价的目的，是要回答武器装备的先进程度如何，其简单的评价指标是先进、一般和落后。但是，由于武器装备技术是不断发展的，进行评价的参照系是随着时间的推移而不断变化的，也就是说，对武器装备的技术水平进行评价不能脱离时间这个因素。

成熟度评估的结果经常有时滞。因此，依据评估结果所做的决策也有时间上的滞后。本文尝试提出一种基于马氏链的成熟度评估方法以解决这个问题。将关键技术成熟度的评估值作为一个马尔可夫链，初始分布和转移概率矩阵可以通过统计计算得到。进一步可以通过计算得到成熟度的稳态分布。将成熟度的数学期望作为成熟度的评估值。基于初始分布、转移概率矩阵和稳态分布的计算可以分别得到成熟度的初始值、演化值和稳态值。可以通过成熟度评估次数的增加，不断更新转移概率矩阵的信息，从而更新成熟度的评估值以保证评估结果的准确性。

1 马尔可夫链与成熟度评估

成熟度评估存在时滞问题，针对时滞可采用马氏链方法进行解决。本节主要对成熟度评估的时滞问题和马氏链的基本理论进行分析。

“来不及了，这篇小说早就发表了，后来我还专门找来读了。该小说跟小说家跟我说过的构思差不多，你说的也不无道理，但对被盯梢者兼男主人公来说是不公平的。”

1.1 时滞问题

本文所谓的时滞有两方面含义：一方面评估结果有时滞；另一方面决策有时滞。

1）时间与空间无时无刻不在变化。基于这样的事实，对传统评估方法而言，当经过一系列复杂的过程得到评估结果时，被评价对象的状态已经发生改变。这样一来，评估结果有了时间上的滞后，即评估的时滞。

2）所谓的决策时滞是指由于评估的时滞，基于评估结果所作出的决策也不可避免地具有时滞和误差。这样一来决策也就有了时间上的滞后。

1.2 解决方案-马尔可夫链

郑英龙（1974-），男，浙江温州人，中国政法大学司法文明协同创新中心博士研究生，浙江工商大学法学院教授。

以服装产业为例，中国的服装行业和产业正发生着前所未有的巨变，社会对专业人才的知识结构和能力结构的要求越来越严格。目前，地方高校的服装设计与工程专业普遍存在着人才培养定位不准确、人才培养缺乏针对性和应用性、毕业生的工作能力不适应地方经济发展和行业需求等问题。

蓄能器Pacc和发动机泵P1的功率，通过使用功率分流因子概念从P2的预计功率确定.这个概念被选择作为模型的控制输入变量,并且可以在μ∈(-∞,1]范围内变化，则

给定σ-域B上的状态空间S，一个由在S中取值的随机变量 ξ_n组成的随机过程{ξ_n，n=0，1，2，…}称为马尔可夫链，如果对每一个非负整数n和集合T∈B，几乎处处成立

当 P{ξ_n+1∈T|ξ_n}不依赖于起始时刻 n 的取值，则称之为齐次马尔可夫链。当S是一个有限集合时，称之为有限状态齐次马尔可夫链（马氏链），式（1）称为马尔可夫性（马氏性）。

1.2.2 平稳分布（不变测度）

事物的发展，是多层次的随机事件与必然事件相互交替和相互作用的过程。虽然个别随机事件在某次试验或观察中可以出现也可以不出现，但在大量试验中它却呈现出明显的规律性-频率稳定性^［1］。下面对相关概念进行介绍^［2-4］。

称概率分布{π_i，i∈S}为马尔可夫链ξ的平稳分布，如果

式（2）可以表述为向量形式

其中，

课堂是学生求知、创新、展示自我、体验成功的平台。激励机制是学生发展的动力，在课堂中，学生时常会提出一些新颖的看法，独到的见解，对于那些发表不同想法的同学也应加以表扬，给予鼓励，提升他们的信心，给学生以无限的思维空间。一个人总是会跌倒的，但关键是要鼓励他能够迅速爬起来,这种激励不是针对一部分人，而是对所有的学生。通过导语激励，目标激励，创新激励和评价激励等，慢慢消除职高学生的自卑心理，帮助他们树立信心，学生才会感动。

1.2.3 基本假设

1.5 疗效判定指标[6] 痊愈：B超检测子宫大小、胎儿大小与孕周基本符合，BPD、FL及S/D比值都有明显改善，双侧RI均<0.8；有效：B超检测子宫大小、胎儿大小与孕周基本符合，BPD、FL、RI及S/D比值都有所改善；无效：B超检测胎儿发育不良或停止发育，BPD、FL、RI及S/D比值均无改善。

当对一个系统进行评估时，第k+1次评估的结果仅与第k次的有关，而与第k次之前的评估结果无关。换句话说，专家对每个指标的打分具有马尔可夫性。如果打分值在整数集中取值且与起始时刻无关，那么每个指标的分值组成的随机过程就是一个有限状态齐次马尔可夫链。

移动车体的结构参数确定为长650 mm，宽L3=500 mm，高985 mm，移动车体距离开关柜L3=400 mm，机械臂底座距离车体侧面L2=400 mm，开关柜高2250 mm。根据机械臂的运动学模型、移动车体及开关柜的结构参数，利用Matlab Robotics Toolbox工具箱建立机械臂、移动车体及开关柜的三维模型，如图13所示，然后基于蒙特卡洛法对机械臂在工作空间中的可操作性进行研究。分析结果如图14所示，该图为机械臂在开关柜正面区域的可操作性图，可以看出在该区域内机械臂普遍具有较高的可操作性，满足了机械臂在开关柜正面及背面检测及操作点处的运动性能要求。

在如上假设条件下，马尔可夫链的相关理论就可以用来解决涉及专家打分评估的相关问题。马氏链中初始分布的概念可以用来刻画指标得分的初始出现频率。转移概率矩阵的概念可以用来描述第k+1次评估和第k次评估之间的关系。对打分结果的统计分析恰好可以描述指标权重的当前状态；借助转移概率矩阵，又可以得到指标权重的演化趋势；最后，平稳分布的概念可以用来确定权重的稳态值。

首先，初始权重可以通过计算指标打分值的频率来获得；其次，指标权重随着时间推进的演化趋势可以通过转移概率矩阵的相关运算得到；最后，各个指标最终的稳态值可以通过平稳分布的相关计算进行确定。

根据初始权重（通常是传统方法所得到的评估结果），决策者可以像往常一样作出他们的决策；根据权重的演化趋势，决策者可以了解评估结果的变化趋势。最终，再结合稳态分步的结果，由于评估结果时滞的减小，并且系统的当前状态及演化趋势已知，决策者的决策误差也就有了相应减小的可能。

基于马氏链的评估方法过程简单、模型通用，下面将用来消除成熟评估过程的时滞，建立基于马氏链的成熟度评估方法

2 基于马尔可夫链的装备技术成熟度评估

建立基于马尔可夫链的成熟度评估方法，需要一些合理的假设和适当的准备工作。本节还将对所提出方法的具体步骤、与经典模型的比较、向技术体系成熟度评估的推广进行分析。

[2]吕峰,杨莉,董如福,杨宏.智能化医院管理信息系统建设分析[J].医学信息学杂志,2015,36(01):19-21.

2.1 基本假设

根据1.2节中的基本假设，得出如下的基本假设：

1）假设第k+1次成熟度评估的结果仅与第k次成熟度评估的结果有关，而与之前的评估结果无关。这其实是比较符合实际的，因为假设某一技术在k时刻处于TRL₂，那么在k+1时刻，正常情况下它不依赖于之前k-1时刻的成熟度等级，最起码他不会回到TRL₁。这样一来，根据相关定义，TRLs和IRLs就具有了马尔可夫性。

特别地，，有

正如恩格斯所指出的，表面上是偶然性在起作用的地方，这种偶然性始终是受内部隐蔽着的规律支配的，而问题只是在于发现这些规律。技术发展规律的分析思路是分析现状、趋势、对比找出不足。在假设1）条件下，考虑到成熟度等级的有限状态（一般为5、7或者9状态）以及它的状态转移与起始时刻无关（齐次性），有限状态马尔可夫链的相关理论，可以用来描述并解决技术成熟度评估、集成成熟度评估和系统成熟度评估的问题。

2.2 准备工作

为了用马尔可夫链的相关理论来解决成熟度评估的时滞问题，需要做一些前期的准备工作。

为了计算稳态的成熟度，首先要得到稳态分步。稳态分步可以通过求解如下方程组得到。

为了方便起见，将成熟度等级记为{m₀，m₁，…，m_n}，这样一来，n级的标度及组成了一个状态空间。

2）获取TRL/IRL的初始分布

If (A is a) and (B is b) then (C is c) and (D is d) and (E is e) and (F is f)

对第k次和第k+1次成熟度评估数据进行统计分析，初始分布可以表示为：

教学督导委员会根据调查问卷和访谈记录进行分析，综合评价各项毕业要求的达成情况，可以按达成、基本达成和未达成得出定性评价结论。

其中，P{ξ（k）=m_l}=n（m_l）/∑n（m_l）（l=0，1，…，n），即P{ξ（k）=m_l}表示的是成熟度 m_l在第 k 次成熟度评估中出现的频率。

3）进行进一步的统计分析以获取TRL/IRL的一步概率转移矩阵

计算成熟度从状态i到状态j转移的概率，作为概率转移矩阵P中的元素p_ij，P=（p_ij）_{（n+1）×（n+1）}。

完成以上工作以后，初始成熟度、演化成熟度就可以进行计算了。实际上，成熟度被认为是他们自身的数学期望。

初始成熟度计算如下：

r步转移后的演化成熟度计算如下：

1）确定技术成熟度和集成成熟度的状态空间

其中，1_n+1=（1，…，1）'；0=（0，…，0）'；P=（p_ij）_{（n+1）×（n+1）}；（·）+表示的是矩阵“·”的 Moore-Penrose广义逆矩阵（因为它的存在性和唯一性）。

这样，稳态成熟度的计算如下：

2.3 具体步骤

模型流程如图1所示。

认识你时，你跟我一样远离家乡来到深圳市，进了一家玩具厂做普工。在我生活遇到困难时，你曾大义地帮助过我。你长得文静秀美，非常爱上网。我们两个人的交情很好，下班后常在一起聊天、上网、听歌。跟你在一起我感到很快乐。

图1 基于马氏链的SRA方法流程图

1）关键技术CTEs获取

CTEs通过QFD/TRIZ/WBS进行获取。通过WBS可以获取初始的CTEs，借助QFD/TRIZ的帮助，可以获取CTEs的权重^［5-6］。综合考虑可以得到系统的关键技术指标CTEs。

2）初始IRL/TRL获取和初始SRL的计算

邀请若干评估小组分别进行成熟度评估。记录初始的TRL/IRL及相应的SRL。为了进行评估的方便，设计如图2（表示CTE_i的技术成熟度是TRL_i，CTE_i和CTE_j的集成成熟度是IRL_ij）所示的表示方式。

图2 TRL/IRL关系示意图

假设有n个CTEs，TRL=（TRL_i）_n×1，IRL=（IRL_ij）_n×n，SRL=（SRL）_1×1，那么 TRL/IRL/SRL 的关系可以用公式表示如下：

2) 额定转速工况下，该传动腔室中的No1轴承发热量最大，达到315.6W，工作中需要特别注意对它的润滑冷却。

2）某项关键技术的发展及其与其他关键技术之间的交互存在可循的规律。

其中，l_i（i=1，2，…，n）表示的是和CTE_i有集成关系的CTEs的个数（包括CTE_i本身在内）；ω_i（i=1，2，…，n）表示的是CTE_i的权重。

根据式（5）和TRL/IRL/SRL的初始状态可以获取：。

3）进行第k次成熟度评估

邀请相同数目（尽量和之前相同）的评估小组进行第k次成熟度评估。将评估数据记为：和。

1.2.1 有限状态齐次马尔可夫链

4）进行第k+1次成熟度评估

邀请相同数目（尽量和之前相同）的评估小组进行第k+1次成熟度评估。将评估数据记为：TRL^k+1=，和。

5）获取一步转移概率矩阵

通过对各成熟度的转移频率进行分析，分别获取TRL和IRL的一步转移概率矩阵。

6）求解稳态分步

通过求解式（3），可以得到TRL和IRL的稳态分布TRLs和IRLs。

7）稳态TRL/IRL获取与SRL计算

根据式（4），计算稳态 SRL。

8）结果对比分析并作为决策支持

将初始的、演化的和稳态的SRL进行对比，得出的相关分析结果可以作为决策支持。

2.4 和经典模型的比较

1）从成熟度获取方面。基于马氏链的系统成熟度分析模型通过初始分布、转移概率矩阵和稳态分步获取3种形式（初始、演化、稳态）的成熟度；经典模型仅能通过简单的统计分析或数学计算获取一个成熟度数值。

2）对时滞的处理方面。基于马氏链的模型通过转移概率矩阵和稳态分布来解决时滞问题；经典模型对时滞问题基本没有考虑。

3）从短板的避免方面。对于每一项CTE在装备发展的不同阶段通过不同的阈值（短板的底线）来进行约束，从而避免单项技术短板的出现。同时，由于IRL的引入，还可以克服集成过程中可能出现的“集成短板”，这一点是传统方法不容易实现的。

综合来讲，由于马氏链模型的采用从某种程度上可以降低系统设计中的风险，决策者的决策变得更加准确和科学。

2.5 向技术体系成熟度评估的推广

引入系统集成成熟度（System Integration Readiness Level，SIRL）的概念，用于衡量不同系统（System）之间集成成熟度，并用于计算技术体系成熟度（Technical System of System Readiness Level，TSoSRL）。

其中，h_i（i=1，2，…，n）表示的是和System_i有集成关系的 System的个数（包括System_i本身在内）；μ_i（i=1，2，…，n）表示的是 System_i的权重；m 表示的是系统集成成熟度的分级；SRL/1中的1是因为SRL是［0，1］之间的一个数。

通过类似技术成熟度到系统成熟度的步骤，就可以通过系统成熟度和系统集成成熟度获取技术体系成熟度。

3 应用示例

为了验证基于马氏链的方法的科学性和可行性，下面给出一个基于马氏链的系统成熟度评估的例子。

1）通过WBS/QFD/TRIZ过程，得到4项关键的装备技术。第k次和第k+1次评估结果（评估小组由4个专家组成）分别如图3和下页图4所示。

将六个沉积物的峰度K的平均值从大到小排列：冰碛物(4.37)、湖泊(1.67)、河流沉积物(1.08)、汉源九襄地区沉积物(1.01)、海洋沉积物(0.97)、黄土(0.87)。根据Folk & Ward(1957)对K的分级规定，黄土峰态相对平坦，然后为海洋沉积物，再是汉源九襄地区沉积物和湖泊沉积物，冰碛物峰度值最高，峰态最尖锐。需要指出的是：汉源九襄沉积物峰态较冰川沉积物平坦，K值相差较大，可能是由于冰川融水搬运作用扰动了原位沉积物的粒度特征，使其峰态降低(见图9)。

图3 第k次评估结果

为了计算的方便，将图3和图4转化为表格形式，如下页表1～表4所示。

2）根据马氏链相关公式进行数据分析，将一些关键结果列在下面。

图4 第k+1次评估结果

表1 第k次评估的技术成熟度

表2 第k次评估的集成成熟度

表3 第k+1次评估的技术成熟度

表4 第k+1次评估的集成成熟度

为了更加形象地展示成熟度在不同转移步数时的状态，分别将技术成熟度、集成成熟度和系统成熟度转移步数的变化规律绘制成散点图，如下页图5～图7所示。

图5 TRL演化规律

其中，☆为 TRL₁；×为 TRL₂；▽为 TRL₃为TRL₄。

图6 IRL演化规律

其中，☆为 IRL₁₂；×为 IRL₁₃；▽为 IRL₁₄为 IRL₂₃；◇为 IRL₂₄；+ 为 IRL₃₄。

图7 SRL变化规律

从图5～图7中可以看出，在初始阶段，TRL₁TRL₂TRL₃=TRL₄并且有SRL=0.555 7；在稳态阶段，TRL₁=TRL₂=TRL₃=TRL₄并且有SRL=0.914 3。从中可以看出，如果仅仅按照传统方法（即初始状态成熟度SRL-1）进行决策，则可能会拒绝这一技术方案；但是如果参考稳态的成熟度数值（SRL-5）及到达稳态所需的转移步数（本例大约9步即可到达0.901 2，如图7所示），则有很大可能会接受这一方案。拒绝一个可行的方案将是十分可惜的。基于马氏链的评估方法减少了时滞对决策的影响，从而降低了决策的风险。以上示例证明，方法是科学可行的。

4 结论

本文针对传统成熟度评估中存在的评估结果时间滞后的问题，采用马尔可夫链的相关理论分别对成熟度的初始值、演化值和稳态值进行分析和计算，并通过示例对方法进行了演示验证。结果表明，所提出的基于马氏链的成熟度评估方法正确、合理、可行，且通用性强，可推广到其他领域中。

参考文献：

［1］李贤平.概率论基础［M］.3版.北京：高等教育出版社，2010.

［2］胡奇英，刘建庸.马尔可夫决策过程引论［M］.西安：西安电子科技大学出版社，2000.

［3］STROOCK D.An introduction to markov process［M］.Berlin Heidelberg：Springer-Verlag，2005.

［4］王梓坤.概率论基础及其应用［M］.3版.北京：北京师范大学出版社，2007.

［5］李亮，郭齐胜，李永，等.基于灰关联分析的质量功能配置方法研究［J］.计算机集成制造系统，2007，13（12）：2469-2472，2486.

［6］牛建钊，耿俊豹，魏曙寰，等.指标融合和隐马尔可夫模型的舰船装备技术状态评估［J］.火力与指挥控制，2016，41（11）：85-89.

Research on Markov Chain-based Weapon Equipment Technical Readiness Assessment Method

LI Liang，SONG Jing-hua，GUO Qi-sheng
（Army Academy of Armored Forces，Beijing 100072，China）

Abstract： In order to solve the time-lag problem in the classical Technical Readiness Assessment（TRA），a new method based on Markov Chain was advanced，including Technical Readiness Assessment （TRA）method，Integral Readiness Assessment （IRA）method and System Readiness Assessment （SRA） method.Through the research on the initial status，evolutionary procedure and steady-state distribution，the decision error caused by time-lags is reduced.The method is demonstrated through an illustrative example.The example shows that the method advanced is scientific，effective，versatility.It is meaningful and generalized.

Key words： markov chain，time-lags，technical readiness，integral readiness，system readiness，assessment method

中图分类号： TJ810.3+1

文献标识码： A

DOI： 10.3969/j.issn.1002-0640.2019.02.022

引用格式：李亮，宋敬华，郭齐胜.基于马尔可夫链的装备技术成熟度评估方法研究［J］.火力与指挥控制，2019，44（2）：108-113.

文章编号： 1002-0640（2019）02-0108-06

收稿日期： 2017-12-10

修回日期： 2018-03-09

作者简介：李亮（1982- ），男，河南荥阳人，博士。研究方向：装备需求论证和试验评估。

Citation format： LI L，SONG J H，GUO Q S.Research on markov chain-based weapon equipment technical readiness assessment method［J］.Fire Control&Command Control，2019，44（2）：108-113.

标签：马尔可夫链论文; 时滞论文; 技术成熟度论文; 集成成熟度论文; 系统成熟度论文; 评估方法论文; 陆军装甲兵学院论文;