Variable Sampling Rate based Active Disturbance Control for a Marine Diesel Engine

In this paper, in order to handle the high nonlinearity and the sophisticated disturbance in marine engines, a variable sampling rate based active disturbance rejection controller is developed for engine speed control. In the proposed method, the Active Disturbance Rejection Control (ADRC) is designed with the consideration of the practical application in engine speed control that is known as the Crank-angle (CA) based or event-based sampling and control, which means the sampling interval varies with the engine speed. Such a problem has not been discussed in any previous study regarding the application of ADRC in engine control. To this end, this paper discusses the convergence of the variable sampling rate based Extended State Observer (ESO), as well as its parameters that guarantee stability. To verify the proposed control scheme more properly, a cycle-detailed hybrid nonlinear engine model is employed. Finally, simulations are carried out on the Hardware-in-the-loop (HIL) system to assess the superiority of the proposed strategy. The comparative results with a Fuzzy-Proportional-Integral-Derivative (PID) controller demonstrate that the proposed control scheme has better adaptation to engine speed, load disturbances, and stronger robustness towards model uncertainties, which indicates a promising reduction of time and burden for calibrating the controller. It also proved that the proposed CA based ADRC by variable sampling rate method outperforms the general fixed sampling rate ADRC, which is widely used in previous works. Moreover, the successful application of the proposed algorithm via CA based strategy in a real Engine Control Unit (ECU) indicates its huge potential in practical engine control.


Introduction
Marine diesel engines are widely used in the domain of ship propulsion [1][2][3].For marine main engines, the most crucial task is to retain engine speed steady under desired engine speed in the presence of the inherent instabilities and disturbances coupled with the unpredictable external environment [4].On the one hand, the diesel engine is highly nonlinear and engine speed is related to various aspects, such as the operating condition, coolant and intake air temperature.On the other hand, as the vessels cruise on the sea, the maritime condition seriously affects the engine speed [2,5].In moderate sea conditions, the ocean current introduces wave disturbance to the propeller, causing the engine with undesired speed fluctuation.While in rough or extreme condition, partial or whole the propeller disk would periodically emerge and submerge from the sea surface, leading to large variations in the engine speed [6].Hence, without an accurate and robust speed governor, the variation of engine speed would lead to abnormal operating conditions [1], which decreases the service life of engines and even results in premature failure of the transmission system in case of severe speed fluctuation [7].Moreover, sustained over-speed will cause irreversible damage in the engine system [4], while sharply decreasing in engine speed may make the engine stall, resulting in serious safe problem.From other aspects, for example, the performance of engine speed control has a significant influence on the fuel economy and exhaust emission [5,8].The design of speed controller for marine engines should under the limitation of the increasingly strict emission regulations (see in [9]), which poses a more tough task for the engine speed control.High performance speed controller can help to reduce the fuel consumption and emission in a certain degree.
As can be found in the literature, vast control techniques have been investigated for marine engine speed control.Some of them range from the model-independent methods, such as classical Proportional-Integral-Derivative (PID) control [10], Sliding Model Control (SMC) [2,5], to the model-based approaches, e.g., H ∞ control [11] and Model Predictive Control (MPC) [12,13].Some others are designed by synthesizing classical PID approach with intelligent algorithms, such as fuzzy technique [14,15].Substantial achievement has been obtained, yet there are still some drawbacks among these approaches.For instance, the classical PID controllers need to be readjusted when external conditions deviate far away from their calibrated situations; the chattering phenomenon and parameters tuning in SMC need to be dealt with well; the implementation of MPC is complex and its cost should be considered; the experience and knowledge base are hard to be effectively gained for fuzzy control system.As mentioned above, the engine speed steady is affected strongly by the internal uncertainties and external load disturbances; therefore, additional attempts have been drawn to the Active Disturbance Rejection Control (ADRC), which has proven to be with strong robustness towards uncertainty and disturbance.
ADRC was first proposed by Han [16] to overcome the drawbacks of the conventional PID.In the ADRC theory, both the internal uncertainty and external disturbance are lumped as the "total disturbance", which is treated as a new state of the system then estimated and compensated by the Extended State Observer (ESO) [17].Its considerable performance has been verified by extensive industrial applications, e.g., [18][19][20], and a summary introduction can be found in [21].Recently, for marine engine speed control, some ADRC or ESO based controllers are developed in simulation environments.In [22], a continuous-time ADRC controller was proposed for a MAN B and W low-speed diesel engine and its control effect was assessed on a simplified transfer function engine model.In [23], a combined discrete-time controller based on Cerebellar Model Articulation controller (CMAC) and ADRC was verified on a simplified engine model where the model of propeller and hull dynamics were considered.In the authors' previous work [24], a controller was obtained by integrating nonlinear and linear ADRC and its control effect was tested on a Crank-angle (CA) based Mean Value Engine Model (MVEM) rather than a simple engine model.Although these ADRC related controllers showed good performance, the ADRC methods investigated therein ignore the crucial characteristics in engine speed control (e.g., the event-based fuel injection control), which are discussed as follow.
First, Internal Combustion (IC) engines are typical event-based system, due to the fact that the in-cylinder process is consisted of four distinct stages, i.e., intake, compression, expand, and exhaust, which are CA dependent [25,26].Such characteristics make it more reasonable for the practical speed controller to be CA based rather than time based.By doing such, the control actions can be synchronized with the CA [27].Numerous researchers have studied the benefits of the CA based control.In [25], for a Spark Ignition (SI) engine, linear quadratic regulator (LQR), fuzzy control, and SMC were designed and compared in the CA domain.More recently, in [27], a gain-schedule speed controller was designed in the CA domain for an ethanol engine.
In order to demonstrate more clearly the speed control issue for IC engines, its basic scheme diagram is illustrated in Figure 1.Generally, for IC engine speed control, the crankshaft and camshaft angular position sensors are necessary for providing the exact fuel injection timing.This can be realized by mounting two specific toothed wheels on the crankshaft and camshaft.As shown in Figure 1, for instance, the toothed wheel for crankshaft should be with at least one missing tooth for Electronics 2019, 8, 370 3 of 23 providing an absolute angular position reference, and for the camshaft, at least one tooth for detecting the Top-dead-center (TDC) of a certain cylinder should be provided.More details can be found in [28].The fuel injection can only be executed at a fixed CA, which means the control input (fuel injection mass) can only be updated once for the individual cylinder within a working cycle (4π rad CA for a 4-stroke engine).Such traits mean that if a controller is structured based on a fixed sampling interval, the calculation of control input cannot be synchronized with the in-cylinder powering process unless it is computed at the correct CA position.It is known as the event-based control scheme for engine speed control, which can be treated as a variable sampling interval control issue as the sampling time ∆t n e is related to engine speed, that is where Str represents the stroke of the engine (for 2-stroke engine, Str = 1, for 4-stroke engine, Str = 2), N cyl is the number of cylinders, n e is engine speed in rpm.
Electronics 2018, 7, x FOR PEER REVIEW 3 of 23 providing an absolute angular position reference, and for the camshaft, at least one tooth for detecting the Top-dead-center (TDC) of a certain cylinder should be provided.More details can be found in [28].The fuel injection can only be executed at a fixed CA, which means the control input (fuel injection mass) can only be updated once for the individual cylinder within a working cycle (4 rad CA for a 4-stroke engine).Such traits mean that if a controller is structured based on a fixed sampling interval, the calculation of control input cannot be synchronized with the in-cylinder powering process unless it is computed at the correct CA position.It is known as the event-based control scheme for engine speed control, which can be treated as a variable sampling interval control issue as the sampling time ∆ is related to engine speed, that is where  represents the stroke of the engine (for 2-stroke engine,  =1, for 4-stroke engine,  =2),  is the number of cylinders,  is engine speed in rpm.As stated above, there comes a problem when, in general, according to [29][30][31], the sampling interval is designed to be constant in discrete-time ADRC and it must be small enough to keep accuracy.However, the sampling period for engine speed control varies with the engine speed as shown in Equation (1).Unfortunately, as the best as the authors know, including the IC engines control domain, there is no literature reported concerning digital implementation of ADRC with variable sampling rate.Most of the reported papers focus on the discrete-time ADRC with a fixed sampling rate.Although in [32] the event-triggered ESO is discussed, the sampling interval is still fixed.We should point out that the variable sampling rate-based control (also known as event-based control, or magnitude-driven sampling-based control, etc.) has drawn considerable attention in past years, especially in the domain of wireless sensor network [33].In such a control method, the sampling is based on event, which makes it hard to obtain effective information from sensors or processes by using fixed sampling time [34].A more detailed introduction and advantage of such As stated above, there comes a problem when, in general, according to [29][30][31], the sampling interval is designed to be constant in discrete-time ADRC and it must be small enough to keep accuracy.However, the sampling period for engine speed control varies with the engine speed as shown in Equation (1).Unfortunately, as the best as the authors know, including the IC engines control domain, there is no literature reported concerning digital implementation of ADRC with variable sampling rate.Most of the reported papers focus on the discrete-time ADRC with a fixed sampling rate.Although in [32] the event-triggered ESO is discussed, the sampling interval is still fixed.We should point out that the variable sampling rate-based control (also known as event-based control, or magnitude-driven sampling-based control, etc.) has drawn considerable attention in past years, especially in the domain of wireless sensor network [33].In such a control method, the sampling is based on event, which makes it hard to obtain effective information from sensors or processes by using fixed sampling time [34].A more detailed introduction and advantage of such control scheme can be found in [33][34][35] and the references therein.Hence, it is of importance to study the variable sampling rate ADRC to improve the control performance for engine speed control.In addition to the IC engine domain, this kind of method also can be considered in speed control for DC motors where the rotation speed is measured via photoelectric encoder or limited speed sensors as the sampling rate also related to the rotation speed, e.g., [36].
Second, when ADRC is implemented via digital controllers where the system outputs are sampled, a discretization approach needs to be considered due to its significant influence on the control effect and stability of the control system [29,37], especially when the sampling time is large [30].In previous works that have focused on the application of ADRC in engine speed control, scant attention has been drawn to this issue.In the authors' previous work [24] and another's work [8,23], the Euler approximation method was directly adopted to design ADRC based engine speed controller with a constant sampling interval.Such an issue has not been studied because the aforementioned previous work designed the controller by using a short sampling interval without considering the practical control process in IC engines, i.e., the CA based control scheme.For the CA based control, we found that when engine speed was low the discretization method was not stable due to the large sampling period in the low speed region.Hence, a more accurate discretization method borrowed from [38] is employed in this study.
To test the variable sampling rate of ADRC before applying to a real engine, a CA based cylinder-by-cylinder mean value engine model (MVEM) was used.Such a model can simulate the discrete torque generation and cylinder-by-cylinder variation in IC engines, which makes the design of controller more practical.In [24], it has been proven that the inherent speed fluctuation caused by the aforementioned traits has a significant influence on the performance of ESO.Thus, it is necessary to verify the control effect of ADRC by such a model.The advantages of designing a controller by using CA based engine models can be found in some other recent works [12,24,27,39].In this study, the engine model not only has the traits in [27,39], but also considers the discrete torque generation and cylinder-by-cylinder variation.Hence, it is more practical and reasonable for the algorithm verification.
Besides, it should be pointed out that, in practical, the engine speed could only be measured when the teeth rim goes through the installed sensors.Moreover, the mean engine speed calculated by tooth-to-tooth engine speed within a specified segment of CA (2π/3 rad CA for a 4-stroke 6-cylinder engine) is employed as the input for controller, which contributes extra delay in the close-loop control system [40].This issue has not been considered in previous works concerning the applications of ADRC in diesel engines.Hence, a Hardware-in-the-loop (HIL) system where the engine signals as in a real engine (such as crankshaft and camshaft signals) can be provided and is used to verify the proposed scheme.The Engine Control Unit (ECU) is developed by automatic code generation tool on an NXP Freescale microprocessor MPC5644 from the company of NXP Semiconductors.
As stated above, addressed in this paper is the variable sampling rate ADRC for marine engine speed control with the consideration of practical control process, including the CA based fuel injection and the delay in measuring engine speed.The main contribution of this study is listed as follows: 1.
The variable sampling rate based ADRC is designed for marine engine speed control.

2.
The stability of the variable sampling rate ESO is analyzed and the parameters tuning is discussed for the variable sampling rate ADRC based engine speed controller.

3.
To demonstrate the control effects of the proposed method, abundant simulations are carried out on the CA based engine model by using the HIL system where a real ECU is used.
The rest paper is structured as follows.In the Section 2, a modified MVEM is introduced.In the Section 3, the variable sampling rate ADRC is designed for marine engine speed control, and the stability of variable sampling rate ESO is analyzed.In the Section 4, the simulations for verifying the proposed algorithm is carried out on the HIL system by comparing with a Fuzzy-PID controller.In Section 5, the analysis of the results and future work are discussed.The conclusion is given in Section 6.

Engine Model Description
The target engine is a four-stroke six-cylinder diesel engine equipped with turbocharger.Figure 2 demonstrates the thermodynamic volumes for the engine system.The engine model is comprised of five parts: intake manifold, exhaust manifold, cylinders, intercooler, and turbocharger.The control input is the fuel injection mass per cylinder per cycle (mg/c.c.), denoted as m f , which is controlled by ECU.We assume that the fuel injection mass is the only controlled variable and otherwise controlled variables such as injection advance angle are optimal.

Engine Model Description
The target engine is a four-stroke six-cylinder diesel engine equipped with turbocharger.Figure 2 demonstrates the thermodynamic volumes for the engine system.The engine model is comprised of five parts: intake manifold, exhaust manifold, cylinders, intercooler, and turbocharger.The control input is the fuel injection mass per cylinder per cycle (mg/c.c.), denoted as  , which is controlled by ECU.We assume that the fuel injection mass is the only controlled variable and otherwise controlled variables such as injection advance angle are optimal.

The General MVEM
Note that the MVEM used in this study was originally used in [41,42] and has been validated well therein.Moreover, this MVEM can be accessed online, based on which plenty of research has been published [43].In the authors' previous study [2], marine load torque and wave disturbance are added to study the engine speed control of marine diesel engine for propulsion.More recently, the authors modified it with the consideration of the discrete indicated torque generation process and the unbalanced power ability among cylinders [12,24].Hence, only the main equations are given.More details can refer to the references mentioned above.The intake and exhaust dynamics can be described as follows, where with  =   +   +  ,  = • 10 . is the temperature,  denotes pressure,  means the volume,  is the mass flow rate, the subscript "1" and "2" denote the intake and exhaust manifold, the subscript "3" means the intercooler, the subscript "e" means the cylinders, and the subscript "c" denotes the turbocharger, the subscript "ij" means the mass flow rate from volume "i"

The General MVEM
Note that the MVEM used in this study was originally used in [41,42] and has been validated well therein.Moreover, this MVEM can be accessed online, based on which plenty of research has been published [43].In the authors' previous study [2], marine load torque and wave disturbance are added to study the engine speed control of marine diesel engine for propulsion.More recently, the authors modified it with the consideration of the discrete indicated torque generation process and the unbalanced power ability among cylinders [12,24].Hence, only the main equations are given.More details can refer to the references mentioned above.The intake and exhaust dynamics can be described as follows, where with T is the temperature, p denotes pressure, V means the volume, W is the mass flow rate, the subscript "1" and "2" denote the intake and exhaust manifold, the subscript "3" means the intercooler, the subscript "e" means the cylinders, and the subscript "c" denotes the turbocharger, the subscript "ij" means the mass flow rate from volume "i" to "j".R a and R e are the ideal-gas constants for intake gas and exhaust gas, respectively, η v is the volumetric efficiency, which relates to intake manifold pressure p 1 , engine speed n e , the coefficients c v1 , c v2 , and c v3 .W f is the mass fuel rate injected to cylinders, which can be calculated by using the control input, i.e., m f .V d is the volume of the total cylinders.
Assuming the inter-cooler is ideal, i.e., there is no pressure loss, no mass accumulation, and the efficiency is perfect, resulting in the following equations. ( The turbocharger dynamics can be given as where P denotes the power, subscript "t" means the turbine, n t is the turbine speed, η m is the efficiency of turbocharger, which needs to be calibrated by experimental data, J t is the inertia of the shaft, τ c is the turbo-lag.Finally, for a diesel engine used for propulsion, its rotational dynamics can be given as with , where J e is the total rotary inertia (including the incidental water effect), M ind is the total indicated torque generated by all the cylinders.M total represents all the resistance torque, including the friction torque (decided by the coefficients C f 1 , C f 2 , C f 3 and engine speed n e ), the friction torque M f , the pumping torque M p , the total bounded noise disturbance M noise which can be modeled by bounded Gaussian noise, the load torque M l , the additional torque from the ocean current M wave .q HV is the fuel low calorific value, η ind is the indicated efficiency which can be given as a function of engine speed n e and Air-to-fuel Ratio (AFR) λ. k l is the load factor which is affected by various aspects, such as the ship size and the load on the ship.

The Modified MVEM by CA Resolved Method
To model the cycle details in reciprocating engines, some innovative methods have been proposed in [26,44] where hybrid models (including both discrete and continuous models) are structured to exhibit the discrete torque generation process.They focused more on the discrete behaviors in torque generation and only limited at idle mode.However, some other important characteristics (such as the working imbalance among cylinders and cycle delays) were neglected.Some other methodologies employ the CA based model to model the variations of transient speed, torque, and in-cylinder process during a cycle.The in-cylinder pressure map is required in some of them, e.g., [45], which makes it higher cost and more difficult to be implemented [46,47].Other methods use the Warson or/and Wiebe combustion functions to describe the in-cylinder combustion process, e.g., [48], which not only need more detailed structural parameters of the engine but also demand more computation burden and time [49].
Considering the deficiencies listed above, we previously introduced the idea to modify the common MVEM [26,44,50], wherein we expected to keep the advantages in the MVEM.Meanwhile, to simulate the inherent speed fluctuation caused by discrete torque generation and cyclic deviations among cylinders without making the engine model more complex and harder to be executed in computation, the aforementioned model we carried out.Hence, a cycle-detailed hybrid nonlinear engine model is simply described in this section.For more details about such an engine model, please refer to the authors' previous works [12,51].
To determine the timing sequence of the four strokes in each cylinder during its working cycle, the CA signal ϕ is introduced by where operator 'mod' means modulus.
Compared with the MVEM above, the proposed model treats each cylinder individually.In the proposed engine model, for each cylinder, the indicted torque is replaced by the mean value during the power stage rather than the mean value during the whole in-cylinder cycle.Hence, for a certain cylinder i, the discrete torque generation process can be given as where i = 1, 2, . . ., N cyl , ∅ d is a block pulse function, which represents the discrete torque generation in IC engines.W i f and η i ind denote the fuel mass flow rate and the indicated efficiency of the cylinder i. ϕ F is the firing duration angle (π rad CA for the target engine), ξ i is is defined as the cylinder-by-cylinder variation of the cylinder i.
Then the total indicated torque generated by all the cylinders is For other dynamic equations in the individual cylinder, there is no any difference with the MVEM, only the whole effects from the individual cylinders should be lumped together at the end.They can be given as follows:

The Comparisons Between the Proposed Engine Model and the MVEM
Under the same operating condition, the indicated torque and speed curves for the proposed engine model and MVEM are compared in Figure 3.The right side of Figure 3 denotes the enlarged figures within one engine working cycle.For the proposed engine model, the discrete torque generation process in individual cylinder is clearly illustrated in Figure 3a. Figure 3b shows that the mean effect of the total indicated torque in the proposed engine model is almost the same as that in MVEM.It can be observed from Figure 3c that the inherent speed fluctuation can be modeled in the proposed model, where the speed fluctuation is significantly larger than that in the MVEM.
mean effect of the total indicated torque in the proposed engine model is almost the same as that in MVEM.It can be observed from Figure 3c that the inherent speed fluctuation can be modeled in the proposed model, where the speed fluctuation is significantly larger than that in the MVEM.
Remark 1.The MVEM has been verified in the works provided by Wahlstrom and Eriksson [41,42].In this study, some modifications are made based on the original engine model, including the torque load, the crankshaft rotational dynamic, and the individual cylinder-based torque generation.The speed responses between the two models are almost the same (see Figure 3c), hence, it is reasonable to conclude that the proposed engine model is equivalent with the original MVEM and it can be used to verify the control algorithm.

The Variable Simpling Rate ADRC for Engine Speed Control
For the control design, the MVEM is enough, so that reconstructing the (5) gives Remark 1.The MVEM has been verified in the works provided by Wahlstrom and Eriksson [41,42].In this study, some modifications are made based on the original engine model, including the torque load, the crankshaft rotational dynamic, and the individual cylinder-based torque generation.The speed responses between the two models are almost the same (see Figure 3c), hence, it is reasonable to conclude that the proposed engine model is equivalent with the original MVEM and it can be used to verify the control algorithm.

The Variable Simpling Rate ADRC for Engine Speed Control
For the control design, the MVEM is enough, so that reconstructing the (5) gives Assuming 30   π J e m f M ind = b 0 + ∆b, b 0 is a constant and ∆b is an unknown value, then (10) can be written as follow: where M total is regarded as the generalized disturbances.
Treating f (t) as an extended state, set x 1 = n e , x 2 = f (t), .
x 2 = . f (t), u = m f , meanwhile, consider the discrete-time CA based speed measurement and control, which means the measured output y(t k ) and the control input u(t k ) are held to be constant before the next sampling.Hence, the augmented state space model for (11) can be structured as . where , t k means the real time at the beginning of the sampling instance k, ∆t k,n e can be calculated by (1), which is a function of the engine speed n e .Note that t k is known, as ∆t k,n e for ∀k is calculable, i.e., t k = ∑ k−1 i=0 ∆t i,n e .In order to match with the practical control condition that the sampling and control of engine speed are CA based, a discrete-time observer is designed for the state space model (12) with the variable sampling interval ∆t k,n e at the sampling instance k.As stated in [52], the discretization method is of importance.In general, most previous works made use of the Euler approximation method [53], which has been proven with a serious problem of stability when the sampling period is large.Miklosovic and Gao compared a couple of typical discretization approaches in [52].More recently, Herbst proposed a more accurate discretization method [38], which is borrowed in this study, due to that the sampling interval ∆t k,n e is big when engine speed is low.For simplicity, denoting x(k) = x(t)| t [t k ,t k +∆t k,ne ) , we can get the discrete-time observer as follow: where Note that the above ESO state equation is gained under the following assumptions.
Assumption 1.The lumped disturbance function f (t) is continuous and differentiable.
which means that, during the controllable engine speed range n min e , n max e , the engine speed n e and the total disturbance f (t) would not change very much between two sampling instances.

Stability Analysis of the Variable Sampling Rate based ESO
The sampling period used in Equation ( 13) varies with the engine speed, which is different with the previous works concerning the convergence analysis on the ESO [30].Here, we need to discuss the stability of the variable sampling rate based ESO.
If the following Equation ( 15) is satisfied, we can conclude that the discrete-time ESO with variable sampling period is convergent [30].
According to [38], the error dynamics can be given as where e = e 1 e 2 , e y is the output error.
Under assumption 3, ( 16) can be rewritten as For Equation ( 17), if pair [A d (∆t k,n e ), CA d (∆t k,n e )] is observable, then there exist L c can stabilize Φ(∆t k,n e ).
As the observable matrix is It is obvious that pair [A d (∆t k,n e ), CA d (∆t k,n e )] is observable if ∆t k,n e = 0.It is always can be satisified as the engine speed cannot be infinite (assumption 3).
The characteristic polynomial of Φ(∆t k,n e ) can be given as according to the condition of ∆, we can discuss the spectral radius r(Φ(∆t k,n e )) of Φ(∆t k,n e ) in the following three cases. If . Then the asymptotically stable condition for the system (17) can be given as follow: If ∆ > 0, D Φ(∆t k,ne ) (z) has two real eigenvalues, they are 2 − (l 1 + l 2 ∆t k,n e ) ± √ ∆ /2 .Then, the r(Φ(∆t k,n e )) = 2 − (l 1 + l 2 ∆t k,n e ) + √ ∆ /2.Hence, the asymptotically stable condition for the system (17) can be given as If ∆ < 0, D Φ(∆t k,ne ) (z) has two imaginary eigenvalues, they are 2 − (l The asymptotically stable condition for the system (17) can be got as which can be further simplified to be Remark 2. From Equations ( 20)-( 22), we get a set of constraints for L c to guarantee the stability of error system (17).But for designing the L c easily, we need efficient approach, such as the bandwidth method proposed in [54].It is interesting that, when we employ such method, it is a special condition of Equation ( 21).

Parameters Design
The control law can be given in the form of state feedback by using the estimated state variables, i.e., where r(k) is the reference speed at the sampling instance k, ω c denotes control bandwidth.Finally, the proposed control scheme can be shown in Figure 4.The calculation of ESO and control law is triggered by the CA signal from speed sensors which guarantee the synchronization between the control input and CA.
Remark 2. From Equations ( 20)-( 22), we get a set of constraints for  to guarantee the stability of error system (17).But for designing the  easily, we need efficient approach, such as the bandwidth method proposed in [54].It is interesting that, when we employ such method, it is a special condition of Equation ( 21).

Parameters Design
The control law can be given in the form of state feedback by using the estimated state variables, i.e., where () is the reference speed at the sampling instance k,  denotes control bandwidth.Finally, the proposed control scheme can be shown in Figure 4.The calculation of ESO and control law is triggered by the CA signal from speed sensors which guarantee the synchronization between the control input and CA.The bandwidth approach means   ∆ , () is designed to be in one location, i.e.   ∆ , () = ( − ) .According to [52], we can get Under such parameters, when  ∆ , ≠ 0 , ( 20) is always true, because The bandwidth approach means D Φ(∆t k,ne ) (z) is designed to be in one location, i.e.D Φ(∆t k,ne ) (z) = (z − β) 2 .According to [52], we can get Under such parameters, when ω o ∆t k,n e = 0, (20) Remark 3. Observing Equation ( 24), we know that, actually, ω c should be designed according to engine speed, as the settling time T settle is varying with engine speed.For instance, when engine speed is low, the response of engine speed is slow, so that it cannot reach a short settling time as when it runs at high speed.Hence, in this study, we modify β(∆t k,n e ) to be with where α is a scale factor according to engine speed at instance k and the nominal speed n e0 , in general, the rating condition, ω o,n e0 and ω c,n e0 are the observer bandwidth and control bandwidth under the nominal speed n e0 , respectively, ω o and ω c mean the modifed observer and control bandwidth, respectively.
Remark 4. The design of the modified parameters does not need extra calibration or tuning.Once the parameters under nominal condition is gained, the parameters in other speed conditions can be calculated directly by a scale parameter α according to the engine speed at the sampling instance k.

Evaluation of the Variable Sampling Rate based ADRC Controller for Engine Speed Control by HIL System
In this section, the proposed control scheme is evaluated based on the proposed engine model by using HIL test platform.Note that the load torque is normalized with k l = 0.004, noise load is banded white noise and the cylinder-by-cylinder variation degree is designed to be [ξ 1 , ξ 2 , ξ 3 , ξ 4 , ξ 5 , ξ 6 ] T = [1.0, 0.95, 1.0, 0.95, 0.95, 1.0] T .
The experimental setup is shown in Figure 5.The engine model that is structured in the Matlab/Simulink environment is downloaded into the NI PXI-8880RT which can simulate the engine in real time.The necessary signals for engine control, such as crankshaft, camshaft, are provided to the ECU by the NI PXI-8880RT hardware (including engine speed simulation module, digital, and analog modules).Note that such crankshaft and camshaft signals are modelled strictly as the signals from a real engine (as illustrated in Figure 1).The control input (fuel injected mass) is acquired from the real fuel injectors by the large current acquisition module in the form of injection timing and injection pulse width then transferred into the engine model.The configuration for input and output signals and the settings for simulations are manipulated by the experimental manager software, NI VeriStand.As a result, the control effect of the proposed method can be validated via the closed-loop simulation.
Additionally, it should be noted that the ECU is established based on the NXP Freescale MPC5644 microcontroller and communicates with host PC by CAN (controller area network) bus using Vector's CANape calibration system.The proposed algorithm and the basic code for engine control are programed in the Matlab/Simulink environment via the automatic code generation technology.

The Comparative Controllers
Generally, as commented in [55], it is unreasonable to evaluate a new designed control strategy by comparing it to a simple PID using rough tuned techniques (i.e., Ziegler-Nichols tuning).Hence, we designed a well-tuned Fuzzy-PID for comparative simulations.The fuzzy rules can refer to [14].
In the Fuzzy-PID, an incremental PID structure with anti-windup scheme is employed as follows.
where  (),  (), and  () gain coefficients at the sampling instance k, which was gained by fuzzy technique according to  (), ∆ (),  , and  * , which are the non-saturated and saturated control inputs, respectively;  and  are the lower and upper bounds of the control input, respectively.
Note that the mentioned fuzzy-PID is also based on variable sampling interval, its sampling and control is triggered by CA signal as the same as in the proposed controller.Hence, for simplicity, it is defined to be VSR-Fuzzy-PID.

The Comparative Controllers
Generally, as commented in [55], it is unreasonable to evaluate a new designed control strategy by comparing it to a simple PID using rough tuned techniques (i.e., Ziegler-Nichols tuning).Hence, we designed a well-tuned Fuzzy-PID for comparative simulations.The fuzzy rules can refer to [14].
In the Fuzzy-PID, an incremental PID structure with anti-windup scheme is employed as follows.
where K p (k), K i (k), and K d (k) gain coefficients at the sampling instance k, which was gained by fuzzy technique according to e PID (k), ∆e PID (k), u PID, and u * PID , which are the non-saturated and saturated control inputs, respectively; u min and u max are the lower and upper bounds of the control input, respectively.
Note that the mentioned fuzzy-PID is also based on variable sampling interval, its sampling and control is triggered by CA signal as the same as in the proposed controller.Hence, for simplicity, it is defined to be VSR-Fuzzy-PID.

The Design of Control Parameters for Controllers
The VSR-Fuzzy-PID controller is tuned via the trail-and-error approach.Finally, the initial coefficients are designed to be K p0 = 0.7, K i0 = 0.06, K d0 = 0.015, the ranges of K p (k), K i (k), and K d (k) are given as follows As for the proposed method, according to Equation ( 25) and Remark 3, if the rated operation (engine runs at 1800 rpm with full load) as a nominal condition is chosen, the control parameters under such condition can be given as ω c,n e0 = 12, ω 0,n e0 = 8ω c,n e0 = 96, , b 0 = 25.Then, in other speed regions, the control parameters are 1800 .We define this control scheme to be a variable sampling rate ADRC, i.e., VSR-ADRC.
In order to show the difference between the VSR-ADRC and the general fixed sampling rate ADRC (FSR-ADRC), a FSR-ADRC was also compared.Its control parameters were designed to be the same as in VSR-ADRC, but its sampling and control interval was fixed at 0.01 s.

The Control Performance
To assess the control effect of the proposed VSR-ADRC, the evaluation criterions in overshoot, settling time, integrated absolute error (IAE) of the output error, and the fuel consumption are used.
According to [56], the IAE of the output error can be calculated by

The Influence of Discretization Methods on the Control Effect
To show clearly the control effects obtained using the Euler approximation and the adopted method, the comparisons of them under different engine speed are plotted in Figure 6.
Electronics 2018, 7, x FOR PEER REVIEW 14 of 23 The VSR-Fuzzy-PID controller is tuned via the trail-and-error approach.Finally, the initial coefficients are designed to be  = 0.7,  = 0.06,  = 0.015, the ranges of  (),  (), and  () are given as follows As for the proposed method, according to Equation ( 25) and Remark 3, if the rated operation (engine runs at 1800 rpm with full load) as a nominal condition is chosen, the control parameters under such condition can be given as  , = 12,  , = 8 , = 96, ,  = 25.Then, in other speed regions, the control parameters are  =  , ,  =  , ,  = ( ) .We define this control scheme to be a variable sampling rate ADRC, i.e., VSR-ADRC.
In order to show the difference between the VSR-ADRC and the general fixed sampling rate ADRC (FSR-ADRC), a FSR-ADRC was also compared.Its control parameters were designed to be the same as in VSR-ADRC, but its sampling and control interval was fixed at 0.01 s.

The Control Performance
To assess the control effect of the proposed VSR-ADRC, the evaluation criterions in overshoot, settling time, integrated absolute error (IAE) of the output error, and the fuel consumption are used.
According to [56], the IAE of the output error can be calculated by

The Influence of Discretization Methods on the Control Effect
To show clearly the control effects obtained using the Euler approximation and the adopted method, the comparisons of them under different engine speed are plotted in Figure 6.As shown in Figure 6, when engine speed is low (600 rpm), the engine speed cannot keep steady in the controller where the ESO is discretized using the Euler approximation.Its control input shakes vigorously and tends to diverge, which should be avoided for preventing the engine from damage.However, for the controller where the ESO is discretized by the employed approach, the engine speed As shown in Figure 6, when engine speed is low (600 rpm), the engine speed cannot keep steady in the controller where the ESO is discretized using the Euler approximation.Its control input shakes vigorously and tends to diverge, which should be avoided for preventing the engine from damage.However, for the controller where the ESO is discretized by the employed approach, the engine speed can retain steady.Moreover, the dynamic performance of the Euler approximation-based controller is inferior to the employed approach, as its settling time is significantly longer.
Note that Figure 6 is got with the same control parameters in both controllers.But no matter how the parameters are tuned in the Euler approximation-based controller, it cannot stabilize the engine speed at low speed.The result indicates that the Euler approximation cannot be employed to design a variable sampling rate based ADRC engine speed controller, at least, during a low speed region.In previous papers that have examined the applications of ADRC for engine speed control, this issue has not been discussed due to the fact that, in such papers [8,23], fixed sampling is used and the sampling period is designed to be smaller.Hence, the discussion above can be treated as a reference for future research.

The Speed Tracking Ability
Under the nominal load, the speed tracking curves for the compared controllers are plotted in Figure 7.For the VSR-ADRC, both in the acceleration and deceleration processes, the settling time is less than that in the VSR-Fuzzy-PID and it can reach the desired speed almost without overshoot.Specifically, compared with the VSR-ADRC, during the deceleration process, the VSR-Fuzzy-PID needs one second more before it settles down the engine speed (see: zoom3).In addition, apparent oscillation can be observed (see: zoom3) in the VSR-Fuzzy-PID controller, which would trigger the turbocharger surge phenomenon that is harmful for the engine system [57].Note that, as shown in Figure 7, for the VSR-Fuzzy-PID, although the tracking speed for ramp reference is faster than the proposed controller, the proposed controller needs less time to reach steady state and keeps less speed fluctuation.In terms of the settling time and overshoot, the control effects of VSR-ADRC and FSR-ADRC are similar.However, in a more detailed comparison, as shown in Table 1, the FSR-ADRC is significantly inferior to the VSR-ADRC, especially in indexes of steady-state speed fluctuation and IAE.It is interesting that the control performance difference between the FSR-ADRC and VSR-ADRC gets bigger as the reference speed deviates from 1800 rpm.One reason is that when reference speed is 1800 rpm, the fixed sampling interval in FSR-ADRC (0.01 s) is almost the same as the variable sampling rate in VSR-ADRC (∆t n e = 60Str N cyl •n e = 20 1800 ≈ 0.011 s).It reveals a fact as mentioned in the Introduction, that it is more reasonable to make control actions that can be synchronized with the CA.
Electronics 2018, 7, x FOR PEER REVIEW 15 of 23 can retain steady.Moreover, the dynamic performance of the Euler approximation-based controller is inferior to the employed approach, as its settling time is significantly longer.Note that Figure 6 is got with the same control parameters in both controllers.But no matter how the parameters are tuned in the Euler approximation-based controller, it cannot stabilize the engine speed at low speed.The result indicates that the Euler approximation cannot be employed to design a variable sampling rate based ADRC engine speed controller, at least, during a low speed region.In previous papers that have examined the applications of ADRC for engine speed control, this issue has not been discussed due to the fact that, in such papers [8,23], fixed sampling is used and the sampling period is designed to be smaller.Hence, the discussion above can be treated as a reference for future research.

The Speed Tracking Ability
Under the nominal load, the speed tracking curves for the compared controllers are plotted in Figure 7.For the VSR-ADRC, both in the acceleration and deceleration processes, the settling time is less than that in the VSR-Fuzzy-PID and it can reach the desired speed almost without overshoot.Specifically, compared with the VSR-ADRC, during the deceleration process, the VSR-Fuzzy-PID needs one second more before it settles down the engine speed (see: zoom3).In addition, apparent oscillation can be observed (see: zoom3) in the VSR-Fuzzy-PID controller, which would trigger the turbocharger surge phenomenon that is harmful for the engine system [57].Note that, as shown in Figure 7, for the VSR-Fuzzy-PID, although the tracking speed for ramp reference is faster than the proposed controller, the proposed controller needs less time to reach steady state and keeps less speed fluctuation.In terms of the settling time and overshoot, the control effects of VSR-ADRC and FSR-ADRC are similar.However, in a more detailed comparison, as shown in Table 1, the FSR-ADRC is significantly inferior to the VSR-ADRC, especially in indexes of steady-state speed fluctuation and IAE.It is interesting that the control performance difference between the FSR-ADRC and VSR-ADRC gets bigger as the reference speed deviates from 1800 rpm.One reason is that when reference speed is 1800 rpm, the fixed sampling interval in FSR-ADRC (0.01s) is almost the same as the variable sampling rate in VSR-ADRC (∆ = • = ≈ 0.011s).It reveals a fact as mentioned in the Introduction, that it is more reasonable to make control actions that can be synchronized with the CA.Besides, for the above processes, the intake pressure, exhaust pressure, and turbine rotational speed of the ADRC and Fuzzy-PID controllers are shown in Figure 8.It can be seen that these three variables are smaller in the ADRC during the acceleration processes, which would reduce the pumping power and less the damage of the engine system [2].
Electronics 2018, 7, x FOR PEER REVIEW 16 of 23 140.4 a , 994.4 b , 3422.5 c 1,2,3 represent the dynamic processes of 600 rpm to 1200 rpm, 1200 rpm to 1800 rpm, 1800 rpm to 600 rpm, respectively; a,b,c denote the steady states between the time from 6 s to 8 s, 10 s to 12 s, 18 s to 20 s, respectively .
Besides, for the above processes, the intake pressure, exhaust pressure, and turbine rotational speed of the ADRC and Fuzzy-PID controllers are shown in Figure 8.It can be seen that these three variables are smaller in the ADRC during the acceleration processes, which would reduce the pumping power and less the damage of the engine system [2].

The Attenuating Disturbance Ability
To verify the attenuating disturbance ability of the proposed method, both mutation and wave load disturbances are applied.Figure 9a shows the curve of  (all the resistance torque).The load torque is removed at the time of 2 s and added in again at the time of 3 s; then the wave load is introduced at the time of 5 s. Figure 9b shows that during the mutation load disturbances, both the variations of engine speed and settling time are significantly less in the proposed controller.Especially with regard to the settling time, it is around a half of that in the Fuzzy-PID controller.From the zoomed-in graph in Figure 9b, we know that during the wave disturbance, the speed deviation in the proposed controller is obviously smaller.Other performance indexes for the wave disturbance are given in Table 2, which indicates that the proposed controller can keep less fuel consumption and smaller speed fluctuation in the situation of wave disturbance.When the ship voyages in the sea, it encounters a wave or ocean current.Hence, although the fuel consumption in the proposed controller

The Attenuating Disturbance Ability
To verify the attenuating disturbance ability of the proposed method, both mutation and wave load disturbances are applied.Figure 9a shows the curve of M total (all the resistance torque).The load torque is removed at the time of 2 s and added in again at the time of 3 s; then the wave load is introduced at the time of 5 s. Figure 9b shows that during the mutation load disturbances, both the variations of engine speed and settling time are significantly less in the proposed controller.Especially with regard to the settling time, it is around a half of that in the Fuzzy-PID controller.From the zoomed-in graph in Figure 9b, we know that during the wave disturbance, the speed deviation in the proposed controller is obviously smaller.Other performance indexes for the wave disturbance are given in Table 2, which indicates that the proposed controller can keep less fuel consumption and smaller speed fluctuation in the situation of wave disturbance.When the ship voyages in the sea, it encounters a wave or ocean current.Hence, although the fuel consumption in the proposed controller is only 0.52% less than that in the Fuzzy-PID controller and so we can conclude that the proposed method has considerable fuel economy.
Electronics 2018, 7, x FOR PEER REVIEW 17 of 23 is only 0.52% less than that in the Fuzzy-PID controller and so we can conclude that the proposed method has considerable fuel economy.The estimate error of state  and the estimate of state  for the VSR-ADRC and FSR-ADRC controllers, we demonstrate our findings in Figure 10.From Figure 10a, we know that the state tracking of  is more accurate in the VSR-ADRC.As shown in Figure 10b, the VSR-ADRC consumes less time to make the estimate of extended state  converge to steady-state.As a result, the VSR-ADRC can estimate and compensate the external disturbance more quickly than the FSR-ADRC.Such characteristics are the reason why the VSR-ADRC outperforms the FSR-ADRC.It is easy to understand that, as the sampling and control in the VSR-ADRC are synchronized with the CA based  The estimate error of state x 1 and the estimate of state x 2 for the VSR-ADRC and FSR-ADRC controllers, we demonstrate our findings in Figure 10.From Figure 10a, we know that the state tracking of x 1 is more accurate in the VSR-ADRC.As shown in Figure 10b, the VSR-ADRC consumes less time to make the estimate of extended state x 2 converge to steady-state.As a result, the VSR-ADRC can estimate and compensate the external disturbance more quickly than the FSR-ADRC.Such characteristics are the reason why the VSR-ADRC outperforms the FSR-ADRC.It is easy to understand that, as the sampling and control in the VSR-ADRC are synchronized with the CA based fuel injection and in-cylinder powering process; hence, the VSR-ADRC can control the engine speed more efficiently and obtain better control performance.
Electronics 2018, 7, x FOR PEER REVIEW 18 of 23 fuel injection and in-cylinder powering process; hence, the VSR-ADRC can control the engine speed more efficiently and obtain better control performance.

The Robustness towards Model Uncertainties
In this study, we primarily consider two kinds of uncertainties: the variations of  and  . is affected by the effect of incidental water, which results in dynamic changing in engine speed response.As shown in Equation ( 5), a smaller  leads to a faster response of the engine speed, causing larger overshoot during speed tracking process and bigger speed variation under disturbance conditions.It should be noted that the over-speed is pretty harmful for the engine system which should be well considered in a speed controller.
Therefore, in order to test the robustness of the proposed controller towards the mentioned uncertainties, for the proposed controller and the VSR-Fuzzy-PID controller, the speed response curves of the engine system with a 20% decrease in  and a 10% increase in  under the load disturbance condition mentioned above are plotted in Figure 11.According to Table 2, the corresponding evaluation indexes of Figure 11 are refreshed in Table 3.As shown in Figure 11 and Table 3, the proposed controller still has superiorities in control performance.Especially, as shown in Figure 11a-zoom1, the engine speed in the proposed controller has a margin of 20 rpm before exceeding the maximum engine speed (2000 rpm).On the contrary, the speed in the VSR-Fuzzy-PID controller almost reaches to the maximum speed during the sudden unload process.Hence, we can conclude that the proposed method has better robustness and adaptation towards system uncertainties.

The Robustness towards Model Uncertainties
In this study, we primarily consider two kinds of uncertainties: the variations of J e and k l .J e is affected by the effect of incidental water, which results in dynamic changing in engine speed response.As shown in Equation ( 5), a smaller J e leads to a faster response of the engine speed, causing larger overshoot during speed tracking process and bigger speed variation under disturbance conditions.It should be noted that the over-speed is pretty harmful for the engine system which should be well considered in a speed controller.
Therefore, in order to test the robustness of the proposed controller towards the mentioned uncertainties, for the proposed controller and the VSR-Fuzzy-PID controller, the speed response curves of the engine system with a 20% decrease in J e and a 10% increase in k l under the load disturbance condition mentioned above are plotted in Figure 11.According to Table 2, the corresponding evaluation indexes of Figure 11 are refreshed in Table 3.As shown in Figure 11 and Table 3, the proposed controller still has superiorities in control performance.Especially, as shown in Figure 11a-zoom1, the engine speed in the proposed controller has a margin of 20 rpm before exceeding the maximum engine speed (2000 rpm).On the contrary, the speed in the VSR-Fuzzy-PID controller almost reaches to the maximum speed during the sudden unload process.Hence, we can conclude that the proposed method has better robustness and adaptation towards system uncertainties.

Discussion
The obtained results demonstrate that, compared with the VSR-Fuzzy-PID, the VSR-ADRC gains better adaptation and robust in engine speed control.Such a concept considers the inherent characteristics in IC engines, i.e., the fuel injection and in-cylinder process are CA or event based.In previous studies, with regard to using ADRC in IC engines speed control, there is not much attention drawn to such traits.However, in the practical control of IC engine speed, it is necessary to design the control strategy in the CA domain and some algorithms can only be implemented in the CA domain, such as the torque balancing control [58,59], which needs to treat the individual cylinder independently in the CA domain.As the proposed ADRC method is CA based, therefore, it can be combined with the torque balancing control strategy to further enhance the speed control performance.
In addition, most the studies focus on investigating the ADRC via engine models under simulation condition [22,23].Only [8] provides bench experiment, but the structure used therein is the fixed sample interval based ESO.Nevertheless, in this study, the controller is verified on the HIL system and the employed engine model is CA resolved which is different from the previous studies.Moreover, the control algorithm is implemented on a real ECU.Hence, the controller investigated in this study is of significant meanings which provides a bridge for making use of ADRC into practical engine speed control.
In this study, the control parameters for ADRC are designed via considering the variable sampling rate, so that it is more understandable.In Equation (1), the sampling interval is varying with the engine speed and the control bandwidth  is also designed to vary with the engine speed, resulting in  changing with engine speed as well.In the previous research [24], where a set of fixed parameters are designed for the marine engine speed controller, lead to a compromise control performance for considering the whole operation conditions.It should be noted that the parameters tuning for the proposed scheme is simple.Actually, only two parameters need to be tuned:  and

Discussion
The obtained results demonstrate that, compared with the VSR-Fuzzy-PID, the VSR-ADRC gains better adaptation and robust in engine speed control.Such a concept considers the inherent characteristics in IC engines, i.e., the fuel injection and in-cylinder process are CA or event based.In previous studies, with regard to using ADRC in IC engines speed control, there is not much attention drawn to such traits.However, in the practical control of IC engine speed, it is necessary to design the control strategy in the CA domain and some algorithms can only be implemented in the CA domain, such as the torque balancing control [58,59], which needs to treat the individual cylinder independently in the CA domain.As the proposed ADRC method is CA based, therefore, it can be combined with the torque balancing control strategy to further enhance the speed control performance.
In addition, most the studies focus on investigating the ADRC via engine models under simulation condition [22,23].Only [8] provides bench experiment, but the structure used therein is the fixed sample interval based ESO.Nevertheless, in this study, the controller is verified on the HIL system and the employed engine model is CA resolved which is different from the previous studies.Moreover, the control algorithm is implemented on a real ECU.Hence, the controller investigated in this study is of significant meanings which provides a bridge for making use of ADRC into practical engine speed control.
In this study, the control parameters for ADRC are designed via considering the variable sampling rate, so that it is more understandable.In Equation (1), the sampling interval is varying with the engine speed and the control bandwidth ω c is also designed to vary with the engine speed, resulting in ω o changing with engine speed as well.In the previous research [24], where a set of fixed parameters are designed for the marine engine speed controller, lead to a compromise control performance for considering the whole operation conditions.It should be noted that the parameters tuning for the proposed scheme is simple.Actually, only two parameters need to be tuned: ω c and b 0 .Whereas in the Fuzzy-PID controller, three basic coefficients (K p0 , K i0 , K d0 ,), the ranges of K p (k), K i (k), K d (k), the fuzzy logic, the universe of discourses for e PID (k) and ∆e PID (k), need to be adjusted well, which is a heavy task.
There are still some issues to be considered to further improve the control performance of the VSR-ADRC in a marine engine speed control.As the experimental results show, the variations of rotational inertia J e and load factor k l would deteriorate the control effect of the VSR-ADRC; it would be effective to correct the parameter b 0 then cancel the influence of such variations. 30  π J e m f M ind = b 0 + ∆b indicates that b 0 is associated with J e , m f , and M ind .Besides, m f and M ind contain the information of enigne load (related to k l ) in a certain degree.In [60], an adaptive ESO is employed to improve the performance for a PMSM motor by identifying the load inertia to modify the parameter b 0 in ESO online.Therefore, it is reasonable and practical to improve the control performance by modifying the parameter b 0 online.In the authors' previous study [12], the steady-state gain between system output (n e ) and the control input (m f ) is identified online to improve the control performance.Such steady-state gain can reflect the engine load; hence, the future work can be focused on improving the control performance of the proposed ADRC with the modification of b 0 by identifying the steady-state gain online.

Conclusions
The engine speed control for marine engine is of great importance.The strong system uncertainties and sophisticated external disturbances make it a challenge to design a high-performance speed controller.In this paper, for a marine engine, a CA based ADRC controller is structured by using the variable sampling rate control.With the discretization method proposed in [38], the discrete-time ESO is obtained with variable sampling interval which depends on engine speed.The convergence of such discrete-time ESO is analyzed, and the parameters that can guarantee its stability is given.In order to verify the proposed control approach, a modified MVEM is employed in which the inherent speed fluctuation caused by discrete indicated torque generation and the cylinder-by-cylinder variations are considered.Finally, by using such an engine model, the proposed controller is tested on the HIL system.Substantial experimental results demonstrate that, compared with the fixed sampling rate ADRC and variable sampling rate Fuzzy-PID, the proposed VSR-ADRC has better adaptation and robustness towards load disturbances (mutation load and wave load) and system uncertainties (the variations of rotational inertia and load factor).

Figure 1 .
Figure 1.The schematic diagram for practical engine speed control system.

Figure 1 .
Figure 1.The schematic diagram for practical engine speed control system.

Figure 2 .
Figure 2. The block diagram of the marine engine model.

Figure 2 .
Figure 2. The block diagram of the marine engine model.

Figure 3 .
Figure 3.The comparisons between two engine models: (a) Individual cylinder indicated torque in the proposed hybrid engine model; (b) Total indicated torque in both engine models; (c) The speed responses in both engine models.(The sub-plots in the right side are the enlarged figures during one working cycle).MVEM: Mean Value Engine Model; CA: Crank-angle.

Figure 3 .
Figure 3.The comparisons between two engine models: (a) Individual cylinder indicated torque in the proposed hybrid engine model; (b) Total indicated torque in both engine models; (c) The speed responses in both engine models.(The sub-plots in the right side are the enlarged figures during one working cycle).MVEM: Mean Value Engine Model; CA: Crank-angle.

Assumption 3 .
In engine speed control, for ∀∆t k,n e , n e ∈ n min e , n max e , 0 < n min e < n max e < +∞, the following equations are satisfied

Figure 4 .
Figure 4.The block diagram of the variable sampling rate Active Disturbance Rejection Control (ADRC) for marine engine speed control.

Figure 4 .
Figure 4.The block diagram of the variable sampling rate Active Disturbance Rejection Control (ADRC) for marine engine speed control.

Figure 5 .
Figure 5.The experimental setup by using the Hardware-in-the-loop (HIL) test platform.

Figure 5 .
Figure 5.The experimental setup by using the Hardware-in-the-loop (HIL) test platform.

Figure 6 .
Figure 6.The control performance comparison of different discretization methods for ADRC with variable sampling rate at low speed region: (a) The speed response; (b) The control input.

Figure 6 .
Figure 6.The control performance comparison of different discretization methods for ADRC with variable sampling rate at low speed region: (a) The speed response; (b) The control input.

Figure 7 .Figure 7 .
Figure 7.The control performance comparison of different controllers during the speed tracking mode under standard engine model: (a) The speed response; (b) the control input.

Figure 8 .
Figure 8.The intake manifold pressure, exhaust manifold pressure, and turbine speed responses of the variable sampling rate VSR-Fuzzy-Proportional-Integral-Derivative (PID) and VSR-ADRC controllers during the speed tracking mode under standard engine model.

Figure 8 .
Figure 8.The intake manifold pressure, exhaust manifold pressure, and turbine speed responses of the variable sampling rate VSR-Fuzzy-Proportional-Integral-Derivative (PID) and VSR-ADRC controllers during the speed tracking mode under standard engine model.

Figure 9 .
Figure 9.The control performance comparison of different controllers during the load disturbance conditions under standard engine model: (a) The total resistance torque; (b) the speed responses; (c) the control input.

Figure 9 .
Figure 9.The control performance comparison of different controllers during the load disturbance conditions under standard engine model: (a) The total resistance torque; (b) the speed responses; (c) the control input.

Figure 10 .
Figure 10.The estimate performance comparison of Extended State Observers (ESOs) for VSR-ADRC and FSR-ADRC during the load disturbance conditions under standard engine model: (a) The estimate error of  ; (b) The estimate of  .

Figure 10 .
Figure 10.The estimate performance comparison of Extended State Observers (ESOs) for VSR-ADRC and FSR-ADRC during the load disturbance conditions under standard engine model: (a) The estimate error of x 1 ; (b) The estimate of x 2 .

Figure 11 .Table 3 .
Figure 11.The control performance comparison between the VSR-Fuzzy-PID and VSR-ADRC controllers during the load disturbance conditions under changed engine model: (a) The speed responses; (b) the control input.

Figure 11 .Table 3 .
Figure 11.The control performance comparison between the VSR-Fuzzy-PID and VSR-ADRC controllers during the load disturbance conditions under changed engine model: (a) The speed responses; (b) the control input.Table 3. The performance comparison of different criterions for the controllers with wave load under changed engine model.Controller Speed Fluctuation (rpm) IAE (×10 −3 ) Fuel Consumption (kg×10 −3 ) VSR-Fuzzy-PID 25 28,489.5 9260.6 VSR-ADRC 22 18,783.99247.2

Table 1 .
The performance comparison of different criterions for the controllers during the speed tracking mode under standard engine model.

Table 1 .
The performance comparison of different criterions for the controllers during the speed tracking mode under standard engine model.

Table 2 .
The performance comparison of different criterions for the controllers with wave load disturbances under standard engine model.

Table 2 .
The performance comparison of different criterions for the controllers with wave load disturbances under standard engine model.