Next Article in Journal
Current Status and Challenges in the Commercial Production of Polyhydroxyalkanoate-Based Bioplastic: A Review
Next Article in Special Issue
Sensorless Position Control in High-Speed Domain of PMSM Based on Improved Adaptive Sliding Mode Observer
Previous Article in Journal
Torque Mapping as a Novel Method for In-Line Characterisation of Spherical Agglomeration Process Kinetics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Identification Method of Respiratory Characteristic Parameters during Mechanical Ventilation

The College of Electrical and Information Engineering, Beihua University, Jilin 132021, China
*
Author to whom correspondence should be addressed.
Processes 2024, 12(8), 1719; https://doi.org/10.3390/pr12081719
Submission received: 30 June 2024 / Revised: 12 August 2024 / Accepted: 13 August 2024 / Published: 15 August 2024
(This article belongs to the Special Issue Design and Analysis of Adaptive Identification and Control)

Abstract

:
In order to enhance the accuracy of ventilator parameter setting, this paper analyzes two identification methods for respiratory characteristic parameters of non-invasive ventilators and invasive ventilators. For non-invasive ventilators, a respiratory characteristic parameter identification method based on a respiration model is established. In this method, the patient’s respiratory sample set is obtained through non-invasive measurements. Experimental results demonstrate that the mean relative error of pulmonary elastance identification was 14.25%, and the mean relative error of intrapulmonary pressure identification was 12.33% using the Romberg integral algorithm. For chronic patients using non-invasive ventilators, the fault-tolerant space for ventilator parameter setting is large; this method meets the requirement of auxiliary setting of non-invasive ventilator parameters. For invasive ventilators, a respiratory characteristic parameter identification method based on the AVOV–BP neural network is established. In this method, the patient’s respiratory sample set is obtained through real-time invasive measurements. Even with small sample datasets, experimental results show that the mean relative error of pulmonary elastance identification and intrapulmonary pressure identification were both 0.22%. For critically ill patients using invasive ventilators, the fault-tolerant space for ventilator parameter setting is small; this method meets the requirement of auxiliary setting of invasive ventilator parameters.

1. Introduction

In clinical practice, ventilators are widely used as an effective means of artificial replacement of autonomic ventilation. Improperly set ventilator parameters can lead to man–machine asynchronous problems, resulting in safety issues such as inadequate oxygen supply, lung damage, difficulty in weaning off the ventilator, and even life-threatening situations for patients. Medical staff typically set ventilator parameters such as ventilation volume and pressure based on the patient’s weight and their own clinical experience. However, the characteristics of the respiratory system vary from person to person, including differences in respiratory rhythm and lung elasticity based on factors such as age, gender, physical condition, and illness. Addressing these individual variations by obtaining each patient’s respiratory system characteristics is crucial for setting precise ventilator parameters and resolving man–machine asynchronous issues.
Respiratory compliance is a parameter used to describe pulmonary respiratory ability and is also known as pulmonary elastance. In numerical terms, respiratory compliance and pulmonary elastance are reciprocal to each other. Pulmonary elastance serves as an important basis for setting the airway pressure limit value of the ventilator. Experimental results from the literature [1] indicate that decreased pulmonary elastance leads to increased air flow amplitude and tidal volume, which potentially results in hyperventilation. Conversely, increased pulmonary elastance may result in hypoventilation [2]. Therefore, it is essential to set ventilator parameters according to the patient’s pulmonary elastance in order to maintain proper mechanical ventilation efficiency.
(1)
Identification of respiratory characteristic parameters of non-invasive ventilators
The structure of the human respiratory system is commonly described by respiration model; the respiratory parameters such as pulmonary elastance and airway resistance can be obtained by parameter identification of the respiration model. The mechanical structure of the human respiratory system is usually equivalent to the electrical form [2,3,4,5,6,7]. At present, there are three main kinds of research on respiration models during mechanical ventilation: respiration models based on pneumatic system theory, respiration models based on fluid mechanics theory, and respiration models based on differential equations.
Non-invasive ventilators are typically used for patients who have clear consciousness and unobstructed respiratory tracts. Patients generally receive respiratory support treatment through mouth and nose masks.
  • Respiration model based on pneumatic system theory
In the literature [1,8,9,10], the mechanical ventilation system is equivalent to a pneumatic system. The established respiration model includes an airflow equation, a pressure equation, and a volume equation. Key parameters in the model include respiratory compliance, airway resistance, and equivalent effective area of throttle. Pulmonary elastance is defined as a constant that is easy to calculate but cannot describe its nonlinear characteristics.
In the literature [1], the air flow dynamic characteristics of the mechanical ventilation of a lung simulator are studied. In the literature [8], an aviation oxygen supply system based on a mechanical ventilation model is studied. In the literature [9], the dynamic characteristics of a mechanical ventilation system with spontaneous breathing are studied. In the above literature, the experimental apparatus comprises a ventilator, a flow sensor, a tube, an artificial simulating lung, a pressure sensor, a data acquisition card, and a computer. The simulation results of the respiration models align with experimental measurement results, verifying the accuracy of the respiration models. The influence of key parameters in the respiration model is studied in the papers.
In the literature [10], an online estimation method for respiratory parameters based on a pneumatic model is studied. The pressure and tidal volume measured by the experimental system are used as input and output data for the respiration model, followed by the use of the recursive least-square method to identify the parameters of the respiratory model.
  • Respiration model based on fluid mechanics theory
In the literature [11,12,13], a time-switching mechanical ventilation piecewise model is established based on fluid mechanics theory. This model, hereinafter referred to as the TS respiration model, includes a flow velocity equation, a flow rate equation, a tidal volume equation, an intrapulmonary pressure equation, and an airway pressure equation. Pulmonary elastance is defined as a constant. By setting ventilator parameters, simulation results for airway pressure, tidal volume, and static pulmonary pressure at the end of the ventilator pipeline can be obtained. The simulation results of the TS respiration model are consistent with the experimental measurement results, verifying its accuracy.
  • Respiration model based on differential equation
In the literature [14,15,16,17], the second-order linear ordinary differential equation models of the respiratory system are presented. This model is hereinafter referred to as the DE respiration model. In the literature, pulmonary elastance is defined as a polynomial [14,15], as a basis function form of the RBF network [16], and as an output function form of the GRNN network [17]. These varying definitions increase the complexity of the respiration model but accurately describe dynamic respiratory processes due to pulmonary elastance, being time-varying variables.
In the literature [14,15], the fuzzy reasoning method is used to identify the parameters of the respiratory model. In the literature [16], the recursive least-square method is used to identify the parameters of the respiratory model.
Non-invasive ventilators are unable to measure the patient’s intrathoracic respiratory pressure, resulting in a crucial absence of respiratory data for identifying respiration model parameters. To address this issue, given that the literature [11,12,13] has confirmed consistency between the simulation data of the TS respiration model and measured data, this paper proposes to use the TS respiration model to simulate patients’ respiratory processes. The respiratory parameters of the respiration simulator are set, and airway pressure, flow rate, tidal volume, and pulmonary pressure data from the TS respiration model are sampled as input data for the DE respiration model. Pulmonary elastance and intrapulmonary pressure are then calculated using numerical integration algorithms and least-squares algorithms.
(2)
Identification of respiratory characteristic parameters of invasive ventilator
Invasive ventilators are commonly used for patients with impaired consciousness or obstructed airways who typically receive respiratory support treatment through tracheal intubation and tracheotomy. These invasive ventilators can provide real-time measurement of patients’ respiratory data, including intrathoracic respiratory pressure.
The neural network is a core technology in machine learning, which is widely applied in various fields such as pattern recognition, intelligent robots, prediction and estimation, biology, medicine, and economics.
In the literature [17], a respiratory sample set is composed of measured respiratory data, a PSO algorithm is used to optimize the smoothing factor of a GRNN network, and the PSO–GRNN network is used to predict pulmonary static pressure.
In the literature [18], the respiratory sample set is composed of oxygenation, ventilation, and acid-base balance respiration data. Based on fuzzy rules and neural networks, a ventilator parameter setting strategy is proposed using expert knowledge.
In the literature [19], a machine learning model is developed for estimating mechanical ventilation parameters for lung health. The model utilizes an inverse mapping of artificial neural networks with the Graded Particle Swarm Optimizer.
In the literature [20], an improved whale optimization algorithm is proposed and applied to optimize the parameters of the LSTM network. The SFE–LSTM model is established for predicting ventilator pressure. Experiments conducted on Kaggle’s open ventilator dataset verify the effectiveness of the model.
In this study, the neural network method was utilized to identify respiratory characteristic parameters of invasive ventilators. Firstly, a neural network model was established and then trained using measured respiratory samples. Subsequently, the model was tested with additional measured respiratory samples. Ultimately, the study achieved the identification of pulmonary elastance and the prediction of intrapulmonary pressure.

2. Materials and Methods

2.1. Identification of Respiratory Characteristic Parameters of Non-Invasive Ventilator

2.1.1. Construction of a Respiratory Sample Set Based on the TS Respiration Model

Mechanical ventilation is a type of fluid movement process that involves changes in the velocity of air flow ( v ( t ) ), flow rate ( F ( t ) ), tidal volume ( V ( t ) ), intrapulmonary pressure ( P A ( t ) ), airway pressure ( P a w ( t ) ), and their interrelationships. The physical process of increasing end-inspiratory breath-holding in intermittent positive pressure ventilation mode can be categorized into five distinct stages: inspiratory stage ( t 0 , t 1 ) , early end-inspiratory platform stage ( t 1 , t 2 ) , end-inspiratory equilibrium stage ( t 2 , t 3 ) , expiratory stage ( t 3 , t 4 ) , and end-expiratory stage ( t 4 , t 5 ) . According to the theory of fluid mechanics [11,12,13], v ( t ) , F ( t ) , V ( t ) , P A ( t ) , and P a w ( t ) in the TS respiration model are as follows:
v ( t ) = K 1 [ P S P A ( t ) ] t ( t 0 , t 1 ) K 3 [ P a w ( t ) P A ( t ) ] t ( t 1 , t 2 ) 0 t ( t 2 , t 3 ) K 4 P A ( t ) t ( t 3 , t 4 ) 0 t ( t 4 , t 5 ) ,
F ( t ) = 6 × 10 4 A v ( t ) ,
V ( t ) = 0 t = t 0 [ V ( t 1 ) + 1000 60 Q ( t ) d t L V ( t 1 ) ] ( 1 C m C + C m ) t 0 < t t 1 V ( t 1 ) + 1000 60 Q ( t ) d t t > t 1 , V ( t 1 ) + 1000 60 Q ( t ) > 0 0 t > t 1 , V ( t 1 ) + 1000 60 Q ( t ) > 0 ,
P A ( t ) = V ( t ) C + P E E P V ( t ) C + P E E P > P E E P P E E P V ( t ) C + P E E P P E E P ,
P a w ( t ) = 0 t = t 0 P A ( t ) + K 2 2 ρ v 2 ( t ) t 0 < t t 1 P a w ( t 1 ) C m 1000 Q ( t ) / 60 C m t 1 < t t 3 , P a w ( t 1 ) C m 1000 Q ( t ) / 60 > P A ( t ) P A ( t ) t 1 < t t 3 , P a w ( t 1 ) C m 1000 Q ( t ) / 60 P A ( t ) P A ( t ) K 5 2 ρ v 2 ( t ) t > t 3 , P A ( t ) K 5 2 ρ v 2 ( t ) > P E E P P E E P t > t 3 , P A ( t ) K 5 2 ρ v 2 ( t ) P E E P ,
where L is the length of tube. K 1 , the inspiratory ventilation resistance coefficient, can be calculated as follows:
K 1 = 16 V T max 2 O p T I 2 d 4 π 2 P s ,
where V T max is the maximum tidal volume, O p is the tidal volume adjustment opening, T I is the inspiratory time, and d is the inner diameter of the endotracheal tube.
K 2 , the airway resistance coefficient during the inspiratory period, can be calculated as follows:
K 2 = 4 200 d ,
K 3 , the resistance coefficient of the anatomical airway and endotracheal tube during the plateau period, can be calculated as follows:
K 3 = 2 ρ K 2 ,
K 4 , the whole resistance coefficient during the expiratory period, can be calculated as follows:
K 4 = 1 K 5 ,
K 5 = 5 300 d ,
K 5 , the resistance coefficient of the anatomical airway and endotracheal tube during the expiratory period, can be calculated as follows:
L V ( t ) in Equation (3) can be calculated as follows:
L V ( t ) = 6.2 × 10 5 ( L 2 ) 2 π 2 p a w ( t 1 ) d t ρ ,
Based on the literature [11,12,13], the parameter values for Equations (1)–(5) are specified in Table 1.
Figure 1 presents the waveform curves of v ( t ) , F ( t ) , V ( t ) , P A ( t ) , and P a w ( t ) in four respiratory cycles.

2.1.2. The DE Respiration Model

The DE respiration model for patients undergoing mechanical ventilation without spontaneous respiration is as follows [14,15,16,17]:
P a w ( t ) + a P · a w ( t ) = k 1 V ( t ) + k 2 V 2 ( t ) + + k n V n ( t ) + r 1 V ˙ ( t ) + r 2 F ( t ) F ( t ) + b F ˙ ( t ) + P e e a + e ( t ) ,
where the physical meanings of F ( t ) , V ( t ) , and P a w ( t ) are consistent with the TS respiration model. P e e a is the alveolar pressure at the end of respiration. e ( t ) represents error values, including equation error and measurement noise. k 1 , k 2 , , k n , a, b, r 1 , and r 2 are constant coefficients [18,19,20].
Pulmonary elastance is defined as a polynomial as follows:
f E ( V ) = k 1 + k 2 V ( t ) + k 3 V 2 ( t ) + + k n V n 1 ( t ) ,
P A ( t ) is intrapulmonary pressure; it can be described as follows:
P A ( t ) = f E ( V ) V ( t ) ,

2.1.3. Matrix Representation of the DE Respiration Model

The matrix form of Equation (12) can be represented as follows:
P a w ( t ) = φ p ( t ) θ p + e ( t ) ,
φ p ( t ) = P ˙ a w ( t ) V ( t ) V n ( t ) V ˙ ( t ) F ( t ) F ( t ) F ˙ ( t ) 1.0 ,
θ p = a k 1 k 2 k n r 1 r 2 b P e e a ,
where φ p ( t ) is the given parameter matrix, and θ p is the parameter matrix to be obtained.

2.1.4. Numerical Integration of the DE Respiration Model

The integral function is f(x), and compound numerical integration divides the integral interval [a,b] into n small intervals: [ x k 1 , x k ] , ( k = 1 , 2 , , n ; x 0 = a , x n = b ) . Calculating the approximate value of the integral on each small interval ( I k = x k 1 x k f ( x ) d x ), the approximate value of the integral on the interval [a, b] is a b f ( x ) d x k = 1 n I k . This paper will compare four types of compound numerical integration formulas, including the trapezoidal formula, Simpson formula, Cotes formula, and Romberg formula. The definitions of the trapezoidal formula, Simpson formula, and Coates formula are shown in Equations (18)–(20).
a b f ( x ) d x = T n h 2 [ f ( a ) + 2 k = 1 n 1 f ( x k ) + f ( b ) ] ,
a b f ( x ) d x = S n h 6 [ f ( a ) + 4 k = 0 n 1 f ( x k + 1 2 ) + 2 k = 1 n 1 f ( x k ) + f ( b ) ] ,
a b f ( x ) d x = C n h 90 [ 7 f ( a ) + 32 k = 0 n 1 f ( x k + 1 4 ) + 12 k = 0 n 1 f ( x k + 1 2 ) + 32 k = 0 n 1 f ( x k + 3 4 ) + 14 k = 1 n 1 f ( x k ) + 7 f ( b ) ] ,
where h = b a n , x k + 1 2 = x k + h 2 , x k + 1 4 = x k + h 4 , and x k + 3 4 = x k + 3 h 4 .
In order to reduce the amount of computation, a recursive algorithm is employed. The equal division number of integral intervals is taken as the value of 2 0 , 2 1 , 2 2 , , 2 n . The recursive formula of the trapezoidal algorithm is as follows:
T 2 0 = b a 2 [ f ( a ) + f ( b ) ] T 2 k = 1 2 T 2 k 1 + b a 2 k i = 1 2 k 1 f [ a + ( 2 i 1 ) b a 2 k ) , k = 1 , 2 , ,
The recursive formula of the Simpson algorithm is
S n = 4 3 T 2 n 1 3 T n ,
The recursive formula of the Coates algorithm is
C n = 16 15 S 2 n 1 15 S n ,
The recursive formula of the Romberg algorithm is
R n = 64 63 C 2 n 1 63 C n ,
The DE respiration model is a continuous time function model. Firstly, each respiratory cycle is sampled at equal intervals (T) to obtain a numerical discrete respiratory sample set, F ( n ) , V ( n ) , and P a w ( n ) , where n = 0 , 1 , 2 , , N 1 . Subsequently, numerical integration is performed with the dichotomy degree of the integral interval being 4, that is, the integral length is 16.
The recursive formula for numerical integration can be uniformly written as follows [14,15,16,17]:
y ( 0 ) ( 0 ) ( n ) = 8 [ y ( n ) + y ( n + 16 ) ] y ( 0 ) ( l ) ( n ) = 1 2 y ( 0 ) ( l 1 ) + 16 2 l i = 1 2 l 1 y [ n + ( 2 i 1 ) 16 2 l ] y ( m ) ( k ) ( n ) = 4 m y ( m 1 ) ( k + 1 ) y ( m 1 ) ( k ) 4 m 1 ,
where n = 0 , 1 , , N 17 and l = 1 , 2 , 3 . When m = 0 , k = 0 , 1 , 2 , 3 , 4 , and y(n) is a trapezoidal sequence. When m = 1 , k = 0 , 1 , 2 , 3 , and y(n) is a Simpson sequence. When m = 2 , k = 0 , 1 , 2 , and y(n) is a Coates sequences. When m = 3 , k = 0 , 1 , and y(n) is a Romberg sequences.

2.1.5. Parameter Identification of the DE Respiration Model Based on Least-Square Method

ϕ p ( n ) is the numerical integral operation result of φ p ( t ) , and e ( n ) is the numerical integration result of e ( t ) . The numerical integration of Equation (15) is as follows [14,15,16,17]:
p a w ( n ) = ϕ p T ( n ) θ p + e ( n ) ,
p is the transposed matrix of p a w ( n ) , and p = [ p a w ( 0 ) , p a w ( 1 ) , , p a w ( N 17 ) ] T . Φ is the transposed matrix of ϕ p ( n ) , and Φ = [ ϕ p ( 0 ) , ϕ p ( 1 ) , , ϕ p ( N 17 ) ] T . e is the transposed matrix of e ( n ) , and e = [ e ( 0 ) , e ( 1 ) , , e ( N 17 ) ] T . Equation (26) can be rewritten as follows:
p = Φ θ p + e ,
θ ^ p is the least-square estimate of θ p ; it can be calculated as follows:
θ ^ p = ( Φ T Φ ) 1 Φ T p ,
where ( θ ^ p ) T = [ a ^ , k ^ 1 , k ^ 2 , , k ^ n , r ^ 1 , r ^ 2 , b ^ , P ^ e e a ] .

2.1.6. Simulation Experiment and Results

F ( t ) , V ( t ) , and P a w ( t ) are all periodic functions. The sample data for one period of these functions was collected with a sampling interval of 0.05 s. As the respiratory cycle was set to be 3 s, there were a total of 60 sampling points (N = 60). These sampling data F ( n T ) , V ( n T ) , and P a w ( n T ) were used as input data for Equation (26), and the parameter matrix θ ^ p of the respiratory equation was then identified according to Equation (28).
When the pulmonary elastance of the DE respiration model is second order, f E ( V ) = k 1 + k 2 V ( t ) . The results of the parameter matrix θ ^ p identified by four numerical integration algorithms are presented in Table 2.
Pulmonary elastance was calculated based on the identified values of k ^ 1 , k ^ 2 , and k ^ 3 , and the resulting f E V curve is shown in Figure 2. The set value of pulmonary elastance in the TS respiration model was 0.7, that is, the pulmonary elastance value was 1 0.7 = 1.428   P a / mL .
The mean absolute error and the mean relative error of pulmonary elastance identification in Figure 2 are presented in Table 3.
The intrapulmonary pressure was calculated based on the f E V curve in Figure 2 and Equation (14). The curve of intrapulmonary pressure fitted for one cycle is shown in Figure 3.
For the four fitting curves of intrapulmonary pressure in Figure 3, the mean absolute error and the mean relative error were calculated. The results of the calculation are presented in Table 4.
The relationship curves between tidal volume ( V ( t ) ) and fitted intrapulmonary pressure ( P A ( t ) ) by four numerical integration algorithms are presented in Figure 4. The value of the TS respiration model curve corresponds to the “tidal volume curve” and the “intrapulmonary pressure curve” in Figure 1.
From the experimental results, it is evident that the final fitted P A ( t ) curve and V P A curve closely align with the TS respiration model. Additionally, the errors of the Romberg numerical integration algorithm for the identification of pulmonary elastance and intrapulmonary pressure are minimal. Although the mean absolute error of the identification results of pulmonary elastance in Table 3 is small, it is evident from Figure 2 that there is significant fluctuation in the identification results of pulmonary elastance within one cycle, and the fitting results exhibit instability.

2.2. Identification of Respiratory Characteristic Parameters of Invasive Ventilator

The BP neural network is widely utilized in pattern recognition and model parameter identification due to its simple structure, as well as its strong approximation and generalization abilities.
However, the weight regulation method of the BP neural network itself tends to be slow and prone to falling into the local minimum. In this paper, the African Vulture algorithm (AVOV) was employed to optimize the weights and thresholds of BP neural networks [21]. The AVOV–BP neural network was established to enhance the convergence speed and identification accuracy of the network [22]. Experimental comparisons were made on the identification results of respiratory characteristic parameters using three types of neural networks: the RNN network, the LSTM network, and the AVOV–BP network.

2.2.1. The AVOV Optimization Algorithm

The principles and steps of the AVOV optimization algorithm are as follows:
(1)
Determine the group optimal:
After initialization, the vultures should be grouped based on their mass. Vultures corresponding to the optimal solution will be placed in the first set, and those corresponding to the suboptimal solution will be placed in the second set.
R ( i ) = Best 1 ,   i f   p i = L 1 Best 2 ,   i f   p i = L 2 ( p ( i ) = F / i = 1 n F i ) ,
where Best1 represents the best vultures and Best2 represents the next best vultures. L1 and L2 are two random numbers in the range [0,1], with their sum being equal to 1. Fi represents the fitness values of the first and second groups of vultures, and n represents the total number of vultures in both groups;
(2)
Calculate the population hunger rate:
t = h ( sin w ( π I i 2 M ) ) + cos ( π I i 2 M 1 ) ,
F = ( 2 k 1 + 1 ) × z × ( 1 I i M ) + t ,
where F indicates that the vulture is full, I i is the current iteration number, M is the maximum iteration number, z is a random number between −1 and 1, h is a random number between −2 and 2, and k 1 is a random number between 0 and 1. When z < 0, the vultures are in a state of starvation, and when z > 0, the vultures are in full state. When F > 1 , the AVOA algorithm is in the exploratory phase. When F < 1 , the AVOV algorithm is in the development phase, and vultures forage around the optimal solution;
(3)
Exploration phase:
P ( i + 1 ) = R ( i ) D ( i ) × F ,   i f   P 1 k P 1 R ( i ) F + k 2 × ( ( u b l b ) × k 3 + l b ) , e l s e ,
D ( i ) = X × R ( i ) P ( i ) ,
where k P 1 is a random number between 0 and 1, P1 is the predetermined exploration parameter, P(i + 1) is the vulture-position vector in the next iteration, F is the satiety rate of vultures in the current iteration, X is a constant, R(i) is one of the best vultures, k2 and k3 are random numbers between 0 and 1, and lb and ub are the upper and lower bounds of optimization, respectively;
(4)
Development phase:
  • Stage I
When F is between 0.5 and 1, AVOA enters stage I of the development phase, implementing two different rotational flight and siege strategies. The selection of the strategies is based on P2, as expressed in Equation (34):
P ( i + 1 ) = D ( i ) × ( F + k 4 ) d ( t ) d ( t ) = R ( i ) P ( i ) ,   if   P 2 k P 2 S 1 = R ( i ) × k 5 × P ( i ) 2 π × cos ( P ( i ) ) S 2 = R ( i ) × k 5 × P ( i ) 2 π × sin ( P ( i ) ) , e l s e R ( i ) ( S 1 + S 2 )
where k P 2 , k 4 , k 5 are all random numbers within the range of [0, 1], F is the satiety rate of vultures in the current iteration, and R(i) is one of the best vultures;
  • Stage II
When F < 0.5 , AVOA enters stage II of the development phase. The movement of the two vultures gathers many types of vultures on the food source, leading to a siege and an aggressive struggle for food among them. Strategies should be selected according to P3, as expressed in Equation (35):
P ( i + 1 ) = A 1 + A 2 2 A 1 = B e s t 1 ( i ) B e s t 1 ( i ) × P ( i ) × F B e s t 1 ( i ) P ( i ) 2 , i f   P 3 k P 3 A 2 = B e s t 2 ( i ) B e s t 2 ( i ) × P ( i ) B e s t 2 ( i ) P ( i ) 2 × F R ( i ) d ( t ) × F × l e v y ( d ) d ( t ) = R ( i ) P ( i ) , e l s e ,
where levy(d) represents the flight operations, and d(t) is the distance between the vultures and the best vultures in both groups.

2.2.2. Implementation of AVOA Algorithm for Optimizing BP Neural Networks

The selection of initial weights and biases in the training process of BP neural networks significantly impacts the model’s performance. This paper adopts a method based on AVOA to optimize the initial weights and biases of BP neural networks, thereby enhancing the training efficiency and prediction performance of the model. The following sections will elaborate on how AVOA optimizes BP neural networks in detail:
(1)
Initialization of network structure and data normalization:
Firstly, the structural parameters of the BP neural network are set, including the number of neurons in the input layer, hidden layer, and output layer. It was assumed that the number of neurons in the input layer, hidden layer, and output layer are n i n , n h i d d e n , and n o u t ι , respectively. Then, the training data is normalized to improve the stability of the training process;
(2)
Random initialization of weights and biases:
The AVOA algorithm performs a global search by simulating the foraging behavior of vultures, where each vulture’s position represents a candidate solution for the BP neural network, i.e., the weights and biases. For the BP neural network, the total number of parameters is set to n t o t a l , including the weights from the input layer to the hidden layer, the weights from the hidden layer to the output layer, and the biases of each layer. The initial positions of the vultures are randomly generated within the search space:
position i = l b + ( u b l b ) r a n d ( n t o t a l )
The lower bound ( I b ) and upper bound ( u b ) represent the parameter limits;
(3)
Definition and calculation of the fitness function:
The fitness function is used to evaluate the quality of each vulture’s position. In this paper, the mean squared error (MSE) is used as the fitness function. For each vulture’s position, the corresponding weights and biases are assigned to the BP neural network, and the MSE on the training set is calculated as follows:
M S E = 1 N i = 1 N ( y i y ^ i ) 2
where N is the number of training samples, y i is the true output of the i -th sample, and y ^ i is the neural network’s predicted output;
(4)
Updating vulture positions:
According to the update rules of the AVOA algorithm, the position update of vultures is divided into two scenarios: random flight and approaching the optimal solution.
Random flight simulates the behavior of vultures randomly moving during the search process:
position i n e w = p o s i t i o n i exp ( r t ) cos ( 2 π r )
where r is a random number, and t is the current iteration number.
Approaching the optimal solution simulates the behavior of vultures gradually approaching the best foraging spot:
p o s i t i o n i n e w = b e s t p o s i t i o n + r a n d n ( 1 , n t o t a l ) b e s t p o s i t i o n p o s i t i o n i
where bestposition is the current optimal position;
(5)
Assigning optimal weights and biases:
After multiple iterations of optimization, the best position found by the AVOA algorithm represents the optimal initial weights and biases for the BP neural network. These weights and biases are assigned to the network, followed by BP algorithm training to further optimize the network parameters. The ultimately optimized BP neural network can significantly enhance training efficiency and prediction accuracy.
By introducing the AVOA algorithm to optimize the initial weights and biases of the BP neural network, the training efficiency and prediction performance of the network can be significantly improved. The random search and gradual approach to the optimal solution characteristics of the AVOA algorithm effectively avoid the common local optimum issues encountered during traditional BP neural network training, thereby achieving better global optimization results. Experimental results demonstrate that the AVOA-optimized BP neural network outperforms traditional methods in complex nonlinear prediction tasks, indicating its broad application prospects.

2.2.3. Simulation Experiment and Results

The 26 groups of measured respiratory data of one respiratory cycle in the literature [14] were used as the sample set, as shown in Figure 5, in which 20 groups of data were used as the training set, and six groups of data were used as the test set.The numbers of neurons in the input layer, hidden layer, and output layer of the AVOV–BP neural network were set to 3, 5, and 1, respectively. The inputs of the AVOV–BP neural network were tidal volume, flow rate, and airway pressure. The output was intrapulmonary pressure. In the AVOV optimization algorithm, the number of vultures was set to 20, and the maximum number of iterations was set to 30.
The numbers of neurons in the input layer, hidden layer, and output layer of the RNN and LSTM neural network were set to 3, 100, and 1, respectively. The inputs of the neural network were tidal volume, flow rate, and airway pressure, and the output was intrapulmonary pressure. Using the Adam algorithm, the weight and bias parameters of the network were dynamically adjusted through a combination of the gradient descent method and the momentum method.
The experimental programs were executed multiple times, and the superior experimental results were recorded. The prediction results and prediction errors of intrapulmonary pressure are presented in Table 5 and Table 6.
According to Equation (14), the pulmonary elastance values are equal to P A ( t ) / V ( t ) . Therefore, by dividing the output of the neural network by the corresponding input tidal volume, the predicted values of pulmonary elastance corresponding to the samples were obtained.
The experimental programs were executed multiple times, and the superior experimental results were recorded. The prediction results and prediction errors of pulmonary elastance are presented in Table 7 and Table 8.
The experimental results indicate that for small samples, the AVOV–BP neural network demonstrated higher prediction accuracy for pulmonary elastance and intrapulmonary pressure compared to both the RNN and LSTM neural networks.

3. Discussion

The main factors that affect the accuracy of identifying respiratory characteristic parameters in non-invasive ventilators are:
(1)
Equations (7) and (10) provide approximate expressions for resistance coefficients K 2 and K 5 , which have been determined through a combination of fluid mechanics principles and practical measurements. The measured value and fitting accuracy directly impact the values of v ( t ) , F ( t ) , V ( t ) , P A ( t ) , P a w ( t ) , and the results of identifying respiratory characteristic parameters and fitting error in P A ( t ) ;
(2)
The values set for parameters in Table 1 directly impact the values of v ( t ) , F ( t ) , V ( t ) , P A ( t ) , P a w ( t ) , and the results of identifying respiratory characteristic parameters and fitting error in P A ( t ) ;
(3)
The TS respiration model and the DE respiration model discussed in this paper are both periodic functions, with the data from each respiratory cycle being consistent. However, there are perturbations in human respiration that cannot be ignored, and the data from each respiratory cycle may not be entirely consistent. As a result, the parameters identified by this method, such as pulmonary elastance, intrapulmonary pressure, and airway resistance, may constrain when compared to the actual situation of patients.
In the literature [10], experimental apparatus measurement data are utilized as the input and output for the respiration model. The parameters of the respiratory model are then identified using the recursive least-square algorithm, with an error rate ranging between 1% and 5%. In the literature [17], measurement data is also used as the input and output for the DE respiration model. The parameters of the DE respiratory model are identified using numerical integration and the recursive least-square algorithm; the average error of intrapulmonary pressure estimation is about 0.004.
In this paper, numerical integration and the least-square method are utilized. However, the results presented in Table 3 and Table 4 indicate large errors in the identification of pulmonary elastance and the prediction of intrapulmonary pressure. These errors arise from the differing definitions of pulmonary elastance in the two models: while it is defined as a constant in the TS respiration model, it is defined as a polynomial in the DE respiration model. Consequently, this disparity directly contributes to substantial errors in the identification of pulmonary elastance and intrapulmonary pressure.
Therefore, it can be inferred that the accuracy of numerical integration and least-squares algorithm for respiration model parameter identification is high. However, due to the use of non-invasive ventilators and functional limitations, it is not possible to directly measure invasive respiratory data. The use of the TS respiration model to simulate respiratory data leads to a significant error in subsequent identification. In future studies, addressing how to unify the consistency of the TS respiration model and the DE respiration model on the definition of pulmonary elastance can solve this problem.
The accuracy of respiratory characteristic parameter identification in invasive ventilators is influenced by several key factors:
(1)
The accuracy of respiratory data samples directly impacts the system parameter identification and data fitting results;
(2)
Increasing the quantity and diversity of respiratory samples can enhance the model’s universality;
(3)
Advancements in neural network and intelligent optimization algorithms have an impact on the results of parameter identification and data fitting.
In the literature [17], the PSO–GRNN network is established, the mean error between the predicted value and the true value of intrapulmonary pressure in the test set is 52.4 Pa, and the mean absolute error between the predicted value and the true value of intrapulmonary pressure in the train set is 8.2 Pa. However, the AVOV–BP network established in this paper demonstrated higher prediction accuracy, with a mean absolute error of only 0.1685 Pa for intrapulmonary pressure prediction.
In the literature [20], in the SFE–LSTM network, the mean absolute error of ventilator pressure prediction is 0.162115 Pa. In the WOA–SFE–LSTM network, the mean absolute error of ventilator pressure prediction is 0.157829 Pa. In the MWOA–SFE–LSTM network, the mean absolute error of ventilator pressure prediction is 0.140571 Pa. The prediction accuracy of ventilator pressure in all three network models is significantly higher than that of the LSTM network in Section 2.2.3 and also higher than that of the AVOV–BP network in this paper. The primary objective of this study is to enhance calculation speed. Therefore, this paper focuses on the parameter identification and intrapulmonary pressure prediction of invasive ventilators using a small sample respiratory dataset. The sample set in this paper only consists of 26 sets of data, while the literature [20] contains a much larger number of datasets, totaling 6,036,000. As a result, the calculation speed of the method in this paper is faster than that in the literature [20]. Considering the speed of calculation and the accuracy of prediction, the AVOV–BP network established in this paper demonstrates superior performance.
The experimental results in Section 2.1.6 and Section 2.2.3 of this paper demonstrate that the accuracy of parameter identification for invasive ventilators is significantly higher than that of non-invasive ventilators. This can be attributed to the fact that the respiratory sample data for invasive ventilators are actual measured values, whereas for non-invasive ventilators, there is a larger fault-tolerant space for parameter setting due to chronic patient usage. In contrast, critically ill patients using invasive ventilators have a smaller fault-tolerant space for parameter setting. Therefore, the two respiratory parameter identification methods meet the specific requirements for auxiliary settings of both invasive and non-invasive ventilators, respectively, in clinical practice.

4. Conclusions

The existing literature only studies one type of respiration model. This paper combined the respiration model based on fluid mechanics theory and the respiration model based on differential equations in the respiratory characteristic parameter identification of non-invasive ventilators. Pulmonary elastance was defined differently under the two respiration models, and the experimental results provided a quantitative relationship of pulmonary elastance in both models, verifying their correlation and consistency. This method simulated respiratory data by setting the resistance coefficient value of the TS respiration model. However, it should be noted that the simulated respiratory data are identical for each respiratory cycle without considering differences in the same patient’s respiratory data during different cycles, leading to potential errors compared to actual patient respiratory situations. Nonetheless, this method does not require the real-time measurement of respiratory data and offers an advantage by avoiding the discomfort caused by the invasive measurement of intrapulmonary pressure in patients. It is particularly suitable for chronic patients using non-invasive ventilators and for settings where high accuracy of ventilator parameter setting is not essential, such as in homes, nursing institutions, and community hospitals.
In this paper, a method based on the AVOV–BP neural network was proposed for the identification of respiratory characteristic parameters of invasive ventilators. The method required real-time monitoring of the patient’s respiratory data. Specifically, if a respiratory cycle lasts 3 s, the respiratory characteristic parameters need to be identified every 3 s using the patient’s actual respiratory data. However, it was found that under Matlab2017b software, the method took about 4 s to complete an operation. It does not meet the requirement for the real-time identification of respiratory parameters. In order to minimize the need for frequent adjustment of ventilator parameters, a designated observation period, such as one hour, can be established in clinical practice. During this time, the patient’s respiratory parameter identification results should be recorded, and the ventilator parameter settings should be adjusted based on both these results and clinical observations.
In this paper, two methods of identifying respiratory characteristic parameters during mechanical ventilation were studied. These methods are suitable for different types of ventilators and patients. Therefore, the superiority of one method over the other cannot be determined solely based on the accuracy of parameter identification. In practical applications, the parameter setting of non-invasive ventilators is simple to operate, allowing even household users to manage after basic training. On the other hand, the parameter setting of invasive ventilators is complex and requires completion by highly skilled medical staff.
In this paper, two respiratory parameter identification algorithms were analyzed and improved, and some results were obtained. However, in clinical applications, there are numerous influencing factors involved in the setting of ventilator parameters, such as individual differences in patients, equipment factors, and medical staff factors. The stability and security of the proposed respiratory characteristic parameter identification algorithm need to be further verified for practical application on a large scale. There is still a distance to go before it can be widely applied.
In our follow-up work, our focus will be on addressing the following issues:
(1)
More real patient data will be collected from clinical practice to enrich respiratory samples of different ages, genders, and disease;
(2)
The study aims to investigate the direct relationship between respiratory characteristic parameters and ventilator parameter settings, as well as to explore the specific operating rules of ventilator settings, so as to provide valuable reference for medical staff and patients’ families;
(3)
For the respiratory characteristic parameter identification algorithm based on the neural network, how to reduce the calculation amount will be studied.
The methods for identifying respiratory characteristic parameters in this paper will serve as a valuable reference for doctors when setting ventilator parameters. Additionally, they offer theoretical support for doctors conducting pathological research, physiological research, and clinical diagnosis.

Author Contributions

Methodology, Y.Z.; writing—original draft preparation, Y.Z.; writing—review and editing, J.B.; investigation, J.B.; software, X.M.; project administration, Y.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Jilin Provincial Department of Science and Technology, China, grant number 20230204093YY.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

  1. Niu, J.; Shi, Y.; Cao, Z.; Cai, M.; Chen, W.; Zhu, J.; Xu, W. Study on air flow dynamic characteristic of mechanical ventilation of a lung simulator. Sci. China Technol. Sci. 2016, 60, 243–250. [Google Scholar] [CrossRef]
  2. Rajagiri, A.; Diong, B.; Goldman, M.; Nazeran, H. Can asthma in children be detected by the estimated parameter values of the augmented RIC mode. In Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, New York, NY, USA, 30 August–3 September 2006. [Google Scholar]
  3. Yi, W.W.; Zhang, Q.; Wang, Y.K. The Study of Constant Pressure Ventilation with no spontaneous respiration. In Proceedings of the International Symposium on Test and Measurement, Beijing, China, 5–8 August 2007. [Google Scholar]
  4. Redlarski, G.; Jaworski, J. A new approach to modeling of selected human respiratory system diseases, directed to computer simulations. Comput. Biol. Med. 2013, 43, 1606–1613. [Google Scholar] [CrossRef] [PubMed]
  5. Shi, Y.; Ren, S.; Cai, M.; Xu, W.; Deng, Q. Pressure dynamic characteristics of pressure controlled ventilation system of a lung simulator. Comput. Math. Methods Med. 2014, 2014, 761712. [Google Scholar] [CrossRef] [PubMed]
  6. Polak, A.G. Analysis of multiple linear regression algorithms used for respiratory mechanics monitoring during artificial ventilation. Comput. Methods Programs Biomed. 2011, 101, 126–134. [Google Scholar] [CrossRef] [PubMed]
  7. Knörzer, A.; Docherty, P.; Chiew, Y.; Chase, J.; Möller, K. An Extension to the First Order Model of Pulmonary Mechanics to Capture a Pressure dependent Elastance in the Human Lung. IFAC Proc. 2014, 47, 1176–1181. [Google Scholar] [CrossRef]
  8. Shi, Y.; Wang, Y.; Cai, M.; Zhang, B.; Zhu, J. An aviation oxygen supply system based on a mechanical ventilation model. Chin. J. Aeronaut. 2018, 31, 197–204. [Google Scholar] [CrossRef]
  9. Hao, L.; Shi, Y.; Cai, M.; Ren, S.; Wang, Y.; Zhang, H.; Yu, Q. Dynamic Characteristics of a Mechanical Ventilation System With Spontaneous Breathing. IEEE Access 2019, 7, 172847–172859. [Google Scholar] [CrossRef]
  10. Shi, Y.; Niu, J.; Cao, Z.; Cai, M.; Zhu, J.; Xu, W. Online Estimation Method for Respiratory Parameters Based on a Pneumatic Model. IEEE/ACM Trans. Comput. Biol. Bioinform. 2016, 13, 939–946. [Google Scholar] [CrossRef] [PubMed]
  11. Wu, D.; Zhao, J. Mathematical model and computers imulation of end-inspiratory plateau ventilation. Int. J. Anesthesiol. Resusc. 2012, 33, 308–313. [Google Scholar]
  12. Gao, J. Analysis the Resistance in the Function of Mechanical Ventilation. Master’s Dissertation, Shanxi Medical University, Shanxi, China, 2008. [Google Scholar]
  13. Wen, Z.F. Modeling and Simulation of Physical Procedures of Time-Cycle Mechanical Ventilation. Master’s Dissertation, North University of China, Shanxi, China, 2012. [Google Scholar]
  14. Kanae, S.; Masanori, M. An ant colony optimization method for fuzzy membership functions and its application to estimate the pulmonary elastance. In Proceedings of the CCC, Xi’an, China, 26–28 July 2013. [Google Scholar]
  15. Nakamichi, M.; Kanae, S.; Yang, Z.J.; Wada, K. Estimation of pulmonary elastance by functional type SIRMs fuzzy reasoning method. In Proceedings of the CCC, Yantai, China, 22–24 July 2011. [Google Scholar]
  16. Kanae, S.; Maeda, K.; Yang, Z.J.; Wada, K. Parameter Estimation of Nonlinear Differential Equation Models of Respiratory System by Using Numerical Integration Technique. In Proceedings of the 14th IFAC Symposium on System Identification, Newcastle, Australia, 29–31 March 2006. [Google Scholar]
  17. Zhang, Y.; Shunshoku, K.; Bai, J.; Zhou, Z. The method of pulmonary static pressure value prediction based on PSO_GRNN network. Chin. J. Sci. Instrum. 2020, 41, 174–184. [Google Scholar]
  18. Wang, Y.; Ma, S.; Hu, T.; Ma, D.; Lian, X.; Wang, S.; Zhang, J. Expert knowledge-based strategies for v entilator parameter setting and stepless adaptive adjustment. J. Biomed. Eng. 2023, 40, 945–952. [Google Scholar]
  19. Venkata, S.S.O.; Koenig, A.; Pidaparti, R.M. Mechanical Ventilator Parameter Estimation for Lung Health through Machine Learning. Bioengineering 2021, 8, 60. [Google Scholar] [CrossRef] [PubMed]
  20. Yang, F. Study of Ventilator Pressure Prediction Method Based on Static Feature Enhancement. Master’s Dissertation, Yanshan University, Heibei, China, 2023. [Google Scholar]
  21. Cui, C.; Dong, Z.; Wang, X. An Improved African Vulture Algorithm. Comput. Simul. 2023, 40, 311–318. [Google Scholar]
  22. Wang, Y.; Li, S.; Deng, D. Prediction algorithm of ultimate bearing capacity of angle steel reinforced by GFRP sheet based on improved IAVOA-BP algorithm. Foreign Electron. Meas. Technol. 2023, 42, 65–73. [Google Scholar]
Figure 1. Simulation waveform curves of the TS respiration model.
Figure 1. Simulation waveform curves of the TS respiration model.
Processes 12 01719 g001
Figure 2. Fitting curves of f E ( V ) by four numerical integration algorithms.
Figure 2. Fitting curves of f E ( V ) by four numerical integration algorithms.
Processes 12 01719 g002
Figure 3. Fitting curves of P A ( t ) by four numerical integration algorithms.
Figure 3. Fitting curves of P A ( t ) by four numerical integration algorithms.
Processes 12 01719 g003
Figure 4. Fitting curves of V P A by four numerical integration algorithms.
Figure 4. Fitting curves of V P A by four numerical integration algorithms.
Processes 12 01719 g004
Figure 5. Measured respiratory data of one respiratory cycle.
Figure 5. Measured respiratory data of one respiratory cycle.
Processes 12 01719 g005
Table 1. Parameter values for the TS respiration model.
Table 1. Parameter values for the TS respiration model.
SymbolQuantityValue
P S gas source driving pressure 2500   P a
Crespiratory compliance 0.7   mL / P a
C m ventilator compliance 0.01   mL / P a
ρ gas density 1.208   kg / m 3
T I inspiratory time1 s
T o expiratory time2 s
O p adjustable opening of tidal volume90%
PEEPpositive end-expiratory pressure0 V
dinner diameter of tracheal tube7 mm
Asectional area of tracheal tube π d / 2 2
Lpipe length (air leakage simulation)0
EIPpercentage of end-inspiration platform time40%
V T max maximum tidal volume0.001
Table 2. The parameter identification results of four numerical integration algorithms.
Table 2. The parameter identification results of four numerical integration algorithms.
ParameterTrapezoidal AlgorithmSimpson AlgorithmCoates AlgorithmRomberg Algorithm
a ^ 0.62190.61960.64990.6589
k ^ 1 0.86581.07141.07311.0706
k ^ 2 0.00148.6052 × 10−48.8286 × 10−48.9825 × 10−4
k ^ 3 −8.9085 × 10−7−5.8165 × 10−7−6.1577 × 10−7−6.3135 × 10−7
r ^ 1 −0.6079−0.2818−0.2475−0.2381
r ^ 2 0.34720.33470.33460.3346
b ^ 10.84715.90585.28525.1619
P ^ e e a −123.8772−57.4888−41.8671−37.6651
Table 3. Identification errors of pulmonary elastance.
Table 3. Identification errors of pulmonary elastance.
Trapezoidal AlgorithmSimpson AlgorithmCoates Algorithm Romberg Algorithm
MAE0.2922 (Pa/mL)0.2023 (Pa/mL)0.2020 (Pa/mL)0.2033 (Pa/mL)
MRE14.18%14.18%14.16%14.25%
Table 4. Identification errors of intrapulmonary pressure.
Table 4. Identification errors of intrapulmonary pressure.
Trapezoidal Algorithm Simpson AlgorithmCoates Algorithm Romberg Algorithm
MAE149.674987.214472.949469.0957
MRE25.67%15.18%12.93%12.33%
Table 5. Prediction results of intrapulmonary pressure.
Table 5. Prediction results of intrapulmonary pressure.
Measured Value384.7000 722.1000 795.3000 500.8000 155.900015.6000
AVOV–BP neural network 384.9361 722.1952795.3560500.5971 156.158415.7622
RNN neural network 384.7000717.522795.3004500.313153.98815.6002
LSTM neural network 262.05596.929721.73527.439303.296153.782
Table 6. Prediction errors of intrapulmonary pressure.
Table 6. Prediction errors of intrapulmonary pressure.
AVOV–BP Neural Network RNN Neural NetworkLSTM Neural Network
MAE0.1685 (Pa)1.1629 (Pa)105.6013 (Pa)
MRE0.22%0.33%174.02%
Table 7. Prediction results of pulmonary elastance.
Table 7. Prediction results of pulmonary elastance.
Measured Value1.42851.42851.42861.42881.42771.4312
AVOV–BP neural network1.4294 1.42871.4287 1.4282 1.43001.4461
RNN neural network 1.42621.41931.43041.42771.41771.4132
LSTM neural network 1.00541.20841.34221.57592.6221.4312
Table 8. Prediction errors of pulmonary elastance.
Table 8. Prediction errors of pulmonary elastance.
AVOV–BP Neural NetworkRNN Neural NetworkLSTM Neural Network
MAE0.1685 (Pa)1.1629 (Pa)105.6013 (Pa)
MRE0.22%0.49%24.17%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Bai, J.; Ma, X.; Xu, Y. Research on the Identification Method of Respiratory Characteristic Parameters during Mechanical Ventilation. Processes 2024, 12, 1719. https://doi.org/10.3390/pr12081719

AMA Style

Zhang Y, Bai J, Ma X, Xu Y. Research on the Identification Method of Respiratory Characteristic Parameters during Mechanical Ventilation. Processes. 2024; 12(8):1719. https://doi.org/10.3390/pr12081719

Chicago/Turabian Style

Zhang, Yuxin, Jing Bai, Xingyi Ma, and Yu Xu. 2024. "Research on the Identification Method of Respiratory Characteristic Parameters during Mechanical Ventilation" Processes 12, no. 8: 1719. https://doi.org/10.3390/pr12081719

APA Style

Zhang, Y., Bai, J., Ma, X., & Xu, Y. (2024). Research on the Identification Method of Respiratory Characteristic Parameters during Mechanical Ventilation. Processes, 12(8), 1719. https://doi.org/10.3390/pr12081719

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop