An Online Data-Driven LPV Modeling Method for Turbo-Shaft Engines

: The linear parameter-varying (LPV) model is widely used in aero engine control system design. The conventional local modeling method is inaccurate and inefﬁcient in the full ﬂying envelope. Hence, a novel online data-driven LPV modeling method based on the online sequential extreme learning machine (OS-ELM) with an additional multiplying layer (MLOS-ELM) was proposed. An extra multiplying layer was inserted between the hidden layer and the output layer, where the hidden layer outputs were multiplied by the input variables and state variables of the LPV model. Additionally, the input layer was set to the LPV model’s scheduling parameter. With the multiplying layer added, the state space equation matrices of the LPV model could be easily calculated using online gathered data. Simulation results showed that the outputs of the MLOS-ELM matched that of the component level model of a turbo-shaft engine precisely. The maximum approximation error was less than 0.18%. The predictive outputs of the proposed online data-driven LPV model after ﬁve samples also matched that of the component level model well, and the maximum predictive error within a large ﬂight envelope was less than 1.1% with measurement noise considered. Thus, the efﬁciency and accuracy of the proposed method were validated.


Introduction
The behavior of an aero engine can be described by a linear parameter-varying (LPV) model, which can approximate a nonlinear system or a time-varying linear system with a combination of linear time-invariant (LTI) models [1][2][3]. The parameters of LPV models are expressed as a function of the scheduling variables or the operating point variables, capturing the dynamic changes of the nonlinear system accurately [4,5]. The main advantage of LPV models is that they allow the application of well-established linear design techniques to complex nonlinear systems based on the linear matrix inequality [6,7] and linear quadric methods [8]. For example, the online optimization of model predictive control (MPC) or performance seeking control (PSC) can be reduced to a linear program or a quadratic program, thereby achieving high real-time property when an LPV model is used as a predictive model [9,10]. LPV models preserve the advantageous properties of LTI models, such as low modeling cost, while being able to represent a large class of nonlinear systems with high accuracy [11]. Therefore, the LPV model has been widely used in convex system controller design [12,13] and model-based approaches [14,15], such as robust control, fault diagnosis and optimal control.
The LPV and quasi-LPV identification method has attracted much attention in recent years [16][17][18]. Quasi-means that the scheduling parameters of the LPV model include system variables [19]. There are two main methods of LPV (or quasi-LPV) modeling which have been researched in the literature: the global method and the local method [20,21].
Compared with existing studies, the main contribution of this paper can be summarized as follows: (1) A special structure OS-ELM named MLOS-ELM was proposed in order to build a state space model directly and conveniently in the form of an artificial neural network. With the multiplying layer added, the neural network model and the LPV model were linked closely in the same output form. (2) A novel DD-LPV online modeling method was proposed which could linearize the dynamic system of the turbo-shaft engine at any working state online with the instantly gathered data. (3) The DD-LPV online modeling method provided a high-accuracy online predictive model for model-based approaches, such as MPC, and fault diagnoses of the individual aero-engine.
The paper is organized as follows: Section 2 introduces the description of the turboshaft engine's LPV model. Section 3 presents the special structure of MLOS-ELM for DD-LPV modeling. Section 4 introduces the online updating algorithm of the MLOS-ELM. Section 5 discusses the simulation results and Section 6 concludes this paper.

Brief Introduction of Turbo-Shaft Engine
The sketch of the bi-rotor turbo-shaft engine researched in this paper is shown in Figure 1 [42]. During the operation, the compressed air is drawn from the inlet and flows through the compressor to the combustor where fuel is introduced, mixed, and ignited. Then, the hot gas flows through the gas turbine which extracts a portion of the energy to drive the compressor. Finally, the gas exiting the gas turbine flows through the power turbine which extracts the remaining energy to drive the main rotor [43,44]. All the components work cooperatively to provide power to the helicopter. Compared with existing studies, the main contribution of this paper can be summarized as follows: (1) A special structure OS-ELM named MLOS-ELM was proposed in order to build a state space model directly and conveniently in the form of an artificial neural network. With the multiplying layer added, the neural network model and the LPV model were linked closely in the same output form. (2) A novel DD-LPV online modeling method was proposed which could linearize the dynamic system of the turbo-shaft engine at any working state online with the instantly gathered data. (3) The DD-LPV online modeling method provided a high-accuracy online predictive model for model-based approaches, such as MPC, and fault diagnoses of the individual aero-engine.
The paper is organized as follows: Section 2 introduces the description of the turboshaft engine's LPV model. Section 3 presents the special structure of MLOS-ELM for DD-LPV modeling. Section 4 introduces the online updating algorithm of the MLOS-ELM. Section 5 discusses the simulation results and Section 6 concludes this paper.

Brief Introduction of Turbo-Shaft Engine
The sketch of the bi-rotor turbo-shaft engine researched in this paper is shown in Figure 1 [42]. During the operation, the compressed air is drawn from the inlet and flows through the compressor to the combustor where fuel is introduced, mixed, and ignited. Then, the hot gas flows through the gas turbine which extracts a portion of the energy to drive the compressor. Finally, the gas exiting the gas turbine flows through the power turbine which extracts the remaining energy to drive the main rotor [43,44]. All the components work cooperatively to provide power to the helicopter. For smooth and reliable operation of the helicopter, the rotor speed of the power turbine was always kept at 100% to maintain the main rotor speed constant [42]. This was achieved by adjusting the fuel flow input to the engine. When a flight condition of the helicopter changed, such as forward speed, the required power of the main rotor of the helicopter varied accordingly. This meant that the power turbine needed to provide more or less power to meet the requirement of the main rotor and reach a newly balanced state.

LPV Model Description
The general discrete state-space equation of LPV model can be described as: (1) where α is the scheduling parameter, u(k) is the input vector of the LPV model, x(k) is the state variable vector, y(k) is the output vector of the LPV model and A(α), B(α), C(α) and D(α) are time-varying matrices scheduled by α. For smooth and reliable operation of the helicopter, the rotor speed of the power turbine was always kept at 100% to maintain the main rotor speed constant [42]. This was achieved by adjusting the fuel flow input to the engine. When a flight condition of the helicopter changed, such as forward speed, the required power of the main rotor of the helicopter varied accordingly. This meant that the power turbine needed to provide more or less power to meet the requirement of the main rotor and reach a newly balanced state.

LPV Model Description
The general discrete state-space equation of LPV model can be described as: where α is the scheduling parameter, u(k) is the input vector of the LPV model, x(k) is the state variable vector, y(k) is the output vector of the LPV model and A(α), B(α), C(α) and D(α) are time-varying matrices scheduled by α. For turbo-shaft engine, its main dynamic is derived from the cooperative working of the rotor systems [44], which can be described as: where n g is the rotor speed of gas turbine, n p is the rotor speed of power turbine, W c , W g , W p , W r are the power values of the compressor, gas turbine, power turbine and main rotor, respectively, η g and η p are the shaft efficiency values of the gas turbine and power turbine, respectively and J g and J p are the rotor inertia values of the gas turbine and power turbine, respectively. Equation (2) is a basic equation for the turbo-shaft engine's dynamic calculation in its mechanism model (for more detail, refer to [42][43][44]). Thus, in research relevant to the state space model of turbo-shaft engines, n g and n p are most frequently selected as the state variables [42][43][44][45][46]. In addition, n g can reflect both the output power and the dynamic property of the engine [47,48]. Therefore, the LPV model of turbo-shaft engine can be described in a quasi-LPV form: where n g is selected as both the state variable and the scheduling parameter of the LPV equation, which has α(k) = n g (k). Additionally, n p is the controlled variable and another state variable; thus, we have x(k) = [n g (k) n p (k)]. The fuel flow W f is the main control input u(k) = W f (k). P 3 is the total pressure at the compressor exit and T 44 is the total temperature at the power turbine inlet. P 3 , T 44 together with n g , are all limit parameters for turbo-shaft engine security which should be taken into consideration in the control system design process to avoid overspeed, overpressure and overtemperature. Therefore, the output vector of the LPV model was set to y(k) = [y 1 (k) y 2 (k) y 3 (k) y 4 (k)] T = [n g (k) n p (k) P 3 (k) T 44 (k)] T . Note that n g and n p were already chosen as the states, so its corresponding entries in C and D had a constant value of 0 and 1. For P 3 and T 44 , corresponding entries in matrix C and D varied with the scheduling parameter [15].
In the proposed method, the LPV model was built based on the online training neural network. However, the order of magnitude of the parameters in Equation (3) differed from each other considerably. Therefore, all the data used in the model should be normalized into the order of magnitude of proximity to enhance the accuracy of the model and avoid data overwhelming. The corresponding normalized LPV model is: where the bar over the variable means that the variable is a normalized parameter, the A, B, C, D are normalized system and output matrices. The normalization was carried out as follows: where d max and d min are the maximum and minimum values of any given variable d.

Special Structure of MLOS-ELM for DD-LPV Modeling
The online sequential extreme learning machine (OS-ELM) is a kind of single-hiddenlayer feedforward neural network (SLFN). It was adopted to train neural networks due to its faster online learning and universal approximation capabilities. In OS-ELM, the training sets of input vector and target vector come one by one or chunk by chunk. For one training data set {(x in , t)}, the OS-ELM model with L additive hidden nodes and activation function f (·) can be formulated as:ŷ whereŷ is OS-ELM model's output, W i is the weight from the input layer to the i th hidden node, b i is the bias of the i th hidden node, β i is the weight from the i th hidden node to the output layer. For DD-LPV modeling, a special layer called multiplying layer was inserted between the hidden layer and the output layer based on common OS-ELM, as shown in Figure 2. The inputs of the multiplying layer were the outputs of the hidden layer, the normalized states, and the normalized inputs of the LPV model. The output of the multiplying layer node was the product of its inputs.
where the bar over the variable means that the variable is a normalized parameter, the A , B , C , D are normalized system and output matrices.
The normalization was carried out as follows: (5) where dmax and dmin are the maximum and minimum values of any given variable d.

Special Structure of MLOS-ELM for DD-LPV Modeling
The online sequential extreme learning machine (OS-ELM) is a kind of single-hiddenlayer feedforward neural network (SLFN). It was adopted to train neural networks due to its faster online learning and universal approximation capabilities. In OS-ELM, the training sets of input vector and target vector come one by one or chunk by chunk. For one training data set {( , )} in x t , the OS-ELM model with L additive hidden nodes and activation function f (·) can be formulated as: where ŷ is OS-ELM model's output, i W is the weight from the input layer to the i th hidden node, bi is the bias of the i th hidden node, βi is the weight from the i th hidden node to the output layer. For DD-LPV modeling, a special layer called multiplying layer was inserted between the hidden layer and the output layer based on common OS-ELM, as shown in Figure 2. The inputs of the multiplying layer were the outputs of the hidden layer, the normalized states, and the normalized inputs of the LPV model. The output of the multiplying layer node was the product of its inputs.  Multiplying layer Hidden layer Input layer There were one-to-one local connections between the hidden layer and the multiplying layer, whose weight was set to 1. The nodes of the hidden layer were divided into 3 groups according to the total dimension of x (2-dimension) and u (1-dimension), and each group had 2 nodes, determined by heuristics. Therefore, the outputs of the MLOS-ELM become: , W 1j is the weight from the input node to the j th hidden node, b j is the bias of the j th hidden node, β jm is the weight from the j th multiplying layer node to the m th output node. (β jm is related to time k because it is updated online, while W 1j and b j are invariant, as detailed in Section 4).
The normalized LPV model of a turbo-shaft engine described by Equation (4), can be written in the form of matrix entries as: According to Equation (7), the entries of A(α), B(α), C(α), D(α) in Equation (8) were obtained: As can be seen from Equation (9), the non-constant elements of A, B, C, D in Equation (4) were derived from the parameters of the MLOS-ELM, varying with the scheduling parameter α and output weight β. By inserting a multiplying layer between the hidden layer and the output layer, the output function of the MLOS-ELM was presented in the mathematical description form of the LPV model. That meant that the weights and the input scheduling parameter of the MSOL-ELM were given the ability to derive the parameters of the LPV model. In other words, Equation (9) linked the neural network model and the LPV model closely, which made the LPV model able to be derived directly and conveniently by the network's inner parameters at every sample time k. Remark 1. The numbers of the hidden nodes and the multiplying nodes were both set to 6 by heuristics in this research, which was sufficient to achieve an accurate MLOS-ELM model. The smallest number of nodes in each group was 2, to ensure the accuracy of the MLOS-ELM in this research. If each group had only 1 node, MLOS-ELM's accuracy would have been negatively affected. However, if each group had more than 2 nodes, its accuracy would not be improved but the computational burden would increase rapidly. Moreover, they could be set to different numbers (no less than 3) separately. Owing to there being 3 independent state and input variables x 1 , x 2 , u in turbo-shaft engine's LPV model, the nodes in the multiplying layer must be divided into 3 groups, evenly or not. Then, the output of the network would always have the expression in form of Equation (7) which in turn had an LPV form. The even structure presented in the paper was simple and easily described.

Remark 2.
MLOS-ELM's input layer had only one input (the 1-dimensional scheduling parameter α(k) = n g (k)) which was enough to express the variation of the matrices of the DD-LPV model. Other factors that could have affected the matrices, such as Mach number Ma or altitude H, were not used as the inputs because the online training method was adopted. Therein, the output weight of the network was adapted to the current working condition, including the Mach number and the altitude change, so the entries in Equation (9) could automatically adjust according to the flying condition. Moreover, taking more parameters as the inputs of MLOS-ELM could have affected the accuracy of the DD-LPV.

Online Updating Algorithm of MLOS-ELM
OS-ELM was utilized to build the network for its high real-time properties and competence in recursive learning. The convergence of OS-ELM was proven in [49]. At each sample time k, the newly gathered data set was used to update OS-ELM's output weight β.
According to the MLOS-ELM described in Figure 2. the output vector of the hidden layer at sample time k is: where W ∈ R 1×6 and b ∈ R 1×6 are randomly generated at the first time step and kept constant, f (Wα + b) = 1 1+e −(Wα+b) . Then the output vector of the multiplying layer is: The outputs of MLOS-ELM model of Equation (7) can be rewritten in a compact form: where If k = 1, the weight β can be calculated by least square method: However, O(1) ∈ R 1×6 has rank 1 (not full column rank), thus O(1) T O(1) is singular. Then the ridge regression theory is used to rewrite Equation (13): where λ is the generalized regularization parameter. If k > 1, the output weight β was updated recursively online every sampling period to adapt to the new working state of the engine. However, in some instances, the outputs of MLOS-ELM already matched the target outputs well with previous output weight, especially at a steady state. Under such conditions, there was no need to update the output weight at every sampling time, which increased the computational burden. Therefore, a criterion was set to judge whether the output weight needed to be updated according to the output error of the network.
The output error was defined by relative percentage error: If one of MLOS-ELM's output errors err i , i = 1,2,3,4 exceeded the update criterion, the recursive least square method was used to update the weight according to the new gathered data. Otherwise, the weight maintained the same value as the last sampling time.
Then we have the output weight β(k) (k > 1) as: where δ i > 0, i = 1, 2, 3, 4 is the maximum allowable output error of MLOS-ELM, P(1) is defined as Combining Sections 3 and 4, the proposed online DD-LPV modeling method of turboshaft engine based on MLOS-ELM can be summarized in the following steps:

Remark 3.
To be consistent with the description of the discrete state space equation in Equation (4), the outputs of MLOS-ELM and the target vector in Equation (12) included states at time k + 1.
Since the data at time k + 1 could not be collected online at time k in real applications, the DD-LPV model at time k was obtained by the data of the states at time k and the outputs at time k − 1.

Simulation and Discussions
To verify the accuracy of the MLOS-ELM network and the online DD-LPV model derived from it, simulations in the flight envelope were carried out based on the component level model (CLM) of the UH-60 helicopter and T700 engine system platform for lack of rig test data. The accuracy of this integrated UH-60 helicopter and T700 engine model was checked and verified in [44,50] by a great number of tests. The test results prove that the integrated model not only had good reliability, but was able to carry out digital flight simulations of routine missions [44]. The data generated by the integrated model had high fidelity. For more detailed information, refer to [43,44,51].
The output weight update criteria of MLOS-ELM were δ i = 0.005, i = 1, 2, 3, 4 and λ = 10 −8 . The anti-normalized outputs were adopted to evaluate the accuracy of the MLOS-ELM and the DD-LPV model with physical value.
The true value of the outputs was attained by anti-normalization according to Equation (5): In Sections 5.1 and 5.2, all comparisons and relative errors were calculated by antinormalized data.

Approximation Ability Verification of the MLOS-ELM
First, the effectiveness of the proposed MLOS-ELM was tested through simulation.
The target values of t were gathered online from the component level model (CLM) of the turbo-shaft engine. The purpose of this section was to validate whether the output of MLOS-ELMŷ could track the output of CLM accurately.
The simulation was first carried out at altitude 0 m, with the trajectory of the forward speed of the helicopter increasing from 0 m/s to 30 m/s in the first 50 s, remaining invariant for 60 s and then decreasing to 0 m/s in the last 70 s, as shown in Figure 3. The turboshaft engine was controlled to keep the power turbine rotor speed at 100% during the whole process.
The true value of the outputs was attained by anti-normalization according to Equation (5): In Sections 5.1 and 5.2, all comparisons and relative errors were calculated by antinormalized data.

Approximation Ability Verification of the MLOS-ELM
First, the effectiveness of the proposed MLOS-ELM was tested through simulation. The target values of t were gathered online from the component level model (CLM) of the turbo-shaft engine. The purpose of this section was to validate whether the output of MLOS-ELM ŷ could track the output of CLM accurately.
The simulation was first carried out at altitude 0 m, with the trajectory of the forward speed of the helicopter increasing from 0 m/s to 30 m/s in the first 50 s, remaining invariant for 60 s and then decreasing to 0 m/s in the last 70 s, as shown in Figure 3. The turbo-shaft engine was controlled to keep the power turbine rotor speed at 100% during the whole process.  The comparisons of the MLOS-ELM network's output and the CLM's output are shown in Figure 4. In the figure, GT represents the gas turbine and PT represents the power turbine.
It can be seen from Figure 4 that the rotor speed response of the MLOS-ELM was able to match that of the CLM accurately during both the steady state and the dynamic process. The relative error of MLOS-ELM outputs to those of CLM was less than 0.009% for GT rotor speed and 0.003% for PT rotor speed. The error of compressor exit pressure and PT inlet temperature was also very small, which was within 0.18% and 0.12%, respectively. When the helicopter switched from acceleration to constant speed at 50 s and from constant speed to deceleration at 110 s, the network could still track the outputs well (though a sudden change of the forward speed trajectory occurred). The simulation results showed that the special OS-ELM, with an additional multiplying layer, achieved a high degree of approximation accuracy, which implied the effectiveness of the online DD-LPV model derived by MLOS-ELM.

Prediction Ability Validation of the Online DD-LPV Model in Flight Envelope
The online DD-LPV model was designed to be used in MPC or PSC, whose model was to predict the future outputs of the next few sample times. Therefore, in this section, the prediction ability of the proposed DD-LPV model was validated.
For the LPV model derived from the MLOS-ELM at time k, its state space equation matrices were A(k), B(k), C(k), D(k). Its prediction outputs at time k + p were obtained with the input series of the CLM over time k to k + p, and compared with the CLM output at time k + p to verify the prediction ability of the LPV model derived at time k. As k moved along the test horizon, the prediction ability of the models derived at each sampling time were tested. The prediction outputs of the LPV model were acquired by the response calculation of the discrete time system: Proceeding forward from Equation (19), we can obtain x(k + p) and then y(k + p): As a recently proposed data-driven modeling method for aero engines, an equilibrium manifolds expansion (EME) model [35][36][37] was also built on the UH-60/T700 system for prediction ability comparison. The EME model for the turbo-shaft engine was built at an altitude of 1000 m with different power requirements of the helicopter and was extended to the entire flight envelope through similarity transformation [37].
The simulations were carried out at altitudes of 1000 m and 3000 m with the same forward speed variation in Figure 3. The prediction results of k + 5 at 1000 m are shown in Figure 5 and the variations of the DD-LPV model parameters during this process are shown in Figure 6. Figure 7 shows the prediction results of k + 5 at 3000 m and the variation of the corresponding DD-LPV model parameters are shown in Figure 8. Random noises (magnitude 0.2%) were added to the data gathered from the CLM at 3000 m to imitate the measurement noise, which meant that the MLOS-ELM model at 3000 m was trained by the data with noise-and so was the DD-LPV model derived from it. Figure 5 shows that future output prediction at time k + 5 of the online DD-LPV model was able to match that of the CLM more accurately than the EME model when the altitude was set to 1000 m. The relative predictive errors at time k + 5 of DD-LPV model rotor speed to that of CLM were less than 0.05% for GT rotor speed and less than 0.01% for PT rotor speed. The maximum errors of compressor exit pressure and PT inlet temperature were less than 0.13% and 0.12% respectively, which meant that the DD-LPV model achieved high accuracy in the prediction of all output parameters. At 50 s and 110 s, when the forward speed switched, the fluctuation ranges of the prediction error were also very small. Meanwhile, as shown in Figure 5, the EME model's relative predictive errors of these four outputs were less than 0.36%, 0.07%, 2.0% and 1.7%, respectively. Although the accuracy of the EME model, whose similarity transformation was computed at the similar inlet pressure and temperature, was still acceptable at 1000 m, its prediction error was still much greater than that of the online identified DD-LPV model.      Figure 6 shows that some of the parameters, such as a 12 (α), a 21 (α) and b 21 (α), changed very slightly, even their corresponding scheduling parameter α(k) = n g (k) changed within a relatively wide range. Thus, the output weight β(k) of the MLOS-ELM changed accordingly because the parameters of the DD-LPV model were changed not only with the scheduling parameter, but also with the online updated output weight. All the parameters of the LPV model varied smoothly during the forward speed increase and decrease process which guaranteed the DD-LPV model the ability of smooth prediction over the test horizon. (c 11 (α)~c 22 (α), d 11 (α)~d 21 (α) are not shown in the figure since they were fixed values of 1 or 0).
It can be seen from Figure 7 that all the future outputs prediction of the online DD-LPV model was able to track that of the CLM at time k + 5 very well with measurement noise considered. Most of the relative errors fluctuated around the magnitude of 0.2% which was consistent with the magnitude of the measurement noise. The rotor speed prediction error of the LPV model at time instant k + 5 was less than 0.35%. Most compressor exit pressure prediction errors of the online DD-LPV model at time instant k + 5 were less than 0.6%, and the maximum error was less than 0.81%. All the GT inlet temperature prediction errors were less than 0.6%. By contrast, the maximum prediction errors of P 3 and T 44 of the EME model increased greatly, up to 2.62% and 3.17%, respectively, when the flying altitude changed to 3000 m.    Figure 5 shows that future output prediction at time k + 5 of the online DD-LPV model was able to match that of the CLM more accurately than the EME model when the altitude was set to 1000 m. The relative predictive errors at time k + 5 of DD-LPV model rotor speed  Figure 8 shows that the parameters of the DD-LPV model varied in a wider range than that at altitude 1000 m with the same forward speed variation and a much similar scheduling parameter variation range. If the measurement noise was not considered, similar parameter variations resulted, as Figure 6 illustrates. This meant that the changes were aroused by the measurement noise added to the online training data of the MLOS-ELM, which in turn affected the parameters of the DD-LPV model. Even so, all the parameters of the DD-LPV model still changed smoothly during the dynamic process.
To further verify the effectiveness of the proposed online DD-LPV model, online prediction comparations of the DD-LPV model and the EME model (at different altitudes with and without measurement noise) were conducted. The results are listed in Table 1.
It can be seen from Table 1 that the online DD-LPV model was able to predict future outputs more accurately than the EME model at different altitudes within the flight envelope.
Without measurement noise considered, the average prediction errors of N g of the DD-LPV model were between 0.0071% and 0.0115%-much smaller than those of the EME model, which were between 0.1258% and 0.2134%. For N p , the range of the average errors of DD-LPV model was between 0.0011% and 0.0025%, and was also smaller than that of EME model which was between 0.0150% and 0.0222%. For P 3 and T 44 , much better prediction results could be achieved by the online DD-LPV model. For P 3 , the average errors of DD-LPV model were between 0.0172% and 0.0896% while those of EME model were between 0.6233% and 1.1080%. For T 44 , the average errors of DD-LPV model were between 0.0165% and 0.0471% while those of EME model were between 0.5714% and 1.4431%. Additionally, the maximum errors of N g , N p , P 3 and T 44 of DD-LPV model were also obviously better than the EME model. Moreover, the performance of output prediction of the DD-LPV model was more stable, with the most standard deviation lower by more than an order of magnitude, compared with the EME model. Furthermore, the DD-LPV model had no obvious prediction error fluctuations when the altitude and the forward flying speed changed. It identified those changes instantly, according to the online gathered data. On the contrary, from 1000 m to 4000 m, the EME model's average prediction errors of P 3 and T 44 showed a clearly increasing trend from 0.7899% and 0.5714% to 1.1080% and 1.4431%, respectively, mainly contributing to the turbo-shaft engine's poor applicability of similarity transformation. With measurement noise considered, the average prediction errors of the online DD-LPV model's four outputs N g , N p , P 3 and T 44 were within 0.1202%, 0.1125%, 0.2320%, and 0.1895%, respectively, and the maximum errors were less than 0.4447%, 0.3825%, 1.0303% and 0.8425%, obviously better than the EME model. However, it was notable that the average prediction errors of the EME model were not greatly changed with or without the measurement noise, while, for the DD-LPV model, the prediction errors increased greatly with noise considered. This was mainly because the original error of the DD-LPV model without noise considered was much smaller than the magnitude of the noise. As an online training method, noise will affect the update of the output weight, which in turn will affect the accuracy of MLOS-ELM and the DD-LPV model derived from it.

Conclusions
This paper proposed an online data-driven LPV modeling method for turbo-shaft engines using the online sequential extreme learning machine with an additional multiplying layer (MLOS-ELM). Through simulations conducted in a wide fight envelope, several conclusions were drawn.
First, the outputs of the proposed MLOS-ELM were able to track the outputs of the CLM online precisely, which demonstrated that the MLOS-ELM was able to model the turbo-shaft engine dynamic accurately. The insertion of the multiplying layer made the output of the network become a state space LPV form without affecting the approximation accuracy. Hence, the MLOS-ELM could represent a new thought for data-driven modeling methods that deserves further research. It has the potential to adapt to the properties of other individual engine or nonlinear systems.
When the traditional state space modeling methods of for aero engines were carried out based on applying small perturbations on the CLM, the DD-LPV modeling method based on MLOS-ELM provided a CLM-free online linear model identification method. It was able to avoid the modeling errors of the CLM and enhance the accuracy of the LPV model, especially at the dynamic state.
The DD-LPV modeling method was very efficient and convenient. It was able to derive an accurate LTI state space model instantly at each sampling time with the individual engine online gathered data, which made the LPV model free from the influence of engine degradation and usable as an adaptive LPV model. Additionally, it could lower the mechanism modeling burden of physical systems for linearization purposes.
Differing from the EME model, which needed the similarity transformation to be applied to the whole flying envelope, the proposed DD-LPV model was derived online, and was able to match variations in flying conditions well. Future output prediction simulation results showed that the online DD-LPV model could be used to predict the output of the aero engine over the next several steps much more precisely than the EME model, with or without the measurement noise considered.
Therefore, the proposed modeling method based on MLOS-ELM represents a new and promising data-driven online modeling method that is worthy of further research in control system designs where LPV models are widely used, such as MPC, fault diagnosis, performance seeking control, and so on.