CN107149480A - 一种基于驾驶人心电rr间期的驾驶疲劳判别方法 - Google Patents

Abstract

本发明涉及一种基于驾驶人心电RR间期的驾驶疲劳判别方法，该方法包括：驾驶人心电信号的采集及心电RR间期的提取；根据基于心电RR间期序列的平稳性模型计算对应序列结构残差；根据一阶自相关模型计算自回归残差序列对应于各个时刻的自回归残差；根据条件方差的GARCH(1,1)模型计算条件方差序列中对应于各个时刻的条件方差；根据对应于当前时刻的条件方差判别驾驶人的疲劳等级是否变化。本发明解决了现有技术中驾驶疲劳判别主观性过强，以及因个体差异导致的判别标准不同的问题，具有较强的客观性强及普遍适用性。

Description

一种基于驾驶人心电RR间期的驾驶疲劳判别方法

技术领域

本发明属于交通安全领域，涉及一种基于驾驶人心电RR间期的驾驶疲劳判别方法。

背景技术

当驾驶员处于疲劳状态时极易出现观察、判断以及操纵上的失误，进而导致交通事故的发生。据美国国家交通安全协会统计，28％的交通事故由驾驶疲劳造成。因此，及时辨识驾驶员疲劳状态并进行预警对于预防交通事故具有重要意义。

为了能够及时判别驾驶疲劳，应利用一些驾驶人自身的生理指标及车辆运动等参数对驾驶疲劳进行检测并预警。但现有的驾驶疲劳判别存在几点不足：

1、驾驶疲劳判别指标阈值界定的主观性强。因为驾驶疲劳是根据驾驶人的主观感受来界定的，如果通过确定判别指标阈值来界定驾驶人是否达到疲劳，主观性太强，不具有普遍适用性。因此疲劳阈值的界定一直是驾驶疲劳判别的难点。

2、因驾驶人个体差异性导致的疲劳判别的标准不同。驾驶疲劳是根据驾驶人的主观感受来界定的，因驾驶人的年龄、性别以及驾龄等因素影响导致驾驶疲劳判别的评判标准难以确定。因此由驾驶人个体差异性导致的疲劳判别误差较大问题一直是驾驶疲劳判别的难点。

发明内容

本发明要解决的技术问题是提供一种基于驾驶人心电RR间期的驾驶疲劳判别方法，该方法能够避免驾驶疲劳判别主观性过强，以及因个体差异导致的判别标准不同的缺陷，具有较强的客观性及普遍适用性。

为了解决上述技术问题，本发明的基于驾驶人心电RR间期的驾驶疲劳判别方法包括下述步骤：

步骤一、利用生理记录仪实时采集驾驶人心电信号并进行预处理获得驾驶人心电RR间期的时间序列{x₁,x₂,......,x_T}，T为初始时刻到当前时刻采集的心率RR间期样本总数；

步骤二、根据基于心电RR间期序列的平稳性模型公式(1)计算对应序列结构残差u_t；其中1≤t≤T；

x_t＝c+u_t (1)

式中：x_t—对应于时刻t的驾驶人心电RR间期；

c—初始时刻到当前时刻T驾驶人心电RR间期时间序列均值，即

步骤三、根据一阶自相关模型公式(2)计算自回归残差序列对应于时刻t的自回归残差e_t；

u_t＝ρu_t-1+e_t (2)

其中u_t-1为对应于时刻t-1的序列结构残差；

步骤四、根据最小二乘法公式(3)估计模型参数ρ，ρ取值为[-1,1]；若u_t和u_t-1正相关，则ρ为正，若u_t和u_t-1负相关则ρ为负；

步骤五、根据条件方差的GARCH(1,1)模型公式(4)计算条件方差序列中对应于时刻t的条件方差h_t；

其中h_t-1是对应于上一时刻t-1的条件方差，h_t的初始值h₁＝0；e_t-1是对应于时刻t-1的自回归残差；c₁、c₂、c₃是模型参数，c₁≥0,c₂≥0,c₃≥0，并且c₁、c₂、c₃的取值应当使得f(z₁)×f(z₂)×......f(z_T)最大；z_T为对应于当前时刻T的标准化残差，f(z_t)为正态分布函数；

步骤六、根据对应于时刻t的条件方差h_t判别驾驶人的疲劳等级是否变化：

设v_a表示条件方差h_t的平均值，σ₀代表平均值v_a的标准差，则

设置初始计数值level＝1；当h_T-v_a＞3σ₀时，将计数值level加1；若level＝1，则表示驾驶人疲劳等级为清醒，level＝2，则表示驾驶人疲劳等级为轻度疲劳，计数值level＝3则表示驾驶人疲劳等级为重度疲劳；计数值level＝4则表示驾驶人疲劳等级为困倦4个等级。

本发明的有益效果：

1、解决了通过界定指标阈值来判别驾驶疲劳的主观性过强的问题

因为驾驶疲劳是根据驾驶人的主观感受来界定的，如果通过确定判别指标阈值来界定驾驶人是否达到疲劳，主观性太强，不具有普遍适用性。本发明利用驾驶人心电RR间期序列的条件方差有效地拟合了序列的集群效应波动特征，基于波动特征判别驾驶疲劳，解决了通过界定指标阈值来判别驾驶疲劳的主观性过强的问题。

2、解决了因驾驶人个体差异性导致的疲劳判别的标准不同的问题

驾驶疲劳是根据驾驶人的主观感受来界定的，因驾驶人的年龄、性别以及驾龄等因素影响导致驾驶疲劳判别的评判标准难以确定。本发明通过建立基于心电RR间期的时间序列判别模型，利用RR间期序列的波动性表征驾驶人自身状态变化进而判别驾驶疲劳，避免了利用同一指标阈值判别不同驾驶人的疲劳。解决了因驾驶人个体差异性导致的疲劳判别标准不同的问题。

附图说明

下面结合附图和具体实施方式对本发明作进一步详细说明。

图1是本发明的基于心电RR间期的驾驶人疲劳判别方法的流程图。

具体实施方式

如图1所示，本发明的基于心电RR间期的驾驶人疲劳判别方法具体包括下述步骤：

步骤一、数据获取及预处理:依次流程包括搭建实车道路试验平台，该平台包括美国Biopac生理记录仪及与其连接的计算机；采集驾驶人心电信号；利用Biopac生理记录仪的软件分析模块中选择hemodynamics-ECG interval Extraction，对驾驶人心电信号进行预处理获取驾驶人心电RR间期的时间序列{x₁,x₂,......,x_T}，T为初始时刻到当前时刻采集的心率RR间期样本总数；

步骤二、设置初始计数值level＝1；根据基于心电RR间期序列的平稳性模型公式(1)计算对应于时刻t的序列结构残差u_t，其中1≤t≤T；

x_t＝c+u_t (1)

式中：x_t—对应于时刻t的驾驶人心电RR间期；

c—初始时刻到当前时刻T驾驶人心电RR间期时间序列均值；心电RR间期时间序列围绕固定值c随机波动；

所述初始时刻到当前时刻驾驶人心电RR间期时间序列均值c可以通过计算均值得到；

x_t是对应于任一时刻t的驾驶人心电RR间期，T-初始时刻到当前时刻心率RR间期样本总数；

步骤三、根据一阶自相关模型公式(2)计算自回归残差序列对应于当前时刻t的自回归残差e_t；

u_t＝ρu_t-1+e_t (2)

其中u_t-1为对应于上一时刻t-1的序列结构残差；

步骤四、根据最小二乘法公式(3)估计模型参数ρ，ρ取值为[-1,1]；若u_t和u_t-1正相关，则ρ为正，若u_t和u_t-1负相关则ρ为负；

步骤五、根据条件方差的GARCH(1,1)模型公式(4)计算条件方差序列中对应于当前时刻t的条件方差h_t；

其中h_t-1是对应于上一时刻t-1的条件方差，h_t的初始值h₁＝0；e_t-1是对应于上一时刻t-1的自回归残差；c₁、c₂、c₃是模型参数，c₁≥0,c₂≥0,c₃≥0,并且c₁、c₂、c₃的取值应当使得f(z₁)×f(z₂)×......f(z_T)最大；z_T为对应于当前时刻T的标准化残差，f(z_t)为正态分布函数；

步骤六、实时或者每隔15min左右根据对应于当前时刻t的条件方差h_t判别驾驶人的疲劳等级是否变化；

设v_a表示条件方差h_t的平均值，σ₀代表平均值v_a的标准差，则

若h_t-v_a＞3σ₀，表明心率R-R间期序列的残差序列发生了巨大的变化，进而说明驾驶人的疲劳程度发生了显著变化，则将计数值level加1；若level＝1，则表示驾驶人疲劳等级为清醒，level＝2，则表示驾驶人疲劳等级为轻度疲劳，计数值level＝3则表示驾驶人疲劳等级为重度疲劳；计数值level＝4则表示驾驶人疲劳等级为困倦4个等级。

Claims (1)

1.一种基于驾驶人心电RR间期的驾驶疲劳判别方法，其特征在于包括下述步骤：

步骤一、利用生理记录仪实时采集驾驶人心电信号并进行预处理获得驾驶人心电RR间期的时间序列{x₁,x₂,......,x_T}，T为初始时刻到当前时刻采集的心率RR间期样本总数；

步骤二、根据基于心电RR间期序列的平稳性模型公式(1)计算对应序列结构残差u_t；其中1≤t≤T；

x_t＝c+u_t (1)

式中：x_t—对应于时刻t的驾驶人心电RR间期；

c—初始时刻到当前时刻T驾驶人心电RR间期时间序列均值，即

步骤三、根据一阶自相关模型公式(2)计算自回归残差序列对应于时刻t的自回归残差e_t；

u_t＝ρu_t-1+e_t (2)

其中u_t-1为对应于时刻t-1的序列结构残差；

步骤四、根据最小二乘法公式(3)估计模型参数ρ，ρ取值为[-1,1]；若u_t和u_t-1正相关，则ρ为正，若u_t和u_t-1负相关则ρ为负；

<mrow> <mi>&amp;rho;</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <msub> <mi>u</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>u</mi> <mi>t</mi> </msub> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <msubsup> <mi>u</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

步骤五、根据条件方差的GARCH(1,1)模型公式(4)计算条件方差序列中对应于时刻t的条件方差h_t；

<mrow> <msub> <mi>h</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msubsup> <mi>e</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <msub> <mi>h</mi> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

其中h_t-1是对应于上一时刻t-1的条件方差，h_t的初始值h₁＝0；e_t-1是对应于时刻t-1的自回归残差；c₁、c₂、c₃是模型参数，c₁≥0,c₂≥0,c₃≥0，并且c₁、c₂、c₃的取值应当使得f(z₁)×f(z₂)×......f(z_T)最大；z_T为对应于当前时刻T的标准化残差，f(z_t)为正态分布函数；

<mrow> <msub> <mi>e</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>z</mi> <mi>t</mi> </msub> <msqrt> <msub> <mi>h</mi> <mi>t</mi> </msub> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

步骤六、根据对应于时刻t的条件方差h_t判别驾驶人的疲劳等级是否变化：

设v_a表示条件方差h_t的平均值，σ₀代表平均值v_a的标准差，则

<mrow> <msub> <mi>v</mi> <mi>a</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <msub> <mi>h</mi> <mi>t</mi> </msub> </mrow> <mi>T</mi> </mfrac> <mo>,</mo> <msub> <mi>&amp;sigma;</mi> <mn>0</mn> </msub> <mo>=</mo> <msqrt> <msub> <mi>v</mi> <mi>a</mi> </msub> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

设置初始计数值level＝1；当h_T-v_a＞3σ₀时，将计数值level加1；若level＝1，则表示驾驶人疲劳等级为清醒，level＝2，则表示驾驶人疲劳等级为轻度疲劳，计数值level＝3则表示驾驶人疲劳等级为重度疲劳；计数值level＝4则表示驾驶人疲劳等级为困倦4个等级。