CN112971791B - 一种个体化的生理状态监测分析方法和设备 - Google Patents
一种个体化的生理状态监测分析方法和设备 Download PDFInfo
- Publication number
- CN112971791B CN112971791B CN202010205666.6A CN202010205666A CN112971791B CN 112971791 B CN112971791 B CN 112971791B CN 202010205666 A CN202010205666 A CN 202010205666A CN 112971791 B CN112971791 B CN 112971791B
- Authority
- CN
- China
- Prior art keywords
- index
- vector
- vectors
- physiological
- difference
- 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
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000035790 physiological processes and functions Effects 0.000 title abstract description 21
- 239000013598 vector Substances 0.000 claims abstract description 173
- 230000008859 change Effects 0.000 claims abstract description 36
- 238000004458 analytical method Methods 0.000 claims description 17
- 230000004962 physiological condition Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000000241 respiratory effect Effects 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 claims 6
- 210000002216 heart Anatomy 0.000 description 20
- 239000011159 matrix material Substances 0.000 description 15
- 230000029058 respiratory gaseous exchange Effects 0.000 description 13
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 6
- 239000008280 blood Substances 0.000 description 6
- 210000004369 blood Anatomy 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 229910052760 oxygen Inorganic materials 0.000 description 6
- 239000001301 oxygen Substances 0.000 description 6
- 206010021143 Hypoxia Diseases 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000007954 hypoxia Effects 0.000 description 5
- 206010003658 Atrial Fibrillation Diseases 0.000 description 4
- 208000008445 altitude sickness Diseases 0.000 description 4
- 230000036387 respiratory rate Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002565 electrocardiography Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
- A61B5/6804—Garments; Clothes
- A61B5/6805—Vests
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Pulmonology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Optics & Photonics (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本申请公开一种个体化的生理状态监测分析方法及设备,其通过在指标数据库中找到与所述输入向量最为接近的预定个数的指标向量,用该预定个数的指标向量计算该输入向量的估计向量;该预定个数的指标向量构成基准向量集;对于基准向量集中的每个指标向量,在指标数据库中找到与其最为接近的预定个数的指标向量,并用该预定个数的指标向量计算与其对应的估计向量;计算该输入向量与其估计向量之间的差值,作为第一差值;计算所述基准向量集中的每个指标向量与其估计向量的差值,作为第二差值;根据第一差值在全部第二差值中的分布情况,计算状态变化指数,通过状态变化指数反映个体当前的状态偏离其正常状态的程度。
Description
技术领域
本发明涉及人体生理参数的监测和分析技术,尤其涉及一种通过可穿戴设备获取建册对象的基本生理信号,对连续生理信号进行学习,进而对生理状态变化进行量化分析的方法和设备。
背景技术
监护技术在上世纪就已经出现了。从早期的床旁监护,到移动监护,进而到现在的穿戴式生理监护。在生理监护领域,有一个重要的概念就是如何刻画生理状态的变化,对这种生理状态的变化进行量化。传统的方法往往是依赖人为经验来观察某一项或某几项生理指标绝对数值的变化,当观察的指标超过一定的阈值之后,观察者就会认为该状态发生了异常,从而做出下一步的应变措施。这种主观经验方法目前在绝大多数工业和产业界都是主要的方法。但随着系统的日趋复杂,这种方法已不再能够满足任务的实际需要,其原因在于很多问题中阈值的确立伴随着很强的个人主观性,对于复杂系统,个人的经验很多时候是不准的;另一方面,由于系统的复杂性,当观察的指标超过实际设定的阈值时,系统已经处于崩溃的状态或奔溃的边缘,这种主观经验方法无法起到一个早期报警的作用,从而限制了其在很多实际应用的效果。对于人体的状态监测来说,更是如此。人体是一个及其复杂的系统,心脏、肝、肺等器官之间有着千丝万缕的联系,很多时候如果只观察某一个器官的指标很难看出人体状态的变化,而且人与人之间的个体化差异十分大,对一个人的定下的阈值对另一个人可能不再适用,这又极大复杂了病理监护人员的工作。
发明内容
鉴于上述问题,本申请旨在提出一种个体化的生理状态监测分析方法,其不是基于传统的阈值方法的生理参数监护分析技术,而是通过被监测对象的生理时间序列纵向对比分析,识别生理状态变化,从而实现更加灵敏的、特异性的、个体化的生理状态监测。
本申请的个体化的生理状态监测分析方法,其包括:
将在同一时间窗的多个生理参数的数值构成的向量作为指标向量;
用该个体在正常状态下的多个指标向量构成指标数据库;
将该个体的当前时间窗的指标向量作为输入向量;
在指标数据库中找到与所述输入向量最为接近的预定个数的指标向量,用该预定个数的指标向量计算该输入向量的估计向量;该预定个数的指标向量构成基准向量集;
对于基准向量集中的每个指标向量,在指标数据库中找到与其最为接近的预定个数的指标向量,并用该预定个数的指标向量计算与其对应的估计向量;
计算该输入向量与其估计向量之间的差值,作为第一差值;
计算所述基准向量集中的每个指标向量与其估计向量的差值,作为第二差值;
根据第一差值在全部第二差值中的分布情况,计算状态变化指数,通过状态变化指数反映个体当前的状态偏离其正常状态的程度。
优选地,所述多个生理参数中至少包括心电信号或呼吸信号。
优选地,所述多个生理参数还可以包括体位信号和体动信号。
优选地,所述指标数据库存储于计算平台;所述预定个数的指标向量的查找、所述估计向量的计算、所述第一差值的计算、所述第二差值的计算、所述第一差值在整个第二差值中的分布的计算、以及状态变化指数的计算均在该计算平台上进行。
优选地,所述计算平台位于本地,穿戴在个体身上;
所述多个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得;
所述多个生理参数传输至所述计算平台;所述计算平台根据接收的所述多个生理参数而获得所述指标向量。
优选地,所述计算平台位于远程;
所述多个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得;
所述多个生理参数通过无线传输发送到所述计算平台;所述计算平台根据接收的所述多个生理参数而获得所述指标向量。
优选地,所述状态变化指数在0-1之间,数值越高表示偏离正常状态越远。
优选地,如果状态变化指数大于0且小于预定值,则将该输入向量作为指标数据库的一个指标向量而纳入指标数据库。
本申请的个体化的生理状态监测分析设备,其包括:计算平台;
计算平台包括指标数据库单元、估计单元、差值单元、状态变化指数单元;
计算平台被配置,用该个体的同一时间窗的多个生理参数的数值构成向量,作为指标向量;
指标数据库单元用该个体在正常状态下的多个指标向量构成指标数据库;
估计单元用于计算与输入到估计单元的指标向量对应的估计向量;
差值单元用于计算两个指标向量之间的差值;
状态变化指数单元用于计算反映个体当前的状态偏离其正常状态的程度的状态变化指数;
计算平台将个体的当前时间窗的指标向量作为输入向量,根据该输入向量自指标数据库单元中选出与所述输入向量最为接近的预定个数的指标向量,估计单元根据该预定个数的指标向量计算出该输入向量的估计向量;该预定个数的指标向量构成基准向量集;对于基准向量集中的每个指标向量,计算平台自指标数据库单元中找到与其最为接近的预定个数的指标向量,估计单元根据该预定个数的指标向量计算与其对应的估计向量;
差值单元计算出该输入向量与其估计向量之间的差值,作为第一差值;差值单元计算出基准向量集中的每个指标向量与其估计向量的差值,作为第二差值;
状态变化指数单元根据第一差值在全部第二差值中的分布情况,计算状态变化指数。
优选地,所述设备被穿戴在该个体身上;所述多个生理参数由该个体穿戴的生理参数传感器感测得到的生理信号经处理后获得。
与以往的主观经验阈值方法相比,本申请从考虑单个个体前后状态的相似性出发,考察多个指标之间的协同相似联系,从而避免了个体之间差异性大的问题。另一方面,即使某些单项指标还处在正常的波动范围里面,也能够更加早期的预测出生理状态变化情况。
附图说明
图1为可穿戴系统采集和发送数据流程图;
图2为可穿戴系统及其采集的生理信号;
图3为心率、呼吸率和三轴加速度信号及30s滑动时间窗;
图4为生理状态监测分析整体流程图;
图5a为实施例一的模拟高原低氧环境实验结果;
图5b为实施例二的模拟高原低氧环境实验结果;
图6为实施例三的突发房颤抢救患者案例。
具体实施方式
下面,结合附图对本申请进行详细说明。
首先通过让被监测对象穿戴胸衣进行数据采集,穿戴胸衣包括带有心电导联接口、呼吸传感器、加速度传感器和血氧饱和度采集器组成,分别可以采集心电信号、呼吸信号、体位体动信号和血氧饱和度,并把采集到的信号集成到信号采集盒。数据发送模块会把采集到的数据中央计算平台,中央计算平台可以是本地的,也可以是远程的。信号流程图如图1所示。
可穿戴系统可以为一背心或者胸带,采集心电、呼吸和三轴加速度传感器信号,心电信号和呼吸信号经过进一步处理,可以获得心率和呼吸率时间序列。可穿戴系统及其采集的原始生理信号(心电、呼吸、三轴加速度)以及心率和呼吸率时间序列如图2所示。
对原始的心率、呼吸率信号,通过平滑滤波依次作用于各项原始信号指标,并从中提取到佩戴者心率指标、呼吸率指标、和三轴加速度指标。采用30s时间窗内(图3中黑色矩形框)的心率、呼吸率中位数组成二维向量,表征测试者当前的生命体征,通常选取被监测者健康、稳定状态下的向量来构建生理指标状态向量库。
当数据传送到中央计算平台后,中央计算平台会为每个被监测者提供单独的计算资源,通过监测一段时间的生理参数,建立单独的生理指标状态库,用G表示。完成G的建立之后,即可开展个体化的生理状态监测分析,设输入指标向量为Vin,给定Vin后,系统会从G中找出与Vin最相似的15个向量组成D矩阵,刻画相似程度可以用欧氏距离或核函数等等,这里采用表示核函数运算,
刻画两个向量xi与xj的相近程度由以下公式给出:
其中核函数可以高斯核、三角函数核等合法核函数。
在寻找最优估计向量Vest时,本申请采用了三种方法。
方法一,得到矩阵D后,运用优化理论,将寻找Vin最优估计向量Vest的问题转换成先求参数w的优化问题:
w≥0
其中,矩阵E使用frobenius范数,D为与Vin最相近的15个向量组成的矩阵,E为形为15*15的噪声矩阵,e为形为15*1的噪声向量,w为15*1的系数向量。
方法二,针对参数w要进行进一步的正则化,比如对于系数w的稀疏化,我们会将寻找最相似向量问题转换为下列的优化问题:
w≥0
其中,矩阵E使用frobenius范数,系数w向量使用L1范数,λ1为惩罚项系数。
方法三,也可以简化噪声部分,仅仅考虑噪声向量e:
w≥0
其中,e为形为15*1的噪声向量,系数w的范数使用L2范数。
在采用上述方法得到系数w后,可以得到与Vin的最优估计向量Vest:
vest=Dw
之后考虑Vest与Vin的差值Vre:
vre=vin-vest
通过Vre可以反映这两个向量之间相似度,来刻画状态改变与否,当输入的向量Vin与计算系统挑出来的向量Vest很相似的时候,认为当下的生理状态没有发生太大的变化,如果输入的向量与计算系统挑出来的向量差别较大时,可以认为当下的生理状态发生了变化。
除了观察Vre数值的大小来反映生理状态的变化,更为稳健的方法是观察Vre的分布。
进一步的,基于多个生理参数的变化,可以状态变化指数(SCI)进行估计。
由于在模板库G样本处理时,为消除异常向量对结果的影响,我们只保留了95%分位点对应的等概率范围内的样本点,所以由此计算的残差不再是随机样本,估计到的协方差矩阵的元素绝对值偏小,所以需要推导出近似的真实协方差矩阵。
通过样本G的残差Gre计算其均值与方差μ1,∑1。
计算其马氏距离平方最大值:
其中,xi为Gre中的样本残差向量。计算自由度为样本特征维度的卡方分布0.95分位点的λ值:
其中f(x)为自由度为样本特征维度的卡方分布密度函数。
计算得到近似真实的协方差矩阵:
∑=∑1·h/λ
根据当前计算的残差Vre在样本Gre残差分布(μ1,∑)中的情况确定状态变化指数的大小,越接近与1表示偏离状态指标库G越远。
首先计算马氏距离:
hre=(Vre-μ1)T∑-1(Vre-μ1)
之后通过自由度为样本特征维度k的卡方分布累积分布函数来计算SCI的具体数值:
这时,可以根据Vre偏离的情况,选择进一步将Vin添加到生理指标状态向量库中,作为系统出现的一种新状态,来丰富状态指标库G。
该计算平台的一大优势在于当新的状态指标数据Vin来临时,可以进行定期更新,从而能够自动长时间监测,并及时作出早期预警。
首先,为了验证算法的准确性与可行性,针对健康人进行了模拟高原实验,实验采用测试者日常生活中的一天的生理参数作为状态向量库G,在模拟高原低氧环境中,测试者会在模拟海拔四千米低氧的环境下分别做静止看书、快步走、高抬腿三个动作,尽量保持与前一天接近的活动强度,验证算法是否能够准确检测测试者当前的生理体征是否偏离状态向量库G。
实施例一
图5a为一例测试者的模拟高原低氧实验结果图,其实施的是第一种方法,图中加的血氧信号为了更好地观察测试者所处的低氧环境。
在实验刚开始时,测试者还处在正常海拔下,当前的状态向量为Vin=[73.2602,25.0000],经过相似度比较,从G中取出15个向量组成矩阵D:
心率 | 呼吸率 |
73.2602 | 25 |
73.7101 | 25 |
72.3765 | 25 |
74.4417 | 25 |
74.4417 | 25 |
72.0289 | 25 |
71.8572 | 25 |
71.8564 | 25 |
71.6847 | 25 |
71.3445 | 25 |
75.2824 | 25 |
73.3496 | 26 |
70.9220 | 25 |
72.9931 | 24 |
73.6206 | 26 |
最终计算出来的最优估计值Vest=[73.9942,25.0637],
差值Vre=[0.3160,0.0637],计算Vre在Gre分布中的SCI指数为0.1584,偏离状态向量库G较小。相反,我们取高原实验开始后的向量:Vin=[92.7359,22.0000],对应的D矩阵为:
心率 | 呼吸率 |
76.9231 | 23 |
79.1557 | 24 |
77.7202 | 23.5 |
78.7402 | 24 |
76.4332 | 23 |
80.0000 | 25 |
77.6208 | 24 |
74.0741 | 22.5 |
72.8160 | 22 |
72.5523 | 22 |
81.3032 | 26 |
76.6295 | 24 |
71.6855 | 22 |
78.3290 | 25 |
73.7101 | 21 |
最终计算出来的最优估计值Vest=[77.7888,24.0335],
差值Vre=[14.9471,2.0335],计算后的SCI指数达到0.98,严重偏离状态向量库G。SCI指数介于0-1之间,数值越高表明偏离日常状态越远。从图5a中可以看出,测试者在低氧环境下生理状态发生改变(血氧饱和度降低),SCI指数也明显增大,能够有效表征被测试者生理状态变化。
实施例二
图5b为一例测试者的模拟高原低氧实验结果图,其中图5b实施的是第二种方法,图中加的血氧信号为了更好地观察测试者所处的低氧环境。
在实验刚开始时,测试者还处在正常海拔下,当前的状态向量为Vin=[68.1818,22.0000],经过相似度比较,从G中取出15个向量组成矩阵D:
心率 | 呼吸率 |
68.3766 | 22 |
67.7966 | 22 |
68.9655 | 22 |
67.2316 | 22 |
67.2274 | 22 |
69.1648 | 22 |
67.0391 | 22 |
67.0391 | 22 |
68.1818 | 23 |
68.1815 | 21 |
68.3766 | 23 |
65.9341 | 22 |
67.7966 | 21 |
67.7966 | 21 |
65.9341 | 22 |
最终计算出来的最优估计值Vest=[67.8312,21.9314],
差值Vre=[0.3505,0.0685],计算Vre在Gre分布中的SCI指数为0.0397,偏离状态向量库G较小。相反,我们取高原实验开始后的向量:Vin=[96.6186,16],对应的D矩阵为:
心率 | 呼吸率 |
69.9714 | 19 |
75.4717 | 22 |
68.9655 | 19 |
72.2892 | 21 |
76.9231 | 23.5 |
73.1707 | 22 |
67.0391 | 19 |
67.0391 | 19 |
71.0059 | 21 |
71.0059 | 21 |
68.5714 | 20 |
72.5082 | 22 |
74.5342 | 23 |
70.3818 | 21 |
72.2918 | 22 |
最终计算出来的最优估计值Vest=[75.1937,22.7126],
差值Vre=[21,4248,6.7126],计算后的SCI指数达到0.99,严重偏离状态向量库G。SCI指数介于0-1之间,数值越高表明偏离日常状态越远。从图5b中可以看出,测试者在低氧环境下生理状态发生改变(血氧饱和度降低),SCI指数也明显增大,能够有效表征被测试者生理状态变化。
实施例三
图6为一突发严重房颤进行抢救的患者案例。图6中实施的是方法三。选取患者处于相对平稳状态下的一天的生理数据构建生理状态向量库G,图6中(a)为经过抢救治疗后生命体征平稳的某个上午的数据,比如上午08点52分左右,当前的状态向量为:Vin=[51.2820,22],经过相似度比较,从G中取出15个向量组成矩阵D:
心率 | 呼吸率 |
51.3921 | 22 |
51.3921 | 22 |
51.3921 | 22 |
51.0638 | 22 |
51.0638 | 22 |
50.9556 | 22 |
50.8475 | 22 |
51.7241 | 22 |
50.6329 | 22 |
51.2821 | 21 |
52.0610 | 22 |
50.4202 | 22 |
50.4202 | 22 |
51.1729 | 23 |
52.1739 | 22 |
最终计算出的最优化估计值Vest=[51.2051,21.9986],差值向量为:Vre=[0.0769,0.0014],计算其在Gre分布:中的SCI值为0.0853,从图中可见其SCI值绝大部分处在正常水平。图6中(b)为患者发生房颤当天采集的数据,大概在下午19点40分时,当前的输入向量为Vin=[72.5,12.5],对应的矩阵D为:
心率 | 吁吸率 |
51.8361 | 14 |
50.6329 | 13 |
51.0638 | 14 |
52.4017 | 16 |
50.0000 | 13 |
52.2878 | 16 |
50.6329 | 14 |
50.6329 | 14 |
49.7925 | 13 |
53.5725 | 18 |
49.5868 | 13 |
51.8361 | 16 |
50.2092 | 14 |
50.2092 | 14 |
50.2092 | 14 |
计算的最优估计向量:Vest=[53.5725,18],差值Vre=[18.9275,5.5],计算后的SCI高达0.99,根据医疗记录显示,该患者确实从下午开始就感觉身体不适,在21点50左右才发生房颤,但SCI与正常水平相比,从19点40分开始,其数值一直处在较高的水平。
本申请中,通过被监测对象的生理时间序列纵向对比分析,识别生理状态变化,从而实现更加灵敏的、特异性的、个体化的生理状态监测。
本申请中计算平台可以以单片机、DSP、计算机等实现,指标数据库单元、估计单元、差值单元、状态变化指数单元可以是计算平台上通过程序而实现的功能模块。
除非另有定义,本申请中使用的所有技术和/或科学术语具有与由本发明所涉及的领域的普通技术人员通常理解的相同含义。本申请中提到的材料、方法和实施例仅为说明性的,而非限制性的。
虽然已结合具体实施方式对本发明进行了描述,在本申请的发明主旨下,本领域的技术人员可以进行适当的替换、修改和变化,这种替换、修改和变化仍属于本申请的保护范围。
Claims (11)
1.一种个体化的生理状态监测分析方法,其包括:
将在同一时间窗的多个生理参数的数值构成的向量作为指标向量;
用该个体在正常状态下的多个指标向量构成指标数据库;
将该个体的当前时间窗的指标向量作为输入向量;
在指标数据库中找到与所述输入向量最为接近的预定个数的指标向量,用该预定个数的指标向量计算该输入向量的估计向量;该预定个数的指标向量构成基准向量集;
对于基准向量集中的每个指标向量,在指标数据库中找到与其最为接近的预定个数的指标向量,并用该预定个数的指标向量计算与其对应的估计向量;
计算该输入向量与其估计向量之间的差值,作为第一差值;
计算所述基准向量集中的每个指标向量与其估计向量的差值,作为第二差值;
根据第一差值在全部第二差值中的卡方分布累积分布情况,计算状态变化指数,通过状态变化指数反映个体当前的状态偏离其正常状态的程度;其中,所述卡方分布累积分布的自由度为所述生理参数的个数,所述卡方分布累积分布的变量为根据第一差值和全部第二差值计算得到的马氏距离。
2.如权利要求1所述的个体化的生理状态监测分析方法,其特征在于:
所述多个生理参数中至少包括心电信号和呼吸信号。
3.如权利要求2所述的个体化的生理状态监测分析方法,其特征在于:
所述多个生理参数还包括体位信号和体动信号。
4.如权利要求1所述的个体化的生理状态监测分析方法,其特征在于:
所述指标数据库存储于计算平台;所述预定个数的指标向量的查找、所述估计向量的计算、所述第一差值的计算、所述第二差值的计算、所述第一差值在全部第二差值中的分布的计算、以及状态变化指数的计算均在该计算平台上进行。
5.如权利要求4所述的个体化的生理状态监测分析方法,其特征在于:
所述计算平台位于本地,穿戴在个体身上;
所述多个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得;
所述多个生理参数传输至所述计算平台;所述计算平台根据接收的所述多个生理参数而获得所述指标向量。
6.如权利要求4所述的个体化的生理状态监测分析方法,其特征在于:
所述计算平台位于远程;
所述多个生理参数由该个体所穿戴的生理参数传感器感测得到的生理信号经处理后获得;
所述多个生理参数通过无线传输发送到所述计算平台;所述计算平台根据接收的所述多个生理参数而获得所述指标向量。
7.如权利要求1所述的个体化的生理状态监测分析方法,其特征在于:
所述状态变化指数在0-1之间,数值越高表示偏离正常状态越远。
8.如权利要求1所述的个体化的生理状态监测分析方法,其特征在于:
如果状态变化指数大于0且小于预定值,则将该输入向量作为指标数据库的一个指标向量而纳入指标数据库。
9.一种个体化的生理状态监测分析设备,其包括:计算平台;
计算平台包括指标数据库单元、估计单元、差值单元、状态变化指数单元;
计算平台被配置,用该个体的同一时间窗的多个生理参数的数值构成向量,作为指标向量;
指标数据库单元用该个体在正常状态下的多个指标向量构成指标数据库;
估计单元用于计算与输入到估计单元的指标向量对应的估计向量;
差值单元用于计算两个指标向量之间的差值;
状态变化指数单元用于计算反映个体当前的状态偏离其正常状态的程度的状态变化指数;
计算平台将个体的当前时间窗的指标向量作为输入向量,根据该输入向量自指标数据库单元中选出与所述输入向量最为接近的预定个数的指标向量,估计单元根据该预定个数的指标向量计算出该输入向量的估计向量;该预定个数的指标向量构成基准向量集;对于基准向量集中的每个指标向量,计算平台自指标数据库单元中找到与其最为接近的预定个数的指标向量,估计单元根据该预定个数的指标向量计算与其对应的估计向量;
差值单元计算出该输入向量与其估计向量之间的差值,作为第一差值;差值单元计算出基准向量集中的每个指标向量与其估计向量的差值,作为第二差值;
状态变化指数单元根据第一差值在全部第二差值中的卡方分布累积分布情况,计算状态变化指数;其中,所述卡方分布累积分布的自由度为所述生理参数的个数,所述卡方分布累积分布的变量为根据第一差值和全部第二差值计算得到的马氏距离。
10.如权利要求9所述的个体化的生理状态监测分析设备,其特征在于:
所述设备被穿戴在该个体身上;所述多个生理参数由该个体穿戴的生理参数传感器感测得到的生理信号经处理后获得。
11.如权利要求9所述的个体化的生理状态监测分析设备,其特征在于:
所述计算平台位于远程;所述多个生理参数由该个体穿戴的生理参数传感器感测得到的生理信号经处理后获得;所述多个生理参数通过无线方式传输给所述计算平台。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010205666.6A CN112971791B (zh) | 2020-03-23 | 2020-03-23 | 一种个体化的生理状态监测分析方法和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010205666.6A CN112971791B (zh) | 2020-03-23 | 2020-03-23 | 一种个体化的生理状态监测分析方法和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112971791A CN112971791A (zh) | 2021-06-18 |
CN112971791B true CN112971791B (zh) | 2024-04-05 |
Family
ID=76344206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010205666.6A Active CN112971791B (zh) | 2020-03-23 | 2020-03-23 | 一种个体化的生理状态监测分析方法和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112971791B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005110029A2 (en) * | 2004-05-07 | 2005-11-24 | Intermed Advisor, Inc. | Method and apparatus for real time predictive modeling for chronically ill patients |
CN101365378A (zh) * | 2005-11-29 | 2009-02-11 | 风险获利有限公司 | 根据残差监测人体健康 |
CN102917661A (zh) * | 2010-01-14 | 2013-02-06 | 风险获利有限公司 | 基于多变量残差的用于人体健康监视的健康指数 |
CN105636061A (zh) * | 2016-01-14 | 2016-06-01 | 南京邮电大学 | 一种在一般性衰落信道上统一精确的能量检测方法 |
CN106510689A (zh) * | 2016-12-30 | 2017-03-22 | 东南大学 | 一种穿戴式多生理参数采集装置 |
CN108601026A (zh) * | 2018-04-02 | 2018-09-28 | 浙江大学 | 基于随机抽样一致性的感知数据错误化攻击检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2916131A1 (en) * | 2015-12-22 | 2017-06-22 | Jannatec Technologies | System for regulating body temperature of a subject |
-
2020
- 2020-03-23 CN CN202010205666.6A patent/CN112971791B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005110029A2 (en) * | 2004-05-07 | 2005-11-24 | Intermed Advisor, Inc. | Method and apparatus for real time predictive modeling for chronically ill patients |
CN101365378A (zh) * | 2005-11-29 | 2009-02-11 | 风险获利有限公司 | 根据残差监测人体健康 |
CN102917661A (zh) * | 2010-01-14 | 2013-02-06 | 风险获利有限公司 | 基于多变量残差的用于人体健康监视的健康指数 |
CN105636061A (zh) * | 2016-01-14 | 2016-06-01 | 南京邮电大学 | 一种在一般性衰落信道上统一精确的能量检测方法 |
CN106510689A (zh) * | 2016-12-30 | 2017-03-22 | 东南大学 | 一种穿戴式多生理参数采集装置 |
CN108601026A (zh) * | 2018-04-02 | 2018-09-28 | 浙江大学 | 基于随机抽样一致性的感知数据错误化攻击检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112971791A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10722179B2 (en) | Residual-based monitoring of human health | |
Radha et al. | Estimating blood pressure trends and the nocturnal dip from photoplethysmography | |
CN109009017B (zh) | 一种智能健康监测系统及其数据处理方法 | |
US20120203511A1 (en) | Method and system for real-time signal classification | |
US20160128638A1 (en) | System and method for detecting and quantifying deviations from physiological signals normality | |
Chen et al. | A new deep learning framework based on blood pressure range constraint for continuous cuffless BP estimation | |
Hung | Detection of central sleep apnea based on a single-lead ECG | |
Hung | Central sleep apnea detection using an accelerometer | |
Cimr et al. | Automatic detection of breathing disorder from ballistocardiography signals | |
CN113171080A (zh) | 一种基于可穿戴传感信息融合的能量代谢评估方法和系统 | |
CN105593860B (zh) | 用于确定复合评分分布的装置和患者健康状态确定器 | |
AU2021236462B2 (en) | Residual-based monitoring of human health | |
US20080246617A1 (en) | Monitor apparatus, system and method | |
CN112998717B (zh) | 一种量化人体生理状态的分析方法和设备 | |
CN111278353A (zh) | 一种生命体征信号噪声的检测方法与系统 | |
CN112971791B (zh) | 一种个体化的生理状态监测分析方法和设备 | |
CN112971792B (zh) | 一种基于连续生理数据的个体化状态监测分析方法和设备 | |
Mitas et al. | Wearable system for activity monitoring of the elderly | |
CN112971742A (zh) | 一种个体化的生理状态监测分析方法和设备 | |
CN112992379A (zh) | 一种个体化的生理状态监测分析方法和设备 | |
Vijayan et al. | Implementing Pattern Recognition and Matching techniques to automatically detect standardized functional tests from wearable technology | |
Nouriani et al. | System-Identification-Based Activity Recognition Algorithms With Inertial Sensors | |
Kaur et al. | Physiological and Inertial Features based Dataset for Falls and Activities: PIF v2 | |
KR102645586B1 (ko) | 생체신호를 이용한 수면 호흡 상태 분류 장치 및 방법 | |
Karako | Integration of wearable devices and deep learning: New possibilities for health management and disease prevention |
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 |