基于形态学的T波交替散点图法的心电信号分析方法.pdf

1、(10)授权公告号 (45)授权公告日 (21)申请号 201310011202.1 (22)申请日 2013.01.11 A61B 5/0452(2006.01) (73)专利权人 山东师范大学 地址 250014 山东省济南市历下区文化东路 88 号 (72)发明人 赵捷 郭慧 李菲 李田田 (74)专利代理机构 济南圣达知识产权代理有限 公司 37221 代理人 王吉勇 CN 101176660 A,2008.05.14, CN 1149440 A,1997.05.14, US 7813791 B1,2010.10.12, US 2006167364 A1,2006.07.27, WO

2、0161550 A2,2001.08.23, 徐舫舟 .T 波交替散点图法 ： 利用非线性动 力学系统表示周期性的散点图 .中国组织工程 研究与临床康复 .2010, 第 14 卷 ( 第 9 期 ), 第 1645-1648 页 . (54) 发明名称 基于形态学的 T 波交替散点图法的心电信号 分析方法 (57) 摘要 本发明公开了一种基于形态学的 T 波交替散 点图法的心电信号分析方法。针对谱分析法检测 T波交替的计算复杂、 检测每个心拍T波交替幅值 的困难等缺点， 提出一种基于形态学的 T 波交替 散点图法的研究， 为国内相关临床研究开展提供 一种新的 T 波交替检测有效量化指标。本发

3、明通 过T波窗口分析法， 对选自美国MIT/BIH标准心律 失常数据库和欧洲ST-T心电数据库的128个连续 心拍的 T 波进行采样， 利用所得采样点作一次差 分散点图， 用散点图的矢量角度指数对 T 波交替 进行定量分析， 得出判断 T 波交替存在与否的矢 量角度指数取值范围。将散点图法的矢量角度指 数和谱分析法的 T 波交替幅值对 T 波交替的检测 结果进行比较， 判断矢量角度指数检测 T 波交替 的可靠性。 (51)Int.Cl. (56)对比文件 审查员 高瑞玲 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要求书2页 说明书8页 附图2页 CN 103006207 B

4、2016.03.02 CN 103006207 B 1/2 页 2 1.一种基于形态学的 T 波交替散点图法的心电信号分析方法， 其特征在于， 包括以下 步骤 ： (1) 选取心电信号并进行心电信号预处理 ； 步骤 (1) 中所述心电信号预处理是指去除基线漂移、 工频干扰及肌电干扰的随机干 扰， 获得具有明显特征的心电信号 ； (2) 获取采样点 ： 采用 T 波窗口分析法， 对 128 个心拍进行采样， 每一个心拍采用 T 波 窗口法选定 7 个采样点， 就得到检测数据共有 128*7 个采样点 ； 采用 T 波窗口分析法， 以便 于更准确地检测微伏 TWA， 由 RR 间期的大小来选择 T

5、 波窗口起点及窗宽， T 波起点按经验选 取， 窗宽取根据 QT 和 RR 间期的关系式， 计算出窗口的终点， 如公式所示 ： 其中， Q 表示 Q 波起点， 近似在 R 峰前 50ms 处 ； T 表示 T 波起点 ； RR 表示间期， k 是常 数 ； (3) 绘制散点图 ： 将 128 个相邻心拍的采样点作一次差分， 绘制一次差分的 T 波散点 图 ； (4) 确定矢量角度指数取值范围， 以判断 T 波交替存在与否 ； (5) 绘制功率谱曲线图 ： 利用快速傅立叶变换将采样点的幅度转化成能量谱并进行叠 加平均得到其功率谱曲线图 ； (6) 计算 T 波交替幅度值 ： 对心电信号中 N 个

6、连续 ST-T 心动周期的一组采样点进行 FFT 快速傅里叶变换分析， 其 功率谱的计算公式如下 ： 其中， 0 l 127， N 128， m 7， bi为第 i 个样点的值。 设定0.460.49周期/心拍内的噪声视为背景噪声， 其平均频率值用Sn表示 ； 0.5周 期 / 心拍处的频率值用 S0.5表示， T 波交替的功率值记作 S TWA， STWA的计算公式如下 ： STWA S 0.5-Sn 规定 ： 当STWA0， T 波交替幅度值为 判定 T 波交替存在 ； 通过观察 0.5cycles/beat 处是否存在一明显的谱线， 即为 T 波交替的标识谱线， 可判 断在心动周期内是否

7、存在 T 波交替 ： 若线谱不明显， 则心跳周期中不存在 T 波交替 ； 若有明 显线谱， 则存在 T 波交替 ； 步骤 (4) 中所述判断 T 波交替存在与否的散点图矢量角度指数取值范围的确定， 包含 以下步骤 ： 41) 根据散点聚集区域的不同形态以及散点图的各种定量参数值， 可得到所需的不同 权 利 要 求 书 CN 103006207 B 2 2/2 页 3 信息 ； 通过计算散点图的矢量角度指数 VAI 来判断 T 波交替的存在与否， 其计算公式为 ： 式中， N 128， xi为第 i 点的横坐标， i为第 i 点到坐标原点的连线与横坐标之间的 夹角度， 即 i arctan(x

8、i+2-xi+1)/(xi+1-xi) ； 用软件进行仿真时， 将 45 度化成弧度制取为 0.7854， VAI 的大小表示相邻 T 波间的幅 值差沿 45线的散开程度 ； 42) 规定 ： 当矢量角度指数 0.9rad VAI 1.03rad 时， 则在心电信号中， 含有 T 波交 替 ； 当 VAI 0.9rad 或 VAI 1.03rad 时， 则在心电信号中， 不含 T 波交替。 2.根据权利要求 1 所述的分析方法， 其特征在于， 所述步骤 (3)、 (4)、 (6) 均通过 MATLAB7.0 软件实现。 3.根据权利要求1所述的分析方法， 其特征在于， 所述步骤(1)中的心电信

9、号是从美国 MIT/BIH 标准心律失常数据库和欧洲 ST-T 心电数据库中获得。 权 利 要 求 书 CN 103006207 B 3 1/8 页 4 基于形态学的 T 波交替散点图法的心电信号分析方法 技术领域 0001 本发明涉及一种心电信号分析方法， 具体涉及一种基于形态学的 T 波交替散点图 法的心电信号分析方法。属于生物医学信号处理技术领域。 背景技术 0002 T 波交替 (T Wave Alternans， TWA) 是在规整的心律时， 一种体现在体表心电图上 的 T 波形态、 极性和振幅逐拍交替变化的心电变异现象。虽然 T 波交替的产生机理与室性 心律失常的关系尚未有一致的定

10、论， 但大量的临床和动物实验研究表明， T 波交替与室性心 律失常、 室速、 室颤、 心源性猝死等有直接的联系， 微伏级 T 波交替已成为预测室性心律失 常的重要指标。 0003 T 波交替是一种非稳态的变异性现象。按照 T 波交替检测统计方法的不同， T 波交 替检测被分为短时傅里叶变换法、 符号变换法和非线性法。频谱分析法 (SM) 是短时傅里叶 变换法中最成熟一种方法， 散点图法是非线性法的一种。随着数字信号处理技术和计算机 技术的快速发展， 微伏级的 T 波交替被检出的精确度越来越高。使用数字信号处理技术检 测体表电路图中是否存在微伏级的 T 波交替现象， 已作为预测室性心律失常、 心

11、源性猝死 风险的独立指标， 并成为一项新的无创心电检测方法。 0004 文献 利用 Poincare 散点图法检测 t 波交替的指标研究 ( 李斐、 赵捷等， 现代生 物医学进展， VOL.12， NO.14MAY2012) 和 基于 Poincare 散点图与谱分析法的 T 波交替检 测 ( 李斐， 山东师范大学， 2012 年 ) 主要是根据散点图形态和 T 波交替幅值变化的特点， 利 用散点到 x+y 0 直线的距离均值作为定量检测指标来判断 T 波交替的存在与否。由于该 方法受 T 波幅值变化影响较大， 对于幅值差别较大的心电， 会出现误检， 降低了检测结果的 准确性。 0005 文献

12、 心率的 Poincare 散点图量化指标 ( 生物医学工程学杂志， 2000 ： 17(4)， P433-436) 只提出了基于散点图的四个参数 ： 长轴、 短轴、 面积和角度， 并没有将矢量角度 指数应用于实际心电进行检测， 也没验证其可靠性。 0006 中国专利申请201210064226公开了一种心电图中的T波交替检测的方法， 采用相 关分析法检测 T 波交替的存在与否， 其缺点是相关分析法的频率分辨性较差， 如当呼吸调 制的振幅较大时， 对相关分析法结果的影响要比谱分析法大得多。 0007 中国专利申请 201110418860.3 公开了一种基于模型的动态心电图 T 波交替定量 分

13、析方法， 采用 T 波的解析函数即维格纳分布 (Winger-Ville Distribution， WVD) 对心电 信号进行定量分析， 而 WVD 容易产生交叉干扰项， 加之其运算量较大， 已超过了现有硬件的 处理能力。 0008 中国专利申请 200410003212.1 公开了一种用于改善 T 波交替的测量的方法和系 统， 利用T波交替是每二次搏动自动重复并已经连接到基本心脏不稳定性的心电图的ST波 段或T波的规则的或搏动对搏动的变化的特点， 将某些搏动从T波交替测量中排除出去， 该 方法只是提取出T波交替信号， 而并没有对T波交替信号做进一步的检测， 缺乏对搏动变率 说 明 书 CN

14、 103006207 B 4 2/8 页 5 可靠性的检验。 发明内容 0009 本发明的目的是为克服上述现有技术的不足， 提供一种基于形态学的 T 波交替散 点图法的心电信号分析方法。针对谱分析法检测 T 波交替的计算复杂、 检测每个心拍 T 波 交替幅值的困难等缺点， 提出一种基于形态学的 T 波交替散点图法的研究， 为国内相关临 床研究开展提供一种新的 T 波交替检测有效量化指标。本发明的技术方案是通过 T 波窗口 分析法， 对选自美国 MIT/BIH 标准心律失常数据库和欧洲 ST-T 心电数据库的 128 个连续心 拍的T波进行采样， 利用所得采样点作一次差分散点图， 用散点图的矢量

15、角度指数对T波交 替进行定量分析， 得出判断 T 波交替存在与否的矢量角度指数取值范围。将散点图法的矢 量角度指数和谱分析法的 T 波交替幅值对 T 波交替的检测结果进行比较， 判断矢量角度指 数检测 T 波交替的可靠性。 0010 为实现上述目的， 本发明采用下述技术方案 ： 0011 一种基于形态学的 T 波交替散点图法的心电信号分析方法， 包括以下步骤 ： 0012 (1) 选取心电信号并进行心电信号预处理 ； 0013 (2)获取采样点 ： 采用T波窗口分析法， 对128个心拍进行采样， 每一个心拍采用T 波窗口法选定 7 个采样点， 就得到检测数据共有 128*7 个采样点 ； 00

16、14 (3) 绘制散点图法 ： 将 128 个相邻心拍的采样点作一次差分， 绘制一次差分的 T 波 散点图 ； 0015 (4) 确定矢量角度指数取值范围， 以判断 T 波交替存在与否 ； 0016 (5) 绘制功率谱曲线图 ： 利用快速傅立叶变换将采样点的幅度转化成能量谱并进 行叠加平均得到其功率谱曲线图 ； 0017 (6) 计算 T 波交替幅度值 (VTWA) ： 0018 对心电信号中 N 个连续 ST-T 心动周期的一组采样点进行 FFT 快速傅里叶变换分 析， 其功率谱的计算公式如下 ： 0019 0020 其中， 0 l 127， N 128， m 7， bi为第 i 个样点的值

17、。 0021 设定 0.46 0.49 周期 / 心拍内的噪声视为背景噪声， 其平均频率值用 Sn 表示 ； 0.5 周期 / 心拍处的频率值用 S0.5表示， T 波交替的功率值记作 S TWA， STWA的计算公式如下 ： 0022 STWA S 0.5-Sn 0023 规定 ： 当STWA0， T 波交替幅度 值为判定 T 波交替存在 ； 0024 通过观察 0.5cycles/beat 处是否存在一明显的谱线， 即为 T 波交替的标识谱线， 可判断在心动周期内是否存在 T 波交替 ： 若线谱不明显， 则心跳周期中不存在 T 波交替 ； 若 有明显线谱， 则存在 T 波交替。 说 明 书

18、 CN 103006207 B 5 3/8 页 6 0025 所述步骤 (3)、 (4)、 (6) 均通过 MATLAB7.0 软件实现。 0026 所述步骤 (1) 中的心电信号是从美国 MIT/BIH 标准心律失常数据库和欧洲 ST-T 心电数据库中获得。 0027 步骤 (1) 中所述心电信号预处理是指去除基线漂移、 工频干扰及肌电干扰的随机 干扰， 获得具有明显特征的心电信号。 0028 步骤 (4) 中所述判断 T 波交替存在与否的散点图矢量角度指数取值范围的确定， 包含以下步骤 ： 0029 41) 根据散点聚集区域的不同形态以及散点图的各种定量参数值， 可得到所需的 不同的信息

19、； 通过计算散点图的矢量角度指数(VAI)来判断T波交替的存在与否， 其计算公 式为 ： 0030 0031 式中， N 128， xi为第 i 点的横坐标， i为第 i 点到坐标原点的连线与横坐标之 间的夹角度， 即 i arctan(x i+2-xi+1)/(xi+1-xi) ； 0032 用软件仿真时， 将 45 度化成弧度制取为 0.7854， VAI 的大小表示相邻 T 波间的幅 值差沿 45线的散开程度 ； 0033 42) 规定 ： 当矢量角度指数 0.9 rad VAI 1.03rad 时， 则在心电信号中， 含有 T 波交替 ； 当 VAI 0.9rad 或 VAI 1.03

20、rad 时， 则在心电信号中， 不含 T 波交替。 0034 最后， 进行了仿真结果分析 ： 首先用 MATLAB7.0 软件和对谱分析法和散点图法所 得矢量角度指数和 T 交替幅值的数据进行对比分析， 检验散点图中矢量角度指数的可靠 性 ； 再用 Kappa 一致性检验进行分析谱分析法与散点图法各个指标检测结果的一致性， 检 验散点图中矢量角度指数的可靠性。 0035 所述仿真结果分析， 具体包含以下步骤 ： 0036 a) 选取美国 MIT/BIH 标准心律失常数据库和欧洲 ST-T 心电数据库的数据， 并连 续取 128 个周期心电信号进行软件仿真检测。MIT/BIH 标准心律失常数据库

21、原始抽样频率 为 360Hz， 而欧洲 ST-T 心电数据库的频率为 250Hz， 故为分析数据方便， 对两组数据库的数 据进行再抽样， 抽样频率为 200Hz ； 0037 b) 用 Curve Fitting Tool 工具对仿真所得的 VTWA和 VAI 的离散数据进行二次曲 线拟合， VTWA和 VAI 之间关系式如 ： 0038 VTWA 4300*VAI 2-7185*VAI+3001 0039 VTWA和 VAI 二次曲线的拟合度为 0.867。V TWA和 VAI 之间的相关系数， 两者之间的 相关系数 0.8601。由于 VT 波交替和 VAI 之间的互相关系数 0.8601

22、， 这表明由 散点图法得到的VAI与非常成熟的谱分析法得到的VT波交替具有很强的相关性， 同时这也 说明矢量角度 VAI 能够很好应用到 T 波交替的检测中。 0040 本发明的优点是 ： 0041 1.本发明采用矢量角度指数， 即散点图每一点到原点的连线与X轴构成的角度Ai 与 45 度之差的和的平均值， 来判断 T 波交替的存在与否， 故该方法与 T 波幅值变化率有关， 进而减小了因 T 波幅值变化而使检测准确率下降的可能性。 说 明 书 CN 103006207 B 6 4/8 页 7 0042 2. 本发明将矢量角度指数应用于实际心电进行检测， 并找出了判断 T 波交替存在 与否的矢量

23、角度指数合适阈值。 0043 3. 本发明采用非线性方法中的散点图法， 非线性法的抗干扰能力较强， 当存在异 常心搏、 急剧变化、 人为干扰以及 “大尾巴” 的生理噪声时， 非线性方法的检测结果比较理 想， 鲁棒性很好。 0044 4. 本发明采用用散点图的矢量角度指数对 T 波交替进行定量分析， 计算比较简 单， 易于硬件操作。 0045 5. 本发明不仅找出了判断 T 波交替存在与否的矢量角度指数合适阈值， 还对其可 靠性作出评估。 0046 6. 具有直观形态且计算简单， 由于是差值作图， 抗干扰性较强， 对信号没有高质量 要求， 不仅可以测得相邻 T 波交替幅值， 而且包含时域信息。

24、附图说明 0047 图 1 为不存在和存在 T 波交替散点图 (a) 存在 T 波交替散点图 (b) 不存在 T 波交 替的散点图 ； 0048 图 2 为不存在 T 波交替交替的功率谱图和存在 T 波交替交替的功率谱图 (a) 不存 在 T 波交替的功率谱图 (b) 存在 T 波交替交替的功率谱图 ； 0049 图 3 为 VTWA和 VAI 的曲线拟合图。 具体实施方式 0050 下面结合附图和实施例对本发明进行进一步的阐述， 应该说明的是， 下述说明仅 是为了解释本发明， 并不对其内容进行限定。 0051 一种基于形态学的 T 波交替散点图法的心电信号分析方法， 包括以下步骤 ： 005

25、2 (1) 选取心电信号并进行心电信号预处理 ： 从美国 MIT/BIH 标准心律失常数据库 和欧洲 ST-T 心电数据库中获得心电信号， 进行去除基线漂移、 工频干扰及肌电干扰的随机 干扰， 获得具有明显特征的心电信号。用简单整系数法去掉 50Hz 工频干扰， 利用零相位数 字滤波器去除心电信号的基线漂移以及采用 bior2.2 小波函数进行小波变换阈值去噪法 去除肌电干扰， 进而获得具有明显特征的心电信号。 0053 (2) 获取采样点 ： 0054 我们采用 T 波窗口分析法， 以便于更准确地检测微伏 TWA(MTWA)。由 RR 间期的大 小来选择 T 波窗口起点及窗宽， T 波起点按

26、经验选取见表 1， 窗宽取 0055 表 1 0056 RR 间期 (s)RR 0.60.6RR 1.1 T波起点(距R峰) 60ms100ms150ms 0057 根据 QT 和 RR 间期的关系式， 计算出窗口的终点， 如公式所示 : 0058 说 明 书 CN 103006207 B 7 5/8 页 8 0059 其中， Q 表示 Q 波起点， 近似在 R 峰前 50ms 处 ； T 表示 T 波起点 ； RR 表示间期。k 是 常数， 一般取 k 0.4。 0060 采用 T 波窗口分析法， 对 128 个心拍进行采样， 每一个心拍采用 T 波窗口法选定 7 个采样点， 就得到检测数据

27、点集合， 共有128*7个采样点 ； 本发明采用T波窗口分析法， 以便 于更准确地检测微伏 T 波交替 (MT 波交替 )， 即选取 L 128 个连续的心跳周期， 每一个周 期采用 T 波窗口法选定一组 m 7 个采样点。 0061 在一个心跳周期内， 将T波分析窗口平均分配， 每两个采样数据间的间距为ID。 采 样点 bi由如下公式得到， 公式为 ： 0062 bi+1 b 1+i*ID (i 1， 26) 0063 利用 QT 和 bi+1的公式， 若每一个心跳周期采样点数为 7 个， 且连续选取 L 个心跳 周期， 可得到7*L个采样点， 进而我们得到一组信号采样点序列Px1,x2,.

28、,xi,.， 取 序列中的采样点的差值组成新的序列 Q x2-x1,x3-x2,.,xi-xi-1,.。 0064 (3) 绘制散点图 ： 将 128 个相邻心拍的采样点作一次差分， 绘制一次差分的 T 波散 点图 ； 得出 (xi+2-xi+1) 和 (xi+1-xi) 的散点图， 其中 i 是识别心跳序数指数。根据散点聚集区 域的不同形态以及散点图的各种定量参数值， 可得到所需的不同信息。 0065 (4) 确定矢量角度指数取值范围， 以判断 T 波交替存在与否的 ； 0066 矢量角度指数 (VAI) 的计算公式为 ： 0067 0068 式中， N 128， xi为第 i 点的横坐标，

29、 i为第 i 点到坐标原点的连线与横坐标之 间的夹角度， 即 i arctan(xi+2-xi+1)/(xi+1-xi)。用软件仿真时， 将 45 度化成弧度制取 为 0.7854。VAI 的大小表示相邻 T 波间的幅值差沿 45线的散开程度。当矢量角度指数 0.9rad VAI 1.03rad 时， 则在心电信号中， 含有 T 波交替 ； 当 VAI 0.9rad 或 VAI 1.03rad 时， 则在心电信号中， 不含 T 波交替。如图 1 不存在 T 波交替散点图和存在 T 波交 替散点图 (a) 存在 T 波交替散点图 (b) 不存在 T 波交替的散点图。 0069 (5) 绘制功率谱

30、曲线图 ： 利用快速傅立叶变换将采样点的幅度转化成能量谱并进 行叠加平均得到其功率谱曲线图 ； 将采样点的能量谱对应叠加法， 以提高信噪比， 从而能够 更加明显地显现出 T 波功率谱图中的 T 波交替。谱分析法是一种趋近于标准化的 T 波交替 频域检测方法， 也是目前应用最广泛的方法。 0070 (6) 计算 T 波交替幅度值 (VTWA) ： 0071 对心电信号中 N 个连续 ST-T 心动周期的一组采样点进行 FFT 快速傅里叶变换分 析， 其功率谱的计算公式如下 ： 0072 0073 其中， 0 l 127， N 128， m 7， bi为第 i 个样点的值。 0074 设定 0.4

31、6 0.49 周期 / 心拍内的噪声视为背景噪声， 其平均频率值用 Sn 表示 ； 说 明 书 CN 103006207 B 8 6/8 页 9 0.5 周期 / 心拍处的频率值用 S0.5表示， T 波交替的功率值记作 S TWA， STWA的计算公式如下 ： 0075 STWA S 0.5-Sn 0076 规定 ： 当STWA0， T 波交替幅度 值为判定 T 波交替存在。 0077 通过观察 0.5cycles/beat 处是否存在一明显的谱线 ( 即为 T 波交替的标识谱 线 )， 可判断在心动周期内是否存在 T 波交替。若线谱不明显， 则心跳周期中不存在 T 波交 替 ； 若有明显线

32、谱， 则存在 T 波交替。以 100 信号的 00:17:13 00:19 ： 00 时间范围内的功 率谱为例， 如图 2(a) 在 0.5cycles/beat 处不存在一明显的谱线， 故在心跳周期中不存在 T 波交替交替。以 e0105 信号的 00 ： 05 ： 20 00 ： 07 ： 40 的谱为例， 如图 2(b) 在 0.5cycles/ beat 处存在一明显的谱线， 故在心跳周期中存在 T 波交替。 0078 仿真结果分析 ： 0079 本发明矢量角度指数 (VAI) 和 T 波交替幅值 (VTWA) 进行比较。在 MATLAB7.0 软件 的环境下， 用Curve Fitt

33、ing Tool工具对实验中VTWA和VAI的离散数据进行二次曲线拟合， VTWA和 VAI 之间关系式如 ： 0080 VTWA 4300*VAI 2-7185*VAI+3001 0081 VTWA和 VAI 的二次曲线拟合如图 3， 该坐标系中纵坐标为 VAI， 横坐标为 V TWA， 其二 次曲线的拟合度为 0.867。VTWA和 VAI 之间的相关系数， 两者之间的相关系数 0.8601。 由于 VTWA和 VAI 之间的互相关系数 0.8601， 这表明由散点图法得到的 VAI 与由已经非 常成熟的谱分析法得到的VT波交替具有很强的相关性， 同时这也说明矢量角度VAI能够很 好应用到

34、 T 波交替的检测中。 0082 为了从定性的角度验证对于散点图和谱分析两种方法在检验 T 波交替时的一致 性关系， 本发明选用 Kappa 一致性检验。用两种方法分别对 N 个观察对象逐个判断其属于 C 类别中的哪一类， 其判断结果常以 CC 列联表的形式表示， 如表 2 ： 0083 表 2 0084 0085 Kappa 系数的计算公式为 ： 0086 K Po-Pe/1-Pe 0087 其中， Po Aii/N， Pe ai*bi， ai Ai/N， bi Bi/N。P0-Pe 为实际一致率， 说 明 书 CN 103006207 B 9 7/8 页 10 1-Pe 为非机遇一致率，

35、Aii 为 CC 联表中的主对角线上的实际值， N 表示总例数， Ai、 Bi 表示 第 i 行、 第 i 列的边缘值， ai、 bi 表示第 i 行、 第 i 列的边缘概率。目前， 对一致性强度的参 考指标划分有细微差异， Landies 和 Koch 将 Kappa 系数划分为六个区段代表不同强弱的一 致性， 如表 3 ： 0088 表 3 0089 Kappa 系数 K00.0 0.20 0.20 0.40 0.40 0.60 0.60 0.80 0.80 1.00 一致性程度极差微弱弱中度高度极强 0090 本发明采用更为细致的六级划分来判断。由于实际资料计算的 K 值只是一个两本 的

36、统计量， 存在抽样误差， 故计算所得 K 值是否来自 K 值为 “0” 的总体即离散数据的一致程 度是由于机遇造成的， 这需要假设检验即 u 检验， 检验公式为 ： 0091 0092 式中， u 为标准正态分位数， Se(K) 为 K 的标准误。用 Kappa 系数的计算原理定性 判断 T 波交替的存在与否， 可将不存在定性为阴性， 用 “-” 表示 ； 将存在定义为阳性， 用 “+” 表示， 就得到了一个两分类， 故可用一个 2*2 联表的形式表示检测结果。通过统计学中的 Kappa一致性检验， 进行分析谱分析法与散点图法各个指标检测结果的一致性， 设H0两种方 法对应的结果无关， H1两

37、种方法对应的结果有关系。谱分析法指标与各个散点图法指标的 联表如表 4 所示 ： 0093 表 4 0094 0095 P0 0.8222， Pe 0.5238， K 0.6266， 参考评价原则， K (0.4,0.75) 一致性 好， 同时， 为 9.2419， 大于 95标准正态分布位数 1.96， 故 P0.05， 拒绝 H0， 接受 H1， 可 认为两种检测方法具有一致性。 0096 以上两种评价方法表明 VAI 与 VTWA具有很强的相关性， 矢量角度 VAI 能够很好应 用到 T 波交替的检测中。 0097 上述虽然结合附图对本发明的具体实施方式进行了描述， 但并非对本发明保护范 围的限制， 在本发明的技术方案的基础上， 本领域技术人员不需要付出创造性劳动即可做 说 明 书 CN 103006207 B 10 8/8 页 11 出的各种修改或变形仍在本发明的保护范围以内。 说 明 书 CN 103006207 B 11 1/2 页 12 图 1 图 2 说 明 书 附 图 CN 103006207 B 12 2/2 页 13 图 3 说 明 书 附 图 CN 103006207 B 13