CN102488521A - 医用呼吸监测仪呼吸阻力值计算装置 - Google Patents
医用呼吸监测仪呼吸阻力值计算装置 Download PDFInfo
- Publication number
- CN102488521A CN102488521A CN2011104087373A CN201110408737A CN102488521A CN 102488521 A CN102488521 A CN 102488521A CN 2011104087373 A CN2011104087373 A CN 2011104087373A CN 201110408737 A CN201110408737 A CN 201110408737A CN 102488521 A CN102488521 A CN 102488521A
- Authority
- CN
- China
- Prior art keywords
- signal
- respiratory resistance
- oscillation
- flow
- component
- 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
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
医用呼吸监测仪呼吸阻力值计算装置,其乃功能模块构架,其所实施的方法对数据采集得到的原始呼吸振荡压力信号和流量信号在时域上进行移动平均滤波得到呼吸压力成分信号和流量成分信号,用原始呼吸振荡压力信号和流量信号减去呼吸压力成分信号和流量成分信号得到经肺组织吸收折射的振荡压力成分信号和流量成分信号;采用振荡成分信号中的基波分量的傅立叶系数近似计算出振荡压力成分信号和流量成分信号的振幅和相位,据此得出其阻抗分量和/或电抗分量。本发明无需进行运算量巨大的傅立叶变换就可得出近似的呼吸阻力值,计算简单,实时性强。
Description
技术领域
本发明涉及基于强迫振荡技术对所采集到的呼吸信号进行数据分析,尤其是用于实施医用呼吸监测仪呼吸阻力值计算方法而建立的功能模块构架。
背景技术
强迫振荡技术(forced oscillation technique,FOT)是无创性评价呼吸***力学特性的方法,具有简便、客观和实时的特点。其工作的原理是,通过振荡发生器产生一个特定频率和振幅的压力振荡,经口施加在受试者的呼吸气流之上,测量经气道及肺组织吸收并折射的振荡压力和振荡流量,利用时域分析技术,计算振荡压力和流量的商即可获得呼吸总阻抗(impedance,Zrs),习惯上称为呼吸阻力。近年来,单频强迫振荡信号用于追踪呼吸阻力随时间的变化,实时监测病情变化及采取相应的治疗措施,是FOT在临床的新应用。呼吸阻力包含了阻抗(resistance,Rrs)和电抗(reactance,Xrs)分量,分别反映了呼吸疾病的胸、肺阻力及顺应性,监测阻抗和电抗的变化有助于阐明疾病发生的机制、监控病情进展以及采取针对性的治疗。
发明内容
本发明的目的是,根据数据采集得到的原始呼吸振荡压力信号和流量信号,以较小的计算量得出呼吸阻力值。
为此给出本发明技术方案一:
医用呼吸监测仪呼吸阻力值计算方法,
对数据采集得到的原始呼吸振荡压力信号和流量信号在时域上进行移动平均滤波得到呼吸压力成分信号和流量成分信号,用原始呼吸振荡压力信号和流量信号减去呼吸压力成分信号和流量成分信号得到经肺组织吸收折射的振荡压力成分信号和流量成分信号;
在振荡成分信号数量为的一个时间窗内,所述的不小于2,记振荡成分信号的离散时间傅立叶正变换中的基波余弦分量的傅立叶系数为,基波正弦分量的傅立叶系数为,用近似作为该振荡成分信号的振幅,用近似作为该振荡成分信号的相位——以此方式近似计算出振荡成分压力信号在时刻的振幅和相位以及振荡成分流量信号在时刻的振幅和相位;
给出本发明优选技术方案二,其在技术方案一的基础上,对呼吸阻力信号进行移动平滑滤波。
给出本发明优选技术方案三,其在技术方案一或二的基础上,在求取呼吸阻力信号的过程中进行时间校正。
给出本发明优选技术方案五,其在技术方案一的基础上,所述时间窗窗宽等于一个振荡周期。
给出本发明优选技术方案六,其在技术方案一或五的基础上,所述的=20或21。
给出本发明优选技术方案七,其在技术方案二的基础上,移动平滑滤波所用移动时间窗窗宽等于两个振荡周期。
给出本发明优选技术方案八,其在技术方案一或七的基础上,移动平均滤波所用移动时间窗窗宽等于一个振荡周期。
本发明给出的方法,其中的全部或部分步骤可以通过建立功能模块构架,由计算机程序指令控制计算机***来完成。这些计算机程序指令存储在计算机可读存储介质中。
本发明的原理是:以最优线性近似法为理论基础,把呼吸看成一个线性***,肺部力学特征相对稳定,***的输入(振荡压力)与输出(振荡流量)具有因果关系,则两者的关系可利用传递函数(transfer function)表达为。数据采集得到的原始呼吸振荡信号(原始呼吸振荡信号包括压力信号和流量信号)包括振荡成分和呼吸成分,振荡成分经移动平均滤波的求和计算后为零,剩下呼吸成分(呼吸成分信号包括压力信号和流量信号)。用原始呼吸振荡信号减去呼吸成分信号就得到了经肺组织吸收折射的振荡成分信号(振荡成分信号包括压力信号和流量信号)。因为该振荡成分信号所含的众多谐波成分中,基波成分的能量是最大的,所以本发明直接取其基波分量的傅立叶系数来近似该振荡成分信号。本发明无需进行运算量巨大的傅立叶变换就可得出近似的呼吸阻力值,计算简单,实时性强。
附图说明
图1是R-I-E模型的实际粘性阻力与阻抗的一致性分析示意图。
具体实施方式
以下通过模型实验详述本发明。
在BiPAP通气模式下,振荡发生器在IPAP为8~24cmH2O,EPAP为4cmH2O的气道压力波动下,均能稳定地输出频率为1~7Hz,振幅为1~4cmH2O的正弦振荡,选5Hz作为振荡频率。振荡波和谐地叠加于气道压力和流量波形之上。***长时间工作稳定,计算获得的呼吸***Zrs、Rrs和Xrs时间曲线光滑规则,受干扰少,相位的变化与呼吸周期一致。
在BiPAP通气模式下,振荡发生器在IPAP为8~24cmH2O,EPAP为4cmH2O的气道压力波动下,均能稳定地输出频率为1~7Hz,振幅为1~4cmH2O的正弦振荡,振荡波和谐地叠加于气道压力和流量波形之上。
进行数据采集得到原始呼吸振荡压力信号和流量信号,包括振荡成分和呼吸成分,振荡频率为5Hz,采样频率为100Hz。对数据采集得到的原始呼吸振荡压力信号和流量信号在时域上进行移动平均滤波,得到呼吸压力成分信号和流量成分信号。移动平均滤波所用移动时间窗窗宽等于一个振荡周期。用原始呼吸振荡压力信号和流量信号减去呼吸压力成分信号和流量成分信号得到经肺组织吸收折射的振荡成分压力信号和流量信号。
今定义一个窗宽等于一个振荡周期的时间窗,因为振荡频率为5Hz,采样频率为100Hz,所以在该时间窗内有=20个振荡成分压力信号。在该时间窗内,记振荡成分压力信号的离散时间傅立叶正变换中的基波余弦分量的傅立叶系数为,基波正弦分量的傅立叶系数为,用近似作为振荡成分压力信号的振幅,用近似作为振荡成分压力信号的相位。据此,用任意时刻代入即可计算出在任意时刻的振荡成分压力信号的振幅和相位。
在同一个时间窗内,记振荡成分流量信号的离散时间傅立叶正变换中的基波余弦分量的傅立叶系数为,基波正弦分量的傅立叶系数为,用近似作为振荡成分流量信号的振幅,用近似作为振荡成分流量信号的相位。据此,用任意时刻代入即可计算出在任意时刻的振荡成分流量信号的振幅和相位。
然后对所得出的呼吸阻力信号(阻抗分量和/或电抗分量)进行移动平滑滤波,移动平滑滤波所用移动时间窗窗宽等于两个振荡周期,从而保留了低频成分,剔除了高频干扰。
在求取呼吸阻力信号的过程中进行时间校正。因为在计算得到的呼吸阻力时间曲线与原始流量或压力曲线进行比较时,时间上滞后,故给予时间校正更为合适。在移动平均滤波的时候,前后各十个点没有数据;在互相干法进行数据处理的时候,序列前后各十个点是没有数据;在进行移动平滑滤波的时候,所得数据前后各二十个点是没有数据的,以上都要进行数据的时间校正。其次,在计算呼吸阻力时,使用的数据点数不同也会造成相位偏移。上述实施例振荡频率为5Hz,采样频率为100Hz,则每个振荡波被采集数据点数是20个,为偶数,通过这20个点(一个完整的正弦波)计算出的呼吸阻力值对应于第10和第11数据点之间的时间段,从而造成相移,频率越高,相移越大。若用于计算的数据点数是奇数,就不会出现相移现象,例如将上例改为用21个数据点计算得到的呼吸阻力值将对应于第11数据点这一时刻。
相位校正后获得的模型阻抗R FOT (即所得出的呼吸阻力的阻抗分量)的均值为4.55±0.23cmH2O·s·L-1,模型电抗X FOT (即所得出的呼吸阻力的阻抗分量)为2.21±0.39cmH2O·s·L-1。R-I-E模型的实际粘性阻力R M 为3.29±0.27cmH2O·s·L-1,实际弹性阻力E M 为14.74±1.19cmH2O/L。故实际粘性阻力R M 与模型电抗R FOT 平均值之差为1.26±0.44cmH2O·s·L-1,一致性范围在0.36~2.14cmH2O·s·L-1,见图1。图1中,采用Bland-Altman图(Bland-Altman PLOT)实现:将同一个个体测量结果的均数作为横坐标,差值作为纵坐标。中间虚线为差值的均值,上下两条虚线为差值的95%上下置信限。然后看置信限是否在临床容许误差范围内(也就是,这样的误差大小在临床上能否接受),从而判断一致性好不好。95%上下置信限构造方法:对于差值正态分布的数据,采用平均相差±1.96*S,S是标准差;对于非正态数据,先正态转化,再用同样方法求解。各数据点来自R-I-E模型。RFOT为强迫振荡测定的模型阻抗,RM为模型的实际粘性阻力。图中实线为RFOT与RM的均差,虚线为2倍标准差。可见,采用本发明方法得到的结果比较接近于实际值,在临床上是可以接受的。
本文给出的方法,其中的全部或部分步骤可以通过建立功能模块构架,由计算机程序指令控制计算机***来完成。这些计算机程序指令存储在计算机可读存储介质中。
以上仅是本专利的较佳实施例,在此基础上的等同技术方案仍落入专利保护范围。
Claims (8)
1.医用呼吸监测仪呼吸阻力值计算装置,其特征是:
包括对数据采集得到的原始呼吸振荡压力信号和流量信号在时域上进行移动平均滤波得到呼吸压力成分信号和流量成分信号的装置,和用原始呼吸振荡压力信号和流量信号减去呼吸压力成分信号和流量成分信号得到经肺组织吸收折射的振荡压力成分信号和流量成分信号的装置;
在振荡成分信号 数量为的一个时间窗内,所述的不小于2,记振荡成分信号的离散时间傅立叶正变换中的基波余弦分量的傅立叶系数为,基波正弦分量的傅立叶系数为,包括以后述方式近似计算出振荡成分压力信号在时刻的振幅和相位以及振荡成分流量信号在时刻的振幅和相位的装置,所述的方式为用近似作为该振荡成分信号的振幅,用近似作为该振荡成分信号的相位;
2.根据权利要求1所述的医用呼吸监测仪呼吸阻力值计算装置,包括对呼吸阻力信号进行移动平滑滤波的装置。
3.根据权利要求1或2所述的医用呼吸监测仪呼吸阻力值计算装置,包括在求取呼吸阻力信号的过程中进行时间校正的装置。
5.根据权利要求1所述的医用呼吸监测仪呼吸阻力值计算装置,所述时间窗窗宽等于一个振荡周期。
6.根据权利要求1或5所述的医用呼吸监测仪呼吸阻力值计算装置,所述的=20或21。
7.根据权利要求2所述的医用呼吸监测仪呼吸阻力值计算装置,移动平滑滤波所用移动时间窗窗宽等于两个振荡周期。
8.根据权利要求1或7所述的医用呼吸监测仪呼吸阻力值计算装置,移动平均滤波所用移动时间窗窗宽等于一个振荡周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104087373A CN102488521A (zh) | 2011-12-09 | 2011-12-09 | 医用呼吸监测仪呼吸阻力值计算装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011104087373A CN102488521A (zh) | 2011-12-09 | 2011-12-09 | 医用呼吸监测仪呼吸阻力值计算装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102488521A true CN102488521A (zh) | 2012-06-13 |
Family
ID=46180432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011104087373A Pending CN102488521A (zh) | 2011-12-09 | 2011-12-09 | 医用呼吸监测仪呼吸阻力值计算装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102488521A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103300861A (zh) * | 2013-05-08 | 2013-09-18 | 深圳市科曼医疗设备有限公司 | 阻抗式呼吸测量*** |
CN115919255A (zh) * | 2022-11-22 | 2023-04-07 | 飞依诺科技股份有限公司 | 识别睡眠呼吸暂停低通气类型的方法、呼吸机及存储介质 |
CN116831558A (zh) * | 2023-06-30 | 2023-10-03 | 浙江柯洛德健康科技有限公司 | 一种基于强迫振荡的呼吸阻抗计算方法及计算装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4022193A (en) * | 1974-03-22 | 1977-05-10 | Siemens Aktiengesellschaft | Apparatus for determination of respiratory passageway resistance |
CN1972631A (zh) * | 2004-05-04 | 2007-05-30 | 达尔豪斯大学 | 评估气道高反应性中气道可变性的方法 |
CN101959455A (zh) * | 2008-03-10 | 2011-01-26 | 捷斯特公司 | 呼吸阻抗测量装置及其测量方法、呼吸阻抗显示方法 |
-
2011
- 2011-12-09 CN CN2011104087373A patent/CN102488521A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4022193A (en) * | 1974-03-22 | 1977-05-10 | Siemens Aktiengesellschaft | Apparatus for determination of respiratory passageway resistance |
CN1972631A (zh) * | 2004-05-04 | 2007-05-30 | 达尔豪斯大学 | 评估气道高反应性中气道可变性的方法 |
CN101959455A (zh) * | 2008-03-10 | 2011-01-26 | 捷斯特公司 | 呼吸阻抗测量装置及其测量方法、呼吸阻抗显示方法 |
Non-Patent Citations (2)
Title |
---|
王华,陈荣昌,何晟,罗群: "无创通气下应用强迫振荡技术检测呼吸阻力", 《中华生物医学工程杂志》 * |
王华: "《中国博士学位论文全文数据库 医药卫生科技辑》", 15 July 2010, 中国学术期刊(光盘版)电子杂志社 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103300861A (zh) * | 2013-05-08 | 2013-09-18 | 深圳市科曼医疗设备有限公司 | 阻抗式呼吸测量*** |
CN103300861B (zh) * | 2013-05-08 | 2015-04-08 | 深圳市科曼医疗设备有限公司 | 阻抗式呼吸测量*** |
CN115919255A (zh) * | 2022-11-22 | 2023-04-07 | 飞依诺科技股份有限公司 | 识别睡眠呼吸暂停低通气类型的方法、呼吸机及存储介质 |
CN116831558A (zh) * | 2023-06-30 | 2023-10-03 | 浙江柯洛德健康科技有限公司 | 一种基于强迫振荡的呼吸阻抗计算方法及计算装置 |
CN116831558B (zh) * | 2023-06-30 | 2024-03-29 | 浙江柯洛德健康科技有限公司 | 一种基于强迫振荡的呼吸阻抗计算方法及计算装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kaczka et al. | Oscillation mechanics of the respiratory system: applications to lung disease | |
Lutchen et al. | Optimal ventilation waveforms for estimating low-frequency respiratory impedance | |
Döllinger et al. | Dynamic vocal fold parameters with changing adduction in ex-vivo hemilarynx experiments | |
CN102488521A (zh) | 医用呼吸监测仪呼吸阻力值计算装置 | |
CN102551712A (zh) | 基于生物电阻抗的非侵入式排尿报警装置及监测方法 | |
Murugappan et al. | Acoustic characteristics of phonation in “wet voice” conditions | |
CN100386052C (zh) | 一种获得声门下压力值及计算发声效率的方法 | |
Lopes de Melo et al. | New impedance spectrometer for scientific and clinical studies of the respiratory system | |
CN105930665A (zh) | 一种基于超声衰减系数的肝脏脂肪定量测量方法 | |
Beaulieu et al. | Measurement of fractional order model parameters of respiratory mechanical impedance in total liquid ventilation | |
Chan et al. | Vocal fold vibration measurements using laser Doppler vibrometry | |
Alipour et al. | On pressure-frequency relations in the excised larynx | |
Semmler et al. | Fluid-structure-acoustic interactions in an ex vivo porcine phonation model | |
CN106840734B (zh) | 一种呼吸机跟随性能的评估方法、装置以及无创呼吸机 | |
CN106859649A (zh) | 呼吸率检测方法及装置 | |
CN106580324A (zh) | 一种呼吸信号提取方法及装置 | |
CN108994833A (zh) | 一种基于肌电活跃度反馈的关节助力控制方法 | |
Alipour et al. | Ventricular pressures in phonating excised larynges | |
CN203539352U (zh) | 振荡式呼吸功能监测仪 | |
Chen et al. | EMG oscillator model-based energy kernel method for characterizing muscle intrinsic property under isometric contraction | |
Alipour et al. | An experimental study of pulsatile flow in canine larynges | |
Hilberg et al. | Acoustic reflections during rhinometry: spatial resolution and sound loss | |
Bizzotto et al. | A portable fan-based device for evaluating lung function in horses by the forced oscillation technique | |
CN108742628A (zh) | 一种呼吸检测方法 | |
Peslin et al. | Partitioning of airway and respiratory tissue mechanical impedances by body plethysmography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20120613 |