一种个体化的生理状态监测分析方法和设备

摘要

本申请一种个体化的生理状态监测分析方法，其包括：在指标数据库G中，找到与输入向量最为接近的n个指标向量，该n个指标向量构成矩阵D；利用找到优化参数w；利用优化参数w求得估计向量Vest，Vest＝Dw；对于矩阵D的每个指标向量xj，找到与其最为接近的n个指标向量并构成矩阵Dj，利用找到优化参数wj；利用优化参数wj求得估计向量Vjest，Vjest＝Djwj；计算该输入向量与其估计向量之间的差值，作为第一差值Vre；计算矩阵D中的每个指标向量xj与其估计向量Vjest的差值，作为第二差值；全部第二差值构成差值库Gre；根据第一差值Vre在差值库Gre中的分布情况，计算状态变化指数。

说明书

技术领域

本发明涉及人体生理参数的监测和分析技术，尤其涉及一种通过可穿戴设备获取建册对象的基本生理信号，对连续生理信号进行学习，进而对生理状态变化进行量化分析的方法和设备。

背景技术

监护技术在上世纪就已经出现了。从早期的床旁监护，到移动监护，进而到现在的穿戴式生理监护。在生理监护领域，有一个重要的概念就是如何刻画生理状态的变化，对这种生理状态的变化进行量化。传统的方法往往是依赖人为经验来观察某一项或某几项生理指标绝对数值的变化，当观察的指标超过一定的阈值之后，观察者就会认为该状态发生了异常，从而做出下一步的应变措施。这种主观经验方法目前在绝大多数工业和产业界都是主要的方法。但随着系统的日趋复杂，这种方法已不再能够满足任务的实际需要，其原因在于很多问题中阈值的确立伴随着很强的个人主观性，对于复杂系统，个人的经验很多时候是不准的；另一方面，由于系统的复杂性，当观察的指标超过实际设定的阈值时，系统已经处于崩溃的状态或奔溃的边缘，这种主观经验方法无法起到一个早期报警的作用，从而限制了其在很多实际应用的效果。对于人体的状态监测来说，更是如此。人体是一个及其复杂的系统，心脏、肝、肺等器官之间有着千丝万缕的联系，很多时候如果只观察某一个器官的指标很难看出人体状态的变化，而且人与人之间的个体化差异十分大，对一个人的定下的阈值对另一个人可能不再适用，这又极大复杂了病理监护人员的工作。

发明内容

鉴于上述问题，本申请旨在提出一种个体化的生理状态监测分析方法，其不是基于传统的阈值方法的生理参数监护分析技术，而是通过被监测对象的生理时间序列纵向对比分析，识别生理状态变化，从而实现更加灵敏的、特异性的、个体化的生理状态监测。

本申请的个体化的生理状态监测分析方法，其包括：

将在同一时间窗的k个生理参数的数值构成的向量作为指标向量，其中k为大于1的自然数；

用该个体在正常状态下的多个指标向量构成指标数据库G；

将该个体的当前时间窗的指标向量作为输入向量V

在指标数据库G中，通过欧氏距离或者核函数来找到与所述输入向量最为接近的n个指标向量，该n个指标向量构成矩阵D；利用

找到优化参数w；其中e为n*1的噪声向量，λ

利用优化参数w求得估计向量V

对于矩阵D的每个指标向量x

找到优化参数w

利用优化参数w

计算该输入向量与其估计向量之间的差值，作为第一差值V

计算矩阵D中的每个指标向量x

根据第一差值V

优选地，计算差值库G

计算差值库G

其中，x

计算自由度为k的卡方分布0.95分位点的λ值：

其中f(x)为自由度为k的卡方分布密度函数；

计算得到近似真实的协方差矩阵：

∑＝∑

计算马氏距离h

h

然后通过自由度为k的卡方分布累积分布函数来计算SCI的具体数值：

优选地，所述指标数据库存储于计算平台；所述预定个数的指标向量的查找、所述估计向量V

优选地，所述计算平台位于本地，穿戴在个体身上；

所述k个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得；

所述k个生理参数传输至所述计算平台；所述计算平台根据接收的所述k个生理参数而获得所述指标向量。

优选地，所述计算平台位于远程；

所述k个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得；

所述k个生理参数通过无线传输发送到所述计算平台；所述计算平台根据接收的所述k个生理参数而获得所述指标向量。

优选地，所述状态变化指数SCI在0-1之间，数值越高表示偏离正常状态越远。

优选地，如果状态变化指数SCI大于0且小于预定值，则将该输入向量作为指标数据库的一个指标向量而纳入指标数据库。

本申请的个体化的生理状态监测分析设备，其包括：计算平台；

计算平台包括指标数据库单元、估计单元、差值单元、状态变化指数单元；

计算平台被配置，用该个体的同一时间窗的k个生理参数的数值构成向量，作为指标向量；

指标数据库单元用该个体在正常状态下的多个指标向量构成指标数据库G；

估计单元用于计算与输入到估计单元的指标向量对应的估计向量；

差值单元用于计算两个指标向量之间的差值；

状态变化指数单元用于计算反映个体当前的状态偏离其正常状态的程度的状态变化指数SCI；

计算平台将个体的当前时间窗的指标向量作为输入向量V

估计单元利用

找到优化参数w；其中e为n*1的噪声向量，λ

估计单元利用优化参数w求得估计向量V

对于矩阵D的每个指标向量x

估计单元利用

找到优化参数w

估计单元利用优化参数w

差值单元计算该输入向量与其估计向量之间的差值，作为第一差值V

差值单元计算矩阵D中的每个指标向量x

状态变化指数单元根据第一差值V

优选地，状态变化指数单元计算差值库G

状态变化指数单元计算差值库G

其中，x

状态变化指数单元计算自由度为k的卡方分布0.95分位点的λ值：

其中f(x)为自由度为k的卡方分布密度函数；

状态变化指数单元计算得到近似真实的协方差矩阵：

∑＝∑

状态变化指数单元计算马氏距离h

h

然后状态变化指数单元通过自由度为k的卡方分布累积分布函数来计算SCI的具体数值：

优选地，所述设备被穿戴在该个体身上；所述k个生理参数由该个体穿戴的生理参数传感器感测得到的生理信号经处理后获得。

优选地，所述计算平台位于远程；所述k个生理参数由该个体穿戴的生理参数传感器感测得到的生理信号经处理后获得；所述k个生理参数通过无线方式传输给所述计算平台。

与以往的主观经验阈值方法相比，本申请从考虑单个个体前后状态的相似性出发，考察多个指标之间的协同相似联系，从而避免了个体之间差异性大的问题。另一方面，即使某些单项指标还处在正常的波动范围里面，也能够更加早期的预测出生理状态变化情况。

附图说明

图1为可穿戴系统采集和发送数据流程图；

图2为可穿戴系统及其采集的生理信号；

图3为心率、呼吸率和三轴加速度信号及30s滑动时间窗；

图4为生理状态监测分析整体流程图；

图5a为实施例一的模拟高原低氧环境实验结果；

图5b为实施例二的模拟高原低氧环境实验结果；

图6为实施例三的突发房颤抢救患者案例。

具体实施方式

下面，结合附图对本申请进行详细说明。

首先通过让被监测对象穿戴胸衣进行数据采集，穿戴胸衣包括带有心电导联接口、呼吸传感器、加速度传感器和血氧饱和度采集器组成，分别可以采集心电信号、呼吸信号、体位体动信号和血氧饱和度，并把采集到的信号集成到信号采集盒。数据发送模块会把采集到的数据中央计算平台，中央计算平台可以是本地的，也可以是远程的。信号流程图如图1所示。

可穿戴系统可以为一背心或者胸带，采集心电、呼吸和三轴加速度传感器信号，心电信号和呼吸信号经过进一步处理，可以获得心率和呼吸率时间序列。可穿戴系统及其采集的原始生理信号(心电、呼吸、三轴加速度)以及心率和呼吸率时间序列如图2所示。

对原始的心率、呼吸率信号，通过平滑滤波依次作用于各项原始信号指标，并从中提取到佩戴者心率指标、呼吸率指标、和三轴加速度指标。采用30s时间窗内(图3中黑色矩形框)的心率、呼吸率中位数组成二维向量，表征测试者当前的生命体征，通常选取被监测者健康、稳定状态下的向量来构建生理指标状态向量库。

当数据传送到中央计算平台后，中央计算平台会为每个被监测者提供单独的计算资源，通过监测一段时间的生理参数，建立单独的生理指标状态库，用G表示。完成G的建立之后，即可开展个体化的生理状态监测分析，设输入指标向量为V

刻画两个向量x

其中核函数

在寻找最优估计向量V

方法一，得到矩阵D后，运用优化理论，将寻找V

w≥0

其中，矩阵E使用frobenius范数，D为与V

方法二，针对参数w要进行进一步的正则化，比如对于系数w的稀疏化，我们会将寻找最相似向量问题转换为下列的优化问题：

w≥0

其中，矩阵E使用frobenius范数，系数w向量使用L1范数，λ

方法三，也可以简化噪声部分，仅仅考虑噪声向量e：

w≥0

其中，e为形为15*1的噪声向量，系数w的范数使用L2范数。

在采用上述方法得到系数w后,可以得到与V

v

之后考虑V

v

通过V

除了观察V

进一步的，基于多个生理参数的变化，可以状态变化指数(SCI)进行估计。

由于在模板库G样本处理时，为消除异常向量对结果的影响，我们只保留了95％分位点对应的等概率范围内的样本点，所以由此计算的残差不再是随机样本，估计到的协方差矩阵的元素绝对值偏小，所以需要推导出近似的真实协方差矩阵。

通过样本G的残差G

计算其马氏距离平方最大值：

其中，x

其中f(x)为自由度为样本特征维度的卡方分布密度函数。

计算得到近似真实的协方差矩阵：

∑＝∑

根据当前计算的残差V

首先计算马氏距离：

h

之后通过自由度为样本特征维度k的卡方分布累积分布函数来计算SCI的具体数值：

这时，可以根据V

该计算平台的一大优势在于当新的状态指标数据V

首先，为了验证算法的准确性与可行性，针对健康人进行了模拟高原实验，实验采用测试者日常生活中的一天的生理参数作为状态向量库G，在模拟高原低氧环境中，测试者会在模拟海拔四千米低氧的环境下分别做静止看书、快步走、高抬腿三个动作，尽量保持与前一天接近的活动强度，验证算法是否能够准确检测测试者当前的生理体征是否偏离状态向量库G。

实施例一

图5a为一例测试者的模拟高原低氧实验结果图，其实施的是第一种方法，图中加的血氧信号为了更好地观察测试者所处的低氧环境。

在实验刚开始时，测试者还处在正常海拔下，当前的状态向量为V

最终计算出来的最优估计值V

差值V

最终计算出来的最优估计值V

差值V

实施例二

图5b为一例测试者的模拟高原低氧实验结果图，其中图5b实施的是第二种方法，图中加的血氧信号为了更好地观察测试者所处的低氧环境。

在实验刚开始时，测试者还处在正常海拔下，当前的状态向量为V

最终计算出来的最优估计值V

差值V

最终计算出来的最优估计值V

差值V

实施例三

图6为一突发严重房颤进行抢救的患者案例。图6中实施的是方法三。选取患者处于相对平稳状态下的一天的生理数据构建生理状态向量库G，图6中(a)为经过抢救治疗后生命体征平稳的某个上午的数据，比如上午08点52分左右，当前的状态向量为：V

最终计算出的最优化估计值V

计算的最优估计向量：V

本申请中，通过被监测对象的生理时间序列纵向对比分析，识别生理状态变化，从而实现更加灵敏的、特异性的、个体化的生理状态监测。

本申请中计算平台可以以单片机、DSP、计算机等实现，指标数据库单元、估计单元、差值单元、状态变化指数单元可以是计算平台上通过程序而实现的功能模块。

除非另有定义，本申请中使用的所有技术和/或科学术语具有与由本发明所涉及的领域的普通技术人员通常理解的相同含义。本申请中提到的材料、方法和实施例仅为说明性的，而非限制性的。

虽然已结合具体实施方式对本发明进行了描述，在本申请的发明主旨下，本领域的技术人员可以进行适当的替换、修改和变化，这种替换、修改和变化仍属于本申请的保护范围。