CN116796656A - Method for estimating lumped parameter respiratory system model parameters - Google Patents
Method for estimating lumped parameter respiratory system model parameters Download PDFInfo
- Publication number
- CN116796656A CN116796656A CN202310553702.1A CN202310553702A CN116796656A CN 116796656 A CN116796656 A CN 116796656A CN 202310553702 A CN202310553702 A CN 202310553702A CN 116796656 A CN116796656 A CN 116796656A
- Authority
- CN
- China
- Prior art keywords
- respiratory system
- system model
- lumped parameter
- estimating
- derivative
- 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.)
- Granted
Links
- 210000002345 respiratory system Anatomy 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000000241 respiratory effect Effects 0.000 abstract description 6
- 238000009423 ventilation Methods 0.000 abstract description 2
- 230000004069 differentiation Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000007429 general method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000005399 mechanical ventilation Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
The application provides a method for estimating parameters of a lumped parameter respiratory system model, which is used for estimating the parameters of the lumped parameter respiratory system model, and based on the real-time acquisition of the pressure value, the flow velocity value and the accumulated gas quantity of a respiratory circuit in any ventilation mode, each parameter of the lumped parameter respiratory system model in any order can be estimated efficiently and simply according to the method provided by the application.
Description
Technical Field
The present application relates to respiratory system models, and more particularly to a method for estimating lumped parameter respiratory system model parameters.
Background
In respiratory mechanical ventilation research, a non-lumped parameter respiratory system model can be asymptotically utilized by a finite-order lumped parameter model, and parameter estimation of the lumped parameter respiratory system model is an important link for researching respiratory system characteristics and applying medical intelligent mechanical ventilation. At present, a general method for efficiently and simply estimating various constant coefficient parameters of any lumped parameter respiratory system model is lacked.
Disclosure of Invention
In view of the above problems, the present application aims to provide a method for estimating parameters of a lumped parameter respiratory system model, which is used for estimating the parameters of the lumped parameter respiratory system model, and based on real-time acquisition of a pressure value, a flow velocity value and an accumulated gas volume of a respiratory circuit in any ventilation mode, each parameter of an arbitrary order lumped parameter respiratory system model can be efficiently and simply estimated according to the method provided by the present application.
The method of the present application for estimating lumped parameter respiratory system model parameters, wherein,
the lumped parameter respiratory system model is expressed as:
P=a 0 V+a 1 V′+a 2 V″+…+a k V (k) +P 0 ;
wherein P is air pressure, V is the 1 st derivative of V' in the input total air formula, and represents the volume flow rate; v' is the 2 nd derivative of V and represents the airflow acceleration; v (V) (k) K-th derivative of V;
P 0 is the initial pressure in the respiratory system;
a 0 、a 1 、a 2 、…、a k coefficient parameters, a, of a lumped parameter respiratory system model 0 Is elastic; a, a 1 Is airway resistance; a, a 2 Is the inertial resistance coefficient;
using matricesObtaining fitting results of various coefficients;
wherein ,m is the number of terms of the lumped parameter respiratory system model;
the matrix X is an n-row m-column matrix;
is an n-dimensional barometric pressure vector, which is composed of P values of n time sampling points of a breathing circuit.
Preferably, the numerical initial value of the derivative of V above 2 nd order is set to 0, and the k-th order derivative term value at the other timing t is calculated as follows:
where h is the time step.
Preferably, when the lumped parameter respiratory system model is a first order linear model, the lumped parameter respiratory system model is expressed as:
P=a 0 V+a 1 V′+P 0 ;
wherein ,
the method for estimating the parameters of the lumped parameter respiratory system model can be suitable for the lumped parameter respiratory system model of each order, is a general method, and can be used for conveniently and rapidly estimating the coefficient parameters.
Drawings
FIG. 1 is a comparison of the parameter estimation of the lumped parameter respiratory system model of equation 1 with the measured data respiration volume by the method of estimating the lumped parameter respiratory system model parameters of the present application;
FIG. 2 is a comparison of the parameter estimation of the lumped parameter respiratory system model of equation 2 with the measured data respiration volume by the method of estimating the lumped parameter respiratory system model parameters of the present application;
Detailed Description
The method for estimating the parameters of the lumped parameter respiratory system model according to the present application will be described in detail with reference to the accompanying drawings.
1. Presetting the number m of lumped parameter respiratory system model items to be fitted and a corresponding form.
Taking m=5 as an example, i.e. there are 5 terms, the lumped parameter respiratory system model has the following form:
P=a 0 V+a 1 V′+a 2 (V) 2 +a 3 (V′) 2 +P 0
wherein the relation between P air pressure and total input air volume V in the respiratory system is described, wherein V 'is the 1 st derivative of V, namely the air flow velocity, V' is the 2 nd derivative, namely the air flow acceleration, P 0 Is the initial pressure within the respiratory system. At this time a 0 I.e. the elasticity (or reciprocal of the compliance) in a typical breathing model, and a 1 Resistance to the airway, a 2 Is the inertial resistance coefficient.
For higher order system cases, the lumped parameter respiratory system model has the following form:
P=a 0 V+a 1 V′+a 2 V″+…+a k V (k) +P 0
in the formula V(k) Representing the k-th derivative.
The preset model can contain a i (V′) 2 、a j e V″ The principle of the constant coefficient nonlinear term estimation method is that the constant coefficient nonlinear term is obtained based on the fact that the coefficient partial derivative of each term is 0, and the fact that the corresponding term numerical sequence in the input matrix in the next step corresponds to the preset model term one by one is guaranteed.
2. And collecting the pressure value, the flow velocity value and the accumulated gas value of the breathing circuit in real time, and sorting the pressure value, the flow velocity value and the accumulated gas value into corresponding matrix form data according to the number m and the corresponding form of a preset model term and a numerical differentiation method to be used as algorithm input.
If the value of the acquired breathing circuit is n time points, the input data are n-dimensional barometric pressure vector P and n rows and m columns of data matrix X:
the medium-high order differential value can be obtained by a numerical differential method with proper precision, such as a three-point formula based on interpolation, a Simpson numerical differential formula or an extrapolation method. The numerical differentiation method of which accuracy is selected is not the key point of the algorithm, optionally, at the beginning and end of the respiratory cycle (differentiation end point), the numerical initial value of the higher-order differentiation (2 nd order derivative and above) in the above formula is optionally set to 0, and the numerical value of the kth order differentiation term at other time sequences t is as follows (midpoint formula):
where h is the time step.
If the flow rate value and the accumulated gas value are limited by the sampling condition, only one of the flow rate value and the accumulated gas value, the other missing value can be obtained by a proper value differentiation and value integration method.
3. And calculating according to the following matrix to obtain fitting results of various coefficients.
Wherein each coefficient vector a is as follows:
in the above formula, the m-dimensional A vector corresponds to each coefficient to be estimated and fitted, and only involves 2 matrix multiplications, 1 inversion and 1 matrix vector multiplication.
The algorithm method for estimating and fitting each constant coefficient parameter by the lumped parameter respiratory system model is derived based on the extreme point with each bias of 0.
In algorithm derivation, for a specified fitting target, for example, the sum of squares SSR (sum of squared residuals) of residuals or other similar important fitting indexes can be selected, when the fitting target takes the minimum value, the partial derivatives of all coefficients in the lumped model are 0, and at this time, all coefficients to be fitted can be solved by a series of 1-time equations. The nonlinear term pair coefficient bias is still the 1 st order term of the constant coefficient to be estimated, so the parameter estimation method is the same. For clearer presentation, the calculation mode of each coefficient can be further written as a matrix solving mode.
Solving each order coefficient by taking extremum bias as 0 for fitting target is characterized by the method, the change of the form of the arithmetic writing method still belongs to the category of the method, for example, the matrix solving form is disassembled into the following estimation mode of each order coefficient (1-order linear model P=a 0 V+a 1 V′+P 0 Equivalent calculations are examples):
in the process of preparing input data, the selection of numerical differentiation and numerical integration methods with different accuracies is not the core of the method, and the whole framework of the same type of method after replacement and change still belongs to the category of the method.
Next, the lumped parameter respiratory system models of formulas 1 and 2 are estimated by the method of the present application, respectively, and compared with the measured respiratory volume.
P=a 0 V+a 1 V′+a 2 V″+a 3 ln(V+1)+P 0 (1)
The parameter estimation is based on animal experimental measured data (obtained by actual measurement of a Delge Savina 300 ventilator), and data of pressure P (unit mbar), flow velocity V' (unit L/s) and respiration volume V (unit L) of each 10 milliseconds in time sequence are recorded as one line (data see data. Txt, 150 seconds data are intercepted).
For the nonlinear equation of formula 1, V "and the corresponding value obtained using the numerical differentiation method, ln (v+1) is substituted into the value of V, and is sorted into the following numerical matrix form in each order as input.
Then, p0= 3.47632212, a0= 57.29251258, a1= 7.17413342, a2= 4.47134159 and a3= -36.22221325 in the formula 1 are calculated by using the method proposed by the patent. Corresponding to the model represented by equation 1, the correlation coefficient 0.7891 between the calculated breathing amount (dotted line) and the measured breathing amount (solid line) is substituted into the parameter calculation, as shown in fig. 1.
P=a 0 V+a 1 V′+a 2 V″+P 0 (2)
For equation 2, the matrix input data matrix can be obtained based on the measured data, and p0= 3.39674027, a0= 26.11755376, a1= 7.17137465, a2= 0.08733642 can be obtained using the calculation method of this patent. The correlation coefficient between the calculated breathing amount (dotted line) and the measured breathing amount (solid line) of the corresponding model (2) is 0.9970 (the preset model 2 in the example is more consistent with the measured data characteristic than the model 1), as shown in fig. 2.
By using the method of the application, various constant coefficient parameters of any lumped respiratory system model can be estimated, the calculation is efficient and simple, and the data acquisition does not require a specific respiratory mode.
Unless defined otherwise, all technical and/or scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application relates. The materials, methods, and examples mentioned herein are illustrative only and not intended to be limiting.
Although the present application has been described in connection with specific embodiments thereof, those skilled in the art will appreciate that various substitutions, modifications and changes may be made without departing from the spirit of the application.
Claims (3)
1. A method for estimating lumped parameter respiratory system model parameters, wherein,
the lumped parameter respiratory system model is expressed as:
P=a 0 V+a 1 V′+a 2 V"+…+a k V (k) +P 0 ;
wherein P is air pressure, V is the 1 st derivative of V' in the input total air formula, and represents the volume flow rate; v' is the 2 nd derivative of V and represents the airflow acceleration; v (V) (k) K-th derivative of V;
P 0 is the initial pressure in the respiratory system;
a 0 、a 1 、a 2 、…、a k coefficient parameters, a, of a lumped parameter respiratory system model 0 Is elastic; a, a 1 Is airway resistance; a, a 2 Is the inertial resistance coefficient;
using matricesObtaining fitting results of various coefficients;
wherein ,m is the number of terms of the lumped parameter respiratory system model;
the matrix X is an n-row m-column matrix;
is an n-dimensional barometric pressure vector, which is sampled by n times of the breathing circuitThe P value of the dot constitutes.
2. The method of estimating lumped parameter respiratory system model parameters as recited in claim 1 wherein:
the initial value of the derivative of V above 2 th order is set to 0, and the value of the kth order derivative term at the other time t is calculated as follows:
where h is the time step.
3. The method of estimating lumped parameter respiratory system model parameters as recited in claim 1 wherein:
when the lumped parameter respiratory system model is a first order linear model, the lumped parameter respiratory system model is expressed as:
P=a 0 V+a 1 V′+P 0 ;
wherein ,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310553702.1A CN116796656B (en) | 2023-05-17 | 2023-05-17 | Method for estimating lumped parameter respiratory system model parameters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310553702.1A CN116796656B (en) | 2023-05-17 | 2023-05-17 | Method for estimating lumped parameter respiratory system model parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116796656A true CN116796656A (en) | 2023-09-22 |
CN116796656B CN116796656B (en) | 2024-03-19 |
Family
ID=88035323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310553702.1A Active CN116796656B (en) | 2023-05-17 | 2023-05-17 | Method for estimating lumped parameter respiratory system model parameters |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116796656B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117174324A (en) * | 2023-10-27 | 2023-12-05 | 中国人民解放军总医院第一医学中心 | Respiratory system modeling method based on hybrid model and electronic equipment |
CN117180573A (en) * | 2023-11-07 | 2023-12-08 | 中国人民解放军总医院第一医学中心 | Respiratory system mechanical ventilation method and system based on double-index model |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180001041A1 (en) * | 2015-02-12 | 2018-01-04 | Koninklijke Philips N.V. | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
CN115317741A (en) * | 2022-08-31 | 2022-11-11 | 深圳市科曼医疗设备有限公司 | Methods, apparatus, devices and media for estimating airway resistance and compliance |
CN115952745A (en) * | 2022-11-01 | 2023-04-11 | 山东大学 | Respiratory system simulation model construction method, simulation platform and storable medium |
-
2023
- 2023-05-17 CN CN202310553702.1A patent/CN116796656B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180001041A1 (en) * | 2015-02-12 | 2018-01-04 | Koninklijke Philips N.V. | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
CN115317741A (en) * | 2022-08-31 | 2022-11-11 | 深圳市科曼医疗设备有限公司 | Methods, apparatus, devices and media for estimating airway resistance and compliance |
CN115952745A (en) * | 2022-11-01 | 2023-04-11 | 山东大学 | Respiratory system simulation model construction method, simulation platform and storable medium |
Non-Patent Citations (2)
Title |
---|
丰继华 等: "集总参数呼吸机模型仿真研究", 系统仿真学报, vol. 18, no. 10, 31 October 2006 (2006-10-31) * |
刘亚杰 等: "基于最小二乘的Mesic呼吸系统模型辨识研究", 航天医学与医学工程, no. 01, 15 February 2009 (2009-02-15) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117174324A (en) * | 2023-10-27 | 2023-12-05 | 中国人民解放军总医院第一医学中心 | Respiratory system modeling method based on hybrid model and electronic equipment |
CN117174324B (en) * | 2023-10-27 | 2024-02-02 | 中国人民解放军总医院第一医学中心 | Respiratory system modeling method based on hybrid model and electronic equipment |
CN117180573A (en) * | 2023-11-07 | 2023-12-08 | 中国人民解放军总医院第一医学中心 | Respiratory system mechanical ventilation method and system based on double-index model |
CN117180573B (en) * | 2023-11-07 | 2024-03-08 | 中国人民解放军总医院第一医学中心 | Respiratory system mechanical ventilation method and system based on double-index model |
Also Published As
Publication number | Publication date |
---|---|
CN116796656B (en) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116796656B (en) | Method for estimating lumped parameter respiratory system model parameters | |
Stinner et al. | SIMULATION OF TEMPERATURE-DEPENDENT DEVELOPMENT IN POPULATION DYNAMICS MODELS12 | |
JPH06510114A (en) | Apparatus and method for calibration of sensor systems | |
Kupper et al. | Age-period-cohort analysis: an illustration of the problems in assessing interaction in one observation per cell data | |
CN116911212B (en) | Respiratory system modeling method based on fractional calculus | |
Crost et al. | Risk and aversion in the integrated assessment of climate change | |
CN102588214A (en) | Method and monitoring apparatus for automated surveillance of a wind turbine and a method for creating a linear model | |
Li et al. | Model selection and structure specification in ultra-high dimensional generalised semi-varying coefficient models | |
CN110781595B (en) | Method, device, terminal and medium for predicting energy use efficiency (PUE) | |
CN108256631A (en) | A kind of user behavior commending system based on attention model | |
CN105078462A (en) | Method and device for estimating air resistant and compliance | |
CN105841847B (en) | A kind of method for estimating Surface latent heat fluxes | |
Van der Werf et al. | Random regression in animal breeding | |
CN108245749A (en) | Respiratory flow method of counting, device, computer equipment and medium | |
CN108121854A (en) | A kind of region Carbon flux evaluation method based on remote sensing data | |
CN115204491A (en) | Production line working condition prediction method and system based on digital twinning and LSTM | |
CN111640507B (en) | Quantum prediction method for human health state | |
Garthwaite et al. | Non-conjugate prior distribution assessment for multivariate normal sampling | |
Kahn et al. | A dynamic model for the simulation of cattle herd production systems: Part 1—General description and the effects of simulation techniques on model results | |
CN117219223A (en) | Medicament monitoring method and system for medical storage cabinet | |
Maggioli et al. | A physically-inspired approach to the simulation of plant wilting | |
CN117180573B (en) | Respiratory system mechanical ventilation method and system based on double-index model | |
CN110555472A (en) | electrocardiosignal data labeling method based on active learning | |
CN117174324B (en) | Respiratory system modeling method based on hybrid model and electronic equipment | |
McNair et al. | The dynamics of age-structured populations with a gestation period: Density-independent growth and egg ratio methods for estimating the birth rate |
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 |