Abstract
This paper considers the problem of autopilot design for surface ship with environmental disturbances and input saturation. Extended state observer (ESO) is used to reconstruct yaw rate and estimate unknown environmental disturbances. An auxiliary dynamic system is employed to handle the input saturation. By combining with the auxiliary dynamic system and ESO, an event-triggered controller (ETC) is designed, and energy is saved by reduction in the rudder rate. Simulation results illustrate the effectiveness and feasibility of the event-triggered proposed controller.
1. Introduction
Autopilot design of marine surface vehicles has witnessed a surge of interest. Autopilot design with sophisticated controllers has a wide application prospect in reduced manpower, adequate economy, sufficient reliability, and optimum performance [1].
A significant amount of research efforts has been focused on autopilot. In the beginning, because of simple structure and preferable antinoise ability, proportional integral differential (PID) algorithm is applied to design autopilot. However, the parameters for the PID control terms are not easy to tune. To solve these problems [2], consider a class of optimal PID control problems, each of these optimal parameter selection problems is solved as a nonlinear optimization problem [3], and present an adaptive PID control algorithm which improves the grey predictive control for temperature control. In [4], a nonlinear function is employed to describe manoeuvring characteristics that reflect the steady-state relation between the rudder deflection and the rate of turn of the hull, and two configurations of nonlinear controllers are applied to ship course.
In practice, the actual yaw angle is measured from the sensors or observer [5]. Das [6] makes an assumption that the yaw and yaw rate are all measurable, which is difficult to implement. Peng et al. [7] present a robust adaptive steering law based on dynamic surface control technique and neural network, and a state observer is employed to reconstruct the yaw rate.
Environmental disturbances and model uncertainties are important challenges. Rigatos and Tzafestas [8] propose an adaptive fuzzy controller for the autopilot system, and unknown part of the system dynamics can be approximated with the indirect adaptive control. Liu et al. [9] propose an output feedback adaptive steering law for an unmanned surface vehicle, and a neural network using iterative updating law is employed to handle the unknown dynamic and environmental disturbances. A sliding mode control strategy for the autopilot system is proposed in [10], where a fuzzy logic is used to handle gain and avoid the chattering effect. Xia et al. [11] address the problem of adaptive neural network controller for the surface vessel, and adaptive technique and neural network system are used to approximate uncertainties. Disturbance observer is employed to estimate the unmodeled external forces in [12, 13].
All of the above articles do not consider the limitation of rudder angle, which is widespread in practice. Input saturation degrades the quality of the control and can even destabilize the system. Hence, it is necessary to consider the limitation of the rudder angle. Tzeng and Lin [14] propose an adaptive ship steering approach which is characterized by a model-based design approach that provides clear connections between the structure of the controller and ship model. Babaei et al. [15] present an efficient strategy to design the altitude hold mode autopilot for unmanned aerial vehicle, and a multiobjective genetic algorithm is used to mechanize the optimal determination of fuzzy logic controller parameters.
In many applications, rudder limits includes magnitude limits and rate limits. The digital computer with sampler and zero-order holder is employed to emulate the continuous-time control scheme, and the control command will be updated periodically [16]. Katayama and Aoki [17] propose the use of a performance-based supervisory switching control system for trajectory tracking unmanned surface vessels without considering the actuator magnitude limits and rate limits. The performance-based simulator includes the effects of actuator rate limits and saturation, and the results show that the orientation cannot stabilize to the tangent velocity. Furthermore, high-frequency action will lead to the mechanical abrasion and increase fuel consumption. Shipping is no different than other industries and is highly affected by fuel prices [18]. Thus, it is practical to lower the acting frequencies of actuators. An event-triggered robust fuzzy control scheme is developed for underactuated ships in [19], where fuzzy logic systems are employed to approximate the uncertainties, and the ideology of ETC is employed to avoid the frequent acting of actuators.
This paper reports an autopilot design for surface vehicle, subject to environmental disturbances and input saturation. Extended state observer (ESO) is employed to estimate the unknown environmental disturbances. An auxiliary dynamic system is designed to handle the input saturation. The control architecture is developed by ESO, auxiliary dynamic system, and event-triggered strategy. By using the event-triggered strategy, the proposed approach is able to reduce the rudder rate, which has great significance to the reduction of mechanical wear and fuel consumption. The stability analysis using Lyapunov function shows that all error signals in the closed-loop system are uniformly ultimately bounded. Simulation result shows the effectiveness of the proposed method.
This paper is organized as follows. Section 2 formulates the surface ship model with input saturation and the problem formulation. Section 3 depicts the design of the extended state observer. Section 4 presents the design of the controller by using the event-trigger strategy and auxiliary dynamic system, gives the stability analysis, and demonstrates the effectiveness of the proposed control scheme. Section 5 provides a simulation to illustrate the theoretical result. Section 6 concludes this paper.
2. Mathematical Model
2.1. Notation
The following notations will be used throughout this paper. denotes the absolute value unless specified. represents the Euclidean norm. represents the -dimensional Euclidean Space. represents minimum eigenvalue of a square matrix. denotes the transpose of a matrix.
2.2. Model of the Ship
As shown in Figure 1, assume that the ship has an -plane of symmetry, and heave, pitch, and roll motions are neglected. The body-fixed frame coordinate origin is set in the center-line of the ship. The mathematical model of the ship moving in a horizontal plane is described as [20]where is the vector denoting the ship position and yaw angle with coordinates in the earth-fixed frame, is the vector denoting surge, sway, and yaw velocities of the ship in the body-fixed frame, is the vector representing environmental disturbances, and is the control vector of the ship which consists of the surge force , sway force , and yaw moment . The matrices is given bywhere is the mass of the ship, is the moment of inertia about z-axis in the b-frame, is the center of gravity in the body-fixed frame, and , , and are the hydrodynamic parameters:where , , , and . represents a skew-symmetric matrix of Coriolis and centripetal term:where , , , , , and , , , , , , and are linear and quadratic drag coefficients.

2.3. Model Decoupling
The mathematical model represented by equation (1) can be simplified by supposing a constant surge speed , The sway-yaw dynamics can be separated from model and is described aswhere is the state vector, is a vector which represents the force and moment in sway and yaw directions, respectively, and
One notes that can be attained bywhere
Assume that the ship has a single rudder, and the following is obtained:where , and are the coefficients of force in sway and moment in yaw, respectively.
In this model, equation (5) can be transformed into Nomoto’s model by eliminating the sway velocity written aswhere is the parameter related to the characteristics of the ship, is the yaw rate, is the rudder angle, and is the gain constant related to the hydrodynamic coefficients. Experimentally, the constants and are nearly equal. Hence, equation (10) can be written aswhere .
Consider input saturation, and it can be attained that , where is the maximum rudder deflection. Define the mismatch function between input without saturation and with saturation as , where is the rudder angle calculated by the controller.
2.4. Control Objective
The desired angle provided by the guidance system, in here, it is given directly. The desired yaw rate is the signals derived from the following second-order filter:where is a positive constant.
The control objective is to force the surface vehicle to track a reference heading signal and lower the acting frequencies of actuators, such thatwhere is a positive constant which can be made small enough.
3. Extended State Observer
The proposed autopilot controller is designed based on ESO which is designed to reconstruct the unmeasurable yaw rate, environmental disturbances, and model uncertainties.
Reshape dynamic (11) with the mismatch function to facilitate the controller design:where is unknown and bounded environmental disturbances.
In order to facilitate the ESO design, equation (14) is rewritten in the state space form aswhere , , , , , and .
Assumption 1. There exists a positive constant satisfying .
Remark 1. Assumption 1 is reasonable because of the environmental disturbances assumed to be bounded. represents the plant dynamics, represents the environmental disturbances, represents the differential of environmental disturbances, and is input signal.
Expression (15) can be rewritten in terms of the matrix as follows:where is the extended state vector,Consider an ESO as follows:where , is the estimation of , and is observer gain matrix with , is positive constant.
Subtracting equation (18) from equation (16), the following can be attained:where .
Define that with . Hence, equation (19) can be rewritten as follows:Select the appropriate matrix to stabilize the matrix . Under Assumption 1, dynamic (20) shows observer (18) ensures that the estimation of the environmental disturbances converges to its real value in a finite time.
4. Controller Design
Auxiliary dynamic system is constructed to handle the input constraint in [21], in here, the auxiliary dynamic system is designed aswhere and are positive constants.
The controller design process is presented.
Step 1. Define the surface error vector asand choose the Lyapunov candidate function asAlong (23) and (21), the time derivative of is represented bywhere is the estimate of .
Define the second surface error aswhere is the signals derived from the following first-order low-pass filters with is a positive constant and is a virtual control defined aswhere is a constant.
Synthesizing with equations (25) and (27),where .
Step 2. Choose a Lyapunov function as follows:From (18), (22), and (26), the following can be attained:Substituting (30) into the time derivative of and synthesizing with equation (28),To stabilize , a controller is designed aswhere is a positive constant.
Substituting (32) into (31), thenAs shown in [22], it can prove that is bounded. Let be the upper bound of . Using Young’s inequality,
Step 3. Consider the total Lyapunov function candidate asDifferentiating along (20)–(22), (34), and (35), the following can be attained:where and is given according to .
Define is upper bound of , and using Young’s inequality,Then, (37) can be rewritten aswhere .
Step 4. In the common continuous control schemes, the design has already been complete in Step 3. Here, the ideology of event-trigger control is employed and the triggering condition will be established in this step. During the flow period between two successive triggering instants, zero-order hold is employed to keep unchanged. The key to a successful ETC design is to select an appropriate trigger condition. The current triggering instant is denoted as and the next triggering instant is determined by satisfying the following triggering condition:where and is the adjustable variable. Once the triggering condition is satisfied, thus , and it rendersand it implies is bounded. Then,from the definition of , it is sure that , , , and are bounded and convergent. From , then and are bounded. Owing to the convergence of proved, then the convergences of and are proved. Therefore, all signals in the closed-loop are bounded.
Thus, (40) can be rewritten aswhere c is a positive constant, and if the triggering condition is satisfied, , otherwise .
The following can be attained:From the definition of , it can conclude that , , , , and are bounded. From (44), exponentially converges to the compact set that can be made arbitrarily small by adjusting . Then, it holds , and the control object is completed.
5. Simulations and Results
In this section, comparison with normal PID controller (PIDC) is given through Matlab simulation, and simulation result is provided to validate the effectiveness and feasibility of the proposed controller. The parameters for ESO are set to , , and . The rudder angle is constrained as . The sampling period of trigger condition is 0.01 s, i.e., . = 0.01, the reference heading signal is set as for , for .
The output responses of the proposed controller are plotted in Figure 2, where it shows that the actual yaw of the surface ship tracks the reference yaw accurately. By comparison with normal PIDC, ETC has the problems of tiny jitter. Figure 3 shows that ETC has a shock in the yaw rate, and the yaw rate is a smooth curve varying with time under PIDC, and the estimation of the yaw rate converges nicely to the actual yaw rate. Figure 4 describes the environmental disturbances and estimation, and it shows that the ESO can accurately estimate the environmental disturbances. Figure 5 shows the rudder angle varying with time under ETC and PIDC, and the local enlarged diagrams of control input around 65 s for the rudder is shown. The rudder angle is updated 30,000 times under PIDC and 8,512 times under ETC. Though the simulation on the system control part carries on the contrast analysis, PID control strategy reflects the good control performance, but it is hard to implement. The ETC is easier to implement and energy is saved by the reduction in the rudder rate. In Figure 3, the yaw rate suffers from sudden jumps. To solve this problem, the actuator magnitude and rate limits can be considered as a part of future works.




6. Conclusions
This paper addresses the autopilot design for a surface vehicle in the presence of environmental disturbances and input saturation. The control architecture is proposed based on ESO and the event-triggered strategy. The ESO is employed to reconstruct the yaw rate. Simulation result shows the designed controller can ensure that the actual yaw of the surface ship tracks the reference yaw accurately. A stability analysis based on Lyapunov direct method proves that all signals are bounded and tracking errors converge to a neighborhood of the origin.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the 7th Generation Ultra DeepWater Drilling Unit Innovation Project.