阅读说明：本技术 一种基于多源数据信息的疲劳状态因果网络方法 (fatigue state causal network method based on multi-source data information ) 是由 任长娥 袁超 杜涛 王岩 于 2019-09-26 设计创作，主要内容包括：本发明公开一种基于多源数据信息的疲劳状态因果网络方法。步骤如下：一、采集与人体疲劳状态相关的呼吸、心率数据；二、针对每一项生理数据,使用外推拟合方法去除野值；三、任选两项生理数据,检验该两项数据之间是否具有相关关系；四、对于具有相关关系的两项生理数据,采用格兰杰因果分析方法检验二者之间是否具有因果关系；五、循环选择两组生理数据,运用第四步中的格兰杰因果分析方法；六、遍历所有生理数据后,建立变量间的因果关系网络。本发明的方法具有精度高、计算量小、兼容性强等优点。可用于汽车、轮船、飞机等载体工作人员的疲劳状态因果分析与检测,提高基于多源数据信息融合的疲劳状态因果分析准确率,具有较高应用价值。(The invention discloses fatigue state causal network methods based on multi-source data information, which comprises the following steps of collecting respiration and heart rate data related to human fatigue states, secondly removing outliers by using an extrapolation fitting method aiming at each physiological data, thirdly selecting two optional physiological data to check whether the two data have a correlation relationship, fourthly checking whether the two physiological data have the causal relationship by using a Glange causal analysis method for the two physiological data with the correlation relationship, fifthly circularly selecting two groups of physiological data, applying the Gegere causal analysis method in the fourth step, and sixthly establishing a causal relationship network among variables after traversing all the physiological data.)

一种基于多源数据信息的疲劳状态因果网络方法

技术领域

本发明涉及一种基于多源数据信息的疲劳状态因果网络方法，可用于汽车、轮船、飞机等载体工作人员的疲劳状态因果分析与检测，提高基于多源数据信息融合的疲劳状态因果分析准确率。

背景技术

在特殊环境下作业，人的某些生理指标值会发生变化，这些生理数据的变化可能会导致某些生理疾病的发生。例如长时间在密闭的环境中作业，可能更容易疲劳、抑郁等。疲劳状况是对某些生理亚健康状态的预警，因此研究特殊作业环境下人的疲劳状态将有助于对人的健康状态进行预测，以便合理改善他们的工作条件，这具有重要意义。

论文《基于深度学习和面部多源动态行为融合的驾驶员疲劳检测方法研究》提出的机器学习方法可以通过学习某些生理指标的数据集，来判断人体是否处于疲劳状态。但是操作过程复杂，并且无法给出这些生理指标之间的因果关系。

格兰杰因果关系分析不仅能够描述不同生理特征数据相互作用的强度，还可以表明它们之间相互作用的因果关系，这与传统的分析方法如相关、相干及互信息等算法相比具有明显优势。因此，结合格兰杰因果关系分析人的生理特征参数之间的关系有着显著的优势。由于因果关系是有方向、客观存在的，所以可以通过“操纵”原因来获得预期的结果，使得人们对自然世界获得较强的主观能动性，这使得人们研究因果关系具有重要意义。

发明内容

本发明的目的在于提供一种基于多源数据信息的疲劳状态因果网络方法，以克服上述机器学习方法的缺陷，将格兰杰因果分析方法应用于影响人体疲劳状态的各项生理指标进行分析。本发明的方法经过数据采集与预处理、相关关系检验、格兰杰因果关系检验、构建因果网络等步骤，可以快速找到各项生理指标之间的因果关系。本发明的方法可用于汽车、轮船、飞机等载体工作人员的疲劳状态因果分析与检测，提高基于多源数据信息融合的疲劳状态因果分析准确率，具有较高的应用价值。

本发明的技术解决方案为：

一种基于多源数据信息的疲劳状态因果网络方法，其特征在于包括以下步骤：

第一步，采集与人体疲劳状态相关的呼吸、心率等生理数据；

具体包括：利用心肺功能测定仪收集呼吸数据，利用心电监测仪收集心率数据，上述呼吸数据包括：呼吸商(R)、呼吸频率(Rf)、潮波式流气量(VT)、通气量(VE)、每分钟氧气消耗量(VO₂)、每分钟二氧化碳消耗量(VCO₂)、呼气末氧分压(PetO₂)、呼气末二氧化碳分压(PetCO₂)每分钟能量消耗量(EEm)；上述心率数据包括：心率(HR)。

第二步，针对每一项项生理数据，去除野值；

具体可使用外推拟合方法去除野值，步骤包括：针对所采集的呼吸和心率数据中的每一单项生理数据，任选从第n时刻至第n+4时刻的五个连续的测量数据，记为z_n,z_n+1,z_n+2,z_n+3,z_n+4，依据前四个时刻数据，计算第n+4时刻的估计值如下：

计算出估计值

后与第n+4时刻测量值z_n+4进行比较，当满足时，则认为z_n+4是正常值，否则认为z_n+4为野值，将其去除。其中σ为均方误差，可以通过查阅传感器说明书获得。

第三步，任选两项生理数据，检验该两项数据之间是否具有相关关系；

具体包括：从第一步所采集的生理数据中任选两项，检验它们之间是否具有相关关系。相关关系计算公式如下：

其中X与Y分别表示两组变量，cov(X,Y)表示X与Y的协方差，

为X的标准差，

为Y的标准差，ρ_XY表示X与Y的相关系数。设定一个阈值e，当|ρ_XY|≥e时，认为变量X与Y之间存在相关关系。对任意两项生理数据，均执行上述相关关系计算，找出所有存在相关关系的生理指标。

第四步，采用格兰杰因果分析方法检验第三步得到的存在相关关系的其中两组数据之间是否具有因果关系；

具体包括：如果两项生理指标之间的相关系数ρ_XY≥0.8，则认为这两项生理指标之间存在相关关系。当两个变量之间存在相关关系时，才进一步检验它们之间的因果关系。因果关系的检验按照如下方式进行：

假定X(t)和Y(t)为两个具有相关关系的随机平稳时间序列，则变量Y(t)可以用它的自回归方程和它与X(t)之间的相互作用来表达，表达式为：

其中，p,q分别为X(t)和Y(t)回归方程中的滞后阶数；i和j为滞后项数，i的取值为1到p之间的整数，j的取值为1到q之间的整数；X(t-i)为X(t)的第i个滞后项，Y(t-j)为Y(t)的第j个滞后项，α_YX(i)和α_YY(i)为Y(t)的系数估计值，e(t)为Y(t)的白噪声。如果α_YX(i)(i＝1,2,...p)在统计学上整体显著不为零，则说明X(t)是Y(t)的格兰杰原因。如果α_YX(i)(i＝1,2,...p)在统计学上整体显著等于零，则说明X(t)和Y(t)不存在因果关系。

第五步，循环选择两组生理数据，运用第四步中的格兰杰因果分析方法；

具体包括：对于第三步中已确定的具有相关关系的生理数据，运用第四步中的格兰杰因果检验方法，检验任意两组数据之间的格兰杰因果关系，直到所有数据之间的因果关系均被检验。

第六步，遍历所有生理数据后，建立变量间的因果关系网络。

具体包括：对于第五步已经检验出的格兰杰因果关系，画出各项数据之间的因果关系网络图。用箭头表示数据之间的因果关系方向，箭头由“原因”变量指向“结果”变量。

本发明与现有技术相比的优点在于：

(1)本发明的方法操作简单，且具有精度高、计算量小的优点。

(2)本发明的方法可以帮助人们快速找到各项生理数据之间的因果关系网络，从而通过控制“原因”变量，得到预期的结果。本发明的方法可用于汽车、轮船、飞机等载体工作人员的疲劳状态因果分析与检测，提高基于多源数据信息融合的疲劳状态因果分析准确率，具有较高的应用价值。

附图说明

图1为本发明的方法的实现流程图；

图2为本发明的一个实施例的因果网络图；

具体实施方式

以下结合实施例及附图对本发明进行清楚、完整地的描述。显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

在本发明中若涉及“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。

请参照图1，图1为本发明的方法的实现流程图。

步骤S1，用K4b²心肺功能测定仪收集9个被测试者的呼吸数据，用心电监测仪收集被测试者的心率数据。上述数据呼吸包括：呼吸商(R)、呼吸频率(Rf)、潮波式流气量(VT)、通气量(VE)、每分钟氧气消耗量(VO₂)、每分钟二氧化碳消耗量(VCO₂)、呼气末氧分压(PetO₂)、呼气末二氧化碳分压(PetCO₂)、每分钟能量消耗量(EEm)；上述心率数据包括：心率(HR)。以上共计10项生理指标，每项生理指标采集369个数据，每次采集的上述10项数据一一对应保存。

步骤S2，针对所采集的呼吸和心率数据中的每一项生理数据，任选从第n时刻至第n+4时刻的五个连续的测量数据，记为z_n,z_n+1,z_n+2,z_n+3,z_n+4，依据前四个时刻数据，计算第n+4时刻的估计值

将估计值

与第n+4时刻测量值z_n+4进行比较，如果

则认为z_n+4是正常值，否则认为z_n+4为野值，将其去除。其中σ为均方误差，可以通过查阅传感器说明书获得。

步骤S3，检验10项生理数据之间是否具有相关关系。相关系数计算公式如下：

其中X与Y分别表示两组变量，cov(X,Y)表示X与Y的协方差，

为X的标准差，

为Y的标准差，ρ_XY表示X与Y的相关系数。设定一个阈值e，当|r_XY|≥e时，认为变量X与Y之间存在相关关系。在本发明的实施例中，设定e＝0.8，当ρ_XY≥0.8时，认为两项生理指标之间存在相关关系。对任意两项生理数据，均执行上述相关关系计算，找出所有存在相关关系的生理指标。在本发明的一个实施例中，经相关关系检验，存在相关关系的生理指标有：VE和VO₂,VE和VCO₂,VO₂和VCO₂,EEm和VCO₂,EEm和VO₂。

步骤S4，采用格兰杰因果分析方法检验第三步得到的存在相关关系的其中两组数据之间是否具有因果关系；实施方式如下：

当两个变量之间存在相关关系时，才进一步检验它们之间的因果关系。因果关系的检验按照如下方式进行：假定X(t)和Y(t)为步骤S3中所得的两个具有相关关系的随机平稳时间序列，则变量Y(t)可以用它的自回归方程和它与X(t)之间的相互作用来表达，表达式为：

其中，p,q分别为X(t)和Y(t)回归方程中的滞后阶数；i和j为滞后项数，i的取值为1到p之间的整数，j的取值为1到q之间的整数；X(t-i)为X(t)的第i个滞后项，Y(t-j)为Y(t)的第j个滞后项，α_YX(i)和α_YY(i)为Y(t)的系数估计值，e(t)为Y(t)的白噪声。如果α_YX(i)(i＝1,2,...p)在统计学上整体显著不为零，则说明X(t)是Y(t)的格兰杰原因。如果α_YX(i)(i＝1,2,...p)在统计学上整体显著等于零，则说明X(t)和Y(t)不存在因果关系。在本发明的一个实施例中，p和q均取4阶滞后，因果关系检验借助Eviews10软件进行。

步骤S5，对于步骤S3中已确定的具有相关关系的数据，循环运用步骤S4中的格兰杰因果检验方法，检验任意两组数据之间是否具有因果关系，直到所有数据之间的因果关系均被检验。

步骤S6，对于步骤S5中已经检验出的存在因果关系的生理指标，绘制各项指标之间的因果关系网络图。用箭头表示生理指标之间的因果关系方向，箭头由“原因”变量指向“结果”变量。本发明的一个实施例得到的因果关系网络图将在图2中介绍。

请参照图2，在本发明的实施例中，经相关关系检验，存在相关关系的生理指标有：VE和VO₂,VE和VCO₂,VO₂和VCO₂,EEm和VCO₂,EEm和VO₂。进一步对上述5对生理指标进行格兰杰因果关系检验，得到图2所示结果。即：VE是VO₂的格兰杰原因，VCO₂是VE的格兰杰原因，VCO₂是VO₂的格兰杰原因，VCO₂是EEm的格兰杰原因，VO₂是EEm的格兰杰原因，EEm也是VO₂的格兰杰原因。

本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。