TY - JOUR AU - Li, Han AB - 1 Introduction With the development of modern intelligent technology, unmanned surface vehicles (USV) are increasingly in demand for civilian applications, including environmental protection, water-quality testing and aquaculture, effectively bridging labor shortages and improving production efficiency [1–3]. However, USV have posed some inherent challenges like poor autonomous navigation performance, large trajectory tracking error, poor adaptive ability, and single application scenario that could impact the future of USV market. First, there is too much manual intervention in the application of USV. To reduce labor costs and improve efficiency, it is necessary to combine autonomous navigation and control algorithms to design trajectory tracking. Besides, prescribed route tracking within USV also needs to be carried out as route error may result in yaw and instability control [4, 5]. Currently, an efficient design algorithm is essential for effective navigation performance of USV, and it has been becoming a hot topic of scientific studies in society now. Second, the prompt development of the integration within control system is crucial to the future development of USV. When the control integration is high, the advantages are more obvious, and the adaptive ability can be improved significantly [6]. Lastly, increasing the demand for multi-scene applications of USV is necessary to track more path points, as it leads to more precise of the driving lines and thereby improves the adaptability of USV in different environment. Hardware structure has been investigated by many researchers, which includes Arduino, PX4, and ARM microcontroller. PX4 is widely adopted because of powerful functionality and open access. Zhao et al. [7] chose Arduino Mega2560 as the controller and built a control system with GPS, a data acquisition module and a motion control module driving four DC motors, which solved the dynamic real-time tracking problem under interference uncertainty and proved the feasibility of this control system. Xiong et al. [8] built Arduino Mega2560, GPS, GPRS, sensors of water quality detection, and motor drive hardware platform to link with a host computer, solving the problem of fixed-point detection for water quality. The literatures [9, 10] selected PX4 as the main controller for the regulation of USV to prove the feasibility of the controller. Following this, we adopted PX4 as the controller, the actuators employed GPS, data transmission module, and drive module as trajectory tracking preparation, including automatic navigation and motion. The trajectory tracking method has been commonly used to determine the effective control and stable operation. For the appropriate method of the trajectory tracking, it is a prerequisite that the control methodology of the trajectory tracking, either with deterministic model parameters or with uncertain model parameters, are known as precisely as possible. Moreover, the trajectory tracking methods of the control system usually include five types (line-of-sight (LOS) guidance algorithm [11], adaptive proportional-integral-differential (PID) algorithm [12], fuzzy control method [13], sliding mode control method [14], neural network control method [15]) during the regulating of control parameters. These methods could significantly improve control performance. Among them, adaptive PID algorithm, is considered with large amount of calculation, which leads to a cumulative error. Next, although fuzzy control method could correct the error, the accuracy will be reduced when fuzzy processing. So, we apply the control algorithm of adaptive fuzzy PID algorithm in this study, which could lower the error and improve the accuracy at the same time. Furthermore, the sliding mode control method and the neural network control method hardly meet the requirement of practical applications due to computational complexity and inability to accurately characterize the mathematical model. Therefore, to date, LOS is still the most effective and independent of the object under control. Here, we adopted the LOS guidance algorithm because of its advantages such as easy calculation, simple model, and easy trajectory tracking instructions. There have been numerous studies in the last 3–5 years on the LOS proposing diverse solutions of route processing. The literatures [16–18] almost all used the LOS guidance algorithms to address the problem of unstable route travel. Concurrently, many efforts have been paid to design the tracking schemes aiming to calculate position error in real-time. It is no doubt that the precision of heading control determines the effectiveness of tracking. Shen et al. [19, 20] designed an adaptive path following control method for the problem of straight path following for a class of underactuated bio-inspired snake robots on ground facing unknown friction coefficients, and verified the effectiveness of the controller. Cao et al. [21] introduced a homogeneous nonlinear extended state observer and a dynamic surface control in their study, the purpose was to adjust the saturation of the controller. They also compared their findings with other studies. It is noted that the controller could accurately track trajectories in different environment. Setiawan et al. [22] proposed a trajectory planning as the ground station by using PID method. Qiu et al. [23] established a neural network to reduce the computation and employed a hyperbolic tangent function to diminish sliding mode chattering based on an adaptive sliding mode control scheme. Iqbal et al. [24] proposed a fuzzy-immune adaptive system that enhanced the anti-jamming capability of a self-balancing robot system and proves the stability of the system. Moreover, Iqbal et al. [25] proposed a novel adaptive PD-type iterative learning control with robustness and effectiveness in handling nonlinearities. Considering the application characteristics of the fuzzy adaptive PID for heading control, we calculated the heading deviation and continuously adjusted the parameters. This work aims to develop an accurately-efficient trajectory tracking control system that could be used for on-site baiting, drug spreading and water quality testing in ponds. Additionally, we take a paddle boat, one of USV, as the research object. The control system does not implement a complex operation process but could capture main trajectory characteristics of the paddle boat. To this goal, we designed the LOS guidance and fuzzy adaptive PID algorithm, and achieved autonomous navigation and transmit relevant data such as heading, position and data transmission speed through constructing the hardware and software system. As a result, these data could be monitored in real-time by the host computer. After that, we carried out several tests using a designed operating platform under given routes in the pond to verify the feasibility of the trajectory tracking control. Overall, our trajectory tracking control system may pave the way to realize unprecedented smart aquaculture. This paper contains 4 sections: in Section 2, the focus is on the methodology. In Section 2.1, a three-degree-of-freedom kinematic model of the paddle boat is built; in Section 2.2, the hardware control device of the trajectory tracking system is constructed; and in Section 2.3, the algorithms and control methods used are highlighted. The experimental results and analysis results are given in Section 3. Finally, some conclusions are presented in Section 4. 2 Methodology 2.1 Model description To establish the kinematics model of the paddle boat, dimensions of the geometric structure used in this study are shown in Fig 1. The basic overall dimensions l×b×h equals 1770mm×1600mm×580mm, the water intrusion depths of the boat and two wheels are d = 110mm and c = 85mm, respectively. Then, the diameter (D) and width (s) of the wheels are 440mm and 110mm, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Geometric dimensions of the paddle boat: (a) rear view; (b)side view. https://doi.org/10.1371/journal.pone.0290246.g001 It is very important to know the acting force between the paddle boat and water. As a general rule, azimuth, pitch and roll angles have little effect on the motion in the horizontal plane. Therefore, we considered only a simplified kinematic model of that is also proposed by Fossen [26] and was characterized by the ease of observing the acting force, which facilitated analytical calculations. The difference between this model and this paper is that the rudder angle is transformed into the steering angle caused by the speed difference of the paddles. Considering the existing longitudinal swing, transverse swing and bow rocking, based on the inertial coordinate system On-XnYnZn fixed on the Earth’s surface and hull coordinate system O-X0Y0Z0 fixed on the paddle boat, we set up a three-degree-of-freedom kinematics model, as shown in Fig 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Schematic diagram of inertial coordinate system and hull coordinate system. https://doi.org/10.1371/journal.pone.0290246.g002 The transformation between the inertial coordinate system and the hull coordinate system can be expressed as: (1) where, x, y, and Ψ represent the x-point, y-point, and the bow angle of the paddle boat in the inertial coordinate system, respectively; u, v, and r represent the forward velocity, traverse velocity, and the swing angle velocity in the hull coordinate system, respectively. In solving the simulation model of the paddle boat, we processed the velocities of the left and right wheels as input parameters and the bow angle as an output parameter. Afterward, we analyzed the propulsive forces of the paddle wheels. According to Newton’s second law of motion, the propulsive forces FL/R generated by the rotations of the left and right wheels in the water are as follows [27]: (2) where, ρ is the water density, kg/m3; B is the number of individual paddle blade; nL/R is the velocities of the left and right wheels, r/s; θ is the angle of blade immersion depth, deg; V0 is the initial velocity, m/s; KT is the thrust coefficient. When the rotation speeds of the left and right wheels are the same, the paddle boat travels in a straight line. In contrast, the rotation speeds of the left and right wheels are asynchronous, the left and right thrusts of the paddles could be inconsistent, then a torque will be generated at this time. The torque NT caused by the rotation speeds difference can be written as: (3) The torque is controlled by the rudder angle during turning. The correlation between the torque NT and the virtual rudder angle δR of the paddle boat can be expressed as: (4) where, NδT is the hydrodynamic coefficient of the rudder angle. According to the descriptions in literatures [28, 29], the main kinematics variables are the forward traverse velocity u, the traverse velocity v, and the swing angle velocity r, respectively. When the paddle boat is in course holding, trajectory holding or maneuvering below moderate intensity. Consequently, the kinematic equations for the paddle boat are defined as follow: (5) where, X and Y are the hydrodynamic forces along the X0-axis and Y0-axis, respectively; N is the dynamic moment along the Z0-axis; m is the weight of the paddle boat, kg; Iz is the moment of inertia of the hull with respect to the Z0-axis, kg·m2; is the navigation acceleration, m/s2; is the transverse drift acceleration,m/s2;ṙ is the bow angle acceleration,°/s2; xG is the coordinate of gravity of the paddle boat, m. The data from the paddle boat are brought into the hydrodynamic expressions in the literature [26] to obtain the new expressions: (6) where, are the hydrodynamic coefficients. Ignoring the viscous force of water, Eq (5) represents the inertia forces and moments along X0, Y0 and Z0-axis. Usually, the forward motion in the X0-axis direction is independent of the motion in the other two-degree-of-freedom. In addition, from the perspective of velocity control, the motion on the X0-direction degree-of-freedom can be considered separately, while there is a strong coupling between the traverse and the swing, which are closely related to the heading and trajectory control of the paddle boat. Combined with Eq (6), the expression for the paddle boat can be written as: (7) The conversion of the state-space model into a linear response model expressed as a transfer function has been widely used in the design of the heading controllers of USV. The paddle boat moves with low-frequency properties and adopts the Nomoto first-order model (Nomoto) on the basis of neglecting the lateral drift velocity, which facilitates the simulation calculation of the motion. Thus, the transfer function is [30]: (8) where, K and T are the manoeuvring coefficients of the paddle boat. 2.2 Hardware aspect and working process This work implemented a series of experimental testings for validation using the hardware control system shown in Fig 3. It shows the overall component hardware, and the experimental modules include a main controller, a power supply system, a communication system, a positioning system, and a driving system. The paddle boat was able to accomplish the mission of autonomous navigation with the cooperation of various subsystems. These experiments employed a main controller (Pixhawk 4) to generate well-controlled autonomous sailing to study the trajectory tracking of the paddle boat. The power supply system consist of a 2S (S stands for cell connection in series) lead-acid battery module and a current sensor module. In fact, the battery module was capable of generating power for the main controller, GPS, remote control, receiver and other devices. In addition, the main functions of the current sensor module were: (i) monitored the values of the current and voltage. (ii) dropped and regulated the voltage. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Hardware control system. https://doi.org/10.1371/journal.pone.0290246.g003 The communication system used include a receiver, a remote control to operate and switch modes, a data transmission module, and a host computer to connect with the computer software through the Mavlink protocol (mainly used for the unmanned carrier communication). The working process of the communication system was described as follows: We built a connection between the ground side and the control side by using the data transmission module, which could collect 50 data per second, resulting in a high data collection effort. Therefore, we chose a mean filtering algorithm in the data acquisition, and extracted the continuous sampled data per second for accumulation, then assigned the obtained mean value to the data at this time. This type of processing could guarantee the stability of the output data. The formula about the assignment data is given by: (9) where, bi is the i-th sample, n is the number. Then, the host computer sent out the information of the waypoints, position, and velocity, and monitored the track. We designed three main navigation modes, namely hold mode, automatic mode, and manual mode. The different modes could be obtained by switching the host computer or remote control, where corresponded to pulse width modulation signals 0 to1230, 1231 to1360, and 1361 to 1490, respectively. We compiled control algorithms to solve for waypoints and positions, and calculated the course and position that the paddle boat would travel. The core of a positioning system is GPS, which is used to find the exact location of the paddle boat by longitude and latitude. A driving system consist two components: two electronic governors and two motors, they provide driving forces and adjusted the rotational speeds for the left and right wheels, respectively. The process of the speed control could be described as follows: first, to calibrate the motors, we connected the signal lines of the electronic governors to the receiver; second, we linked the signal lines to the pulse position modulation decoder board, aiming to receive the motion command signal and spin the motors. All components were assembled together using the hardware control system, the details of which are shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Hardware model selection. https://doi.org/10.1371/journal.pone.0290246.t001 2.3 Algorithms and control As stated in section 2.1, the velocity and turning of the paddle boat are adjusted by the rotational speed and rotational speed difference between the left and right wheels, respectively. At this point, the turning strategy is to rotate one motor and stop the other one. Additionally, to enable the paddle boat to navigate under the prescribed route, we further develop the trajectory tracking control system in this study, which is mainly composed by two parts, namely LOS guidance system and fuzzy adaptive PID control system, and Fig 4 shows the diagram of the control strategy. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Control strategy of the trajectory tracking. https://doi.org/10.1371/journal.pone.0290246.g004 The trajectory tracking control operation includes the host computer monitoring system, the tracking control system and the execution system, where the tracking control system includes the LOS guidance and the fuzzy adaptive PID control algorithms, as shown in Fig 5. The host computer displays the current heading and tracks the target heading, and the two motors make the relevant commands after being solved by the previously mentioned control algorithm. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Control operation flow. https://doi.org/10.1371/journal.pone.0290246.g005 2.3.1 LOS implementation. The LOS guidance algorithm is able to select line-of-sight points on this section of the track formed by the starting point and the target point when the paddle boat is underway, and controls the course of the paddle boat to navigate in the desired direction. We set the starting point at Pk (Xk,Yk) and the target point Pk+1(Xk+1,Yk+1) in the inertial coordinate system On-XnYnZn, and straight line PkPk+1 was the preset trajectory of the paddle boat, as illustrated in Fig 6. The real-time position of the paddle boat can be expressed as: (10) where, X0 and Y0 are the horizontal and vertical coordinates of the initial position of the paddle boat, respectively; RLOS is the sailing radius, which is a variable, RLOS = nl, and n is the captain ratio. From Fig 6, it could be seen that when the turning radius was set to the variable RLOS, the left and right limit radii could form a sector region. Then we marked PLOS(XLOS,YLOS) as the intersection of this sector region with the preset trajectory.Assuming that the positive direction of the Yn-axis is north, then we record the target heading angle in real-time between the north direction and the straight line OPLOS as Ψr, and Ψr can be written as: (11) Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. LOS guidance principle. https://doi.org/10.1371/journal.pone.0290246.g006 At the same time, we used GPS positioning system to determine the actual heading angle Ψ. Then by calculation we got the deviation angle ΔΨ, which was the difference between Ψr and Ψ. We could also observe from Fig 6 that when the paddle boat entered the circular area at point Pk+1, the target point would be shifted to Pk+2 and at this time, Pk+1 would be associated with Pk+2 by iteration. In this way, the calculation continued in this iterative manner until the computation reached the point Pk+n point. Based on the above LOS algorithm analysis, it was understandable that the position of the PLOS point could be deduced from the real-time position, the starting point and the target point of the paddle boat [31]. 2.3.2 Fuzzy adaptive PID control system. PID control algorithm is a control algorithm that combines proportional, integral and differential links in one, which can be expressed as a function of the deviation as [32]: (12) where, u(t) represents the control output; e(t) represents the deviation between the actual heading angle and the target heading angle; Kpi represents the proportional gain; Ti and Td represent the integration time constant and the differential time constant, respectively. It is no doubt that the three parts of PID, proportional, integral and differential, are linear, thus resulting in a PID controller being a linear controller. However, the paddle boat exhibits a nonlinear characteristic during navigation. If we apply the linear PID controller to the trajectory tracking of the paddle boat, it may reduce the regulation accuracy. At the same time, it needs to adjust the parameters artificially continuously during the control process, which may not applicable to practical situations. Compared with general PID linear controller, the fuzzy adaptive PID control has a better nonlinear characteristic, a higher stability and reliability, can achieve system stability in the fastest time and the least amount of overshoot. Therefore, it is necessary to combine the experience of manual operation with the fuzzy control on the basis of general PID, thus forming the fuzzy PID control algorithm, which can yield more precise results. The fuzzy controller consists of three main parts: the fuzzy input, the fuzzy inference rule, and the fuzzy output. The inputs we set were e and ec, which were the heading error between the given heading and the actual heading as well as the heading error rate, respectively. Additionally, we used the adjusted values of the PID parameters Kp1, Ki1, and Kd1, i.e. ΔKp, ΔKi, and ΔKd, as outputs. We summed the ΔKp, ΔKi, and ΔKd of the fuzzy control output with Kp1, Ki1 and Kd1 of the PID, respectively, to obtain the adjusted Kp, Ki, and Kd, as shown in Eq (13) and Fig 7. (13) where, Kp1, Ki1, and Kd1 are the proportional, integral and differential coefficients of the PID, respectively, which are constantly adjusted by the experimental process of referencing; ΔKp, ΔKi, and ΔKd are the proportional, integral and differential coefficients of the fuzzy control, respectively; Kp, Ki, and Kd are the proportional, integral and differential coefficients of the final action on the paddle boat, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Fuzzy adaptive PID control process. https://doi.org/10.1371/journal.pone.0290246.g007 In the fuzzy adaptive PID control, we designed the fuzzy inference rules using a 7-segment fuzzy method, setting the theoretical domain ranges of both e and ec to [-6, 6], ΔKp, ΔKi, and ΔKd to [-6, 6], [-0.6, 0.6], and [-10, 10], respectively. And, we set the subsets of the fuzzy to {NB, NM, NS, ZO, PS, PM, PB}, the element within the subsets represents {negative large, negative medium, negative small, zero, positive small, positive medium, positive large}. The main rules of the regulation for each period in fuzzy adaptive PID were: (i) in the early stage, the proportional coefficient of the fuzzy control ΔKp should be taken as a larger value to improve the response speed. At the same time, to prevent the integral saturation, the integral coefficient ΔKi should be as small as possible, or even can take zero. In addition, the differential coefficient ΔKd should be taken as a small value to avoid overshooting. (ii) in the mid-regulation, a smaller ΔKp should be taken in order to make the overshoot of the system smaller and ensure a high response rate. Thus, to avoid affecting the stability of the system, ΔKi was taken as an intermediate value. Since the regulation characteristics are sensitive to the change of ΔKd value, the ΔKd value should be appropriately smaller and should be kept constant. (iii) in the post-regulation, larger ΔKp, ΔKi and smaller ΔKd were chosen to reduce the static error, improve the control accuracy and reduce braking effects of the controlled process, respectively. We chose the input and output variables as triangular affiliation functions and determined the corresponding fuzzy rules using the center of gravity method, the details of which are shown in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Fuzzy rules for ΔKp, ΔKi, and ΔKd. https://doi.org/10.1371/journal.pone.0290246.t002 2.1 Model description To establish the kinematics model of the paddle boat, dimensions of the geometric structure used in this study are shown in Fig 1. The basic overall dimensions l×b×h equals 1770mm×1600mm×580mm, the water intrusion depths of the boat and two wheels are d = 110mm and c = 85mm, respectively. Then, the diameter (D) and width (s) of the wheels are 440mm and 110mm, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Geometric dimensions of the paddle boat: (a) rear view; (b)side view. https://doi.org/10.1371/journal.pone.0290246.g001 It is very important to know the acting force between the paddle boat and water. As a general rule, azimuth, pitch and roll angles have little effect on the motion in the horizontal plane. Therefore, we considered only a simplified kinematic model of that is also proposed by Fossen [26] and was characterized by the ease of observing the acting force, which facilitated analytical calculations. The difference between this model and this paper is that the rudder angle is transformed into the steering angle caused by the speed difference of the paddles. Considering the existing longitudinal swing, transverse swing and bow rocking, based on the inertial coordinate system On-XnYnZn fixed on the Earth’s surface and hull coordinate system O-X0Y0Z0 fixed on the paddle boat, we set up a three-degree-of-freedom kinematics model, as shown in Fig 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Schematic diagram of inertial coordinate system and hull coordinate system. https://doi.org/10.1371/journal.pone.0290246.g002 The transformation between the inertial coordinate system and the hull coordinate system can be expressed as: (1) where, x, y, and Ψ represent the x-point, y-point, and the bow angle of the paddle boat in the inertial coordinate system, respectively; u, v, and r represent the forward velocity, traverse velocity, and the swing angle velocity in the hull coordinate system, respectively. In solving the simulation model of the paddle boat, we processed the velocities of the left and right wheels as input parameters and the bow angle as an output parameter. Afterward, we analyzed the propulsive forces of the paddle wheels. According to Newton’s second law of motion, the propulsive forces FL/R generated by the rotations of the left and right wheels in the water are as follows [27]: (2) where, ρ is the water density, kg/m3; B is the number of individual paddle blade; nL/R is the velocities of the left and right wheels, r/s; θ is the angle of blade immersion depth, deg; V0 is the initial velocity, m/s; KT is the thrust coefficient. When the rotation speeds of the left and right wheels are the same, the paddle boat travels in a straight line. In contrast, the rotation speeds of the left and right wheels are asynchronous, the left and right thrusts of the paddles could be inconsistent, then a torque will be generated at this time. The torque NT caused by the rotation speeds difference can be written as: (3) The torque is controlled by the rudder angle during turning. The correlation between the torque NT and the virtual rudder angle δR of the paddle boat can be expressed as: (4) where, NδT is the hydrodynamic coefficient of the rudder angle. According to the descriptions in literatures [28, 29], the main kinematics variables are the forward traverse velocity u, the traverse velocity v, and the swing angle velocity r, respectively. When the paddle boat is in course holding, trajectory holding or maneuvering below moderate intensity. Consequently, the kinematic equations for the paddle boat are defined as follow: (5) where, X and Y are the hydrodynamic forces along the X0-axis and Y0-axis, respectively; N is the dynamic moment along the Z0-axis; m is the weight of the paddle boat, kg; Iz is the moment of inertia of the hull with respect to the Z0-axis, kg·m2; is the navigation acceleration, m/s2; is the transverse drift acceleration,m/s2;ṙ is the bow angle acceleration,°/s2; xG is the coordinate of gravity of the paddle boat, m. The data from the paddle boat are brought into the hydrodynamic expressions in the literature [26] to obtain the new expressions: (6) where, are the hydrodynamic coefficients. Ignoring the viscous force of water, Eq (5) represents the inertia forces and moments along X0, Y0 and Z0-axis. Usually, the forward motion in the X0-axis direction is independent of the motion in the other two-degree-of-freedom. In addition, from the perspective of velocity control, the motion on the X0-direction degree-of-freedom can be considered separately, while there is a strong coupling between the traverse and the swing, which are closely related to the heading and trajectory control of the paddle boat. Combined with Eq (6), the expression for the paddle boat can be written as: (7) The conversion of the state-space model into a linear response model expressed as a transfer function has been widely used in the design of the heading controllers of USV. The paddle boat moves with low-frequency properties and adopts the Nomoto first-order model (Nomoto) on the basis of neglecting the lateral drift velocity, which facilitates the simulation calculation of the motion. Thus, the transfer function is [30]: (8) where, K and T are the manoeuvring coefficients of the paddle boat. 2.2 Hardware aspect and working process This work implemented a series of experimental testings for validation using the hardware control system shown in Fig 3. It shows the overall component hardware, and the experimental modules include a main controller, a power supply system, a communication system, a positioning system, and a driving system. The paddle boat was able to accomplish the mission of autonomous navigation with the cooperation of various subsystems. These experiments employed a main controller (Pixhawk 4) to generate well-controlled autonomous sailing to study the trajectory tracking of the paddle boat. The power supply system consist of a 2S (S stands for cell connection in series) lead-acid battery module and a current sensor module. In fact, the battery module was capable of generating power for the main controller, GPS, remote control, receiver and other devices. In addition, the main functions of the current sensor module were: (i) monitored the values of the current and voltage. (ii) dropped and regulated the voltage. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Hardware control system. https://doi.org/10.1371/journal.pone.0290246.g003 The communication system used include a receiver, a remote control to operate and switch modes, a data transmission module, and a host computer to connect with the computer software through the Mavlink protocol (mainly used for the unmanned carrier communication). The working process of the communication system was described as follows: We built a connection between the ground side and the control side by using the data transmission module, which could collect 50 data per second, resulting in a high data collection effort. Therefore, we chose a mean filtering algorithm in the data acquisition, and extracted the continuous sampled data per second for accumulation, then assigned the obtained mean value to the data at this time. This type of processing could guarantee the stability of the output data. The formula about the assignment data is given by: (9) where, bi is the i-th sample, n is the number. Then, the host computer sent out the information of the waypoints, position, and velocity, and monitored the track. We designed three main navigation modes, namely hold mode, automatic mode, and manual mode. The different modes could be obtained by switching the host computer or remote control, where corresponded to pulse width modulation signals 0 to1230, 1231 to1360, and 1361 to 1490, respectively. We compiled control algorithms to solve for waypoints and positions, and calculated the course and position that the paddle boat would travel. The core of a positioning system is GPS, which is used to find the exact location of the paddle boat by longitude and latitude. A driving system consist two components: two electronic governors and two motors, they provide driving forces and adjusted the rotational speeds for the left and right wheels, respectively. The process of the speed control could be described as follows: first, to calibrate the motors, we connected the signal lines of the electronic governors to the receiver; second, we linked the signal lines to the pulse position modulation decoder board, aiming to receive the motion command signal and spin the motors. All components were assembled together using the hardware control system, the details of which are shown in Table 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Hardware model selection. https://doi.org/10.1371/journal.pone.0290246.t001 2.3 Algorithms and control As stated in section 2.1, the velocity and turning of the paddle boat are adjusted by the rotational speed and rotational speed difference between the left and right wheels, respectively. At this point, the turning strategy is to rotate one motor and stop the other one. Additionally, to enable the paddle boat to navigate under the prescribed route, we further develop the trajectory tracking control system in this study, which is mainly composed by two parts, namely LOS guidance system and fuzzy adaptive PID control system, and Fig 4 shows the diagram of the control strategy. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Control strategy of the trajectory tracking. https://doi.org/10.1371/journal.pone.0290246.g004 The trajectory tracking control operation includes the host computer monitoring system, the tracking control system and the execution system, where the tracking control system includes the LOS guidance and the fuzzy adaptive PID control algorithms, as shown in Fig 5. The host computer displays the current heading and tracks the target heading, and the two motors make the relevant commands after being solved by the previously mentioned control algorithm. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Control operation flow. https://doi.org/10.1371/journal.pone.0290246.g005 2.3.1 LOS implementation. The LOS guidance algorithm is able to select line-of-sight points on this section of the track formed by the starting point and the target point when the paddle boat is underway, and controls the course of the paddle boat to navigate in the desired direction. We set the starting point at Pk (Xk,Yk) and the target point Pk+1(Xk+1,Yk+1) in the inertial coordinate system On-XnYnZn, and straight line PkPk+1 was the preset trajectory of the paddle boat, as illustrated in Fig 6. The real-time position of the paddle boat can be expressed as: (10) where, X0 and Y0 are the horizontal and vertical coordinates of the initial position of the paddle boat, respectively; RLOS is the sailing radius, which is a variable, RLOS = nl, and n is the captain ratio. From Fig 6, it could be seen that when the turning radius was set to the variable RLOS, the left and right limit radii could form a sector region. Then we marked PLOS(XLOS,YLOS) as the intersection of this sector region with the preset trajectory.Assuming that the positive direction of the Yn-axis is north, then we record the target heading angle in real-time between the north direction and the straight line OPLOS as Ψr, and Ψr can be written as: (11) Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. LOS guidance principle. https://doi.org/10.1371/journal.pone.0290246.g006 At the same time, we used GPS positioning system to determine the actual heading angle Ψ. Then by calculation we got the deviation angle ΔΨ, which was the difference between Ψr and Ψ. We could also observe from Fig 6 that when the paddle boat entered the circular area at point Pk+1, the target point would be shifted to Pk+2 and at this time, Pk+1 would be associated with Pk+2 by iteration. In this way, the calculation continued in this iterative manner until the computation reached the point Pk+n point. Based on the above LOS algorithm analysis, it was understandable that the position of the PLOS point could be deduced from the real-time position, the starting point and the target point of the paddle boat [31]. 2.3.2 Fuzzy adaptive PID control system. PID control algorithm is a control algorithm that combines proportional, integral and differential links in one, which can be expressed as a function of the deviation as [32]: (12) where, u(t) represents the control output; e(t) represents the deviation between the actual heading angle and the target heading angle; Kpi represents the proportional gain; Ti and Td represent the integration time constant and the differential time constant, respectively. It is no doubt that the three parts of PID, proportional, integral and differential, are linear, thus resulting in a PID controller being a linear controller. However, the paddle boat exhibits a nonlinear characteristic during navigation. If we apply the linear PID controller to the trajectory tracking of the paddle boat, it may reduce the regulation accuracy. At the same time, it needs to adjust the parameters artificially continuously during the control process, which may not applicable to practical situations. Compared with general PID linear controller, the fuzzy adaptive PID control has a better nonlinear characteristic, a higher stability and reliability, can achieve system stability in the fastest time and the least amount of overshoot. Therefore, it is necessary to combine the experience of manual operation with the fuzzy control on the basis of general PID, thus forming the fuzzy PID control algorithm, which can yield more precise results. The fuzzy controller consists of three main parts: the fuzzy input, the fuzzy inference rule, and the fuzzy output. The inputs we set were e and ec, which were the heading error between the given heading and the actual heading as well as the heading error rate, respectively. Additionally, we used the adjusted values of the PID parameters Kp1, Ki1, and Kd1, i.e. ΔKp, ΔKi, and ΔKd, as outputs. We summed the ΔKp, ΔKi, and ΔKd of the fuzzy control output with Kp1, Ki1 and Kd1 of the PID, respectively, to obtain the adjusted Kp, Ki, and Kd, as shown in Eq (13) and Fig 7. (13) where, Kp1, Ki1, and Kd1 are the proportional, integral and differential coefficients of the PID, respectively, which are constantly adjusted by the experimental process of referencing; ΔKp, ΔKi, and ΔKd are the proportional, integral and differential coefficients of the fuzzy control, respectively; Kp, Ki, and Kd are the proportional, integral and differential coefficients of the final action on the paddle boat, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Fuzzy adaptive PID control process. https://doi.org/10.1371/journal.pone.0290246.g007 In the fuzzy adaptive PID control, we designed the fuzzy inference rules using a 7-segment fuzzy method, setting the theoretical domain ranges of both e and ec to [-6, 6], ΔKp, ΔKi, and ΔKd to [-6, 6], [-0.6, 0.6], and [-10, 10], respectively. And, we set the subsets of the fuzzy to {NB, NM, NS, ZO, PS, PM, PB}, the element within the subsets represents {negative large, negative medium, negative small, zero, positive small, positive medium, positive large}. The main rules of the regulation for each period in fuzzy adaptive PID were: (i) in the early stage, the proportional coefficient of the fuzzy control ΔKp should be taken as a larger value to improve the response speed. At the same time, to prevent the integral saturation, the integral coefficient ΔKi should be as small as possible, or even can take zero. In addition, the differential coefficient ΔKd should be taken as a small value to avoid overshooting. (ii) in the mid-regulation, a smaller ΔKp should be taken in order to make the overshoot of the system smaller and ensure a high response rate. Thus, to avoid affecting the stability of the system, ΔKi was taken as an intermediate value. Since the regulation characteristics are sensitive to the change of ΔKd value, the ΔKd value should be appropriately smaller and should be kept constant. (iii) in the post-regulation, larger ΔKp, ΔKi and smaller ΔKd were chosen to reduce the static error, improve the control accuracy and reduce braking effects of the controlled process, respectively. We chose the input and output variables as triangular affiliation functions and determined the corresponding fuzzy rules using the center of gravity method, the details of which are shown in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Fuzzy rules for ΔKp, ΔKi, and ΔKd. https://doi.org/10.1371/journal.pone.0290246.t002 2.3.1 LOS implementation. The LOS guidance algorithm is able to select line-of-sight points on this section of the track formed by the starting point and the target point when the paddle boat is underway, and controls the course of the paddle boat to navigate in the desired direction. We set the starting point at Pk (Xk,Yk) and the target point Pk+1(Xk+1,Yk+1) in the inertial coordinate system On-XnYnZn, and straight line PkPk+1 was the preset trajectory of the paddle boat, as illustrated in Fig 6. The real-time position of the paddle boat can be expressed as: (10) where, X0 and Y0 are the horizontal and vertical coordinates of the initial position of the paddle boat, respectively; RLOS is the sailing radius, which is a variable, RLOS = nl, and n is the captain ratio. From Fig 6, it could be seen that when the turning radius was set to the variable RLOS, the left and right limit radii could form a sector region. Then we marked PLOS(XLOS,YLOS) as the intersection of this sector region with the preset trajectory.Assuming that the positive direction of the Yn-axis is north, then we record the target heading angle in real-time between the north direction and the straight line OPLOS as Ψr, and Ψr can be written as: (11) Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. LOS guidance principle. https://doi.org/10.1371/journal.pone.0290246.g006 At the same time, we used GPS positioning system to determine the actual heading angle Ψ. Then by calculation we got the deviation angle ΔΨ, which was the difference between Ψr and Ψ. We could also observe from Fig 6 that when the paddle boat entered the circular area at point Pk+1, the target point would be shifted to Pk+2 and at this time, Pk+1 would be associated with Pk+2 by iteration. In this way, the calculation continued in this iterative manner until the computation reached the point Pk+n point. Based on the above LOS algorithm analysis, it was understandable that the position of the PLOS point could be deduced from the real-time position, the starting point and the target point of the paddle boat [31]. 2.3.2 Fuzzy adaptive PID control system. PID control algorithm is a control algorithm that combines proportional, integral and differential links in one, which can be expressed as a function of the deviation as [32]: (12) where, u(t) represents the control output; e(t) represents the deviation between the actual heading angle and the target heading angle; Kpi represents the proportional gain; Ti and Td represent the integration time constant and the differential time constant, respectively. It is no doubt that the three parts of PID, proportional, integral and differential, are linear, thus resulting in a PID controller being a linear controller. However, the paddle boat exhibits a nonlinear characteristic during navigation. If we apply the linear PID controller to the trajectory tracking of the paddle boat, it may reduce the regulation accuracy. At the same time, it needs to adjust the parameters artificially continuously during the control process, which may not applicable to practical situations. Compared with general PID linear controller, the fuzzy adaptive PID control has a better nonlinear characteristic, a higher stability and reliability, can achieve system stability in the fastest time and the least amount of overshoot. Therefore, it is necessary to combine the experience of manual operation with the fuzzy control on the basis of general PID, thus forming the fuzzy PID control algorithm, which can yield more precise results. The fuzzy controller consists of three main parts: the fuzzy input, the fuzzy inference rule, and the fuzzy output. The inputs we set were e and ec, which were the heading error between the given heading and the actual heading as well as the heading error rate, respectively. Additionally, we used the adjusted values of the PID parameters Kp1, Ki1, and Kd1, i.e. ΔKp, ΔKi, and ΔKd, as outputs. We summed the ΔKp, ΔKi, and ΔKd of the fuzzy control output with Kp1, Ki1 and Kd1 of the PID, respectively, to obtain the adjusted Kp, Ki, and Kd, as shown in Eq (13) and Fig 7. (13) where, Kp1, Ki1, and Kd1 are the proportional, integral and differential coefficients of the PID, respectively, which are constantly adjusted by the experimental process of referencing; ΔKp, ΔKi, and ΔKd are the proportional, integral and differential coefficients of the fuzzy control, respectively; Kp, Ki, and Kd are the proportional, integral and differential coefficients of the final action on the paddle boat, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Fuzzy adaptive PID control process. https://doi.org/10.1371/journal.pone.0290246.g007 In the fuzzy adaptive PID control, we designed the fuzzy inference rules using a 7-segment fuzzy method, setting the theoretical domain ranges of both e and ec to [-6, 6], ΔKp, ΔKi, and ΔKd to [-6, 6], [-0.6, 0.6], and [-10, 10], respectively. And, we set the subsets of the fuzzy to {NB, NM, NS, ZO, PS, PM, PB}, the element within the subsets represents {negative large, negative medium, negative small, zero, positive small, positive medium, positive large}. The main rules of the regulation for each period in fuzzy adaptive PID were: (i) in the early stage, the proportional coefficient of the fuzzy control ΔKp should be taken as a larger value to improve the response speed. At the same time, to prevent the integral saturation, the integral coefficient ΔKi should be as small as possible, or even can take zero. In addition, the differential coefficient ΔKd should be taken as a small value to avoid overshooting. (ii) in the mid-regulation, a smaller ΔKp should be taken in order to make the overshoot of the system smaller and ensure a high response rate. Thus, to avoid affecting the stability of the system, ΔKi was taken as an intermediate value. Since the regulation characteristics are sensitive to the change of ΔKd value, the ΔKd value should be appropriately smaller and should be kept constant. (iii) in the post-regulation, larger ΔKp, ΔKi and smaller ΔKd were chosen to reduce the static error, improve the control accuracy and reduce braking effects of the controlled process, respectively. We chose the input and output variables as triangular affiliation functions and determined the corresponding fuzzy rules using the center of gravity method, the details of which are shown in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Fuzzy rules for ΔKp, ΔKi, and ΔKd. https://doi.org/10.1371/journal.pone.0290246.t002 3 Results and discussion Fig 8 shows a test operating platform. We tested the trajectory tracking of the paddle boat in a pond, with a size of 9,000mm×100,000mm, and the weather on the day of the test was sunny, with the temperature ranging from 10 to 25°C and wind force 2 to 3. According to the software and hardware platform built in the previous section, we also created the host monitoring for the Mission Planner platform, a ground station software of open access developed in C#. Under these conditions, we set the radius of the waypoint to 2,000mm and the visible radius of the LOS guidance algorithm to 3,000mm. To further verify the applicability of the trajectory tracking control system, we carried out the tests in two ways. One was to compare the trajectory tracking at different velocities on the same route to select the optimal velocity, the other was to compare different turning angles under various preset trajectories to observe the effect of the trajectory tracking. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Test operating platform. https://doi.org/10.1371/journal.pone.0290246.g008 3.1 Trajectories at different velocities As you know, when the velocity of the paddle boat is relatively fast, it is easy to be disturbed by wind and waves, even cause deviation of the course. In contrast, when the velocity of the paddle boat is relatively slow, it may affect the efficiency of the work. Therefore, we selected the velocities of 0.5m/s, 1m/s and 1.5m/s for the comparison tests of the trajectory tracking on the same route. Fig 9 shows the comparison of the trajectory tracking under different velocities. In Fig 9(A), the position of the paddle boat was the origin, i.e. the starting point. The green marker indicated the first few target points, the red line was the current heading of the paddle boat, the black line was the GPS tracking line, and the orange line expressed the next waypoint index. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Comparison of trajectory tracking at different velocities on the same route:(a) A straight line reference trajectory; (b) Actual trajectory at a speed of 0.5m/s; (c) Actual trajectory at a speed of 1m/s; (d) Actual trajectory at a speed of 1.5m/s; (e) Comparison between the reference trajectory and the actual tracks. https://doi.org/10.1371/journal.pone.0290246.g009 As shown in Fig 9, the overall trend of trajectory were consistent obviously within the error allowed. Apparently, when the velocity increased, the tracking error also increased gradually. However, we do not consider the effect of wind and waves on the heading in our algorithm, increasing boat velocity will lead to the adaptive PID control algorithm to fail to keep up with the convergence rate. And yet, when the heading error e between the given heading and the actual heading reached the set value, the LOS guidance algorithm could assist in turning and correcting the deflection to make the paddle boat towards the target waypoint, proving the stability of this control system. By comparing the cases above, it can be seen that the actual trajectory of the paddle boat at 0.5m/s matches better with the reference trajectory. Therefore, 0.5m/s was chosen for the subsequent tests. 3.2 Trajectories at different turning angles In the actual application scenario, the paddle boat is required to adapt to diverse shapes and areas of ponds, resulting in accurate tracking under different turning angles. To be able to perform tasks in different scenarios, it is necessary to verify the turning angle for the accurate trajectory tracking, see Fig 10.Through the actual measurement we found that the small boat turning was very flexible and its turning radius was about 2 m. Here, we set the turning angles as acute (about 30°), right angle (about 90°) and obtuse angle (about 120°), respectively, and studied the tracking of the reference trajectory by the paddle boat when the reference routes were triangle, square and pentagon, respectively. The test results were shown in Figs 11–13. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Left turn of the experimental boat: (a) Actual trajectory when testing turning radius; (b) Extracted left-turn trajectory. https://doi.org/10.1371/journal.pone.0290246.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Trajectory of a triangle with an approximate 60° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Trajectory of a square with an approximate 90° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g012 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Trajectory of a pentagon with an approximate 120° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g013 From the Figs 11–13, it could be observed that with the LOS guidance algorithm and fuzzy adaptive PID control system, the paddle boat could achieve the trajectory tracking under different turning angles with a fixed velocity of 0.5m/s. By comparing the trajectory tracking errors with those under different literatures, it is found that the average error in this paper tracks the best and the maximum error is relatively good, as shown in Table 3. We expand on this with pentagon trajectory tracking and make an experimental analysis of trajectory tracking under different shapes. The actual trajectory and the reference trajectory were in good agreement, but had slightly difference at different turning angles. The reasons are mainly the following two parts: (1) the distance between points Pk and PK+1 may be shorter than the turning radius. Therefore, even if the minimum turning radius is maintained, the actual trajectory may still deviate from the reference trajectory. (2) the fluid and wind in the actual test will affect the experimental results and thus lead to deviation from the reference trajectory. However, when the distance is longer than the turning radius, the LOS guidance algorithm will automatically solve the next tracking point and use fuzzy adaptive PID control to correct the heading, and thereby to achieve trajectory tracking. Overall, the navigation error is within the allowed range, which proves that the control system has high reliability. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Comparison of trajectory tracking under different literatures. https://doi.org/10.1371/journal.pone.0290246.t003 3.1 Trajectories at different velocities As you know, when the velocity of the paddle boat is relatively fast, it is easy to be disturbed by wind and waves, even cause deviation of the course. In contrast, when the velocity of the paddle boat is relatively slow, it may affect the efficiency of the work. Therefore, we selected the velocities of 0.5m/s, 1m/s and 1.5m/s for the comparison tests of the trajectory tracking on the same route. Fig 9 shows the comparison of the trajectory tracking under different velocities. In Fig 9(A), the position of the paddle boat was the origin, i.e. the starting point. The green marker indicated the first few target points, the red line was the current heading of the paddle boat, the black line was the GPS tracking line, and the orange line expressed the next waypoint index. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Comparison of trajectory tracking at different velocities on the same route:(a) A straight line reference trajectory; (b) Actual trajectory at a speed of 0.5m/s; (c) Actual trajectory at a speed of 1m/s; (d) Actual trajectory at a speed of 1.5m/s; (e) Comparison between the reference trajectory and the actual tracks. https://doi.org/10.1371/journal.pone.0290246.g009 As shown in Fig 9, the overall trend of trajectory were consistent obviously within the error allowed. Apparently, when the velocity increased, the tracking error also increased gradually. However, we do not consider the effect of wind and waves on the heading in our algorithm, increasing boat velocity will lead to the adaptive PID control algorithm to fail to keep up with the convergence rate. And yet, when the heading error e between the given heading and the actual heading reached the set value, the LOS guidance algorithm could assist in turning and correcting the deflection to make the paddle boat towards the target waypoint, proving the stability of this control system. By comparing the cases above, it can be seen that the actual trajectory of the paddle boat at 0.5m/s matches better with the reference trajectory. Therefore, 0.5m/s was chosen for the subsequent tests. 3.2 Trajectories at different turning angles In the actual application scenario, the paddle boat is required to adapt to diverse shapes and areas of ponds, resulting in accurate tracking under different turning angles. To be able to perform tasks in different scenarios, it is necessary to verify the turning angle for the accurate trajectory tracking, see Fig 10.Through the actual measurement we found that the small boat turning was very flexible and its turning radius was about 2 m. Here, we set the turning angles as acute (about 30°), right angle (about 90°) and obtuse angle (about 120°), respectively, and studied the tracking of the reference trajectory by the paddle boat when the reference routes were triangle, square and pentagon, respectively. The test results were shown in Figs 11–13. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Left turn of the experimental boat: (a) Actual trajectory when testing turning radius; (b) Extracted left-turn trajectory. https://doi.org/10.1371/journal.pone.0290246.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Trajectory of a triangle with an approximate 60° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Trajectory of a square with an approximate 90° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g012 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Trajectory of a pentagon with an approximate 120° turning angle: (a) Reference trajectory; (b) Actual trajectory; (c) Comparison of the reference trajectory and the actual trajectory. https://doi.org/10.1371/journal.pone.0290246.g013 From the Figs 11–13, it could be observed that with the LOS guidance algorithm and fuzzy adaptive PID control system, the paddle boat could achieve the trajectory tracking under different turning angles with a fixed velocity of 0.5m/s. By comparing the trajectory tracking errors with those under different literatures, it is found that the average error in this paper tracks the best and the maximum error is relatively good, as shown in Table 3. We expand on this with pentagon trajectory tracking and make an experimental analysis of trajectory tracking under different shapes. The actual trajectory and the reference trajectory were in good agreement, but had slightly difference at different turning angles. The reasons are mainly the following two parts: (1) the distance between points Pk and PK+1 may be shorter than the turning radius. Therefore, even if the minimum turning radius is maintained, the actual trajectory may still deviate from the reference trajectory. (2) the fluid and wind in the actual test will affect the experimental results and thus lead to deviation from the reference trajectory. However, when the distance is longer than the turning radius, the LOS guidance algorithm will automatically solve the next tracking point and use fuzzy adaptive PID control to correct the heading, and thereby to achieve trajectory tracking. Overall, the navigation error is within the allowed range, which proves that the control system has high reliability. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Comparison of trajectory tracking under different literatures. https://doi.org/10.1371/journal.pone.0290246.t003 4 Conclusions The present paper reports a trajectory tracking control system for a paddle boat. In addition, we successfully develops a test operating platform based on Mission Planner for monitoring navigation track in the paddle boat and providing control method recommendations using the LOS guidance and fuzzy adaptive PID algorithm. We also test the trajectories at different velocities and turning angles successfully using the software and hardware control system and evaluate using preciseness. The significant observations based on the present investigation are listed below: A three-degree-of-freedom kinematics model is set up to study the navigation of the paddle boat in the water. Moreover, the simplified model can describe the relationship between the velocities, forces and moments of the paddle boat. The overall hardware includes a main controller, a power supply system, a communication system, a positioning system, and a driving system. All components are assembled together and work with each other to perform the task of autonomous navigation. The LOS guidance algorithm is used to calculate the LOS angle of the preset track, and the waypoint is tracked and analyzed. Simultaneously, the heading is controlled by fuzzy adaptive PID, and the parameters of the PID are constantly corrected through the transformation of the paddle boat’s heading to make the trajectory tracking more accurate. By testing the trajectory tracking effect of different velocities and turning angles under the same route, it is found that the system has good tracking effect under 0.5m/s, and can execute each waypoint stably according to the preset trajectory with good control performance. In addition, the paddle boat in this study has a good turning effect. Furthermore, these results are expected to provide some significant physical insights for the design and trajectory tracking optimization of such paddle boat. As a perspective for this work, the disturbance of the system by noise such as wind, fluid velocity and water waves cannot be neglected and there is a strong coupling with the motion of the paddle boat, such system will be investigated and analyzed in the near future based on the developed model. Supporting information S1 Data. https://doi.org/10.1371/journal.pone.0290246.s001 (XLSX) S1 Raw images. https://doi.org/10.1371/journal.pone.0290246.s002 (PDF) TI - A trajectory tracking control system for paddle boat in intelligent aquaculture JF - PLoS ONE DO - 10.1371/journal.pone.0290246 DA - 2023-08-17 UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/a-trajectory-tracking-control-system-for-paddle-boat-in-intelligent-WQuwgERY0V SP - e0290246 VL - 18 IS - 8 DP - DeepDyve ER -