CN109241583B - Human-computer interaction system reliability solving method based on Markov - Google Patents

Human-computer interaction system reliability solving method based on Markov Download PDF

Info

Publication number
CN109241583B
CN109241583B CN201810941561.XA CN201810941561A CN109241583B CN 109241583 B CN109241583 B CN 109241583B CN 201810941561 A CN201810941561 A CN 201810941561A CN 109241583 B CN109241583 B CN 109241583B
Authority
CN
China
Prior art keywords
state
fault
human
states
event
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810941561.XA
Other languages
Chinese (zh)
Other versions
CN109241583A (en
Inventor
尤启东
曾声奎
郭健彬
吕红红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201810941561.XA priority Critical patent/CN109241583B/en
Publication of CN109241583A publication Critical patent/CN109241583A/en
Application granted granted Critical
Publication of CN109241583B publication Critical patent/CN109241583B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Test And Diagnosis Of Digital Computers (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a human-computer interaction system reliability solving method based on Markov, which comprises the following steps: the method comprises the following steps: analyzing a human-computer interaction system, an environment and a task scene, and analyzing fault logic to determine a bottom event set causing a human-computer system fault; step two: determining a complete set of state space, and drawing a system state transition diagram according to system fault logic; step three: calculating to obtain the transfer rate between states, and constructing a transfer rate matrix; step four: listing a state equation, and solving to obtain the instantaneous reliability of the human-computer system; through the steps, the fault logic in the human-computer interaction process is analyzed, the set of system fault bottom events is determined, and the instantaneous reliability of the human-computer system is finally solved, so that the effects of accurately solving the reliability of the human-computer system and searching system defects are achieved, and the problem that the multi-factor coupling characteristic of the human-computer system is difficult to describe in the conventional human-computer system reliability modeling is solved.

Description

基于马尔科夫的人机交互系统可靠度求解方法A Markov-based Reliability Solution for Human-Computer Interaction System

技术领域technical field

本发明提供一种基于马尔科夫的人机交互系统可靠度求解方法,能够深入考虑人的认知特性、场景任务特性和人机耦合特性,属于人机交互可靠性定量建模与预计领域。The invention provides a Markov-based reliability solution method for a human-computer interaction system, which can deeply consider human cognitive characteristics, scene task characteristics and human-computer coupling characteristics, and belongs to the field of human-computer interaction reliability quantitative modeling and prediction.

背景技术Background technique

人机交互建模的目的是在充分考虑存在环境扰动、系统故障和人为失误的情况下,将人、机、环作为整体,对任务需求下的人机交互动态过程进行描述,研究人机之间的信息传递,分析可能导致事故的潜在风险场景。人机交互建模依据实现方式进行分类,可以分为基于逻辑的方法和基于仿真的方法两类。第一类方法是通过逻辑分析对系统内的人机环要素以及相互之间的影响关系进行建模,得到元素组合或序列;第二类方法是通过仿真技术来实际模拟人机交互过程,将环境扰动、系统故障和人为失误加入人机系统,从而分析异常态下的人机交互行为。The purpose of human-computer interaction modeling is to describe the dynamic process of human-computer interaction under task requirements by taking human, machine, and environment as a whole while fully considering the existence of environmental disturbances, system failures and human errors. information transfer between, and analysis of potential risk scenarios that may lead to accidents. Human-computer interaction modeling is classified according to the implementation, and can be divided into two categories: logic-based methods and simulation-based methods. The first type of method is to model the elements of the human-machine loop in the system and the relationship between them through logical analysis to obtain element combinations or sequences; the second type of method is to use simulation technology to actually simulate the human-computer interaction process. Environmental disturbances, system failures and human errors are added to the human-machine system to analyze the human-machine interaction behavior under abnormal conditions.

在逻辑建模方面,美国贝尔实验室于1961年首次提出故障树分析(Fault TreeAnalysis,FTA)方法,Dugan引入了功能相关门、优先与门、备件门等动态逻辑门,由此形成了动态故障树分析(Dynamic Fault Tree Analysis,DFTA)方法,并提出使用马尔科夫模型来进行定量分析,众多学者对其不足进行了改进,得到了广泛引用。但该类方法仅将人误视为底事件而没有深入考虑人的认知特性、场景任务特性和人机耦合特性。事件树分析(Event Tree Analysis,ETA)方法是一种逻辑演绎的,自下而上的分析方法,能从因果逻辑关系对人机交互过程进行描述,建模过程简单,但是仅能以线性链条的方式进行描述而无法对组合事件(人机环)进行分析,描述能力不足。在仿真建模方面,马里兰大学风险与可靠性中心开发了班组信息、决策和行为响应(Information,Decision,and Action in CrewContext,IDAC)仿真分析系统;塔菲研究所、麻省理工大学、芝加哥大学、英国伯明翰大学等国外研究机构开展了多主体(Multi Agent,MA)仿真技术研究并在开发相应平台。基于仿真的建模方法需要针对不同任务场景中的多种认知失误模型和故障机理模型分别进行仿真模型的构建,而遍历事件组合或时序十分耗时,因而工作量浩大且通用性差,同时仿真结果难以进行校核。In terms of logic modeling, Bell Labs first proposed the Fault Tree Analysis (FTA) method in 1961. Dugan introduced dynamic logic gates such as function-related gates, priority AND gates, and spare parts gates, thus forming dynamic faults. Tree analysis (Dynamic Fault Tree Analysis, DFTA) method, and proposed the use of Markov model for quantitative analysis, many scholars have improved its shortcomings, and has been widely cited. However, this kind of method only treats people as the bottom event without taking into account the cognitive characteristics, scene task characteristics and human-machine coupling characteristics. The Event Tree Analysis (ETA) method is a logical deductive, bottom-up analysis method, which can describe the human-computer interaction process from the causal logical relationship. The modeling process is simple, but it can only be used in a linear chain. However, the combined events (human-machine loop) cannot be analyzed in the way of description, and the description ability is insufficient. In terms of simulation modeling, the Risk and Reliability Center of the University of Maryland has developed a team information, decision, and action in CrewContext (IDAC) simulation analysis system; Taffy Institute, Massachusetts Institute of Technology, University of Chicago , University of Birmingham and other foreign research institutions have carried out research on multi-agent (MA) simulation technology and are developing corresponding platforms. Simulation-based modeling methods need to construct simulation models for various cognitive failure models and failure mechanism models in different task scenarios, and traversing event combinations or time sequences is very time-consuming, so the workload is huge and the versatility is poor. Simultaneous simulation The results are difficult to check.

目前,自动化复杂化的人机系统已经广泛地应用于众多重要领域。在一定程度上,人已经从操作者转换成为普通场景中的监视者和控制者、重要场景中的处理者和紧急场景的应急处理者。在各个场景中,人往往具有不同的信息获取渠道、以不同的感知方式从高度密集的信息群中不间断的获取信息、理解信息、筛选整理信息、分析加工信息并做出决策和相应操作,同时完成任务和确保安全。一定时间内复杂的认知处理活动会导致人的认知负荷的增加,降低人的认知行为能力和作业绩效水平,当认知负荷过高、飞行面临复杂气象条件时,信息的获取与分析都可能产生错误,也可能导致任务的失败或系统的失效并引发事故。因此从信息认知层面分析复杂系统人机交互可靠性至关重要。At present, the man-machine system with complex automation has been widely used in many important fields. To a certain extent, people have been transformed from operators to monitors and controllers in ordinary scenarios, handlers in important scenarios, and emergency handlers in emergency scenarios. In various scenarios, people often have different information acquisition channels, continuously acquire information from highly dense information groups in different ways of perception, understand information, filter and organize information, analyze and process information, and make decisions and corresponding operations. Complete missions and stay safe at the same time. Complex cognitive processing activities within a certain period of time will lead to an increase in human cognitive load and reduce human cognitive behavioral ability and operational performance. When the cognitive load is too high and the flight faces complex meteorological conditions, information acquisition and analysis Errors may occur, and may also lead to the failure of tasks or system failures and lead to accidents. Therefore, it is very important to analyze the reliability of human-computer interaction of complex systems from the level of information cognition.

认知过载和模式混淆是两种典型认知层面的人误行为。认知过载是指在时间有限的前提下,人由于认知资源有限而未能感知到所需的全部信息的状态。也就是说,认知过载故障与时间是否充足、认知资源限度和所需感知信息相关,导致的结果是有选择性的舍弃完成某一任务和感知与该任务相关的信息,因此未能感知到某些信息或忘记执行某项操作。模式混淆是指人由于感知信息缺失(由于认知过载引起的放弃感知故障诊断所需信息或由于设备故障/恶劣环境导致的故障诊断所需信息全部缺失)、感知信息不完全(部分缺失即交互信息不完整)或感知信息错误(即交互信息错误,包括虚警)而形成对系统故障模式的错误判断的状态。这两种人误模式能解释绝大部分的人机交互故障,因此在对人机交互系统进行建模分析时若能完整的刻画这两种故障逻辑,就能更加精确地分析系统的可靠性。Cognitive overload and pattern confusion are two typical cognitive-level human errors. Cognitive overload refers to the state in which people fail to perceive all the information they need due to limited cognitive resources under the premise of limited time. That is to say, the cognitive overload failure is related to the adequacy of time, the limit of cognitive resources and the required perception information, resulting in the selective abandonment of completing a certain task and perceiving information related to the task, thus failing to perceive to some information or forget to perform an action. Pattern confusion refers to the lack of perceptual information (the information required for perceptual fault diagnosis is abandoned due to cognitive overload, or the information required for fault diagnosis due to equipment failure/harsh environment is completely missing), and the perceptual information is incomplete (partial missing means interaction. Incomplete information) or perceptual information errors (that is, interactive information errors, including false alarms), resulting in a state of erroneous judgment on the system failure mode. These two human error modes can explain most of the human-computer interaction faults. Therefore, if these two fault logics can be completely described when modeling and analyzing the human-computer interaction system, the reliability of the system can be analyzed more accurately. .

发明内容SUMMARY OF THE INVENTION

(1)目的:(1. Purpose:

本发明提供一种基于马尔科夫的人机交互系统可靠度的求解方法,它从人机交互过程中人对信息的认知层面出发,充分考虑人的认知特性、场景任务特性和人机耦合特性,完善了系统故障逻辑,从而更加准确地度量了复杂人机系统的可靠性。通过分析能够发现容易诱发人为错误的场景和设计缺陷,对于提高人员绩效能力水平和增强人机系统可靠性和安全性有着重要作用。The invention provides a Markov-based method for solving reliability of human-computer interaction system, which starts from the cognitive level of human to information in the process of human-computer interaction, and fully considers human cognitive characteristics, scene task characteristics and human-computer interaction. The coupling characteristics improve the system fault logic, so as to measure the reliability of the complex human-machine system more accurately. Through analysis, it is possible to find scenarios and design defects that are easy to induce human error, which plays an important role in improving the level of human performance capabilities and enhancing the reliability and safety of human-machine systems.

(2)技术方案:(2) Technical solution:

本发明是一种基于马尔科夫的人机系统可靠度的求解方法,包括确定系统故障底事件集合、画出系统状态转移图、确定状态转移率矩阵和求解系统瞬时可靠度四步。该方法能够更加准确地求解出复杂人机系统的可靠度,并能发现容易诱发人为错误的场景和设计缺陷,对于提高人员绩效能力水平和增强人机系统可靠性和安全性有着重要作用。The invention is a Markov-based man-machine system reliability solution method, which includes four steps of determining system fault bottom event set, drawing system state transition diagram, determining state transition rate matrix and solving system instantaneous reliability. This method can more accurately solve the reliability of complex human-machine systems, and can find scenarios and design defects that are easy to induce human errors, which plays an important role in improving the level of human performance capabilities and enhancing the reliability and safety of human-machine systems.

本发明一种基于马尔科夫的人机系统可靠度的求解方法,其具体步骤介绍如下:A method for solving the reliability of a Markov-based man-machine system of the present invention, and its specific steps are introduced as follows:

步骤一:分析人机交互系统、环境和任务场景,分析故障逻辑确定导致人机系统故障的底事件集合;Step 1: Analyze the human-computer interaction system, environment and task scenarios, analyze the fault logic and determine the bottom event set that causes the human-computer system failure;

步骤二:确定状态空间的全集,根据系统故障逻辑画出系统状态转移图;Step 2: Determine the complete set of the state space, and draw a system state transition diagram according to the system fault logic;

步骤三:计算得出状态间转移率,构建转移率矩阵;Step 3: Calculate the transition rate between states and construct a transition rate matrix;

步骤四:列出状态方程,求解得出人机系统的瞬时可靠度;Step 4: List the state equation and solve to obtain the instantaneous reliability of the man-machine system;

通过以上步骤,分析了人机交互过程中的故障逻辑,确定了系统故障底事件的集合,并最终求解出人机系统的瞬时可靠度,达到了精确求解人机系统可靠度和查找系统缺陷的效果,解决了以往人机系统可靠性建模人机系统多因素耦合特性描述难的问题。Through the above steps, the fault logic in the human-computer interaction process is analyzed, the set of system fault bottom events is determined, and the instantaneous reliability of the human-computer system is finally solved, which achieves the accuracy of solving the reliability of the human-computer system and finding system defects. As a result, it solves the problem that the multi-factor coupling characteristics of the human-machine system is difficult to describe in the past reliability modeling of the human-machine system.

其中,在步骤一中所述的“确定导致人机系统故障的底事件集合”是整个分析和计算过程的基础;在分析系统故障逻辑时,应分为独立失效和人机交互故障两部分进行,并分别确定这两部分的底事件集合,通过求并集即可获得人机系统故障的底事件集合,包括以下步骤:Among them, the "determination of the bottom event set that causes the human-machine system failure" described in step 1 is the basis of the entire analysis and calculation process; when analyzing the system failure logic, it should be divided into two parts: independent failure and human-machine interaction failure. , and determine the bottom event sets of the two parts respectively, and obtain the bottom event set of the human-machine system failure by summing the union set, including the following steps:

步骤1):确定独立失效部分底事件集合Step 1): Determine the set of independent failure partial bottom events

独立失效部分底事件集合的确定采用的是故障树即FTA方法,该方法在工程实际中已经得到了广泛的应用,是一种较为成熟的方法;结合当前人机系统、环境和任务情景,对系统故障进行从上到下的逻辑分析,构建事故树即可获得独立失效部分的底事件集合;将此集合中人机交互故障传递至下一步骤,作为分析的基础;The fault tree or FTA method is used to determine the bottom event set of the independent failure part, which has been widely used in engineering practice and is a relatively mature method. The system fault is logically analyzed from the top to the bottom, and the bottom event set of the independent failure part can be obtained by constructing the fault tree; the human-computer interaction fault in this set is passed to the next step as the basis of the analysis;

步骤2):确定人机交互故障部分底事件集合Step 2): Determine the bottom event set of the human-computer interaction fault part

这部分的分析建立在故障树即FTA分析的基础之上,根据认知过载和模式混淆两种人机交互故障的定义和故障逻辑分析当前人机交互系统、环境和任务情景下发生的人机交互故障,从而确定人机交互故障部分的底事件集合;首先确定其触发事件即与当前情境下人所承担的识别系统机械故障的任务,以此作为认知过载故障的底事件;其次,分析导致人感知信息缺失、错误或不完全的因素,将其与需识别的设备故障作为模式混淆故障的底事件;二者并集即为人机交互故障部分的底事件集合;最后通过求独立失效部分和人机交互故障部分底事件集合的并集得到人机系统故障的底事件集合。This part of the analysis is based on the fault tree (FTA) analysis. According to the definition and fault logic of two types of human-computer interaction faults, cognitive overload and mode confusion, the human-computer interaction system, environment and task scenarios that occur in the current human-computer interaction system are analyzed. Interaction fault, so as to determine the bottom event set of the human-computer interaction fault part; first, determine the triggering event, which is the task of identifying the mechanical fault of the system under the current situation, as the bottom event of cognitive overload fault; secondly, analyze The factors that cause people to perceive the lack, error or incompleteness of information, and the equipment failure to be identified are regarded as the bottom event of the mode confusion failure; the combination of the two is the bottom event set of the human-computer interaction failure part; finally, by finding the independent failure part The union of the bottom event set of the human-computer interaction fault part obtains the bottom event set of the human-computer system fault.

其中,在步骤二中所述的“画出系统状态转移图”包括以下步骤:Wherein, "drawing a system state transition diagram" described in step 2 includes the following steps:

步骤1):定义系统状态Step 1): Define System State

所有底事件各种状态的随机组合构成了系统的状态,因此首先应定义底事件的状态;一般将底事件的发生状态定义为正常(用0表示)和故障(用1表示),则所有底事件的状态的零一组合构成了系统的不同状态,所有系统状态的集合称为状态空间;The random combination of all the states of the bottom event constitutes the state of the system, so the state of the bottom event should be defined first; generally, the occurrence state of the bottom event is defined as normal (represented by 0) and fault (represented by 1), then all the bottom events The zero-one combination of the state of the event constitutes the different states of the system, and the collection of all system states is called the state space;

步骤2):画出系统状态转移图Step 2): Draw a system state transition diagram

以底事件故障数作为横坐标,从上到下依次罗列系统的所有状态;接着根据系统的状态转移逻辑用有向线条连接状态点,该线条表示状态之间的转移;Taking the bottom event failure number as the abscissa, list all the states of the system from top to bottom; then connect the state points with a directed line according to the state transition logic of the system, and the line represents the transition between states;

认知过载的触发条件是输入事件数量a满足a≥2,而输入事件不同对应的输出事件的发生概率不同,因此在马尔科夫状态转移图中的同一条链中认知过载可能被触发多次;例如,在某一情境下人机系统中人承担A、B、C三项任务发生了认知过载故障,导致放弃了某些任务或者放弃感知与这些任务相关的信息D、E,这种情况就发生了触发了两次认知过载故障;The trigger condition of cognitive overload is that the number of input events a satisfies a ≥ 2, and the probability of occurrence of output events corresponding to different input events is different. Therefore, cognitive overload may be triggered more than once in the same chain in the Markov state transition diagram. For example, in a certain situation, in the human-machine system, a cognitive overload failure occurs when a person undertakes three tasks, A, B, and C, resulting in giving up some tasks or giving up perception of information D and E related to these tasks. In this case, two cognitive overload faults are triggered;

模式混淆的触发条件是条件判定事件发生,因此在状态转移图中触发条件的发生应和模式混淆合并;如当发生故障M时,由于事件P导致人将故障M错误识别为故障N;以此最终形成二维网状结构图;The triggering condition of the pattern confusion is the occurrence of a conditional judgment event, so the occurrence of the triggering condition in the state transition diagram should be combined with the pattern confusion; for example, when a fault M occurs, the fault M is mistakenly identified as a fault N due to the event P; Finally, a two-dimensional network structure diagram is formed;

步骤3):模型简化Step 3): Model Simplification

根据步骤一中故障逻辑的分析确定系统的故障状态点和非故障状态点;将系统故障状态集合F和非故障状态集合NF中的多个状态进行合并,简化模型同时,减少系统状态数;简化后的状态空间从0至n标注序号,并将其输入到下一步骤中从而达到简化求解过程的目的。Determine the fault state point and non-fault state point of the system according to the analysis of the fault logic in step 1; combine multiple states in the system fault state set F and non-fault state set NF to simplify the model and reduce the number of system states; simplify The latter state space is numbered from 0 to n, and it is input to the next step to simplify the solution process.

其中,在步骤三中所述的“计算得出状态间转移率,构建转移率矩阵”,其作法如下:计算状态间转移率应根据状态间转移形式进行分类,包括独立失效部分、认知过载故障和模式混淆故障三类;因此对于两不同状态间的转移率应进行分类讨论;状态间独立失效部分的转移率由专家根据实际经验和理论基础给出,人机交互部分的状态转移率由计算得出;最后将状态间转移率用矩阵的形式表示出来构建转移率矩阵;状态间转移率计算步骤如下:Among them, the method of "calculating the transition rate between states and constructing a transition rate matrix" described in step 3 is as follows: calculating the transition rate between states should be classified according to the form of transition between states, including independent failure parts, cognitive overload There are three types of faults and mode confusion faults; therefore, the transition rate between two different states should be classified and discussed; the transition rate of the independent failure part between states is given by experts based on actual experience and theoretical basis, and the state transition rate of the human-computer interaction part is given by Calculated; finally, the transition rate between states is expressed in the form of a matrix to construct a transition rate matrix; the calculation steps of the transition rate between states are as follows:

步骤1)状态i到状态j(i≠j)转移率求解Step 1) Solution of transition rate from state i to state j (i≠j)

若状态i转移到状态j与人机交互故障逻辑不相关且i≠j,则状态之间的转移仅包含独立失效部分,状态间转移率为:If the transition from state i to state j is not related to the human-computer interaction fault logic and i≠j, the transition between states only includes the independent failure part, and the transition rate between states is:

qij=λij (1)q ijij (1)

式中:λij为状态i至状态j(i≠j)独立失效部分状态转移率;In the formula: λ ij is the state transition rate of the independent failure part from state i to state j (i≠j);

若状态i转移到状态j与认知过载的输出事件相关且i≠j,则状态间转移包含独立失效部分和认知过载相关部分,状态间转移率为:If the transition from state i to state j is related to the output event of cognitive overload and i≠j, the transition between states includes the independent failure part and the cognitive overload-related part, and the transition rate between states is:

Figure BDA0001769200080000051
Figure BDA0001769200080000051

式中:Pc为认知过载输出事件在输入事件发生时发生的条件概率;where P c is the conditional probability that the cognitive overload output event occurs when the input event occurs;

若状态i转移到状态j与模式混淆的输出事件相关且i≠j,即该输出事件对应的输入事件发生状态一定导致输出事件发生,则将状态i和状态j合并,依据公式(1)(4.5)进行计算;If the transition from state i to state j is related to the output event of mode confusion and i≠j, that is, the occurrence state of the input event corresponding to the output event must lead to the occurrence of the output event, then the state i and state j are combined, according to formula (1)( 4.5) make calculations;

步骤2)状态i到状态j(i=j)的转移率的求解Step 2) Solution of the transition rate from state i to state j (i=j)

状态i转移到其余n个状态时,可能存在多个与认知过载相关的状态转移过程,则状态i到状态j(i=j)的转移率为:When state i transitions to the remaining n states, there may be multiple state transition processes related to cognitive overload, then the transition rate from state i to state j (i=j) is:

Figure BDA0001769200080000061
Figure BDA0001769200080000061

式中:m(m≤n)为状态i转移到其余n个状态时与认知过载相关的过程数;Pcb为第b(i=1,2,3Λm)个与认知过载相关的过程中认知过载输出事件在输入事件发生时发生的条件概率;where m (m≤n) is the number of processes related to cognitive overload when state i transitions to the other n states; P cb is the bth (i=1, 2, 3Λm) process related to cognitive overload The conditional probability that the output event of cognitive overload occurs when the input event occurs;

由此求得所有状态间转移率,将其以矩阵形式表示得转移率矩阵Q,形如:From this, the transition rates between all states are obtained, and they are expressed in matrix form to obtain the transition rate matrix Q, as follows:

Figure BDA0001769200080000062
Figure BDA0001769200080000062
;

其中,在步骤四中所述的“列出状态方程,求解得出人机系统的瞬时可靠度”,其作法如下:系统瞬时可靠度的求解实际上是求解线性微分方程,根据其求解过程包括如下步骤:Among them, the method of "listing the state equation and solving to obtain the instantaneous reliability of the man-machine system" described in step 4 is as follows: the solution of the instantaneous reliability of the system is actually solving the linear differential equation. According to the solution process, it includes: Follow the steps below:

步骤1):确定系统初始状态分布Step 1): Determine the initial state distribution of the system

将系统在t时刻时处于状态i的概率记为Pi(t),令P(t)=(P0(t),P1(t),Λ,Pn(t)),则根据状态方程有P(t)=P(0)·eQt;转移率矩阵Q已由步骤二中确定,因此应首先确定系统的初始状态分布,即t=0时刻P0(t),P1(t),Λ,Pn(t)的值,从而确定P(0),然后列出系统状态方程;Denote the probability that the system is in state i at time t as P i (t), let P(t)=(P 0 (t), P 1 (t), Λ, P n (t)), then according to the state The equation has P(t)=P(0)·e Qt ; the transition rate matrix Q has been determined in step 2, so the initial state distribution of the system should be determined first, that is, at t=0 time P 0 (t), P 1 ( t), Λ, P n (t) values to determine P(0), and then list the system state equation;

步骤2):求解系统状态方程Step 2): Solve the system state equation

首先将转移率矩阵Q对角化如式(4),为简化表达将X中各列记作n+1个n+1阶列向量A0,A1,ΛAn,将X-1中各行记作n+1个n+1阶行向量B0,B1,ΛBnFirst, the transition rate matrix Q is diagonalized as in formula (4). To simplify the expression, each column in X is denoted as n+1 n+1-order column vectors A 0 , A 1 , ΛA n , and each row in X -1 Denoted as n+1 n+1 order row vectors B 0 , B 1 , ΛB n ;

Figure BDA0001769200080000063
Figure BDA0001769200080000063

式中:Q为转移率矩阵;In the formula: Q is the transfer rate matrix;

X为变换矩阵;X is the transformation matrix;

X-1为X的逆矩阵;X -1 is the inverse matrix of X;

D为与Q相似的对角阵;D is a diagonal matrix similar to Q;

A0,A1,ΛAn为矩阵X中的n+1行;A 0 , A 1 , ΛA n are n+1 rows in matrix X;

B0,B1,ΛBn为矩阵X-1中的n+1列;B 0 , B 1 , ΛB n are the n+1 columns in the matrix X -1 ;

n为矩阵Q的阶数减去1;n is the order of the matrix Q minus 1;

d0,d1……dn为Q的n+1个特征根;d 0 , d 1 ......d n is the n+1 characteristic root of Q;

其次计算转移率矩阵的指数eQtNext, calculate the exponent e Qt of the transition rate matrix:

Figure BDA0001769200080000071
Figure BDA0001769200080000071

则各状态瞬时概率为:Then the instantaneous probability of each state is:

Figure BDA0001769200080000072
Figure BDA0001769200080000072

步骤3):确定系统瞬时可靠度Step 3): Determine the instantaneous reliability of the system

通过步骤2)的求解得出系统所有状态的瞬时概率值;将系统的失效状态集合F中的所有状态概率求和即可得到系统的瞬时失效概率;因此系统在t时刻故障的概率为:The instantaneous probability value of all states of the system is obtained through the solution of step 2); the instantaneous failure probability of the system can be obtained by summing all the state probabilities in the failure state set F of the system; therefore, the probability of the system failure at time t is:

Figure BDA0001769200080000073
Figure BDA0001769200080000073

其中,

Figure BDA0001769200080000074
in,
Figure BDA0001769200080000074

式中:P(t)为系统瞬时故障率;In the formula: P(t) is the instantaneous failure rate of the system;

t为时间,单位s;t is time, in s;

F为故障集合;F is the fault set;

则可知系统的瞬时可靠度为:Then the instantaneous reliability of the system can be known as:

R(t)=1-P(t) (8)R(t)=1-P(t) (8)

式中:R(t)为瞬时可靠度;In the formula: R(t) is the instantaneous reliability;

此步骤为求解方程的过程,其求解过程清晰简洁,可通过MATLAB、VB、C语言或C++编程替代人工求解。This step is the process of solving the equation, and the solving process is clear and concise, and manual solution can be replaced by MATLAB, VB, C language or C++ programming.

(3)功效、优点(3) Efficacy and advantages

本发明提出了一种基于马尔可夫的人机交互系统可靠性计算方法,考虑人的认知特性、场景任务特性和人机耦合特性,完整地描述了系统故障逻辑,能够适应复杂人机系统的可靠性定量建模分析的要求。同时提出了系统状态间转移率的计算方法,能够更加准确地计算出人机系统的可靠度。The invention proposes a Markov-based human-computer interaction system reliability calculation method, which considers human cognitive characteristics, scene task characteristics and human-computer coupling characteristics, completely describes the system fault logic, and can adapt to complex human-computer systems. Quantitative modeling analysis of reliability requirements. At the same time, a calculation method of the transition rate between system states is proposed, which can more accurately calculate the reliability of the man-machine system.

附图说明Description of drawings

图1本发明所述方法流程图。Fig. 1 is a flow chart of the method of the present invention.

图2双发舰载直升机坠毁事故树分析图。Figure 2. Analysis of the crash accident tree of a twin-engine carrier-based helicopter.

图3双发舰载直升机坠毁事故树简化图。Figure 3. A simplified diagram of the crash accident tree of a twin-engine carrier-based helicopter.

图4马尔科夫状态转移图。Figure 4. Markov state transition diagram.

图5马尔科夫状态转移简化图。Figure 5 Simplified diagram of Markov state transition.

图6事故分析结果对比图。Figure 6. Comparison of accident analysis results.

图中序号、符号、代号说明如下:The serial numbers, symbols and codes in the figure are explained as follows:

图2图3中英文字母加数字为事件代号,M代表中间事件、X代表底事件,FDEP为故障树中功能相关逻辑门;Fig. 2 and Fig. 3 Chinese and English letters plus numbers are the event codes, M represents the middle event, X represents the bottom event, and FDEP is the function-related logic gate in the fault tree;

图4图5中X加数字代表底事件代号,椭圆框代表底事件构成的系统的状态,数字0、1代表底事件的发生状态,0代表事件发生、1代表事件不发生,箭头代表系统的转移方向,实线代表设备独立失效引起的状态转移,虚线代表人误引起的状态转移。In Figure 4 and Figure 5, X plus a number represents the code of the bottom event, the ellipse represents the state of the system composed of the bottom event, the numbers 0 and 1 represent the occurrence state of the bottom event, 0 represents the event occurs, 1 represents the event does not occur, and the arrow represents the system In the direction of transition, the solid line represents the state transition caused by the independent failure of the equipment, and the dotted line represents the state transition caused by human error.

具体实施方式Detailed ways

本发明提供一种基于马尔科夫的人机系统可靠度求解方法,如图1所示;该方法依次按照下述四个步骤进行;本发明采用双发舰载直升机作为案例,在本案例中飞机处于着舰阶段,高度较低,当自动油门系统处于接通状态,直升机其中一个发动机发生油门卡阻故障时,飞机为保持进近速度而自动补偿引起的推力不对称情况,自动收回另一发动机的油门手柄,同时该发动机的参数急速下降。机组成员易将发动机的油门卡阻故障错误判断为另一个发动机故障,人工关停正常的发动机,导致飞机丧失动力而坠毁。The present invention provides a Markov-based man-machine system reliability solution method, as shown in FIG. 1 ; the method is carried out in sequence according to the following four steps; the present invention uses a twin-engine ship-based helicopter as a case, in this case The aircraft is in the landing stage and the altitude is low. When the auto-throttle system is turned on, and one of the helicopter's engines has a throttle jam failure, the aircraft automatically compensates for the asymmetric thrust caused by maintaining the approach speed, and automatically retracts the other. The throttle handle of the engine, while the parameters of the engine drop rapidly. The crew members easily misjudged the throttle jamming failure of the engine as another engine failure, and manually shut down the normal engine, causing the aircraft to lose power and crash.

本发明一种基于马尔科夫的人机系统可靠度的求解方法,见图1所示,具体实施方式详述如下:A method for solving the reliability of a Markov-based man-machine system of the present invention is shown in Figure 1, and the specific implementation is described in detail as follows:

步骤一:确定导致人机系统故障的底事件集合Step 1: Identify the bottom set of events that lead to the failure of the human-machine system

分为以下两个步骤:Divided into the following two steps:

步骤1):确定独立失效部分底事件集合Step 1): Determine the set of independent failure partial bottom events

利用FTA方法对系统进行建模分析,即可得到独立失效部分底事件集合。在本例中以双发舰载直升机坠毁为顶事件,通过自上而下的逻辑分析构建事故树,如图2,其中不包含认知过载和模式混淆两种故障逻辑。为简化计算过程将中间事件(M4,M5,M9,M11)省略,简化的模型如图3。则可得独立失效部分底事件见下表:Using the FTA method to model and analyze the system, the set of independent failure partial bottom events can be obtained. In this example, the crash of a dual-engine carrier-based helicopter is taken as the top event, and an accident tree is constructed through top-down logic analysis, as shown in Figure 2, which does not include the two fault logics of cognitive overload and mode confusion. In order to simplify the calculation process, the intermediate events (M4, M5, M9, M11) are omitted, and the simplified model is shown in Figure 3. Then the bottom events of the independent failure part can be obtained as shown in the following table:

表1独立失效部分底事件表Table 1 Bottom event table of independent failure parts

编号Numbering 事件event X1X1 油门卡阻Throttle stuck X2X2 人工错误关停2号发动机Manual error shutting down the No. 2 engine X3X3 未查看发动机参数Engine parameters not checked X4X4 2号发动机油门手柄收回No. 2 engine throttle handle retracted

步骤2):确定人机交互故障部分底事件集合Step 2): Determine the bottom event set of the human-computer interaction fault part

本案例中既包括认知过载故障和由认知过载导致的模式混淆故障。首先应分析认知分析过载故障逻辑,在本案例中人承担正常飞行任务、通讯型任务以及识别识别机械故障任务,此故障表征为2号发动机手柄收回。这三项任务均需占用人的认知资源导致认知过载故障的发生从而导致人放弃识别发动机参数。结合故障表征,从而发生模式混淆导致人将1号发动机故障错认为是2号发动机故障而错误关停2号发动机。根据认知过载和模式混淆底事件的定义将二者底事件列于表2:This case includes both cognitive overload failure and pattern confusion failure caused by cognitive overload. First of all, the logic of cognitive analysis overload fault should be analyzed. In this case, people undertake normal flight tasks, communication tasks, and tasks of identifying and identifying mechanical faults. This fault is characterized by the retraction of the No. 2 engine handle. These three tasks all need to occupy human cognitive resources, leading to cognitive overload failure, which leads people to give up identifying engine parameters. Combined with the fault characterization, the resulting pattern confusion led people to mistake the No. 1 engine failure as the No. 2 engine fault and erroneously shut down the No. 2 engine. According to the definition of cognitive overload and mode confusion bottom events, the two bottom events are listed in Table 2:

表2人机交互故障部分底事件集合Table 2 The bottom event set of the human-computer interaction fault part

编号Numbering 事件event X4X4 2号发动机油门手柄收回No. 2 engine throttle handle retracted

从而可知系统故障底事件集合为{X1,X4,X2,X3}。Thus it can be known that the set of system fault bottom events is {X1, X4, X2, X3}.

步骤二:画出系统状态转移图包括以下步骤Step 2: Draw the system state transition diagram including the following steps

步骤1):定义系统状态Step 1): Define System State

对系统状态的定义应首先定义底事件的状态,将底事件发生记为状态1,底事件未发生记为状态0,即可明确区分底事件各状态。本案例中共有四个底事件,这四个事件状态的零一组合构成系统的状态空间的全集,该集合中共有16个状态点。The definition of the system state should first define the state of the bottom event, record the occurrence of the bottom event as state 1, and record the occurrence of the bottom event as state 0, so that the states of the bottom event can be clearly distinguished. There are four bottom events in this case, and the zero-one combination of these four event states constitutes the complete set of the state space of the system, and there are 16 state points in this set.

步骤2):画出系统状态转移图Step 2): Draw a system state transition diagram

将各底事件均未发生设为初始状态点,以底事件发生数为横坐标,从上到下以此罗列系统状态,根据步骤一中故障逻辑的分析依次用有向线条连接各系统状态点,形成网状逻辑图,如图4。Take the number of bottom events as the abscissa, and list the system states from top to bottom. According to the analysis of the fault logic in step 1, connect the system state points with directional lines in turn. , forming a network logic diagram, as shown in Figure 4.

步骤3):模型简化Step 3): Model Simplification

根据步骤一中故障逻辑的分析确定系统的故障状态点和非故障状态点。将系统故障状态集合F和非故障状态集合NF中的多个状态进行合并,剩余12个状态空间,如图5。简化后的状态空间从0至11标注序号,并将其输入到下一步骤中从而达到简化求解过程的目的。According to the analysis of the fault logic in step 1, the fault state point and the non-fault state point of the system are determined. The multiple states in the system fault state set F and the non-fault state set NF are merged, and 12 state spaces remain, as shown in Figure 5. The simplified state space is numbered from 0 to 11 and input to the next step to simplify the solution process.

步骤三:计算状态间转移率Step 3: Calculate the transition rate between states

对于系统间状态转移率的计算应首先确定独立失效部分的状态转移率。而在系统状态转移时其独立失效部分是由于底事件的发生引起的,因此独立失效部分的状态转移率就等于引起系统状态转移的底事件的发生率。做出如下假设:For the calculation of the state transition rate between systems, the state transition rate of the independent failure part should be determined first. The independent failure part of the system state transition is caused by the occurrence of the bottom event, so the state transition rate of the independent failure part is equal to the occurrence rate of the bottom event that causes the system state transition. Make the following assumptions:

λX1=λX2=λX3=λX4=1×10-4 (9)λ X1X2X3X4 =1×10 −4 (9)

式中λX1、λX2、λX3、λX4为底事件X1、X2、X3、X4的发生率。where λ X1 , λ X2 , λ X3 , and λ X4 are the occurrence rates of the bottom events X1, X2, X3, and X4.

步骤1)状态i到状态j(i≠j)转移率求解Step 1) Solution of transition rate from state i to state j (i≠j)

首先应分析系统状态间转移的故障分类,在本案例中包括独立失效部分和人机交互故障部分,在人机交互故障部分模式混淆由认知过载引起可将两种故障合为一个整体分析。经分析将系统间状态转移关系列于表3:First of all, the fault classification of system state transition should be analyzed. In this case, it includes the independent failure part and the human-computer interaction fault part. In the human-computer interaction fault part, the mode confusion is caused by cognitive overload, and the two faults can be combined into one overall analysis. After analysis, the state transition relationship between systems is listed in Table 3:

表3系统间状态转移分类表Table 3 Classification of state transitions between systems

Figure BDA0001769200080000111
Figure BDA0001769200080000111

则有仅包含独立失效部分的状态间转移率为:Then there is a transition rate between states containing only independent failure parts:

q01=q21=q39=q45=q65=q79=q8,11=q10,11=λX1 q 01 =q 21 =q 39 =q 45 =q 65 =q 79 =q 8,11 =q 10,11X1

q02=q37=q46=q8,11=λX4 q 02 =q 37 =q 46 =q 8,11X4

q03=q19=q27=q48=q5,11=q6,10=λX2 q 03 =q 19 =q 27 =q 48 =q 5,11 =q 6,10X2

q04=q15=q26=q38=q7,10=q9,11=λX2 q 04 =q 15 =q 26 =q 38 =q 7,10 =q 9,11X2

对于包含人机交互故障部分,应首先确定认知过载输出事件在输入事件发生是的条件概率。在本案例中,认知过载输入事件为人所承担的正常飞行任务、通讯型任务以及识别识别机械故障任务,输出事件为放弃识别发动机参数。假设此条件概率Pc=0.8,则人机交互故障部分状态间转移率为:For the part containing human-computer interaction faults, the conditional probability that the cognitive overload output event occurs at the input event should be determined first. In this case, the input events of cognitive overload are normal flight tasks, communication tasks and tasks of identifying and identifying mechanical faults undertaken by humans, and the output events are abandoning the identification of engine parameters. Assuming this conditional probability P c =0.8, the transition rate between states of the human-computer interaction fault part is:

q1,11=(λX2X3)Pc q 1,11 =(λ X2X3 )P c

q2,10=(λX1X2X3)Pc q 2,10 =(λ X1X2X3 )P c

步骤2)状态i到状态j(i=j)的转移率的求解Step 2) Solution of the transition rate from state i to state j (i=j)

已知在系统全部12个状态中仅有状态1与状态2在转移至其他状态过程中存在认知过载故障。在两过程均中认知过载的输入事件和输出事件均相同,其条件概率为Pc=0.8。则有状态间转移率:It is known that among all 12 states of the system, only state 1 and state 2 have cognitive overload failures in the process of transitioning to other states. The input and output events of cognitive overload were the same in both processes, with a conditional probability of P c =0.8. Then there is the transition rate between states:

q00=-(λX1X4X2X3)q 00 =-(λ X1X4X2X3 )

q11=-(λX2X3)(1-pc)q 11 =-(λ X2X3 )(1-p c )

q22=-(λX1X2X3)(1-pc)q 22 =-(λ X1X2X3 )(1-p c )

q33=-(λX1X4X3)q 33 =-(λ X1X4X3 )

q44=-(λX1X4X2)q 44 =-(λ X1X4X2 )

q55=-λX2 q 55 =-λ X2

q66=-λX1X2 q 66 = -λ X1X2

q77=-λX1X3 q 77 = -λ X1X3

q88=-λX1X4 q 88 = -λ X1X4

q99=-λX3 q 99 = -λ X3

q10,10=-λX1 q 10,10 = -λ X1

其余系统间状态转移率均为0,则系统间状态转移率矩阵为:The other inter-system state transition rates are all 0, then the inter-system state transition rate matrix is:

Figure BDA0001769200080000121
Figure BDA0001769200080000121

步骤四:求解系统瞬时可靠度Step 4: Solve the instantaneous reliability of the system

首先,应确定系统初始时刻状态分布,在本案例中将底事件均为发生的状态定义为系统的初始状态即状态点(0000),则可知系统的初始状态分布为P(0)=(1,0,Λ,0)12,之后可利用编程的手段对系统可靠度进行求解。为对比方法的效果,本案例求解出仅考虑独立失效部分的系统瞬时可靠度如下:First of all, the state distribution at the initial time of the system should be determined. In this case, the state where the bottom events are all occurrences is defined as the initial state of the system, that is, the state point (0000), then it can be known that the initial state distribution of the system is P(0)=(1 ,0,Λ,0) 12 , and then the system reliability can be solved by means of programming. In order to compare the effect of the method, the instantaneous reliability of the system considering only the independent failure part is solved in this case as follows:

Figure BDA0001769200080000122
Figure BDA0001769200080000122

画出同时考虑独立失效和人机交互故障的系统瞬时可靠度和仅考虑独立失效时的系统瞬时可靠度对比图,如图6,可以看出前者显著大于后者。Draw a comparison chart of the instantaneous reliability of the system considering both independent failures and human-computer interaction failures and the instantaneous reliability of the system when only considering independent failures, as shown in Figure 6. It can be seen that the former is significantly greater than the latter.

Claims (2)

1. A human-computer interaction system reliability solving method based on Markov is characterized in that: the method comprises the following steps:
the method comprises the following steps: analyzing a human-computer interaction system, an environment and a task scene, and analyzing fault logic to determine a bottom event set causing a human-computer system fault;
step two: determining a complete set of state space, and drawing a system state transition diagram according to system fault logic;
step three: calculating to obtain the transfer rate between states, and constructing a transfer rate matrix;
step four: listing a state equation, and solving to obtain the instantaneous reliability of the human-computer system;
the step one of determining a bottom event set causing the man-machine system fault is the basis of the whole analysis and calculation process; when analyzing system fault logic, the method is divided into two parts of independent failure and man-machine interaction fault, bottom event sets of the two parts are respectively determined, and the bottom event set of the man-machine system fault can be obtained by solving a union set, and the method comprises the following steps:
step 1.1): determining a set of independent failure partial bottom events
Determining the independent failure part bottom event set by adopting a Fault Tree (FTA) method, and performing top-to-bottom logic analysis on system faults by combining the current man-machine system, environment and task situation to construct a fault tree, namely obtaining the independent failure part bottom event set; transmitting the man-machine interaction faults in the set to the next step to serve as a basis for analysis;
step 1.2): determining a set of partial bottom events of human-computer interaction faults
The analysis of the part is established on the basis of Fault Tree (FTA) analysis, and human-computer interaction faults occurring under the current human-computer interaction system, environment and task situations are analyzed according to the definition and fault logic of cognitive overload and mode confusion of two human-computer interaction faults, so that a bottom event set of the human-computer interaction fault part is determined; firstly, determining a trigger event, namely a task of identifying a system mechanical fault borne by a person under the current situation, and taking the trigger event as a bottom event of the cognitive overload fault; secondly, analyzing factors causing human perception information missing, errors and incompleteness, and taking the factors and equipment faults needing to be identified as bottom events of mode confusion faults; the union set of the two events is a bottom event set of the human-computer interaction fault part; finally, obtaining a bottom event set of the man-machine system fault by solving a union set of bottom event sets of the independent failure part and the man-machine interaction fault part;
the step two of drawing the system state transition diagram comprises the following steps:
step 2.1): defining system states
The random combination of all the states of the bottom event constitutes the state of the system, so the state of the bottom event should be defined first; the occurrence state of the bottom event is represented as normal by '0', and the occurrence state of the bottom event is represented as fault by '1', so that zero combination of the states of all the bottom events forms different states of the system, and a set of all the states of the system is called as a state space;
step 2.2): drawing a system state transition diagram
Taking the number of faults of the bottom event as an abscissa, and sequentially listing all states of the system from top to bottom; then, according to the state transition logic of the system, connecting state points by using a directed line, wherein the line represents the transition between states;
the cognitive overload triggering condition is that the number a of input events is more than or equal to 2, and the occurrence probability of output events corresponding to different input events is different, so that the cognitive overload is triggered for multiple times in the same chain in the Markov state transition diagram;
the trigger condition of the mode confusion is that a condition judgment event occurs, so the occurrence of the trigger condition in the state transition diagram should be merged with the mode confusion; if a fault M occurs, a person mistakenly identifies the fault M as a fault N due to an event P; thereby finally forming a two-dimensional mesh structure diagram;
step 2.3): model simplification
Determining fault state points and non-fault state points of the system according to the analysis of the fault logic in the step one; combining a plurality of states in a system fault state set F and a non-fault state set NF, simplifying a model and reducing the number of system states; the simplified state space is marked with serial numbers from 0 to n and is input into the next step, so that the aim of simplifying the solving process is fulfilled;
the "calculating inter-state transition rates, constructing the transition rate matrix" described in step three is performed as follows: calculating the transfer rate among the states, and classifying the transfer rate among the states according to the transfer form among the states, wherein the transfer rate includes three types, namely an independent failure part, a cognitive overload fault and a mode confusion fault; therefore, the transition rate between two different states should be discussed in a classification way; the transfer rate of the independent failure part between the states is given by experts according to actual experience and theoretical basis, and the state transfer rate of the human-computer interaction part is obtained by calculation; finally, expressing the transition rate between states in a matrix form to construct a transition rate matrix; the calculation steps of the transition rate between states are as follows:
step 3.1) solving the transition rate from the state i to the state j;
if the state i is transferred to the state j which is not related to the man-machine interaction fault logic and i is not equal to j, the transfer between the states only comprises an independent failure part, and the transfer rate between the states is as follows:
q ij =λ ij (1)
in the formula: lambda ij The state transition rate of the independent failure part from the state i to the state j;
if the state i is transferred to the state j and is related to the output event of the cognitive overload, and i is not equal to j, the inter-state transfer comprises an independent failure part and a cognitive overload related part, and the inter-state transfer rate is as follows:
Figure FDA0003762561320000031
in the formula: p is c Outputting a conditional probability of an event occurring at the time of an input event for cognitive overload;
if the state i is transferred to the state j and is related to the mode-confused output event and i is not equal to j, namely the output event is generated because the input event generating state corresponding to the output event definitely causes the output event, combining the state i and the state j, and calculating according to a formula (1);
step 3.2) solving the transfer rate from the state i to the state j;
when the state i is transferred to the rest n states, a plurality of state transfer processes related to cognitive overload exist, and the transfer rate from the state i to the state j is as follows:
Figure FDA0003762561320000032
in the formula: m is the number of processes related to cognitive overload when the state i is transferred to the other n states, wherein m is less than or equal to n; i =1,2,3 … m; p cb The conditional probability of the cognitive overload output event occurring when the input event occurs in the process related to the cognitive overload is the b < th > process;
therefore, the transition rate among all states is obtained, and the transition rate matrix Q is expressed in a matrix form and is in the form of:
Figure FDA0003762561320000033
2. the markov-based human-computer interaction system reliability solution method of claim 1, wherein: "list the state equations and solve to get the instantaneous reliability of the human-machine system" described in step four, which is done as follows: the solution of the system instantaneous reliability is actually to solve a linear differential equation, and the solution process comprises the following steps:
step 4.1): determining system initial state distribution
Let P denote the probability that the system is in state i at time t i (t), let P (t) = (P) 0 (t),P 1 (t),…,P n (t)), then there is P (t) = P (0) · e) according to the equation of state Qt (ii) a The transfer rate matrix Q has been determined in step two, so the initial state distribution of the system should be determined first, i.e. at time P when t =0 0 (t),P 1 (t),…,P n (t) to determine P (0), and then listing the system equation of state;
step 4.2): solving system state equation
First, the transfer rate matrix Q is diagonalized as shown in formula (4), and each column in X is recorded as n +1 order column vectors A for simplifying expression 0 ,A 1 ,…A n Is mixing X -1 The middle lines are marked as n +1 order line vectors B 0 ,B 1 ,…B n
Figure FDA0003762561320000041
In the formula: q is a transfer rate matrix;
x is a transformation matrix;
X -1 an inverse matrix of X;
d is a diagonal matrix similar to Q;
A 0 ,A 1 ,…A n is n +1 rows in the matrix X;
B 0 ,B 1 ,…B n is a matrix X -1 N +1 columns of (1);
n is the order of the matrix Q minus 1;
d 0 ,d 1 ……d n n +1 characteristic roots for Q;
second, the index e of the transfer rate matrix is calculated Qt
Figure FDA0003762561320000042
Then the transient probability of each state is:
Figure FDA0003762561320000043
step 4.3): determining system instantaneous reliability
Obtaining instantaneous probability values of all states of the system through the solution of the step 4.2); summing all the state probabilities in the failure state set F of the system to obtain the instantaneous failure probability of the system; the probability of the system failing at time t is therefore:
Figure FDA0003762561320000044
wherein,
Figure FDA0003762561320000051
in the formula: p (t) is the instantaneous failure rate of the system;
t is time, unit s;
f is a fault set;
then the instantaneous reliability of the system is known as:
R(t)=1-P(t) (8)
in the formula: r (t) is the instantaneous reliability.
CN201810941561.XA 2018-08-17 2018-08-17 Human-computer interaction system reliability solving method based on Markov Active CN109241583B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810941561.XA CN109241583B (en) 2018-08-17 2018-08-17 Human-computer interaction system reliability solving method based on Markov

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810941561.XA CN109241583B (en) 2018-08-17 2018-08-17 Human-computer interaction system reliability solving method based on Markov

Publications (2)

Publication Number Publication Date
CN109241583A CN109241583A (en) 2019-01-18
CN109241583B true CN109241583B (en) 2022-10-18

Family

ID=65071467

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810941561.XA Active CN109241583B (en) 2018-08-17 2018-08-17 Human-computer interaction system reliability solving method based on Markov

Country Status (1)

Country Link
CN (1) CN109241583B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110322977B (en) * 2019-07-10 2021-02-09 河北工业大学 A reliability analysis method of nuclear power reactor core water level monitoring system
CN112183755A (en) * 2020-10-21 2021-01-05 中国核动力研究设计院 Markov model construction simplification method and system applied to complex system
CN112433609B (en) * 2020-11-19 2022-03-15 北京航空航天大学 A Multi-agent-based Information Hierarchical Human-Computer Interaction Security Modeling Method
CN113836741B (en) * 2021-09-30 2024-01-26 中国工程物理研究院研究生院 Reconstruction and reliability evaluation method based on multi-functional system degradation process
CN114329911B (en) * 2021-12-07 2024-04-16 中国航空综合技术研究所 Flight control system function fault analysis method based on scene model
CN115114804A (en) * 2022-07-28 2022-09-27 哈尔滨理工大学 Markov state transfer-based water surface unmanned patrol boat multi-state system modeling method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663408A (en) * 2011-12-31 2012-09-12 电子科技大学 Backup structure-oriented fault tree analysis method
CN104392072A (en) * 2014-12-12 2015-03-04 北京交通大学 Equipment integration system reliability analysis method based on dynamic fault tree
CN104898636A (en) * 2015-03-15 2015-09-09 国家电网公司 Safety and stability control device reliability analysis method in consideration of multistate operation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663408A (en) * 2011-12-31 2012-09-12 电子科技大学 Backup structure-oriented fault tree analysis method
CN104392072A (en) * 2014-12-12 2015-03-04 北京交通大学 Equipment integration system reliability analysis method based on dynamic fault tree
CN104898636A (en) * 2015-03-15 2015-09-09 国家电网公司 Safety and stability control device reliability analysis method in consideration of multistate operation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An enhanced CREAM with stakeholder-graded protocols for tanker shipping safety application;Qingji Zhou 等;《Safety Science》;20170227;第95卷;140-147 *
动态系统故障的混杂传播特征及建模方法;郭健彬等;《系统工程与电子技术》;20150131;第37卷(第01期);224-228 *

Also Published As

Publication number Publication date
CN109241583A (en) 2019-01-18

Similar Documents

Publication Publication Date Title
CN109241583B (en) Human-computer interaction system reliability solving method based on Markov
CN112100865A (en) Method for predicting remaining life of aircraft engine based on parallel CNN model
CN106250631B (en) Fault diagnosis method based on fault-test correlation matrix
CN104134010B (en) Satellite fault diagnosis method for discrete type data based on Naive Bayes
CN108334907A (en) A kind of complex equipment point method for detecting abnormality and system based on deep learning
CN111260261B (en) Automatic evaluation method for probability risk of system
US8560279B2 (en) Method of determining the influence of a variable in a phenomenon
JP7527702B1 (en) Bridge deformation monitoring method, computer device, and computer-readable storage medium
CN103810383A (en) Method for identifying and analyzing single-point failure mode of launch vehicle
CN105512483A (en) Remaining Life Prediction of Aeroengine Based on Damage Propagation Model and Data Analysis
CN114925536A (en) Airborne system PHM testability modeling and diagnosis strategy optimization method and device
CN108519769A (en) A Rule-Based Method for Fault Diagnosis of Flight Control System
De Felice et al. A hybrid model for human error probability analysis
Jaw et al. Engine diagnostics in the eyes of machine learning
Guo et al. Deep-learning-based model for accident-type prediction during approach and landing
CN106777745B (en) Markov-based security assessment method
Zibaei et al. Diagnosis of safety incidents for cyber-physical systems: A uav example
Zhao et al. Remote intelligent expert system for operation state of marine gas turbine engine
CN117708741A (en) A method for anomaly detection of UAV multi-dimensional flight data for parameter selection
Rabenoro et al. A methodology for the diagnostic of aircraft engine based on indicators aggregation
Vesely Probabilistic risk assessment
Mack et al. Using tree augmented Naıve Bayes classifiers to improve engine fault models
Ding et al. A method of error mode effect analysis for a human-computer interaction system in aviation
CN106528401B (en) The tripping test method and device of control system based on redundancy structure
Yang et al. An expert system reasoning machine based on the combination of fault tree and generalized regression neural network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant