Next Article in Journal
Monocular Real Time Full Resolution Depth Estimation Arrangement with a Tunable Lens
Next Article in Special Issue
Automatic Bug Triaging via Deep Reinforcement Learning
Previous Article in Journal
Vibration and Wave Analyses in the Functionally Graded Graphene-Reinforced Composite Plates Based on the First-Order Shear Deformation Plate Theory
Previous Article in Special Issue
Reliability Analysis of Military Vehicles Based on Censored Failures Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order Controller for Course-Keeping of Underactuated Surface Vessels Based on Frequency Domain Specification and Improved Particle Swarm Optimization Algorithm

1
School of Software, Dalian Jiaotong University, Dalian 116026, China
2
School of Computer Science, China West Normal University, Nanchong 637002, China
3
School of Electronic Information and Automation, Civil Aviation University of China, Tianjin 300300, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(6), 3139; https://doi.org/10.3390/app12063139
Submission received: 16 February 2022 / Revised: 14 March 2022 / Accepted: 16 March 2022 / Published: 18 March 2022
(This article belongs to the Special Issue Soft Computing Application to Engineering Design)

Abstract

:
In this paper, a new fractional-order (FO) PIλDµ controller is designed with the desired gain and phase margin for the automatic rudder of underactuated surface vessels (USVs). The integral order λ and the differential order μ are introduced in the controller, and the two additional adjustable factors make the FO PIλDµ controller have better accuracy and robustness. Simulations are carried out for comparison with a ship’s digital PID autopilot. The results show that the FO PIλDµ controller has the advantages of a small overshoot, short adjustment time, and precise control. Due to the uncertainty of the model parameters of USVs and two extra parameters, it is difficult to compute the parameters of an FO PIλDµ controller. Secondly, this paper proposes a novel particle swarm optimization (PSO) algorithm for dynamic adjustment of the FO PIλDµ controller parameters. By dynamically changing the learning factor, the particles carefully search in their own neighborhoods at the early stage of the algorithm to prevent them from missing the global optimum and converging on the local optimum, while at the later stage of evolution, the particles converge on the global optimal solution quickly and accurately to speed up PSO convergence. Finally, comparative experiments of four different controllers under different sailing conditions are carried out, and the results show that the FO PIλDµ controller based on the IPSO algorithm has the advantages of a small overshoot, short adjustment time, precise control, and strong anti-disturbance control.

1. Introduction

In recent years, the maneuverability of underactuated surface vessels (USVs) has become a hot topic of ship control [1,2,3,4,5]. An underactuated system is one in which the number of degrees of freedom of the system is greater than the number of control inputs. An underactuated system refers to a situation where the dimension of the system control input vector space is smaller than the dimension of the system’s generalized coordinate vector space [6]. In the study of ship maneuverability, the effectiveness of ship course control has a direct impact on ship maneuverability and the safety of navigation at sea. For vessels at sea, to keep the ship on the expected course, the rudder must be controlled by a control system to correct for the yaw caused by disturbances. Digital PID rudders are widely used as course controllers for ships at sea, but they are too sensitive to disturbances such as wind and waves, which can lead to frequent steering of the rudder. A great deal of research has been carried out on automatic rudder steering for ships in order to achieve satisfactory maneuvering performance, and many results have been obtained [7,8,9,10,11,12]. In [7], Le proposed to combine fuzzy logic control with linear control theory for the design of a ship autopilot controller and used the controller for ship course control, achieving good control results. In [8], optimal control based on the linear quadratic optimal was proposed and applied to a ship’s autopilot. In [9], the PSO method for an optimal autopilot PID controller was designed. In [10], the H∞ controller was solved by the “2-Riccati” equation, and it was used for the rudder and flap of a ship‘s course under system uncertainty. In [11], a new robust adaptive nonlinear feedback algorithm based on error-driven function was proposed, and it was used for ship course-keeping. In [12], an indirect adaptive robust controller was proposed. First, the parameter estimation law under system uncertainty was studied for a ship model. Secondly, an indirect adaptive robust controller was applied for a ship’s heading.
With the continuous improvement of precision control theory, fractional order calculus theory [13,14,15] and fractional order PIλDµ controllers [16,17,18] have been proposed one after another. In 1999, Podlubny proposed fractional order PIλDµ controllers on the basis of fractional order calculus theory [16]. Compared with the traditional PID controller, the introduction of integral order λ and differential order μ makes the PID controller have two more control parameters, expanding the original “point” control to “surface” control, further improving the control accuracy and making the controller more flexible [19,20,21,22,23,24,25,26,27]. In [19], a two-degrees-of-freedom planar robot was controlled by an FO PID controller tuned with particle swarm optimization (PSO) and a genetic algorithm (GA). In [20], an FO PIλDµ controller was applied to the speed control of subway trains, and optimization studies were carried out to make the train achieve better speed control and stability. In [21], a new FO PIλDµ controller was proposed, and a structural transformation method was designed to obtain the controller parameters, which were used to improve the anti-interference of the controller. In [22], a special case of fractional-order systems was first investigated, and on this basis, the FOPD controller and the FOPID controller were designed, which were used in the control of fractional-order systems to achieve a closed-loop system with a monotonic amplitude-frequency response. In [23], for a three-link rigid manipulator system with coupled multiple-input multiple-output characteristics and high nonlinearity, a fractional-order self-tuned fuzzy PID controller was proposed, which was applied to a robotic manipulator system and achieved good control results. In [24], an adaptive fractional order PID controller was proposed, which was applied to a photovoltaic grid system. In [25], an FOPID controller was proposed in a magnetic ball suspension system, and a modified opposition-based learning technique and the hunger games search algorithm were used to tune the parameters of the FOPID controller. In [26], an FO PIλDµ controller was applied to an automatic voltage regulator (AVR), and a self-regulated off-line optimal tuning method based on the gradient-based optimization algorithm was proposed to tune the five parameters of the FO PIλDµ controller. Davut Izci proposed an FOPID controller for a buck converter system and hybrid Lévy flight distribution and simulated annealing algorithms (LFDSAs), which involved a balanced structure in terms of explorative and exploitative phases being used to tune the parameters of the FOPID controller [27]. In addition, some new algorithms were also proposed to optimize the PID [28,29,30,31,32,33,34].
In this paper, an FO PIλDµ controller is applied to a ship’s course control. A modified Oustaloup approximation method is adopted. Given the values of the gain margin Am and phase margin ϕm, the parameters of the FO PIλDµ controller are calculated. The effectiveness of the algorithm is demonstrated by simulations comparing the FO PIλDµ controller with a PID controller. However, as the FO PIλDµ controller is designed for ship course control at a rated speed, if the ship speed changes dynamically, the non-linear model of the ship will also change accordingly, which will result in the original controller parameters no longer meeting the needs of the ship’s course control and needing to be recalculated. Therefore, this paper proposes an improved particle swarm algorithm (IPSO) for dynamic adjustment of the FO PIλDµ controller parameters. By dynamically changing the learning factor, the particles carefully search in their own neighborhoods at the early stage of the algorithm to prevent them from skipping the global optimum and converging on a local optimum, while at the later stage of evolution, the particles converge on the global optimal solution quickly and accurately to improve the convergence speed. Comparative experiments are carried out without wind and wave disturbances and under wind and wave disturbances. The results show that the FO PIλDµ controller based on the IPSO algorithm has the advantages of a small overshoot, short adjustment time, precise control, and strong anti-disturbance control.

2. USV Model

2.1. Ship Motion Model

The movement of a ship has six degrees of freedom (DOF), which represent the position and direction of the ship’s movement, and therefore six separate coordinates are required to describe these six degrees of freedom. (x, y, z) represents the position of a ship in the inertial frame, and the first-order differentiation of (x, y, z) represents the translational motion along x, y, and z. The first three DOF are surge, sway, and heave. (φ, θ, ψ) represents the direction of a ship, and the first-order differentiation of (φ, θ, ψ) represents the rotation of the ship. The last three DOF are roll, pitch, and yaw [35]. The ship motion description is shown in Figure 1.
Ship maneuvering usually considers only the movement of the ship in the horizontal plane. For most ship movements and their control problems, it is possible to ignore heave, pitch, and roll motions, and thus only the sway, surge, and yaw are discussed. Therefore, the ship motion can be regarded as a plane motion, and the plane motion has only three degrees of freedom, which further simplifies the research object. The maneuvering mathematical model of a USV [36] can be written as
x ˙ = u   c o s ψ v   s i n ψ y ˙ = u   s i n ψ + v   c o s ψ ψ ˙ = r
u ˙ = m 22 m 11 v r d 11 m 11 u + 1 m 11 τ u v ˙ = m 11 m 22 u r d 22 m 22 v r ˙ = m 11 m 22 m 33 u v d 33 m 33 r + 1 m 33 τ r
where x and y denote the surge and sway, ψ denotes the yaw angle, u denotes the surge velocity, v denotes the sway velocity, r denotes the yaw velocity, m11, m22, and m33 are uncertainty parameters representing the intrinsic and additional masses, d11, d22, and d33 are also uncertainty parameters, d11 is the surge hydrodynamic damping factor, d22 is the sway hydrodynamic damping factor, and d33 is the yaw angle hydrodynamic damping factor. For simplicity, the high-order nonlinear damping term is ignored. τu is the longitudinal force of the propeller, and τr is the torque of the propeller.

2.2. Nonlinear Ship Model

The nonlinear ship model is a dynamic system. The rudder angle δ is the input, and the yaw angle ψ is the output. By capturing the main threads of the ship’s dynamics from δ ψ ˙ ψ, the resulting differential equation retains the nonlinear influences and even allows the effects of wind and wave disturbances to be converted into a disturbed rudder angle δD as an input signal that enters the ship’s model together with the actual rudder angle δ. The responding nonlinear ship model is shown in Figure 2.
The Norrbin nonlinear model for course control [37] can be derived from the ship motion maneuvering model in cases where steering is not very frequent:
ψ = a 1 ψ ˙ a 2 ψ ˙ 3 + b δ
where ψ and δ are the yaw angle and rudder angle. a1 = αb, a2 = βb, and b = K/T, where K, T, α, and β are all related to the navigation speed and their values have a large impact on the accuracy of the simulation model.

3. FO PIλDµ Controller

Fractional order (FO) calculus is essentially an arbitrary-order calculus whose order can be real or even complex, which is a further extension of integer order calculus. To improve the accuracy of the PID controller, the differential order µ and the integral order λ are extended to the real domain to obtain the FO PIλDµ controller, whose differential equation [16] is expressed as
u t = K p e t   +   K i D t λ e t   +   K d D t μ e t
where D t μ D t α a c is defined by Caputo and λ > 0 and µ > 0 are any real number and the order of the FO controller. We can then find the Laplace transform of the FO calculus defined by Caputo:
L D t α a c f t = s α F s k = 0 n 1 s a k 1 f k 0
From Equations (4) and (5), the FO PIλDµ controller transfer function can be obtained:
G c s = K p + K i s λ + K d s μ
Since the differential µ and the integral order λ are arbitrary positive real numbers, the PI controller is λ = 1 and µ = 0, the PD controller is λ = 0 and µ = 1, and the PID controller is λ = 1 and µ = 1.
Since λ and µ are allowed to take non-integer values, the adjustment range of the controller parameters is changed from the “point” adjustment of the PID controller in Figure 3a to the “plane” adjustment in Figure 3b, thus enabling the FO PIλDµ controller to control the controlled object more flexibly and to obtain better control results.
In order to better study the FO PIλDµ controllers, researchers wish to extend the design methods of integer-order controllers to FO controllers, which requires approximating the FO calculus to integer-order calculus. This paper uses an improved Oustaloup approximation for the FO calculus operator s α to convert it to integer-order calculus [38]. We set the order to N , where ( ω b , ω h ) is the fitting frequency range and s α is transformed into an integer-order transfer function as follows:
s α = d ω b b α d s 2 + b s ω h d 1 α s 2 + b s ω h + d α 1 + s d ω b / b 1 + s d ω h / b α
where 0 < α < 1 , b > 0 , d > 0 , and d = j ω . We express the FO part of Equation (7) (i.e., K s ) as a rational transfer function in the form of zeros and poles:
K ( s ) = lim N K N ( s ) = lim N k = N N 1 + s / ω k 1 + s / ω k
The k zero pole is
ω k = b d 2 k α 2 N + 1 ω h N + k + 1 2 1 α 2 N + 1 ω b N k + 1 2 1 + α 2 N + 1
ω k = b d 2 k + α 2 N + 1 ω h N + k + 1 2 1 + α 2 N + 1 ω b N k + 1 2 1 α 2 N + 1
An FO calculus operator is obtained by constructing a continuous rational transfer function model:
G s = K d s 2 + b ω h s d 1 α s 2 + b ω h s + d α k = N N 1 + s / ω k 1 + s / ω k
where K = ( ω b ω h ) α .

4. The FO PIλDµ Controller of a Ship’s Course

4.1. FO PIλDµ Controller Design

The working principle of the FO PIλDµ controller system of a ship’s course is shown in Figure 4.
In Figure 4, the actual course is determined with the compass and compared with the given course to generate the course deviation, which enters the FO PIλDµ controller. It then calculates the rudder angle signal to be generated, and the rudder is turned to the required angle when the rudder receives the command signal to make the ship start to change its course. The rudder then returns to its zero position, completing the course revision.
In this paper, given the values of the gain margin A m and phase margin ϕ m , the parameters of the FO PIλDµ controller parameters can be calculated [39]. Starting from the basic definitions of A m and ϕ m , the dynamically controlled object G p s and the controller G c s should satisfy the following relationships:
ϕ m = a r g G c j ω g G p j ω g   +   π
A m = 1 G c j ω p G p j ω p
where ω g meets
G c j ω g G p j ω g   = 1
and where ω p meets
a r g G c j ω p G p j ω p   = π
By substituting G c s with Equation (6), the following relationship is obtained:
K p + K i c o s π λ 2 ω p λ + K d c o s π μ 2 ω p μ = R m p
K p + K i c o s π λ 2 ω g λ + K d c o s π μ 2 ω g μ = R m g
K i s i n π λ 2 ω p λ + K d s i n π μ 2 ω p μ = I m p
K i s i n π λ 2 ω g λ + K d s i n π μ 2 ω g μ = I m g
where
1 A m G p j ω p = R m p + j I m p
c o s φ m j s i n φ m G p j ω g = R m g + j I m g
When the controller is designed, as the controlled object G p s and A m and ϕ m are known from Equations (16)–(19), the error squared minimization equation is introduced into the parameter calculation:
J = 0 e 2 t d t
Here, ω p , ω g , λ , and μ need to meet the following constraints:
ω g λ + μ ω p λ + μ R m p R m g + ω p λ + μ I m p + ω g λ + μ I m g c o t π μ 2 + ω p λ + μ I m g + ω g λ + μ I m p c o t π λ 2 ω p λ ω g μ I m p + ω g λ ω p μ I m g c o t π λ 2 + c o t π μ 2 = 0
If ω p , ω g , λ , and μ are known, then the parameters K p , K i , and K d are expressed as follows:
K p = ω p λ R m p ω g λ R m g ω p λ I m p ω g λ I m g c o t π μ 2 / ω p λ ω g λ
K i = ω p λ ω g λ ω g μ I m p ω p μ I m g ω p λ + μ ω g λ + μ s i n π λ 2
K d = ω p λ I m p ω g λ I m g ω p λ + μ ω g λ + μ s i n π μ 2

4.2. Autotuning of the FO PIλDµ Controller Based on the IPSO Algorithm

However, as the FO PIλDµ controller is designed for ship course control at a rated speed, if the ship speed changes dynamically, the non-linear model of the ship will also change accordingly, which will result in the original controller parameters no longer meeting the needs of the ship’s course control and needing to be recalculated. Therefore, many novel optimization algorithms [40,41,42,43,44,45,46,47,48,49,50,51,52,53,54] are used to tune the parameters of the controller. This paper proposes an IPSO algorithm for dynamic adjustment of the FO PIλDµ controller parameters.
The PSO algorithm was originally designed to optimize the unpredictable movement of simulated flocks of birds. It is an iterative random search algorithm that is easy to implement, is robust, has good parallel processing capabilities, and is likely to obtain a global optimal solution [55,56]. However, in the later stages of optimization, the convergence time is long, the accuracy of the optimization is not high, and it easily falls into local extremes [57,58,59]. The mathematical description of the PSO algorithm [55,56] is as follows. We set the population size to N . At the iteration time, the position vector of each particle in the D -dimensional space is expressed as x ¯ i t = x i 1 , x i 2 , , x i d , , x i D . The velocity of each particle in the D -dimensional space is v ¯ i t = v i 1 , v i 2 , , v i d , , v i D . The optimal position for the particle i to search thus far is p i ¯ t = p i 1 , p i 2 , , p i d , , p i D . The optimal position for all of the particles searched thus far is p g ¯ t = p g 1 , p g 2 , , p g d , , p g D . The velocity v ¯ i t and the position x ¯ i t of particle i are adjusted at time t + 1 as follows:
v ¯ i t + 1   =   ω v ¯ i t   +   c 1 r 1 p ¯ i t     x ¯ i t   + c 2 r 2 p ¯ g t     x ¯ i t
x ¯ i t + 1   =   x ¯ i t   +   v ¯ i t + 1
v i d = v m a x , i f   v i d > v m a x           v i d = v m a x , i f   v i d < v m a x
where i = 1 , N , ω is the inertia weight and non-negative, learning factors c 1 and c 2 are non-negative, and r 1 and r 2 are pseudo-random numbers that obey a uniform distribution on 0 ,   1 . v i d v m a x , v m a x , where v m a x is a constant. If the search space is in x m a x , x m a x , one can set v m a x = k v m a x , where 0.1 k 1.0 .
From Equation (27), it can be seen that in the PSO algorithm, the particle seeking performance is related to the learning factors c 1 and c 2 , which are used to regulate the maximum step size in the direction of motion of the particles p i ¯ t and p g ¯ t respectively. A small learning factor can make the particles oscillate in the region away from the target, while a large learning factor can make the particles move rapidly toward the target region or even move away from the target. The learning factors c 1 and c 2 reflect the exchange of information between particle swarms and determine the influence of a single particle and swarm experience on the particle trajectories. The effect of the values of c 1 and c 2 on the search for an optimum is as follows:
  • If c1 and c2 are both set to 0, the particle flight velocity will remain constant, making it impossible to search the entire space;
  • If c1 = 0, the particle loses its cognitive ability and easily falls into a local optimum solution;
  • If c2 = 0, which leads to a lack of communication and cooperation between particles, the algorithm will fall into a local minimum, and the probability of obtaining an optimal solution will become smaller;
  • Smaller values for both c1 and c2 result in particles flying away from the target area and oscillating in that area;
  • With larger values for both c1 and c2, particles fly faster toward the target area, and this may also cause particles to fly away from the target area.
In general, when using PSO algorithms to solve engineering optimization problems, the goal is to allow the particle in the swarm to perform optimization in the entire search space at the early stage to avoid the problem of lost local optima. In the later stages of the search, the goal is to increase the convergence accuracy and algorithm speed, which can lead to finding the global optimal solution quickly and efficiently. Therefore, the disadvantages of traditional methods can be overcome by introducing the dynamic accelerated learning factors c 1 and c 2 in the PSO algorithms. The maximum initial value of the learning factor is set to 2, while c 1 and c 2 are constructed as functions (i.e., monotonically decreasing and monotonically increasing functions, respectively) with the expressions shown in Equations (30) and (31):
c 1 t = 2 1 e t / t m a x
c 2 t = 2 e t t m a x / t m a x
Therefore, the particle swarm velocity update equation is
v ¯ i t + 1   =   ω v ¯ i t   +   c 1 t r 1 p ¯ i t     x ¯ i t +   c 2 t r 2 p ¯ g t     x ¯ i t
In the early stages of evolution, if c 1 takes a larger value and c 2 takes a smaller value, this will allow the particles to move throughout the search space. In the later stages, if c 1 takes a smaller value and c 2 takes a larger value, this will increase the convergence rate toward the optimal solution.
In this paper, two classical benchmark functions (Griewank and Schaffer) are selected to conduct comparative simulations on the PSO algorithm [31] and the IPSO algorithm with the dynamic acceleration learning factor, which are used to verify the correctness and superiority of the IPSO algorithm.
The Griewank function is as follows:
m i n f x i   = i = 1 N x i 2 4000 i = 1 N c o s x i i   +   1
where x i 600 , 600 . The function has a local minimal value, the number of which is related to the dimensions of the problem, with a global minimum of 0, and x 1 , x 2 , , x n   =   0 , 0 , , 0 can be obtained. The function is typically a nonlinear multi-modal function with a wide search space and is often considered a complex multi-modal problem that is difficult for optimization algorithms to handle.
The Schaffer function is as follows:
m i n f x 1 , x 2 = 0.5 + s i n x 1 + 2 x 2 2 2 0.5 1 + 0.01 x 1 + 2 x 2 2 2
where x 1 , x 2 10 , 10 . The function is a two-dimensional complex function with an infinite number of minimal value points, achieving a minimum value of 0 at (0,0). Since the function has a strongly oscillatory nature, it is difficult to find a globally optimal solution. The curves of the Griewank and Schaffer functions when their independent variables are two-dimensional are shown in Figure 5.
In this paper, the Griewank and Schaffer functions are used to test the performance of the IPSO algorithm. The parameter settings are shown in Table 1, and the test results are shown in Figure 6.
From the simulations, it can be seen that the test results using the Griewank function show that the PSO and IPSO algorithms could converge to the optimal solution, but the number of iterations used to find the optimal IPSO algorithm was significantly smaller than that of the PSO algorithm. The results of the Schaffer function test show that the PSO and IPSO algorithms could converge to the optimal solution, but the number of iterations used to find the optimal solution was significantly smaller for the IPSO algorithm than for the PSO algorithm. Therefore, it is reasonable to believe that the addition of dynamic learning factors to the PSO algorithm in the process of finding the best search faster and with higher convergence accuracy and higher stability to solve the PSO algorithm easily fell into the local optimum problem.
In this part, the IPSO algorithm is used to autotune the five parameters ( K p , K i , K d , λ , and μ ) of the FO PIλDµ controller in Equation (6), which is used for course-keeping of the USV. The working principal diagram and model of the FO PIλDµ controller with the IPSO algorithm are shown in Figure 7 and Figure 8.
The IPSO algorithm is used to autotune the five parameters ( K p , K i , K d , λ , and μ ) of the FO PIλDµ control, which is a five-dimensional solution space. Considering the system overshoot and the square integral of the control loop error, the equation of the optimization objective function is as follows:
J = ϖ 1 0 + e 2 t d t + ϖ 2 σ
where e t is the systematic error, σ is the system overshoot, and ϖ 1 and ϖ 2 are weights.

5. Simulations

5.1. FO P I λ D μ Controller

The research object used the actual ship data of COSCO’s large container vessel of 5446 TEU, whose parameters are shown in Table 2. Let the ship’s rated speed be V = 19.8 knots = 10.186 m/s (1 knot = 1.852 km/s) and the ship’s nonlinear course model adopt the Norrbin nonlinear model. Since the ship could be regarded as a large mass object in movement, the ship’s course control could be regarded as a kind of slow turning movement when the steering was not very frequent, so Equation (3) could be simplified. The simplified equation [37] is shown here:
T ψ · + ψ = K δ
According to the calculation, the model coefficients K = 0.1955 and T = 255.8837 were obtained. Am = 1.3, m = π/4, λ = 0.1, and μ = 0.85 for the FO PIλDµ controller were determined by the system characteristics and practical experience. Then, Kp = 1.4516, Ki = 1.1355, and Kd = 95.0636 were calculated from Equations (22)–(26), respectively. The fitting frequency of the modified Oustaloup filtering method (ωb, ωh) chose [0.001, 1000] and the order N = 5. According to the design method of the PID ship autopilot [37], the Kp, Ki, and Kd are shown here:
K p = T ω n 2 K , K i = T ω n 3 10 K , K d = 2 T ε ω n 1 K
where ω n is the system’s natural frequency and ε is the system’s relative attenuation coefficient. Usually, one would have 0.8 ≤ ε ≤ 1.0, ωn = 0.06, and ε = 0.9. From Equation (37), we could obtain Kp = 4.712, Ki = 0.028, and Kd = 136.243. The FO PIλDµ controller was applied to the ship’s course, where the maximum rudder angle was restricted from −35° to +35°, and the desired course angle was 30° when time was 0–1000 s. The FO PIλDµ controller and PID controller were compared. The comparison results can be seen in Figure 9, Figure 10 and Figure 11 and Table 3.
In the next step, we set the ship’s course angle to change in the following order: (1) 0 ≤ t ≤ 1000 while setting the desired course angle to ψr(t) = 20° and (2) 1000 ≤ t ≤ 2000, while setting the desired course angle to ψr(t) = 30°. The course-tracking curve and rudder angle curve with the FO PIλDµ controller and the PID controller are shown in Figure 12 and Figure 13.
By comparing the simulation experiments, it can be seen that the FO PIλDµ controller had a significantly better control effect than the PID automatic rudder controller proposed in [37]. The FO PIλDµ controller had the advantages of a small overshoot, short adjustment time, and accurate control.

5.2. FO PIλDµ Controller Based on the IPSO Algorithm

Since the ship’s nonlinear model parameters changed with the speed, the calculated FO PIλDµ controller parameters could not meet the dynamic changes of the model’s parameters. The corresponding nonlinear model’s parameter values for ship course control at different speeds are shown in Table 4.
In the next stage, the objective was to check the self-adjustment of the controller based on IPSO. The simulation used the Norrbin nonlinear model in Equation (3). The five parameters of the FO PIλDµ controller were adjusted using the IPSO algorithm for the ship’s course, and the effectiveness of the algorithm was verified by simulation experiments. Following the design step of the FO PIλDµ controller, the fitting frequency of the Oustaloup filtering method ( ω b , ω k ) was [0.001,1000] and N = 5 . The population size was set to 50. The maximum number of the birds’ steps was 50. The particle dimensions were 5, meansing [Kp, Ki, Kd, λ, µ], and the initial values of c1 and c2 were set at 2. The inertia weight was set to 0.9. Four different control strategies were used for the ship’s course, which were the PID controller [37], the PID controller based on the PSO algorithm [9], the FO PIλDµ controller based on the GA algorithm [19], and the FO PIλDµ controller based on the IPSO algorithm. The ship’s speed increased from V = 19.8 knots to V = 24.5 knots at a constant speed, the maximum rudder angle was restricted from −35° to +35°, and the desired course angle was 30° when the time was 0–1000 s. The comparison results can be seen in Figure 14 and Figure 15 and Table 5.
Meanwhile, wind and wave disturbances have been the main causes of yawing of ships when they are sailing at sea. In this design, wind and wave disturbances were introduced into the system, which verified the robustness. The disturbances used white noise to drive a typical second-order oscillation. The wave model [60] obtained under the effect of a force-six wind is that of Equation (38), and the simulation curve is shown in Figure 16.
h ( s ) = 0.4198 s s 2 + 0.3638 s + 0.3675
Four different control strategies were used for a ship’s course under wind and waves, and the results can be seen in Figure 17 and Figure 18.
It can be seen from Figure 17 and Figure 18 that the FO PIλDµ controller based on the IPSO algorithm had a significantly better control effect than the automatic rudder PID controller, the PID controller based on the PSO algorithm, and the FO PIλDµ controller based on a GA algorithm. The FO PIλDµ controller based on the IPSO algorithm had the advantages of a small overshoot, short adjustment time, accurate control, and strong wind and wave resistance.

6. Conclusions

This paper proposed two methods for an FO PIλDµ controller for USV course control. First, this paper proposed an FO PIλDµ controller based on frequency domain specification, which used a desired gain margin and phase margin to meet the performance requirements of the ship’s course control. The simulation results showed that the FO PIλDµ controller had the advantages of a small overshoot, short adjustment time, and precise control over the automatic rudder PID controller. However, if the ship’s speed changed, the ship’s nonlinear course model would change. The parameters of the FO PIλDµ controller no longer met the needs of course control. Secondly, the FO PIλDµ controller based on the IPSO algorithm was proposed for USV course control. To address the influence of the learning factor on the search speed of the particles in the solution space, this paper constructed a particle swarm optimization algorithm with a time-varying learning factor. In the initial stage of optimization, the particles had a stronger self-learning ability and weaker social learning ability to strengthen the global search ability. In the later stage, the particles had a stronger social learning ability and weaker self-learning ability to facilitate convergence at the global optimum. Comparative experiments were carried out without wind and wave disturbances and with wind and wave disturbances. The results show that the FO PIλDµ controller based on the IPSO algorithm had the advantages of a small overshoot, short adjustment time, precise control, and strong anti-disturbance control. The controller output was subjected to wind and wave interference, and the control object itself changed due to certain factors. It could also overcome these effects and had a strong robustness. The control performance was better than that of the three other controllers, and it was easy to modify the existing PID-type automatic rudder while having a high promotion value.
In light of the effectiveness of the FO PIλDµ controller based on the IPSO algorithm, more improvements need to be studied in future. An interesting aspect for future study is how to combine an FO PIλDµ controller with other intelligent algorithms (e.g., neural networks, Lévy flight distribution, the Nelder–Mead algorithm, and fuzzy control) and apply them to a ship’s course control to further verify the superiority of the FO PIλDµ controller in ship course control. Another aspect which needs to be studied involves the use of FO PIλDµ controllers for a ship’s trajectory control.

Author Contributions

Conceptualization, G.L. and Y.L.; methodology, G.L. and Y.L.; software, H.C.; validation, G.L. and H.C.; data curation, G.L. and Y.L.; writing—original draft preparation, G.L. and Y.L.; writing—review and editing, W.D. and H.C.; project administration, H.C.; funding acquisition, G.L. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science Researching Plans of Liaoning Province Education Department under Grant JDL2019003.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Park, B.S.; Kwon, J.W.; Kim, H. Neural network-based output feedback control for reference tracking of underactuated surface vessels. Automatica 2017, 77, 353–359. [Google Scholar] [CrossRef]
  2. Lu, Y.; Zhang, G.; Sun, Z.; Zhang, W. Robust adaptive formation control of underactuated autonomous surface vessels based on MLP and DOB. Nonlinear Dyn. 2018, 94, 503–519. [Google Scholar] [CrossRef]
  3. Liu, C.; Negenborn, R.R.; Chu, X.; Zheng, H. Predictive path following based on adaptive line-of-sight for underactuated autonomous surface vessels. J. Mar. Sci. Technol. 2018, 23, 483–494. [Google Scholar] [CrossRef]
  4. Sahu, B.K.; Subudhi, B.; Gupta, M.M. Stability Analysis of an Underactuated Autonomous Underwater Vehicle Using Extended-Routh’s Stability Method. Int. J. Autom. Comput. 2018, 15, 299–309. [Google Scholar] [CrossRef]
  5. Xie, W.; Reis, J.; Cabecinhas, D.; Silvestre, C. Design and experimental validation of a nonlinear controller for underactuated surface vessels. Nonlinear Dyn. 2020, 102, 2563–2581. [Google Scholar] [CrossRef]
  6. Kolmanovsky, I.; McClamroch, N.H. Developments in nonholonomic control problems. IEEE Control Syst. Mag. 1995, 15, 20–36. [Google Scholar]
  7. Le, M.D.; Nguyen, S.H.; Nguyen, L.A. Study on a new and effective fuzzy PID ship autopilot. Artif. Life Robot. 2004, 8, 197–201. [Google Scholar] [CrossRef]
  8. Gao, C.; Guo, J. On the existence of an optimal control of ship automatic steering instruments. J. Ocean Univ. China 2005, 4, 185–188. [Google Scholar] [CrossRef]
  9. Dlabač, T.; Ćalasan, M.; Krčum, M.; Marvučić, N. PSO-BASED PID controller design for ship course-keeping autopilot. Brodogradnja 2019, 70, 1–15. [Google Scholar] [CrossRef]
  10. Liu, S.; Fang, L. Application of H∞ control in rudder/flap vector robust control of a ship’s course. J. Mar. Sci. Appl. 2007, 6, 62–67. [Google Scholar] [CrossRef]
  11. Zhang, Q.; Zhang, M.; Hu, Y.; Zhu, G. Error-driven-based adaptive nonlinear feedback control of course-keeping for ships. J. Mar. Sci. Technol. 2021, 26, 357–367. [Google Scholar] [CrossRef]
  12. Wu, J.; Zeng, C.; Hu, Y. Indirect Adaptive Robust Control Design for Course Tracking of Ships Subject to Unknown Control Coefficient and Disturbances. Int. J. Control Autom. Syst. 2021, 19, 2059–2067. [Google Scholar] [CrossRef]
  13. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  14. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  15. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  16. Podlubny, I. Fractional-order System and Fractional-order Control. IEEE Trans. Autom. Control 1999, 44, 208–214. [Google Scholar] [CrossRef]
  17. Podlubny, I. Realization of Fractional-Order Control. Acta Montan. Slovaca 2003, 8, 233–235. [Google Scholar]
  18. Podlubny, I.; Petráš, I.; Vinagre, B.M. Analogue Realizations of Fractional-Order Controllers. Nonlinear Dyn. 2002, 29, 281–296. [Google Scholar] [CrossRef]
  19. Bingül, Z.A.F.E.R.; Karahan, O.Ğ.U.Z.H.A.N. Fractional PID controllers tuned by evolutionary algorithms for robot trajectory control. Turk. J. Electr. Eng. Comput. Sci. 2012, 20, 1123–1136. [Google Scholar]
  20. Zhang, C.; Tan, N.; Zhou, T.; Liu, M.; Shan, H. Research on Optimal Control of Subway Train Based on Fractional Order PID Controller. J. China Railw. Soc. 2018, 40, 8–14. [Google Scholar]
  21. Qi, Z.D.; Bian, H.J.; Leng, B.Y.; Shan, L. Design of a new fractional order PIλ–PDµ controller based on SQP. Control Decis. 2016, 31, 2275–2279. [Google Scholar]
  22. Safikhani Mohammadzadeh, H.; Tabatabaei, M. Design of Non-overshooting Fractional-Order PD and PID Controllers for Special Case of Fractional-Order Plants. J. Control Autom. Electr. Syst. 2019, 30, 611–621. [Google Scholar] [CrossRef]
  23. Kumar, J.; Kumar, V.; Rana, K.P.S. Fractional-order self-tuned fuzzy PID controller for three-link robotic manipulator system. Neural Comput. Appl. 2020, 32, 7235–7257. [Google Scholar] [CrossRef]
  24. Nasir, A.; Rasool, I.; Sibtain, D.; Kamran, R. Adaptive Fractional Order PID Controller Based MPPT for PV Connected Grid System under Changing Weather Conditions. J. Electr. Eng. Technol. 2021, 16, 2599–2610. [Google Scholar] [CrossRef]
  25. Izci, D.; Ekinci, S.; Eker, E.; Kayri, M. A novel modified opposition-based hunger games search algorithm to design fractional order proportional-integral-derivative controller for magnetic ball suspension system. Adv. Control Appl. Eng. Ind. Syst. 2022, 4, e96. [Google Scholar] [CrossRef]
  26. Altbawi, S.M.A.; Mokhtar, A.S.B.; Jumani, T.A. Optimal design of Fractional order PID controller based Automatic voltage regulator system using gradient-based optimization algorithm. J. King Saud Univ. Eng. Sci. 2021, in press. [Google Scholar] [CrossRef]
  27. Izci, D.; Ekinci, S.; Hekimoğlu, B. Fractional-Order PID Controller Design for Buck Converter System via Hybrid Lévy Flight Distribution and Simulated Annealing Algorithm. Arab. J. Sci. Eng. 2022, in press. [Google Scholar] [CrossRef]
  28. Cui, H.; Guan, Y.; Chen, H. Rolling element fault diagnosis based on VMD and sensitivity MCKD. IEEE Access 2021, 9, 120297–120308. [Google Scholar] [CrossRef]
  29. Deng, W.; Xu, J.; Gao, X.; Zhao, H. An enhanced MSIQDE algorithm with novel multiple strategies for global optimization problems. IEEE Trans. Syst. Man. Cybern. Syst. 2022, 52, 578–1587. [Google Scholar] [CrossRef]
  30. Ran, X.; Zhou, X.; Lei, M.; Tepsan, W.; Deng, W. A novel k-means clustering algorithm with a noise algorithm for capturing urban hotspots. Appl. Sci. 2021, 11, 11202. [Google Scholar] [CrossRef]
  31. Deng, W.; Zhang, X.X.; Zhou, Y.Q.; Liu, Y.; Zhou, X.B.; Chen, H.L.; Zhao, H.M. An enhanced fast non-dominated solution sorting genetic algorithm for multi-objective problems. Inform. Sci. 2022, 585, 441–453. [Google Scholar] [CrossRef]
  32. Li, T.Y.; Qian, Z.J.; Deng, W.; Zhang, D.Z.; Lu, H.H.; Wang, S.H. Forecasting crude oil prices based on variational mode decomposition and random sparse Bayesian learning. Appl. Soft Comput. 2021, 113, 108032. [Google Scholar] [CrossRef]
  33. Shao, H.D.; Lin, J.; Zhang, L.W.; Galar, D.; Kumar, U. A novel approach of multisensory fusion to collaborative fault diagnosis in maintenance. Inf. Fusion 2021, 74, 65–76. [Google Scholar] [CrossRef]
  34. Deng, W.; Xu, J.; Zhao, H.; Song, Y. A novel gate resource allocation method using improved PSO-based QEA. IEEE Tran. Intell. Transp. Syst. 2022, 23, 1737–1745. [Google Scholar] [CrossRef]
  35. Wu, X.H. Ship Maneuverability and Seakeeping; China Communications Press: Beijing, China, 1999. [Google Scholar]
  36. Wei, M.; Guo, C.; Liu, Y. Robust Adaptive path following for underactuated surface vessels with uncertain dynamics. J. Mar. Sci. 2012, 11, 244–250. [Google Scholar]
  37. Fossen, T.I. Guidance and Control of Ocean Vehicle; John Wiley & Sons Ltd.: Chichester, UK, 1994. [Google Scholar]
  38. Oustaloup, A.; Levron, F.; Mathieu, B.; Nanot, F.M. Frequency-band complex non-integer differentiator: Characterization and synthesis. IEEE Trans. Circuits Syst.-I Fundam. Theory Appl. 2000, 47, 25–39. [Google Scholar] [CrossRef]
  39. Xue, D.Y.; Zhao, C.N. Fractional order PID controller design for fractional order system. Control Theory Appl. 2007, 24, 771–776. [Google Scholar]
  40. Li, T.Y.; Shi, J.Y.; Zhang, D.Z. Color image encryption based on joint permutation and diffusion. J. Electron. Imaging 2021, 30, 013008. [Google Scholar] [CrossRef]
  41. Zhang, Z.H.; Min, F.; Chen, G.S.; Shen, S.P.; Wen, Z.C.; Zhou, X.B. Tri-partition state alphabet-based sequential pattern for multivariate time series. Cogn. Comput. 2021, 1–19. [Google Scholar] [CrossRef]
  42. Chen, H.; Zhang, Q.; Luo, J. An enhanced Bacterial Foraging Optimization and its application for training kernel extreme learning machine. Appl. Soft Comput. 2020, 86, 105884. [Google Scholar] [CrossRef]
  43. Deng, W.; Shang, S.; Cai, X.; Zhao, H.; Zhou, Y.; Chen, H.; Deng, W. Quantum differential evolution with cooperative coevolution framework and hybrid mutation strategy for large scale optimization. Knowl.-Based Syst. 2021, 224, 107080. [Google Scholar] [CrossRef]
  44. Wang, X.; Wang, H.Y.; Du, C.Z.; Fan, X.; Cui, L.; Chen, H.; Deng, F.; Tong, Q.; He, M.; Yang, M.; et al. Custom-molded offloading footwear effectively prevents recurrence and amputation, and lowers mortality rates in high-risk diabetic foot patients: A multicenter, prospective observational study. Diabetes Metab. Syndr. Obes. Targets Ther. 2022, 15, 103–109. [Google Scholar]
  45. Wei, Y.Y.; Zhou, Y.Q.; Luo, Q.F.; Deng, W. Optimal reactive power dispatch using an improved slime mould algorithm. Energy Rep. 2021, 7, 8742–8759. [Google Scholar] [CrossRef]
  46. Wang, L.Z.; Zhang, J.B.; Liu, P.; Choo, K.K.; Huang, F. Spectral-spatial multi-feature-based deep learning for hyperspectral remote sensing image classification. Appl. Soft Comput. 2017, 21, 213–221. [Google Scholar] [CrossRef]
  47. Luo, X.; Shen, Z.; Xue, R. Unsupervised band selection method based on importance-assisted column subset selection. IEEE Access 2018, 7, 517–527. [Google Scholar] [CrossRef]
  48. Chang, C.I.; Kuo, Y.M.; Chen, S. Self-mutual information-based band selection for hyperspectral image classification. IEEE Trans. Geosci. Remote Sens. 2021, 59, 5979–5997. [Google Scholar] [CrossRef]
  49. Lin, Z.; Yan, L. A support vector machine classifier based on a new kernel function model for hyperspectral data. Mapp. Sci. Remote Sens. 2015, 53, 85–101. [Google Scholar] [CrossRef]
  50. He, Z.Y.; Shao, H.D.; Zhong, X.; Zhao, X.Z. Ensemble transfer CNNs driven by multi-channel signals for fault diagnosis of rotating machinery cross working conditions. Knowl.-Based Syst. 2020, 207, 106396. [Google Scholar] [CrossRef]
  51. Li, T.Y.; Shi, J.Y.; Deng, W.; Hu, Z.D. Pyramid particle swarm optimization with novel strategies of competition and cooperation. Appl. Soft Comput. 2021, 121, 108731. [Google Scholar]
  52. Deng, W.; Li, Z.; Li, X.; Chen, H.; Zhao, H. Compound fault diagnosis using optimized MCKD and sparse representation for rolling bearings. IEEE Trans. Instrum. Meas. 2022, in press. [Google Scholar] [CrossRef]
  53. Izci, D. A novel improved atom search optimization algorithm for designing power system stabilizer. Evol. Intell. 2021, in press. [Google Scholar] [CrossRef]
  54. Izci, D. Design and application of an optimally tuned PID controller for DC motor speed regulation via a novel hybrid Lévy flight distribution and Nelder–Mead algorithm. Trans. Inst. Meas. Control 2021, 43, 3195–3211. [Google Scholar] [CrossRef]
  55. Eberhart, R.C.; Kennedy, J. A New Optimizer Using Particle Swarm Theory. In Proceedings of the 6th International on Symposium of Micromachine Human Science, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  56. Kennedy, J.; Eberhart, R.C. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  57. Angeline, P.J. Using selection to improve particle swarm optimization. In Proceedings of the IEEE International Conference on Evolutionary Computation and IEEE World Congress on Computational Intelligence, Anchorage, AK, USA, 4–9 May 1998; pp. 84–89. [Google Scholar]
  58. Naka, S.; Genji, T.; Yura, T.; Fukuyama, Y. A hybrid particle swarm optimization for distribution state estimation. IEEE Trans. Power Syst. 2003, 18, 60–68. [Google Scholar] [CrossRef]
  59. Wonohadidjojo, D.M.; Kothapalli, G.; Hassan, M.Y. Position control of electro-hydraulic actuator system using fuzzy logic controller optimized by particle swarm optimization. Int. J. Autom. Comput. 2013, 10, 181–193. [Google Scholar] [CrossRef] [Green Version]
  60. Jia, X.L.; Zhang, X.K. Intelligent Control of Ship Motion with H∞ Robust Control; Dalian Maritime University Press: Dalian, China, 2002. [Google Scholar]
Figure 1. Ship motion model.
Figure 1. Ship motion model.
Applsci 12 03139 g001
Figure 2. Responding nonlinear ship model.
Figure 2. Responding nonlinear ship model.
Applsci 12 03139 g002
Figure 3. PID controller and FO PIλDµ controller: (a) PID controller and (b) FO PIλDµ controller.
Figure 3. PID controller and FO PIλDµ controller: (a) PID controller and (b) FO PIλDµ controller.
Applsci 12 03139 g003
Figure 4. FO PIλDµ control system of a ship’s course.
Figure 4. FO PIλDµ control system of a ship’s course.
Applsci 12 03139 g004
Figure 5. Distribution chart of test functions.
Figure 5. Distribution chart of test functions.
Applsci 12 03139 g005
Figure 6. Performance comparison of PSO and IPSO algorithms.
Figure 6. Performance comparison of PSO and IPSO algorithms.
Applsci 12 03139 g006
Figure 7. FO PIλDµ control system of ship’s course based on IPSO.
Figure 7. FO PIλDµ control system of ship’s course based on IPSO.
Applsci 12 03139 g007
Figure 8. Process of IPSO optimizing the FO PIλDµ controller.
Figure 8. Process of IPSO optimizing the FO PIλDµ controller.
Applsci 12 03139 g008
Figure 9. Course-keeping curve with FO PIλDµ controller and the PID controller.
Figure 9. Course-keeping curve with FO PIλDµ controller and the PID controller.
Applsci 12 03139 g009
Figure 10. Rudder angle curve with FO PIλDµ controller and the PID controller.
Figure 10. Rudder angle curve with FO PIλDµ controller and the PID controller.
Applsci 12 03139 g010
Figure 11. Steady state error curve with FO PIλDµ controller and the PID controller.
Figure 11. Steady state error curve with FO PIλDµ controller and the PID controller.
Applsci 12 03139 g011
Figure 12. Course-tracking curves with FO PIλDµ controller and PID controller.
Figure 12. Course-tracking curves with FO PIλDµ controller and PID controller.
Applsci 12 03139 g012
Figure 13. Rudder angle curves with FO PIλDµ controller and PID controller.
Figure 13. Rudder angle curves with FO PIλDµ controller and PID controller.
Applsci 12 03139 g013
Figure 14. Course-keeping curves without disturbance.
Figure 14. Course-keeping curves without disturbance.
Applsci 12 03139 g014
Figure 15. Rudder angle curves without disturbance.
Figure 15. Rudder angle curves without disturbance.
Applsci 12 03139 g015
Figure 16. The output curve of wind and wave disturbances.
Figure 16. The output curve of wind and wave disturbances.
Applsci 12 03139 g016
Figure 17. Course-keeping curves under wind and wave disturbances.
Figure 17. Course-keeping curves under wind and wave disturbances.
Applsci 12 03139 g017
Figure 18. Rudder angle curves under wind and wave disturbances.
Figure 18. Rudder angle curves under wind and wave disturbances.
Applsci 12 03139 g018
Table 1. Parameter settings for different PSO algorithms.
Table 1. Parameter settings for different PSO algorithms.
AlgorithmParticles
m
Weight
ω
Learning Factors
c1, c2
Maximum Speed VmaxMaximum Iterations Tmax
PSO2012, 21100
IPSO201Equations (30) and (31)1100
Table 2. Structure data of 5446 TEU container vessel.
Table 2. Structure data of 5446 TEU container vessel.
StructureValueStructureValue
length of ship280 mlength of two columns267 m
width of ship39.8 mdistance of gravity center2.64 m
area of Rudder61.0 m2square coefficient0.67
no-load weight3.5453 M tonsfull load weight6.5531 M tons
draft12.532 mdraft with fully load14.023 m
Table 3. Performance of PID and FO PIλDµ controller.
Table 3. Performance of PID and FO PIλDµ controller.
Controller TypeAdjust Time tsRise Time trOvershoot MEss
FO PIλDµ114754.42%0.002
PID controller5197114.25%0.007
Table 4. Corresponding model parameters under different speeds.
Table 4. Corresponding model parameters under different speeds.
V (knots)KTαβ
27.10.2676186.955610.49157.5343
26.50.2617191.188510.72918.0578
24.50.2419206.795811.604910.1966
22.50.2222225.177612.636613.1644
19.80.1955255.883714.360119.3172
Table 5. Performance of four different controllers.
Table 5. Performance of four different controllers.
Controller TypeAdjust Time tsRise Time trOvershoot MEss
IPSO-FOPID (proposed)1031030.03%0.005
PSO-PID4166731.16%0.024
PID5197114.25%0.007
GA-FOPID163659.53%0.086
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, G.; Li, Y.; Chen, H.; Deng, W. Fractional-Order Controller for Course-Keeping of Underactuated Surface Vessels Based on Frequency Domain Specification and Improved Particle Swarm Optimization Algorithm. Appl. Sci. 2022, 12, 3139. https://doi.org/10.3390/app12063139

AMA Style

Li G, Li Y, Chen H, Deng W. Fractional-Order Controller for Course-Keeping of Underactuated Surface Vessels Based on Frequency Domain Specification and Improved Particle Swarm Optimization Algorithm. Applied Sciences. 2022; 12(6):3139. https://doi.org/10.3390/app12063139

Chicago/Turabian Style

Li, Guangyu, Yanxin Li, Huayue Chen, and Wu Deng. 2022. "Fractional-Order Controller for Course-Keeping of Underactuated Surface Vessels Based on Frequency Domain Specification and Improved Particle Swarm Optimization Algorithm" Applied Sciences 12, no. 6: 3139. https://doi.org/10.3390/app12063139

APA Style

Li, G., Li, Y., Chen, H., & Deng, W. (2022). Fractional-Order Controller for Course-Keeping of Underactuated Surface Vessels Based on Frequency Domain Specification and Improved Particle Swarm Optimization Algorithm. Applied Sciences, 12(6), 3139. https://doi.org/10.3390/app12063139

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

Article Metrics

Back to TopTop