一种心率变异性的C*模式度指标检测分析方法及仪器技术

本发明专利技术为一种心率变异性的Ｃ↓［０］模式度指标检测分析方法及仪器。本发明专利技术采用Ｃ↓［０］模式度ＰＭ作为心率变异性的非线性分析指标，其检测分析方法包括利用动态心电仪器记录２４小时心电信号，在人工辅助下自动识别动态心电ＲＲ间期，构成时间序列ＲＲｉ，运用记录的时间序列，给出了非线性Ｃ↓［０］模式度ＰＭ的计算方法。本发明专利技术提出的检测指标能够较好地反应心率变化的动力学行为特征和生理、病理特征，该指标区分能力强，显著性高，具有广泛的临床应用价值。（*该技术在2024年保护过期，可自由使用*）

【技术实现步骤摘要】

本专利技术是一种动态心电图中心率变异性的检测分析方法以及测量仪器。
技术介绍
心率变异的临床指标是一种定量的、可重复应用的指标，它既是对窦性心律不齐的程度的定量指标，也是评价自主神经调节功能的重要定量指标。正确测定心率变异的种种指标，无论是对生理、病理研究，还是对临床的诊断和防治，都具有指导性的意义。特别对分析和观察交感神经和迷走神经的调控功能、对界定心肌梗死病人预后的危险因子、对糖尿病、阵发性心律失常、高血压病、睡眠呼吸暂停综合症等的诊断具有重要的临床应用价值。目前最为常用或者常见的方法是(1)线性的时域指标，主要包括检测动态心电全程心电RR间期的标准差SDNN、短程每5分钟间期的SDANN/SDANN的均值，相邻间期的参数SDSD/NN50/PNN90、微分/对数指数等；(2)线性的频域指标，主要包括短程(每5分钟)的总功率、极低频、低频、高频、低频高频正态以及低频高频比值；长程(24小时)的超低频、极低频、低频、高频以及对数坐标中的频谱的线性内插值斜率；(3)非线性分析方法，目前运用最多的是RR间期、RR间期差值的庞加莱(Poincaré)散点图，主要依靠肉眼直观图形特征来进行区分。对于现行时域、频域线性方法的进一步临床运用，其主要症结是难以准确刻画和分析心率变化——这一受到血流动力学，电生理及激素的变化及自主神经和中枢系统等多项因素影响的非常复杂的过程。所以现行的线性方法过于粗糙，对外界噪声干扰比较敏感，其数学模型并不符合心率调控系统本身是非线性系统的这一基本事实。对于现行的非线性散点图分析方法，主要还是凭借肉眼的观察和一些简单数值定量指标的分析，而这是远远不够的，也不十分符合心率变化这一高维非线性系统的基本事实。所以无论是现行的那一类方法，都有其显著的局限性，对心率变异性指标的临床运用和推广带来了阻碍。
技术实现思路
本专利技术的目的在于提出一种能够反映心率变异过程和非线性特征的心率变异性非线性指标及其检测方法，并提供一种相应于该方法的心率变异性非线性指标的检测仪器。本专利技术中，心率变异性的非线性分析指标称作为C0模式度PM(PatternMeasurement)。C0模式度PM的检测分析方法如下(1)利用动态心电仪器记录24小时心电信号ECG；(2)将采集的心电信号ECG在计算机中回放，其中考虑滤波、数字信号压缩等环节，从而显示24小时动态心电波形；(3)基于人工辅助下的自动识别动态心电RR间期，构成时间序列RRi＝{ri1，ri2，…，rij，…rin}，其中i＝1，2，…，24；(4)基于人工辅助下的自动识别动态心电R峰的跃度(以水平基线为准)，构成时间序列RLi＝{li1，li2，…，lij，…lin}，其中i＝1，2，…，24；(5)通过上述记录的时间序列，用计算机分析测量相应的心率变异性的非线性C0模式度PM的值。上述这种非线性C0模式度PM指标完全是基于建立合适的理论模型以及严格的计算方法。下面我们给出这种指标的相关模型和计算方法。考虑一小时的RR间期时间序列RRi＝{ri1，ri2，…，rij，…rin}(i＝1，2，…，24)。为了给出能够反映自主神经调节功能强弱正常与否的一个非线性心率变异指标，则需要引入离散信号傅立叶变换和其逆变换。为了记号方便，我们重新记RR间期时间序列rik＝ri(k)，i＝1，2，…，24，k＝1，2，…，n。于是，F(j)=1nΣk=1nri(k)·e-2π-1kjn,j=1,2,···,n,--(2.1)]]>为时间序列rik＝ri(k)的傅立叶变换，生成新的时间序列DFT＝{F(j)，j＝1，2，…，n}。另一方面，新生成时间序列DFT的逆变换为，ri(k)=F~(k)=Σj=1nF(j)·e2π-1kjn,k=1,2,···,n,--(2.2)]]>下面给出心率变异非线性指标C0模式度的计算步骤(1)将每一个小时的RR间期心电时间序列RRi＝{ri1，ri1，…，rij，…rin}＝{ri(1)，ri(2)，…，ri(j)，…ri(n)}根据公式(2.1)作傅立叶变换得到新的序列DFTi＝{Fi(1)，Fi(2)，…，Fi(j)，…Fi(n)}；(2)记Gin=1n(Σj=1n|Fi(j)|2),]]>如下构造新的序列 于是，得到新的序列集合 DF^Ti={F^i(1),F^i(2),···,F^i(j),···F^i(n)};]]>(3)根据公式(2.2)，对于序列集合 作傅立叶逆变换，从而得到r.i={r.i(1),r.i(2),···,r.i(j),···,r.i(n)};]]>(4)于是，我们就可以如下定义第i小时的RR间期序列模式度PMi的运算公式PMi=Σk=1n|ri(k)-ri.(k)|2Σk=1n|ri(k)|2.--(2.4)]]>对于如上(2.4)定义的第i小时的RR间期序列模式度模式度(见图1)，具有如下若干性质(a)如果心跳RR间期为常值，即ri(k)＝const；那么PMi＝0；(b)如果心跳RR间期为周期的，即存在某个自然数к，使得ri(k)＝ri(k+к)成立，那么PMi趋向于0；(c)如果心跳间期RRi＝{ri1，ri2，…，rij，…rin}是一个独立的随机时间序列，均服从相同的概率分布，并且存在有限的四阶矩，μ是均值，σ是方差，那么PMi依概率1趋向于常数σ2μ2+σ2,]]>特别地，当μ＝0时，PMi依概率1趋向于1。以上即给出了定量评价自主神经功能的心率变异性非线性指标——C0模式度——的检测方法和模型构造的基本思路。当然对于给出R峰跃度时间序列的相应非线性指标，可以用以上类似方法来说明医学临床生理机制。根据上述心率变异非线性C0模式度指标的检测方法，本专利技术特别设计了相应的检测装置，该装置包括动态心电记录器、心电信号预处理系统和非线性C0模式度指标检测系统。动态心电记录器分为两种，一种为磁带式记录器，一种为闪光卡式记录器。均具备24-30小时动态心电采集存储功能。在心电信号预处理系统中，对于磁带式记录器提供的数据采用以单片机为核心，配有读写存储器(RAM)、只读存储器(ROM)、A/D转换器以及计算机连接并行通讯接口；对于闪光卡式记录器提供的数据可以直接通过Flash卡USB读取器及配以相应的驱动程序、解码软件即可。在预处理本文档来自技高网...

【技术保护点】
一种心率变异性的Ｃ↓［０］模式度ＰＭ指标检测分析方法，其特征在于步骤如下（１）利用动态心电仪器记录２４小时心电信号ＥＣＧ；（２）将采集的心电信号ＥＣＧ在计算机中回放，显示２４小时动态心电波形；（３）基于人工辅助下的自动识别动态心电ＲＲ间期，构成时间序列ＲＲ↓［ｉ］＝｛ｒｉ↓［１］，ｒｉ↓［２］，…，ｒｉ↓［ｊ］，…ｒｉ↓［ｎ］｝，其中ｉ＝１，２，…，２４；（４）基于人工辅助下的自动识别动态心电Ｒ峰的跃度，构成时间序列ＲＬ↓［ｉ］＝｛ｌｉ↓［１］，ｌｉ↓［２］，…，ｌｉ↓［ｊ］，…ｌｉ↓［ｎ］｝，其中ｉ＝１，２，…，２４；（５）通过上述记录的时间序列，用计算机分析测量相应的心率变异性的非线性Ｃ↓［０］模式度ＰＭ的值。

【技术特征摘要】

【专利技术属性】
技术研发人员：阮炯，林伟，顾凡及，蔡志杰，
申请(专利权)人：阮炯，林伟，顾凡及，蔡志杰，
类型：发明
国别省市：31[中国|上海]