Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Multi-AUV Cooperative Underwater Multi-Target Tracking Based on Dynamic-Switching-enabled Multi-Agent Reinforcement Learning

Shengbo Wang, Chuan Lin, , Guangjie Han, , Shengchao Zhu, Zhixian Li, Zhenyu Wang Corresponding author: Chuan Lin & Guangjie HanShengbo Wang and Shengchao Zhu contributed equally to this research.Shengbo Wang, Zhixian Li and Zhenyu Wang are with Software College, Northeastern University, Shenyang, China (e-mails: puppytagge@gmail.com; liaozhaixian@gmail.com; 15063371639@163.com).Chuan Lin is with Software College, Northeastern University, Shenyang, China and is also with Key Laboratory of Data Analytics and Optimization for Smart Industry (Northeastern University), Ministry of Education, China (e-mails: chuanlin1988@gmail.com).Guangjie Han and Shengchao Zhu are with the Department of Internet of Things Engineering, Hohai University, Changzhou, 213022, China. (e-mails: hanguangjie@gmail.com; zhushengchao77@gmail.com).
Abstract

With the rapid development of underwater communication, sensing, automation, robot technologies, autonomous underwater vehicle (AUV) swarms are gradually becoming popular and have been widely promoted in ocean exploration and underwater tracking or surveillance, etc. However, the complex underwater environment poses significant challenges for AUV swarm-based accurate tracking for the underwater moving targets. In this paper, we aim at proposing a multi-AUV cooperative underwater multi-target tracking algorithm especially when the real underwater factors are taken into account. We first give normally modelling approach for the underwater sonar-based detection and the ocean current interference on the target tracking process. Then, we regard the AUV swarm as a underwater ad-hoc network and propose a novel Multi-Agent Reinforcement Learning (MARL) architecture towards the AUV swarm based on Software-Defined Networking (SDN). It enhances the flexibility and scalability of the AUV swarm through centralized management and distributed operations. Based on the proposed MARL architecture, we propose the ”dynamic-attention switching” and ”dynamic-resampling switching” mechanisms, to enhance the efficiency and accuracy of AUV swarm cooperation during task execution. Finally, based on a proposed AUV classification method, we propose an efficient cooperative tracking algorithm called ASMA. Evaluation results demonstrate that our proposed tracking algorithm can perform precise underwater multi-target tracking, comparing with many of recent research products in terms of convergence speed and tracking accuracy.

Index Terms:
AUV swarm network, multi-target tracking, dynamic-switching attention, dynamic-switching resampling, software-defined networking.

I Introduction

Covering more than 70% of Earth’s surface, the ocean is a vital life-support system with abundant biological and mineral resources [1, 2, 3]. With the advancement of technology and our increasing knowledge of the oceans, the need for efficient and intelligent exploration technologies, particularly in ocean environment surveillance [4, 5], disaster prevention, and rescue operations [6, 7], becomes increasingly pressing. In these years, the rapid development of underwater robots, underwater communication, underwater sensing technologies, has given birth to the Autonomous Underwater Vehicle (AUV). As of today, AUV has become indispensable tools in a wide range of marine applications especially in underwater target tracking [8, 9], due to its ability to autonomously navigate underwater environments.

With the continuous improvement and development of swarm intelligence theory, the control technology of AUV swarms has been greatly propelled and enhanced [10, 11]. The AUV swarms have demonstrated unique advantages in executing tasks in complex marine environments, e.g., the AUV-swarm-based cooperative underwater data collection [12], cooperative underwater target tracking [13]. Through cooperation manner, the AUVs in AUV swarm can cover broader areas, enhancing the efficiency and accuracy of underwater operatiion. Particularly in underwater cooperative multi-target tracking, the AUV swarms have been demonstrated exceptional capabilities in cooperatively tracking various underwater targets such as underwater biological communities [14] or environmental pollutants [15]. This is crucial not only for underwater scientific research such as ecological [16] surveys but also for various categories of civil [17] and military applications [18], such as supporting maritime search, rescue operations and maritime boundary surveillance, etc.

However, effectively managing and cooperating such a complex multi-agent system, especially in dynamic and unpredictable marine environments, has to face many challenges. Particularly in underwater multi-target tracking, these challenges mainly manifest in several aspects: The complexity and dynamism of the marine environment, such as ocean currents interference [19, 20], this demands AUV swarms to possess high adaptability and flexibility to cope with sudden environmental changes and unknown underwater obstacles; Effective communication is crucial for guaranteeing the AUV swarm-based underwater cooperative operations, but the high latency and low bandwidth characteristics of underwater communication especially underwater acoustic-based communication limit the efficiency of the information exchange among the AUVs in AUV swarm and cannot support flexible swarm cooperation [21, 22, 23]. Thus, it requires a scalable network architecture to re-organize and re-define the information exchange architecture for AUV swarm system; Meanwhile, in reality, due to productivity constraints, AUV swarms often comprise various heterogeneous AUVs, making AUV swarm-based cooperative underwater operation more challenging and impeding support for cooperative swarm intelligence [24]. Further, traditional MARL algorithms suffer from long runtimes, slow convergence speeds, and high resource consumption, making MARL challenging to deploy in large-scale AUV swarm.

As a cutting-edge networking approach, Software-Defined Networking (SDN) separates the control plane from the data plane in traditional network architectures [25, 26, 27, 28]. By this architecture, the network control functions are centralized into one or more controllers. This allows the network resources together with the network information to be flexibly managed and utilized. Motivated by this, in this work, we display how to utilize SDN technology to re-define the AUV swarm as a AUV swarm network with self-learning ability. For dedicated MARL architecture towards AUV swarm network, we propose a MARL algorithm incorporating a dynamic-switching mechanism, which demonstrates faster convergence and better performance in multi-target tracking tasks under ocean current interference. Totally, this paper mainly makes the following contributions:

  1. 1.

    We regard the AUV swarm as a AUV swarm network, utilize SDN to optimize AUV swarm network and enhance the flexibility and scalability of the AUV swarm;

  2. 2.

    We propose a multi-agent reinforcement learning algorithm incorporating dynamic-switching mechanism (including both ”dynamic-attention switching” and ”dynamic-resampling switching”) to enhance AUV swarm network’s self-learning efficiency and accuracy;

  3. 3.

    We introduces an AUV formation classification algorithm based on fuzzy logic and rule-based expert systems, to enable the AUV swarm network to perform precise tracking of multiple targets in the interference of ocean currents.

The rest of this paper is organized as follows: Related works on the paper’s theme is surveyed in Sec. II; Preliminary materials are introduced in Sec. III; The proposed FATHOM-Net is detailed in Sec. IV; The proposed cooperative tracking algorithm is presented in Sec. V; Evaluation results are shown in Sec. VI; The paper concludes in Sec. VII with a discussion of some future research directions derived from this paper.

II Related Works

In this section, we investigate recent advances related to the main research topics, specifically categorized into the following three areas: 1) SDN-based underwater multi-agent system, 2) Single AUV-based target tracking, and 3) AUV swarm-based target tracking.

II-A SDN-based underwater multi-agent system

SDN technology, characterized by its network architecture that separates the control plane from the data forwarding plane, provides strong support for centralized control and flexible management of underwater multi-agent systems. This technology is currently widely applied in fields such as underwater multi-agent networks, opening new possibilities for underwater operations and research.

For instance, in [29], the authors propose the DRSIR model, a deep reinforcement learning approach for routing selection in SDN. This method utilizes path state indicators and demonstrates higher efficiency and intelligence compared to traditional Dijkstra and RSIR algorithms. It can dynamically adjust routing strategies based on changes in network traffic, showcasing practical feasibility and superior performance in SDN routing.

In [30], an SDN-based architecture for AUV underwater wireless networks is introduced to support multi-AUV cooperative search. This architecture achieves network information synchronization, node localization, multi-AUV cooperative control, and intelligent data transmission scheduling through software-defined beaconing, hierarchical localization, cooperative control, and software-defined hybrid data transmission frameworks. Building upon these studies, we leverage SDN technology to enhance the underwater vehicle swarm network system.

II-B Single AUV-based target tracking

In Single AUV-based target tracking, the focus is on enhancing the perception [31], decision-making [32], and control capabilities [33] of an individual AUV to autonomously perform identification and tracking of specific targets. This approach necessitates highly integrated sensor systems and advanced data processing algorithms to adapt to complex underwater environments and achieve high-precision tracking.

In [33], the authors develop a DRL-based control strategy for X-rudder AUVs, utilizing the DDPG algorithm to enable precise posture control and efficient target tracking. This approach achieves precise three-degree-of-freedom posture control and rapid DRL algorithm deployment, and significantly enhances AUV maneuverability for dynamic underwater target tracking.

In [32], the authors introduce an optimized AUV path planning method (PPM-BBD) employing L-SHADE, focusing on energy-efficient diving through decision-making under motion constraints. This approach minimizes energy consumption by up to 9% compared to conventional methods, demonstrating the importance of strategic decision-making in solitary AUV target tracking amidst environmental challenges.

II-C AUV swarm-based target tracking

AUV swarm-based target tracking leverages the coordinated efforts of multiple AUVs, utilizing principles of swarm intelligence to improve tracking range, efficiency, and robustness. This method not only extends beyond the limitations of a single AUV but also emphasizes distributed decision-making and information sharing among the AUVs to effectively track multiple or wide-area targets in dynamic environments [34, 35].

In [34], the authors propose a novel multi-AUV cooperative tracking scheme named CTDE, which employs a central training and distributed execution approach to achieve safe and effective tracking of moving targets using MADDPG. By conducting central training within a designed secure private network, information sharing is not required during task execution, thereby enhancing the security of the entire system.

In [35], the authors discuss trajectory tracking control of multiple AUVs using discrete-time control methods in weak communication environments. The study focuses on considering communication delays (bounded and unbounded) and packet loss scenarios, ensuring stable tracking of given trajectories by designing methods based on leader-follower and virtual leader approaches.

Totally, multi-AUV-based underwater target tracking technology has received widespread research attention. However, this research field still faces several challenges, such as adaptability to dynamic and complex marine environments, flexibility and robustness of AUV swarm networks, and how to stably and accurately track multiple targets. In response to these issues, our work employ the SDN technology to improve the scalability and functionability of the AUV swarm system. On account of the above system architecture, we propose a dynamic-switching-enabled MARL architecture which is more feasible for AUV swarm system’s self-learning. And, we propose as well a novel AUV classification method to improve the scalability of the AUV swarm-based tracking system. In the following, we will prove our work demonstrates precise multi-target tracking in dynamic and complex marine environments, compared to many existing studies.

III Preliminary Materials

In order to accurately track underwater targets and adapt to the dynamic factors in real underwater environments, we consider the effects of sonar detection and ocean current interference. Furthermore, we transform the target tracking scenario based on AUV swarm network into a Markov Decision Process (MDP).

III-A Sonar detection modeling

In real complex underwater environments, traditional electromagnetic-based detection methods are limited in their effectiveness. Therefore, in this paper, we utilize sonar technology to detect AUVs’ position together with their tracking targets. Equipped with sonar, AUVs emit directional sound waves to scan their surrounding environment. This process involves using a sector array sonar to capture echo signals from different directions, and the variations in signal intensity allow us to estimate the position of targets. In this study, a specific active sonar equation model is employed to simulate and optimize the target detection process, which can be expressed by Eq. 1.

EM=SL2TL+TS(NLDI)DT𝐸𝑀𝑆𝐿2𝑇𝐿𝑇𝑆𝑁𝐿𝐷𝐼𝐷𝑇EM=SL-2TL+TS-(NL-DI)-DTitalic_E italic_M = italic_S italic_L - 2 italic_T italic_L + italic_T italic_S - ( italic_N italic_L - italic_D italic_I ) - italic_D italic_T (1)

where SL𝑆𝐿SLitalic_S italic_L is the strength of the sonar signal emitted, TL𝑇𝐿TLitalic_T italic_L is the attenuation of sound propagation, TS𝑇𝑆TSitalic_T italic_S is the target echo strength, NL𝑁𝐿NLitalic_N italic_L is the ambient marine noise, DI𝐷𝐼DIitalic_D italic_I is the system’s noise suppression capability, and DT𝐷𝑇DTitalic_D italic_T is the minimum signal-to-noise ratio required for the equipment to function properly. The units for the notations mentioned are all in decibels (dB).

III-B Ocean current modeling

The underwater environments’ uncertainty, especially the influence of ocean currents, poses significant challenges to the navigation and positioning capabilities of AUVs. It demands AUVs to possess a high level of adaptability and intelligence.

The Navier-Stokes equations are the fundamental equations for describing fluid dynamics behavior. They can simulate the characteristics of fluid flow. Additionally, they provide a computational environment for calculating the forces exerted by ocean currents (Eq. 2).

ρ(𝐮t+𝐮𝐮)=p+μ2𝐮+𝐅𝜌𝐮𝑡𝐮𝐮𝑝𝜇superscript2𝐮𝐅\rho\left(\frac{\partial\mathbf{u}}{\partial t}+\mathbf{u}\cdot\nabla\mathbf{u% }\right)=-\nabla p+\mu\nabla^{2}\mathbf{u}+\mathbf{F}italic_ρ ( divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG + bold_u ⋅ ∇ bold_u ) = - ∇ italic_p + italic_μ ∇ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT bold_u + bold_F (2)

where ρ𝜌\rhoitalic_ρ is fluid density, 𝐮𝐮\mathbf{u}bold_u is velocity, 𝐮t𝐮𝑡\frac{\partial\mathbf{u}}{\partial t}divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG is velocity change over time, 𝐮𝐮𝐮𝐮\mathbf{u}\cdot\nabla\mathbf{u}bold_u ⋅ ∇ bold_u is self-advection, p𝑝\nabla p∇ italic_p is pressure change, μ𝜇\muitalic_μ is viscosity, 2𝐮superscript2𝐮\nabla^{2}\mathbf{u}∇ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT bold_u is the diffusion term, and 𝐅𝐅\mathbf{F}bold_F is external forces.

However, the Navier-Stokes equations involve a large number of partial differential calculations, resulting in high computational costs and long processing times. To reduce the complexity of the model and obtain preliminary fluid dynamic characteristics, we assume that fluid flow is laminar and relatively uniform. Neglecting the nonlinear convective term 𝐮𝐮𝐮𝐮\mathbf{u}\cdot\nabla\mathbf{u}bold_u ⋅ ∇ bold_u in the Navier-Stokes equations simplifies it to a linear problem that can be expressed by the following.

ρ𝐮t=p+μ2𝐮+𝐅𝜌𝐮𝑡𝑝𝜇superscript2𝐮𝐅\rho\frac{\partial\mathbf{u}}{\partial t}=-\nabla p+\mu\nabla^{2}\mathbf{u}+% \mathbf{F}italic_ρ divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG = - ∇ italic_p + italic_μ ∇ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT bold_u + bold_F (3)

After calculating 𝐮𝐮\mathbf{u}bold_u and 𝐮t𝐮𝑡\frac{\partial\mathbf{u}}{\partial t}divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG through the Navier-Stokes equations, we use the drag equation (Eq. 4), lift equation (Eq. 5), and virtual mass equation (Eq. 6) to calculate the forces acting on the structure under the influence of ocean currents:

FD=12ρu2CDAsubscript𝐹𝐷12𝜌superscript𝑢2subscript𝐶𝐷𝐴F_{D}=\frac{1}{2}\rho u^{2}C_{D}Aitalic_F start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ρ italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT italic_A (4)

where FDsubscript𝐹𝐷F_{D}italic_F start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT is the drag force, ρ𝜌\rhoitalic_ρ is fluid density, u𝑢uitalic_u is relative fluid velocity, CDsubscript𝐶𝐷C_{D}italic_C start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT is drag coefficient, and A𝐴Aitalic_A is object’s frontal area.

FL=12ρu2CLAsubscript𝐹𝐿12𝜌superscript𝑢2subscript𝐶𝐿𝐴F_{L}=\frac{1}{2}\rho u^{2}C_{L}Aitalic_F start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ρ italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_A (5)

where FLsubscript𝐹𝐿F_{L}italic_F start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT is the lift force, CLsubscript𝐶𝐿C_{L}italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT is lift coefficient.

FVM=ρCVMV𝐮tsubscript𝐹𝑉𝑀𝜌subscript𝐶𝑉𝑀𝑉𝐮𝑡F_{VM}=\rho C_{VM}V\frac{\partial\mathbf{u}}{\partial t}italic_F start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT = italic_ρ italic_C start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT italic_V divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG (6)

where FVMsubscript𝐹𝑉𝑀F_{VM}italic_F start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT is the virtual mass force, CVMsubscript𝐶𝑉𝑀C_{VM}italic_C start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT is the virtual mass coefficient, V𝑉Vitalic_V is AUV’s volume, and ut𝑢𝑡\frac{\partial u}{\partial t}divide start_ARG ∂ italic_u end_ARG start_ARG ∂ italic_t end_ARG is the rate of change of velocity.

Through the above equations (Eq. 4-Eq. 6), we can compute the complex and variable forces exerted by ocean currents on the AUV. The three forces are shown in Fig. 1. The overall calculation for the ocean current can be summarized in Algorithm 1.

Refer to caption
Figure 1: Schematic diagram of ocean current
Algorithm 1 Simplified fluid dynamics calculation algorithm
1:Initialize the Set NEsubscript𝑁𝐸N_{E}italic_N start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT, T𝑇Titalic_T respectively.
2:Initialize AUV’s frontal area A𝐴Aitalic_A, volume V, drag coefficient CDsubscript𝐶𝐷C_{D}italic_C start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT, lift coefficient CLsubscript𝐶𝐿C_{L}italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT, and inertia coefficient CVMsubscript𝐶𝑉𝑀C_{VM}italic_C start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT.
3:for Episode=1Episode1\text{Episode}=1Episode = 1 to NEsubscript𝑁𝐸N_{E}italic_N start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT do
4:     Initialize fluid density ρ𝜌\rhoitalic_ρ, dynamic viscosity μ𝜇\muitalic_μ, and external forces 𝐅𝐅\mathbf{F}bold_F.
5:     Initialize environment and fluid velocity 𝐮𝐮\mathbf{u}bold_u and pressure gradient p𝑝\nabla p∇ italic_p.
6:     for t=1𝑡1t=1italic_t = 1 to T𝑇Titalic_T do
7:         Calculate the velocity time derivative 𝐮t𝐮𝑡\frac{\partial\mathbf{u}}{\partial t}divide start_ARG ∂ bold_u end_ARG start_ARG ∂ italic_t end_ARG by simplified Navier-Stokes equation.
8:         Update fluid velocity 𝐮𝐮\mathbf{u}bold_u and pressure gradient p𝑝\nabla p∇ italic_p.
9:         Calculate the drag force FDsubscript𝐹𝐷F_{D}italic_F start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT, lift force FLsubscript𝐹𝐿F_{L}italic_F start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT, and virtual mass force FVMsubscript𝐹𝑉𝑀F_{VM}italic_F start_POSTSUBSCRIPT italic_V italic_M end_POSTSUBSCRIPT by Eq. 4 to Eq. 6.
10:         Combine forces to determine net force on the AUV.      

III-C Markov decision process modeling

In the decision-making process of AUV swarm network, the environment is defined as a MDP (Eq. 7), consisting of five main components: the state set 𝒮𝒮\mathcal{S}caligraphic_S, the action set 𝒜𝒜\mathcal{A}caligraphic_A, the state transition probabilities 𝒫𝒫\mathcal{P}caligraphic_P, the reward function \mathcal{R}caligraphic_R, and the discount factor γ𝛾\gammaitalic_γ:

=(𝒮,𝒜,𝒫,,γ)𝒮𝒜𝒫𝛾\mathcal{M}=(\mathcal{S},\mathcal{A},\mathcal{P},\mathcal{R},\gamma)caligraphic_M = ( caligraphic_S , caligraphic_A , caligraphic_P , caligraphic_R , italic_γ ) (7)

The state set 𝒮𝒮\mathcal{S}caligraphic_S encompasses the sonar detection results in different coordinate systems and their location information. Specifically, it is a set composed of the target and the AUV’s echo intensity along with their location information. 𝒮𝒮\mathcal{S}caligraphic_S can be represented as 𝒮={s1,,sN}𝒮subscript𝑠1subscript𝑠𝑁\mathcal{S}=\{s_{1},\ldots,s_{N}\}caligraphic_S = { italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_s start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT }, where si=(ηi,ϕi,oi)subscript𝑠𝑖subscript𝜂𝑖subscriptitalic-ϕ𝑖subscript𝑜𝑖s_{i}=(\eta_{i},\phi_{i},o_{i})italic_s start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_η start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_ϕ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_o start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ), oi=(κi,σi)Nκ+Nσsubscript𝑜𝑖subscript𝜅𝑖subscript𝜎𝑖superscriptsubscript𝑁𝜅subscript𝑁𝜎o_{i}=(\kappa_{i},\sigma_{i})\in\mathbb{R}^{N_{\kappa}+N_{\sigma}}italic_o start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ∈ blackboard_R start_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_κ end_POSTSUBSCRIPT + italic_N start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT end_POSTSUPERSCRIPT, κisubscript𝜅𝑖\kappa_{i}italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT represents the echo intensity from the target, σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT represents the echo intensity from other AUVs, Nκsubscript𝑁𝜅N_{\kappa}italic_N start_POSTSUBSCRIPT italic_κ end_POSTSUBSCRIPT is the number of targets being tracked, and Nσsubscript𝑁𝜎N_{\sigma}italic_N start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT is the number of other AUVs. The action set 𝒜𝒜\mathcal{A}caligraphic_A in this study represents the discrete action set of all AUVs, denoted as 𝒜={a1,,aN}𝒜subscript𝑎1subscript𝑎𝑁\mathcal{A}=\{a_{1},\ldots,a_{N}\}caligraphic_A = { italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_a start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT }, where each action aisubscript𝑎𝑖a_{i}italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is a vector representing seven possible actions (e.g., up, down, left, right, forward, backward, and stay in place). The state transition probabilities 𝒫𝒫\mathcal{P}caligraphic_P describe the dynamic transition relationship between states, which is key for the algorithm to understand how different states transition into each other. The discount factor γ𝛾\gammaitalic_γ, where the value of γ𝛾\gammaitalic_γ ranges between 0 and 1, together determines the calculation method for long-term benefits in policy evaluation. Further, the reward function \mathcal{R}caligraphic_R will be explained in Sec. V, where we consider tracking accuracy, collision avoidance, ocean current stability, and energy efficiency to enhance multi-AUV cooperative tracking effects.

IV Hierarchical Software-Defined Multi-AUV Reinforcement Learning Architecture

In this section, we introduce Flexible Autonomous Training Hierarchy for Operational Marine Network (FATHOM-Net), which is a software-defined MARL architecture based on Centralized Training and Decentralized Execution.

Refer to caption
Figure 2: Proposed FATHOM-Net

IV-A Overview of FATHOM-Net

FATHOM-Net framework consists of three functional layers: global-control training layer, local-control training layer, and application execution tracking layer, which can be detailed as the following.

Global-control training layer: The global-control training layer comprises a global controller based on Unmanned Surface Vessel-based Global Controller (USV-GC). It aims to coordinate Local Controller-based AUV (LC-AUV)s in various regions through a dedicated northbound interface, ensuring efficient classification of large-scale tasks to suitable local-control training layers by deconstructing and optimizing tasks. The strategy updates and task scheduling mechanism in this layer not only support the efficient execution of complex tasks but also continuously optimize the global network view by periodically obtaining status information from the lower layers. This ensures the flexibility and efficiency of network operations, especially when facing environmental challenges and changing task requirements.

Local-control training layer: It focuses on optimizing tasks issued by the global-control training layer into specific strategies through training, and conveying them to each managed Execution Tracking-based AUV (ET-AUV). Upon assignment of subregion tasks by the global-control training layer, each LC-AUV begins centralized local training. This process uses the unique characteristics of the ET-AUV subset it manages, following the parameters and tasks provided by the global-control training layer. This further refines tasks and formulates specific execution strategies to guide its managed ET-AUVs. Additionally, it periodically maintains the local network view through a dedicated southbound interface. As LC-AUVs require a certain amount of computation, they are typically equipped with high-performance batteries and computing units. This training mode ensures that all ET-AUVs can execute tasks according to uniform standards and optimal strategies, while also allowing for rapid adjustments tailored to local environments and specific task requirements.

Application execution tracking layer: It is directly responsible for executing specific strategies. This layer consists of ordinary ET-AUVs swarms. Each ET-AUV independently executes tasks and strategies received from the local-control training layer. Therefore, communication among ET-AUVs is not imposed. This idea (the idea of ”separate control from operation” in SDN) reduces communication costs and enhances the flexibility and robustness of the AUV swarm network. Even if individual ET-AUVs fails or are destroyed, it does not affect the overall efficiency of the network.

Overall, FATHOM-Net reduces the coupling between individual AUVs, addressing the scalability issues present in traditional architectures. The FATHOM-Net framework is summarized in Fig. 2. The entire task classification processing is summarized in Fig. 3.

IV-B Beacon framework for FATHOM-Net

In order to ensure precise target tracking based on the proposed FATHOM-Net, it is crucial to establish both local and global network views for LC-AUVs and USV-GC. Therefore, we propose an SDN beacon framework to synchronize and share information among ET-AUVs, LC-AUVs, and USV-GC. The proposed SDN beacon framework can be carried out by the following:

Phase 1: To update local and global network views, the beaconing is activated between USV-GC/LC-AUVs and their subordinate levels within the region requesting synchronization (SYN) information. Upon receiving a synchronization request, each LC-AUV/ET-AUV replies with key information such as its AUV ID and status to update the network views.

Phase 2: To allocate specific-tasks/strategies to LC-AUVs/ET-AUVs through centralized training. In each training round, USV-GC/LC-AUVs send requests to LC-AUVs/ET-AUVs they manage after training completion, and send specific-tasks/strategies to them through operation requests. Each LC-AUV/ET-AUV then replies with key information.

In summary, Phase 1 describes the network view updating procedure, while Phase 2 displays the task classification processing. To improve the information exchange efficiency, the reply information in Phase 2 also includes information for updating the network views. Therefore, these replies are also considered as an update to the network views.

V Proposed Cooperative Tracking Algorithm

In this section, we showcase our proposed cooperative tracking algorithm. Especially, we firstly propose a classification algorithm called AUV Scalable Management and Allocation (ASMA) to classify/allocate ET-AUVs to each formation for performing cooperative multi-target tracking. To enhance the tracking accuracy and the convergence speed of the AUV swarm network, we employ the concept of ”Dynamic-Switching”, integrate the ”Dynamic-Switching Attention” mechanism into the Actor-Critic model and develop the Dymiac-Switching-Based MARL (DSBM) architecture. Additionally, to improve the efficiency of better sample utilization in DSBM, we introduce a novel ”Dynamic-Switching Resampling” mechanism based on experience replay buffer mechanism.

V-A AUV classification method

In underwater multi-target tracking problems, efficiently and reasonably classifying tracking formations (consisting of a series of ET-AUVs) for each underwater target is an important issue. Totally, the proposed ASMA achieves efficient and accurate formation assignments for ET-AUVs by utilizing a combination of fuzzy logic and rule-based expert systems, integrating various performance indicators.

The evaluation process initially employs fuzzy logic scoring to map specific performance indicators of ET-AUVs onto a membership degree. The mathematical expression for this function is as follows (Eq. 8):

Refer to caption
Figure 3: Task classification based on FATHOM-Net
TABLE I: Special notations in Sec. V
Notation Description
αisubscript𝛼𝑖\alpha_{i}italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Acceleration score
βisubscript𝛽𝑖\beta_{i}italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Remaining battery score
κisubscript𝜅𝑖\kappa_{i}italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Carrying capacity score
ϵisubscriptitalic-ϵ𝑖\epsilon_{i}italic_ϵ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Energy consumption score
di,jsubscript𝑑𝑖𝑗d_{i,j}italic_d start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT Distance between ET-AUV i𝑖iitalic_i to the target j𝑗jitalic_j
Disubscript𝐷𝑖D_{i}italic_D start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Perception range
visubscript𝑣𝑖v_{i}italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT Speed score
μ(x)={0if x<range_min1if x>range_maxxrange_minrange_maxrange_minotherwise𝜇𝑥cases0if 𝑥range_min1if 𝑥range_max𝑥range_minrange_maxrange_minotherwise\mu(x)=\begin{cases}0&\text{if }x<\text{range\_min}\\ 1&\text{if }x>\text{range\_max}\\ \frac{x-\text{range\_min}}{\text{range\_max}-\text{range\_min}}&\text{% otherwise}\end{cases}italic_μ ( italic_x ) = { start_ROW start_CELL 0 end_CELL start_CELL if italic_x < range_min end_CELL end_ROW start_ROW start_CELL 1 end_CELL start_CELL if italic_x > range_max end_CELL end_ROW start_ROW start_CELL divide start_ARG italic_x - range_min end_ARG start_ARG range_max - range_min end_ARG end_CELL start_CELL otherwise end_CELL end_ROW (8)

where μ(x)𝜇𝑥\mu(x)italic_μ ( italic_x ) is the membership degree of notation x𝑥xitalic_x normalized between 0 and 1, x𝑥xitalic_x is the actual value of the notation, and range_min and range_max set the normalization range for x𝑥xitalic_x, specifying the minimum and maximum values, respectively.

The rule-based expert system takes the fuzzy membership degrees as input. It then calculates the overall score Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT acquired by agent i𝑖iitalic_i and target j𝑗jitalic_j. The notations and their meanings are shown in Table I.

Thus, based on ASMA, the score Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT acquired by ET-AUV i𝑖iitalic_i and underwater tracking target j𝑗jitalic_j can be computed by the following steps:

Step 1: Energy regulation score:

Si,j+=(vi3+2αi2)×11+eβi×(1di,jDi)limit-fromsubscript𝑆𝑖𝑗superscriptsubscript𝑣𝑖32superscriptsubscript𝛼𝑖211superscript𝑒subscript𝛽𝑖1subscript𝑑𝑖𝑗subscript𝐷𝑖\mathit{S}_{i,j}+=\left(v_{i}^{3}+2\alpha_{i}^{2}\right)\times\frac{1}{1+e^{-% \beta_{i}}}\times\left(1-\frac{d_{i,j}}{D_{i}}\right)italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = ( italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) × divide start_ARG 1 end_ARG start_ARG 1 + italic_e start_POSTSUPERSCRIPT - italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG × ( 1 - divide start_ARG italic_d start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT end_ARG start_ARG italic_D start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ) (9)

Significant increases in total Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT occur when ET-AUVs have fast speeds, stable accelerations, and sufficient power. However, when the remaining power is low, even if the ET-AUV accelerates quickly, the Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT will be lower to avoid that ET-AUVs quickly deplete their power. Additionally, as the distance between the ET-AUV and the target approaches the maximum perception range, the score decreases to prevent losing track of the target.

Step 2: Efficiency load balancing score:

Si,j+=cosh(ϵiκi)×|κiϵi|limit-fromsubscript𝑆𝑖𝑗subscriptitalic-ϵ𝑖subscript𝜅𝑖subscript𝜅𝑖subscriptitalic-ϵ𝑖\mathit{S}_{i,j}+=\cosh(\epsilon_{i}-\kappa_{i})\times\sqrt{\lvert\kappa_{i}-% \epsilon_{i}\rvert}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = roman_cosh ( italic_ϵ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) × square-root start_ARG | italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_ϵ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | end_ARG (10)

By comparing energy consumption with payload capacity, the balance between energy efficiency and stability of ET-AUVs during task execution is emphasized.

Step 3: Comprehensive performance score:

Si,j+=vi2+αi2+ϵi23limit-fromsubscript𝑆𝑖𝑗superscriptsubscript𝑣𝑖2superscriptsubscript𝛼𝑖2superscriptsubscriptitalic-ϵ𝑖23\mathit{S}_{i,j}+=\frac{v_{i}^{2}+\alpha_{i}^{2}+\epsilon_{i}^{2}}{3}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = divide start_ARG italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ϵ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 end_ARG (11)

This rule considers the comprehensive impact of speed, acceleration, and energy consumption, reflecting the stability of ET-AUV’s overall performance under different operating conditions.

Step 4: Dynamic counter-reaction score:

Si,j+=11+eαivilimit-fromsubscript𝑆𝑖𝑗11superscript𝑒subscript𝛼𝑖subscript𝑣𝑖\mathit{S}_{i,j}+=\frac{1}{1+e^{\alpha_{i}-v_{i}}}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = divide start_ARG 1 end_ARG start_ARG 1 + italic_e start_POSTSUPERSCRIPT italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG (12)

If the ET-AUV can smoothly adjust its speed, it indicates high-performance control capability under dynamic conditions. This is crucial for EC-AUVs to perform complex underwater evasive action, e.g., collision avoidance, precise positioning, efficient navigation, etc.

Step 5: Fine operational expertise score:

Si,j+=if v0.3 then v×κi else 0limit-fromsubscript𝑆𝑖𝑗if 𝑣0.3 then 𝑣subscript𝜅𝑖 else 0\mathit{S}_{i,j}+=\text{if }v\leq 0.3\text{ then }v\times\kappa_{i}\text{ else% }0italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = if italic_v ≤ 0.3 then italic_v × italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT else 0 (13)

In underwater tasks requiring precise control, ET-AUVs with low speed and high payload capacity will receive additional Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT due to higher operational accuracy.

Step 6: Remote rapid response score:

Si,j+=if v0.7 then v×βi×di,jDi else 0limit-fromsubscript𝑆𝑖𝑗if 𝑣0.7 then 𝑣subscript𝛽𝑖subscript𝑑𝑖𝑗subscript𝐷𝑖 else 0\mathit{S}_{i,j}+=\text{if }v\geq 0.7\text{ then }v\times\beta_{i}\times\frac{% d_{i,j}}{D_{i}}\text{ else }0italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT + = if italic_v ≥ 0.7 then italic_v × italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × divide start_ARG italic_d start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT end_ARG start_ARG italic_D start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG else 0 (14)

In underwater urgent tasks requiring the ET-AUV to fleetly arrive at the distant targets, high-speed ET-AUVs at greater distances from the target receive additional Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT, reflecting their rapid response capability to distant targets.

Totally, the proposed ASMA is carried out based on the Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT of each ET-AUV for each target. This ensures that every ET-AUV tracks a target and each target is tracked by at least one ET-AUV. The optimization goal is to maximize the cumulative score.

In summary, ASMA can be summarized as Algorithm 2, which considers various special factors in different scenarios, classifying rational and efficient tracking formations.

Algorithm 2 ET-AUV performance classification algorithm
1:Initialize the number k𝑘kitalic_k of ET-AUVs and the number T𝑇Titalic_T of targets.
2:for i=1𝑖1i=1italic_i = 1 to k𝑘kitalic_k do
3:     Initialize ET-AUV i𝑖iitalic_i’s αisubscript𝛼𝑖\alpha_{i}italic_α start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, βisubscript𝛽𝑖\beta_{i}italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, κisubscript𝜅𝑖\kappa_{i}italic_κ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, ϵisubscriptitalic-ϵ𝑖\epsilon_{i}italic_ϵ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, di,jsubscript𝑑𝑖𝑗d_{i,j}italic_d start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT, Disubscript𝐷𝑖D_{i}italic_D start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and visubscript𝑣𝑖v_{i}italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT.
4:for i=1𝑖1i=1italic_i = 1 to k𝑘kitalic_k do
5:     Calculate membership values for ET-AUV i𝑖iitalic_i by Eq. 8.
6:for i=1𝑖1i=1italic_i = 1 to k𝑘kitalic_k do
7:     for j=1𝑗1j=1italic_j = 1 to T𝑇Titalic_T do
8:         Update the score Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT of ET-AUV i𝑖iitalic_i for target j𝑗jitalic_j by Eq. 9 to Eq. 14.      
9:Assign the best formation plan based on score Si,jsubscript𝑆𝑖𝑗\mathit{S}_{i,j}italic_S start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT.

V-B Proposed dynamic-switching-based MARL

Normally, in MARL, there are two major challenges: the high variance problem of traditional strategy estimation methods and the slow convergence problem caused by the high dimensionality of Critic networks in multi-agent systems. To address these issues, in this paper, we introduce Temporal Difference techniques to reduce evaluation variance and propose DSBM, which adopts the ”Dynamic-Switching” mechanism to improve learning efficiency and accelerate policy convergence. The ”Dynamic-Switching” mechanism focuses on optimizing the handling of high-value samples, as shown on the right side of Fig. 2.

Dynamic-switching attention mechanism: During the process of tracking targets, there are a total of N𝑁Nitalic_N ET-AUVs. For a specific ET-AUV i𝑖iitalic_i, we need to calculate its Q-value function Qψi(o,a)subscript𝑄subscript𝜓𝑖𝑜𝑎Q_{\psi_{i}}(o,a)italic_Q start_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_o , italic_a ), where o𝑜oitalic_o represents the current state and a𝑎aitalic_a represents the action. The Qψi(o,a)subscript𝑄subscript𝜓𝑖𝑜𝑎Q_{\psi_{i}}(o,a)italic_Q start_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_o , italic_a ), which represents the value of the action taken by ET-AUV i𝑖iitalic_i given a particular observation and action. It can be obtained through the following steps.

For each ET-AUV, first identify the ET-AUV with the highest reward in the current training round (AUVbe𝐴𝑈𝑉𝑏𝑒AUVbeitalic_A italic_U italic_V italic_b italic_e) (Eq. 15).

AUVbe=argmax(Rf)𝐴𝑈subscript𝑉beargmaxsubscript𝑅𝑓AUV_{\text{be}}=\text{argmax}(R_{f})italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT = argmax ( italic_R start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) (15)

where Rfsubscript𝑅𝑓R_{f}italic_R start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT is the set of reward value for the AUVs at round f𝑓fitalic_f.

After that, it randomly select the experience from the remaining ET-AUVs (Eq. 16).

AUVad={AUVberan(R,2)if i=AUVbeAUVberan(R,1)if iAUVbe𝐴𝑈subscript𝑉adcases𝐴𝑈subscript𝑉beransuperscript𝑅2if 𝑖𝐴𝑈subscript𝑉be𝐴𝑈subscript𝑉beransuperscript𝑅1if 𝑖𝐴𝑈subscript𝑉beAUV_{\text{ad}}=\begin{cases}AUV_{\text{be}}\cup\text{ran}(R^{*},2)&\text{if }% i=AUV_{\text{be}}\\ AUV_{\text{be}}\cup\text{ran}(R^{*},1)&\text{if }i\neq AUV_{\text{be}}\end{cases}italic_A italic_U italic_V start_POSTSUBSCRIPT ad end_POSTSUBSCRIPT = { start_ROW start_CELL italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT ∪ ran ( italic_R start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , 2 ) end_CELL start_CELL if italic_i = italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT ∪ ran ( italic_R start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , 1 ) end_CELL start_CELL if italic_i ≠ italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT end_CELL end_ROW (16)

where AUVad𝐴𝑈subscript𝑉adAUV_{\text{ad}}italic_A italic_U italic_V start_POSTSUBSCRIPT ad end_POSTSUBSCRIPT is the chosen ET-AUVs for update, and ran(X,n)ran𝑋𝑛\text{ran}(X,n)ran ( italic_X , italic_n ) selects n𝑛nitalic_n experience from X𝑋Xitalic_X randomly. The set Rsuperscript𝑅R^{*}italic_R start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT denotes all ET-AUVs excluding AUVbe𝐴𝑈subscript𝑉beAUV_{\text{be}}italic_A italic_U italic_V start_POSTSUBSCRIPT be end_POSTSUBSCRIPT.

The specific selection process is as follows:

  1. 1.

    If the current ET-AUV has the highest reward, then randomly select two from the remaining ET-AUVs to avoid the algorithm getting stuck in local optima.

  2. 2.

    If it doesn’t, then only randomly select one from the other ET-AUVs to increase the proportion of experiences from the best ET-AUV.

After obtaining AUVad𝐴𝑈𝑉𝑎𝑑AUVaditalic_A italic_U italic_V italic_a italic_d, exclude AUVad𝐴𝑈𝑉𝑎𝑑AUVaditalic_A italic_U italic_V italic_a italic_d from the ET-AUV set to obtain AUVord𝐴𝑈𝑉𝑜𝑟𝑑AUVorditalic_A italic_U italic_V italic_o italic_r italic_d (Eq. 17).

AUVord=Rf\AUVad𝐴𝑈subscript𝑉ord\subscript𝑅𝑓𝐴𝑈subscript𝑉adAUV_{\text{ord}}=R_{f}\backslash AUV_{\text{ad}}italic_A italic_U italic_V start_POSTSUBSCRIPT ord end_POSTSUBSCRIPT = italic_R start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT \ italic_A italic_U italic_V start_POSTSUBSCRIPT ad end_POSTSUBSCRIPT (17)

Then, concatenate the state and action information of ET-AUVs in AUVad𝐴𝑈𝑉𝑎𝑑AUVaditalic_A italic_U italic_V italic_a italic_d to form m1𝑚1m1italic_m 1 (Eq. 18).

m1=catkAUVad(ok,ak)subscript𝑚1𝑐𝑎subscript𝑡𝑘𝐴𝑈subscript𝑉𝑎𝑑subscript𝑜𝑘subscript𝑎𝑘m_{1}=cat_{k\in AUV_{ad}}(o_{k},a_{k})italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_c italic_a italic_t start_POSTSUBSCRIPT italic_k ∈ italic_A italic_U italic_V start_POSTSUBSCRIPT italic_a italic_d end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_o start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) (18)

Apply neural networks and dynamic-switching attention weight processing to ET-AUVs in AUVord𝐴𝑈𝑉𝑜𝑟𝑑AUVorditalic_A italic_U italic_V italic_o italic_r italic_d to form m2𝑚2m2italic_m 2 (Eq. 19).

m2=kAUVordωkf(ok,ak)subscript𝑚2subscript𝑘𝐴𝑈subscript𝑉𝑜𝑟𝑑subscript𝜔𝑘𝑓subscript𝑜𝑘subscript𝑎𝑘m_{2}=\sum_{k\in AUV_{ord}}\omega_{k}f(o_{k},a_{k})italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_k ∈ italic_A italic_U italic_V start_POSTSUBSCRIPT italic_o italic_r italic_d end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_f ( italic_o start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) (19)

where the dynamic-switching attention weight ωksubscript𝜔𝑘\omega_{k}italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT, is to evaluate the importance of the information of the other AUVs, and f(ok,ak)𝑓subscript𝑜𝑘subscript𝑎𝑘f(o_{k},a_{k})italic_f ( italic_o start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is a function to extract key features from oksubscript𝑜𝑘o_{k}italic_o start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT and aksubscript𝑎𝑘a_{k}italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT of the k𝑘kitalic_k-th ET-AUV.

ωk=rfkmAUVordrfmsubscript𝜔𝑘subscript𝑟𝑓𝑘subscript𝑚𝐴𝑈subscript𝑉𝑜𝑟𝑑subscript𝑟𝑓𝑚\omega_{k}=\frac{r_{fk}}{\sum_{m\in AUV_{ord}}r_{fm}}italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = divide start_ARG italic_r start_POSTSUBSCRIPT italic_f italic_k end_POSTSUBSCRIPT end_ARG start_ARG ∑ start_POSTSUBSCRIPT italic_m ∈ italic_A italic_U italic_V start_POSTSUBSCRIPT italic_o italic_r italic_d end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_f italic_m end_POSTSUBSCRIPT end_ARG (20)

where rfksubscript𝑟𝑓𝑘r_{fk}italic_r start_POSTSUBSCRIPT italic_f italic_k end_POSTSUBSCRIPT is the reward of k𝑘kitalic_k-th ET-AUV at f𝑓fitalic_f-th round.

Finally, input m1𝑚1m1italic_m 1 and m2𝑚2m2italic_m 2 into a neural network (MLP) to calculate the Q-value function (Eq. 21).

Qiψ(o,a)=fi(m1,m2)subscriptsuperscript𝑄𝜓𝑖𝑜𝑎subscript𝑓𝑖subscript𝑚1subscript𝑚2Q^{\psi}_{i}(o,a)=f_{i}(m_{1},m_{2})italic_Q start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_o , italic_a ) = italic_f start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) (21)

where fisubscript𝑓𝑖f_{i}italic_f start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is a three-layer MLP.

Dynamic-switching resampling mechanism: In our study, to enhance the training efficiency by using better samples, we propose an improved sampling strategy based on the experience replay buffer D𝐷Ditalic_D. The buffer D𝐷Ditalic_D accumulates data from multiple training rounds. This data is stored in the form of a series of tuples Disubscript𝐷𝑖D_{i}italic_D start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, where each tuple records the states, actions, and rewards of all ET-AUVs in a specific round.

Our approach involves extracting a subset D1subscript𝐷1D_{1}italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT comprising n𝑛nitalic_n samples from D𝐷Ditalic_D in each training round, which is then used to update the Critic network. Subsequently, the update process selects the top 50% highest-reward samples from D1subscript𝐷1D_{1}italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. These selected samples form a new set D2subscript𝐷2D_{2}italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, whose samples are then used for a second update for the Critic network. This mechanism ensures that the focus of network learning is on the most valuable experiences.

Network updates: Further, we employ a progressive strategy to update network parameters using the soft update method (Eq. 22), gradually adjusting the target network parameters to stabilize the training process.

ω^τω+(1τ)ω^^𝜔𝜏𝜔1𝜏^𝜔\hat{\omega}\leftarrow\tau\omega+(1-\tau)\hat{\omega}over^ start_ARG italic_ω end_ARG ← italic_τ italic_ω + ( 1 - italic_τ ) over^ start_ARG italic_ω end_ARG (22)

where ω𝜔\omegaitalic_ω represents the parameters of the original network and ω¯¯𝜔\bar{\omega}over¯ start_ARG italic_ω end_ARG represents the parameters of the target work.

Since all Critic networks share common features, we establish a shared loss function to centrally optimize all the Critics (Eq. 23). This loss function (regarded as a measure of training effectiveness) is based on the sum of differences between expected and actual Q-values.

Q(ψ)=i=1NE(o,a,r,o)[(Qiψ(o,a)yi)2],subscript𝑄𝜓superscriptsubscript𝑖1𝑁subscript𝐸𝑜𝑎𝑟superscript𝑜delimited-[]superscriptsuperscriptsubscript𝑄𝑖𝜓𝑜𝑎subscript𝑦𝑖2\begin{array}[]{c}\mathcal{L}_{Q}(\psi)=\sum_{i=1}^{N}E_{\left(o,a,r,o^{\prime% }\right)}\left[\left(Q_{i}^{\psi}(o,a)-y_{i}\right)^{2}\right],\end{array}start_ARRAY start_ROW start_CELL caligraphic_L start_POSTSUBSCRIPT italic_Q end_POSTSUBSCRIPT ( italic_ψ ) = ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT italic_E start_POSTSUBSCRIPT ( italic_o , italic_a , italic_r , italic_o start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT [ ( italic_Q start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT ( italic_o , italic_a ) - italic_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , end_CELL end_ROW end_ARRAY (23)

where

yi=ri+γQiψ¯(o,a)subscript𝑦𝑖subscript𝑟𝑖𝛾superscriptsubscript𝑄𝑖¯𝜓superscriptosuperscript𝑎y_{i}=r_{i}+\gamma Q_{i}^{\bar{\psi}}\left(\mathrm{o}^{\prime},a^{\prime}\right)italic_y start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_γ italic_Q start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_ψ end_ARG end_POSTSUPERSCRIPT ( roman_o start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) (24)

To slightly adjust the policy of each agent, we utilize gradient ascent to update their policy parameters (Eq. 25).

θiJ(πθ)=Eo,a𝒟[θiπθi(aioi)aiQiψ(o,a)]subscriptsubscript𝜃𝑖𝐽subscript𝜋𝜃subscript𝐸similar-to𝑜𝑎𝒟delimited-[]subscriptsubscript𝜃𝑖subscript𝜋subscript𝜃𝑖conditionalsubscript𝑎𝑖subscript𝑜𝑖subscriptsubscript𝑎𝑖superscriptsubscript𝑄𝑖𝜓𝑜𝑎\nabla_{\theta_{i}}J\left(\pi_{\theta}\right)=E_{o,a\sim\mathcal{D}}\left[% \nabla_{\theta_{i}}\pi_{\theta_{i}}\left(a_{i}\mid o_{i}\right)\nabla_{a_{i}}Q% _{i}^{\psi}(o,a)\right]∇ start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_J ( italic_π start_POSTSUBSCRIPT italic_θ end_POSTSUBSCRIPT ) = italic_E start_POSTSUBSCRIPT italic_o , italic_a ∼ caligraphic_D end_POSTSUBSCRIPT [ ∇ start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_π start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∣ italic_o start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ∇ start_POSTSUBSCRIPT italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT ( italic_o , italic_a ) ] (25)

where ai=πθi(oi)subscript𝑎𝑖subscript𝜋subscript𝜃𝑖subscript𝑜𝑖a_{i}=\pi_{\theta_{i}}(o_{i})italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_π start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_o start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ).

V-C Proposed cooperative tracking policy

In this section, we introduce the proposed reward function. In reinforcement learning, by modeling the problem as MDP, the problem can be transformed into a corresponding maximization of the cumulative expected reward. Therefore, we define the reward function as follows (Eq. 26):

Ri(s,a)=αriδ+βriζ+δriμ+γriσsubscript𝑅𝑖𝑠𝑎𝛼superscriptsubscript𝑟𝑖𝛿𝛽superscriptsubscript𝑟𝑖𝜁𝛿superscriptsubscript𝑟𝑖𝜇𝛾superscriptsubscript𝑟𝑖𝜎R_{i}(s,a)=\alpha r_{i}^{\delta}+\beta r_{i}^{\zeta}+\delta r_{i}^{\mu}+\gamma r% _{i}^{\sigma}italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_s , italic_a ) = italic_α italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_δ end_POSTSUPERSCRIPT + italic_β italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ζ end_POSTSUPERSCRIPT + italic_δ italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT + italic_γ italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_σ end_POSTSUPERSCRIPT (26)

where riδsuperscriptsubscript𝑟𝑖𝛿r_{i}^{\delta}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_δ end_POSTSUPERSCRIPT is the tracking reward, riζsuperscriptsubscript𝑟𝑖𝜁r_{i}^{\zeta}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ζ end_POSTSUPERSCRIPT is the collision-avoidance reward, riμsuperscriptsubscript𝑟𝑖𝜇r_{i}^{\mu}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT is the ocean current stability reward, and riσsuperscriptsubscript𝑟𝑖𝜎r_{i}^{\sigma}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_σ end_POSTSUPERSCRIPT is the velocity matching reward.

For ET-AUV i𝑖iitalic_i, the tracking reward riδsuperscriptsubscript𝑟𝑖𝛿r_{i}^{\delta}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_δ end_POSTSUPERSCRIPT is designed to control the distance between EC-AUV and the tracking target, thereby maintaining the agent and the target within an best tracking distance range (Eq. 27).

riδ=ω1|d(i,j)ϕdbeϕ|dbeϕsuperscriptsubscript𝑟𝑖𝛿subscript𝜔1subscriptsuperscript𝑑italic-ϕ𝑖𝑗subscriptsuperscript𝑑italic-ϕbesubscriptsuperscript𝑑italic-ϕber_{i}^{\delta}=-\omega_{1}\frac{|d^{\phi}_{(i,j)}-d^{\phi}_{\text{be}}|}{d^{% \phi}_{\text{be}}}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_δ end_POSTSUPERSCRIPT = - italic_ω start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT divide start_ARG | italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ( italic_i , italic_j ) end_POSTSUBSCRIPT - italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT | end_ARG start_ARG italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT end_ARG (27)

where d(i,j)ϕsubscriptsuperscript𝑑italic-ϕ𝑖𝑗d^{\phi}_{(i,j)}italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ( italic_i , italic_j ) end_POSTSUBSCRIPT is the distance between ET-AUV i𝑖iitalic_i and target j𝑗jitalic_j, and dbeϕsubscriptsuperscript𝑑italic-ϕbed^{\phi}_{\text{be}}italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT is the best tracking distance.

At the same time, to ensure a safe distance between ET-AUVs during the tracking process, we introduce a collision-avoidance reward riζsuperscriptsubscript𝑟𝑖𝜁r_{i}^{\zeta}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ζ end_POSTSUPERSCRIPT (Eq. 28):

riζ=ω2j=1jiN|d(i,j)ψdbeψ|dbeψsuperscriptsubscript𝑟𝑖𝜁subscript𝜔2superscriptsubscript𝑗1𝑗𝑖𝑁subscriptsuperscript𝑑𝜓𝑖𝑗subscriptsuperscript𝑑𝜓besubscriptsuperscript𝑑𝜓ber_{i}^{\zeta}=-\omega_{2}\sum_{\begin{subarray}{c}j=1\\ j\neq i\end{subarray}}^{N}\frac{|d^{\psi}_{(i,j)}-d^{\psi}_{\text{be}}|}{d^{% \psi}_{\text{be}}}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ζ end_POSTSUPERSCRIPT = - italic_ω start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT start_ARG start_ROW start_CELL italic_j = 1 end_CELL end_ROW start_ROW start_CELL italic_j ≠ italic_i end_CELL end_ROW end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT divide start_ARG | italic_d start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ( italic_i , italic_j ) end_POSTSUBSCRIPT - italic_d start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT | end_ARG start_ARG italic_d start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT end_ARG (28)

where d(i,j)ψsubscriptsuperscript𝑑𝜓𝑖𝑗d^{\psi}_{(i,j)}italic_d start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ( italic_i , italic_j ) end_POSTSUBSCRIPT is the distance between ET-AUV i𝑖iitalic_i and j𝑗jitalic_j, and dbeψsubscriptsuperscript𝑑𝜓bed^{\psi}_{\text{be}}italic_d start_POSTSUPERSCRIPT italic_ψ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT be end_POSTSUBSCRIPT is the best distance between ET-AUVs.

To promote the smooth movement of ET-AUVs under the influence of ocean currents, we introduce an ocean current stability reward riμsuperscriptsubscript𝑟𝑖𝜇r_{i}^{\mu}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT (Eq. 29). This reward helps mitigate the effects on ET-AUV’s velocity caused by ocean currents.

riμ=ω3vcur(i)vpre(i)vpre(i)superscriptsubscript𝑟𝑖𝜇subscript𝜔3delimited-∥∥subscript𝑣cur𝑖subscript𝑣pre𝑖delimited-∥∥subscript𝑣pre𝑖r_{i}^{\mu}=-\omega_{3}\frac{\lVert v_{\text{cur}}(i)-v_{\text{pre}}(i)\rVert}% {\lVert v_{\text{pre}}(i)\rVert}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = - italic_ω start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT divide start_ARG ∥ italic_v start_POSTSUBSCRIPT cur end_POSTSUBSCRIPT ( italic_i ) - italic_v start_POSTSUBSCRIPT pre end_POSTSUBSCRIPT ( italic_i ) ∥ end_ARG start_ARG ∥ italic_v start_POSTSUBSCRIPT pre end_POSTSUBSCRIPT ( italic_i ) ∥ end_ARG (29)

where vcur(i)subscript𝑣cur𝑖v_{\text{cur}}(i)italic_v start_POSTSUBSCRIPT cur end_POSTSUBSCRIPT ( italic_i ) is ET-AUV i𝑖iitalic_i’s current velocity, vpre(i)subscript𝑣pre𝑖v_{\text{pre}}(i)italic_v start_POSTSUBSCRIPT pre end_POSTSUBSCRIPT ( italic_i ) is ET-AUV i𝑖iitalic_i’s previous velocity.

To avoid excessive energy consumption by ET-AUVs due to frequent acceleration and deceleration, we introduce a velocity matching reward riσsuperscriptsubscript𝑟𝑖𝜎r_{i}^{\sigma}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_σ end_POSTSUPERSCRIPT (Eq. 30):

riσ=ω4vage(i)vtar(j)vtar(j)superscriptsubscript𝑟𝑖𝜎subscript𝜔4delimited-∥∥subscript𝑣age𝑖subscript𝑣tar𝑗delimited-∥∥subscript𝑣tar𝑗r_{i}^{\sigma}=-\omega_{4}\frac{\lVert v_{\text{age}}(i)-v_{\text{tar}}(j)% \rVert}{\lVert v_{\text{tar}}(j)\rVert}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_σ end_POSTSUPERSCRIPT = - italic_ω start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT divide start_ARG ∥ italic_v start_POSTSUBSCRIPT age end_POSTSUBSCRIPT ( italic_i ) - italic_v start_POSTSUBSCRIPT tar end_POSTSUBSCRIPT ( italic_j ) ∥ end_ARG start_ARG ∥ italic_v start_POSTSUBSCRIPT tar end_POSTSUBSCRIPT ( italic_j ) ∥ end_ARG (30)

where vage(i)subscript𝑣age𝑖v_{\text{age}}(i)italic_v start_POSTSUBSCRIPT age end_POSTSUBSCRIPT ( italic_i ) is the velocity of the agent, vtar(j)subscript𝑣tar𝑗v_{\text{tar}}(j)italic_v start_POSTSUBSCRIPT tar end_POSTSUBSCRIPT ( italic_j ) is the velocity of the target.

Based on the AUV swarm network, the cooperative underwater target tracking algorithm can be summarized as shown in Algorithm 3.

The proposed Algorithm 3 combines the proposed DSBM architecture with the proposed tracking policy. This integration enables accurate and cooperative performance in underwater cooperative target tracking, while the ET-AUV’s endurance and resistance to ocean current interference is taken into account.

Algorithm 3 FATHOM-Net-based DSBM tracking algorithm
1:Initialize the Set NEsubscript𝑁𝐸N_{E}italic_N start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT, T𝑇Titalic_T, N𝑁Nitalic_N, respectively.
2:Initialize the weights of Actor and Critic networks corresponding to each ET-AUV.
3:Initialize a random process 𝒩𝒩\mathcal{N}caligraphic_N for exploration of action.
4:Initialize replay buffer, environment.
5:for Episode=1Episode1\text{Episode}=1Episode = 1 to NEsubscript𝑁𝐸N_{E}italic_N start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT do
6:     Reset environment.
7:     Assign tracking formations according to Algorithm 2.
8:     for t=1𝑡1t=1italic_t = 1 to T𝑇Titalic_T do
9:         Each ET-AUV takes action ai=πθi(oi)+𝒩tsubscript𝑎𝑖subscript𝜋subscript𝜃𝑖subscript𝑜𝑖subscript𝒩𝑡a_{i}=\pi_{\theta_{i}}(o_{i})+\mathcal{N}_{t}italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_π start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_o start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) + caligraphic_N start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT based on the policy and exploration.
10:         Calculate current force according to Algorithm 1.
11:         Update environment, obtain next state and reward.
12:         Add experience to the replay buffer.
13:         for i=1𝑖1i=1italic_i = 1 to N𝑁Nitalic_N do
14:              Select ET-AUVs based on dynamic-switching attention mechanism.
15:              Calculate Qi(ψ)subscript𝑄𝑖𝜓Q_{i}(\psi)italic_Q start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ψ ) based on selected ET-AUVs.
16:              Update the actor according to Eq. 25.
17:              Update the Critic according to Eq. 23.               
18:     Perform dynamic-switching resampling mechanism from (xt,at,rt,xt+1)subscript𝑥𝑡subscript𝑎𝑡subscript𝑟𝑡subscript𝑥𝑡1(x_{t},a_{t},r_{t},x_{t+1})( italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_a start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_r start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_x start_POSTSUBSCRIPT italic_t + 1 end_POSTSUBSCRIPT ).
19:     Update the target network of each ET-AUV according to Eq. 22.

VI Evaluations

This section presents the evaluation results, where we conduct an in-depth analysis of the newly proposed tracking algorithm from two perspectives: the effectiveness of the algorithm and the performance of target tracking. We demonstrate our proposed algorithm by comparing it with some of current mainstream MARL algorithms.

VI-A Simulation setup

The test is conducted on a computer equipped with an Intel(R) Core(TM) i9-12900H 2.50GHz processor and 32GB of RAM.

We construct a model using coordinate gridization, treating ET-AUVs as moving particles in a 3D underwater environment. The targets initially locate near the origin and move with a random velocity, while the ET-AUVs are uniformly distributed within a region one kilometer away from the targets.

The relevant evaluation notations is detailed in Table II.

TABLE II: Simulation notation
Notation Description Value
NAsubscript𝑁𝐴N_{A}italic_N start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT Number of ET-AUVs [4,6,12]
NTsubscript𝑁𝑇N_{T}italic_N start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT Number of Targets [2,3,4]
lrsubscript𝑙𝑟l_{r}italic_l start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT Learning rate 3e33𝑒33e-33 italic_e - 3
NEsubscript𝑁𝐸N_{E}italic_N start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT Number of training rounds 5000
Nhsubscript𝑁N_{h}italic_N start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT Number of neurons in the hidden layer 64
γ𝛾\gammaitalic_γ Discount factor 0.95
τ𝜏\tauitalic_τ Network update coefficient 1e21𝑒21e-21 italic_e - 2
dminϕsubscriptsuperscript𝑑italic-ϕmind^{\phi}_{\text{min}}italic_d start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT min end_POSTSUBSCRIPT Minimum tracking distance 80 m
dminκsubscriptsuperscript𝑑𝜅mind^{\kappa}_{\text{min}}italic_d start_POSTSUPERSCRIPT italic_κ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT min end_POSTSUBSCRIPT Minimum distance among ET-AUVs 80 m
Refer to caption
Figure 4: Convergence speed comparison under stable ocean currents (4 ET-AUVs tracking 2 targets)
Refer to caption
Figure 5: Convergence speed comparison under ocean current interference (4 ET-AUVs tracking 2 targets)
Refer to caption
Figure 6: Convergence speed comparison under stable ocean currents (12 ET-AUVs tracking 4 targets)
Refer to caption
Figure 7: Convergence speed comparison under ocean current interference (12 ET-AUVs tracking 4 targets)

VI-B Results

First, to measure the effectiveness of the proposed DSBM, we compare DSBM, with MADDPG [36], MAAC [37], and MADDPG-SAC [38]. Note that MADDPG-SAC is based on MADDPG where the entropy value (that is adaptively adjusted by SAC), is incorporated into the optimization goal to change the update of the Critic network. We compare DSBM with MADDPG-SAC to evaluate the significance of entropy value in the algorithm comparison. We conduct the comparison by respectively utilizing different algorithms to tracking the underwater targets in terms of convergence speed, tracking distance maintenance, difference in velocity from target, and intra-formation strategy consistency.

Convergence speed: To evaluate the stability of algorithm convergence, we repeat multiple rounds of experiments (under different tracking scenarios) and select the samples in a 95% confidence interval. Totally, as shown in Fig. 4 to Fig. 7, our proposed DBSM performs the best than the other three compared objects in convergence speed.

Refer to caption
Figure 8: Distribution of Agent-Target distance with stable ocean currents
Refer to caption
Figure 9: Distribution of Agent-Target distance with ocean current interference
Refer to caption
Figure 10: Agent-Target velocity difference with stable ocean currents
Refer to caption
Figure 11: Agent-Target velocity difference with ocean current interference

In details, in Fig. 4, the proposed DSBM showcases the best convergence performance, even if the other algorithms eventually show a trend of convergence. It’s worth noting, MADDPG-SAC converges the fastest in the early stages of training, but the final convergence effect is not satisfactory. This is because the entropy value-based mechanism in MADDPG-SAC improves the exploration scope in the early stage of training, but it also results in that the training to be unstable in the later stages. Different with the case, as shown in Fig. 5, where the ocean current interference is employed, only the DSBM algorithm is able to converge. This demonstrates that our proposed DSBM’s superior adaptability to ocean currents.

Refer to caption
Figure 12: Strategic consistency with stable ocean currents
Refer to caption
Figure 13: Strategic consistency with ocean current interference

Furthermore, as shown in Fig. 6 and Fig. 7, DSBM also showcases the fastest convergence speed and the best performance with/without ocean current interference.

The results in Fig. 4, Fig. 5, Fig. 6 and Fig. 7 demonstrate that: DSBM pays more attention to the information of better-performing ET-AUVs and conducts more meaningful exploration, thus the time the ineffective exploration can be significantly reduced.

Tracking distance maintenance: As it is aforementioned, maintaining the distance between the agent and the target, i.e., the ”Agent-Target distance” during the tracking process is crucial. Thus, we test as well the distribution of Agent-Target distance with/without ocean current interference.

As the results shown in Fig. 8, with stable ocean current, DSBM can maintain the distance between ET-AUV and its targets within a small (pre-defined) distance change range. In contrast, the other algorithms all display significant deviations. Further, in Fig. 9, the similar results are also achieved in the scenario with ocean current interference.

Velocity difference maintenance: In AUV swarm network-based underwater cooperative tracking, frequent acceleration and deceleration for ET-AUV can rapidly consume the energy, making against the long-duration tracking tasks especially in underwater. Therefore, the velocity difference between the ET-AUV and target, i.e., the ”Agent-Target velocity difference” can be used to measure the endurance capacity and tracking efficiency of the tracking policy.

As shown in Fig. 10, with stable ocean current conditions, the Agent-Target velocity differences of DSBM are within a very small range. In contrast, the other three algorithms showcase a wider range of Agent-Target velocity differences, indicating that the ET-AUVs have to frequently accelerate and decelerate throughout the entire tracking process. This will significantly reduce the endurance time under the same state conditions.

As shown in Fig. 11, with ocean current interference, the velocity of ET-AUVs guided by DSBM almost maintains consistency with the velocity of the tracking target. The results in both Fig. 10 and Fig. 11 demonstrate that the proposed DSBM can perform stable and high-efficient underwater cooperative tracking with complex underwater environment (with/without ocean current interference).

Refer to caption
(a) Scenario 1_Phase 1
Refer to caption
(b) Scenario 1_Phase 2
Refer to caption
(c) Scenario 1_Phase 3
Refer to caption
(d) Scenario 2_Phase 1
Refer to caption
(e) Scenario 2_Phase 2
Refer to caption
(f) Scenario 2_Phase 3
Refer to caption
(g) Scenario 3_Phase 1
Refer to caption
(h) Scenario 3_Phase 2
Refer to caption
(i) Scenario 3_Phase 3
Figure 14: Test for the availability of the proposed approaches in 3D environment.

Intra-formation strategy consistency: During the tracking process, a consistent strategy for the AUV formation helps improve efficiency in target identification and tracking, especially in complex underwater environments. This is crucial to prevent tracking errors or relief target loss due to environmental factors. Hence, we take the experimental scenario where twelve ET-AUVs tracking four targets as an example. And the ET-AUVs in each formation each choose actions from seven possible actions, to evaluate the consistency of these actions.

As the results shown in both Fig. 12 and 13, we showcase the probability distribution of the number of the consistent ET-AUVs in the formation. As shown in Fig. 12, before the ocean current interference is taken into account, the actions of all the ET-AUVs under our proposed DSBM are more likely consistent. Similarly, the proposed DSBM obviously performs the best than the other compared objects, while the ocean current interference is considered. This benefits from the proposed DSBM is based on our proposed FATHOM-Net framework, by which the AUV formation within a region undergoes local training, focusing more on the tracking targets within the region rather than the entire tracking task.

Finally, to make the entire tracking process directed by our proposed DSBM more clear, we use Unity to simulate the underwater environment for visualizing the tracking process. As shown in Fig. 14, the cubes represent the targets, and the spheres represent ET-AUVs. The red solid line represents the movement trajectory of the target, while other colored solid lines represent the tracking trajectories of the ET-AUVs. Specifically, we simulate three scenarios of multi-target tracking. Fig. 14(a)-Fig. 14(c) represent four ET-AUVs tracking two targets, Fig. 14(d)-Fig. 14(f) represent six ET-AUVs tracking three targets, and Fig. 14(g)-Fig. 14(i) represent twelve ET-AUVs tracking four targets. The visualization in Fig. 14 also proves the effectiveness of the proposed DSBM in tracking underwater targets in a 3D environment.

VII Conclusion

In this study, we propose a dynamic-switching-enabled multi-agent reinforcement learning framework for performing AUV swarm-network-based cooperative multi-target tracking in underwater environment. We first give modelling approach for the underwater sonar-based detection and the ocean current interference. Then, we propose FATHOM-Net, a novel MARL architecture based on SDN, which enhances the flexibility and robustness of the AUV swarm. The proposed DSBM includes the ”dynamic-switching attention” and ”dynamic-switching resampling” mechanisms, both based on an experience replay buffer. It enhances AUV swarm network’s self-learning efficiency and accuracy. Finally, to achieve efficient formation classification in the interference of ocean currents, we propose an algorithm named ASMA. It efficiently and accurately assigns formations for ET-AUV through fuzzy logic and rule-based expert systems. Evaluation results demonstrate that our proposed tracking algorithm achieves the fastest convergence speed and tracking performance compared to various popular research products. Future research directions derived from this work can be summarized as the following:

  1. 1.

    Optimizing underwater obstacle avoidance mechanisms for ET-AUVs to mitigate potential damage;

  2. 2.

    Balancing energy consumption among the AUVs to enhance the endurance of the AUV swarm network-based cooperative tracking systems;

  3. 3.

    Designing AUV swarm robustness control framework when the unstable underwater communication (e.g., the underwater acoustic-based communication) has to been taken into account.

References

  • [1] F. A. Setiawan and P. Rahmadi, “Indoalgae: The database of indonesian native strains of potential marine algae,” in 2021 Sixth International Conference on Informatics and Computing (ICIC), 2021, pp. 1–5.
  • [2] Z. Yan, “Revolution system of spatial-temporal pattern of marine resources development and management based on deep learning algorithms,” in 2023 International Conference on Data Science and Network Security (ICDSNS), 2023, pp. 1–6.
  • [3] A. Rumson, A. B. M. Tvedt, and C. Aas, “An efficient coring technology for deep sea mineral exploration,” in OCEANS 2023 - Limerick, 2023, pp. 1–11.
  • [4] L. Zong, H. Wang, and G. Luo, “Transmission control over satellite network for marine environmental monitoring system,” IEEE Transactions on Intelligent Transportation Systems, vol. 23, no. 10, pp. 19 668–19 675, 2022.
  • [5] S. Yuan, Y. Li, F. Bao, H. Xu, Y. Yang, Q. Yan, S. Zhong, H. Yin, J. Xu, Z. Huang, and J. Lin, “Marine environmental monitoring with unmanned vehicle platforms: Present applications and future prospects,” Science of The Total Environment, vol. 858, no. Part 1, p. 159741, 2023.
  • [6] S. Hou, W. Li, T. Liu, S. Zhou, J. Guan, R. Qin, and Z. Wang, “D2cl: A dense dilated convolutional lstm model for sea surface temperature prediction,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 14, pp. 12 514–12 523, 2021.
  • [7] Q. Ma, D. Zhang, C. Wan, J. Zhang, and N. Lyu, “Multi-objective emergency resources allocation optimization for maritime search and rescue considering accident black-spots,” Ocean Engineering, vol. 261, p. 112178, 2022.
  • [8] K. Shojaei, “Neural network feedback linearization target tracking control of underactuated autonomous underwater vehicles with a guaranteed performance,” Ocean Engineering, vol. 258, p. 111827, 2022.
  • [9] X. Wang, B. Xu, and Y. Guo, “Fuzzy logic system-based robust adaptive control of auv with target tracking,” International Journal of Fuzzy Systems, vol. 25, no. 1, pp. 338–346, 2023.
  • [10] G. Han, A. Gong, H. Wang, M. Martínez-García, and Y. Peng, “Multi-auv collaborative data collection algorithm based on q-learning in underwater acoustic sensor networks,” IEEE Transactions on Vehicular Technology, vol. 70, no. 9, pp. 9294–9305, 2021.
  • [11] Z. Yan, K. Zhang, L. Qiao, Y. Hu, and B. Song, “A multiload wireless power transfer system with concentrated magnetic field for auv cluster system,” IEEE Transactions on Industry Applications, vol. 58, no. 1, pp. 1307–1314, 2022.
  • [12] B. Jiang, J. Du, C. Jiang, Z. Han, and M. Debbah, “Underwater searching and multiround data collection via auv swarms: An energy-efficient aoi-aware mappo approach,” IEEE Internet of Things Journal, vol. 11, no. 7, pp. 12 768–12 782, 2024.
  • [13] N. T. Hung, F. F. C. Rego, and A. M. Pascoal, “Cooperative distributed estimation and control of multiple autonomous vehicles for range-based underwater target localization and pursuit,” IEEE Transactions on Control Systems Technology, vol. 30, no. 4, pp. 1433–1447, 2022.
  • [14] T. Zhou, Y. Wang, L. Zhang, B. Chen, and X. Yu, “Underwater multitarget tracking method based on threshold segmentation,” IEEE Journal of Oceanic Engineering, vol. 48, no. 4, pp. 1255–1269, 2023.
  • [15] H. Flores, N. H. Motlagh, A. Zuniga, M. Liyanage, M. Passananti, S. Tarkoma, M. Youssef, and P. Nurmi, “Toward large-scale autonomous marine pollution monitoring,” IEEE Internet of Things Magazine, vol. 4, no. 1, pp. 40–45, 2021.
  • [16] Y. Girdhar, N. McGuire, L. Cai, S. Jamieson, S. McCammon, B. Claus, J. E. S. Soucie, J. E. Todd, and T. A. Mooney, “Curee: A curious underwater robot for ecosystem exploration,” in 2023 IEEE International Conference on Robotics and Automation (ICRA), 2023, pp. 11 411–11 417.
  • [17] S. Rani, H. Babbar, P. Kaur, M. D. Alshehri, and S. H. Shah, “An optimized approach of dynamic target nodes in wireless sensor network using bio inspired algorithms for maritime rescue,” IEEE Transactions on Intelligent Transportation Systems, vol. 24, no. 2, pp. 2548–2555, 2023.
  • [18] G. Soldi, D. Gaglione, N. Forti, A. D. Simone, F. C. Daffinà, G. Bottini, D. Quattrociocchi, L. M. Millefiori, P. Braca, S. Carniel, P. Willett, A. Iodice, D. Riccio, and A. Farina, “Space-based global maritime surveillance. part i: Satellite technologies,” IEEE Aerospace and Electronic Systems Magazine, vol. 36, no. 9, pp. 8–28, 2021.
  • [19] Z. Chu, F. Wang, T. Lei, and C. Luo, “Path planning based on deep reinforcement learning for autonomous underwater vehicles under ocean current disturbance,” IEEE Transactions on Intelligent Vehicles, vol. 8, no. 1, pp. 108–120, 2023.
  • [20] W. Lan, X. Jin, T. Wang, and H. Zhou, “Improved rrt algorithms to solve path planning of multi-glider in time-varying ocean currents,” IEEE Access, vol. 9, pp. 158 098–158 115, 2021.
  • [21] Y. Yang, Y. Xiao, and T. Li, “A survey of autonomous underwater vehicle formation: Performance, formation control, and communication capability,” IEEE Communications Surveys & Tutorials, vol. 23, no. 2, pp. 815–841, 2021.
  • [22] Z. Fang, J. Wang, J. Du, X. Hou, Y. Ren, and Z. Han, “Stochastic optimization-aided energy-efficient information collection in internet of underwater things networks,” IEEE Internet of Things Journal, vol. 9, no. 3, pp. 1775–1789, 2022.
  • [23] H. Guo, Z. Sun, and P. Wang, “Joint design of communication, wireless energy transfer, and control for swarm autonomous underwater vehicles,” IEEE Transactions on Vehicular Technology, vol. 70, no. 2, pp. 1821–1835, 2021.
  • [24] J. Tang, G. Liu, and Q. Pan, “A review on representative swarm intelligence algorithms for solving optimization problems: Applications and trends,” IEEE/CAA Journal of Automatica Sinica, vol. 8, no. 10, pp. 1627–1643, 2021.
  • [25] C. Lin, G. Han, J. Jiang, C. Li, S. B. H. Shah, and Q. Liu, “Underwater pollution tracking based on software-defined multi-tier edge computing in 6g-based underwater wireless networks,” IEEE Journal on Selected Areas in Communications, vol. 41, no. 2, pp. 491–503, 2023.
  • [26] C. Lin, G. Han, T. Wang, Y. Bi, J. Du, and B. Zhang, “Fast node clustering based on an improved birch algorithm for data collection towards software-defined underwater acoustic sensor networks,” IEEE Sensors Journal, vol. 21, no. 22, pp. 25 480–25 488, 2021.
  • [27] C. Lin, G. Han, T. Zhang, S. B. H. Shah, and Y. Peng, “Smart underwater pollution detection based on graph-based multi-agent reinforcement learning towards auv-based network its,” IEEE Transactions on Intelligent Transportation Systems, vol. 24, no. 7, pp. 7494–7505, 2023.
  • [28] F. Li, H. Yao, J. Du, C. Jiang, Z. Han, and Y. Liu, “Auction design for edge computation offloading in sdn-based ultra dense networks,” IEEE Transactions on Mobile Computing, vol. 21, no. 5, pp. 1580–1595, 2022.
  • [29] D. M. Casas-Velasco, O. M. C. Rendon, and N. L. S. da Fonseca, “Drsir: A deep reinforcement learning approach for routing in software-defined networking,” IEEE Transactions on Network and Service Management, vol. 19, no. 4, pp. 4807–4820, 2022.
  • [30] C. Lin, G. Han, M. Guizani, Y. Bi, J. Du, and L. Shu, “An sdn architecture for auv-based underwater wireless networks to enable cooperative underwater search,” IEEE Wireless Communications, vol. 27, no. 3, pp. 132–139, 2020.
  • [31] D. Q. Huy, N. Sadjoli, A. B. Azam, B. Elhadidi, Y. Cai, and G. Seet, “Object perception in underwater environments: a survey on sensors and sensing methodologies,” Ocean Engineering, vol. 267, p. 113202, 2023.
  • [32] G. Chen, Y. Shen, N. Qu, and B. He, “Path planning of auv during diving process based on behavioral decision-making,” Ocean Engineering, vol. 234, p. 109073, 2021.
  • [33] Y. Fang, Z. Huang, J. Pu, and J. Zhang, “Auv position tracking and trajectory control based on fast-deployed deep reinforcement learning method,” Ocean Engineering, vol. 245, p. 110452, 2022.
  • [34] Z. Yang, J. Du, Z. Xia, C. Jiang, A. Benslimane, and Y. Ren, “Secure and cooperative target tracking via auv swarm: A reinforcement learning approach,” in 2021 IEEE Global Communications Conference (GLOBECOM), 2021, pp. 1–6.
  • [35] Z. Yan, C. Zhang, W. Tian, and M. Zhang, “Formation trajectory tracking control of discrete-time multi-auv in a weak communication environment,” Ocean Engineering, vol. 245, p. 110495, 2022.
  • [36] R. Lowe, Y. Wu, A. Tamar, J. Harb, P. Abbeel, and I. Mordatch, “Multi-agent actor-critic for mixed cooperative-competitive environments,” Neural Information Processing Systems (NIPS), 2017.
  • [37] S. Iqbal and F. Sha, “Actor-attention-critic for multi-agent reinforcement learning,” in International conference on machine learning, vol. 97.   PMLR, 2019, pp. 2961–2970.
  • [38] H. Yong, J. Seo, J. Kim, M. Kim, and J. Choi, “Suspension control strategies using switched soft actor-critic models for real roads,” IEEE Transactions on Industrial Electronics, vol. 70, no. 1, pp. 824–832, Jan 2023.
[Uncaptioned image] Shengbo Wang is currently pursuing a Bachelor’s degree at the Software College, Northeastern University, Shenyang, China. His current research interests include reinforcement learning, computer vision, and machine learning.
[Uncaptioned image] Chuan Lin [S’17, M’20] is currently an associate professor with the Software College, Northeastern University, Shenyang, China. He received the B.S. degree in Computer Science and Technology from Liaoning University, Shenyang, China in 2011, the M.S. degree in Computer Science and Technology from Northeastern University, Shenyang, China in 2013, and the Ph.D. degree in computer architecture in 2018. From Nove. 2018 to Nove. 2020, he is a Postdoctoral Researcher with the School of Software, Dalian University of Technology, Dalian, China. His research interests include UWSNs, industrial IoT, software-defined networking.
[Uncaptioned image] Guangjie Han [S’03-M’05-SM’18-F’22] is currently a Professor with the Department of Internet of Things Engineering, Hohai University, Changzhou, China. He received his Ph.D. degree from Northeastern University, Shenyang, China, in 2004. In February 2008, he finished his work as a Postdoctoral Researcher with the Department of Computer Science, Chonnam National University, Gwangju, Korea. From October 2010 to October 2011, he was a Visiting Research Scholar with Osaka University, Suita, Japan. From January 2017 to February 2017, he was a Visiting Professor with City University of Hong Kong, China. From July 2017 to July 2020, he was a Distinguished Professor with Dalian University of Technology, China. His current research interests include Internet of Things, Industrial Internet, Machine Learning and Artificial Intelligence, Mobile Computing, Security and Privacy. Dr. Han has over 500 peer-reviewed journal and conference papers, in addition to 160 granted and pending patents. Currently, his H-index is 65 and i10-index is 282 in Google Citation (Google Scholar). The total citation count of his papers raises above 15500+ times. Dr. Han is a Fellow of the UK Institution of Engineering and Technology (FIET). He has served on the Editorial Boards of up to 10 international journals, including the IEEE TII, IEEE TCCN, IEEE TVT, IEEE Systems, etc. He has guest-edited several special issues in IEEE Journals and Magazines, including the IEEE JSAC, IEEE Communications, IEEE Wireless Communications, Computer Networks, etc. Dr. Han has also served as chair of organizing and technical committees in many international conferences. He has been awarded 2020 IEEE Systems Journal Annual Best Paper Award and the 2017-2019 IEEE ACCESS Outstanding Associate Editor Award. He is a Fellow of IEEE.
[Uncaptioned image] Shengchao Zhu received his B.S. degree in Internet of Things Engineering from Hohai University,Changzhou, China, in 2023. He is currently pursuing the Ph.D. degree with the Department of Computer Science and Technology at Hohai University, Nan- jing, China. His current research interests include swarm intelligence, swarm ocean, Multi-Agent Reinforcement Learning.
[Uncaptioned image] Zhixian Li is currently pursuing a bachelor’s degree at the Software College, Northeastern University, Shenyang, China. His research interests include natural language processing, computer vision and distributed computing.
[Uncaptioned image] Zhenyu Wang is currently pursuing a Bachelor’s degree at the Software College, Northeastern University, Shenyang, China. His research interests include machine learning, software architecture and data mining.