Hongliang Li

Hao Chen
email: chenhao@shmtu.edu.cn

Ning Gao
email: ngao@shmtu.edu.cn

Nadia Aït-Ahmed
email: nadia.ait-ahmed@univ-nantes.fr

Jean-Frederic Charpentier
email: jean-frederic.charpentier@ecole-navale.fr

Mohamed Benbouzid
email: mohamed.benbouzid@univ-brest.fr

Ship Dynamic Positioning Control Based on Active Disturbance Rejection Control

Keywords: Active Disturbance Rejection Control, fixed-point control, straight track control, fal function filter, phase prediction, DP control

Nowadays, Dynamic Positioning (DP) is applied to various tasks such as subsea pipeline laying and the requirements for the positioning performance in marine operations are higher and higher. The main objective of this paper is to design a DP controller based on the Active Disturbance Rejection Control (ADRC) to solve the problems of long response time, large overshoot and low positioning accuracy in ship positioning. Firstly, the mathematical models of the ship and environmental disturbances are established. Secondly, the basic principle of ADRC is described. Meanwhile, stability analysis of the control system is introduced. Thirdly, ADRC is improved by the fal function filter and phase prediction method which solve the problem of dither and phase delay in tracking differentiators. Finally, simulations are carried out to verify the performance of the designed ADRC and the improved ADRC. Several simulation results show that the designed ADRC can realize the fixed-point control of the ship, which effectively solves the problems of long response time, overshoot and positioning accuracy, and compared with the traditional ADRC, the improved ADRC can reduce the error of straight track control, which indicates that ADRC can meet the requirements for the positioning performance and has a strong application value.

Introduction

With the development of modernization and the continuous improvement of living standards, people's demand for coal and oil is increasing. Due to resource shortages caused by the over-exploitation of terrestrial resources, people began to turn their attention to the ocean whose rich resources are of great significance for people to improve economic development and social progress. Dynamic Positioning (DP) is one of the indispensable functions of modern ships, which refers to a technology in which the ship uses its propulsion device to generate thrust to resist external interference and keep the ship in a certain position or sail along a fixed trajectory [START_REF] Wang | Summary of Research on Related Technologies of Ship Dynamic Positioning System[END_REF]. It can be applied to lots of maritime activities, such as convoy, subsea pipeline laying, exploitation of marine resources and hydrological investigation. Research on DP can make ships' control position, heading and trajectory more accurate so that ships can realize safe and flexible operation in crowded offshore production fields without risk of damage to seabed infrastructure and risers [START_REF] Veksler | Dynamic Positioning with Model Predictive Control[END_REF]. A cablelaying vessel with DP system is shown in Figure 1 [3]. Since the 1960s, many scholars all over the world have already conducted studies on the DP. The challenges facing DP control methods are nonlinearity of the model, uncertain disturbances and complexity of the algorithm [START_REF] Fossen | A survey on Nonlinear Ship Control: From Theory to Practice[END_REF][START_REF] Yuh | Design and Control of Autonomous Underwater Robots: A survey[END_REF][START_REF] Smallwood | Model-based dynamic positioning of underwater robotic vehicles: Theory and experiment[END_REF][START_REF] Sørensen | A survey of dynamic positioning control systems[END_REF][START_REF] Mehrzadi | Review of Dynamic Positioning Control in Maritime Microgrid Systems[END_REF]. Therefore, traditional control methods may be not applicable, leading to many advanced methods being proposed to solve the above problems, such as fuzzy control, intelligent algorithm, robust control, Model Predictive Control (MPC), improved Proportional-Integral-Derivative (PID), Sliding Mode Control (SMC) and adaptive control.

Fuzzy Control

Fuzzy control can be used in combination with other algorithms, which plays a great role in parameter problems. In [START_REF] Xu | A fuzzy rule-based PID controller for dynamic positioning of vessels in variable environmental disturbances[END_REF], Xu et al. proposed a fuzzy-PID controller to solve the problems of difficult parameter adjustment and nonlinear ship dynamics. The designed controller considered the positioning error and the low frequency velocity as inputs and the Proportional-Derivative (PD) parameters as outputs of the fuzzy inference. The integral control parameters remained unchanged. The simulation results showed that the proposed fuzzy PID controller could automatically tune the control coefficients according to the positioning accuracy to reduce the percent maximum overshoot, rise time and settling time. In [START_REF] Zhang | Dynamic modeling and adaptive controlling in GPS-intelligent buoy (GIB) systems based on neural-fuzzy networks[END_REF], Zhang et al. proposed an adaptive control based on neural-fuzzy networks to overcome negative effects caused by unavailable velocities and unidentified control parameters. The simulations showed that the designed controllers were capable of making the ship converge to the desired point.

Intelligent Algorithm

Intelligent algorithms include many novel algorithms, such as deep reinforcement learning or Particle Swarm Optimization (PSO). In [START_REF] Øvereng | Dynamic Positioning using Deep Reinforcement Learning[END_REF], Ø vereng et al. introduced a DP control method based on deep reinforcement learning to solve modeling inaccuracies and computational complexity, which were numerically trained through neural networks and used an optimized learning algorithm called Proximal Policy Optimization (PPO) algorithm. The simulation results showed that the proposed method could eliminate body frame errors compared with traditional methods. In [START_REF] Ryalat | Particle Swarm Optimization of a Passivity-Based Controller for Dynamic Positioning of Ships[END_REF], Ryalat et al. proposed an Interconnection and Damping Assignment-Passivity-Based Control (IAD-PBC) controller and used PSO to obtain the gains of the IDA-PBC controller which provides an effective way to tune these gains, thus improving the control performance of DP system. The simulation showed that the states with PSO converged to the desired position and heading angle and had less overshooting.

Robust Control

Robust control means that the system can maintain certain performance under parameter perturbation. In [START_REF] Zhang | Anti-disturbance control for dynamic positioning system of ships with disturbances[END_REF], Zhang et al. studied a robust anti-disturbance strategy to Since the 1960s, many scholars all over the world have already conducted studies on the DP. The challenges facing DP control methods are nonlinearity of the model, uncertain disturbances and complexity of the algorithm [START_REF] Fossen | A survey on Nonlinear Ship Control: From Theory to Practice[END_REF][START_REF] Yuh | Design and Control of Autonomous Underwater Robots: A survey[END_REF][START_REF] Smallwood | Model-based dynamic positioning of underwater robotic vehicles: Theory and experiment[END_REF][START_REF] Sørensen | A survey of dynamic positioning control systems[END_REF][START_REF] Mehrzadi | Review of Dynamic Positioning Control in Maritime Microgrid Systems[END_REF]. Therefore, traditional control methods may be not applicable, leading to many advanced methods being proposed to solve the above problems, such as fuzzy control, intelligent algorithm, robust control, Model Predictive Control (MPC), improved Proportional-Integral-Derivative (PID), Sliding Mode Control (SMC) and adaptive control.

Fuzzy Control

Fuzzy control can be used in combination with other algorithms, which plays a great role in parameter problems. In [START_REF] Xu | A fuzzy rule-based PID controller for dynamic positioning of vessels in variable environmental disturbances[END_REF], Xu et al. proposed a fuzzy-PID controller to solve the problems of difficult parameter adjustment and nonlinear ship dynamics. The designed controller considered the positioning error and the low frequency velocity as inputs and the Proportional-Derivative (PD) parameters as outputs of the fuzzy inference. The integral control parameters remained unchanged. The simulation results showed that the proposed fuzzy PID controller could automatically tune the control coefficients according to the positioning accuracy to reduce the percent maximum overshoot, rise time and settling time. In [START_REF] Zhang | Dynamic modeling and adaptive controlling in GPS-intelligent buoy (GIB) systems based on neural-fuzzy networks[END_REF], Zhang et al. proposed an adaptive control based on neural-fuzzy networks to overcome negative effects caused by unavailable velocities and unidentified control parameters. The simulations showed that the designed controllers were capable of making the ship converge to the desired point.

Intelligent Algorithm

Intelligent algorithms include many novel algorithms, such as deep reinforcement learning or Particle Swarm Optimization (PSO). In [START_REF] Øvereng | Dynamic Positioning using Deep Reinforcement Learning[END_REF], Øvereng et al. introduced a DP control method based on deep reinforcement learning to solve modeling inaccuracies and computational complexity, which were numerically trained through neural networks and used an optimized learning algorithm called Proximal Policy Optimization (PPO) algorithm. The simulation results showed that the proposed method could eliminate body frame errors compared with traditional methods. In [START_REF] Ryalat | Particle Swarm Optimization of a Passivity-Based Controller for Dynamic Positioning of Ships[END_REF], Ryalat et al. proposed an Interconnection and Damping Assignment-Passivity-Based Control (IAD-PBC) controller and used PSO to obtain the gains of the IDA-PBC controller which provides an effective way to tune these gains, thus improving the control performance of DP system. The simulation showed that the states with PSO converged to the desired position and heading angle and had less overshooting.

Robust Control

Robust control means that the system can maintain certain performance under parameter perturbation. In [START_REF] Zhang | Anti-disturbance control for dynamic positioning system of ships with disturbances[END_REF], Zhang et al. studied a robust anti-disturbance strategy to deal with unknown time-varying disturbances of DP system and established a disturbance observer to estimate the disturbances online. The simulation results showed that the strategy could realize the fixed point and yaw control. To solve the problems of unmeasurable speed and disturbances estimation, a Finite-Time State Observer (FTSO) and Finite-Time Feedback Control (FTFC) system based on finite-time control law were proposed in [START_REF] Xia | Finite Time Output Feedback Control for Ship Dynamic Positioning Assisted Mooring Positioning System with Disturbances[END_REF] by Xia et al. The simulation results demonstrated that the proposed FTSO and FTFC could provide a faster convergence speed, better disturbance and parameter perturbation rejection and higher accuracy.

MPC

MPC is a special control related to the model. In [START_REF] Li | Model predictive controller design for ship dynamic positioning system based on state-space equations[END_REF], Li et al. applied the MPC algorithm with a state estimator to ship DP system, aiming at solving the problem of unmeasurable states and measurement noises to improve the control accuracy. The simulation results showed that the MPC controller with the state estimator had the characteristics of fast, stable and accurate responses. In [START_REF] Zheng | Trajectory tracking of autonomous vessels using model predictive control[END_REF], Zheng et al. introduced two different MPC approaches to address the nonlinear horizontal trajectory tracking problem of surface vehicles. The simulation results showed that Nonlinear MPC (NMPC) was better than Linearized MPC (LMPC) in tracking error through a comparison of tracking error of NMPC and LMPC, while NMPC took a much longer time to run the simulation.

5.

Improved PID Due to the limitations of classical PID control, many scholars began to improve PID. A Multi-Variable (MV) PID controller was proposed by Tiwari and Krishnankutty in [START_REF] Tiwari | Comparison of PID and LQR Controllers for Dynamic Positioning of an Oceanographic Research Vessel[END_REF] to overcome the coupled nature of the ship system. The results showed that the force requirement of MV PID controller was minimal only with a bit of sluggishness and a slight overshoot compared with parallel PID and Linear Quadratic Regulator (LQR) controllers. Reference [START_REF] Xu | A fuzzy rule-based PID controller for dynamic positioning of vessels in variable environmental disturbances[END_REF] also introduced a PID controller with fuzzy control, which reduced the percentage maximum overshoot, rise time and settling time.

SMC

SMC is also called variable structure control, which is essentially a special nonlinear control. In [START_REF] Ashrafiuon | Sliding-mode tracking control of surface vessels[END_REF][START_REF] Ashrafiuon | Sliding mode tracking control of surface vessels[END_REF], Ashrafiuon proposed a sliding mode control method using the first-order sliding mode surface and the second-order sliding mode surface to overcome the uncertainty of the system and enhance robustness to interference. The simulation results showed that the proposed methods could successfully make the small boat follow straight-line and circular trajectories. Tannuri et al. proposed a sliding mode control to solve the problem that the positioning performance varied with environmental or loading conditions due to the nonlinearity of the ship in [START_REF] Tannuri | Dynamic positionning systems: An experimental analysis of sliding mode control[END_REF]. This paper made use of systematic experimental tests with a scale model to verify its robustness to variations in its displacement and environmental conditions. In [START_REF] Alattas | Design of a Non-Singular Adaptive Integral-Type Finite Time Tracking Control for Nonlinear Systems with External Disturbances[END_REF], Alattas et al. proposed an adaptive non-singular Fast Terminal Sliding Mode Control (FTSMC) with integral surface for the finite time tracking control of nonlinear systems with external disturbances. The method derived an adaptive parameter-tuning law to tackle the unknown bounded disturbances and alleviate the undesired chattering problem. In [START_REF] Vu | Station-keeping control of a hovering overactuated autonomous underwater vehicle under ocean current effects and model uncertainties in the horizontal plane[END_REF][START_REF] Vu | Robust Position Control of an Over-actuated Underwater Vehicle under Model Uncertainties and Ocean Current Effects Using Dynamic Sliding Mode Surface and Optimal Allocation Control[END_REF], Vu et al. adopted SMC to design controllers and proved the stability of the system based on Lyapunov criteria. The former reference proposed a robust Station-Keeping (SK) control algorithm based on SMC to tackle model uncertainties and ocean current disturbance. Besides, an optimal allocation control was also designed to keep the linear position and Euler angles. The latter reference proposed a motion control based on Dynamic Sliding Mode Control (DSMC) to improve the system robustness with the model uncertainties. Compared with the Least Square (LS) and Quadratic Programming (QP), the QP method could distribute a proper thrust and provide higher stability with smaller steady-state error and stronger robustness. In [START_REF] Piao | Adaptive Backstepping Sliding Mode Dynamic Positioning System for Pod Driven Unmanned Surface Vessel Based on Cerebellar Model Articulation Controller[END_REF], Piao et al. designed a sliding mode controller based on approach law and adaptive backstepping control and added a tracking differentiator to eliminate the large chattering. Compared with traditional sliding mode control, the designed controller had a better performance on dynamic positioning under certain disturbances.

7.

Adaptive Control Hu et al. introduced an adaptive control approach based on Dynamic Surface Control (DSC) and Minimal Learning Parameter (MLP) techniques in [START_REF] Hu | Robust Adaptive Fuzzy Design for Ship Linear-tracking Control with Input Saturation[END_REF]. The proposed complete system was developed to overcome the negative effect of input saturation constraints and the singularity problem of the controller. The simulation results showed that the tracking error could converge to a small interval. In [START_REF] Zhao | Adaptive nonlinear controller design for course-tracking of ships based on dynamic surface control and backstepping[END_REF], Zhao et al. designed an adaptive nonlinear course controller based on the backstepping control method and dynamic surface control to solve the problems of parameter uncertainties and completely unknown control gain, but also to overcome the "explosion of terms" of the traditional backstepping control method. The simulation results showed that the method could make the error converge approximately to 0, which was proved effective. In [START_REF] Li | Robust adaptive neural network control for dynamic positioning of marine vessels with prescribed performance under model uncertainties and input saturation[END_REF], Li et al. proposed a Robust Adaptive Neural Network Control (RANNC) scheme for DP of marine vessels to overcome model uncertainties, external disturbances and input saturation. This scheme combined adaptive control and neural network, which tackled the changes of inertial matrix and damping matrix.

In the above studies, although the requirement of DP ships can be met, there are common problems of long response time, large overshoot or low positioning accuracy which may cause position deviation in laying subsea cables or oil extraction, especially when the environmental interference is large; it may cause the motion of ships to be out of control. At the same time, the ship modeling is complex and needs nonlinear approaches. Therefore, it is difficult to define their accurate kinematic and kinetic models. The Active Disturbance Rejection Control (ADRC) can overcome the dependence on the precise model and has the advantage of fast speed, and its estimation and compensation of the total disturbance can avoid the side effect of error integral feedback, which makes the signal achieve complete tracking to reduce the error and increase positioning accuracy without large overshooting. Therefore, this paper uses ADRC to realize the fixed-point control and the straight track control of the ship. For the setting of the trajectory, this paper uses a parametric equation to simplify complex problems. The objective of the fixed-point control is to keep the ships constantly near of the set position and stay berthed after reaching the set position. In the fixed-point control, through comparison between ADRC and backstepping method, it can be seen that ADRC does have a fast response speed. The objective of the straight track control is to make the ship sail along the set trajectory and minimize the error, while keeping the ship speed at the set speed. Due to the possible problems of dither and phase delay of classical ADRC, this paper is to use fal function filter and phase prediction to improve these limitations. The fal function has fast convergence and can filter out dithering well. This paper applies the fal function filter to the tracking differentiator to solve the dither problem. The phase prediction based on mathematical differential allows good phase compensation and is easy to realize. This paper combines the two methods with ADRC into dynamic positioning to improve the control performance. The simulation results demonstrate the effectiveness and practicality of the proposed controllers and their robustness to external disturbances.

The paper is constructed as follows. Firstly, Section 2 introduces the ship model including coordinate systems and mathematical models of the ship. Next, Section 3 describes the models of environmental disturbances including wind and waves. Then, Section 4 introduces the DP controller based on the basic principle of ADRC and describes the stability analysis of the control system. After that, Section 5 analyzes the problems of the controller and improves the controller by fal function filter and phase prediction. Subsequently, Section 6 presents the simulations that the designed controller is used to realize the fixed-point control and yaw control compared with backstepping method, and the effectiveness of the Improved Active Disturbance Rejection Control (IADRC) is verified under the condition of straight track control. Finally, Section 7 provides conclusions of this paper and introduces suggestions for future research.

Ship Model 2.1. Assumptions

To simplify the problem, this work is based on the following assumptions:

(1) The motion of the ship in roll, pitch and heave is neglected.

(2) The motion of the ship is regarded as a plane motion.

(3) The water area of the ship is wide enough, and the hull draft is unchanged during the movement.

Coordinate System

The first step in solving the DP problem is to describe the motion of the ship in an appropriate coordinate system. Generally speaking, two coordinate systems are used to describe the motion of the ship, including the North-East-Down (NED) coordinate system (X 0 , Y 0 , Z 0 ) and the body-fixed coordinate system (X, Y, Z). The two coordinate systems are shown in Figure 2 under the condition of straight track control. Finally, Section 7 provides conclusions of this paper and introduces suggestions for future research.

Ship Model

Assumptions

To simplify the problem, this work is based on the following assumptions:

(1) The motion of the ship in roll, pitch and heave is neglected.

(2) The motion of the ship is regarded as a plane motion.

(3) The water area of the ship is wide enough, and the hull draft is unchanged during the movement.

Coordinate System

The first step in solving the DP problem is to describe the motion of the ship in an appropriate coordinate system. Generally speaking, two coordinate systems are used to describe the motion of the ship, including the North-East-Down (NED) coordinate system (X 0 , Y 0 , Z 0 ) and the body-fixed coordinate system (X, Y, Z). The two coordinate systems are shown in Figure 2 The NED coordinate system is a coordinate system based on the earth. Its origin O can be selected at any point on the surface of the earth, and its three axes point to north, east and the direction of gravity, respectively. Once determined, the three-axis directions are always maintained constant.

The origin of the body-fixed coordinate system is located at the center of gravity of the ship: its horizontal axis X is parallel to the hull and points to the forward direction; the longitudinal axis Y is perpendicular to the hull and points to the starboard side of the hull; the axis Z is same as the direction of gravity and points to the direction of the keel.

Mathematical Model of the Ship

The mathematical models of ship motion include the kinematic model and the dynamic model. The actual motion of the ship is extremely complex, and generally has six Degrees of Freedom (DOF). The mathematical model proposed by T.I. Fossen in [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF] is used to model the ship motion in this study.

Kinematic Model

For marine ships, six DOF models are often used to describe the motion of ships. These six DOF are named surge, sway, heave, roll, pitch and yaw. The vector η consists The NED coordinate system is a coordinate system based on the earth. Its origin O can be selected at any point on the surface of the earth, and its three axes point to north, east and the direction of gravity, respectively. Once determined, the three-axis directions are always maintained constant.

The origin of the body-fixed coordinate system is located at the center of gravity of the ship: its horizontal axis X is parallel to the hull and points to the forward direction; the longitudinal axis Y is perpendicular to the hull and points to the starboard side of the hull; the axis Z is same as the direction of gravity and points to the direction of the keel.

Mathematical Model of the Ship

The mathematical models of ship motion include the kinematic model and the dynamic model. The actual motion of the ship is extremely complex, and generally has six Degrees of Freedom (DOF). The mathematical model proposed by T.I. Fossen in [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF] is used to model the ship motion in this study.

Kinematic Model

For marine ships, six DOF models are often used to describe the motion of ships. These six DOF are named surge, sway, heave, roll, pitch and yaw. The vector η consists of the position components and angle component of the marine ships. The vector ν consists of the linear speed components and angular speed component of the marine ships. The vectors can be described as follows:

η = [x y z ϕ θ ψ] T ,η 1 = [x y z] T , η 2 = [ϕ θ ψ] T ; ν = [u v w p q r] T , ν 1 = [u v w] T , ν 2 = [p q r] T .
The kinematic model of ships is shown in Equation (1):

.

η 1 = J 1 (η)ν 1 . η 2 = J 2 (η)ν 2 (1)
The J 1 (η) and J 2 (η) are rotation matrices:

J 1 (η) =   cos θ cos ψ cos ψ sin θ sin ϕ -cos ϕ sin ψ sin ψ sin ϕ + cos ψ cos ϕ sin θ cos θ sin ψ sin ψ sin θ sin ϕ + cos ϕ cos ψ -cos ψ sin ϕ + sin ψ cos ϕ sin θ -sin θ cos θ sin ϕ cos θ cos ϕ   , J 2 (η) =   1 sin ϕ tan θ cos ϕ tan θ 0 cos ϕ -sin ϕ 0 sin ϕ/ cos θ cos ϕ/ cos θ   .
So Equation ( 1) can be expressed as Equation ( 2

): . η 1 . η 2 = J 1 (η) 0 3×3 0 3×3 J 2 (η) ν 1 ν 2 (2)
For general surface vessels, the six DOF can be simplified to three DOF including surge, sway and yaw because vessels always travel at a low speed. In this case, ϕ and θ are very small and can be neglected. It is a good approximate expression for most conventional ships and offshore platforms [START_REF] Du | A robust adaptive neural networks controller for maritime dynamic positioning system[END_REF][START_REF] Hu | Adaptive fuzzy controller design for dynamic positioning system of vessels[END_REF]. Therefore, the rotation matrices can be expressed as:

J 1 (η) ≈   cos ψ -sin ψ 0 sin ψ cos ψ 0 0 0 1   , J 2 (η) ≈ I 3×3
Therefore, the kinematic model of ships with three DOF can be rewritten as:

  . x . y . ψ   =   cos ψ -sin ψ 0 sin ψ cos ψ 0 0 0 1     u v r   (3)
That is:

. η = J(ψ)ν (4)
where

J(ψ) =   cos ψ -sin ψ 0 sin ψ cos ψ 0 0 0 1   .
In the kinematic model ( 4), the simplified vector η = x y ψ T denotes the ship position (x, y) and yaw component ψ in the earth fixed coordinate system. The simplified vector ν = u v r T denotes the ship linear speed (surge: u, sway: v) and angular speed r (yaw component). J(ψ) is the rotation matrix between earth fixed frame and body-fixed frame, which is nonsingular for all ψ and has the following property: J -1 (ψ) = J T (ψ) [START_REF] Peng | A robust unscented Kalman filter and its application in estimating dynamic positioning ship motion states[END_REF].

Dynamic Model

The rigid body dynamic equation can be expressed in the following form:

M . ν + C(ν)ν + D(ν)ν = τ + a (5)
In the dynamic model ( 5), M is the total inertia matrix including additional mass, which can be expressed as:

M =   m -X . u 0 0 0 m -Y . v mx g -Y . r 0 mx g -N . v I zz -N . r   (6)
The matrix M can be expressed as M = M A +M RB . M A is the inertia matrix of the hydrodynamic system and M RB is the inertia matrix of the rigid body system (platform or vessel).

C(ν) is the Coriolis-centripetal force matrix. Generally, in the case of low-speed motion, the Coriolis-centripetal force matrix C(ν) is very small and can be neglected, namely C(ν) ≈ 0. D(ν) is the damping matrix which the surge mode can be decoupled from steering modes (sway and yaw). Therefore, the linearized damping matrix can be written:

D(ν) =   -X u 0 0 0 -Y v -Y r 0 -N v -N r   (7)
The damping coefficients can be considered linear in the case of low-speed operation which means that N v = Y r such that D T = D [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF].

The vector τ = [Fx Fy Fn] represents the total thrust forces and moment generated by thrusters. The term a is used to represent the external environmental disturbance forces and moments.

In summary, the three DOF motion mathematical model of the ship can be expressed as a combination of the kinematic model and the dynamic model, as shown in Equation ( 8) [START_REF] Ngongi | Design of generalised predictive controller for dynamic positioning system of surface ships[END_REF]:

. η = J(ψ)ν M . ν + D(ν)ν = τ + a (8)

Environmental Disturbances including Wind and Wave

Some experts and scholars have conducted studies on environmental disturbances. Reference [START_REF] Gao | RANS Prediction of Wave-Induced Ship Motions, and Steady Wave Forces and Moments in Regular Waves[END_REF] predicted for KVLCC2 performing straight ahead motion in head regular waves by using the expanded RANS solver based on OpenFOAM. Reference [START_REF] Borkowski | Numerical Modeling of Wave Disturbances in the Process of Ship Movement Control[END_REF] proposed a numerical model of sea wave generation and a forming filter based on the wave spectrum. This paper adopted external disturbances models from Thor I. Fossen. For control system design it is common to assume the principle of superposition when considering wind and wave disturbances. For most marine control applications this is a good approximation [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF]. Environmental disturbances vector a can be expressed as: a = a wind + a wave [START_REF] Xu | A fuzzy rule-based PID controller for dynamic positioning of vessels in variable environmental disturbances[END_REF] a wind means the wind forces and moment vector. a wave means the waves forces and moment vector.

Wind Disturbance

The wind speed V T at any height can be expressed as:

V T = V * ( h 10 ) 1/7 ( 10
)
where V is the average wind speed and h is height.

Relative wind speed V r can be expressed as:

V r = V 2 T +V 2 S + 2V T V S cos(ϕ wind -ϕ ship + γ) (11)
where V s is the speed of ship, ϕ wind is the wind direction, ϕ ship is the ship direction and γ is the drift angle. The generalized wind force vector is a wind = [X wind , Y wind , N wind ]. The wind distur- bance models in this paper proposed by Isherwood [START_REF] Isherwood | Wind resistance of merchant ships[END_REF] are expressed in Equation ( 12):

   X wind = 0.5C X ρ α V 2 r A T Y wind = 0.5C Y ρ α V 2 r A L N wind = 0.5C N ρ α V 2 r LA L ( 12
)
where C X and C Y are the empirical force coefficients, C N is a moment coefficient, which can be obtained by regression analysis. ρ α is the density of air, A T and A L are the transverse and lateral projected areas, and L is the overall length of the ship.

Wave Disturbance

The wave force can be divided into first-order wave force and second-order wave force. The first-order wave force is often ignored due to its small effect. The second-order wave force can be considered to be generated by the wave spectra. This paper adopts Pierson-Moskowitz Spectrum (PM Spectrum) to model the wave, which is shown as:

S(ω i ) = 0.0081g 2 ω -5 i exp(- 3.11 H 2 s ω -4 i ) ( 13
)
where g is gravitational acceleration, ω i is the wave frequency of the i-th wave and H s is wave height. The waves force vector is a wave = [X wave , Y wave , N wave ]. The wave disturbance can be expressed as follows:

   X wave = -0.5ρLµ 2 cos χC XW Y wave = 0.5ρLµ 2 sin χC YW N wave = 0.5ρL 2 µ 2 sin χC NW [START_REF] Xia | Finite Time Output Feedback Control for Ship Dynamic Positioning Assisted Mooring Positioning System with Disturbances[END_REF] where χ is encounter angle and µ is wave amplitude. C XW , C YW and C NW are coefficients obtained by regression analysis. In this paper, the fourth-level sea state is taken into account as environmental disturbances, which have 2 m wave height and 10 m/s average wind speed.

ADRC Ship DP Controller

The control objective of the system modeling is to find a control law such that the actual position and trajectory can converge to a desired position and trajectory.

Basic Principle of ADRC

Deriving Equation ( 8) allows to establish Equation (15):

.. η = . Jν -JM -1 Dν + JM -1 a + JM -1 τ (15)
By setting x 1 = η and x 2 = .

η, a state space model of the ship can be obtained from Equation [START_REF] Li | Model predictive controller design for ship dynamic positioning system based on state-space equations[END_REF] as follows:

   . x 1 = x 2 . x 2 = f(x 1 ,x 2 ,t) + ω(t) + bu y = x 1 (16)
where u = JM -1 τ is the control quantity, ω(t) = JM -1 a is the unknown disturbance and f(x 1 ,x 2 ,t) = .

Jν-JM -1 Dν is the known disturbance. y is the output and b is the control gain. Due to the good decoupling performance of ADRC, a single dimension is considered in the following part to simplify analysis. The control analysis of this system is described as follows:

Step 1: Arranging the transition process.

In classical control, . η = η(t)-η(t-T)

T is generally used to generate differential signals. The smaller T is, the closer the output is to the differential signal. However, if T is too small, this process may result in the amplification of noise effects. Therefore, it is necessary to take a suitable method to extract the differential signal.

For the system shown below, it is hoped that η 1 can be asymptotically stable at the target η 0 and η 2 can be asymptotically stable at the differential of the target. The algorithm formula is as follows:

.

η 1 = η 1 + h * η 2 . η 2 = η 2 + h * Γ (17)
Γ is a control law which can be taken as a sign function or a sat function (a saturation function). However, the two functions will produce high-frequency noise when the system enters a steady state. Hence, this paper adopts fst function to design the fastest discrete feedback system. The specific expression of fst function is [START_REF] Wang | Design of auto disturbance rejection controller for train traction control system based on artificial bee colony algorithm[END_REF]:

                           d = r 0 * h 0 d 0 = d * h 0 y = η 1 +h 0 * η 2 α 0 = d 2 + 8r 0 |y| α = η 2 + (α 0 -d)sign(y) 2 , |y| > d 0 η 2 + y h 0 , |y | ≤ d 0 fst = -r 0 * α d , |α | ≤ d -r 0 * sign(α) , |α| > d (18)
Equation ( 18) is the fastest synthesis function, which can be obtained by the 'Isochronic Region method'. Due to the complexity of derivation process, the specific derivation process and variable meanings can be learned in reference [START_REF] Han | Active Disturbance Rejection Control Technique-the Technique for Estimating and Compensating the Uncertainties[END_REF] written by Prof. Han. In the Equation ( 18), r 0 is the factor of tracking speed and h 0 is filter factor that represents filter effect. For the speed factor r 0 , the value should not be too large or too small. If the value is too small, the tracking speed is too slow, while if the value is too large, the differential signal is poor when the given signal will be noisy.

Therefore, transition process can be designed as Equation [START_REF] Ashrafiuon | Sliding mode tracking control of surface vessels[END_REF] and then the tracking signal η 1 and differential signal η 2 can be obtained according to the given input signal η 0 . This module can be called Tracking Differentiator (TD). The specific algorithm formula is shown in Equation ( 19):

.

η 1 = η 1 + h * η 2 . η 2 = η 2 + h * fst(η 1 -η 0 , η 2 , r 0 , h 0 (19)
where fst(•) is the modified fastest synthesis function. η 1 is the tracking signal and the η 2 is the differential signal. h is the sample step time. For the step h, the value should not be too small, otherwise the differential performance will not be guaranteed. The value of the step h can generally be chosen in the (0.001, 0.1) range.

Step 2: Estimating the state of each order of the controlled system. The total disturbance f(x 1 , x 2 , t)+ω(t) in the second-order ship system ( 16) is expanded into a new state variable f. Then state equation of the system can be expressed as:

       . x 1 = x 2 . x 2 = f + bu . f = ω 0 (t) y = x 1 (20)
For such a system, the state observer can be expressed as:

         e = η1 -η . η1 = η2 -β 01 * e . η2 = f -β 02 * |e| 0.5 sign(e) + b * u . f = -β 03 * |e| 0.25 sign(e) (21)
where η1 , η2 , f are the estimations of x 1 , x 2 , f, respectively. β 01 , β 02 , β 03 are adjustable parameters whose values should be greater than 0. The three parameters determine the effect of estimating the disturbances and have a mutually restrictive relationship, the ideal estimation effect of which can be obtained by coordination action. Hence, as long as the appropriate β 01 , β 02 and β 03 are selected, the system can estimate the state variables x 1 , x 2 and f, namely η1 → x 1 , η2 → x 2 and f → f . In order to avoid the phenomenon of high-frequency flutter, |e| 0.5 sign(e) and |e| 0.25 sign(e) need to be transformed into a nonlinear function which is continuous around the origin. The fal function was therefore proposed by Han in [START_REF] Han | From PID to Active Disturbance Rejectio Control[END_REF]. Its expression is shown as follows:

fal(e, α, δ) = |e| α sign(e) , |e|> δ e δ 1-α , |e| ≤ δ ( 22
)
Fal function is a nonlinear function that has fast convergence which can be referred to in Section 5.2.1. δ is a small number greater than 0 and α is a parameter that determines the shape of the fal function. This parameter is chosen in (0, 1) range.

According to the input u and output η of the system, the observer can be designed as shown in Equation ( 23):

         e = η1 -η . η1 = η2 -β 01 * e . η2 = f -β 02 * fal(e, α 1 , δ)+b * u . f = -β 03 * fal(e, α 2 , δ) (23)
where b is control coefficient. β 01 , β 02 , β 03 are adjustable parameters greater than 0.

This observer is called Extended State Observer (ESO) and the control effect of the system depends on the observation of the total disturbances by ESO.

Step 3: Non-linear combination and the compensation of total disturbance. The control variable u of classical control is usually a linear combination of the error related to present (e(t)), past ( t 0 e(t)dt) and future ( de(t) dt ). However, such a combination is not necessarily applicable in a large number of engineering practices. Therefore, ADRC uses the fal function mentioned before to change the linear combination into the nonlinear combination. The non-linear combination is called Non-linear State Error Feedback (NLSEF) and its expression is shown in Equation ( 24):

       e 1 = η 1 -η1 e 2 = η 2 -η2 u 0 = k 1 * fal(e 1 , α n1 , δ)+k 2 * fal(e 2 , α n2 , δ) u = u 0 -f/b ( 24
)
where k 1 and k 2 are adjustable parameters. As long as the parameters of fal function and k 1 and k 2 are selected appropriately, it can ensure that the controller may have good robustness and adaptability.

Stability Analysis of Control System

Rewrite the state space form of the second-order controlled ship system in this paper as follows:

   . x 1 = x 2 . x 2 = f + bu y = x 1 (25)
To simplify the analysis of the problem, this paper uses the following assumptions:

(1) Set z 1 = η1 , z 2 = η2 and z 3 = f. The three variables are the estimations of tracking signal, differential signal and total disturbances, respectively. (2) Let the input be 0, then the output of TD is also 0.

(3) Change the Nonlinear State Error Feedback (NLSEF) into Linear State Error Feedback (LSEF):

u = k 1 z 1 + k 2 z 2 - z 3 b (26)
(4) Select nonlinear function ϕ(x) which satisfies that if x = 0, then xϕ(x)> 0. Define ϕ ∈ F(ς 1 , ς 2 ) to satisfy the following conditions: 1 ϕ(0)= 0; 2 ς 1 x 2 < xϕ(x)< ς 2 x 2 (∀x = 0) and then xϕ(x)> 0, namely ϕ ∈ F(0, ∞).

   . z 1 = z 2 -β 01 ϕ(z 1 -y) . z 2 = z 3 -β 02 ϕ(z 1 -y) + bu . z 3 = -ϕ(z 1 -y) (27)
(5) Let the control object be a linear time-invariant object:

..

x = -a 2 xa 1 .

x

+ bu y = x (28)
Substituting Equation [START_REF] Zhao | Adaptive nonlinear controller design for course-tracking of ships based on dynamic surface control and backstepping[END_REF] and Equation [START_REF] Du | Adaptive Robust Output Feedback Control for a Marine Dynamic Positioning System Based on a High-Gain Observer[END_REF] into Equation ( 25) can get:

. X = A 11 X + A 12 Z + a 13 z 3 y = x 1 [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF] where

X = [x 1 , x 2 ] T , Z = [z 1 , z 2 ] T , a 13 = [0, -1] T , A 11 = 0 1 -a 2 -a 1 , A 12 = 0 0 bk 1 bk 2 .
Substitute Equation [START_REF] Zhao | Adaptive nonlinear controller design for course-tracking of ships based on dynamic surface control and backstepping[END_REF] into Equation ( 27) to get:

   . Z = A 22 Z + b 2 ζ . z 3 = β 03 ζ ζ = -ϕ(z 1 -y) (30) where b 2 = [β 01 , β 02 ] T , A 22 = 0 1 bk 1 bk 2 .
Combine Equation [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF] and Equation ( 30) can get:

             . X = A 11 X + A 12 Z + a 13 z 3 . Z = A 22 Z + b 2 ζ . z 3 = β 03 ζ ζ = -ϕ(θ) θ = c T 1 X + c T 2 Z (31)
where

c 1 = [ -1, 0] T , c 2 = [1, 0] T .
Let Y = A 11 X + a 13 z 3 , z 3 = z 3 /β 03 and Equation ( 32) can be available after arranging the equation:

             . Y = A 11 Y + A 11 A 12 Z + a 13 β 03 ζ . Z = A 22 Z + b 2 ζ . z 3 = β 03 ζ ζ = -ϕ(θ) θ = c T 1 A -1 11 Y + c T 2 Z -c T 1 A -1 11 a 13 β 03 z 3 (32)
That is:

                 . Y . Z = A 11 A 11 A 12 0 A 22 Y Z + a 13 β 03 b 2 ζ . z 3 = β 03 ζ ζ = -ϕ(θ) θ = [c T 1 A -1 11 c T 2 ] Y Z -c T 1 A -1 11 a 13 z 3 (33)
Equation ( 33) conforms to a standard form of the first critical case of absolute stability, which is:

       . x = Ax + bu . ξ = u u = -ϕ(y) y = c T x + ρξ (34)
where

A = A 11 A 11 A 12 0 A 22 , b = a 13 β 03 b 2 , c T = [c T 1 A -1 11 c T 2 , ρ = -c T 1 A -1 11 a 13 .
Theorem 1. The necessary condition for the absolute stability of the zero solution of F(0, ∞) in System [START_REF] Gao | RANS Prediction of Wave-Induced Ship Motions, and Steady Wave Forces and Moments in Regular Waves[END_REF] is Re[λ(A)] < 0, ρ > 0, [START_REF] Xie | Absolute Stability Theory and Application[END_REF].

We can construct a Lyapunov function:

V(x, y) = x T Qx+ (y -c T x) 2 + β y 0 ϕ(y)dy (35)
If Q is a positive definite matrix and , β ≥ 0, V(x, y) is a positive definite function.

Find the derivative of V(x, y):

- . V(x, y) = x T -QA -A T Q x + 2b T Q -2 ρc T -βc T A xϕ(y) +(βρ + βc T b)ϕ 2 (y) + 2 ρyϕ(y) (36)
Since yϕ(y) > 0, ρ > 0, as long as ≥ 0 and the first three terms of Equation ( 36) are the positive definite quadratic form, then .

V must be negative definite, thus ensuring the global asymptotic stability of the system. Therefore, as long as the parameters are properly selected to satisfy the relevant conditions, the stability can be ensured.

Structure of ADRC

According to the above analysis, ADRC method of the DP ship can be obtained based on Equations ( 19), ( 23) and [START_REF] Piao | Adaptive Backstepping Sliding Mode Dynamic Positioning System for Pod Driven Unmanned Surface Vessel Based on Cerebellar Model Articulation Controller[END_REF]. ADRC is composed of three parts [START_REF] Ye | A design of dredger cutter motor synchronous speed control system based on ADRC[END_REF]: TD, ESO and NLSEF. The structure of ADRC is shown in Figure 3 [START_REF] Hu | Research on the course control of USV based on improved ADRC[END_REF]. ADRC is a nonlinear control that does not rely on precise mathematical models [START_REF] Wu | On transitioning from PID to ADRC in thermal power plants[END_REF]. It promotes and enriches the error-based idea of classical control and has strong advantages, especially in uncertain models and the ocean environment [START_REF] Wang | Improved ADRC control strategy in FPSO dynamic positioning control application[END_REF], which has fast speed, strong anti-disturbance ability and good nonlinear control performance.

According to the above analysis, ADRC method of the DP ship can be obtained based on Equations ( 19), ( 23) and [START_REF] Piao | Adaptive Backstepping Sliding Mode Dynamic Positioning System for Pod Driven Unmanned Surface Vessel Based on Cerebellar Model Articulation Controller[END_REF]. ADRC is composed of three parts [START_REF] Ye | A design of dredger cutter motor synchronous speed control system based on ADRC[END_REF]: TD, ESO and NLSEF. The structure of ADRC is shown in Figure 3 [START_REF] Hu | Research on the course control of USV based on improved ADRC[END_REF]. ADRC is a nonlinear control that does not rely on precise mathematical models [START_REF] Wu | On transitioning from PID to ADRC in thermal power plants[END_REF]. It promotes and enriches the errorbased idea of classical control and has strong advantages, especially in uncertain models and the ocean environment [START_REF] Wang | Improved ADRC control strategy in FPSO dynamic positioning control application[END_REF], which has fast speed, strong anti-disturbance ability and good nonlinear control performance.

Improvement of ADRC Controller

Analysis of TD Behavior and Limitations

Since ADRC has good decoupling performance and different working conditions have different requirements on the ship's trajectory, this section only considers the control effect of single control channel to conduct simplified theoretical analysis, which is shown in Figure 4. This paper assumes that the input reference signal is a sinusoidal signal ϕ(t)=sin(0.1t) to illustrate the existing problems. The simulation result on TD is shown in Figure 5. The red line is the input signal ϕ(t), and the blue line is the tracking signal ϕ 1 (t). It is found that the tracking signal ϕ 1 (t) has a phase delay and dither, which has a great impact on the track control and may cause the tracking error to be non-zero all the time. This error may lead to a decrease in the accuracy of the ship's DP control, possibly resulting in position deviations in some marine operations such as oil extraction and cable laying. Therefore, this paper adopts an effective method to solve the problems of dither and phase delay to reduce the error of straight track control.

Improvement of ADRC Controller

Analysis of TD Behavior and Limitations

Since ADRC has good decoupling performance and different working conditions have different requirements on the ship's trajectory, this section only considers the control effect of single control channel to conduct simplified theoretical analysis, which is shown in Figure 4.

NLSEF. The structure of ADRC is shown in Figure 3 [START_REF] Hu | Research on the course control of USV based on improved ADRC[END_REF]. ADRC is a nonlinear control that does not rely on precise mathematical models [START_REF] Wu | On transitioning from PID to ADRC in thermal power plants[END_REF]. It promotes and enriches the errorbased idea of classical control and has strong advantages, especially in uncertain models and the ocean environment [START_REF] Wang | Improved ADRC control strategy in FPSO dynamic positioning control application[END_REF], which has fast speed, strong anti-disturbance ability and good nonlinear control performance.

Improvement of ADRC Controller

Analysis of TD Behavior and Limitations

Since ADRC has good decoupling performance and different working conditions have different requirements on the ship's trajectory, this section only considers the control effect of single control channel to conduct simplified theoretical analysis, which is shown in Figure 4. This paper assumes that the input reference signal is a sinusoidal signal ϕ(t)=sin(0.1t) to illustrate the existing problems. The simulation result on TD is shown in Figure 5. The red line is the input signal ϕ(t), and the blue line is the tracking signal ϕ 1 (t). It is found that the tracking signal ϕ 1 (t) has a phase delay and dither, which has a great impact on the track control and may cause the tracking error to be non-zero all the time. This error may lead to a decrease in the accuracy of the ship's DP control, possibly resulting in position deviations in some marine operations such as oil extraction and cable laying. Therefore, this paper adopts an effective method to solve the problems of dither and phase delay to reduce the error of straight track control. This paper assumes that the input reference signal is a sinusoidal signal φ(t) = sin(0 .1t) to illustrate the existing problems. The simulation result on TD is shown in Figure 5. The red line is the input signal φ(t), and the blue line is the tracking signal φ 1 (t). It is found that the tracking signal φ 1 (t) has a phase delay and dither, which has a great impact on the track control and may cause the tracking error to be non-zero all the time. This error may lead to a decrease in the accuracy of the ship's DP control, possibly resulting in position deviations in some marine operations such as oil extraction and cable laying. Therefore, this paper adopts an effective method to solve the problems of dither and phase delay to reduce the error of straight track control.

Improvement Methods

Fal Function Filter

To solve the dithering problem, this paper adopts the fal function filter to eliminate dither. The fal function has good filtering performance and can also solve the dithering problem well. The characteristic analysis of the fal function filter is given below.

Improvement Methods

Fal Function Filter

To solve the dithering problem, this paper adopts the fal function filter to eliminate dither. The fal function has good filtering performance and can also solve the dithering problem well. The characteristic analysis of the fal function filter is given below.

The following system is available:

. x = w(x, t) + s [START_REF] Wang | Design of auto disturbance rejection controller for train traction control system based on artificial bee colony algorithm[END_REF] where w(x, t) represents environmental disturbances and is less than w d which is maximum value of external disturbances. The form s is the control variable. This paper takes the feedback s of state x as s = s(x). Let s(x) = -k * fal(x, α, δ) and then the fal function can be rewritten as follows:

fal(x, α, δ) = |x| α sign(x) , |x|> δ x δ 1-α , |x| ≤ δ , δ > 0 ( 38
)
When |x|> δ , system (37) becomes:

. x = w(x, t) -k|x| α sign(x)

Multiply both sides by 2x and Equation ( 40) can be obtained as follows:

d x 2 dt ≤ -2k|x|( |x| α - w d k ) < 0 (40)
x = w(x, t) -k x δ 1-α (41)
Multiply both sides by 2x and Equation ( 42) can be obtained as follows:

d x 2 dt ≤ -2x 2 k δ 1-α - w d x < 0 (42) When |x (t)| > |w d |δ 1-α k
, control variable s = -k|x| α sign(x) can control the steady error caused by w(x, t) within the range of |x| ≤ |w d |δ 1-α k . The feedback structure of fal function filter can be shown in Figure 6 [START_REF] Lei | Noise suppression by ADRC for ship positioning and yaw control[END_REF]. The formula of this filter can be expressed as [START_REF] Dong | An active disturbance rejection controller design for ball and plate system[END_REF]:

   . q = k * fal(e, α, δ) e = p -q p 0 = q ( 43
)
where k is the scale factor. p is the input signal of the filter and p 0 is the output of the filter. The feedback structure of fal function filter can be shown in Figure 6 [START_REF] Lei | Noise suppression by ADRC for ship positioning and yaw control[END_REF]. The formula of this filter can be expressed as [START_REF] Dong | An active disturbance rejection controller design for ball and plate system[END_REF]:

{ q̇=k * fal(e,α,δ) e=p -q p 0 =q ( 43
)
where k is the scale factor. p is the input signal of the filter and p 0 is the output of the filter. This paper analyzes the filter shown in Equation ( 43). When |e|>δ, nonlinear feedback k * |e| α sign(e) makes the system state rapidly approach the input signal p, so that the error e approaches δ.

When |e|≤δ, q̇=k * fal(e,α,δ)=k * e/δ 1-α . Set k f =k/δ 1-α . In this case, a transfer func- This paper analyzes the filter shown in Equation ( 43).

When |e| > δ , nonlinear feedback k * |e| α sign(e) makes the system state rapidly approach the input signal p, so that the error e approaches δ.

When |e| ≤ δ, . q = k * fal(e, α, δ) = k * e/δ 1-α . Set k f = k/δ 1-α . In this case, a transfer function from input to output can be obtained in Equation ( 44):

p 0 p = k f /s 1 + k f /s = k f s + k f (44)
Equation ( 44) is actually a low-pass filter. When k decreases, the bandwidth becomes narrower and the filtering effect may be better.

According to the above analysis, a fal function filter is used at the output of the TD. As shown in Figure 7, the blue signal is the filtered signal. It can be clearly seen that compared with Figure 5, the original dither signal is turned into a smooth signal by the operation of the fal function filter, which solves the dithering problem. This paper analyzes the filter shown in Equation [START_REF] Wu | On transitioning from PID to ADRC in thermal power plants[END_REF]. When |e|>δ, nonlinear feedback k * |e| α sign(e) makes the system state rapidly approach the input signal p, so that the error e approaches δ.

When |e|≤δ, q̇=k * fal(e,α,δ)=k * e/δ 1-α . Set k f =k/δ 1-α . In this case, a transfer function from input to output can be obtained in Equation ( 44):

p 0 p = k f /s 1+k f /s = k f s+k f (44)
Equation ( 44) is actually a low-pass filter. When k decreases, the bandwidth becomes narrower and the filtering effect may be better.

According to the above analysis, a fal function filter is used at the output of the TD. As shown in Figure 7, the blue signal is the filtered signal. It can be clearly seen that compared with Figure 5, the original dither signal is turned into a smooth signal by the operation of the fal function filter, which solves the dithering problem.

Phase Prediction

As is shown in Figure 8, the phase prediction method is used in this paper to solve the phase delay problem.

Phase Prediction

As is shown in Figure 8, the phase prediction method is used in this paper to solve the phase delay problem. In Figure 8, the red signal is the input curve ϕ(t) and the blue one is the tracking curve ϕ 1 (t).

According to Figure 8, the new tracking curve ϕ 1 ' (t) can be expressed as Equation ( 45):

ϕ 1 ' (t)=ϕ 1 (t)+Δϕ (45)
where Δϕ=ϕ(t 2 )ϕ(t 1 ). In Figure 8, the red signal is the input curve φ(t) and the blue one is the tracking curve φ 1 (t).

According to Figure 8, the new tracking curve φ 1 (t) can be expressed as Equation ( 45):

φ 1 (t) = φ 1 (t) + ∆φ (45)
where ∆φ = φ(t 2 )φ(t 1 ).

Considering that φ 2 (t) is the differential signal of the input signal φ(t), Equation ( 46) can be obtained according to the definition of differential:

φ 2 (t) = lim t 2 →t 1 φ(t 2 ) -φ(t 1 ) t 2 -t 1 ( 46
)
When the time interval is small, Equation ( 46) can be approximately rewritten as:

∆φ = (t 2 -t 1 )φ 2 (t) (47)
Then Equation ( 45) can be rewritten as:

φ 1 (t) = φ 1 (t) + (t 2 -t 1 φ 2 (t) ( 48
)
where t 2t 1 = λ * h. h is the step time in the tracking differentiator parameters, λ = 2 * r 1 = 2 * h 0 h is the compensation coefficient and h 0 is the filter factor which represents the filter effect.

Then the new tracking signal φ 1 (t) can be expressed as:

φ 1 (t) =φ 1 (t) + λ * h * φ 2 (t) (49)
The meaning of Equation ( 49) is that the new tracking signal of the input signal is taken as the original tracking signal plus the product of the differential signal and the predicted step size.

The input signal is already defined above. The simulation result is shown in Figure 9. It can be seen that the error between the input signal and the tracking signal is significantly reduced compared with Figure 7, indicating that the phase prediction method can well overcome the phase delay problem of the tracking differentiator. The two improvements are applied to design an IADRC to the straight track control. The principle of IADRC is shown in Figure 10. The two improvements are applied to design an IADRC to the straight track control. The principle of IADRC is shown in Figure 10. The two improvements are applied to design an IADRC to the straight track control. The principle of IADRC is shown in Figure 10.

Simulation and Analysis

During the dynamic positioning and tracking process of the ship, the controller calculates the corresponding thrust according to the deviation between the set position and the current position and transmits it to the ship to drive the ship [START_REF] Ye | Design and Implementation of Pseudo-Inverse Thrust Allocation Algorithm for Ship Dynamic Positioning[END_REF]. The ship then feeds the position information back to the controller through the navigation system to form a closed control loop. The control structure of DP system is shown in Figure 11 [START_REF] Wang | ROV Dynamic Positioning[END_REF]. Since the ship system in this paper is a coupled MIMO system, ADRC has good decoupling performance. It means that ADRC can be designed separately for each control channel [START_REF] Zheng | A practical approach to disturbance decoupling control[END_REF]. Therefore, the surge, sway and yaw can be controlled separately in the following study in this paper.

To verify the performance of the controller designed above, a case study ship is used. Its parameters are shown in Table 1 and the mass matrix and damping matrix of the ship follow [START_REF] Zhang | Ship dynamic positioning system based on backstepping control[END_REF].

Simulation and Analysis

During the dynamic positioning and tracking process of ship, the controller calculates the corresponding thrust according to the deviation between the set position and the current position and transmits it to the ship to drive the ship [START_REF] Ye | Design and Implementation of Pseudo-Inverse Thrust Allocation Algorithm for Ship Dynamic Positioning[END_REF]. The ship then feeds the position information back to the controller through the navigation system to form a closed control loop. The control structure of DP system is shown in Figure 11 [START_REF] Wang | ROV Dynamic Positioning[END_REF]. The two improvements are applied to design an IADRC to the straight track control. The principle of IADRC is shown in Figure 10.

Simulation and Analysis

During the dynamic positioning and tracking process of the ship, the controller calculates the corresponding thrust according to the deviation between the set position and the current position and transmits it to the ship to drive the ship [START_REF] Ye | Design and Implementation of Pseudo-Inverse Thrust Allocation Algorithm for Ship Dynamic Positioning[END_REF]. The ship then feeds the position information back to the controller through the navigation system to form a closed control loop. The control structure of DP system is shown in Figure 11 [START_REF] Wang | ROV Dynamic Positioning[END_REF]. Since the ship system in this paper is a coupled MIMO system, ADRC has good decoupling performance. It means that ADRC can be designed separately for each control channel [START_REF] Zheng | A practical approach to disturbance decoupling control[END_REF]. Therefore, the surge, sway and yaw can be controlled separately in the following study in this paper.

To verify the performance of the controller designed above, a case study ship is used. Its parameters are shown in Table 1 and the mass matrix and damping matrix of the ship follow [START_REF] Zhang | Ship dynamic positioning system based on backstepping control[END_REF]. Since the ship system in this paper is a coupled MIMO system, ADRC has good decoupling performance. It means that ADRC can be designed separately for each control channel [START_REF] Zheng | A practical approach to disturbance decoupling control[END_REF]. Therefore, the surge, sway and yaw can be controlled separately in the following study in this paper.

To verify the performance of the controller designed above, a case study ship is used. Its parameters are shown in Table 1 and the mass matrix and damping matrix of the ship follow [START_REF] Zhang | Ship dynamic positioning system based on backstepping control[END_REF]. The mass matrix and damping matrix of the ship are as follows:

M =   1.1274 0 0 0 1.8902 -0.0744 0 -0.0744 0.1278   , D =   0.0358 0 0 0 0.1183 -0.0124 0 -0.0041 0.0308   .
The focus of this paper is to improve the structure of ADRC to achieve the control effect. Therefore, parameter tuning of ADRC is not involved in this paper and the parameters of the controllers in this paper are empirical parameters. The automatic adjustment of parameters will be studied in detail in the future.

Fixed-Point Control and Yaw Control Based on Classical ADRC

To solve the problems of large overshoot, long response time and low positioning accuracy, this paper uses firstly classical ADRC to verify the fixed-point control and yaw control of the ship compared with classical backstepping control.

In the late 1980s, the idea of backstepping was initiated. Some scholars proposed a backstepping design method based on Lyapunov stability for partially linear strict feedback systems [START_REF] Jin | Command-Filtered Backstepping Integral Sliding Mode Control with Prescribed Performance for Ship Roll Stabilization[END_REF]. The idea of backstepping control is to divide the entire design process into multiple steps and each step is designed with virtual control and stabilization functions. Then, the control starts from the lowest-order equation of the system, gradually recurses to higher-order equations, and finally realizes the global adjustment of the system. The control structure composed of multiple subsystems can effectively control the n-order nonlinear system. Therefore, the backstepping method has a certain representativeness. The idea of dynamic positioning backstepping controller can be found in the Reference [START_REF] Song | Study on Ship Dynamic Positioning with Stochastic Backstepping Nonlinear Controller[END_REF].

Simulation under Ideal Sea Condition

Firstly, the fixed-point control and yaw control under ideal sea condition are considered. In the simulation, the initial position signal and yaw signal of the ship are set as η 0 = [0 m 0 m 0 • ] T , and the desired target position and yaw signal are set as η d = [20 m 30 m 45 • ] T . The simulation time is 200 s. The surge, sway, yaw and position under ideal sea condition can be seen in Figures 12 and13, respectively. The speed of the ship is shown in Figure 14. The forces and moment of thrusters are shown in Figure 15. It is obvious from Figures 12 and 13 that both controllers can bring the ship to the set position without overshooting. After reaching the desired position, the ship remains berthed, and the yaw angle reaches the predetermined angle and remains constant. However, ADRC takes about 30 s, while backstepping control takes about 50 s, indicating that the designed ADRC controller can effectively solve the problem of long response time compared with backstepping control. It can be seen from Figure 14 that the linear speed in the surge and sway directions and the angular speed in the yaw direction of the ship with ADRC have large accelerations at the start, and then gradually decelerate. The linear speed in the surge and sway directions of the ship with backstepping control fluctuates greatly in the early stage, and the largest fluctuation is about 2 m/s. The angular speed in the yaw direction of that is consistent with the trend of ADRC, but its maximum speed is only about 2.5 m/s and the deceleration is more stable. After the ship reaches the specified position, the speed decreases to 0 and stays parked. Figure 15 shows the forces and moments of the two controllers. It can be seen that the oscillations of the forces in surge and sway direction in the early stage of backstepping method are more obvious which may cause damage to thrusters and do not exceed about 2 kN and 5 kN. The yaw moment of the backstepping method is much smaller than that of the ADRC, but there is still a small oscillation. The above results prove that the two controllers can both realize fixed-point control and yaw control under ideal sea condition, but according to the comprehensive analysis of position, speed and forces, the performance of ADRC is better than backstepping.

Simulation under Environmental Disturbances

For control system design, it is common to assume the principle of superposition when considering wind and wave disturbances. For most marine control applications, this is a good approximation. In this paper, the fourth-level sea state has 2 m wave height and 10 m/s average wind speed [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF][START_REF] Fossen | Guidance and Control of Ocean Vehicles[END_REF]. It is obvious from Figures 12 and13 that both controllers can bring the ship to the set position without overshooting. After reaching the desired position, the ship remains berthed, and the yaw angle reaches the predetermined angle and remains constant. However, ADRC takes about 30 s, while backstepping control takes about 50 s, indicating that the designed ADRC controller can effectively solve the problem of long response time compared with backstepping control. It can be seen from Figure 14 that the linear speed in the surge and sway directions and the angular speed in the yaw direction of the ship with ADRC have large accelerations at the start, and then gradually decelerate. The linear speed in the surge and sway directions of the ship with backstepping control fluctuates greatly in the early stage, and the largest fluctuation is about 2 m/s. The angular speed in the yaw direction of that is consistent with the trend of ADRC, but its maximum speed is only about 2.5 m/s and the deceleration is more stable. After the ship reaches the specified position, the speed decreases to 0 and stays parked. Figure 15 shows the forces and moments of the two controllers. It can be seen that the oscillations of the forces in surge and sway direction in the early stage of backstepping method are more obvious which may cause damage to thrusters and do not exceed about 2 kN and 5 kN. The yaw moment of the backstepping method is much smaller than that of the ADRC, but there is still a small oscillation. The above results prove that the two controllers can both realize fixed-point control and yaw control under ideal sea condition, but according to the comprehensive analysis of position, speed and forces, the performance of ADRC is better than backstepping.

Simulation under Environmental Disturbances

For control system design, it is common to assume the principle of superposition when considering wind and wave disturbances. For most marine control applications, this is a good approximation. In this paper, the fourth-level sea state has 2 m wave height and 10 m/s average wind speed [START_REF] Fossen | Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles[END_REF][START_REF] Fossen | Guidance and Control of Ocean Vehicles[END_REF].

The initial signal and desired signal are still the same as the simulation under ideal sea conditions. The surge, sway, yaw and position under environmental disturbances can be seen in Figures 16 and17, respectively. The speed of the ship is shown in Figure 18. Figure 19 shows the forces and moment of thrusters under environmental disturbances. It can be seen from Figure 16 to Figure 17 that after adding environmental disturbances, the ship with the two controllers can still reach the desired position in a state of oscillation due to the influence of wind and wave interference. However, it can be observed that ADRC reveals more stable performance with a smoother trajectory than backstepping and the ADRC is still about 20 s faster than backstepping control to reach the set position. Figure 18 indicates the comparison of the speed of three DOF with ADRC and backstepping control. As can be seen from Figure 18, the speed change trend of ADRC is basically the same as that without disturbances, which oscillates much more slightly than backstepping. It is obvious from Figure 19 that there is a huge fluctuation in the forces and moment of thrusters with disturbances. The maximum abrupt values of forces and moment in the three directions are 10 kN, 20 kN and 40 kN•m, respectively. From the positioning results, it proves that the designed controllers can still realize the positioning function to reach the desired position without large overshooting. ADRC has a stronger antiinterference ability than backstepping control, but the forces and moment both have a fluctuation, which may be harmful to the use of thrusters.

In summary, the use of the classical ADRC can realize the fixed-point positioning control and yaw control of the ship, which proves that ADRC does make the ship reach the desired reference quickly without overshooting and have higher positioning accuracy than backstepping control. However, from the results of forces and moment, there is much room for improvement in the future research.

Straight Track Control Based on IADRC

According to the above simulations and analysis of the fixed-point control and yaw control, the classical ADRC has a better effect than backstepping control. However, according to the analysis in Section 5.1, due to the complexity of track control, the problems of dither and phase delay in the classical ADRC may have an impact on the straight track control, and the classical ADRC may be not suitable in this case. Therefore, the IADRC which has been presented in Section 5.2 is applied and compared with the classic ADRC.

To compare the errors before and after improvement more intuitively, the mean value of the square value of the tracking error is used to compare the size of the tracking error, as shown in Equation [START_REF] Zhang | Ship dynamic positioning system based on backstepping control[END_REF]:

ε= ∑ error 2 n ( 50
)
where error is the distance between the desired position and the actual position and n is the number of samples. ε represents the degree to which the actual value deviates from It can be seen from Figure 16 to Figure 17 that after adding environmental disturbances, the ship with the two controllers can still reach the desired position in a state of oscillation due to the influence of wind and wave interference. However, it can be observed that ADRC reveals more stable performance with a smoother trajectory than backstepping and the ADRC is still about 20 s faster than backstepping control to reach the set position. Figure 18 indicates the comparison of the speed of three DOF with ADRC and backstepping control. As can be seen from Figure 18, the speed change trend of ADRC is basically the same as that without disturbances, which oscillates much more slightly than backstepping. It is obvious from Figure 19 that there is a huge fluctuation in the forces and moment of thrusters with disturbances. The maximum abrupt values of forces and moment in the three directions are 10 kN, 20 kN and 40 kN•m, respectively. From the positioning results, it proves that the designed controllers can still realize the positioning function to reach the desired position without large overshooting. ADRC has a stronger anti-interference ability than backstepping control, but the forces and moment both have a fluctuation, which may be harmful to the use of thrusters.

In summary, the use of the classical ADRC can realize the fixed-point positioning control and yaw control of the ship, which proves that ADRC does make the ship reach the desired reference quickly without overshooting and have higher positioning accuracy than backstepping control. However, from the results of forces and moment, there is much room for improvement in the future research.

Straight Track Control Based on IADRC

According to the above simulations and analysis of the fixed-point control and yaw control, the classical ADRC has a better effect than backstepping control. However, according to the analysis in Section 5.1, due to the complexity of track control, the problems of dither and phase delay in the classical ADRC may have an impact on the straight track control, and the classical ADRC may be not suitable in this case. Therefore, the IADRC which has been presented in Section 5.2 is applied and compared with the classic ADRC.

To compare the errors before and after improvement more intuitively, the mean value of the square value of the tracking error is used to compare the size of the tracking error, as shown in Equation ( 50):

ε = ∑ error 2 n ( 50
)
where error is the distance between the desired position and the actual position and n is the number of samples. ε represents the degree to which the actual value deviates from the set value, which is similar to Mean Square Error (MSE) in statistics. MSE is the averaged squared difference between the outputs and the targets [START_REF] Ahmed | Comparative analysis of DGPS predicted corrections using dynamic neural networks[END_REF]. The smaller the MSE is, the less the actual value deviates from the set value and the more stable the system is.

Simulation under Ideal Sea Condition

The initial state of the ship is that the initial position is at the origin of the coordinates. The straight reference track is set as y = x + 5. The set yaw angle is 45 • and the set speed is 1 m/s. In the case of low speed, the ship needs to keep the set yaw angle unchanged when sailing along the set trajectory.

First, the environment is ideal sea condition. The simulation time is 200 s. Figure 20 shows the straight track of the ship. Figure 21 shows the tracking error and Figure 22 shows the speed of the ship. Figures 232425show the forces and moment.

the set value, which is similar to Mean Square Error (MSE) in statistics. MSE is the averaged squared difference between the outputs and the targets [START_REF] Ahmed | Comparative analysis of DGPS predicted corrections using dynamic neural networks[END_REF]. The smaller the MSE is, the less the actual value deviates from the set value and the more stable the system is.

Simulation under Ideal Sea Condition

The initial state of the ship is that the initial position is at the origin of the coordinates. The straight reference track is set as y=x+5. The set yaw angle is 45° and the set speed is 1 m/s. In the case of low speed, the ship needs to keep the set yaw angle unchanged when sailing along the set trajectory.

First, the environment is ideal sea condition. The simulation time is 200 s. Figure 20 shows the straight track of the ship. Figure 21 shows the tracking error and Figure 22 shows the speed of the ship. the set value, which is similar to Mean Square Error (MSE) in statistics. MSE is the averaged squared difference between the outputs and the targets [START_REF] Ahmed | Comparative analysis of DGPS predicted corrections using dynamic neural networks[END_REF]. The smaller the MSE is, the less the actual value deviates from the set value and the more stable the system is.

Simulation under Ideal Sea Condition

The initial state of the ship is that the initial position is at the origin of the coordinates. The straight reference track is set as y=x+5. The set yaw angle is 45° and the set speed is 1 m/s. In the case of low speed, the ship needs to keep the set yaw angle unchanged when sailing along the set trajectory.

First, the environment is ideal sea condition. The simulation time is 200 s. Figure 20 shows the straight track of the ship. Figure 21 shows the tracking error and Figure 22 shows the speed of the ship. It can be seen from Figures 20 and21 that both ADRC and IADRC can effectively make the ship track the given trajectory. As is clearly shown in the two figures, the tracking error of the two controllers gradually decreases from 5 m and the tracking error of them is less than 1 m after the sails steadily. However, the tracking error of IADRC is smaller than that of classical ADRC, which proves that the DP performance of IADRC is better than that of ADRC under ideal sea condition. According to Equation ( 50), the error dispersion before and after improvement can be calculated as ε 1 =0.1644 and ε 2 =0.1260, respectively, in which ε 2 is about 23.4% less than ε 1 . It indicates that the improved controller brings the trajectory closer to the given track. Figure 22 shows the speed of the ship with the two controllers. It can be seen that the ship starts from the origin, It can be seen from Figures 20 and21 that both ADRC and IADRC can effectively make the ship track the given trajectory. As is clearly shown in the two figures, the tracking error of the two controllers gradually decreases from 5 m and the tracking error of them is less than 1 m after the ship sails steadily. However, the tracking error of IADRC is smaller than that of classical ADRC, which proves that the DP performance of IADRC is better than that of ADRC under ideal sea condition. According to Equation ( 50), the error dispersion before and after improvement can be calculated as ε 1 = 0.1644 and ε 2 = 0.1260, respectively, in which ε 2 is about 23.4% less than ε 1 . It indicates that the improved controller brings the trajectory closer to the given track. Figure 22 shows the speed of the ship with the two controllers. It can be seen that the ship starts from the origin, accelerates with a certain acceleration, decelerates after basically synchronizing with the given track, and finally maintains a low speed of about 1 m/s. The maximum speed of the two controllers is around 3.1 m/s. Figures 232425show the comparison of the forces and moment of three DOF with ADRC and IADRC. In the initial stage of the tracking process, in order to enable the ship to achieve fast tracking, the control forces and moment of thrusters change more rapidly, and the maximum abrupt values of forces and moment in the three directions are 8 kN, 150 kN and 100 kN•m, respectively. Although the forces in surge and sway direction after improvement is smaller than that before improvement, both of the two have severe fluctuation. This analysis shows that the IADRC can reduce the error of straight track control compared with the classical ADRC so that the actual trajectory of the ship is more synchronized with the set trajectory, which ensures the stability of sailing.

Simulation under Environmental Disturbances

This case considers the straight track control of the ship under the environmental disturbances the same as in the previous case (fixed-point control) and the initial state is the same as Section 6.2.1. Figure 26 shows the straight track of the ship. Figure 27 shows the tracking error and Figure 28 shows the speed of the ship. Figures 293031show the forces and moment. It can be seen from Figure 26 to Figure 27 that the two controllers can still realize the straight track control of the ship under the environmental disturbances. According to Equation ( 50), the error dispersion before and after improvement can be calculated as ε 1 =0.1095 and ε 2 =0.0736, respectively, in which ε 2 is about 32.8% less than ε 1 . It indicates that the improved controller brings the track closer to the given track and the error is smaller under the environmental disturbances. The speed of the ship can be stabilized It can be seen from Figure 26 to Figure 27 that the two controllers can still realize the straight track control of the ship under the environmental disturbances. According to Equation ( 50), the error dispersion before and after improvement can be calculated as ε 1 = 0.1095 and ε 2 = 0.0736, respectively, in which ε 2 is about 32.8% less than ε 1 . It indicates that the improved controller brings the track closer to the given track and the error is smaller under the environmental disturbances. The speed of the ship can be stabilized at about 1 m/s and the maximum speed is slightly larger than that under ideal sea condition due to the environmental disturbances. It can be seen in Figures 293031that the forces and moment of thrusters in the three directions with ADRC or IADRC all have different degrees of oscillation due to external disturbances. The maximum fluctuation amplitude in the surge and sway directions does not exceed about 10 kN and that in the yaw direction is not more than about 20 kN•m. In conclusion, the IADRC can better realize the straight track control of the ship under environmental disturbances according to the positioning results.

From the above analysis, it can be seen that ADRC can achieve fast and accurate tracking both for the fixed-point control and the straight track control, and the IADRC can significantly reduce the error, indicating that the control methods used in this paper have good dynamic control performance. Although the above controllers can all have good results from a position perspective, the results of forces and moment are not perfect, which are also the focuses of future work.

Conclusions

This paper focuses on DP control of ships based on ADRC. This paper proposes a new ADRC controller which is significantly improved by the uses of fal function filter and phase prediction. The designed ADRC can realize the fixed-point and yaw control in a short time of about 30 s compared with backstepping control. It can reach the set value without overshoot, and it also has better positioning performance under environmental disturbances. With the fal function filter and phase prediction, the phase delay and dither problems of the tracking differentiator are improved. Compared with the classical ADRC, the IADRC can not only realize the straight track control of the ship, but also effectively reduce the error, among which the degree of the error ε decreases by 23.4% under ideal sea condition and by 32.8% under external environmental disturbances. The simulations indicate that the proposed methods can meet the requirements of DP and have better positioning performance than more classical methods. However, although the proposed methods can achieve dynamic positioning, they still have some limitations. The forces and moment oscillations are large, which are harmful to the stable operation of the thrusters and may damage the equipment for a long time. Besides, ADRC has many parameters that need to be tuned. Hence, the future work will focus on improving control output and parameter optimization. Artificial intelligence algorithms, conjugate gradient algorithm, sliding mode control or fuzzy control will be considered to optimize control parameters and adjust control output to improve the control effect.

η

The position and angle vector η 1 ,η 2

The position vector and angle vector respectively x, y, z

The positions (m) ϕ, θ, ψ The angles (deg) ν

The speed vector ν 1 ,ν 2

The linear speed vector and angular speed vector respectively u, v, w

The linear speed (m/s) p, q, r

The angular speed (deg/s) J, J 1 ,J 2

The rotation matrices τ

The total thrust and moment generated by thrusters Fx, Fy, Fn

The force (kN) and moment (kN•m) of thrusters a

The external environmental disturbance force and moment a wind , a wave

The force and moment of wind and wave M

The total inertia matrix M A , M RB The inertia matrix of the hydrodynamic system and the inertia matrix of the rigid body system D

The damping matrix C

The Coriolis-centripetal force matrix

V T , V, V r , V s
The absolute wind speed, the average wind speed, the relative wind speed and the speed of ship (m/s) ϕ wind , ϕ ship

The wind direction and the ship direction (deg) γ

The drift angle (deg) X wind , Y wind , Z wind The force and moment of wind in three directions C X , C Y , C N

The empirical force and moment coefficients ρ α

The density of air (kg/m 3 ) A T , A L

The transverse and lateral projected areas (m 2 ) L

The overall length of the ship (m) g

The gravitational acceleration (m/s 2 ) ω i

The wave frequency of the i-th wave (Hz) H s

The wave height (m) ρ

The density of sea (kg/m 3 ) χ

The encounter angle (deg) µ

The wave amplitude (m) C XW , C YW , C NW The coefficient obtained by regression analysis

Figure 1 .

 1 Figure 1. A cable laying vessel with DP system [3].

Figure 1 .

 1 Figure 1. A cable laying vessel with DP system [3].

  [START_REF] Du | Adaptive Robust Output Feedback Control for a Marine Dynamic Positioning System Based on a High-Gain Observer[END_REF].

Figure 2 .

 2 Figure 2. (a) The NED coordinate and (b) Body-fixed coordinate [28].

Figure 2 .

 2 Figure 2. (a) The NED coordinate and (b) Body-fixed coordinate [28].

Figure 3 .

 3 Figure 3. The structure of ADRC [42].

Figure 4 .

 4 Figure 4. Single control of TD.

Figure 3 .

 3 Figure 3. The structure of ADRC [42].

Figure 3 .

 3 Figure 3. The structure of ADRC [42].

Figure 4 .

 4 Figure 4. Single control of TD.

Figure 4 .

 4 Figure 4. Single control of TD.

Figure 5 .

 5 Figure 5. The output of TD.

Figure 5 .

 5 Figure 5. The output of TD.

When |x (t)| > ( w d k 1 α 1 α

 11 , control variable s = -k|x| α sign(x) can control the steady error caused by w(x, t) within the range of |x| ≤ w d k .When |x| ≤ δ, system (37) becomes:

.

Figure 6 .

 6 Figure 6. The feedback structure of fal function filter [45].

Figure 6 .

 6 Figure 6. The feedback structure of fal function filter [45].

Figure 6 .

 6 Figure 6. The feedback structure of fal function filter [45].

Figure 7 .

 7 Figure 7. The tracking signal after fal function filter.

Figure 7 .

 7 Figure 7. The tracking signal after fal function filter.

31 Figure 8 .

 318 Figure 8. The principle figure of phase prediction.

Figure 8 .

 8 Figure 8. The principle figure of phase prediction.

Figure 9 .

 9 Figure 9. The output of IADRC.

Figure 9 .

 9 Figure 9. The output of IADRC.

Figure 9 .

 9 Figure 9. The output of IADRC.

Figure 10 .

 10 Figure 10. The principle of IADRC.

Figure 11 .

 11 Figure 11. The structure of DP system [48].

Figure 10 .

 10 Figure 10. The principle of IADRC.

Figure 9 .

 9 Figure 9. The output of IADRC.

Figure 10 .

 10 Figure 10. The principle of IADRC.

Figure 11 .

 11 Figure 11. The structure of DP system [48].

Figure 11 .

 11 Figure 11. The structure of DP system [48].

31 Figure 12 .

 3112 Figure 12. The position of surge, sway and yaw under ideal sea condition.

Figure 13 .

 13 Figure 13. The trajectory of the ship under ideal sea condition.

Figure 12 . 31 Figure 12 .

 123112 Figure 12. The position of surge, sway and yaw under ideal sea condition.

Figure 13 .

 13 Figure 13. The trajectory of the ship under ideal sea condition.Figure 13. The trajectory of the ship under ideal sea condition.

Figure 13 .

 13 Figure 13. The trajectory of the ship under ideal sea condition.Figure 13. The trajectory of the ship under ideal sea condition.

Figure 13 .

 13 Figure 13. The trajectory of the ship under ideal sea condition.

Figure 14 .

 14 Figure 14. The speed of the ship under ideal sea condition.Figure 14. The speed of the ship under ideal sea condition.

Figure 14 .

 14 Figure 14. The speed of the ship under ideal sea condition.Figure 14. The speed of the ship under ideal sea condition.

Figure 15 .

 15 Figure 15. The forces and moment under ideal sea condition.

Figure 15 .

 15 Figure 15. The forces and moment under ideal sea condition.

J 31 Figure 16 .

 3116 Figure 16. The position of surge, sway and yaw under environmental disturbances.

Figure 17 .

 17 Figure 17. The trajectory of the ship under environmental disturbances.

Figure 18 .

 18 Figure 18. The speed of the ship under environmental disturbances.

Figure 16 . 31 Figure 16 .

 163116 Figure 16. The position of surge, sway and yaw under environmental disturbances.

Figure 17 .

 17 Figure 17. The trajectory of the ship under environmental disturbances.

Figure 18 .

 18 Figure 18. The speed of the ship under environmental disturbances.

Figure 17 . 31 Figure 16 .

 173116 Figure 17. The trajectory of the ship under environmental disturbances.

Figure 17 .

 17 Figure 17. The trajectory of the ship under environmental disturbances.

Figure 18 .

 18 Figure 18. The speed of the ship under environmental disturbances.Figure 18. The speed of the ship under environmental disturbances.

Figure 18 .

 18 Figure 18. The speed of the ship under environmental disturbances.Figure 18. The speed of the ship under environmental disturbances.

Figure 19 .

 19 Figure 19. The forces and moment under environmental disturbances.

Figure 19 .

 19 Figure 19. The forces and moment under environmental disturbances.

   show the forces and moment.

Figure 20 .

 20 Figure 20. The straight track of the ship under ideal sea condition.

Figure 21 .

 21 Figure 21. Tracking error under ideal sea condition.

Figure 20 .

 20 Figure 20. The straight track of the ship under ideal sea condition.

   show the forces and moment.

Figure 20 .

 20 Figure 20. The straight track of the ship under ideal sea condition.

Figure 21 .

 21 Figure 21. Tracking error under ideal sea condition.Figure 21. Tracking error under ideal sea condition.

Figure 21 .

 21 Figure 21. Tracking error under ideal sea condition.Figure 21. Tracking error under ideal sea condition.

Figure 22 .Figure 22 .

 2222 Figure 22. The speed of the ship under ideal sea condition. Figure 22. The speed of the ship under ideal sea condition.

Figure 22 .

 22 Figure 22. The speed of the ship under ideal sea condition.

Figure 23 .

 23 Figure 23. The Force in the surge direction under ideal sea condition.

Figure 24 .

 24 Figure 24. The Force in the sway direction under ideal sea condition.

Figure 23 .

 23 Figure 23. The Force in the surge direction under ideal sea condition.

Figure 22 .

 22 Figure 22. The speed of the ship under ideal sea condition.

Figure 23 .

 23 Figure 23. The Force in the surge direction under ideal sea condition.

Figure 24 .

 24 Figure 24. The Force in the sway direction under ideal sea condition.Figure 24. The Force in the sway direction under ideal sea condition.

Figure 24 .

 24 Figure 24. The Force in the sway direction under ideal sea condition.Figure 24. The Force in the sway direction under ideal sea condition.

Figure 25 .

 25 Figure 25. The Moment in the yaw direction under ideal sea condition.

Figure 25 .

 25 Figure 25. The Moment in the yaw direction under ideal sea condition.

31 Figure 26 .

 3126 Figure 26. The straight track of the ship under environmental disturbances.

Figure 27 .

 27 Figure 27. Tracking error under environmental disturbances.

Figure 26 . 31 Figure 26 .

 263126 Figure 26. The straight track of the ship under environmental disturbances.

Figure 27 .

 27 Figure 27. Tracking error under environmental disturbances.Figure 27. Tracking error under environmental disturbances.

Figure 27 .

 27 Figure 27. Tracking error under environmental disturbances.Figure 27. Tracking error under environmental disturbances.

Figure 27 .

 27 Figure 27. Tracking error under environmental disturbances.

Figure 28 .

 28 Figure 28. The speed of the ship under environmental disturbances.Figure 28. The speed of the ship under environmental disturbances.

Figure 28 . 31 Figure 29 .

 283129 Figure 28. The speed of the ship under environmental disturbances.Figure 28. The speed of the ship under environmental disturbances. J. Mar. Sci. Eng. 2022, 10, x FOR PEER REVIEW 27 of 31

Figure 30 .

 30 Figure 30. The force in the sway direction under environmental disturbances.

Figure 29 . 31 Figure 29 .

 293129 Figure 29. The force in the surge direction under environmental disturbances.

Figure 30 .

 30 Figure 30. The force in the sway direction under environmental disturbances.Figure 30. The force in the sway direction under environmental disturbances.

Figure 30 .

 30 Figure 30. The force in the sway direction under environmental disturbances.Figure 30. The force in the sway direction under environmental disturbances.

Figure 30 .

 30 Figure 30. The force in the sway direction under environmental disturbances.

Figure 31 .

 31 Figure 31. The moment in the yaw direction under environmental disturbances.

Figure 31 .

 31 Figure 31. The moment in the yaw direction under environmental disturbances.

  [START_REF] Du | Adaptive Robust Output Feedback Control for a Marine Dynamic Positioning System Based on a High-Gain Observer[END_REF].

	J. Mar. Sci. Eng. 2022, 10, x FOR PEER REVIEW	5 of 31

Table 1 .

 1 Parameters of the example ship.

	Entry	Data
	Ship Length	76.2 m
	Ship Width	18.8 m
	Ship Height	82.5 m
	Draft	6.25 m
	Displacement	4200 t
	Power	3533 kW

J. Mar. Sci. Eng. 2022, 10, 865. https://doi.org/10.3390/jmse10070865 https://www.mdpi.com/journal/jmse

J. Mar. Sci. Eng. 2022, 10, x FOR PEER REVIEW

Funding: This research was supported by: (a) National Natural Science Foundation of China, China (Grant No. 61503242, 61673260); (b) Shanghai Science and Technology Commission project (Grant No. 20040501200); (c) Science Technology Commission of Shanghai Municipality and Shanghai Engineering Research Center of Ship Intelligent Maintenance and Energy Efficiency under Grant 20DZ2252300. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable.

Data Availability Statement:

The datas that support the findings of this study are available within the article.

Conflicts of Interest:

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature