CN110755077B - 一种基于呼吸装置的肺弹性系数测量方法及系统 - Google Patents

一种基于呼吸装置的肺弹性系数测量方法及系统 Download PDF

Info

Publication number
CN110755077B
CN110755077B CN201911127403.1A CN201911127403A CN110755077B CN 110755077 B CN110755077 B CN 110755077B CN 201911127403 A CN201911127403 A CN 201911127403A CN 110755077 B CN110755077 B CN 110755077B
Authority
CN
China
Prior art keywords
lung
coefficient
pipeline
pressure
volume
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.)
Active
Application number
CN201911127403.1A
Other languages
English (en)
Other versions
CN110755077A (zh
Inventor
张玉欣
金江春植
白晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihua University
Original Assignee
Beihua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihua University filed Critical Beihua University
Priority to CN201911127403.1A priority Critical patent/CN110755077B/zh
Publication of CN110755077A publication Critical patent/CN110755077A/zh
Application granted granted Critical
Publication of CN110755077B publication Critical patent/CN110755077B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/085—Measuring impedance of respiratory organs or lung elasticity
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4836—Diagnosis combined with treatment in closed-loop systems or 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/0003—Accessories therefor, e.g. sensors, vibrators, negative pressure
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2230/00—Measuring parameters of the user
    • A61M2230/40—Respiratory characteristics
    • A61M2230/46—Resistance or compliance of the lungs

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Pulmonology (AREA)
  • Pathology (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Physiology (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Hematology (AREA)
  • Emergency Medicine (AREA)
  • Operations Research (AREA)
  • Anesthesiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Computing Systems (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开一种基于呼吸装置的肺弹性系数测量方法及系统。肺弹性系数为肺内部压力与肺体积的商。该测量方法包括：获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；根据管道末端压力和肺内气体体积利用基于广义回归神经网络计算肺弹性变量系数；建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；将管道末端压力、肺内气体体积、管道气体流量和肺弹性变量系数代入呼吸方程，利用最小二乘法求呼吸方程的解，得到肺弹性常数系数；将肺弹性常数系数乘以肺弹性变量系数得到肺弹性系数。本发明能够适用于呼吸装置的实时调节和自动调节。

Description

一种基于呼吸装置的肺弹性系数测量方法及系统
技术领域
本发明涉及呼吸装置领域，特别是涉及一种基于呼吸装置的肺弹性系数测量方法及系统。
背景技术
当处于特殊环境下进行作业(例如水中营救溺水人员)或为临床病人供氧时通常需要配备呼吸装置。现有的呼吸装置大多具有自动调节通气量和/或通气压力的功能。
目前自动调节通气量的方式一般按照分钟通气量＝潮气量×通气频率来估算的。其中成人的潮气量根据体重按照6～8ml/kg估算，儿童的潮气量根据体重按照15～23ml/kg估算，成人通气频率为14～20次/分，儿童通气频率为18～40次/分。
自动调节通气压力的方式一般为为保证通气安全设置压力上下限。成人压力上限一般为50～60cmH2O，儿童压力上限一般为20～40cmH2O。
通气量和通气压力的双重调节方式一般为以固定压力通气，通气量不足时以固定容量的通气量来补充。该调节方式需要监测肺功能指标从而自动调节。但所监测的功能指标局限于潮气量的测定，而潮气量的测定的任何误差都会导致呼吸装置自动调控的失误。
对于呼吸装置的自动调节，比较理想的状态是根据人体的差异进行自适应调节。该差异体现在肺弹性系数等特性上。虽然肺弹性系数能够很好的适用于呼吸装置的自动调节，但肺弹性系数的测量需要暂时阻断人的呼吸气流，给人带来一定的痛苦，不适用于呼吸装置的实时调节和自动调节。
发明内容
本发明的目的是提供一种基于呼吸装置的肺弹性系数测量方法及系统，适用于呼吸装置的实时调节和自动调节。
为实现上述目的，本发明提供了如下方案：
一种基于呼吸装置的肺弹性系数测量方法，所述肺弹性系数为肺内部压力与肺体积的商；所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积，该测量方法包括：
获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；
根据所述管道末端压力和所述肺内气体体积利用基于广义回归神经网络计算所述肺弹性变量系数；
建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力。
可选的，所述根据所述管道末端压力和所述肺内气体体积利用基于广义回归神经网络计算所述肺弹性变量系数，具体包括：
将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
按照公式
Figure BDA0002277290530000021
计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子。
可选的，所述建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程，具体包括：
建立基本呼吸方程：
Figure BDA0002277290530000031
其中Pao(t)为管道末端压力，
Figure BDA0002277290530000032
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure BDA0002277290530000033
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值；
将所述基本呼吸方程进行关系转换，得到转换后的呼吸方程：
Figure BDA0002277290530000034
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure BDA0002277290530000035
为V(t)的一阶导数，F(t)为管道气体流量，
Figure BDA0002277290530000036
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数；
将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure BDA0002277290530000037
其中，
Figure BDA0002277290530000038
为已知量矩阵，
Figure BDA0002277290530000039
θg为待求量矩阵，
Figure BDA00022772905300000310
可选的，所述将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数，具体包括：
对所述矩阵表达式的两侧进行积分运算，得到积分表达式；
将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
一种基于呼吸装置的肺弹性系数测量系统，所述肺弹性系数为肺内部压力与肺体积的商；所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积，该测量系统包括：
获取模块，用于获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；
肺弹性变量系数计算模块，用于根据所述管道末端压力和所述肺内气体体积利用基于广义回归神经网络计算所述肺弹性变量系数；
呼吸方程建立模块，用于建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
最小二乘求解模块，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
肺弹性系数确定模块，用于将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力。
可选的，所述肺弹性变量系数计算模块包括：
中心点选取单元，用于将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
肺弹性变量系数计算单元，用于按照公式
Figure BDA0002277290530000041
计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子。
可选的，所述呼吸方程建立模块包括：
基本呼吸方程建立单元，用于建立基本呼吸方程：
Figure BDA0002277290530000051
其中Pao(t)为管道末端压力，
Figure BDA0002277290530000052
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure BDA0002277290530000053
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值；
关系转换单元，用于将所述基本呼吸方程进行关系转换，得到转换后的呼吸方程：
Figure BDA0002277290530000054
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure BDA0002277290530000055
为V(t)的一阶导数，F(t)为管道气体流量，
Figure BDA0002277290530000056
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数；
矩阵转换单元，用于将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure BDA0002277290530000057
其中，
Figure BDA0002277290530000058
为已知量矩阵，
Figure BDA0002277290530000059
θg为待求量矩阵，
Figure BDA00022772905300000510
可选的，所述最小二乘求解模块包括：
积分单元，用于对所述矩阵表达式的两侧进行积分运算，得到积分表达式；
最小二乘求解单元，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明的基于呼吸装置的肺弹性系数测量方法及系统，在获取到呼吸装置直接采集到的数据的基础上利用广义回归神经网络计算肺弹性变量系数，在得到肺弹性变量系数的基础上利用呼吸装置直接采集到的数据对呼吸方程进行求解得到肺弹性常数系数，从而实现肺弹性系数的测量。本发明依托广义回归神经网络计算肺弹性变量系数，从而直接得到肺弹性变量系数，并为肺弹性常数系数的计算建立基础，实现了直接利用呼吸装置采集的数据进行肺弹性系数的测量，避免阻断人的呼吸气流，适用于呼吸装置的实时调节和自动调节。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例1的基于呼吸装置的肺弹性系数测量方法的方法流程图；
图2为婴儿使用呼吸装置时呼吸装置采集到的管道末端压力曲线图；
图3为婴儿使用呼吸装置时呼吸装置采集到的肺内气体体积曲线图；
图4为婴儿使用呼吸装置时呼吸装置采集到的管道气体流量曲线图；
图5为婴儿吸气时管道末端压力与肺内气体体积的关系曲线图；
图6为婴儿呼气时管道末端压力与肺内气体体积的关系曲线图；
图7为吸气时肺内部静态压力的真实值与本发明的估计值的对比结果图；
图8为呼气时肺内部静态压力的真实值与本发明的估计值的对比结果图；
图9为本发明实施例2的基于呼吸装置的肺弹性系数测量系统的系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1：
该实施例1提供一种基于呼吸装置的肺弹性系数测量方法，所述肺弹性系数为肺内部压力与肺体积的商；所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积。
图1为本发明实施例1的基于呼吸装置的肺弹性系数测量方法的方法流程图。
参见图1，该测量方法包括：
步骤101：获取呼吸装置检测到的一个呼吸周期内的管道末端压力(管道末端压力表示的是用户气道压力)Pao(t)、肺内气体体积V(t)和管道气体流量F(t)。
以处于呼吸窘迫状态的某个婴儿为例。采样间隔T＝0.005秒，即采样频率200Hz，采样点数为N0＝229，其中吸气点数N1是116，呼气点数N2是113。
图2为婴儿使用呼吸装置时呼吸装置采集到的管道末端压力曲线图。
图3为婴儿使用呼吸装置时呼吸装置采集到的肺内气体体积曲线图。
图4为婴儿使用呼吸装置时呼吸装置采集到的管道气体流量曲线图。
图5为婴儿吸气时管道末端压力与肺内气体体积的关系曲线图。
图6为婴儿呼气时管道末端压力与肺内气体体积的关系曲线图。
获取的样本数据参见图2～图6。
步骤102：根据所述管道末端压力和所述肺内气体体积利用基于广义回归神经网络计算所述肺弹性变量系数。
广义回归神经网络(General Regression Neural Network，GRNN)是径向基网络(Radial basis function network，RBF)的变形形式，具有良好的非线性逼近能力，训练更为方便。
该步骤102具体包括：
1、将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
2、按照公式
Figure BDA0002277290530000081
(通过将fg(V)定义为GRNN网络输出函数形式得到)计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子。
根据呼吸样本值分布，设定GRNN网络隐含层神经元个数N＝7。Pi和Vi分别是隐含层第i个中心点对应的呼吸样本的Pao(t)值和V(t)值，每一个隐含层神经元即为一个数据中心点，在吸气样本中选择7个数据点做为数据中心点，对全部吸气样本进行GRNN网络输出值计算；在呼气样本中选择7个数据点做为数据中心点，对全部呼气样本进行GRNN网络输出值计算。
隐含层数据中心点的选择，中心点的选择要尽可能平均分布在整个呼吸周期中，不能过分集中。
步骤103：建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
该步骤103，具体包括：
1、建立基本呼吸方程：
Figure BDA0002277290530000082
其中Pao(t)为管道末端压力，
Figure BDA0002277290530000083
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure BDA0002277290530000084
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值，包含方程误差和测量噪声；a和b为常系数。
2、将所述基本呼吸方程进行关系转换
肺内部压力表达式Pl(t)＝cfg(V)V(t)
其中cfg(V)是肺弹性系数。肺内部压力表达式描述的肺内部压力Pl(t)与肺内气体体积V的关系，就是肺的静态特性。c是待求的肺弹性常数系数。
呼吸装置管道压力损耗表达式Pr(t)＝fr(F)F(t)
其中fr(F)是气道阻力。
气体流量和气体体积的关系式
Figure BDA0002277290530000091
所以
Figure BDA0002277290530000092
如果考虑到气体流量的二阶作用，
Figure BDA0002277290530000093
此时
Figure BDA0002277290530000094
将Pl(t)、Pr(t)代人到呼吸系统方程，可得：
Figure BDA0002277290530000095
利用气体流量和气体体积的关系式
Figure BDA0002277290530000096
对上式进行改写，可得转换后的呼吸方程：
Figure BDA0002277290530000097
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure BDA0002277290530000098
为V(t)的一阶导数，F(t)为管道气体流量，
Figure BDA0002277290530000099
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数。
3、将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure BDA00022772905300000910
其中，
Figure BDA00022772905300000911
为已知量矩阵，
Figure BDA00022772905300000912
θg为待求量矩阵，
Figure BDA00022772905300000913
步骤104：将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
该步骤104具体包括：
1、对所述矩阵表达式的两侧进行积分运算，得到积分表达式。
矩阵表达式是连续时间函数模型，而呼吸样本Pao(t)、V(t)、F(t)均为等时间间隔T的采样数据，因此，对矩阵表达式采用数值积分运算的方法进行处理。对矩阵表达式左右两端进行[kT,(k-l)T]的积分，其中k＝1,2,3…,N0，N0是呼吸样本个数，l＝20。
矩阵表达式左侧进行数值积分运算结果定义为p(k)，则：
Figure BDA0002277290530000101
其中gj是利用复合梯形公式法进行数值积分时的系数
Figure BDA0002277290530000102
定义
Figure BDA0002277290530000103
对应的数值积分运算结果为φg(k)，定义e(t)对应的数值积分运算结果为e(k)，则：
Figure BDA0002277290530000104
Figure BDA0002277290530000105
呼吸方程的数值积分运算方程为：
Figure BDA0002277290530000106
2、将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
将呼吸样本数据Pao(t)代入到积分公式p(k)可以求出p(k)，将呼吸样本数据V(t)、F(t)代入积分公式φg(k)可以求出φg(k)。k的取值范围是[l+1,N0]，如果l＝20，则k的取值范围是[21,N0]。
设由p(k)，k＝l+1,…,N0，组成的矩阵的转置矩阵为P，即P＝[p(l+1),p(l+2),…,p(N0)]T；φg(k)，k＝l+1,…,N0，组成的矩阵的转置矩阵为Φ，即Φ＝[φg(l+1),φg(l+2),…,φg(N0)]T；e(k)，k＝l+1,…,N0，组成的矩阵的转置矩阵为e，即e＝[e(l+1),e(l+1),…，e(N0)]T。
因此呼吸方程的数值积分运算方程可以写为
p＝Φθg+e
定义
Figure BDA0002277290530000111
是利用最小二乘法求出的θg的估计值，
Figure BDA0002277290530000112
矩阵中的参数为θg矩阵参数的估计值，定义为
Figure BDA0002277290530000113
Figure BDA0002277290530000114
从而计算得到矩阵θg中的各个参数。
定义最小二乘法估计误差是J，则：
J＝(P-Φθg)T(P-Φθg)
步骤105：将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力。
将计算出来的参数
Figure BDA0002277290530000115
和步骤102出来的fg(V)相乘得到肺弹性系数估计值
Figure BDA0002277290530000116
当需要调节通气量或通气压力时，可以基于肺弹性系数的计算结果进行进一步计算：
将实时测量的V(t)值乘以肺弹性系数估计值求出此时肺内部静态压力Pl(t)的估计值
Figure BDA0002277290530000117
将矩阵
Figure BDA0002277290530000118
和步骤102计算出来的fg(V)代入呼吸方程，可以得到符合用户呼吸特性的呼吸方程。
根据参数
Figure BDA0002277290530000119
以及实时测量的V(t)求出时呼吸装置管道丢失压力Pr(t)的估计值。
下面通过静态数据拟合实验对本实施例1的方案进行详细说明：
针对步骤101中的呼吸数据样本，利用本发明的方案计算吸气样本的肺内部静态压力Pl(t)估算值，从而构建图7中曲线；计算呼气样本的肺内部静态压力Pl(t)估算值，从而构建图8中曲线。吸气样本的最小二乘法估计误差是J＝0.0038；呼气样本的最小二乘法估计误差是J＝0.0025。
图7、图8中矩形数据点，是利用特殊的实验方法测出的病人在吸气、呼气过程中肺内部Pl(t)的真实值，即肺静态值。
从图7和图8可以看出，本方案计算的Pl(t)值曲线对肺静态值拟合较好。
肺呼吸时和吸气时的静态特性并不完全相同，图7、图8的仿真实验结果显示的Pl(t)值曲线不是完全重合的曲线，也印证了这个特点。
误差分析：
①对图7中吸气样本的Pl(t)估算值进行数据拟合。由于吸气样本估算值共116个，数据点过多，全部用来进行曲线拟合计算量过大，且准确性降低。采用分段数据点参与的方式进行曲线拟合。以图7中的8个静态数据的Pao(t)值为中心，分别选择6个Pl(t)估算值点进行曲线拟合。
以图8中静态点A为例，A点实测的Pl＝7.4，V＝3.6。在116个吸气样本的Pl(t)估算值中选择第25到第30个数据点进行拟合，这6个点的Pl(t)值分别为[6.2420，6.7217，7.1942，7.6519，8.0867，8.4980]，V值分别为[3.2124,3.4127,3.6100,3.8011,3.9827,4.1544]。
采用二阶多项式进行曲线拟合，即V＝a0+a1Pl+a2Pl 2，代人这6个点数据得到的拟合多项式为V＝0.609+0.4167Pl+0.0001Pl 2。当Pl＝7.4时，拟合得出的V＝3.6959。与A点的V实测值之间误差为0.0959。
采用同样的方法，计算出的图7中8个吸气静态数据点的误差值如表1所示，平均绝对值误差为0.1693；图8中的8个呼气静态数据点的误差值如表2所示，平均绝对值误差为：0.3728。误差值符合实际应用需求。
综上，依据本发明的方法，对肺弹性量值估算得出肺内压力的方法，误差值较小。
表1图7中8个吸气静态数据点的误差值
Figure BDA0002277290530000121
表2图8中的8个呼气静态数据点的误差值
Figure BDA0002277290530000131
实施例2：
该实施例2提供一种基于呼吸装置的肺弹性系数测量系统。所述肺弹性系数为肺内部压力与肺体积的商；所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积。
图9为本发明实施例2的基于呼吸装置的肺弹性系数测量系统的系统结构图。
参见图9，该测量系统包括：
获取模块201，用于获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；
肺弹性变量系数计算模块202，用于根据所述管道末端压力和所述肺内气体体积利用基于广义回归神经网络计算所述肺弹性变量系数；
呼吸方程建立模块203，用于建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
最小二乘求解模块204，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
肺弹性系数确定模块205，用于将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力。
可选的，所述肺弹性变量系数计算模块202包括：
中心点选取单元，用于将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
肺弹性变量系数计算单元，用于按照公式
Figure BDA0002277290530000141
计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子。
可选的，所述呼吸方程建立模块203包括：
基本呼吸方程建立单元，用于建立基本呼吸方程：
Figure BDA0002277290530000142
其中Pao(t)为管道末端压力，
Figure BDA0002277290530000143
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure BDA0002277290530000144
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值；
关系转换单元，用于将所述基本呼吸方程进行关系转换，得到转换后的呼吸方程：
Figure BDA0002277290530000145
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure BDA0002277290530000146
为V(t)的一阶导数，F(t)为管道气体流量，
Figure BDA0002277290530000147
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数；
矩阵转换单元，用于将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure BDA0002277290530000148
其中，
Figure BDA0002277290530000149
为已知量矩阵，
Figure BDA00022772905300001410
θg为待求量矩阵，
Figure BDA00022772905300001411
可选的，所述最小二乘求解模块204包括：
积分单元，用于对所述矩阵表达式的两侧进行积分运算，得到积分表达式；
最小二乘求解单元，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明的基于呼吸装置的肺弹性系数测量方法及系统，在获取到呼吸装置直接采集到的数据的基础上利用广义回归神经网络计算肺弹性变量系数，在得到肺弹性变量系数的基础上利用呼吸装置直接采集到的数据对呼吸方程进行求解得到肺弹性常数系数，从而实现肺弹性系数的测量。本发明依托广义回归神经网络计算肺弹性变量系数，从而直接得到肺弹性变量系数，并为肺弹性常数系数的计算建立基础，实现了直接利用呼吸装置采集的数据进行肺弹性系数的测量，避免阻断人的呼吸气流，适用于呼吸装置的实时调节和自动调节。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (6)

1.一种基于呼吸装置的肺弹性系数测量方法，所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积，其特征在于，该测量方法包括：
获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；
根据所述管道末端压力Pao(t)和所述肺内气体体积V(t)利用广义回归神经网络GRNN计算所述肺弹性变量系数；
建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力；
所述根据所述管道末端压力Pao(t)和所述肺内气体体积V(t)利用广义回归神经网络GRNN计算所述肺弹性变量系数，具体包括：
将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
按照公式
Figure FDA0003458606740000011
计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子；
根据呼吸样本值分布，设定GRNN网络隐含层神经元个数N＝7；Pi和Vi分别是隐含层第i个中心点对应的呼吸样本的Pao(t)值和V(t)值，每一个隐含层神经元即为一个数据中心点，在吸气样本中选择7个数据点做为数据中心点，对全部吸气样本进行GRNN网络输出值计算；在呼气样本中选择7个数据点做为数据中心点，对全部呼气样本进行GRNN网络输出值计算；隐含层数据中心点的选择，中心点的选择要尽可能平均分布在整个呼吸周期中，不能过分集中。
2.根据权利要求1所述的基于呼吸装置的肺弹性系数测量方法，其特征在于，所述建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程，具体包括：
建立基本呼吸方程：
Figure FDA0003458606740000021
其中Pao(t)为管道末端压力，
Figure FDA0003458606740000022
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure FDA0003458606740000023
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值；
将所述基本呼吸方程进行关系转换，得到转换后的呼吸方程：
Figure FDA0003458606740000024
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure FDA0003458606740000025
为V(t)的一阶导数，F(t)为管道气体流量，
Figure FDA0003458606740000026
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数；
将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure FDA0003458606740000027
其中，
Figure FDA0003458606740000028
为已知量矩阵，
Figure FDA0003458606740000029
θg为待求量矩阵，
Figure FDA00034586067400000210
3.根据权利要求2所述的基于呼吸装置的肺弹性系数测量方法，其特征在于，所述将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数，具体包括：
对所述矩阵表达式的两侧进行积分运算，得到积分表达式；
将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
4.一种基于呼吸装置的肺弹性系数测量系统，所述肺弹性系数为肺弹性常数系数和肺弹性变量系数相乘的积，其特征在于，该测量系统包括：
获取模块，用于获取呼吸装置检测到的一个呼吸周期内的管道末端压力、肺内气体体积和管道气体流量；
肺弹性变量系数计算模块，用于根据所述管道末端压力Pao(t)和所述肺内气体体积V(t)利用广义回归神经网络GRNN计算所述肺弹性变量系数；
呼吸方程建立模块，用于建立表达管道末端压力、肺内气体体积、管道气体流量、肺弹性常数系数和肺弹性变量系数之间关系的呼吸方程；
最小二乘求解模块，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述呼吸方程，利用最小二乘法求所述呼吸方程的解，得到所述肺弹性常数系数；
肺弹性系数确定模块，用于将所述肺弹性常数系数乘以所述肺弹性变量系数得到所述肺弹性系数；所述肺弹性系数用于调节所述呼吸装置的通气量和/或通气压力；
所述肺弹性变量系数计算模块包括：
中心点选取单元，用于将一个呼吸周期内的管道末端压力和肺内气体体积数据作为样本数据，从呼气阶段的样本数据中选取N个数据中心点，从吸气阶段的样本数据中选取N个数据中心点；
肺弹性变量系数计算单元，用于按照公式
Figure FDA0003458606740000041
计算所述肺弹性变量系数，其中fg(V)为肺弹性变量系数，i表示各数据中心点，Pi为第i个数据中心点的管道末端压力，V为除数据中心点以外的其他样本数据的肺内气体体积，Vi为第i个数据中心点的肺内气体体积，σ为光滑因子；
根据呼吸样本值分布，设定GRNN网络隐含层神经元个数N＝7；Pi和Vi分别是隐含层第i个中心点对应的呼吸样本的Pao(t)值和V(t)值，每一个隐含层神经元即为一个数据中心点，在吸气样本中选择7个数据点做为数据中心点，对全部吸气样本进行GRNN网络输出值计算；在呼气样本中选择7个数据点做为数据中心点，对全部呼气样本进行GRNN网络输出值计算；隐含层数据中心点的选择，中心点的选择要尽可能平均分布在整个呼吸周期中，不能过分集中。
5.根据权利要求4所述的基于呼吸装置的肺弹性系数测量系统，其特征在于，所述呼吸方程建立模块包括：
基本呼吸方程建立单元，用于建立基本呼吸方程：
Figure FDA0003458606740000042
其中Pao(t)为管道末端压力，
Figure FDA0003458606740000043
为Pao(t)的一阶导数，Pl(t)为肺内部压力，Pr(t)为呼吸装置管道压力损耗，V(t)为肺内气体体积，
Figure FDA0003458606740000044
为V(t)的二阶导数，Peea为呼吸末端肺泡压力，e(t)为误差值；
关系转换单元，用于将所述基本呼吸方程进行关系转换，得到转换后的呼吸方程：
Figure FDA0003458606740000045
其中cfg构成所述肺弹性系数，c为肺弹性常数系数，fg为肺弹性变量系数，
Figure FDA0003458606740000046
为V(t)的一阶导数，F(t)为管道气体流量，
Figure FDA0003458606740000047
为F(t)的一阶导数，a、c、fg、r1、r2、b均为系数；
矩阵转换单元，用于将所述转换后的呼吸方程转换为矩阵形式，得到矩阵表达式：
Figure FDA0003458606740000051
其中，
Figure FDA0003458606740000052
为已知量矩阵，
Figure FDA0003458606740000053
θg为待求量矩阵，
Figure FDA0003458606740000054
6.根据权利要求5所述的基于呼吸装置的肺弹性系数测量系统，其特征在于，所述最小二乘求解模块包括：
积分单元，用于对所述矩阵表达式的两侧进行积分运算，得到积分表达式；
最小二乘求解单元，用于将所述管道末端压力、所述肺内气体体积、所述管道气体流量和所述肺弹性变量系数代入所述积分表达式，利用最小二乘法求解，得到所述待求量矩阵内各元素的值。
CN201911127403.1A 2019-11-18 2019-11-18 一种基于呼吸装置的肺弹性系数测量方法及系统 Active CN110755077B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911127403.1A CN110755077B (zh) 2019-11-18 2019-11-18 一种基于呼吸装置的肺弹性系数测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911127403.1A CN110755077B (zh) 2019-11-18 2019-11-18 一种基于呼吸装置的肺弹性系数测量方法及系统

Publications (2)

Publication Number Publication Date
CN110755077A CN110755077A (zh) 2020-02-07
CN110755077B true CN110755077B (zh) 2022-02-22

Family

ID=69338103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911127403.1A Active CN110755077B (zh) 2019-11-18 2019-11-18 一种基于呼吸装置的肺弹性系数测量方法及系统

Country Status (1)

Country Link
CN (1) CN110755077B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115297916A (zh) * 2020-03-23 2022-11-04 深圳迈瑞生物医疗电子股份有限公司 一种呼吸监测装置及方法
CN111973188A (zh) * 2020-08-26 2020-11-24 北京航空航天大学 一种基于神经网络估算呼吸力学参数的方法
CN113769217B (zh) * 2021-09-17 2023-07-07 北华大学 一种呼吸系统潮气量控制方法、装置及设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6068602A (en) * 1997-09-26 2000-05-30 Ohmeda Inc. Method and apparatus for determining airway resistance and lung compliance
CN104050517A (zh) * 2014-06-27 2014-09-17 哈尔滨工业大学 基于grnn神经网络的光伏发电预测方法
CN105899249A (zh) * 2013-10-29 2016-08-24 通用电气公司 用于确定患者特定的通气需求的方法和布置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109803708B (zh) * 2016-10-07 2022-03-01 皇家飞利浦有限公司 呼吸设备和通气机设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6068602A (en) * 1997-09-26 2000-05-30 Ohmeda Inc. Method and apparatus for determining airway resistance and lung compliance
CN105899249A (zh) * 2013-10-29 2016-08-24 通用电气公司 用于确定患者特定的通气需求的方法和布置
CN104050517A (zh) * 2014-06-27 2014-09-17 哈尔滨工业大学 基于grnn神经网络的光伏发电预测方法

Also Published As

Publication number Publication date
CN110755077A (zh) 2020-02-07

Similar Documents

Publication Publication Date Title
CN110755077B (zh) 一种基于呼吸装置的肺弹性系数测量方法及系统
CN102770070B (zh) 用于实时评估肺力学的系统和方法
US7460959B2 (en) System and method for estimating oxygen concentration in a mixed gas experiencing pressure fluctuations
US20110146683A1 (en) Sensor Model
CN109803708A (zh) 使用压力控制呼吸来估计肺顺应性和肺阻力以允许所有呼吸肌反冲生成的压力消失
JP2005537068A (ja) 呼吸仕事量を予測するための方法および装置
JP2012503527A (ja) 医療用人工呼吸器における患者呼吸努力力学のモデル予測オンライン同定
WO2011027243A1 (en) System and method for quantifying lung compliance in a self-ventilating subject
WO2011145014A1 (en) System and method for estimating upper airway resistance and lung compliance employing induced central apneas
CN111728298A (zh) 一种高同步供氧呼吸装置、呼吸监测系统及方法
CN111565780A (zh) 麻醉机及麻药输出浓度监测方法、系统、设备、存储介质
EP3397330B1 (en) Device and method for controlling ventilatory assist
CN107690310B (zh) 用于经由连续参数估计来监测患者呼吸状态的无创方法
Stout et al. Pulmonary blood flow determined by continuous analysis of pulmonary N2O exchange
CN112754465A (zh) 一种压力控制机械通气下肺部准静态顺应性估测方法
WO2020054360A1 (ja) 呼吸パラメータ推定装置および呼吸パラメータ推定方法
Kanae et al. Modeling of respiration and estimation of pulmonary elastance
Miyamoto et al. Online computer for assessing respiratory and metabolic function during exercise
CN111818848A (zh) 经由气流阻滞器调节的呼气气流受限检测
US20230191067A1 (en) Expiratory filter with embedded detectors
CN215653198U (zh) 球囊复苏器
CN213046880U (zh) 一种新型手持式鼻阻力测试装置
US20220167870A1 (en) Device for measuring the concentration of gases in exhaled air
CN115844371A (zh) 基于肺内气体空间分布变化的肺功能参数表征方法
CN116763289A (zh) 基于单一呼出气体的能量代谢评估方法及系统

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
GR01 Patent grant
GR01 Patent grant