CN112891098A - 用于健康监测仪的体重测量方法 - Google Patents
用于健康监测仪的体重测量方法 Download PDFInfo
- Publication number
- CN112891098A CN112891098A CN202110072174.9A CN202110072174A CN112891098A CN 112891098 A CN112891098 A CN 112891098A CN 202110072174 A CN202110072174 A CN 202110072174A CN 112891098 A CN112891098 A CN 112891098A
- Authority
- CN
- China
- Prior art keywords
- weight
- steady state
- body weight
- state
- measurement
- 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.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61G—TRANSPORT, PERSONAL CONVEYANCES, OR ACCOMMODATION SPECIALLY ADAPTED FOR PATIENTS OR DISABLED PERSONS; OPERATING TABLES OR CHAIRS; CHAIRS FOR DENTISTRY; FUNERAL DEVICES
- A61G7/00—Beds specially adapted for nursing; Devices for lifting patients or disabled persons
- A61G7/05—Parts, details or accessories of beds
- A61G7/0527—Weighing devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
- G01G19/44—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing persons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
- G01G19/44—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing persons
- G01G19/445—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing persons in a horizontal position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
- G01G19/52—Weighing apparatus combined with other objects, e.g. furniture
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nursing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开了一种用于健康监测仪的体重测量方法,在床下设置多个压力传感器,并确定一个测量时段的时长;定义了长时间窗,并通过长时间窗判断稳态;当在床状态的稳态结束时,根据该稳态期间某一测量时刻或多个测量时刻测量的瞬间体重的值确定使用者在该稳态期间的稳态体重;当一个测量时段结束后,根据该测量时段的一个或多个稳态期间的稳态体重确定使用者在该测量时段的时段体重。本发明中,通过判断长时间窗是否为稳定状态的方法,能够计算出在床状态的稳态的时长,并得出在床状态的各稳态的稳态体重;通过评估各稳态体重的可信度,根据可信度高的稳态体重得到能够表征使用者实际体重的体重值,为基础健康监测提供准确的体重数据。
Description
技术领域
本发明涉及基础健康监测领域,特别涉及一种用于健康监测仪的体重测量方法。
背景技术
传统的电子体重秤精度一般是0.1kg,在这一精度下,读数比较容易稳定下来,一旦获得稳定,电子体重秤就可将其输出到显示屏,看起来就好像只读出了一个值。而实际上,使用任何电子称量设备称量体重,获得的读数都是一个序列,而不是单独的一个值。由于基础健康监测需要获取使用者的高精度体重数据,以便于分析监测健康状态,一般使用10g以下的高精度压力传感器,而压力传感器的测量是个动态过程,只要达到一定精度,轻微的扰动就会产生读数差异,假设每秒进行一次测量,则每秒都会产生不同的读数,这样体重数据就很难稳定下来,例如使用者在床时长为8小时,即8×60×60=28800秒,就会产生28800个称重读数;因而难以判断哪一个称重数据能够更准确的反映使用者的实际体重。
另外,人在床上的动作,比如:正在上床、正在下床、正在翻身等,都会对压力传感器读数产生影响,出现较大的瞬时的波动;使用者的一些行为习惯,比如,是否将手机、眼镜、床品、衣物之类拿上拿下等也会对读数造成扰动。如果不剔除这些影响,则会使测量的体重值产生较大的失真。
发明内容
本发明要解决的技术问题是提供了一种能够提供用于健康监测仪的体重测量方法。
本发明的技术方案如下:
一种用于健康监测仪的体重测量方法,包括以下步骤:
步骤S1、在床下设置多个压力传感器,各压力传感器每间隔第一预设时间T0测量一次压力,并将压力转换为重量,定义总读数A表示各压力传感器读数值之和,预定义空床读数B表示空床状态下的总读数;定义每一测量时刻测量的瞬间体重W=A-B;
步骤S2、确定一个测量时段的时长,并确定使用者在该测量时段内的参考体重Wr;
步骤S3、在每一次测量后判断该测量时刻是否为空床状态,如果为空床状态,则返回继续执行步骤S3;否则,执行步骤S4;
步骤S4、定义起始时刻为压力传感器的某一测量时刻,且长度为uT0的时间段为一个长时间窗,其中u为自然数,判断以当前测量时刻为结束时刻的长时间窗内的瞬间体重W是否为稳定状态,如果该长时间窗内的瞬间体重W为稳定状态,则判定当前处于稳态,并执行步骤S5;否则,返回执行步骤S3;
步骤S5、判断稳态是否结束,如果稳态结束则执行步骤S6;否则,返回继续执行步骤S5;
步骤S6、定义稳态体重Wc表示使用者在稳态期间的体重值,根据该稳态期间某一测量时刻或多个测量时刻测量的瞬间体重W的值确定使用者在该稳态期间的稳态体重Wc;返回执行步骤S3;
定义时段体重Wz表示使用者在一个测量时段测量得到的体重值,在测量时段结束后,根据该测量时段的一个或多个稳态期间的稳态体重Wc确定使用者在该测量时段的时段体重Wz。
进一步的,在所述步骤S2中,确定使用者在测量时段内的参考体重Wr的方法为:对于第一个测量时段,预先设置一个重量值作为参考体重Wr;对于第二个测量时段,参考体重Wr为第一个测量时段计算得到的使用者的时段体重Wz,对于第三个及以后的测量时段,根据该测量时段之前两个测量时段的时段体重Wz形成的趋势线确定参考体重Wr。
进一步的,在所述步骤S2中,确定使用者在测量时段内的参考体重Wr的方法为:对于第一个测量时段,预先设置一个重量值作为参考体重Wr;对于第二个至第m个测量时段,参考体重Wr为该测量时段之前所有测量时段的时段体重Wz的平均值或中位数,对于第(m+1)个及以后的测量时段,参考体重Wr为该测量时段之前m个测量时段的时段体重Wz的平均值或中位数;m为大于2的自然数。
进一步的,在所述步骤S3中,判断是否为空床状态的方法为:定义起始时刻为压力传感器的某一测量时刻,且长度为vT0的时间段为一个短时间窗,其中v为小于u的自然数,判断以当前测量时刻为结束时刻的短时间窗是否为稳定状态,如果该短时间窗为稳定状态,则将该短时间窗内各瞬间体重W的平均值或中位数做为该测量时刻的测量体重W1,比较该测量体重W1的值是否小于空床阈值,如果该测量体重W1的值小于空床阈值则判定该测量时刻为空床状态。
进一步的,判断短时间窗是否为稳定状态的方法为:定义一个短时间窗内记录的所有瞬间体重W的值的标准差为σTWD,设定一个短稳态标准差阈值δ0,当一个短时间窗结束时,若该短时间窗的σTWD≤δ0,则判定该短时间窗内的瞬间体重W处于稳定状态;若σTWD>δ0,则判定该短时间窗内的瞬间体重W处于非稳定状态。
进一步的,在所述步骤S4中,判断长时间窗内的瞬间体重W是否为稳定状态的方法为:定义一个长时间窗内记录的所有瞬间体重W的值的标准差为σTWC,设定一个长稳态标准差阈值δ1,当一个长时间窗结束时,若该长时间窗的σTWC≤δ1,则判定该长时间窗内的瞬间体重W处于稳定状态;若σTWC>δ1,则判定该长时间窗内的瞬间体重W处于非稳定状态。
进一步的,在每一稳态结束时还计算该稳态的置信度,其中稳态的置信度表示稳态中测试的稳态体重Wc的可信度,记为稳态置信度C;在测量时段结束后,确定使用者在该测量时段的时段体重Wz的方法为,将该测量时段中的各稳态按稳态置信度进行排序,根据排序情况取一个或多个稳态期间的稳态体重Wc确定使用者在该测量时段的时段体重Wz。
进一步的,计算稳态置信度C的方法为:定义参考体重置信度Cr表示根据该稳态的稳态体重Wc与参考体重Wr的差值得出的置信度,差值越小,参考体重置信度Cr越高;定义偏心置信度Cb表示根据使用者上床后总载荷的受力作用点偏离床的几何中心的偏心距得出的置信度,偏心距越小,偏心置信度Cb越高;定义稳态时长置信度Ch表示根据稳态时长得出的置信度,稳态时长越长,稳态时长置信度Ch越高;使用Cr、Cb、Ch、Cr×Ch、Cr×Cb、Cb×Ch或Cr×Cb×Ch作为稳态置信度C。
进一步的,所述参考体重置信度Cr的计算方法为:
预先指定使用者的体重标准差σr的值;
Wc<Wr时,参考体重置信度Cr表示为:
Wc≥Wr时,参考体重置信度Cr表示为:
其中,x=Wc,μ=Wr,σ=σr。
有益效果:本发明中,通过划分长时间窗并判断长时间窗是否为稳定状态的方法,能够计算出在床状态的稳态的时长,并得出在床状态的各稳态的稳态体重;通过评估各稳态体重的可信度,根据可信度高的稳态体重得到能够表征使用者实际体重的体重值,为基础健康监测提供准确的体重数据。
附图说明
图1为本发明用于健康监测仪的体重测量方法的一个优选实施例的流程图;
图2为参考体重为80Kg且使用正态分布时参考体重置信度的分别图;
图3为使用四个压力传感器时床板受力情况的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明实施例中的技术方案,并使本发明实施例的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明实施例中技术方案作进一步详细的说明。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。另外,为叙述方便,本申请中“重量”、“力”和“荷载”的单位均以质量计,单位为g。
如图1所示,本发明用于健康监测仪的体重测量方法的一个优选实施例包括以下步骤:
步骤S1、在床下设置多个压力传感器,各压力传感器每间隔第一预设时间T0测试一次压力,并将压力转换为重量,优选为在床的四个床脚下分别设置一个压力传感器,第一预设时间T0优选为1s,当然,第一预设时间T0也可设置为其他值;定义总读数A表示各压力传感器读数值之和,预定义空床读数B表示空床状态下的总读数,空床读数B的值可在测量过程中进行修改和校正;定义每一测量时刻测量的瞬间体重W=A-B。
步骤S2、确定一个测量时段的时长,并确定使用者在该测量时段内的参考体重Wr。为便于统计,一般设置一个测量时段的时长为24小时,优选为选取每一天的中午12:00:00至次日中午的11:59:59为一个测量时段。
定义时段体重Wz表示使用者在该测量时段测量得到的体重值,确定使用者在测量时段内的参考体重Wr的方法为:对于第一个测量时段,预先设置一个重量值作为参考体重Wr;对于第二个测量时段,参考体重Wr为第一个测量时段计算得到的使用者的时段体重Wz,对于第三个及以后的测量时段,根据之前两个测量时段的时段体重Wz形成的趋势线确定参考体重Wr;例如,当该测量时段之前第二个测量时段的时段体重Wz为60.1KG,该测量时段之前第一个测量时段的时段体重Wz为60KG时,则设置该测量时段的参考体重Wr=59.9KG,使三个体重值在一条直线上。当然,对于第二个及以后的测量时段,也可根据之前测量时段的平均值或中位数来确定参考体重Wr;即,对于第二个至第m个测量时段,参考体重Wr为该测量时段之前所有测量时段的时段体重Wz的平均值或中位数,对于第(m+1)个及以后的测量时段,参考体重Wr为该测量时段之前m个测量时段的时段体重Wz的平均值或中位数。m为大于2的自然数,例如,当m=5时,将当测量时段之前的5个测量时段(当测量时段之前不到5个测量时段则为当测量时段之前的所有测量时段)的时段体重Wz的平均值或中位数为该测量时段的参考体重Wr,以减小某一测量时段的测量值出现异常时对之后测量时段的影响。
步骤S3、在每一次测量后根据该测量时刻的测量体重W1的值判断该测量时刻是否为空床状态,如果为空床状态,则返回继续执行步骤S3;如果为在床状态,则执行步骤S4。
本实施方式中,判断是否为空床状态的方法为:定义起始时刻为压力传感器的某一测量时刻,且长度为vT0的时间段为一个短时间窗,其中v为小于u的自然数,判断以当前测量时刻为结束时刻的短时间窗是否为稳定状态,如果该短时间窗为稳定状态,则将该短时间窗内各瞬间体重W的平均值或中位数做为该测量时刻的测量体重W1,比较该测量体重W1的值是否小于空床阈值,如果该测量体重W1的值小于空床阈值则判定该测量时刻为空床状态;否则,判定为在床状态。由于压力传感器的测量精度在10g以下(约为7g),因此可将空床阈值设置为10g,当然,也可设置为15g、20g等其他数值。
判断短时间窗是否为稳定状态的方法为:定义一个短时间窗内记录的所有瞬间体重W的值的标准差为σTWD,设定一个短稳态标准差阈值δ0,当一个短时间窗结束时,若该短时间窗的σTWD≤δ0,则判定该短时间窗内的瞬间体重W处于稳定状态;若σTWD>δ0,则判定该短时间窗内的瞬间体重W处于非稳定状态。
在空床状态的稳态期间,每间隔第二预设时间对各压力传感器进行一次校准,第二预设时间优选为30分钟,当然,第二预设时间也可设置为其他值。本实施方式中,校准的方法为以当前测量的总读数A的值(即以当前测量时刻为结束时刻的短时间窗内各测量时刻的总读数的平均值或中位数)作为空床读数B的值,从而使该测量时刻的测量体重W1的值为0。通过校准可以减小温度变化对压力传感器的影响,提高测量精度。
步骤S4、定义起始时刻为压力传感器的某一测量时刻,且长度为uT0的时间段为一个长时间窗,其中u为自然数,优选为u=60(即uT0为1分钟),判断以当前测量时刻为结束时刻的长时间窗内的瞬间体重W是否为稳定状态,如果该长时间窗内的瞬间体重W为稳定状态,则判定当前处于稳态,并执行步骤S5;如果该长时间窗内的瞬间体重W为非稳定状态,则执行步骤S3。
判断长时间窗内的瞬间体重W是否为稳定状态的方法为:定义一个长时间窗内记录的所有瞬间体重W的值的标准差为σTWC,设定一个长稳态标准差阈值δ1,当一个长时间窗结束时,若该长时间窗的σTWC≤δ1,则判定该长时间窗内的瞬间体重W处于稳定状态;若σTWC>δ1,则判定该长时间窗内的瞬间体重W处于非稳定状态。
稳态的定义为:如果相邻的两个非稳定状态的长时间窗之间包括了连续的稳定状态的长时间窗,则将该相邻两个非稳定状态的长时间窗之间稳定状态的长时间窗的持续时期定义为一个稳态。
本实施方式中,定义第i次测量的时刻为ti,i为自然数;定义起始时刻为ti的长时间窗为TWi,定义起始时刻为ti+1的长时间窗为TWi+1,如果长时间窗TWi为非稳定状态,长时间窗TWi+1为稳定状态,则认为从长时间窗TWi+1的起始时刻ti+1开始进入了稳态,将ti+1定义为稳态的开始时刻;如果长时间窗TWi~TWi+k均为稳定状态,k为自然数,长时间窗TWi+k+1为非稳定状态,则认为在长时间窗TWi+k的结束时刻(ti+k+uT0)稳态结束,将(ti+k+uT0)定义为稳态的结束时刻。
步骤S5、判断稳态是否结束,如果稳态结束则执行步骤S6;如果稳态未结束,则返回继续执行步骤S5。
为了更全面地对使用者的健康状况进行检测,有时还会判断使用者在睡眠过程中有无发生异常安静事件。即预先设定一个异常安静标准差阈值δ2,在执行步骤S5的结果为稳态未结束时先执行以下步骤:
步骤S501、将该长时间窗的σTWC与δ2进行比较,若该长时间窗的σTWC≤δ2,则判断存在异常安静事件,执行步骤S502;否则返回执行步骤S5。异常安静事件也可作为判断使用者健康状况的依据,如果发生了异常安静事件则说明使用者的健康状况存在隐患,还可根据异常安静事件发生的频率及次数对使用者的健康状态进一步进行评估。
步骤S502、上报异常安静事件,返回执行步骤S5。
步骤S6、定义稳态体重Wc表示使用者在稳态期间的体重值,根据该稳态期间某一测量时刻或多个测量时刻测量的瞬间体重W的值确定使用者在该稳态期间的稳态体重Wc,例如,可将该稳态结束时刻的瞬间体重W的值作为该稳态期间的稳态体重Wc,也可将该稳态期间所有测量时刻的瞬间体重W的平均值作为该稳态期间的稳态体重Wc;计算该稳态的置信度,并返回执行步骤S3。
定义稳态置信度C表示稳态中测试的稳态体重Wc的可信度,计算稳态置信度C的方法为:定义参考体重置信度Cr表示根据该稳态的稳态体重Wc与参考体重Wr的差值得出的置信度,差值越小,参考体重置信度Cr越高;定义偏心置信度Cb表示根据使用者上床后总载荷的受力作用点偏离床的几何中心的偏心距得出的置信度,偏心距越小,偏心置信度Cb越高;定义稳态时长置信度Ch表示根据稳态时长得出的置信度,稳态时长越长,稳态时长置信度Ch越高;使用Cr、Cb、Ch、Cr×Ch、Cr×Cb、Cb×Ch或Cr×Cb×Ch作为稳态置信度C,优选为采用Cr×Cb×Ch作为稳态置信度C。
1、稳态时长置信度的计算方法
所述稳态时长置信度由稳态的时长确定,具体来说,先将稳态的结束时刻减去稳态的开始时刻计算出稳态的时长,例如,稳态中第一个稳定状态的长时间窗为TWi+1,最后一个稳定状态的长时间窗为TWi+K,则稳态时长为:ti+K+uT0-ti+1=(u+k-1)T0。然后人为设定一个期望的稳态时长的标准值TA,此值可经过计算不断迭代更新,因此理论上其初始值可随意指定,例如,可指定TA=1h;选取一个以TA为期望值的统计分布,这是一个基于时间间隔的统计,因此我们暂时按照一般经验选取指数分布,当然,也可选取其他分布。
对于一个随机变量X,若服从指数分布,则写作X~Exp(λ),其累积分布函数可表示为:
其中,自变量x为某一段稳态的时长,λ表示单位时间里非稳态事件发生的频率,定义TA为稳态时长的期望值,即期望的稳态时长的标准值,则λ=1/TA。
在本系统中,自变量x是某一段稳态的时长TN,λ是单位时间里非稳态事件发生的频率。使用此函数的目的是,计算小于一定时长的所有时长里出现非稳态事件的概率累积的和,即累积分布函数,作为我们的稳态时长置信度函数,这样就能给出任意一个稳态的时长TN在某一期望下的稳态时长置信度,记为Ch,相当于给出了一个置信度评分,满分为1分,分数越高,Ch越高。Ch的计算方法与累积分布函数相同,即:
其中,λ=1/TA,x=TN。
在计算前需要预先指定期望值TA的值,例如可指定此分布的期望值TA=1h(即1个小时之内只发生1次非稳态事件),则根据函数表达式,可以得出稳态时长TN和稳态时长置信度Ch的取值对应关系如表1所示:
表1稳态时长和稳态时长置信度对应表(TA=1h)
x=T<sub>N</sub> | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.0952 | 0.1813 | 0.2592 | 0.3297 | 0.3935 | 0.4512 |
x=T<sub>N</sub> | 0.7 | 0.8 | 0.9 | 1.0 | 1.1 | 1.2 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.5034 | 0.5507 | 0.5934 | 0.6321 | 0.6671 | 0.6988 |
x=T<sub>N</sub> | 1.3 | 1.4 | 1.5 | 1.6 | 1.7 | 1.8 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.7275 | 0.7534 | 0.7769 | 0.7981 | 0.8173 | 0.8347 |
x=T<sub>N</sub> | 1.9 | 2.0 | 2.1 | 2.2 | 2.3 | 2.4 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.8504 | 0.8647 | 0.8775 | 0.8892 | 0.8997 | 0.9093 |
x=T<sub>N</sub> | 2.5 | 2.6 | 2.7 | 2.8 | 2.9 | 3.0 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.9179 | 0.9257 | 0.9328 | 0.9392 | 0.9450 | 0.9502 |
查表可知,当稳态时长TN为1h时,稳态时长置信度Ch为0.6321,当稳态时长TN为0.5h时,稳态时长置信度Ch为0.3935。
稳态时长置信度Ch主要用于与参考体重置信度和偏心置信度相结合,对同一次睡眠的各个稳态下获得的稳态体重Wc与实际体重的接近程度进行比较,因此系统实际运行中,若发现设定的期望值TA过高(人在睡眠中很难获得长达1h的稳态),从而造成得分太低,尤其三种置信度结合之后,稳态置信度数C的值太小,小数点后0太多不便于观察,可适当调低期望值TA,例如,指定期望值TA=0.5h时,稳态时长TN和稳态时长置信度Ch的关系如表2所示:
表2稳态时长和稳态时长置信度对应表(TA=0.5h)
x=T<sub>N</sub> | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.1813 | 0.3297 | 0.4512 | 0.5507 | 0.6321 | 0.6988 |
x=T<sub>N</sub> | 0.7 | 0.8 | 0.9 | 1 | 1.1 | 1.2 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.7534 | 0.7981 | 0.8347 | 0.8647 | 0.8892 | 0.9093 |
x=T<sub>N</sub> | 1.3 | 1.4 | 1.5 | 1.6 | 1.7 | 1.8 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.9257 | 0.9392 | 0.9502 | 0.9592 | 0.9666 | 0.9727 |
x=T<sub>N</sub> | 1.9 | 2 | 2.1 | 2.2 | 2.3 | 2.4 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.9776 | 0.9817 | 0.985 | 0.9877 | 0.9899 | 0.9918 |
x=T<sub>N</sub> | 2.5 | 2.6 | 2.7 | 2.8 | 2.9 | 3 |
F<sub>X</sub>(x)=C<sub>h</sub> | 0.9933 | 0.9945 | 0.9955 | 0.9963 | 0.997 | 0.9975 |
如此,则稳态时长TN为0.5h时,稳态时长置信度Ch为0.6321,稳态时长TN为1h时,稳态时长置信度Ch达到0.8647。
2、参考体重置信度的计算方法
所述参考体重置信度由该稳态的稳态体重Wc与参考体重Wr的差值确定,假如某人上一次称量的体重值为Wi,在不附加任何其它干扰条件(如是否进食排便)的情况下,显然,得出的结果就可信度而言,下一次称量的体重值等于Wi时的可信度肯定优于体重值为(Wi±1kg)时,下一次称量的体重值为(Wi±1kg)时的可信度肯定也优于体重值为(Wi±2kg)时,以此类推,也就是说,下一次称量的体重值越接近Wi时越可信,否则越不可信。在实际应用中,我们以参考体重Wr作为计算的基准,来评估稳态体重Wc的可信度,即参考体重置信度,记为Cr。
本实施方式中,我们使用正态分布函数非置信区间的概率来表征参考体重置信度Cr。具体计算方法如下:
若随机变量X服从一个位置参数(均值)为μ、尺度参数(标准差)为σ的正态分布,则其概率密度可表示为:
其中,x是概率密度分布函数的自变量。定义随机事件X偏离μ的程度小于或者等于x的概率为P(X≤x),则其累积分布函数可表示为:
令μ=Wr,即以参考体重Wr作为正态分布的均值μ;标准差σ为体重标准差σr,此时,需先指定体重标准差σr的值,首次测量时可以凭经验指定σr的初始值;σr的值指定的不合适也没有关系,可在后续测量时根据历史数据的统计情况对体重标准差σr的值进行迭代更新。
对于任何一个x,考虑随机事件X偏离μ的程度大于x的概率,也就是,当x<μ时,X<x这一概率为P(X<x,和X≥2μ-x);当x≥μ时,这一概率为P(X<2μ-x,和X≥x)。括号内X的范围,就是X的非置信区间(与平时常用的置信区间相反,如图示,本系统关注的是深颜色阴影部分),而相应的P,就是X落在此非置信区间的概率,通称置信水平。
x<μ时,P(X<x,和X≥2μ-x)=2F(x);
x≥μ时,P(X<2μ-x,和X≥x)=2[1-F(x)]。
同理,对于任一稳态的稳态体重Wc,任一测量时刻的瞬间体重W偏离Wr的程度超过Wc偏离Wr的程度的概率(也就是W落在由Wc与Wr所定义的非置信区间的概率),即反映了Wc这一读数的置信水平,这就是参考体重置信度Cr的严格定义。即:
Wc<Wr时,参考体重置信度Cr表示为:
Wc≥Wr时,参考体重置信度Cr表示为:
其中,x=Wc,μ=Wr,σ=σr。
举例如下:成年人一天之内的体重波动,是他全部体重的±1%非常正常,±2%也不少见,但是超过±3%的情况会比较少。根据正态分布定义,μ±2σr之内事件出现的概率约为0.9545,我们可暂时指定±3%这个偏移水平为±2σr,那么对一个参考体重Wr=80kg的成年人来说,2σr=80×3%=2.4,即σr=1.2kg。据此计算的结果如表3所示:
表3测量体重与参考体重置信度对应表
如图2所示,深色阴影部分即为稳态体重Wc超出(Wr±2σr)范围(即Wc>82.4kg或Wc<77.6kg)时的参考体重置信度。使用正态分布计算,二者偏离程度相同,方向相反,置信度是相同。当然,由于Wc≥Wr时的置信度要高于Wc<Wr时的置信度,因此也可将正态分布调整为偏态分布来计算参考体重置信度Cr。
3、偏心置信度的计算方法
所述偏心置信度由根据使用者上床后总载荷的受力作用点偏离床的几何中心的偏心距确定,下面以在床的四个床脚各设置1个压力传感器为例进行说明;定义四个压力传感器的测量值,即绝对输出分别为A1、A2、A3、A4;定义四个压力传感器在空床状态下的测量值,即基准输出分别为B1、B2、B3、B4;将单只压力传感器在在床状态下相对于空床状态下所增加的压力值定义为该压力传感器的读数,分别记为I1、I2、I3、I4,则有I1=A1-B1,I2=A2-B2,I3=A3-B3,I4=A4-B4;床板受力情况如图3所示。
假定此床板是理想矩形,床板的重心即为其几何中心,压力传感器准确安装在其四角,则可设定如图3所示的平面直角坐标系,以I3的作用点(x3,y3)为原点,I3的作用点与I4的作用点(x4,y4)连线为x轴,I3的作用点与I1的作用点(x1,y1)连线为y轴,I2的作用点的坐标为(x2,y2),则x3=x1=0,y3=y4=0,(xc,yc)为床板的几何中心。
与参考体重置信度的计算方法相类似,按照正态分布的累积分布函数计算公式,取稳态任一时刻的总荷载F作为该稳态的总荷载F0,
优选为取稳态结束时刻的总荷载F作为F0,假设其受力作用点的坐标为(x0,y0),令概率密度分布函数的自变量x=x0,则任意测量时刻的总荷载F的受力作用点的横坐标偏离xc的程度超过x0偏离xc的程度的概率(也就是总荷载F的作用点的横坐标分量落在由x0与xc所定义的非置信区间的概率),即反映了F0由于其偏心程度而造成的置信水平,这就是偏心置信度Cb的严格定义。
根据受力平衡、分别以x轴和y轴为转轴的力矩平衡,可得联立式:
由于x2=x4,y1=y2,容易得到:
坐标点(x0,y0)偏离床板几何中心(xc,yc)的距离,可称为偏心距。但一般而言,人在床上的行为以顺着短边方向,或者说以近似平行于长边的一条线为轴的侧向翻转为主,因此一般情况下我们只需要关心顺着短边方向的偏心,可以忽略顺着长边方向的偏心,以及整个偏心距。那么,若横坐标x轴方向是床的短边,我们只需关心这一方向的偏心即可。(注:实例图中横坐标x轴方向更像长边,是由于透视原因。)
偏心置信度从直观上看,就是人越靠床中间,称重读数越可靠,越靠近床边,称重读数越不可靠。
当然,在实际应用中,压力可能通过床腿传递,床腿可能并不严格处于床的四角,四个受力点可能并非理想矩形,但这些因素造成的干扰可以忽略不计,对于一个四只床脚分布接近矩形的系统,我们只要知道床脚间的距离,就可根据上式计算。
对于一个n只床脚架在n只压力传感器上的情况,n为大于3的自然数,只要建立平面直角坐标系并定出每只压力传感器的横坐标和纵坐标,可同理进行计算,受力表达式为:
计算可得到:
其中,Ii表示第i个压力传感器在在床状态下相对于空床状态下所增加的压力值,xi表示第i个压力传感器的横坐标,yi表示第i个压力传感器的纵坐标。
只考虑横坐标x轴方向的偏心距时,可以得到偏心置信度Cb的计算方式如下:
x0<xc时,偏心置信度Cb表示为
x0≥xc时,偏心置信度Cb表示为
其中,σb为横坐标x轴方向的偏心标准差,首次测量时可以凭经验指定σb的初始值;σb的值指定的不合适也没有关系,可在后续测量时根据历史数据的统计情况对偏心标准差σb进行迭代更新。
在测量时段结束后,将该测量时段中的各稳态按稳态置信度进行排序,选取在床状态下置信度最高的稳态的稳态体重Wc或排名前几位的稳态的稳态体重Wc的平均值作为使用者的时段体重Wz,稳态体重Wc可以取该稳态在任一测量时刻测量的瞬间体重W的值,也可以取该稳态内所有测量时刻测量的瞬间体重W的平均值,优选为将在该稳态结束时刻测量的瞬间体重W的值作为稳态体重Wc。另外,还可在该测量时段内选取一次主睡眠作为目标睡眠,如果不存在主睡眠,则选取该测量时段内稳态时长最长的一次睡眠,作为目标睡眠,使用目标睡眠中稳态的稳态体重Wc进行每日参考体重Wr的计算。
主睡眠的判定方法为预先设置主睡眠时间阈值,将总睡眠时长大于或等于主睡眠时间阈值的睡眠判定为主睡眠,若使用者在睡眠中途离床,但在上下床时间阈值内重新上床,则前后算作一次睡眠,不算作两次,对于普通使用者,可以将总在床时长(即床上为在床状态的时长)作为总睡眠时长;主睡眠时间阈值和上下床时间阈值的取值可根据实际情况进行设置和调整;例如,可设置主睡眠时间阈值为3h,上下床时间阈值为30min。
可能存在如下特殊情况:
(1)横跨中午12:00:00这一时间点的睡眠(即横跨两个测量时段的睡眠),在计算整个睡眠时长时,先将该睡眠分为12:00:00之前的部分和12:00:00之后的部分,若前一部分的时长多于后一部分的时长,则将此睡眠计入前一天的统计,若前一部分的时长少于或等于后一部分的时长,则将此睡眠计入后一天的统计;因此,只要后一部分的时长超过前一部分,就可以开始进行前一天的睡眠报告统计,而不必等到这一次睡眠结束。
(2)对于长时间在床的使用者,由于其离床时间经常小于半小时,使得超过我们规定的24小时的时间里,无法分割出一次睡眠的情况,则将此使用者标记为“长期卧床”,然后基于稳态时长而不是在床时长、入睡时长的方式来标记上述的主睡眠及目标睡眠,并减小用于判断主睡眠的主睡眠时间阈值和上下床时间阈值的取值;例如,可设置主睡眠时间阈值为30min,上下床时间阈值为1min;当然,其具体取值可根据实际情况进行调整。
选定目标睡眠后,如果目标睡眠包括多段稳态,将目标睡眠中所有的稳态按置信度C从高到低排序,取置信度C最高的一段稳态的稳态体重Wc作为使用者在上一测量时段的时段体重Wz和当前测量时段的参考体重Wr;当然也可计算置信度C排序最靠前的几段稳态的稳态体重Wc的平均值,以此作为使用者在上一测量时段的时段体重Wz和当前测量时段的参考体重Wr。
在测量时段结束后,可根据该时段内的睡眠次数、睡眠时长和睡眠过程中发生的异常安静事件等指标对使用者在该测量时段的睡眠质量进行评估;并可根据之前每一测量时段的时段体重Wz的变化情况和睡眠质量情况对使用者的健康状态进行评估。
本发明未描述部分与现有技术一致,在此不做赘述。
以上仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构,直接或间接运用在其他相关的技术领域,均同理在本发明的专利保护范围之内。
Claims (9)
1.一种用于健康监测仪的体重测量方法,其特征在于,包括以下步骤:
步骤S1、在床下设置多个压力传感器,各压力传感器每间隔第一预设时间T0测量一次压力,并将压力转换为重量,定义总读数A表示各压力传感器读数值之和,预定义空床读数B表示空床状态下的总读数;定义每一测量时刻测量的瞬间体重W=A-B;
步骤S2、确定一个测量时段的时长,并确定使用者在该测量时段内的参考体重Wr;
步骤S3、在每一次测量后判断该测量时刻是否为空床状态,如果为空床状态,则返回继续执行步骤S3;否则,执行步骤S4;
步骤S4、定义起始时刻为压力传感器的某一测量时刻,且长度为uT0的时间段为一个长时间窗,其中u为自然数,判断以当前测量时刻为结束时刻的长时间窗内的瞬间体重W是否为稳定状态,如果该长时间窗内的瞬间体重W为稳定状态,则判定当前处于稳态,并执行步骤S5;否则,返回执行步骤S3;
步骤S5、判断稳态是否结束,如果稳态结束则执行步骤S6;否则,返回继续执行步骤S5;
步骤S6、定义稳态体重WC表示使用者在稳态期间的体重值,根据该稳态期间某一测量时刻或多个测量时刻测量的瞬间体重W的值确定使用者在该稳态期间的稳态体重WC;返回执行步骤S3;
定义时段体重WZ表示使用者在一个测量时段测量得到的体重值,在测量时段结束后,根据该测量时段的一个或多个稳态期间的稳态体重WC确定使用者在该测量时段的时段体重WZ。
2.根据权利要求1所述的用于健康监测仪的体重测量方法,其特征在于,在所述步骤S2中,确定使用者在测量时段内的参考体重Wr的方法为:对于第一个测量时段,预先设置一个重量值作为参考体重Wr;对于第二个测量时段,参考体重Wr为第一个测量时段计算得到的使用者的时段体重WZ,对于第三个及以后的测量时段,根据该测量时段之前两个测量时段的时段体重WZ形成的趋势线确定参考体重Wr。
3.根据权利要求1所述的用于健康监测仪的体重测量方法,其特征在于,在所述步骤S2中,确定使用者在测量时段内的参考体重Wr的方法为:对于第一个测量时段,预先设置一个重量值作为参考体重Wr;对于第二个至第m个测量时段,参考体重Wr为该测量时段之前所有测量时段的时段体重WZ的平均值或中位数,对于第(m+1)个及以后的测量时段,参考体重Wr为该测量时段之前m个测量时段的时段体重WZ的平均值或中位数;m为大于2的自然数。
4.根据权利要求1所述的用于健康监测仪的体重测量方法,其特征在于,在所述步骤S3中,判断是否为空床状态的方法为:定义起始时刻为压力传感器的某一测量时刻,且长度为vT0的时间段为一个短时间窗,其中v为小于u的自然数,判断以当前测量时刻为结束时刻的短时间窗是否为稳定状态,如果该短时间窗为稳定状态,则将该短时间窗内各瞬间体重W的平均值或中位数做为该测量时刻的测量体重W1,比较该测量体重W1的值是否小于空床阈值,如果该测量体重W1的值小于空床阈值则判定该测量时刻为空床状态。
5.根据权利要求4所述的用于健康监测仪的体重测量方法,其特征在于,判断短时间窗是否为稳定状态的方法为:定义一个短时间窗内记录的所有瞬间体重W的值的标准差为σTWD,设定一个短稳态标准差阈值δ0,当一个短时间窗结束时,若该短时间窗的σTWD≤δ0,则判定该短时间窗内的瞬间体重W处于稳定状态;若σTWD>δ0,则判定该短时间窗内的瞬间体重W处于非稳定状态。
6.根据权利要求1所述的用于健康监测仪的体重测量方法,其特征在于,在所述步骤S4中,判断长时间窗内的瞬间体重W是否为稳定状态的方法为:定义一个长时间窗内记录的所有瞬间体重W的值的标准差为σTWC,设定一个长稳态标准差阈值δ1,当一个长时间窗结束时,若该长时间窗的σTWC≤δ1,则判定该长时间窗内的瞬间体重W处于稳定状态;若σTWC>δ1,则判定该长时间窗内的瞬间体重W处于非稳定状态。
7.根据权利要求1所述的用于健康监测仪的体重测量方法,其特征在于,在每一稳态结束时还计算该稳态的置信度,其中稳态的置信度表示稳态中测试的稳态体重WC的可信度,记为稳态置信度C;在测量时段结束后,确定使用者在该测量时段的时段体重WZ的方法为,将该测量时段中的各稳态按稳态置信度进行排序,根据排序情况取一个或多个稳态期间的稳态体重WC确定使用者在该测量时段的时段体重WZ。
8.根据权利要求7所述的用于健康监测仪的体重测量方法,其特征在于,计算稳态置信度C的方法为:定义参考体重置信度Cr表示根据该稳态的稳态体重WC与参考体重Wr的差值得出的置信度,差值越小,参考体重置信度Cr越高;定义偏心置信度Cb表示根据使用者上床后总载荷的受力作用点偏离床的几何中心的偏心距得出的置信度,偏心距越小,偏心置信度Cb越高;定义稳态时长置信度Ch表示根据稳态时长得出的置信度,稳态时长越长,稳态时长置信度Ch越高;使用Cr、Cb、Ch、Cr×Ch、Cr×Cb、Cb×Ch或Cr×Cb×Ch作为稳态置信度C。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110072174.9A CN112891098B (zh) | 2021-01-19 | 2021-01-19 | 用于健康监测仪的体重测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110072174.9A CN112891098B (zh) | 2021-01-19 | 2021-01-19 | 用于健康监测仪的体重测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112891098A true CN112891098A (zh) | 2021-06-04 |
CN112891098B CN112891098B (zh) | 2022-01-25 |
Family
ID=76116182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110072174.9A Active CN112891098B (zh) | 2021-01-19 | 2021-01-19 | 用于健康监测仪的体重测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112891098B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113974567A (zh) * | 2021-11-09 | 2022-01-28 | 重庆火后草科技有限公司 | 睡眠过程的代谢率的计算方法 |
CN114027667A (zh) * | 2021-12-01 | 2022-02-11 | 慕思健康睡眠股份有限公司 | 一种在离床状态判定方法、装置、智能床垫及介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104287703A (zh) * | 2013-06-03 | 2015-01-21 | 飞比特公司 | 个人健身跟踪装置中的陀螺仪的使用 |
CN107854126A (zh) * | 2017-11-28 | 2018-03-30 | 广东乐心医疗电子股份有限公司 | 睡眠监测方法和睡眠监测设备 |
CN108366741A (zh) * | 2015-12-07 | 2018-08-03 | 通用电气公司 | 人体生理参数的确定 |
CN108836264A (zh) * | 2018-04-23 | 2018-11-20 | 中冶赛迪技术研究中心有限公司 | 生命指数计算方法及体征指数仪 |
CN110400622A (zh) * | 2019-07-17 | 2019-11-01 | 北京市海淀区疾病预防控制中心 | 一种体重健康管理系统及方法 |
WO2020023621A1 (en) * | 2018-07-24 | 2020-01-30 | Puneeth Iyengar | Improved weight loss therapy |
CN112147557A (zh) * | 2019-06-26 | 2020-12-29 | 贵州汉能移动能源研究院有限公司 | 重复性和再现性的测试方法 |
CN112204616A (zh) * | 2018-04-23 | 2021-01-08 | 皇家飞利浦有限公司 | 用于医学成像的自动化对象监测 |
-
2021
- 2021-01-19 CN CN202110072174.9A patent/CN112891098B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104287703A (zh) * | 2013-06-03 | 2015-01-21 | 飞比特公司 | 个人健身跟踪装置中的陀螺仪的使用 |
CN108366741A (zh) * | 2015-12-07 | 2018-08-03 | 通用电气公司 | 人体生理参数的确定 |
CN107854126A (zh) * | 2017-11-28 | 2018-03-30 | 广东乐心医疗电子股份有限公司 | 睡眠监测方法和睡眠监测设备 |
CN108836264A (zh) * | 2018-04-23 | 2018-11-20 | 中冶赛迪技术研究中心有限公司 | 生命指数计算方法及体征指数仪 |
CN112204616A (zh) * | 2018-04-23 | 2021-01-08 | 皇家飞利浦有限公司 | 用于医学成像的自动化对象监测 |
WO2020023621A1 (en) * | 2018-07-24 | 2020-01-30 | Puneeth Iyengar | Improved weight loss therapy |
CN112147557A (zh) * | 2019-06-26 | 2020-12-29 | 贵州汉能移动能源研究院有限公司 | 重复性和再现性的测试方法 |
CN110400622A (zh) * | 2019-07-17 | 2019-11-01 | 北京市海淀区疾病预防控制中心 | 一种体重健康管理系统及方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113974567A (zh) * | 2021-11-09 | 2022-01-28 | 重庆火后草科技有限公司 | 睡眠过程的代谢率的计算方法 |
CN113974567B (zh) * | 2021-11-09 | 2024-03-26 | 重庆火后草科技有限公司 | 睡眠过程的代谢率的计算方法 |
CN114027667A (zh) * | 2021-12-01 | 2022-02-11 | 慕思健康睡眠股份有限公司 | 一种在离床状态判定方法、装置、智能床垫及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112891098B (zh) | 2022-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112924007B (zh) | 基于目标睡眠的体重测量方法 | |
CN112891098B (zh) | 用于健康监测仪的体重测量方法 | |
CN112914884B (zh) | 通过稳态时长置信度测量睡眠状态下的体重值的方法 | |
CN112914883B (zh) | 通过偏心置信度测量睡眠状态下的体重值的方法 | |
US20060224726A1 (en) | Monitoring system | |
CN107677393B (zh) | 一种芯片动态热管理中热传感器温度实时校准的方法 | |
Adami et al. | Unobtrusive monitoring of sleep patterns | |
CN110849459B (zh) | 一种称重传感器的蠕变修正方法 | |
CN1104818A (zh) | 采用高速数字频差计数器的同步电路 | |
WO2014052232A1 (en) | Process variable measurement noise diagnostic | |
JPH10274555A (ja) | 電子秤 | |
CN117100237B (zh) | 一种居家养老用智能监测系统 | |
CN114027792B (zh) | 基于线性相关系数去干扰的睡眠过程的代谢率检测方法 | |
JPH08261854A (ja) | 天気予報の気圧読み取り方法 | |
CN112415141B (zh) | 一种甲醛测量浓度显示值的补偿方法及补偿装置 | |
WO2021189560A1 (zh) | 体温检测的矫正方法、装置、计算机设备和存储介质 | |
CN113974568B (zh) | 基于斜率去干扰的睡眠过程的代谢率的计算方法 | |
CN113974567B (zh) | 睡眠过程的代谢率的计算方法 | |
CN113552997A (zh) | 参数变化趋势确定方法、装置、设备及存储介质 | |
US6757641B1 (en) | Multi sensor transducer and weight factor | |
JP6733415B2 (ja) | 周波数カウンタおよびフィールド機器 | |
US11796381B1 (en) | Weight verification and tare process for scale | |
TWI614510B (zh) | 濕度傳感器的校正方法 | |
CN115932711B (zh) | 一种电能表计量可信度的数据处理方法、装置及介质 | |
JP6626543B2 (ja) | 呼吸障害判定装置及び呼吸障害判定方法 |
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 |