CN108175406A - 一种治疗中枢性睡眠呼吸暂停症的方法及装置 - Google Patents
一种治疗中枢性睡眠呼吸暂停症的方法及装置 Download PDFInfo
- Publication number
- CN108175406A CN108175406A CN201711361057.4A CN201711361057A CN108175406A CN 108175406 A CN108175406 A CN 108175406A CN 201711361057 A CN201711361057 A CN 201711361057A CN 108175406 A CN108175406 A CN 108175406A
- Authority
- CN
- China
- Prior art keywords
- signal
- respiratory muscle
- envelope
- value
- sleep apnea
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 208000003417 Central Sleep Apnea Diseases 0.000 title abstract description 8
- 210000003019 respiratory muscle Anatomy 0.000 claims abstract description 56
- 238000001914 filtration Methods 0.000 claims abstract description 20
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 19
- 239000001301 oxygen Substances 0.000 claims abstract description 19
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 19
- 239000008280 blood Substances 0.000 claims abstract description 17
- 210000004369 blood Anatomy 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 201000002859 sleep apnea Diseases 0.000 claims description 27
- 238000002567 electromyography Methods 0.000 claims description 24
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 17
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 14
- 210000003205 muscle Anatomy 0.000 claims description 13
- 210000000876 intercostal muscle Anatomy 0.000 claims description 11
- 210000003489 abdominal muscle Anatomy 0.000 claims description 8
- 238000000338 in vitro Methods 0.000 claims description 8
- 238000000205 computational method Methods 0.000 claims description 6
- 230000033001 locomotion Effects 0.000 claims description 5
- 210000000062 pectoralis major Anatomy 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 2
- 230000000747 cardiac effect Effects 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000003183 myoelectrical effect Effects 0.000 abstract description 3
- 230000000241 respiratory effect Effects 0.000 abstract description 2
- 108010064719 Oxyhemoglobins Proteins 0.000 abstract 2
- 239000000523 sample Substances 0.000 abstract 1
- 210000000038 chest Anatomy 0.000 description 6
- 208000014486 central sleep apnea syndrome Diseases 0.000 description 5
- 230000003434 inspiratory effect Effects 0.000 description 5
- 206010019280 Heart failures Diseases 0.000 description 4
- 208000008784 apnea Diseases 0.000 description 4
- 208000011580 syndromic disease Diseases 0.000 description 4
- 238000002560 therapeutic procedure Methods 0.000 description 4
- 230000008602 contraction Effects 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 206010020591 Hypercapnia Diseases 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000024891 symptom Diseases 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000000115 thoracic cavity Anatomy 0.000 description 2
- 238000009423 ventilation Methods 0.000 description 2
- 206010008501 Cheyne-Stokes respiration Diseases 0.000 description 1
- 206010012289 Dementia Diseases 0.000 description 1
- 208000007590 Disorders of Excessive Somnolence Diseases 0.000 description 1
- 206010041349 Somnolence Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000004070 electrodeposition Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 210000002216 heart Anatomy 0.000 description 1
- 230000003601 intercostal effect Effects 0.000 description 1
- 210000000936 intestine Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 235000013372 meat Nutrition 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 201000001119 neuropathy Diseases 0.000 description 1
- 230000007823 neuropathy Effects 0.000 description 1
- 230000000414 obstructive effect Effects 0.000 description 1
- 208000001797 obstructive sleep apnea Diseases 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000008289 pathophysiological mechanism Effects 0.000 description 1
- 208000033808 peripheral neuropathy Diseases 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000003519 ventilatory effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/3601—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation of respiratory organs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/36014—External stimulators, e.g. with patch electrodes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pulmonology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种治疗中枢性睡眠呼吸暂停症的方法及装置，该方法包括以下步骤：1)、采集呼吸肌的表面肌电信号；2)、滤波滤除；3)、进行表面肌电信号的包络计算；4)产生呼吸肌电标示信号；5)、检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头；6)、如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过3％，则对呼吸肌进行体外起搏，以及应用该方法的装置，本发明精确有效，使用效果好。
Description
技术领域
本发明涉及医疗技术领域，尤其涉及一种治疗中枢性睡眠呼吸暂停症的方法和装置。
背景技术
中枢性睡眠呼吸暂停(CSA)特点是睡眠时因中枢驱动功能受损而引起的反复气流中断。中枢性呼吸暂停是指呼吸暂停时间≥10秒，无呼吸运动。通常这种情况每小时＞5次就认为是异常的。CSA综合征是指睡眠期间平均每小时出现5 次以上中枢性呼吸暂停，并且出现睡眠片段(频繁觉醒)相关症状和/或白天过度嗜睡。阻塞性呼吸暂停患者也可发生中枢性呼吸暂停，因此应着重辨别CSA，尽管在这个问题上还没有绝对标准，通常对CSA患者的研究要求50％甚至80％以上的症状要符合中枢性。CSA通常非单一因素诱发，因此出现了大量的综合征，且各自存在一些潜在的病理生理学机制。潮式呼吸(CSR)与神经病变，如脑血管疾病、痴呆等有关，亦常发生在心力衰竭患者中，特点是呼吸呈渐强- 渐弱改变，在呼吸运动终末出现中枢性呼吸暂停或低通气。心力衰竭患者发生 CSR是由于循环时间延长以及通气控制系统高增益对高碳酸血症反应增强的联合作用所致。这种联合作用导致通气控制的不稳定性以及特定形式的周期性呼吸。原发性CSA特点是因高碳酸血症引起深度通气反应导致的通气不稳。
中枢性睡眠呼吸暂停综合征是一种严重影响工作和生活的疾病，目前的治疗手段极其有限，主要为完善心力衰竭治疗、气道正压通气、夜间吸氧等。完善心力衰竭治疗、夜间吸氧对于中枢性睡眠呼吸暂停综合征的确切疗效还没有明确的结论。气道正压通气治疗要求患者夜间使用睡眠呼吸机，但大多数患者对于使用睡眠呼吸机感觉舒适性很差，不能长期坚持。
人体的呼吸相关肌肉主要有膈肌、肋间外肌、胸大肌、腹肌、胸锁乳突肌等，其中膈肌和肋间外肌是最主要的吸气肌。吸气时，膈肌收缩，膈顶下将，胸腔增大；呼气时，膈肌舒张，膈顶上升，胸腔缩小。吸气时，肋间外肌收缩，肋骨向上向外运动，体积增大；呼气时，肋间外肌舒张，肋骨向下向内运动，体积缩小。用力吸气时，除了膈肌、肋间外肌的收缩，胸大肌、胸锁乳突肌等发生收缩，参与扩张胸廓。用力呼气时，除了膈肌、肋间外肌的舒张，肋间内肌、腹肌等发生收缩，参与收缩胸廓。
研究表明，吸气期间吸气肌的肌电信号出现，随着吸气努力的增加，吸气肌的肌电信号也增强；而吸气肌的肌电信号随吸气的结束和呼气的开始而消失。呼气期间呼气肌的肌电信号出现；而呼气肌的肌电信号随呼气的结束和吸气的开始而消失。如果能够采集到呼吸肌的肌电信号并对肌电活动进行计时，就可以初步判断有无呼吸暂停。
膈肌是最主要的呼吸肌，但表面膈肌肌电信号非常微弱，健康个体在平静呼吸时表面电极所测膈肌肌电幅值约为10～100μV，极易受多种噪声的影响。同时，由于解剖的差异性，膈肌肌电受电极与膈肌的距离以及膈肌肌纤维分布密度等的影响，不同个体之间、同一个体不同时间段的膈肌肌电幅值存在较大差别。为保证高精度测量膈肌肌电信号并预留一定裕量，所需测量范围将高达数千倍以上。
膈肌肌电信号微弱并受多种因素的影响(如心脏、胃肠等肌性器官及胸腹肌肉的电活动、电极位置移动等的干扰)，获取稳定的膈肌肌电信号难度较大。肋间外肌和胸大肌、腹肌的表面肌电信号采集时也受心电信号的很大干扰。
由于表面肌电图采集呼吸相关肌肉的肌电信号时受心电等信号的干扰很严重，难于提取到纯洁的肌电信号，使得呼吸肌表面肌电图应用受限。
中国专利申请号CN 2574596A公开了一种睡眠呼吸障碍快捷诊断监护装置，旨在安全、快捷、方便地获得睡眠呼吸暂停综合症特征病理参数，为实现睡眠呼吸暂停综合症普查提供一种有效的技术手段，为睡眠呼吸暂停综合症患者提供必要的监护。它包括两个探测呼吸运动的体位传感器、呼吸信号处理器、双信号相关分析器、记录显示呼吸参数的终端、可控电路开关和报警救助装置。患者在睡眠中出现呼吸障碍时，实时记录显示关键病理参数，并自动提供必要的救助。这种装置不扰动自然呼吸过程，能够大幅度降低睡眠呼吸暂停综合症诊断装置的制造成本，还可以成为家用健康监测装置。
CN 202505987U公开一种一种基于体外膈肌起搏的睡眠呼吸暂停治疗装置，由传感器、膈肌起搏器和电极组成，其中，传感器与膈肌起搏器连接，膈肌起搏器与电极连接。使用时，当气流传感器、声音传感器、腹部运动传感器或/和氧饱和度传感器检测到睡眠呼吸暂停后，输出控制信号至膈肌起搏器，脉冲信号产生器按照预置的程序发出脉冲信号，通过体表的膈神经附近的电极传送电脉冲，膈神经受到脉冲刺激后传递至膈肌，使膈肌收缩，呼吸恢复正常。本实用新型对中枢性和阻塞性睡眠呼吸暂停均有良好的效果，对中枢性睡眠呼吸暂停尤为明显。
CN 103458781A公开了一种诊断患者的中枢性睡眠呼吸暂停的方法和系统，采用几乎不结合呼吸参数的相对短的测试。与通常使用的耗时且复杂的多参数睡眠研究相比，所述测试时间短。所述患者使用气道正压通气设备，并且遵从呼吸提示。在传送呼吸提示期间，中枢性呼吸暂停的发生可以支持对中枢性睡眠呼吸暂停的诊断。
目前这些技术来相对来说都比较简单，没有实现与人体自然呼吸的同步，有可能与人体自然呼吸冲突，使治疗效果变差。表面膈肌肌电信号非常微弱，健康个体在平静呼吸时表面电极所测膈肌肌电幅值约为10～100μV，极易受多种噪声的影响。没有实现表面肌电信号应用于中枢性睡眠呼吸暂停综合征的治疗。
发明内容
为了克服现有技术的不足，本发明的目的之一在于提供一种治疗中枢性睡眠呼吸暂停症的方法，其能够解决现有技术中不能够将表面肌电信号应用到中枢性睡眠呼吸暂停综合征的治疗中。
本发明的目的之二在于提供一种治疗中枢性睡眠呼吸暂停症的装置，其能够解决现有技术中不能够将表面肌电信号应用到中枢性睡眠呼吸暂停综合征的治疗中。
本发明的目的之一采用如下技术方案实现：
一种治疗中枢性睡眠呼吸暂停症的方法，该方法包括以下步骤：
1)、采集呼吸肌，包括但不限于膈肌、肋间外肌、胸大肌、腹肌等的表面肌电信号，使用仪表放大电路和右腿驱动提高采集放大电路的抗共模干扰能力，进行多级放大；
2)、低通滤波滤除高频噪音，高通滤波滤除极化电压、伪迹和不稳定成分，陷波和梳状滤波去除工频及其谐波干扰；
3)、进行表面肌电信号的包络计算，使用均方根算法实现，包络值大于阈值表示肌肉收缩，阈值根据一定时段包络值的均值确定；
4)、产生呼吸肌电标示信号，呼吸肌电标示信号用来标记呼吸肌的收缩运动，由于心电信号幅值通常大于呼吸肌的肌电信号，但心电的短时过零率远小于肌电的短时过零率，可以用短时过零率来区分心电还是呼吸肌的肌电信号；
5)、检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头；
6)、如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过 3％，则对呼吸肌进行体外起搏。
进一步地，所述呼吸肌电标示信号产生的条件为某一小段时间内的平均包络值大于设定的阈值而且短时过零率大于某个设定值。
一种治疗中枢性睡眠呼吸暂停症的方法，该方法包括以下步骤：
1)、功能需求及参数取值
a、对体表膈肌肌电测量电路的输出信号进行采样和AD转换，精度为12bit，采样频率初选4K/s；
b、对信号作四阶低通滤波和工频陷波处理；
c、计算肌电包络，并形成呼吸肌电标示信号；
2)、滤波、包络计算和同步信号产生算法
a、四阶低通滤波与工频陷波
低通滤波和工频陷波，其脉冲传递函数为：
其中G(z)表示系统输出变量的脉冲序列的z变换，记滤波器输入和输出为 u和y，则其差分方程表示为：
其中，参数可根据巴特沃斯低通滤波器或连续信号滤波器离散化等设计方法来确定，工频陷波仅滤掉50/60Hz附近很窄的频段，基本不影响有效信号；
b、包络计算
包络计算借鉴了均方根算法实现，其计算方法为：
取En(k)作为当前时刻的包络值，包络值大于阈值表示呼吸肌收缩，阈值根据一定时段包络值的均值确定；
c、呼吸肌电标示信号产生
程序自动统计当前若干毫秒内包络的均值，以此作为判断的阈值；程序按帧处理数据，单帧平均包络计算方法为：
其中，En(k-i)为公式(4)所得的包络值；进而根据(6)确定阈值Threshold；
此外，呼吸肌电标示信号产生的条件为：Avg(k)＞Threshold且cnt＞4；
3)、对呼吸肌进行体外起搏的条件
检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头；
如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过3％，则对呼吸肌进行体外起搏，此处所述呼吸肌包括但不限于膈肌、肋间外肌、腹肌。
进一步地，所述步骤2)c中，考虑到心电信号幅值通常大于呼吸肌的肌电信号，但其短时过零率远小于呼吸肌的肌电信号，因此，当信号单帧平均包络超过阈值时，程序进一步运用短时过零率来区分其是由心电还是呼吸肌收缩所引起。
进一步地，所述步骤2)c中，短时过零率cnt计算方式为：单帧内峰值幅度超过阈值c1Threshold而相邻谷值小于0的峰值点数，其中，比例系数c1可取0.3～0.6。
进一步地，所述工频陷波可为选50Hz。
进一步地，所述工频陷波可为选60Hz。
本发明的目的之二采用如下技术方案实现：
一种应用如前所述的治疗中枢性睡眠呼吸暂停症的方法的装置。
进一步地，所述装置还加设有显示驱动电路、用于显示工作状态的显示器和微处理器，所述微处理器具有显示信号输出端，该显示信号输出端经显示驱动电路与显示器相连。
进一步地，所述装置还具有电源管理模块。
相比现有技术，本发明的有益效果在于：
本发明成功地实现了表面肌电信号应用于中枢性睡眠呼吸暂停综合征的治疗。
具体实施方式
下面，结合具体实施方式，对本发明做进一步描述，需要说明的是，在不相冲突的前提下，以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。
一个应用本发明具体的实施例：
一、功能需求及参数取值
1)对体表呼吸肌肌电测量电路的输出信号进行采样和AD转换，例如：精度12bit，采样频率为4K/s。
2)对信号作四阶低通滤波(例如：截止频率为600Hz)和工频陷波(可选 50Hz和60Hz)处理。
3)高通滤波滤除极化电压、伪迹和不稳定成分，例如：截止频率为10Hz。
4)计算肌电包络，并形成呼吸肌电标示信号。
二、滤波、包络计算和呼吸肌电标示信号产生算法
2.1、四阶低通滤波与工频陷波
低通滤波和工频陷波，其脉冲传递函数为：
分别记滤波器输入和输出为u和y，则其差分方程表示为：
其中，参数可根据巴特沃斯低通滤波器或连续信号滤波器离散化等设计方法来确定。工频陷波仅滤掉50/60Hz附近很窄的频段，基本不影响有效信号。
2.2、包络计算
包络计算借鉴了均方根算法实现，其计算方法为：
取En(k)作为当前时刻的包络值。包络值大于阈值表示呼吸肌收缩。阈值根据一定时段包络值的均值确定。
2.3、呼吸肌电标示信号产生
程序自动统计当前若干毫秒内包络的均值，以此作为判断的阈值。程序按帧处理数据，单帧平均包络计算方法为：
其中，En(k-i)为式(4)所得包络值。进而根据(6)确定阈值Threshold。
此外，考虑到心电信号幅值通常大于呼吸肌的肌电信号，但其短时过零率远小于呼吸肌的肌电信号。因此，当信号单帧平均包络超过阈值时，程序进一步运用短时过零率来区分其是由心电还是呼吸肌收缩所引起。短时过零率cnt 计算方式为：单帧内峰值幅度超过阈值c1Threshold而相邻谷值小于0的峰值点数，其中，比例系数c1可取0.3～0.6。
呼吸肌电标示信号产生的条件为：Avg(k)＞Threshold且cnt＞4。
三、对呼吸肌进行体外起搏的条件
检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头。
如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过3％，则对呼吸肌进行体外起搏。此处所述呼吸肌包括但不限于膈肌、肋间外肌、腹肌。
上述实施方式仅为本发明的优选实施方式，不能以此来限定本发明保护的范围，本领域的技术人员在本发明的基础上所做的任何非实质性的变化及替换均属于本发明所要求保护的范围。
Claims (10)
1.一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：该方法包括以下步骤：
1)、采集呼吸肌，包括但不限于膈肌、肋间外肌、胸大肌、腹肌等的表面肌电信号，使用仪表放大电路和右腿驱动提高采集放大电路的抗共模干扰能力，进行多级放大；
2)、低通滤波滤除高频噪音，高通滤波滤除极化电压、伪迹和不稳定成分，陷波和梳状滤波去除工频及其谐波干扰；
3)、进行表面肌电信号的包络计算，使用均方根算法实现，包络值大于阈值表示肌肉收缩，阈值根据一定时段包络值的均值确定；
4)、产生呼吸肌电标示信号，呼吸肌电标示信号用来标记呼吸肌的收缩运动，由于心电信号幅值通常大于呼吸肌的肌电信号，但心电的短时过零率远小于肌电的短时过零率，可以用短时过零率来区分心电还是呼吸肌的肌电信号；
5)、检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头；
6)、如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过3％，则对呼吸肌进行体外起搏。
2.如权利要求1所述的治疗中枢性睡眠呼吸暂停症的方法，其特征在于：所述呼吸肌电标示信号产生的条件为某一小段时间内的平均包络值大于设定的阈值而且短时过零率大于某个设定值。
3.一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：该方法包括以下步骤：
1)、功能需求及参数取值
a、对体表膈肌肌电测量电路的输出信号进行采样和AD转换，精度为12bit，采样频率初选4K/s；
b、对信号作四阶低通滤波和工频陷波处理；
c、计算肌电包络，并形成呼吸肌电标示信号；
2)、滤波、包络计算和同步信号产生算法
a、四阶低通滤波与工频陷波
低通滤波和工频陷波，其脉冲传递函数为：
其中G(z)表示系统输出变量的脉冲序列的z变换，记滤波器输入和输出为u和y，则其差分方程表示为：
其中，参数可根据巴特沃斯低通滤波器或连续信号滤波器离散化等设计方法来确定，工频陷波仅滤掉50/60Hz附近很窄的频段，基本不影响有效信号；
b、包络计算
包络计算借鉴了均方根算法实现，其计算方法为：
取En(k)作为当前时刻的包络值，包络值大于阈值表示呼吸肌收缩，阈值根据一定时段包络值的均值确定；
c、呼吸肌电标示信号产生
程序自动统计当前若干毫秒内包络的均值，以此作为判断的阈值；程序按帧处理数据，单帧平均包络计算方法为：
其中，En(k-i)为公式(4)所得的包络值；进而根据(6)确定阈值Threshold；
此外，呼吸肌电标示信号产生的条件为：Avg(k)>Threshold且cnt>4；
3)、对呼吸肌进行体外起搏的条件
检测患者的血氧饱和度，可以使用夹于手指的血氧饱和度探头；
如果连续10秒没有采样到呼吸肌电标示信号且血氧饱和度下降超过3％，则对呼吸肌进行体外起搏，此处所述呼吸肌包括但不限于膈肌、肋间外肌、腹肌。
4.如权利要求1-3所述的一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：所述步骤2)c中，考虑到心电信号幅值通常大于呼吸肌的肌电信号，但其短时过零率远小于呼吸肌的肌电信号，因此，当信号单帧平均包络超过阈值时，程序进一步运用短时过零率来区分其是由心电还是呼吸肌收缩所引起。
5.如权利要求1-4所述的一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：所述步骤2)c中，短时过零率cnt计算方式为：单帧内峰值幅度超过阈值c1Threshold而相邻谷值小于0的峰值点数，其中，比例系数c1可取0.3～0.6。
6.如权利要求1-5所述的一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：所述工频陷波可为选50Hz。
7.如权利要求3所述的一种治疗中枢性睡眠呼吸暂停症的方法，其特征在于：所述工频陷波可为选60Hz。
8.一种应用权利要求1-7所述的治疗中枢性睡眠呼吸暂停症的方法的装置。
9.一种如权利要求8所述的装置，其特征在于：所述装置还加设有显示驱动电路、用于显示工作状态的显示器和微处理器，所述微处理器具有显示信号输出端，该显示信号输出端经显示驱动电路与显示器相连。
10.一种如权利要求8所述的装置，其特征在于：所述装置还具有电源管理模块。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711361057.4A CN108175406A (zh) | 2017-12-14 | 2017-12-14 | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711361057.4A CN108175406A (zh) | 2017-12-14 | 2017-12-14 | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108175406A true CN108175406A (zh) | 2018-06-19 |
Family
ID=62546197
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711361057.4A Pending CN108175406A (zh) | 2017-12-14 | 2017-12-14 | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108175406A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110464944A (zh) * | 2019-08-26 | 2019-11-19 | 深圳创达云睿智能科技有限公司 | 一种呼吸机的控制方法、装置及呼吸机 |
| CN111012341A (zh) * | 2020-01-08 | 2020-04-17 | 东南大学 | 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202505987U (zh) * | 2012-02-13 | 2012-10-31 | 山东大学 | 一种基于体外膈肌起搏的睡眠呼吸暂停治疗装置 |
| CN105879223A (zh) * | 2016-04-22 | 2016-08-24 | 广州雪利昂生物科技有限公司 | 表面肌电信号作为同步信号触发体外膈肌起搏器的方法及装置 |
| CN105943040A (zh) * | 2016-06-03 | 2016-09-21 | 广州雪利昂生物科技有限公司 | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 |
| CN106362282A (zh) * | 2016-09-27 | 2017-02-01 | 广州雪利昂生物科技有限公司 | 一种治疗中枢性睡眠呼吸暂停症的方法和装置 |
-
2017
- 2017-12-14 CN CN201711361057.4A patent/CN108175406A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202505987U (zh) * | 2012-02-13 | 2012-10-31 | 山东大学 | 一种基于体外膈肌起搏的睡眠呼吸暂停治疗装置 |
| CN105879223A (zh) * | 2016-04-22 | 2016-08-24 | 广州雪利昂生物科技有限公司 | 表面肌电信号作为同步信号触发体外膈肌起搏器的方法及装置 |
| CN105943040A (zh) * | 2016-06-03 | 2016-09-21 | 广州雪利昂生物科技有限公司 | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 |
| CN106362282A (zh) * | 2016-09-27 | 2017-02-01 | 广州雪利昂生物科技有限公司 | 一种治疗中枢性睡眠呼吸暂停症的方法和装置 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110464944A (zh) * | 2019-08-26 | 2019-11-19 | 深圳创达云睿智能科技有限公司 | 一种呼吸机的控制方法、装置及呼吸机 |
| CN111012341A (zh) * | 2020-01-08 | 2020-04-17 | 东南大学 | 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 |
| CN111012341B (zh) * | 2020-01-08 | 2022-04-22 | 东南大学 | 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101656611B1 (ko) | 무구속적으로 측정한 생체신호를 이용하여 산소탈포화지수를 획득하는 방법 | |
| CN108135535B (zh) | 用于处理与呼吸活动有关的肌电图信号的装置和方法 | |
| Estrada et al. | Improvement in neural respiratory drive estimation from diaphragm electromyographic signals using fixed sample entropy | |
| Helfenbein et al. | Development of three methods for extracting respiration from the surface ECG: A review | |
| KR101864642B1 (ko) | 수면 무호흡 측정장치 및 그 방법 | |
| CN109414204A (zh) | 用于确定针对对象的呼吸信息的方法和装置 | |
| JP6387352B2 (ja) | 脳系生体信号を検出するシステム | |
| CN105769122B (zh) | 一种睡眠呼吸暂停低通气指数的测量设备 | |
| Nabavi et al. | A robust fusion method for motion artifacts reduction in photoplethysmography signal | |
| US11986314B2 (en) | Apparatus and methods for screening, diagnosis and monitoring of respiratory disorders | |
| CN109770892A (zh) | 一种基于心电信号的睡眠分期方法 | |
| CA2559857A1 (en) | Non-invasive method and device for detecting inspiratory effort | |
| WO2015001188A1 (en) | Method for respiratory measurement | |
| Młyńczak et al. | Impedance pneumography: Is it possible? | |
| Lee et al. | Respiratory rate extraction from pulse oximeter and electrocardiographic recordings | |
| JP6709116B2 (ja) | 呼吸検出装置、呼吸検出方法および呼吸検出用プログラム | |
| Przystup et al. | A multisensor detector of a sleep apnea for using at home | |
| CN105943040A (zh) | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 | |
| CN106362282A (zh) | 一种治疗中枢性睡眠呼吸暂停症的方法和装置 | |
| CN108175406A (zh) | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 | |
| KR101755791B1 (ko) | 공진호흡주파수 및 위상 검출을 통한 최적의 호흡패턴 결정 방법, 시스템 및 컴퓨터 프로그램 | |
| CN114027824B (zh) | 普适性肺通气量与经胸电阻抗的线性模型构建方法及应用 | |
| Seppä | Development and clinical application of impedance pneumography technique | |
| CN207679453U (zh) | 一种治疗中枢性睡眠呼吸暂停症的装置 | |
| Vowles et al. | The Use of Conductive Lycra Fabric in the Prototype Design of a Wearable Device to Monitor Physiological Signals |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180619 |
|
| RJ01 | Rejection of invention patent application after publication |