一种基于单通道ECG信号的心肺耦合分析方法技术

本发明专利技术公开了一种基于单通道ECG信号的心肺耦合分析方法，具体包括：从体表单导联ECG信号中提取RR间期序列和EDR信号；使用三次样条插值将RR间期序列和EDR序列插值至相同采样率；计算RR间期序列和EDR序列的自相关函数和互相关函数；进行傅里叶变换分别得到自相关功率和互相关功率；通过互相关功率计算心肺耦合结果，得到心肺耦合图谱；根据心肺耦合图谱对睡眠质量进行分析。本发明专利技术通过互相关功率计算心肺耦合结果，进而得到心肺耦合图谱对睡眠质量进行分析，提高了睡眠质量分析准确度，可用于居家检测反应真实睡眠状况，便于对使用者睡眠好坏进行准确判断，有利于制定个性化调理方案。案。案。

【技术实现步骤摘要】
一种基于单通道ECG信号的心肺耦合分析方法


[0001]本专利技术涉及医疗健康检测
，具体涉及一种基于单通道ECG信号的心肺耦合分析方法。

技术介绍

[0002]上个世纪脑电波的发现让睡眠研究人员十分兴奋，并把脑电波定为评价睡眠的核心讯号，后来又增加了眼电、肌电、呼吸、心电、血氧等等多种生理信号的采集，形成了目前医院里使用的多导睡眠监测仪(PSG)。多导睡眠仪虽有它的优点，但同时它也有很多的缺点，比如，被试者需要在身上连接各种导线，并要在医院的睡眠实验室睡一个晚上。这样的环境和条件往往无法反映被试者平时真实的睡眠；另一大问题是必须依赖有经验的技术人员花几个小时去人工识别出脑电图中的特征波形，对睡眠进行分期。
[0003]睡眠是全身整体的生理过程，大脑放松的同时，肌肉也是在放松的状态，心率会减慢、呼吸也会减慢并变得更平稳。也就是说，在睡眠中，不仅只有大脑进入特殊的状态，身体的其他系统也会表现出睡眠的特征。经过对大量临床数据的分析，医学院专家们发现，在睡眠的时候，心电和呼吸还存在一种奇妙的耦合关系。当人在熟睡的时候，这种耦合性会增强。而在清醒、浅睡、或者一些疾病状态下，这种耦合性也会呈现出不同的特点。基于此发现，医学院的睡眠专家团队在2005年创新性地提出了Cardiopulmonary Coupling，也就是心肺耦合的概念，简称CPC，并且发展出成熟的技术，可以准确地给出睡眠的分期。通过临床验证，这些结果与传统基于脑电波方法判断的睡眠分期结果有很好的一致性。
[0004]在一整夜的睡眠过程中，人体会经历从清醒，浅睡，到熟睡的一个周期性循环，其中还会穿插若干段与梦境相关的快速动眼睡眠(REM)，每一个睡眠阶段都应该在一个正常的时间范围内。CPC可以判断每一个睡眠阶段发生的时间和比例，准确的刻画出睡眠的结构，提示用户是否存在睡眠的问题。目前已有的心肺耦合检测方法较为基础，还需进一步完善。
[0005]如中国专利CN104757967A，公开日2015年7月8日，一种心肺耦合反馈方法及其装置，其包括：第一步：通过将受试者的心血管信号和呼吸信号进行耦合分析，得到受试者的心肺耦合强度；第二步：根据所述受试者的心肺耦合强度向受试者输出运动指示信号，以提高受试者的心肺耦合强度，使用者在心肺耦合强度的指导下作特定运动，不断提高心肺耦合强度。其存在检测手段较为简单，睡眠质量分析准确度不高的问题。

技术实现思路

[0006]本专利技术要解决的技术问题是：目前的心肺耦合分析方法存在睡眠质量分析准确度不高的技术问题。提出了一种能够通过心肺耦合图谱提高睡眠质量分析准确度的基于单通道ECG信号的心肺耦合分析方法。
[0007]为解决上述技术问题，本专利技术所采取的技术方案为：一种基于单通道ECG信号的心肺耦合分析方法，包括如下步骤：
S1：从体表单导联ECG信号中提取RR间期序列和EDR信号；S2：使用三次样条插值将RR间期序列和EDR序列插值至相同采样率；S3：计算RR间期序列和EDR序列的自相关函数和互相关函数；S4：进行傅里叶变换分别得到自相关功率和互相关功率；S5：通过互相关功率计算心肺耦合结果，得到心肺耦合图谱；S6：根据心肺耦合图谱对睡眠质量进行分析。
[0008]RR间期为心电图上QRS波群中R波和R波之间的距离，信号为采样率500Hz的ECG信号时，需将ECG信号降采样至100Hz；再对降采样之后的信号进行10Hz～30Hz的带通滤波；之后对ECG信号进行标准化操作。
[0009]作为优选，步骤S1中提取RR间期序列的过程包括自动识别ECG信号中所有R波的位置的过程，识别ECG信号中所有R波的位置包括如下步骤：A1：判断ECG信号是否倒置；A2：计算ECG
3000
的一阶差分ECG
‑
diff
3000
，将ECG
‑
diff
3000
划分为长度为1秒的片段，求出每个片段的最小值的均值mean
diff
，将R波的阈值设为μ
r
＝0.4*mean
diff
；A3：计算ECG信号的一阶差分ECG
‑
diff，找出ECG
‑
diff的所有极值点，得到R波位置序列I
R
。
[0010]ECG信号倒置为R波峰朝下，可能由于电极贴反导致。
[0011]作为优选，步骤A1包括如下步骤：A11：取ECG信号中前3000秒信号ECG
3000
，将其分为3000个长度为1秒的片段；A12：找出每1秒片段的最大值索引i
max
以及最小值索引i
min
；A13：将3000个片段中满足i
max
＜i
min
和的片段的总数记为n,若则ECG信号方向正确，否则倒置ECG信号。
[0012]fz为ECG信号的采样率，正常1秒心电信号中最大值对应R波峰，最小值对应S波，R波应在S波之前，且间距一般小于0.3秒。
[0013]作为优选，所述步骤S1包括如下步骤：”B1：计算I
R
中相邻两个R波的中点得到I
R
，对I
R
进行插值，使得RR间期序列与EDR序列保持相同采样率；B2：计算I
R
相邻R波的距离得到间期序列RR；B3：根据每个R波位置对应的振幅值得到心电图派生呼吸序列EDR。
[0014]I
R
'表示ECG信号的基线，ECG'＝ECG
‑
I
r
'为去除基线漂移的心电信号，R波振幅的为ECG'(R
i
)
‑
ECG'(S
i
)，即ECG'中R波位置对应的值减S波位置对应的值。
[0015]作为优选，所述步骤S3包括如下步骤：S31：计算RR序列和EDR序列对应的自相关函数f
xx
，f
yy
，以及RR序列与EDR序列的互相关函数f
xy
；S32：计算f
xx
、f
yy
及f
xy
的傅里叶变换，得到相应的自功率谱S
xx
，S
yy
和互功率谱S
xy
。
[0016]以每300秒为单位计算RR序列和EDR序列的自功率谱S
xx
、S
yy
以及RR与EDR的互功率
谱S
xy
。
[0017]作为优选，步骤S31中的自相关函数公式为：其中N为信号长度互相关函数公式为：其中N为信号长度。
[0018]作为优选，步骤S32中采用的离散傅里叶变换公式为：其中N为信号长度(0≤k≤N
‑
1)。
[0019]作为优选，步骤S5中计算心肺耦合结果的计算式为：
[0020]通过本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于单通道ECG信号的心肺耦合分析方法，其特征在于，包括如下步骤：S1：从体表单导联ECG信号中提取RR间期序列和EDR信号；S2：使用三次样条插值将RR间期序列和EDR序列插值至相同采样率；S3：计算RR间期序列和EDR序列的自相关函数和互相关函数；S4：进行傅里叶变换分别得到自相关功率和互相关功率；S5：通过互相关功率计算心肺耦合结果，得到心肺耦合图谱；S6：根据心肺耦合图谱对睡眠质量进行分析。2.根据权利要求1所述的一种基于单通道ECG信号的心肺耦合分析方法，其特征在于，步骤S1中提取RR间期序列的过程包括自动识别ECG信号中所有R波的位置的过程，识别ECG信号中所有R波的位置包括如下步骤：A1：判断ECG信号是否倒置；A2：计算ECG
3000
的一阶差分ECG
‑
diff
3000
，将ECG
‑
diff
3000
划分为长度为1秒的片段，求出每个片段的最小值的均值mean
diff
，将R波的阈值设为μ
r
＝0.4*mean
diff
；A3：计算ECG信号的一阶差分ECG
‑
diff，找出ECG
‑
diff的所有极值点，得到R波位置序列I
R
。3.根据权利要求2所述的一种基于单通道ECG信号的心肺耦合分析方法，其特征在于，步骤A1包括如下步骤：A11：取ECG信号中前3000秒信号ECG
3000
，将其分为3000个长度为1秒的片段；A12：找出每1秒片段的最大值索引i
max
以及最小值索引i
min
；A13：将3000个片段中满足...

【专利技术属性】
技术研发人员：袁靖嘉，戴松青，王凯平，罗睿，
申请(专利权)人：杭州菲诗奥医疗科技有限公司，
类型：发明
国别省市：