Next Article in Journal
Estimation and Reliability Research of Post-Earthquake Traffic Travel Time Distribution Based on Floating Car Data
Previous Article in Journal
Upgraded Combined Inject-and-Transfer System for Serial Femtosecond Crystallography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pursuit–Evasion Problem of Unmanned Surface Vehicles in a Complex Marine Environment

1
College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China
2
School of Marine Engineering, Jimei University, Xiamen 361021, China
3
China Ship Scientific Research Center, Shanghai 200011, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(18), 9120; https://doi.org/10.3390/app12189120
Submission received: 17 August 2022 / Revised: 6 September 2022 / Accepted: 8 September 2022 / Published: 11 September 2022
(This article belongs to the Section Marine Science and Engineering)

Abstract

:
In order to solve the interception problem of escaping targets under environmental interference and inaccurate target perception, this paper introduces the guidance interception algorithm and introduces a Kalman filter to predict the trajectory of escaping targets. In order to make a USV pursuit under a disturbance in a time-varying real marine environment, a disturbance observer is proposed to estimate the external disturbance and model uncertainty. The simulation of pursuit–evasion results shows that the interception algorithm designed enables the pursuer to intercept the escaping target in a time-varying, disturbed marine environment. Both the simulation and experimental results are provided to validate the effectiveness of the method.

1. Introduction

In recent years, with the continuous development of intelligent control theory and with continuous improvements to actual demand, unmanned surface vehicles (USVs) have been considered for a wide range of potential applications in ocean exploration, coastal defense, and intelligent cruises [1,2]. Under an efficient communication network technology [3], there has been a lot of work and focus on the path tracking control of USVs, such as environmental monitoring, hydrological measurement, and path tracking [4,5]. Since the pursuit–evasion problem involves many hot issues in control fields, this issue has always been the focus of intelligent systems [6,7]. Nevertheless, the pursuit–evasion of USVs is still an open problem that needs to be overcome.
In the implementation of target interception, Ref. [8] intercepted an escaping target based on the initial position and relative orientation of the tracker and the escaping target. The authors of [9] achieved the capture of an escaping target, but the motion model assumed in that article only allows the evader and the hunter to move in the horizontal or vertical directions, which is obviously inconsistent with actual motion. In addition, a method based on Reinforcement Learning [10] was gradually applied to solve the pursuit–evasion problem and allowed multiple robots to track a single target. Ref. [11] encircled escaping targets with multiple hunters in an environment of obstacles. By establishing multiple encirclement points around the target, multiple pursuers blocked the escape route of the evader in all directions [12]. Additionally, the encirclement strategy is often used to encircle static targets [13]. However, once the trajectory of the dynamic target becomes complex, the applicability of this algorithm becomes questionable. Although more pursuers can ensure a higher success rate of hunting, the strategy of multiple pursuers is based on tracking and interception by a single pursuer.
An accurate target location is of great significance for intercepting the escaping target. In article [14], a target position estimator was proposed, and a corresponding control law was designed. However, it cannot be ignored that the escaping target has high mobility. Moreover, the interception direction of the pursuer should be time-varying; that is, it should be changed according to the movement of the target. The Kalman filter is widely used because it can predict the motion and position of the target simply and conveniently. To intercept the escaping target more accurately, the Kalman filter algorithm is adopted in this paper to track and predict an escaping target’s trajectory. For target tracking, the algorithm based on line of sight (LOS) is generally used in the marine field [15]. However, the LOS cannot predict the target trajectory in advance for interception. As a simple and effective interception algorithm, the missile guidance law is widely used in the aviation field. Therefore, the missile guidance law is adopted to intercept escaping vessels in the marine field.
When performing tasks in an actual environment, an USV will be affected by environmental interference by the wind, waves, and currents in a marine environment [16] and by threats from moving obstacles such as sailing ships [17]. A small USV will be more sensitive to external disturbances due to its small main scale and displacement [18]. Although Ref. [19] carried out sea trials of an USV and avoided obstacles and rescued people from drowning, it did not consider interference from the external environment. How to deal with uncertain nonlinearity and interference is an important problem in target tracking and interception [20]. However, the above article does not consider the uncertainty of the environment, which brings challenges to the interception of dynamic escaping targets.
Previous studies on the pursuit–evasion problem had several shortcomings. First, many studies set the escaping target as static, which is not consistent with pursuit–evasion movement in an actual scenario. Second, in order to validate the effectiveness of the interception strategy, only the path of the interceptor is optimized, and the motion of the escaping target is simplified to linear motion, which is also inconsistent with an actual scene. Finally, at present, there are few experiments on the pursuit–evasion problem and the interception strategy is only validated by simulation.
In view of the above challenges, we aim to design a strategy that can intercept escaping targets. Compared with previous work, the contribution of this paper is mainly reflected in the following three aspects. Firstly, considering the dynamic curve motion of the target, a surface interception strategy is constructed to intercept fast maneuvering targets. Secondly, this paper considers predictability based on the Kalman filter, which can improve the efficiency of interception. Thirdly, a disturbance observer is proposed, for which the accuracy of estimating a disturbance is independent from the designed pursuit algorithm.
The rest of this article is organized as follows: in the Section 2, the basic assumptions of the pursuit–evasion problem are introduced, and the guidance algorithm is adopted to solve the interception problem of evaders. The kinematics and dynamics of USV and the Kalman filter algorithm are also given. In the Section 3, the controller and disturbance observer are designed, and a stability analysis of the system is presented. A simulation of the pursuit algorithm is presented in the Section 4. The Section 5 provides the experimental results. Finally, some conclusions are drawn in the Section 6.

2. Problem Formulation and Tracking Strategy

2.1. Pursuit–Evasion Problem Statement

Suppose that there are pursuit USVs and escaping target USVs in a sea area. When the escaping target ship enters the detection range of the pursuit USV, the interceptor will exhibit an interception behavior. The target also has the ability to perceive the pursuer; when the target USV senses the surrounding captors, it will change its route. The pursuer cannot determine the escape strategy of the target and can only adjust its interception direction in real time according to the perceived position and speed of the target.

2.2. Guidance Method in Pursuit–Evasion Games

Suppose at time k that the pursuer P ( x k , y k ) has perceived the escaping target E ( x E k , y E k ) , and the line l formed by the pursuer pointing to the escaping target E is called the line of sight. The angle formed by the line of sight and the x-axis is the line of sight angle. The angle between the target speed V E and l is ψ E , and the angle between the pursuer speed and l  is ψ P . The angle between the line connecting P  and the virtual hunting point E , and the x-axis is the desired heading angle ψ d . A schematic diagram of the guidance method is shown in Figure 1.
The relative motion equation of the pursuer and the evader is as follows:
{ l ˙ = V P cos ψ p V E cos ψ E l θ ˙ = V P sin ψ p V E sin ψ E
The normal speed of the pursuer needs to be equal to the escaping target, that is θ ˙ = 0 ; substituting it into the Formula (1), we have the following:
V P sin ψ p = V E sin ψ E
According to the geometric relationship, the desired heading angle is as follows:
ψ d = sin 1 ( V E V P sin ψ E ) + θ

2.3. Kinematics and Dynamics of USV

A reasonable mathematical model of motion is the basis of ship motion simulation and control algorithm design. The motion model of USV is simplified to a three degree of freedom (DOF) motion model of the water surface. The motion needed to be taken into consideration is in the horizontal plane, including the surge, the sway, and the yaw. The vessel motion model in the horizontal plane is shown in Figure 2.
The model of a USV with three DOFs can be described as follows:
{ x ˙ = u cos ψ v sin ψ y ˙ = u sin ψ + v cos ψ ψ ˙ = r
In the above expressions, x , y , and ψ denote the position and orientation in the earth-fixed coordination system, and u , v , and r are the velocity vectors of the USV with respect to the body-fixed frame, which represent the surge, the sway, and the yaw.
The dynamic equation of the USV when encountering an external disturbance is as follows:
{ u ˙ = m 22 m 11 v r d 11 m 11 u + 1 m 11 ( τ u + τ w u ) v ˙ = m 11 m 22 u r d 22 m 22 v + 1 m 22 ( τ v + τ w v ) r ˙ = m 11 m 22 m 33 u v d 33 m 33 r + 1 m 33 ( τ r + τ w r )
In the above expressions, m i i is the positive definite inertia matrix and d i i is the damping matrix. τ u i and τ r i represent the longitudinal force and yaw moment of the USV, and τ w u , τ w v , and τ w r represent the influence of an external time-varying uncertain disturbance on the USV motion.

2.4. Kalman Filter Algorithm

In this paper, the Kalman filter algorithm is used to estimate and predict the position of the target in real time, and the target state is constantly updated during a pursuit. The dynamic equation and measurement equation of the discrete-time system of the Kalman filter can be expressed as follows:
{ x ^ k = F x ^ k 1 + μ η Z k = H X k + ω
where x ^ k refers to the state estimation at time k 1 and F represents the state transition matrix. η is a white noise series subject to Gaussian distribution with a mean value of zero and a positive definite covariance matrix of Q . μ is the process noise gain coefficient. Z k is the observed measurement of the system, H is the observation coefficient matrix, and ω is the white noise series subject to Gaussian distribution with zero mean and positive definite covariance matrix. R , η , and ω are independent from each other.
The equation for the covariance matrix of the state vector at the predicted state time is as follows:
P k = F P k 1 F T + Q
where P k represents the error covariance matrix at time k . The updated prediction state and noise covariance matrix with observations is as follows:
K k = P k H T ( H P k H T + R ) 1
where K k represents the gain of the filter at time k , which is the intermediate calculation result of the filter. R is the measurement noise, and its value can be directly set according to the accuracy of the observation equipment.
The updated prediction state and noise covariance matrix with observations is:
{ x ^ k = x ^ k + K k ( Z k H x ^ k ) P k = ( I 4 × 4 K k H ) P k
where ( Z k H x ^ k ) represents the residual error of the actual observation and the predicted observation, which will correct the a priori value together with the Kalman gain. P k represents the covariance of a priori estimation of x ^ k at time k, and P k represents the covariance of the posterior estimate of x ^ k . P represents the error covariance, and a smaller value represents a more accurate prediction of the state. The initial value only affects the initial convergence speed, and the value of P will be updated with the iteration of the Kalman filter. Set initial value P 0 = 0 .
The motion equation of the escaping target object can be expressed as follows:
{ x k = x k 1 + v k 1 x T v k 1 x = v k 1 x + α k 1 x T y k = y k 1 + v k 1 y T v k 1 y = v k 1 y + α k 1 y T
where x k 1 x , v k 1 x , and α k 1 x are the position, velocity, and acceleration of the target in the x-axis direction at time t k 1 , and y k 1 , v k 1 y , and α k 1 y are the position, velocity, and acceleration of the target in the y-axis direction at time t k 1 . T is the time interval between two adjacent sampling points. When the sampling frequency is approximately 10 Hz, the value of T can be taken as 0.1 s. The above formula in matrix form can be written as follows:
[ x k v k x y k v k y ] = [ 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ] [ x k 1 v k 1 x y k 1 v k 1 y ] + [ 0.5 T 2 0 T 0 0 0.5 T 2 0 T ] [ α k 1 x α k 1 x α k 1 y α k 1 y ]
The dynamic noise vector of the system is as follows:
η = [ 0.5 α k 1 x T 2 , α k 1 x T , 0.5 α k 1 y T 2 , α k 1 y T ] T
After the above system state equation and observation equation are established, the Kalman filter equation can be used to continuously predict the position of the target at the next time.
H = [ 1 0 0 0 0 0 1 0 ]

3. Controller Design of USV

3.1. Disturbance Observer Design

The external environmental disturbance is treated as the combined effect of the uncertainty of the model and disturbances by the wind, waves, and the current in the external marine environment. The environmental disturbance force can be expressed as follows:
τ E = [ τ w u   τ w v   τ w r ] T
The disturbance is regarded as an uncertain time-varying disturbance with upper bound and slow variation. τ w u , τ w v , and τ w r represent the influence of ocean disturbances and model uncertainty acting on the surge, sway, and yaw directions, respectively.
In this paper, the disturbance in the surge and yaw directions was estimated online, and the trajectory tracking controller with disturbance compensation was designed in combination with the output of the disturbance observer:
{ τ ^ e = μ + K e M υ μ ˙ = K e ( C ( υ ) υ + D ( υ ) υ τ K e M ν ) K e μ
where τ ^ e = [ τ ^ e u , τ ^ e v , τ ^ e r ] Τ is introduced to represent the external estimation and K e = [ K e u , K e v , K e r ] Τ is the disturbance gain.
The estimation error of each degree of freedom disturbance variable is defined as follows:
{ τ ˜ w u = τ w u τ ^ w u τ ˜ v = τ v τ ^ v τ ˜ w r = τ w r τ ^ w r
The Lyapunov function is designed as follows:
V w u = 1 2 τ ˜ w u 2
Taking the derivation of both sides of Equation (17) and substituting Equation (15) into it, we obtain the following:
V ˙ w u = τ ˜ w u τ ˜ ˙ w u = τ ˜ w u ( τ ˙ w u τ ^ ˙ w u ) = τ ˜ w u τ ˙ w u + k w u τ ˜ w u { C ( υ ) υ + D ( υ ) υ τ k w u M υ + μ + M υ ˙ } = k w u τ ˜ w u 2 0
Therefore, the longitudinal disturbance observation error τ ^ w u is globally exponentially stable, and similarly, the disturbance observations τ ^ w v and τ ^ w r can also accurately estimate the actual disturbance. The disturbance observer and control algorithm designed are independent from each other; the accuracy of disturbance estimation does not affect the controller design.

3.2. Heading Control

In the control strategy, define the angle tracking error as follows:
ψ ˜ = ψ ˜ ψ d
In order to adopt the first-order sliding mode to converge the heading tracking angle error, it is necessary to obtain the yaw rate error of the USV. The yaw rate error r ˜ can be obtained by the time derivative of Equation (19):
r ˜ = r r d
Considering the fast convergence of sliding mode control, the first-order sliding surface is described by the following formula:
S φ r = ψ ˜ ˙ + λ φ ψ ˜ = r r d + λ φ φ ˜
where λ φ r is a positive constant and it represents the sliding mode design variable.
The Lyapunov function considering heading angle and external disturbance can be expressed as follows:
V φ r = 1 2 S φ 2 + 1 2 τ ˜ w r 2
By differentiating the above formula, we obtain the following:
V ˙ φ r = S φ S ˙ φ + τ ˜ r 2 τ ˜ ˙ r = S φ ( r ˙ r ˙ d + λ φ ( r r d ) ) + τ ˜ r 2 τ ˜ ˙ r = S φ ( m 11 m 22 m 33 u v d 33 m 33 r + 1 m 33 ( τ r + τ w r ) r ˙ d ) + S φ ( λ φ ( r r d ) ) + τ ˜ r 2 τ ˜ ˙ r
To ensure that the pursuer can hunt the target more effectively, the heading angle controller is designed as follows:
τ r = m 33 ( m 11 m 22 m 33 u v d 33 m 33 r + τ ^ e r ) + r ˙ d λ φ ( r r d ) - k τ r S φ r
where k τ r is a positive constant and it represents the heading angle gain parameter.
Substituting the above formula into Formula (24), we obtain the following:
V ˙ φ r = k τ r S φ 2 k w r τ ˜ r 2 τ ˜ ˙ r 0
It can be known that the sliding surface S φ r will converge to zero in finite time and that the heading angle error will also converge to zero, which means lim t ( ψ ˜ , r ˜ ) = ( 0 , 0 ) .

3.3. Surge Control

Longitudinal velocity tracking error variables can be defined as follows:
lim t u u d
In the pursuit problem, u d is the pursuer’s maximum speed. The surge following error u e is described as follows:
u e = u u d
The Lyapunov function considering longitudinal thrust and external disturbance is designed as follows:
V u e = 1 2 u e 2 + 1 2 τ ˜ w u 2
Taking the derivative of the above formula, we obtain the following:
V ˙ u e = u e ( m 22 m 11 v r d 11 m 11 u + 1 m 11 ( τ u + τ e u ) u ˙ d ) + τ ˜ w u τ ˜ ˙ w u
In order to ensure the convergence of the velocity term, the design surge control law for τ u is chosen as follows:
τ u = m 11 ( k u u e m 22 m 11 v r + d 11 m 11 u + u ˙ d ) τ ^ e u
where k u is a positive constant. If the surge force designed in Formula (31) can ensure the function value is always negative, then the actual velocity of the USV will be consistent with the desired value. Substituting the above expression into Formula (29), we have the following:
V ˙ u e = k u u e 2 k w u τ ˜ w u 2 < 0
According to Formula (31), the velocity error will gradually converge to zero, and the velocity of the pursuer will reach the desired rate. When the direction and velocity of the pursuit USV reach the desired value and converge continuously, the escaping target can be intercepted by the pursuer.

4. Simulation and Analysis

In order to validate the effectiveness of intercepting escaping targets under the guidance strategy proposed in this paper, a simulation experiment was carried out for the algorithm. The ship model used in the simulation is the scale model of CyberShip II (CS2) [21], which is the test tool of the Marine Cybernetics Laboratory (MCLab).
The hydrodynamic parameters are also taken from the parameters of the well-known surface vehicle CS2. In order to test the robustness of the proposed control algorithm to external disturbances, it is assumed that there is an error of 10% between the nominal value and the actual value of the motion parameters. indicates the nominal value and the specific nominal value of hydrodynamic parameters, which can be found in the literature [21]. The actual values of the motion parameters are shown in Table 1.
The parameter model of USV was taken as a constant disturbance and summed up in the disturbance model for external wind, waves, and currents. The environmental disturbance model is defined as follows:
{ τ w u = ( h ( s ) w u ( s ) + ξ w u ) m 11 τ w v = ( h ( s ) w v ( s ) + ξ w v ) m 22 τ w r = ( h ( s ) w r ( s ) + ξ w r ) m 33
with
h ( s ) = K ω s s 2 + 2 ζ ω 0 s + ω 0 2
where w u , w v , and w r are Gaussian white noise; ξ w u , ξ w v , and ξ w r are the uncertainty error coefficient of the model. ζ = 0.3 , K ω = 0.25 , σ m = 0.4 , ω 0 = 0.8 , and k w u = k w v = k w r = 0.1 . The maximum values of τ e u / m 11 , τ e v / m 22 , and τ e r / m 33 are 0.1 m/s2, 0.1 m/s2, and 0.5 m/s2, respectively.
The recommended values of parameters in the controller in the simulation are shown in Table 2.
The initial state of the pursuer is cruising. When the escaping target enters its sensing range, the pursuer adjusts its state and begins its pursuit. The escaping target moves along a curved path. The curved path can be parameterized as follows:
{ ( x E , y E ) = ( 0 . 4 ω + 10   ,   0 . 6 ω + 10 )                                                     , 0 ω < 100 ( x E , y E ) = ( 0 . 6 ω 10 ,   20 sin ( 0 . 2 ( ω / 25 4 ) ) + 70 )   ,   ω 100
where the path parameter is ω = V E / x E 2 + y E 2 . Within the range of simulation time t = 0–170 s, the trajectories of the pursuer and evader are shown in Figure 3.
Set the parameters in the filter as the parameters in the actual observer to maintain the consistent performance for the actual observer. H is the measurement vector of the sensor, which is the coordinate value obtained through the UWB in this system. As the source of the systematic observation data, the UWB positioning equipment has a coordinate accuracy of 0.1 m, so the element value in matrix R was taken as 0.1. The element value in matrix Q was taken as 0.1.
The observation trajectory obtained by the observation station is obviously oscillating, but the observation noise value of the observer is relatively small, so it is closer to the value of the real trajectory. This measurement value can be directly used for target-locking and interception in the simulation. However, since the actual observer may have abnormal values, it is necessary to determine and discard abnormal values. The position errors before and after filtering are shown in Figure 4.
The detection range of the pursuer is a circular area with a radius of 20 m. When the escaping target enters the pursuer’s perceptual range around t = 23 s, the pursuer adjusts the direction of movement and performs the pursuit. However, due to the limitation of the minimum turning radius of the pursuit USV and the influence of the evader’s rapid adjustment of the direction of motion, there will be a short deviation in the route when the pursuer tracks the evader in real time.
Figure 5 shows the angle changes in the pursuer and evader. When the evader enters the perceptual range of the pursuer around t = 23 s, the pursuer adjusts the motion angle to intercept the target.
From Figure 6, it can be seen that in the process of performing the pursuit, although the evader adopts a variety of motion modes, such as linear motion and curve motion, the pursuer always maintains a good tracking effect. It can be seen that in about 150 s, the distance between the pursuer and the escaping target is less than 1 m. This is enough to ensure successful interception of the escaping target by the pursuit USV.
In order to test the estimation effect of the disturbance observer on the uncertainty of the system, the actual value of the total disturbance momentum of the system in this simulation is compared with the estimated value. Because the USV belongs to the underactuated control system, only the designed thrust disturbance and rudder force disturbance are compared. The observation results of system disturbance are shown in Figure 7 and Figure 8.
From the disturbance observation results, it can be observed that the estimated values of system disturbances τ w u and τ w r fit the expected values more accurately. The designed disturbance observation algorithm has a good estimation effect on unknown system parameter perturbation and the time-varying external disturbance force.

5. Experimental Results

In order to further validate the effectiveness of the interception algorithm proposed in this paper, the pursuit-evasion experiment ispresented in this section. Two USVs with the same principal dimensions are defined as the evader (yellow) and the pursuer (orange). The mass of the USV is 12 kg, the length is 0.8 m, and the width is 0.24 m. The maximum speed of the evader and the pursuer are 0.5 m/s and 0.6 m/s, respectively. The position is measured by the ultra-wideband (UWB), and the positioning accuracy is 0.1 m. In addition, the pursuer is equipped with cameras to hunt the evader. The composition of the USV is shown in Figure 9.
The UWB position measurement system is composed of three base station modules at fixed positions, and the final positioning data are generated only on the “base station” module connected to the controller. The “label” module is fixed on the moving ship. The “tag” can be easily added or reduced. Figure 10 shows the structure of UWB position.
UWB is used to obtain and predict the position of the pursuit vehicle and the target vehicle through the Kalman filtering algorithm. Since the movement speed of the pursuer and target in this test changes slowly, and the sampling frequency of UWB is greater than 10 Hz, the ship can be regarded as moving at a uniform speed within the sampling frequency. Due to the noise and error of the data directly obtained by the UWB sensor in real time, the data need to be processed by the filtering algorithm. First, obviously wrong UWB data are eliminated. Data points beyond the scope of the tested water area or beyond the radius of 2 m of the current position obtained by the least squares method will be regarded as abnormal measurement data points. Then, the least squares method is used to perform second-order fitting on 20 data points that were recently saved. The fitting result is taken as the prediction value, and the Kalman filter algorithm is used to correct the latest measurement value to obtain the optimal target motion information.
In the test, the evader always moves at a heading angle of 10°. The pursuer senses the target and executes an interception in order to show the trajectories of the evader and the pursuer more clearly. Figure 11 shows their trajectory curves according to the UWB data. Although the escaping target moves along a curve, it can still be seen that the proposed encirclement algorithm can intercept the escaping target.

6. Conclusions

Combining the guidance algorithm and the underactuated characteristics of a USV, this paper proposes an interception algorithm for escaping targets. For the escaping target with a flexible escape route, the algorithm can still achieve rapid interception. According to the Lyapunov stability theory, it is validated that the designed controller can ensure the convergence of the motion error of USV in a complex marine environment. Finally, the simulation and experimental results validate the effectiveness of the interception algorithm. In this paper, the uncertainty of the environment and the flexibility of the escaping target are comprehensively considered. In addition, the missile guidance method is applied to the interception field, and the trajectory of the escaping target is further predicted by combining the Kalman algorithm.
However, this study still has several shortcomings. First, although the controller considering uncertain external disturbance is designed in the algorithm, due to the lack of equipment, this function is not realized in the experiment. Second, the route of the escaping target should be adjusted in real time according to the state of the pursuer. In this paper, the escaping target moves according to the predefined route in advance. In addition, obstacles are an important factor affecting interception. In future research, we plan to consider the above shortcomings. Future work will further consider the game strategy between the evader and the pursuer, and the obstacle avoidance algorithm will be applied to the pursuit algorithm in more complex environments.

Author Contributions

Conceptualization, Z.S. and H.S.; Formal analysis, Z.S. and X.L.; Funding acquisition, H.S.; Methodology, Z.S. and H.S.; Software, L.D. and P.L.; Writing—original draft, Z.S.; Writing—review and editing, X.L. and P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 52071100), by the National Natural Science Foundation of China (Grant No. 51909103), and by The Fundamental Research Funds for the central Universities (No. 3072022QBZ0102).

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. Tang, P.; Zhang, R.; Liu, D.; Huang, L.; Liu, G.; Deng, T. Local reactive obstacle avoidance approach for high-speed unmanned surface vehicle. Ocean Eng. 2015, 106, 128–140. [Google Scholar] [CrossRef]
  2. Sun, Z.; Sun, H.; Li, P.; Zou, J.; Zhuang, J.; Tan, G. An innovative distributed self-organizing control of unmanned surface vehicle swarm with collision avoidance. Ocean Eng. 2022, 254, 111342. [Google Scholar] [CrossRef]
  3. Zhang, S.; Xiang, X.; Yang, S. Communication Network and Qos Evaluation for Formation Control of Unmanned Surface Vehicles. Brodogradnja 2020, 71, 19–37. [Google Scholar] [CrossRef]
  4. Van, M. An enhanced tracking control of marine surface vessels based on adaptive integral sliding mode control and disturbance observer. ISA Trans. 2019, 90, 30–40. [Google Scholar] [CrossRef]
  5. Lc, A.; Cl, A.; Bing, X.B.; Yg, A. Formation-containment control of networked Euler-Lagrange systems: An event-triggered framework. ISA Trans. 2019, 86, 87–89. [Google Scholar]
  6. Dan, X.; Tan, A.H. Cooperative reinforcement learning in topology-based multi-agent systems. Auton. Agents Multi-Agent Syst. 2013, 26, 86–119. [Google Scholar]
  7. Li, W. The confinement-escape problem of a defender against an evader escaping from a circular region. IEEE Trans. Cybern. 2016, 46, 1028–1039. [Google Scholar] [CrossRef]
  8. Ramana, M.V.; Kothari, M. Pursuit-Evasion Games of High Speed Evader. J. Intell. Robot. Syst. 2016, 85, 293–306. [Google Scholar] [CrossRef]
  9. Ge, H.; Chen, G.; Xu, G. Multi-AUV cooperative target hunting based on Improved potential field in a surface-water environment. Appl. Sci. 2018, 8, 973. [Google Scholar] [CrossRef]
  10. Zhou, L.; Pei, Y.; Chen, C.; Yang, G. Multi-agent reinforcement learning with sparse interactions by negotiation and knowledge transfer. IEEE Trans. Cybern. 2017, 47, 1238–1250. [Google Scholar] [CrossRef]
  11. Sun, Z.; Sun, H.; Li, P.; Zou, J. Self-Organizing cooperative pursuit strategy for multi-USV with dynamic obstacle ships. J. Mar. Sci. Eng. 2022, 10, 562. [Google Scholar] [CrossRef]
  12. Chen, M.; Zhu, D. A Novel Cooperative Hunting Algorithm for Inhomogeneous Multiple Autonomous Underwater Vehicles. IEEE Access 2018, 6, 7818–7828. [Google Scholar] [CrossRef]
  13. Huang, Z.; Zhu, D.; Sun, B. A multi-AUV cooperative hunting method in 3-D underwater environment with obstacle. Eng. Appl. Artif. Intell. 2016, 50, 192–200. [Google Scholar] [CrossRef]
  14. Miao, Z.; Wang, Y.; Fierro, R. Cooperative circumnavigation of a moving target with multiple nonholonomic robots using backstepping design. Syst. Control Lett. 2017, 103, 58–65. [Google Scholar] [CrossRef]
  15. Gao, Z.; Guo, G. Velocity free leader-follower formation control for autonomous underwater vehicles with line-of-sight range and angle constraints. Inf. Sci. 2019, 486, 359–378. [Google Scholar] [CrossRef]
  16. Tan, G.; Zhuang, J.; Zou, J.; Wan, L. Coordination control for multiple unmanned surface vehicles using hybrid behavior-based method. Ocean Eng. 2021, 232, 109147. [Google Scholar] [CrossRef]
  17. Wang, X.; Feng, K.; Wang, G.; Wang, Q. Local path optimization method for unmanned ship based on particle swarm acceleration calculation and dynamic optimal control. Appl. Ocean Res. 2021, 110, 102588. [Google Scholar] [CrossRef]
  18. Sun, Z.; Sun, H.; Li, P.; Zou, J. Formation Control of Multiple Underactuated Surface Vessels with a Disturbance Observer. J. Mar. Sci. Eng. 2022, 10, 1016. [Google Scholar] [CrossRef]
  19. Liu, C.; Xiang, X.; Huang, J.; Yang, S.; Shaoze, Z.; Su, X.; Zhang, Y. Development of Usv Autonomy: Architecture, Implementation and Sea Trials. Brodogradnja 2022, 73, 89–107. [Google Scholar] [CrossRef]
  20. Zhang, G.; Huang, C.; Li, J.; Zhang, X. Constrained coordinated path-following control for underactuated surface vessels with the disturbance rejection mechanism. Ocean Eng. 2020, 196, 106725. [Google Scholar] [CrossRef]
  21. Skjetne, R.; Fossen, T.I.; Kokotović, P.V. Adaptive maneuvering, with experiments, for a model ship in a marine control laboratory. Automatica 2005, 41, 289–298. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of guidance method.
Figure 1. Schematic diagram of guidance method.
Applsci 12 09120 g001
Figure 2. Coordinate systems of USV.
Figure 2. Coordinate systems of USV.
Applsci 12 09120 g002
Figure 3. The trajectories of pursuer and evader.
Figure 3. The trajectories of pursuer and evader.
Applsci 12 09120 g003
Figure 4. Filter performance.
Figure 4. Filter performance.
Applsci 12 09120 g004
Figure 5. Heading angles of pursuer and evader.
Figure 5. Heading angles of pursuer and evader.
Applsci 12 09120 g005
Figure 6. Distances between pursuer and evader.
Figure 6. Distances between pursuer and evader.
Applsci 12 09120 g006
Figure 7. The estimated value of uncertain unknown τ w u .
Figure 7. The estimated value of uncertain unknown τ w u .
Applsci 12 09120 g007
Figure 8. The estimated value of uncertain unknown τ w r .
Figure 8. The estimated value of uncertain unknown τ w r .
Applsci 12 09120 g008
Figure 9. Composition diagram of the USV.
Figure 9. Composition diagram of the USV.
Applsci 12 09120 g009
Figure 10. Structure of UWB position measurement system.
Figure 10. Structure of UWB position measurement system.
Applsci 12 09120 g010
Figure 11. Experimental results of the pursuit trajectory. (a) The pursuer sensed the evader and began to intercept it. (b) The pursuer gradually approached the evader. (c) The pursuer successfully intercepted the evader.
Figure 11. Experimental results of the pursuit trajectory. (a) The pursuer sensed the evader and began to intercept it. (b) The pursuer gradually approached the evader. (c) The pursuer successfully intercepted the evader.
Applsci 12 09120 g011
Table 1. Actual values of CS2.
Table 1. Actual values of CS2.
ParametersValueParametersValue
m 11 1.1 m ^ 11 d 11 1.1 d ^ 11 + 0.2 d ^ 11 | u |
m 22 1.1 m ^ 22 d 22 1.1 d ^ 22 + 0.2 d ^ 22 | v |
m 23 1.1 m ^ 23 d 23 1.1 d ^ 23
m 33 1.1 m ^ 33 d 33 1.1 d ^ 33 + 0.2 d ^ 33 | r |
Table 2. USV control parameters.
Table 2. USV control parameters.
ParametersValueParametersValue
k e u 10 λ φ 0.5
k e v 10 k τ r 2
k e r 10 ξ w u 0.1
k w r 2 ξ w v 0.1
k w u 2 ξ w r 0.1
k u 10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Z.; Sun, H.; Li, P.; Li, X.; Du, L. Pursuit–Evasion Problem of Unmanned Surface Vehicles in a Complex Marine Environment. Appl. Sci. 2022, 12, 9120. https://doi.org/10.3390/app12189120

AMA Style

Sun Z, Sun H, Li P, Li X, Du L. Pursuit–Evasion Problem of Unmanned Surface Vehicles in a Complex Marine Environment. Applied Sciences. 2022; 12(18):9120. https://doi.org/10.3390/app12189120

Chicago/Turabian Style

Sun, Zhiyuan, Huawei Sun, Ping Li, Xiaowen Li, and Lei Du. 2022. "Pursuit–Evasion Problem of Unmanned Surface Vehicles in a Complex Marine Environment" Applied Sciences 12, no. 18: 9120. https://doi.org/10.3390/app12189120

APA Style

Sun, Z., Sun, H., Li, P., Li, X., & Du, L. (2022). Pursuit–Evasion Problem of Unmanned Surface Vehicles in a Complex Marine Environment. Applied Sciences, 12(18), 9120. https://doi.org/10.3390/app12189120

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