判定对象昼夜节律特征的方法.pdf

1、(10)申请公布号 CN 103168306 A(43)申请公布日 2013.06.19CN103168306A*CN103168306A*(21)申请号 201180052087.1(22)申请日 2011.10.0510189084.6 2010.10.27 EPG06F 19/00(2006.01)(71)申请人皇家飞利浦电子股份有限公司地址荷兰艾恩德霍芬(72)发明人 A.M.巴罗索 B.萨劳克L.J.M.施兰根(74)专利代理机构中国专利代理(香港)有限公司 72001代理人董宁 汪扬(54) 发明名称判定对象昼夜节律特征的方法(57) 摘要本发明涉及确定对象昼夜节律特征的方法，其中

2、包括测量指示对象的心功能的第一输入信号，测量至少一个指示对象活动的第二输入信号，将第一输入信号和第二输入信号组合为表示对象昼夜节律的周期性输出信号，并且确定该周期性输出信号的至少一个特征。(30)优先权数据(85)PCT申请进入国家阶段日2013.04.26(86)PCT申请的申请数据PCT/IB2011/054379 2011.10.05(87)PCT申请的公布数据WO2012/056342 EN 2012.05.03(51)Int.Cl.权利要求书1页 说明书4页 附图1页(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书1页 说明书4页 附图1页(10)申请公布号 CN

3、103168306 ACN 103168306 A1/1页21. 一种用于确定对象（10）昼夜节律特征的方法，包括以下步骤：- 测量指示对象（10）在时段T内心功能的第一输入信号xn；- 测量至少一个指示对象（10）在与所述时段T重叠的时段内的活动的第二输入信号xn；- 将所述第一输入信号xn和所述第二输入信号xn组合为表示所述对象的所述昼夜节律的周期性输出信号yn，以及- 确定该周期性输出信号yn的至少一个特征。2. 根据权利要求1的方法，其中所述特征为所述周期性输出信号yn的峰值期。3. 根据权利要求1或2的方法，其中所述特征为所述周期性输出信号yn的振幅。4. 根据权利要求1至3的方法，

4、其中所述第一输入信号和所述第二输入信号分别通过采样获取，从而产生大量第一输入信号值xn和第二输入信号值xn，以及根据采样的第一输入信号值xn和第二输入信号值xn计算大量输出信号值yn。5. 根据权利要求4的方法，其中所述周期性输出信号通过回归分析从计算出的大量输出信号值yn导出。6. 根据权利要求4或5的方法，其中所述大量输出信号值yn根据以所述第一输入信号值xn和第二输入信号值xn作为外生输入的自动回归滑动平均模型（ARMAX），从所述第一输入信号值xn和第二输入信号值xn计算得出。7. 根据权利要求4至6的方法，其中所述第一输入信号值xn和所述第二输入信号值xn通过在高于1 Hz的频率上采

5、样获取。8. 根据权利要求4至7的方法，其中所述第一输入信号值xn或所述第二输入信号值xn通过在第一频率上对最初信号值进行采样获取，所述最初信号值通过在高于所述第一频率的第二频率上采样获取。9. 根据上述权利要求之一的方法，其中所述第一输入信号为所述对象（10）的心跳间隔（IBI）。10. 根据权利要求9结合权利要求8的方法，其中原始ECG信号值通过在第二频率512 Hz上对原始ECG信号进行采样获取，IBI信号通过所获取的ECG信号值导出，以及表示所述第一输入信号值的IBI信号值通过在第一频率60 Hz上对所导出的IBI信号进行采样获取。11. 根据上述权利要求之一的方法，其中所述第二输入信

6、号为所述对象（10）的腕动计信号。12. 根据权利要求1至10之一的方法，其中所述第二输入信号为指示所述对象（10）的光照量的信号。权 利 要 求 书CN 103168306 A1/4页3判定对象昼夜节律特征的方法技术领域0001 本发明涉及判定对象昼夜节律（circadian rhythm）特征的方法。背景技术0002 众所周知，心血管和代谢功能存在昼夜变化。但是，只是在最近才有显示，心血管和代谢过程不仅受醒睡周期行为影响，而且还部分地受主生理节奏（circadian）起搏器控制。主生理节奏起搏器与心输出量之间的关系表明，可使用后者可靠地导出前者的振荡参数，例如振幅和峰值期（即，振荡峰值时间

7、）。这些参数在提升人类睡眠和身体或精神表现上具有重要应用。0003 了解个人的生理节奏相位（phase）有助于优化可改变昼夜节律的光照干预时间。这种昼夜节律改变在治疗时差综合症，帮助轮班工人适应夜班工作安排，以及引导具有睡眠障碍（例如，睡眠相位后移症候群或睡眠相位前移症侯群）的人进入更正常的醒睡模式时间的方面很有用。0004 优化个人的生理节奏相位可调整个人的生物节律，以便在期望的时间上出现最大警觉和性能。日常生活、工作绩效、运动成绩等均可预期从此中获益。0005 使用心输出量作为主生理节奏起搏器动态变化估计手段的主要优点是可使用相对不明显的传感器可靠地测量变化。这些传感器产生的数据可用于收集

8、包含细微内生昼夜摆动（circadian variation）的信号。假设这些摆动的振幅很小，那么它们通常会被竞争过程的影响所掩蔽，例如，身体或精神压力对心脏活动的影响。很多这样的掩蔽效应会产生易于过滤的失调，因为它们引入的频谱分量不在正常生理节奏范围之内。但是，其他掩蔽效应本质上也具有生理节奏性，并且需要更精细的过滤法。最近研究表明，心脏活动的生理节奏调节受先前觉醒时期的影响。醒睡周期遵循生理节奏周期，该生理节奏周期影响类似频率范围内心脏活动的方式不同于主生理节奏起搏器。发明内容0006 因此，本发明的目标是提供一种确定对象昼夜节律特征的方法，其允许通过分析对象的心功能，消除所述心功能中存在

9、的内生昼夜节律的掩蔽效应来估计主起搏器振荡参数。0007 此目标通过包括权利要求1的特征的方法实现。0008 根据此方法，测量有关特定时段内对象的心功能的第一输入信号。此外，测量指示所述对象在与测量所述第一输入信号的所述第一时段重叠的时段内的生理活动的第二输入信号。0009 所述第一输入信号和所述第二输入信号通过组合产生表示所述对象昼夜节律的周期性输出信号。根据此输出信号，导出所需的所述昼夜节律特征。0010 组合所述第一输入信号和所述第二输入信号背后的理念是：后者包括有关掩蔽效应的信息，所述效应掩蔽有关所述第一信号中存在的内生昼夜节律的信息。所述组合允许说 明 书CN 103168306 A

10、2/4页4消除所述第一输入信号（表示所述心功能）中的无用特征，以便描述所述昼夜节律动态变化的振荡参数可以直接根据所述周期性输出信号导出。0011 所述第一输入信号和所述第二输入信号的所述组合可通过自回归、消色、独立分量（component）分析或（非线性）主分量分析执行，从而产生由所述周期性输出信号实现的增强昼夜节律表示。该增强信号然后可用于导出所述振荡参数的估计，例如通过应用谐波回归。0012 尽管所述第一输入信号可以是任何心功能标志，但是第二信号可通过来自对象的生理、行为或环境数据表示，这些数据提供有关她/他在与测量所述心功能时段重叠的时段内的先前不眠期的指示，例如，腕动计、睡眠锁（sle

11、ep locks）、问卷回答、皮肤温度、光照量等。0013 优选地，待确定的所述昼夜节律的特征为所述周期性输出信号的峰值期。0014 依据根据本发明的方法的另一优选实施例，所述特征为所述周期性输出信号的振幅。0015 在另一优选实施例中，所述第一输入信号和所述第二输入信号通过采样获取，从而分别产生大量第一输入信号值和第二输入信号值，并且根据采样的第一输入信号值和第二输入信号值计算大量输出信号值。0016 在该实施例中，所述第一和第二输入信号为通过对所述对象的心功能进行采样以及监视所述对象的活动分别获取的离散输入信号值。这样产生第一输入信号值的第一集合和第二输入信号值的第二集合。根据这些集合，例

12、如，可通过对以所述第一和第二输入信号值作为输入变量的方程组求解来获得作为解的所述输出信号值，来导出所述输出信号值的集合。0017 优选地，所述周期性输出信号通过回归分析从计算出的大量输出信号值导出。0018 根据本发明的另一优选实施例，根据以所述第一输入信号值和第二输入信号值作为外生输入的自动回归滑动平均模型（ARMAX），从所述第一输入信号值和第二输入信号值计算出大量输出信号值。0019 自动回归滑动平均模型在信号处理中公知地用于处理自相关时间序列数据。此模型由两部分组成：即自动回归部分和滑动平均部分。其数学基础是将两个时间序列设为相互间关系的线性方程组，本文用所述第一输入信号值和第二输入信

13、号值表示时间序列。对此线性方程求解得出一系列表示所述对象的所述昼夜节律的输出信号值。0020 优选地，所述第一输入信号值和所述第二输入信号值通过在高于1 Hz的频率上采样获取。0021 根据另一优选实施例，所述第一输入信号值或所述第二输入信号值通过在第一频率上对原始信号值进行采样获取，所述原始信号值通过在高于所述第一频率的第二频率上采样获取。0022 在实践中，这意味着由所述第一输入信号值或所述第二输入值表示的时间序列的其中之一通过对最初（原始）信号值进行采样获取，反过来，所述原始信号值通过在相对较高频率上对测量参数进行采样获取。仅作为一个实例，所述原始信号值是通过在频率512 Hz上对心脏参

14、数进行采样获取的原始值。这些原始信号值不能直接与通过在60 Hz上采样获取的第二输入信号值的第二时间序列一起计算。因此，在相对较低的频率60 Hz上对所说 明 书CN 103168306 A3/4页5述原始信号值进行采样以获取所述第一输入信号值，所述第一输入信号值然后可与第二输入信号值进行组合，例如，通过上述ARMAX模型。0023 优选地，所述第一输入信号为所述对象的心跳间隔（IBI）。0024 根据另一优选实施例，原始ECG信号值通过在第二频率512 Hz上对原始ECG信号进行采样获取，IBI信号通过所获取的ECG信号值导出，并且表示所述第一输入信号值的IBI信号值通过在第一频率60 Hz

15、上对所导出的IBI信号进行采样获取。0025 在本发明的另一优选实施例中，所述第二输入信号为所述对象的腕动计信号。0026 根据又一优选实施例，所述第二输入信号为指示所述对象的光照量的信号。附图说明0027 本发明的上述和其他方面将参考下面描述的实施例进行阐述并变得显而易见。0028 在附图中，唯一的一个图是示意性流程图，其显示一个根据本发明确定对象昼夜节律特征的方法的实施例。具体实施方式0029 在附图中的流程图所示的方法中，收集对象10的生理信号。通过采取扣带或类似形式的可佩戴ECG设备中集成的适当传感器获取ECG信号。原始ECG信号值通过在频率512 Hz上采样获取（步骤12）。此采样操

16、作在预定时段T期间执行来产生原始ECG信号值集合。第二，通过另一适当的传感器收集对象10的其他生理数据。这些第二数据指示对象10在与上述原始ECG信号采样时段T重叠的时段内的生理活动。这种对活动相关信号进行采样的操作在较低频率60 Hz上执行（步骤14）。可在袖口或类似位置内集成适当的传感器来测量对象10的手腕活动信号。需要指出，这些指示对象10活动的第二生理数据也可以是基于睡眠日志、皮肤温度、光照量等的其他数据。假设这些数据包含有关对象10的先前不眠期的信息。在下文中，将显示，这些数据可与有关对象10心功能的数据组合以消除心功能信号的掩蔽效应，从而产生表示对象10生理节奏活动的信号函数。00

17、30 根据原始ECG信号，通过计算R-R间隔导出IBI信号。最终的IBI信号通过步骤16在频率60 Hz上进一步被采样。步骤16的结果是时间间隔T内有关对象10 IBI的输入信号值的集合。另一方面，步骤14的结果是具有相同数量第二输入信号值的集合，这些第二输入信号值指示对象10在与时段T重叠的时段内的活动。第一和第二输入信号去趋于消除零频率分量。0031 在下面的步骤18，使用自动回归滑动平均模型（ARMAX）对第一输入信号值和第二输入信号值进行数学组合，在此模型中，第一输入信号值和第二输入信号值表示外生输入。令xn为第二输入信号值，xn为第一输入信号值，则可在ARMAX模型内计算一系列表示x

18、n和xn的线性函数的yn输出信号值。此函数yn是增强昼夜节律表示，包含有关昼夜节律峰值期的信息。对于任何n，均可通过过去值的线性组合，针对xn和xn计算yn。0032 更具体地说明，第一和第二输入信号值xn和xn可通过以下线性模型进行组合以用于获取增强昼夜节律表示yn：（1）说 明 书CN 103168306 A4/4页6当n为非负整数时，q为后移算子，即，qyn=yn-1。Wn为高斯白噪声过程，其中均值为0，标准差为，并且A(q)、B(q)和C(q)为多项式，定义如下：（2）（3）（4）在这些表达式中，值r、s1 s2T、u ai、bij、ci和是模型的给定常数。然后增强节律表示yn服从于余

19、弦法回归，形式如下：(5)系数的回归是所需的昼夜节律估计量。0033 表示峰值期，即，可根据上述方程(5)的结构导出的昼夜节律函数的振荡峰值时间（步骤22）。一旦导出系数，便可使用它控制可改变对象10昼夜节律的光照干预时间并帮助优化生理节奏活动时间。例如，越接近中心体温度最低的时间（休息日自发觉醒之前的大约两小时），光照干预的效果就越好效并且所用的治疗时间也越短。0034 需要指出，表示昼夜节律函数的峰值期的系数只是由周期性输出信号xn表示的一个昼夜节律特征实例。例如，该信号的振幅可以是另一可用于控制光照干预，从而优化对象10生理节奏相位的特征。一般而言，可将主体10的心跳间隔（IBI）以外的

20、另一心功能用作第一输入信号，而使用不用于腕动计信号的另一第二输入信号表示对象10的活动，如上所述。0035 尽管在附图和上面的描述中详细地示出和介绍了本发明，但此类图示和说明应被视为是说明性和示例性的，而非限制性的；本发明不限于所公开的实施例。所公开实施例的其他变形可被实施本发明的所属领域技术人员通过对附图，说明书和所附权利要求书的研究而理解和实现。在权利要求中，单词“包括”不排除其他元素或步骤，不定冠词“一”或“一个”不排除多个。某些特征被记载在相互不同的从属权利要求中并不表示这些特征的组合不能被有利地使用。在权利要求中的任何参考标号不应被解释为对范围的限制。说 明 书CN 103168306 A1/1页7图 1说 明 书 附 图CN 103168306 A