Flight Test of Autonomous Formation Management for Multiple Fixed-Wing UAVs Based on Missile Parallel Method

: This paper reports on the formation and transformation of multiple ﬁxed-wing unmanned aerial vehicles (UAVs) in three-dimensional space. A cooperative guidance law based on the classic missile-type parallel-approach method is designed for the multi-UAV formation control problem. Additionally, formation transformation strategies for multi-UAV autonomous assembly, disbandment, and special circumstances are formed, effective for managing and controlling the formation. When formulating the management strategy for formation establishment, its process is divided into three steps: (i) selecting and allocating target points, (ii) forming loose formations, and (iii) forming short-range formations. The management of disbanding the formation is formulated through reverse thinking: the assembly process is split and recombined in reverse, and a formation disbanding strategy that can achieve a smooth transition from close to lose formation is proposed. Additionally, a strategy is given for adjusting the formation transformation in special cases, and the formation adjustment is completed using the adjacency matrix. Finally, a hardware-in-the-loop simulation and measured ﬂight veriﬁcation using a simulator show the practicality of the guidance law in meeting the control requirements of UAV formation ﬂight for speciﬁc ﬂight tasks.


Introduction
Form formation management is an important research topic for multiple unmanned aerial vehicles (UAVs) flying in cooperative formation [1,2]. Because of the limited control of fixed-wing UAVs, their formation management is very different from that of quadrotor UAVs. By contrast, static and fixed-wing UAVs cannot wait for other UAVs during formation [3], which makes stable and reliable formation challenging.
In recent decades, many scholars have researched the formation of fixed-wing UAVs. Zhang [4] proposed a proportional integral-derivative (PID) integrated control method based on robust control; this comprehensive control method improves the hit rate and flight stability of UAV formation, reduces the dynamic response of the steady-state error, and shortens the convergence time, but the influence of the coupling effect of the integrated controller on the entire control system was not considered in the corresponding experiment. Kada [5] proposed (i) a smooth distributed cooperative control method for multiple aircraft (such as UAVs) based on multi-agent system (MAS) consistency and (ii) a smooth distributed consistency algorithm, as well as designing a formation control model for three-dimensional (3D) geometry tracking; however, the disadvantage of this approach is that it fails to consider (i) formation tracking in the case of external interference and (ii) obstacle avoidance among flying agents. Wang [6] analyzed the optimality of formation configuration, provided an optimal formation design strategy for multi-UAV patrol tasks, and developed a model prediction trajectory replanning algorithm. However, the disadvantages of this approach are that it (i) considers only the mission and motion planning level and (ii) does not involve low-level flight-control problems. Reference [7] proposed a robust decentralized tracking control scheme for a large-scale UAV-formation network control system and verified it in such a system.
Research on formation maintenance and transformation remains sparse. Proposed herein is a formation management method based on cooperative waypoint following. Additionally, formation reconfiguration strategy is designed considering how some typical unexpected situations influence cooperative formations. Combining the guidance law of the parallel-approach method [8,9] and the calculation method of virtual dynamic tracking points, a method is proposed for cooperative path-point following in formation control [10][11][12] to realize multi-machine formation flight of fixed-wing UAVs. It has the advantages of fast tracking response in the low-altitude formation flight, and rapid formation reconstruction in unexpected situations. Finally, the stability and agility of the guidance law in formation cooperative control are verified by hardware-in-the-loop [13][14][15] simulation and flight-test measurements.
Based on the flight control of a single fixed-wing aircraft [13][14][15], this paper designs a layered control framework, as shown in Figure 1. The control framework includes four submodules: the intermachine communication module, the cooperative waypointfollowing control module, the cooperative formation module, and the stabilization loop control module. After receiving the instruction to start the formation, the cooperative formation module generates the desired waypoints for the UAVs. Compared with trajectory tracking, waypoint tracking is more widely used in aviation, and a series of static waypoints often represent common single-plane flight tasks. Waypoint tracking is used more widely in aviation, and a series of static waypoints often represent common single-plane flight tasks.
Drones 2022, 6, x FOR PEER REVIEW 2 of 26 patrol tasks, and developed a model prediction trajectory replanning algorithm. However, the disadvantages of this approach are that it (i) considers only the mission and motion planning level and (ii) does not involve low-level flight-control problems. Reference [7] proposed a robust decentralized tracking control scheme for a large-scale UAV-formation network control system and verified it in such a system. Research on formation maintenance and transformation remains sparse. Proposed herein is a formation management method based on cooperative waypoint following. Additionally, formation reconfiguration strategy is designed considering how some typical unexpected situations influence cooperative formations. Combining the guidance law of the parallel-approach method [8,9] and the calculation method of virtual dynamic tracking points, a method is proposed for cooperative path-point following in formation control [10][11][12] to realize multi-machine formation flight of fixed-wing UAVs. It has the advantages of fast tracking response in the low-altitude formation flight, and rapid formation reconstruction in unexpected situations. Finally, the stability and agility of the guidance law in formation cooperative control are verified by hardware-in-the-loop [13][14][15] simulation and flight-test measurements.
Based on the flight control of a single fixed-wing aircraft [13][14][15], this paper designs a layered control framework, as shown in Figure 1. The control framework includes four submodules: the intermachine communication module, the cooperative waypoint-following control module, the cooperative formation module, and the stabilization loop control module. After receiving the instruction to start the formation, the cooperative formation module generates the desired waypoints for the UAVs. Compared with trajectory tracking, waypoint tracking is more widely used in aviation, and a series of static waypoints often represent common single-plane flight tasks. Waypoint tracking is used more widely in aviation, and a series of static waypoints often represent common single-plane flight tasks.

Design of Guidance Law
The required lateral movement is determined by calculating virtual dynamic tracking points according to the formation relative distance of the formation, and the virtual dynamic tracking points are those of the wingmen.
Assuming that the UAVs are all mass points when flying in formation, as shown in Figure 2, in the case of two UAVs, let point A be where the leader is located, let point B be where the virtual dynamic tracking point is located, and let point C be where the wingman is located. Define the coordinate formation system fixed to the lead plane

Design of Guidance Law
The required lateral movement is determined by calculating virtual dynamic tracking points according to the formation relative distance of the formation, and the virtual dynamic tracking points are those of the wingmen.
Assuming that the UAVs are all mass points when flying in formation, as shown in Figure 2, in the case of two UAVs, let point A be where the leader is located, let point B be where the virtual dynamic tracking point is located, and let point C be where the wingman is located. Define the coordinate formation system O f X f Y f Z f fixed to the lead plane L:O f X f is parallel to the direction of the speed vector of the lead plane and points forward; the axis O f Y f is perpendicular to O f X f and points to the right of the flight direction; the axis O f Z f is perpendicular to O f X f Y f and points downward. Points A and B are on concentric Drones 2022, 6, 99 3 of 24 circular arcs whose lengths are ∆x and arc, respectively, and the concentric circular arcs are on O f X f Y f ; ∆y is the length of the straight line of the leader aircraft to the arc, R L is the flight radius of the lead aircraft, R is the radius of the concentric circular arcs, and ϕ is the included angle of point A and B. In this coordinate system, point B is relative to point A and its positional differences are ∆x err and ∆y err , whose lengths can be obtained geometrically.
Drones 2022, 6, x FOR PEER REVIEW 3 of 26 L: is parallel to the direction of the speed vector of the lead plane and points forward; the axis is perpendicular to and points to the right of the flight direction; the axis is perpendicular to and points downward. Points A and B are on concentric circular arcs whose lengths are x  and arc, respectively, and the concentric circular arcs are on ; y  is the length of the straight line of the leader aircraft to the arc, is the flight radius of the lead aircraft, R is the radius of the concentric circular arcs, and φ is the included angle of point A and B. In this coordinate system, point B is relative to point A and its positional differences are err x  and err y  , whose lengths can be obtained geometrically.
where L  is the angle between the two vector axes f Ox and g Ox .
Then, the coordinate difference between points A and B in the ground coordinate system is When ϕ is small, we have ∆x err = arc ∆y err = ∆y When ϕ is larger, we have where arc is calculated as Therefore, the distance difference [x err y err ] T between the virtual dynamic tracking point B and the leader A can be obtained through the distance [∆x ∆y] T set for the formation.
The transformation matrix for converting from the formation coordinate system x f O f y f to the ground coordinate system x g O g y g is where ψ L is the angle between the two vector axes Ox f and Ox g . Then, the coordinate difference between points A and B in the ground coordinate system is x err y err = L g f x err y err As can be seen, the coordinates of point B are Letting the distance difference between points B and C be [P x P y ] T , we obtain The key to designing the guidance law for the wingmen in the formation is to convert the relative distances between the UAVs and the target tracking points, the velocity vector deviations, the heading deviations, and other information into guidance commands. The algorithm for the guidance law can be designed using either geometric methods or controltheory methods. The former use the geometric relationship between UAV and target point to design the guidance law, characterized by convenient analysis and simulation; geometric methods are the ones used most commonly in engineering practice.
The geometric parallel-approach method is used as the guidance law for cooperative formation. The parallel-approach method was designed originally for missile interception and attacking targets, and the resulting trajectories are relatively straight. The specific implementation of the parallel-approach method herein is to keep the line of sight from UAV to target point unchanged during flight, and the line-of-sight angle forms a set of parallel lines in space, as shown in Figure 3.
As can be seen, the coordinates of point B are Letting the distance difference between points B and C be [ ] T x y P P , we obtain The key to designing the guidance law for the wingmen in the formation is to convert the relative distances between the UAVs and the target tracking points, the velocity vector deviations, the heading deviations, and other information into guidance commands. The algorithm for the guidance law can be designed using either geometric methods or control-theory methods. The former use the geometric relationship between UAV and target point to design the guidance law, characterized by convenient analysis and simulation; geometric methods are the ones used most commonly in engineering practice.
The geometric parallel-approach method is used as the guidance law for cooperative formation. The parallel-approach method was designed originally for missile interception and attacking targets, and the resulting trajectories are relatively straight. The specific implementation of the parallel-approach method herein is to keep the line of sight from UAV to target point unchanged during flight, and the line-of-sight angle forms a set of parallel lines in space, as shown in Figure 3.
In the formation of the lead-aircraft-wingman control mode selected herein, the lead aircraft sends real-time status information through an airborne data link, and the wingmen in the formation obtain m v and m  of the lead aircraft through calculation.
The realization of the parallel-approach method in this paper is to keep the line-ofsight angle q between the wingman and the dynamic tracking point unchanged during the flight. In the ground coordinate system, the line-of-sight angle between the wingman In Figure 3, v m is the velocity vector of the UAV tracking the target point whose azimuth angle is η m . The velocity vector V d of the UAV has components v i that are equal to v m , i.e., v i = v m (8) In the formation of the lead-aircraft-wingman control mode selected herein, the lead aircraft sends real-time status information through an airborne data link, and the wingmen in the formation obtain v m and η m of the lead aircraft through calculation.
The realization of the parallel-approach method in this paper is to keep the line-ofsight angle q between the wingman and the dynamic tracking point unchanged during the flight. In the ground coordinate system, the line-of-sight angle between the wingman and the dynamic tracking point is ψ e , and the angle between the velocity vector and the line-of-sight direction is defined as the lead angle, as shown in Figure 4. Then, the flight speed of the wingman is V F , the speed azimuth is ψ F , and the lead angle is η F ; the flight speed of the dynamic tracking point is V L , the speed azimuth is ψ F , and the lead angle is η F . The horizontal straight-line distance between the wingman and the dynamic tracking point is d horiz .

Simulation Verification of Guidance Law
Here, we introduce the method of using software modeling and simulation to simulate the built formation model, which is an effective way to test the algorithm's robustness and provide data support for the subsequent actual flight. The specific method is to model the stand-alone control models of the lead plane and the wingmen under MATLAB/Simulink and perform modular encapsulation processing to facilitate the subsequent expansion of the formation scale. At the same time, the guidance-law part is also modeled to simulate the control of the aircraft by the guidance law and the control law during actual flight. Finally, the aircraft's control, trajectory, and attitude parameters are output for subsequent data analysis and visual trajectory analysis.

Construction of Simulation Model
The guidance-law simulation is mainly to verify the performance of the guidancelaw module. Here, we establish a particle model of a sample UAV. The channel dynamics of the leader and wingman control models are simulated using transfer functions, and their response is the same as the current design control. The law remains the same. Each control channel of the UAV in the particle model can respond to the command signal given by the guidance law in time and at the same time directly output basic information such as altitude, speed, and attitude angle. The output form is convenient for subsequent analysis of the simulated flight trajectory.
The lateral-heading control channel's modeling involves the roll-angle-rate control loop as a first-order inertial link because bank-to-turn is used in control herein to control the heading channel. Therefore, the modeling of the heading channel can be realized based on the roll channel and on the complex number field; the physical relationship is The following geometric relationships exist between the lead angle, velocity azimuth angle, and line-of-sight angle of the wingman and the dynamic tracking point: The relative motion speed . d horiz of the wingman and the dynamic tracking point is The change law of the line-of-sight angle between the wingman and the dynamic tracking point is These two formulas are the guiding equations of the parallel-approach method. The constraints of the equation are: In the formula, ψ e can be calculated from the position of the wingman and the position of the dynamic tracking point. As shown in Figure 2, the coordinate difference between points B and C is [P x P y ] T , and the positional relationship is shown in Figure 4, then the calculation of ψ e is ψ e = arctan P y P x (13) V L sin η L is the projection of the velocity of the dynamic tracking point in the normal direction of the line of sight, v normal , shown in Figure 4, using the calculation method of vector dot product. Let the position difference vector P = [P x P y ] T between the wingman and the virtual dynamic point, the speed vector of the virtual dynamic point is v = [v n v e ] T , then Pv = |P||v| cos η L , where |v| cos η L is the speed of the dynamic virtual point V L on the vector P Tangential projection of v horiz . v horiz = |v| cos η L = (P x v n + P y v e ) |P| (14) Drones 2022, 6, 99 6 of 24 The vector of velocity v horiz is expressed as The normal projection velocity vector of the velocity V L of the dynamic virtual point on the vector → P is denoted as v normal = v − v horiz . The modulus of v normal is recorded as |v normal |, which is the numerical value of |v normal |, and the expected heading control command of the wingman can be obtained as

Simulation Verification of Guidance Law
Here, we introduce the method of using software modeling and simulation to simulate the built formation model, which is an effective way to test the algorithm's robustness and provide data support for the subsequent actual flight. The specific method is to model the stand-alone control models of the lead plane and the wingmen under MATLAB/Simulink and perform modular encapsulation processing to facilitate the subsequent expansion of the formation scale. At the same time, the guidance-law part is also modeled to simulate the control of the aircraft by the guidance law and the control law during actual flight. Finally, the aircraft's control, trajectory, and attitude parameters are output for subsequent data analysis and visual trajectory analysis.

Construction of Simulation Model
The guidance-law simulation is mainly to verify the performance of the guidance-law module. Here, we establish a particle model of a sample UAV. The channel dynamics of the leader and wingman control models are simulated using transfer functions, and their response is the same as the current design control. The law remains the same. Each control channel of the UAV in the particle model can respond to the command signal given by the guidance law in time and at the same time directly output basic information such as altitude, speed, and attitude angle. The output form is convenient for subsequent analysis of the simulated flight trajectory.
The lateral-heading control channel's modeling involves the roll-angle-rate control loop as a first-order inertial link because bank-to-turn is used in control herein to control the heading channel. Therefore, the modeling of the heading channel can be realized based on the roll channel and on the complex number field; the physical relationship is where ψ is the heading angle, φ is the roll angle, v g is the target speed, and g is the acceleration of gravity. As shown in Figure 5, the model input is the roll-angle command φ c of the formation guidance-law module, and the outputs are the UAV model's roll angle and heading angle. The intermediate variable of the control channel is the control of the angular roll rate, and we use a first-order transfer function with a bandwidth of ω n = 6 to simulate the dynamics of the angular-rate channel. where  is the heading angle,  is the roll angle, g v is the target speed, and g is the acceleration of gravity. As shown in Figure 5, the model input is the roll-angle command c  of the formation guidance-law module, and the outputs are the UAV model's roll angle and heading angle. The intermediate variable of the control channel is the control of the angular roll rate, and we use a first-order transfer function with a bandwidth of 6 n   to simulate the dynamics of the angular-rate channel. The modeling of the longitudinal pitch control channel involves regarding the pitchangle-rate control loop as a first-order inertial link, as shown in Figure 6. The input is the pitch-angle command c  of the formation guidance-law module, and the output is the pitch angle of the UAV model. The intermediate variable of the control channel is the pitch rate, and the modeling uses a first-order transfer function with a bandwidth of 7 n   to The modeling of the longitudinal pitch control channel involves regarding the pitchangle-rate control loop as a first-order inertial link, as shown in Figure 6. The input is the pitch-angle command θ c of the formation guidance-law module, and the output is the pitch angle of the UAV model. The intermediate variable of the control channel is the pitch rate, and the modeling uses a first-order transfer function with a bandwidth of ω n = 7 to simulate the dynamics of the angular-rate channel. n late the dynamics of the angular-rate channel. The modeling of the longitudinal pitch control channel involves regarding the pitchangle-rate control loop as a first-order inertial link, as shown in Figure 6. The input is the pitch-angle command c  of the formation guidance-law module, and the output is the pitch angle of the UAV model. The intermediate variable of the control channel is the pitch rate, and the modeling uses a first-order transfer function with a bandwidth of 7 n   to simulate the dynamics of the angular-rate channel. The longitudinal velocity control channel is modeled by considering acceleration as a first-order inertial link, as shown in Figure 7. The input is the target-speed command gc v of the formation guidance-law module, and the output is the indicated airspeed of the UAV model currently flying. The acceleration is modeled using a first-order transfer function with a bandwidth of 3 n   to simulate its dynamics.

Results of Simulation of Guidance Law
The formation guidance law can be simulated by building the leader-follower formation simulation model. The intermediate variables of the flight trajectory and guidance law are stored temporarily in the MATLAB workspace, and various output data can be visualized and analyzed by graph-drawing scripts.
The core index of the guidance law of the parallel-approach method is the tracking effect of the heading channel. Figure 8 shows the dynamic change in the heading angles and velocities of the leader and the three wingmen with time. When the heading of the leader changes rapidly, the heading channels of the three wingmen respond and track it in time. When the lead plane's heading becomes stable, the wingmen follow it smoothly after a small amount of overshoot in time. The longitudinal velocity control channel is modeled by considering acceleration as a first-order inertial link, as shown in Figure 7. The input is the target-speed command v gc of the formation guidance-law module, and the output is the indicated airspeed of the UAV model currently flying. The acceleration is modeled using a first-order transfer function with a bandwidth of ω n = 3 to simulate its dynamics.
roll rate, and we use a first-order transfer function with a bandwidth of 6 n   to simulate the dynamics of the angular-rate channel. The modeling of the longitudinal pitch control channel involves regarding the pitchangle-rate control loop as a first-order inertial link, as shown in Figure 6. The input is the pitch-angle command c  of the formation guidance-law module, and the output is the pitch angle of the UAV model. The intermediate variable of the control channel is the pitch rate, and the modeling uses a first-order transfer function with a bandwidth of 7 n   to simulate the dynamics of the angular-rate channel. The longitudinal velocity control channel is modeled by considering acceleration as a first-order inertial link, as shown in Figure 7. The input is the target-speed command gc v of the formation guidance-law module, and the output is the indicated airspeed of the UAV model currently flying. The acceleration is modeled using a first-order transfer function with a bandwidth of 3 n   to simulate its dynamics.

Results of Simulation of Guidance Law
The formation guidance law can be simulated by building the leader-follower formation simulation model. The intermediate variables of the flight trajectory and guidance law are stored temporarily in the MATLAB workspace, and various output data can be visualized and analyzed by graph-drawing scripts.
The core index of the guidance law of the parallel-approach method is the tracking effect of the heading channel. Figure 8 shows the dynamic change in the heading angles and velocities of the leader and the three wingmen with time. When the heading of the leader changes rapidly, the heading channels of the three wingmen respond and track it in time. When the lead plane's heading becomes stable, the wingmen follow it smoothly after a small amount of overshoot in time.

Results of Simulation of Guidance Law
The formation guidance law can be simulated by building the leader-follower formation simulation model. The intermediate variables of the flight trajectory and guidance law are stored temporarily in the MATLAB workspace, and various output data can be visualized and analyzed by graph-drawing scripts.
The core index of the guidance law of the parallel-approach method is the tracking effect of the heading channel. Figure 8 shows the dynamic change in the heading angles and velocities of the leader and the three wingmen with time. When the heading of the leader changes rapidly, the heading channels of the three wingmen respond and track it in time. When the lead plane's heading becomes stable, the wingmen follow it smoothly after a small amount of overshoot in time. In summary, the guidance law designed herein leads to a good formation and meets the requirements of fast response and stable tracking. Therefore, it can be used in the subsequent semiphysical simulation and actual flight.  In summary, the guidance law designed herein leads to a good formation and meets the requirements of fast response and stable tracking. Therefore, it can be used in the subsequent semiphysical simulation and actual flight.

Formation Management and Control
The flight mission of the leader in the cooperative formation [16][17][18] is described by a set of waypoints known as a waypoint list. After switching to fixed-wing mode, the lead pilot executes a mission route consisting of waypoints, flying to each waypoint in the list in turn. Before the wingmen enter formation mode, there are also waypoints to maintain normal fixed-wing flight. The waypoints in the lead waypoint list are represented using the position vector and those of the wingman is represented by the position vector where d z is the cross-track error of the wingman waypoint relative to the leader in the vertical direction, thereby allowing calculation of the list of waypoints that each wingman flies before entering formation mode. The cooperative formation-management method [19][20][21][22] herein involves the lead plane flying in formation with the wingmen in the link by managing a set of cooperative waypoints. The lead plane coordinates and moves the set of waypoints when formation changes are required. As shown in Figure 9, a group of parallel cooperative path points is generated by the wingmen following the leader, and a straight formation is formed.

Formation Keeping Based on Coordinated Waypoints
The formation maintenance of cooperative formations is studied in the formation coordinate system. As the coordinator, the lead aircraft generates a list of cooperating waypoints by adjusting their positions, and it sends it to the corresponding wingmen in the link through the intermachine communication module. The wingmen then extract the corresponding waypoints from the list of cooperative waypoints assigned by the lead plane according to its own, thereby realizing coordinated formation flight. Based on this method, we design the following three types of formation.
1. Linear formation. In the formation coordinate system, the position coordinates of the lead aircraft are (0, 0), those of the first cooperative path point are (0, xoffset), and those of the other i cooperative path points are (0, xoffset × i). We then construct the list of cooperative waypoints for a straight-line formation. Figure 10 shows the relationship of coordinated waypoints in straight horizontal formation.

Formation Keeping Based on Coordinated Waypoints
The formation maintenance of cooperative formations is studied in the formation coordinate system. As the coordinator, the lead aircraft generates a list of cooperating waypoints by adjusting their positions, and it sends it to the corresponding wingmen in the link through the intermachine communication module. The wingmen then extract the corresponding waypoints from the list of cooperative waypoints assigned by the lead plane according to its own, thereby realizing coordinated formation flight. Based on this method, we design the following three types of formation.

1.
Linear formation. In the formation coordinate system, the position coordinates of the lead aircraft are (0, 0), those of the first cooperative path point are (0, x offset ), and those of the other i cooperative path points are (0, x offset × i). We then construct the list of cooperative waypoints for a straight-line formation. Figure 10 shows the relationship of coordinated waypoints in straight horizontal formation.
according to its own, thereby realizing coordinated formation flight. Based on this method, we design the following three types of formation.
1. Linear formation. In the formation coordinate system, the position coordinates of the lead aircraft are (0, 0), those of the first cooperative path point are (0, xoffset), and those of the other i cooperative path points are (0, xoffset × i). We then construct the list of cooperative waypoints for a straight-line formation. Figure 10 shows the relationship of coordinated waypoints in straight horizontal formation. 2. Triangular formation. In the formation coordinate system, the position coordinates of the lead aircraft are (0, 0), and there are at most (2n−1) cooperative waypoints in n lines, i.e., one line has one cooperative waypoint, and two lines have three cooperative waypoints. The y-axis coordinate interval between every two rows of collaborative path points is yoffset, and the x-axis coordinate interval between every two columns of collaborative path points is xoffset. A list of cooperating waypoints for a triangular formation can be constructed from this. Figure 11 shows the relationships of cooperating waypoints in horizontal triangular formations. 2. Triangular formation. In the formation coordinate system, the position coordinates of the lead aircraft are (0, 0), and there are at most (2n−1) cooperative waypoints in n lines, i.e., one line has one cooperative waypoint, and two lines have three cooperative waypoints. The y-axis coordinate interval between every two rows of collaborative path points is y offset , and the x-axis coordinate interval between every two columns of collaborative path points is x offset . A list of cooperating waypoints for a triangular formation can be constructed from this. Figure 11 shows the relationships of cooperating waypoints in horizontal triangular formations. 3. Stepped formation. A stepped formation is 3D, and in the formation coordinate system, the position coordinates of the leader are (0, 0) and each line has only one collaborative path point. As shown in the plan view of the stepped formation in Figure 12, the y-axis coordinate interval between every two rows of cooperative waypoints in the horizontal direction is yoffset, and the x-axis coordinate interval between every two columns of cooperative waypoints is xoffset. The side view of the stepped formation in Figure 13 shows that in the vertical direction, the z-axis coordinate interval between every two rows of cooperative path points is zoffset, and the z-axis coordinate interval between every two columns of cooperative path points is xoffset. A list of cooperating waypoints for the stepped formation can be constructed from this. 3.
Stepped formation. A stepped formation is 3D, and in the formation coordinate system, the position coordinates of the leader are (0, 0) and each line has only one collaborative path point. As shown in the plan view of the stepped formation in Figure 12, the y-axis coordinate interval between every two rows of cooperative waypoints in the horizontal direction is y offset , and the x-axis coordinate interval between every two columns of cooperative waypoints is x offset . The side view of the stepped formation in Figure 13 shows that in the vertical direction, the z-axis coordinate interval between every two rows of cooperative path points is z offset , and the z-axis coordinate interval between every two columns of cooperative path points is x offset . A list of cooperating waypoints for the stepped formation can be constructed from this.
waypoints in the horizontal direction is yoffset, and the x-axis coordinate interval between every two columns of cooperative waypoints is xoffset. The side view of the stepped formation in Figure 13 shows that in the vertical direction, the z-axis coordinate interval between every two rows of cooperative path points is zoffset, and the z-axis coordinate interval between every two columns of cooperative path points is xoffset. A list of cooperating waypoints for the stepped formation can be constructed from this.  The formations mentioned above are three common ones in which the target tracking position can be determined by moving several waypoints on the spatial geometric position. Each wingman invokes the guidance-law algorithm according to its position and the positional relationship with the reference cooperative waypoint, and it controls the speed and attitude of its UAV to reach the desired position and form a specific formation.

Research on Formation Transformation Method of Cooperative Formation
In many mission scenarios, UAVs must form different formations to adapt to changes in the mission, which requires research on changes of collaborative formations. After creating a coordinated formation, the UAVs will continue to fly in the existing formation until a new mission requirement arises. Herein, a command from the ground-station software initiates a change in formation. For the coordinated formation of fixed-wing UAVs, the under-actuation and control limitations of the fixed-wing make it impossible to form trajectories from all initial positions to the desired formation; the only possible trajectories involve movement that maintains normal airspeed while approaching the desired location. Additionally, the intermediate transition process of cooperative formation change is very important [23][24][25][26], and it is necessary to consider possible collisions in the process of short-range formation change. According to the control strategy of formation keeping, we design a method for changing formation while also preventing collisions.
Before discussing the formation transformation method, we define the standard for forming a close-range formation. The current position of the center of mass of each wingman is i p , which includes the 3D position information, i.e., The location of each collaborative waypoint is The formations mentioned above are three common ones in which the target tracking position can be determined by moving several waypoints on the spatial geometric position. Each wingman invokes the guidance-law algorithm according to its position and the positional relationship with the reference cooperative waypoint, and it controls the speed and attitude of its UAV to reach the desired position and form a specific formation.

Research on Formation Transformation Method of Cooperative Formation
In many mission scenarios, UAVs must form different formations to adapt to changes in the mission, which requires research on changes of collaborative formations. After creating a coordinated formation, the UAVs will continue to fly in the existing formation until a new mission requirement arises. Herein, a command from the ground-station software initiates a change in formation. For the coordinated formation of fixed-wing UAVs, the under-actuation and control limitations of the fixed-wing make it impossible to form trajectories from all initial positions to the desired formation; the only possible trajectories involve movement that maintains normal airspeed while approaching the desired location. Additionally, the intermediate transition process of cooperative formation change is very important [23][24][25][26], and it is necessary to consider possible collisions in the process of short-range formation change. According to the control strategy of formation keeping, we design a method for changing formation while also preventing collisions.
Before discussing the formation transformation method, we define the standard for forming a close-range formation. The current position of the center of mass of each wingman is p i , which includes the 3D position information, i.e., The location of each collaborative waypoint is and the position error distance d err of a wingman and a cooperative waypoint in the horizontal direction determines whether the formation forms a tight formation on the We define the variable k to indicate whether the cooperative formation has formed a tight formation; k = 0 means that a wingman has not reached the cooperative target waypoint in the horizontal direction, whereas k = 1 implies that the wingman has reached the cooperative target waypoint in the horizontal direction and is now partnered with the leader. The value of k is judged by the relay method, with the judgment condition depending on the position error distance d err between the wingman position and the cooperative path point; we have where r max is the distance difference between judging whether the wingman reaches the cooperative path point, which is related to the UAV's wingspan and formation spacing; for the actual flight, we use r max = 2m. The formation transformation method [27][28][29][30] for cooperative formation proposed herein involves performing formation transformation in the horizontal direction first and then in the vertical direction when the wingmen in the horizontal direction reach the constraint condition of close formation. The specific method is when it means that the wingman has not yet reached the cooperative target waypoint in the horizontal direction, and the control of the altitude channel is the height of the wingman waypoint, which means that the formation has completed the close formation in the horizontal direction. At this time, the wingmen set the control target of the altitude channel. To correspond to the heights of the cooperative waypoints, the transformation process is as shown in Figure 14.
As an example to illustrate the process of formation transformation, we take the transformation from a linear formation to a stepped formation. As shown in Figure 15, the formation keeps flying in a straight formation until time t 1 , and it begins to adjust the flight altitude of each wingman to the corresponding waypoint altitude. At time t 2 , the lead plane starts to redistribute the list information of the cooperative waypoints, transforming from a linear formation to a stepped one. At time t 3 , the positions of the cooperative waypoints are reassigned, and the wingmen begin to track the new cooperative waypoints. At time t 4 , the horizontal direction satisfies the constraints for a tight formation. After that, the vertical formation height control is carried out to meet the corresponding vertical spacing requirements of the stepped formation. After this, the formation has completed a transformation. The advantage of this method is that in the formation transformation process, it effectively reduces the probability of collision between the UAVs in short-range formation.
herein involves performing formation transformation in the horizontal direction first and then in the vertical direction when the wingmen in the horizontal direction reach the constraint condition of close formation. The specific method is when it means that the wingman has not yet reached the cooperative target waypoint in the horizontal direction, and the control of the altitude channel is the height of the wingman waypoint, which means that the formation has completed the close formation in the horizontal direction. At this time, the wingmen set the control target of the altitude channel. To correspond to the heights of the cooperative waypoints, the transformation process is as shown in Figure 14. As an example to illustrate the process of formation transformation, we take the transformation from a linear formation to a stepped formation. As shown in Figure 15, the formation keeps flying in a straight formation until time t1, and it begins to adjust the flight altitude of each wingman to the corresponding waypoint altitude. At time t2, the lead plane starts to redistribute the list information of the cooperative waypoints, transforming from a linear formation to a stepped one. At time t3, the positions of the cooperative waypoints are reassigned, and the wingmen begin to track the new cooperative waypoints. At time t4, the horizontal direction satisfies the constraints for a tight formation. After that, the vertical formation height control is carried out to meet the corresponding vertical spacing requirements of the stepped formation. After this, the formation has completed a transformation. The advantage of this method is that in the formation transformation process, it effectively reduces the probability of collision between the UAVs in short-range formation.

Formation Reconfiguration Strategy
The coordinated formation of UAVs involves some complex flight environments; there is a certain probability that some UAVs will suffer structural damage or power failure during the flight mission. Some complex electromagnetic environments may also cause the loss of communication links. These contingencies must be incorporated into the formation control strategy to complete the established formation flight mission.
To this end, we design a formation reconstruction strategy for some typical unexpected situations, which we divide into two types: (i) data-link topology damage and (ii) UAV self-failure. The data links use the mesh networking mode with antidestruction characteristics for this topic. In cases where there is no direct path between the transmitting node and the target node, automatic path searching can be enabled to find and provide path information to ensure that data are delivered to the destination node. By discovering unknown paths automatically and using the best one to transfer the data to the expected destination effectively, the redundancy of the data channel can be maintained when the data-link topology is damaged. As shown in Figure 16, there are multiple redundant channels for the communication between each UAV and the others. Taking the communication connection from vehicle A to vehicle C as an example, there are three paths in the link that can reach vehicle C, i.e., ABC, AFC, and AEDC; therefore, if the communication of a node is damaged, it will not affect the data transmission from vehicle A to vehicle C.

Formation Reconfiguration Strategy
The coordinated formation of UAVs involves some complex flight environments; there is a certain probability that some UAVs will suffer structural damage or power failure during the flight mission. Some complex electromagnetic environments may also cause the loss of communication links. These contingencies must be incorporated into the formation control strategy to complete the established formation flight mission.
To this end, we design a formation reconstruction strategy for some typical unexpected situations, which we divide into two types: (i) data-link topology damage and (ii) UAV selffailure. The data links use the mesh networking mode with antidestruction characteristics for this topic. In cases where there is no direct path between the transmitting node and the target node, automatic path searching can be enabled to find and provide path information to ensure that data are delivered to the destination node. By discovering unknown paths automatically and using the best one to transfer the data to the expected destination effectively, the redundancy of the data channel can be maintained when the data-link topology is damaged. As shown in Figure 16, there are multiple redundant channels for the communication between each UAV and the others. Taking the communication connection from vehicle A to vehicle C as an example, there are three paths in the link that can reach vehicle C, i.e., ABC, AFC, and AEDC; therefore, if the communication of a node is damaged, it will not affect the data transmission from vehicle A to vehicle C. unknown paths automatically and using the best one to transfer the data to the expected destination effectively, the redundancy of the data channel can be maintained when the data-link topology is damaged. As shown in Figure 16, there are multiple redundant channels for the communication between each UAV and the others. Taking the communication connection from vehicle A to vehicle C as an example, there are three paths in the link that can reach vehicle C, i.e., ABC, AFC, and AEDC; therefore, if the communication of a node is damaged, it will not affect the data transmission from vehicle A to vehicle C.  Failure of a UAV itself is divided into (i) lead-aircraft failure and (ii) wingman failure. In the coordinated formation of leader-follower mode, all the wingmen receive the position information of the leader in real-time. When more than λ communication cycles, the wingman still does not receive the leader's data, and the leader is considered to be damaged. The follower with the minimum ID number in the link then acts as the leader; the other wingmen begin to receive the control information from the new leader, who also begins to output the list of cooperative waypoints.
In coordinated formation flight in leader-wingman mode, the wingmen far outnumber the leader, so the probability of a wingman being damaged is higher. When flying in a coordinated formation, the lead plane maintains a status table of whether each wingman is online. This status table is calculated based on the heartbeat packets sent by the wingmen to the lead plane. The specific logic is as follows: if the lead plane continues to receive a wingman's heartbeat packet for time T min , then add this wingman to the online status table; if the lead plane does not receive a wingman's heartbeat packet for time T max , then remove this wingman from the online status table. Here, T min is the shortest time for the lead plane and a wingman to establish a stable connection, and T max is the longest waiting time for the lead plane and a wingman to lose the connection; these two times are set according to the specific indicators of the data links. Therefore, when a wingman loses connection, the leader must remove it from the online status table, and the formation's coordinated waypoint list is also adjusted accordingly. As shown in Figure 17, if vehicle 2 is damaged, then vehicle 3 moves to take its place. cles, the wingman still does not receive the leader's data, and the leader is considered to be damaged. The follower with the minimum ID number in the link then acts as the leader; the other wingmen begin to receive the control information from the new leader, who also begins to output the list of cooperative waypoints. In coordinated formation flight in leader-wingman mode, the wingmen far outnumber the leader, so the probability of a wingman being damaged is higher. When flying in a coordinated formation, the lead plane maintains a status table of whether each wingman is online. This status table is calculated based on the heartbeat packets sent by the wingmen to the lead plane. The specific logic is as follows: if the lead plane continues to receive a wingman's heartbeat packet for time Tmin, then add this wingman to the online status table; if the lead plane does not receive a wingman's heartbeat packet for time Tmax, then remove this wingman from the online status table. Here, Tmin is the shortest time for the lead plane and a wingman to establish a stable connection, and Tmax is the longest waiting time for the lead plane and a wingman to lose the connection; these two times are set according to the specific indicators of the data links. Therefore, when a wingman loses connection, the leader must remove it from the online status table, and the formation's coordinated waypoint list is also adjusted accordingly. As shown in Figure 17, if vehicle 2 is damaged, then vehicle 3 moves to take its place.

System Composition
The semiphysical simulation system consisted of a simulation machine, a flight-control system, ground-station telemetry software, trajectory display software, and data-link communication equipment. The relationship among the different approaches is shown in Figure 18. The nonlinear six-degrees-of-freedom dynamic model of the UAV ran in the simulator, and the flight-control computer received the model data of the simulator for calculation. The command signal was output to the simulator to control the model. The data links were responsible for the communication between the flight-control computers and the ground-station software. The latter was divided into a display area and a command area: the display area displayed status information about the UAV flight, while the

System Composition
The semiphysical simulation system consisted of a simulation machine, a flight-control system, ground-station telemetry software, trajectory display software, and data-link communication equipment. The relationship among the different approaches is shown in Figure 18. The nonlinear six-degrees-of-freedom dynamic model of the UAV ran in the simulator, and the flight-control computer received the model data of the simulator for calculation. The command signal was output to the simulator to control the model. The data links were responsible for the communication between the flight-control computers and the ground-station software. The latter was divided into a display area and a command area: the display area displayed status information about the UAV flight, while the command area contained common remote-control commands to control the UAV formation. Finally, the Tacview trajectory-display software can display the flight trajectories of all the UAVs in the formation, thereby showing the formation more intuitively. The equipment connection relationship of the constructed semiphysical simulation system is shown in Figure 19. The simulation machine used the Speedgoat real-time simulation platform. The ground-station computer was responsible for starting the simulator, running the ground-station software, and establishing communication with the simulator through network protocols. The ground-station computer was connected to a communication node of the data communication link to establish wireless communication with all flight-control computers participating in the simulation; communication between the simulator and the flight-control computers was carried out through an RS422 serial port. The Speedgoat real-time simulation system is the official hardware platform of MATLAB/Simulink, and its operating system and software were developed by Math-Works and integrates MATLAB/Simulink seamlessly. The simulation system used an Intel processor and had powerful computing ability. All the logic was written in Simulink, and The equipment connection relationship of the constructed semiphysical simulation system is shown in Figure 19. The simulation machine used the Speedgoat real-time simulation platform. The ground-station computer was responsible for starting the simulator, running the ground-station software, and establishing communication with the simulator through network protocols. The ground-station computer was connected to a communication node of the data communication link to establish wireless communication with all flight-control computers participating in the simulation; communication between the simulator and the flight-control computers was carried out through an RS422 serial port. The equipment connection relationship of the constructed semiphysical simulation system is shown in Figure 19. The simulation machine used the Speedgoat real-time simulation platform. The ground-station computer was responsible for starting the simulator, running the ground-station software, and establishing communication with the simulator through network protocols. The ground-station computer was connected to a communication node of the data communication link to establish wireless communication with all flight-control computers participating in the simulation; communication between the simulator and the flight-control computers was carried out through an RS422 serial port. The Speedgoat real-time simulation system is the official hardware platform of MATLAB/Simulink, and its operating system and software were developed by Math-Works and integrates MATLAB/Simulink seamlessly. The simulation system used an Intel processor and had powerful computing ability. All the logic was written in Simulink, and there was no need to write much code; the model could also be converted directly into real-time code for execution. This project used the Baseline-S basic-version simulator from The Speedgoat real-time simulation system is the official hardware platform of MAT-LAB/Simulink, and its operating system and software were developed by MathWorks and integrates MATLAB/Simulink seamlessly. The simulation system used an Intel processor and had powerful computing ability. All the logic was written in Simulink, and there was no need to write much code; the model could also be converted directly into real-time code for execution. This project used the Baseline-S basic-version simulator from Speedgoat, which is small and easy to move, weighing only 2.56 kg. It has two RS232 ports and three USB ports, which can be connected to multiple UAVs simultaneously.
The flight-control computer and the data communication links constituted the flightcontrol system. We used four sets of flight-control systems to simulate the formation of four aircraft. Each control system was connected to the Speedgoat real-time simulation system, receiving simulation data and controlling the dynamic model running in the simulator. The data communication links click the four sets of flight-control computers with the ground-station computers for networking communication.

Semiphysical Simulation Process
The flight-control system had to be started first to perform semiphysical simulation. After establishing a connection with the ground-station software, we set the flight-control system to the hardware-in-the-loop simulation mode, we then started the Speedgoat simulator, whereupon the UAV model in the simulator started to run and output data to the flight-control computer. After the ground-station software had checked that the flight-control computer system was in a normal state, each UAV model took off in sequence, and then a formation start command was sent in fixed-wing mode to enter the formation simulation verification.
The semiphysical simulation was a desktop simulation verification link before the actual flight test in the field, and it was necessary to test the scenes used in the flight accordingly. Herein, we report the semiphysical simulation of the following scenarios.

• Case 1: Simulation of formation-mode flight
This simulation was aimed at multiple UAVs flying in cooperative formation, with a two-aircraft formation and a four-aircraft formation used as the research objects to conduct semiphysical simulations. The specific operation was for each UAV to takeoff according to the process and enter the fixed. Start the coordinated formation flight mode after the wings are in the state. Figure 20 shows the simulation results from the perspectives of the ground station (Figure 20a) and the Tacview trajectory-display software (Figure 20b) in cooperative formation mode. As can be seen, each wingman maintained a fixed distance from the leader to fly cooperatively. Speedgoat, which is small and easy to move, weighing only 2.56 kg. It has two RS232 ports and three USB ports, which can be connected to multiple UAVs simultaneously. The flight-control computer and the data communication links constituted the flightcontrol system. We used four sets of flight-control systems to simulate the formation of four aircraft. Each control system was connected to the Speedgoat real-time simulation system, receiving simulation data and controlling the dynamic model running in the simulator. The data communication links click the four sets of flight-control computers with the ground-station computers for networking communication.

Semiphysical Simulation Process
The flight-control system had to be started first to perform semiphysical simulation. After establishing a connection with the ground-station software, we set the flight-control system to the hardware-in-the-loop simulation mode, we then started the Speedgoat simulator, whereupon the UAV model in the simulator started to run and output data to the flight-control computer. After the ground-station software had checked that the flightcontrol computer system was in a normal state, each UAV model took off in sequence, and then a formation start command was sent in fixed-wing mode to enter the formation simulation verification.
The semiphysical simulation was a desktop simulation verification link before the actual flight test in the field, and it was necessary to test the scenes used in the flight accordingly. Herein, we report the semiphysical simulation of the following scenarios.  Case 1: Simulation of formation-mode flight This simulation was aimed at multiple UAVs flying in cooperative formation, with a two-aircraft formation and a four-aircraft formation used as the research objects to conduct semiphysical simulations. The specific operation was for each UAV to takeoff according to the process and enter the fixed. Start the coordinated formation flight mode after the wings are in the state. Figure 20 shows the simulation results from the perspectives of the ground station ( Figure 20a) and the Tacview trajectory-display software (Figure 20b) in cooperative formation mode. As can be seen, each wingman maintained a fixed distance from the leader to fly cooperatively. The simulation of formation transformation verified the proposed method for transforming the cooperative formation. First, the UAVs were operated to form a stable formation, whereupon the ground-station software issued the command for formation trans- The simulation of formation transformation verified the proposed method for transforming the cooperative formation. First, the UAVs were operated to form a stable formation, whereupon the ground-station software issued the command for formation transformation. The transformation of the four-machine formation from the triangular formation to the stepped formation as the research object shows the transformation results in Figure 21. The position changes in the formation transformation simulation are related to the formation before and after the transformation, as shown in Figure 22, where the ordinate is the position difference between wingman and leader during formation flight. During triangular formation flight, the three wingmen were distributed behind the leader: follower 2 was the closest to the leader and was located directly behind it, and followers 1 and 3 were on either side of follower 2. After the formation had changed to the stepped formation, the three followers were distributed evenly on one side of the leader plane in turn and flew at fixed distances; follower 3 was on the outermost side of the formation and farthest from the lead plane. Figure 22 shows that the trajectories of position movement during the formation transformation were relatively stable with no rapid jumps. The formation transformation simulation also verifies the guidance law's response characteristics to changes regarding the guidance target. The most important parameter to change is the heading angle. In this semiphysical simulation, the lead plane and the wingmen kept the same heading and flew parallel in the triangular formation. The formation management module generated a new path tracking point after receiving the formation change command, and the wingman guidance law guided the UAVs to the unique The position changes in the formation transformation simulation are related to the formation before and after the transformation, as shown in Figure 22, where the ordinate is the position difference between wingman and leader during formation flight. During triangular formation flight, the three wingmen were distributed behind the leader: follower 2 was the closest to the leader and was located directly behind it, and followers 1 and 3 were on either side of follower 2. After the formation had changed to the stepped formation, the three followers were distributed evenly on one side of the leader plane in turn and flew at fixed distances; follower 3 was on the outermost side of the formation and farthest from the lead plane. Figure 22 shows that the trajectories of position movement during the formation transformation were relatively stable with no rapid jumps. The position changes in the formation transformation simulation are related to the formation before and after the transformation, as shown in Figure 22, where the ordinate is the position difference between wingman and leader during formation flight. During triangular formation flight, the three wingmen were distributed behind the leader: follower 2 was the closest to the leader and was located directly behind it, and followers 1 and 3 were on either side of follower 2. After the formation had changed to the stepped formation, the three followers were distributed evenly on one side of the leader plane in turn and flew at fixed distances; follower 3 was on the outermost side of the formation and farthest from the lead plane. Figure 22 shows that the trajectories of position movement during the formation transformation were relatively stable with no rapid jumps. The formation transformation simulation also verifies the guidance law's response characteristics to changes regarding the guidance target. The most important parameter to change is the heading angle. In this semiphysical simulation, the lead plane and the wingmen kept the same heading and flew parallel in the triangular formation. The formation management module generated a new path tracking point after receiving the formation change command, and the wingman guidance law guided the UAVs to the unique The formation transformation simulation also verifies the guidance law's response characteristics to changes regarding the guidance target. The most important parameter to change is the heading angle. In this semiphysical simulation, the lead plane and the wingmen kept the same heading and flew parallel in the triangular formation. The formation management module generated a new path tracking point after receiving the formation change command, and the wingman guidance law guided the UAVs to the unique path tracking point, which mainly produced changes in the heading angle. Having completed the formation change, the lead plane and the wingmen flew in the new stepped formation and kept the same heading while flying in parallel. Figure 23 shows the evolution of heading angle in the formation change. As can be seen, the wingman heading angle changed rapidly after the command to change the formation, and the wingman heading converged quickly upon completion of the formation change.
Drones 2022, 6, x FOR PEER REVIEW 18 of 26 path tracking point, which mainly produced changes in the heading angle. Having completed the formation change, the lead plane and the wingmen flew in the new stepped formation and kept the same heading while flying in parallel. Figure 23 shows the evolution of heading angle in the formation change. As can be seen, the wingman heading angle changed rapidly after the command to change the formation, and the wingman heading converged quickly upon completion of the formation change.

Flight Platform and Ground Measurement and Control System
The UAVs used in the test are electric vertical takeoff and landing aircraft (eVTOL). The aircraft in this layout can takeoff and land vertically with the rotor, thereby extending the total flight time, which can fly in the air for a long time to verify the control method proposed in this paper. As shown in Figure 24, it adopts the layout of large wingspan and high flat tail. The wingspan is 2.2 m and the takeoff weight is 6 kg. The main parameters of the UAV are shown in Table 1.

Flight Platform and Ground Measurement and Control System
The UAVs used in the test are electric vertical takeoff and landing aircraft (eVTOL). The aircraft in this layout can takeoff and land vertically with the rotor, thereby extending the total flight time, which can fly in the air for a long time to verify the control method proposed in this paper. As shown in Figure 24, it adopts the layout of large wingspan and high flat tail. The wingspan is 2.2 m and the takeoff weight is 6 kg. The main parameters of the UAV are shown in Table 1.
Drones 2022, 6, x FOR PEER REVIEW 18 of 26 path tracking point, which mainly produced changes in the heading angle. Having completed the formation change, the lead plane and the wingmen flew in the new stepped formation and kept the same heading while flying in parallel. Figure 23 shows the evolution of heading angle in the formation change. As can be seen, the wingman heading angle changed rapidly after the command to change the formation, and the wingman heading converged quickly upon completion of the formation change.

Flight Platform and Ground Measurement and Control System
The UAVs used in the test are electric vertical takeoff and landing aircraft (eVTOL). The aircraft in this layout can takeoff and land vertically with the rotor, thereby extending the total flight time, which can fly in the air for a long time to verify the control method proposed in this paper. As shown in Figure 24, it adopts the layout of large wingspan and high flat tail. The wingspan is 2.2 m and the takeoff weight is 6 kg. The main parameters of the UAV are shown in Table 1.    As the best way to show the practicality of the designed control method, the actual flight test in the field was the final verification link of the fixed-wing UAV cooperative formation control method intended herein. The main equipment involved in the real flight test was the UAV flight platform and the ground measurement and control system, described separately below.
The UAV flight platform comprised two categories of equipment: (i) the flight-control system comprising the flight-control computer and the data-link radio and (ii) the actuating equipment, such as motors, electric speed controllers, and steering gears. The complete equipment of the flight platform is shown in Figure 25.

Flight Platform and Ground Measurement and Control System
As the best way to show the practicality of the designed control method, the actual flight test in the field was the final verification link of the fixed-wing UAV cooperative formation control method intended herein. The main equipment involved in the real flight test was the UAV flight platform and the ground measurement and control system, described separately below.
The UAV flight platform comprised two categories of equipment: (i) the flight-control system comprising the flight-control computer and the data-link radio and (ii) the actuating equipment, such as motors, electric speed controllers, and steering gears. The complete equipment of the flight platform is shown in Figure 25. The main parts of the ground measurement and control system were the groundstation computer, the real-time kinematic (RTK) base station, and the data-link desktop terminal, as shown in Figure 26. Among these, the ground-station computer was the main equipment for remote control and telemetry between the ground station and the UAVs; it was responsible for checking the real-time status of the UAVs and sending control instructions. The RTK base station generated differential positioning data; coordinated formation flight requires high positioning accuracy. The differential data generated by the base station were transmitted to the flight-control computer through the data link to enable the drones to enter positioning mode with centimeter-level accuracy. The desktop end of the data link was a node of the data link, and this node was connected to the ground-station computer, securing the ground-station software to the entire communication link. The main parts of the ground measurement and control system were the groundstation computer, the real-time kinematic (RTK) base station, and the data-link desktop terminal, as shown in Figure 26. Among these, the ground-station computer was the main equipment for remote control and telemetry between the ground station and the UAVs; it was responsible for checking the real-time status of the UAVs and sending control instructions. The RTK base station generated differential positioning data; coordinated formation flight requires high positioning accuracy. The differential data generated by the base station were transmitted to the flight-control computer through the data link to enable the drones to enter positioning mode with centimeter-level accuracy. The desktop end of the data link was a node of the data link, and this node was connected to the ground-station computer, securing the ground-station software to the entire communication link. The UAV flight platform and the ground measurement and control system were the core equipment of the formation flight test. Figure 27 shows the equipment connection for the four-machine formation flight test, where a dotted line is a wireless connection between data links. The ground-station system monitored the status of each UAV flight platform in real-time during flight and also initiated some control commands, such as open- The UAV flight platform and the ground measurement and control system were the core equipment of the formation flight test. Figure 27 shows the equipment connection for the four-machine formation flight test, where a dotted line is a wireless connection between data links. The ground-station system monitored the status of each UAV flight platform in real-time during flight and also initiated some control commands, such as opening, changing, and closing the formation. The UAV flight platform and the ground measurement and control system were the core equipment of the formation flight test. Figure 27 shows the equipment connection for the four-machine formation flight test, where a dotted line is a wireless connection between data links. The ground-station system monitored the status of each UAV flight platform in real-time during flight and also initiated some control commands, such as opening, changing, and closing the formation.

Formation Flight Test
The outdoor flight test reported herein involved a four-aircraft formation. Figure 28 shows a photograph of the UAVs before takeoff.

Formation Flight Test
The outdoor flight test reported herein involved a four-aircraft formation. Figure 28 shows a photograph of the UAVs before takeoff. First, each UAV was placed at a fixed takeoff point and a power-on operation was performed. After the flight-control computer system was initialized, the ground station began to detect the states of the UAVs. If any of the drones did not meet the takeoff status, the status of all the drones was rechecked until the takeoff conditions were met. Send the takeoff command, then all drones enter the fixed-wing mode, when all drones are flying to the right place. The station sends the "open formation" command, and the drones enter the cooperative formation flight mode. At this time, the command was sent to change the cooperative formation, and the flight modes of the wingmen in the tight formation were changed according to the flight requirements. After the formation test was completed, the ground station sent a "close formation" command, whereupon the leader and wingmen exited the coordinated formation mode, flew along their respective routes, and finally returned to their takeoff points for recovery ( Figure 29). First, each UAV was placed at a fixed takeoff point and a power-on operation was performed. After the flight-control computer system was initialized, the ground station began to detect the states of the UAVs. If any of the drones did not meet the takeoff status, the status of all the drones was rechecked until the takeoff conditions were met. Send the takeoff command, then all drones enter the fixed-wing mode, when all drones are flying to the right place. The station sends the "open formation" command, and the drones enter the cooperative formation flight mode. At this time, the command was sent to change the cooperative formation, and the flight modes of the wingmen in the tight formation were changed according to the flight requirements. After the formation test was completed, the ground station sent a "close formation" command, whereupon the leader and wingmen exited the coordinated formation mode, flew along their respective routes, and finally returned to their takeoff points for recovery ( Figure 29). takeoff command, then all drones enter the fixed-wing mode, when all drones are flying to the right place. The station sends the "open formation" command, and the drones enter the cooperative formation flight mode. At this time, the command was sent to change the cooperative formation, and the flight modes of the wingmen in the tight formation were changed according to the flight requirements. After the formation test was completed, the ground station sent a "close formation" command, whereupon the leader and wingmen exited the coordinated formation mode, flew along their respective routes, and finally returned to their takeoff points for recovery ( Figure 29). The ground-station software was connected to the communication links of the UAVs through the ground data-link node, and it was responsible for real-time communication with all the UAVs in the formation to (i) monitor their status during the flight, (ii) send the takeoff and landing instructions to each aircraft, and (iii) send the mission instructions, among other duties. The flight-control interface of the ground-station software is shown in Figure 30. The ground-station software was connected to the communication links of the UAVs through the ground data-link node, and it was responsible for real-time communication with all the UAVs in the formation to (i) monitor their status during the flight, (ii) send the takeoff and landing instructions to each aircraft, and (iii) send the mission instructions, among other duties. The flight-control interface of the ground-station software is shown in Figure 30. In the figure, the corresponding commands can be sent to the drone through the control panel. The status of the current drone can be seen in the attitude instrument and flight parameter information in the upper-right corner.
The command panel of the ground-station software was the main functional area for controlling the flight of the UAVs, including the command areas related to (i) single-aircraft takeoff and landing and (ii) formation. The instructions for single-aircraft takeoff and landing are to complete the takeoff and landing of a single UAV, including unlocking, locking, takeoff, and landing. The formation commands control the opening and closing of the formation mode and set the formation parameters after the UAVs entered the fixedwing flight state; they included instructions such as turning the formation mode on and off and setting the formation shape and spacing parameters.
In this test, two types of coordinated formation flight were carried out, i.e., the stepped formation and the straight-line formation. The two formation arrangements were relatively In the figure, the corresponding commands can be sent to the drone through the control panel. The status of the current drone can be seen in the attitude instrument and flight parameter information in the upper-right corner.
The command panel of the ground-station software was the main functional area for controlling the flight of the UAVs, including the command areas related to (i) single-aircraft takeoff and landing and (ii) formation. The instructions for single-aircraft takeoff and landing are to complete the takeoff and landing of a single UAV, including unlocking, locking, takeoff, and landing. The formation commands control the opening and closing of the formation mode and set the formation parameters after the UAVs entered the fixedwing flight state; they included instructions such as turning the formation mode on and off and setting the formation shape and spacing parameters.
In this test, two types of coordinated formation flight were carried out, i.e., the stepped formation and the straight-line formation. The two formation arrangements were relatively simple but could also verify the control method of formation maintenance and change proposed herein. This outdoor flight began with a coordinated flight in a straight formation, and then the ground-station software issued the "set stepped formation" command, and the formation changed to a stepped formation and continued to fly. Figure 31 shows the actual outdoor flight conditions of the linear collaborative formation (Figure 31a) and the stepped collaborative formation (Figure 31b).
In the figure, the corresponding commands can be sent to the drone through the control panel. The status of the current drone can be seen in the attitude instrument and flight parameter information in the upper-right corner.
The command panel of the ground-station software was the main functional area for controlling the flight of the UAVs, including the command areas related to (i) single-aircraft takeoff and landing and (ii) formation. The instructions for single-aircraft takeoff and landing are to complete the takeoff and landing of a single UAV, including unlocking, locking, takeoff, and landing. The formation commands control the opening and closing of the formation mode and set the formation parameters after the UAVs entered the fixedwing flight state; they included instructions such as turning the formation mode on and off and setting the formation shape and spacing parameters.
In this test, two types of coordinated formation flight were carried out, i.e., the stepped formation and the straight-line formation. The two formation arrangements were relatively simple but could also verify the control method of formation maintenance and change proposed herein. This outdoor flight began with a coordinated flight in a straight formation, and then the ground-station software issued the "set stepped formation" command, and the formation changed to a stepped formation and continued to fly. Figure 31 shows the actual outdoor flight conditions of the linear collaborative formation (Figure 31a) and the stepped collaborative formation (Figure 31b).  When it was necessary for the UAVs to return, the ground station sent the command to dissolve the formation, and all the UAVs participating in the formation exited the formation mode and automatically returned to the mission route of single-machine flight ( Figure 32). Finally, the ground station sent the software command to return to the takeoff points and complete the return landing and recovery one by one. When it was necessary for the UAVs to return, the ground station sent the command to dissolve the formation, and all the UAVs participating in the formation exited the formation mode and automatically returned to the mission route of single-machine flight ( Figure 32). Finally, the ground station sent the software command to return to the takeoff points and complete the return landing and recovery one by one.

Flight Data Analysis
Analyzing the flight outcome provides direct and effective data to support the results of the outdoor flight. After the outdoor flight test was completed, the track file recorded in the flight-control computer was imported for data analysis. The data were processed in different drawing and visualization software to analyze the flight results from other dimensions.
The direct indicator of the guidance-law effect in the data is the wingman heading angle's tracking development after the guidance target change. Figure 33 shows the change curve of the heading angle of each UAV in the formation for times between 634 and 654 s in the track record file. As can be seen, when the heading angle of the leader changed, those of the wingmen responded and tracked it in time. Therefore, the response characteristics of the guidance law used herein are verified, and it is shown that the guidance-law algorithm can meet practical needs.

Flight Data Analysis
Analyzing the flight outcome provides direct and effective data to support the results of the outdoor flight. After the outdoor flight test was completed, the track file recorded in the flight-control computer was imported for data analysis. The data were processed in different drawing and visualization software to analyze the flight results from other dimensions.
The direct indicator of the guidance-law effect in the data is the wingman heading angle's tracking development after the guidance target change. Figure 33 shows the change curve of the heading angle of each UAV in the formation for times between 634 and 654 s in the track record file. As can be seen, when the heading angle of the leader changed, those of the wingmen responded and tracked it in time. Therefore, the response characteristics of the guidance law used herein are verified, and it is shown that the guidance-law algorithm can meet practical needs.
Analyzing the flight outcome provides direct and effective data to support the results of the outdoor flight. After the outdoor flight test was completed, the track file recorded in the flight-control computer was imported for data analysis. The data were processed in different drawing and visualization software to analyze the flight results from other dimensions.
The direct indicator of the guidance-law effect in the data is the wingman heading angle's tracking development after the guidance target change. Figure 33 shows the change curve of the heading angle of each UAV in the formation for times between 634 and 654 s in the track record file. As can be seen, when the heading angle of the leader changed, those of the wingmen responded and tracked it in time. Therefore, the response characteristics of the guidance law used herein are verified, and it is shown that the guidance-law algorithm can meet practical needs.  Furthermore, the flight trajectory of each UAV in the formation reflects intuitively the effect of cooperative formation from the data. Figure 34a shows the 3D trajectories of the actual flight route as plotted from the data for a stable formation section. The trajectories comprise the three geographical coordinates of longitude, latitude, and altitude; the horizontal coordinates are longitude and latitude, and the vertical coordinate is altitude. As can be seen, each UAV kept a fixed distance within the allowable error range, and the UAVs flew together. Figure 34b shows the entire formation trajectory, which intuitively reflects the flying effect of the trajectory in the real scene. Furthermore, the flight trajectory of each UAV in the formation reflects intuitively the effect of cooperative formation from the data. Figure 34a shows the 3D trajectories of the actual flight route as plotted from the data for a stable formation section. The trajectories comprise the three geographical coordinates of longitude, latitude, and altitude; the horizontal coordinates are longitude and latitude, and the vertical coordinate is altitude. As can be seen, each UAV kept a fixed distance within the allowable error range, and the UAVs flew together. Figure 34b shows the entire formation trajectory, which intuitively reflects the flying effect of the trajectory in the real scene. Within a certain allowable range of control error, this paper reports research on a cooperative formation control method for fixed-wing UAVs. The effect of the actual flight of the formation meets the design requirements and the requirements of formation control for some specific tasks.

Conclusions
The management of UAVs entering and departing from formation is an indispensable part of intelligent and information-based UAV operations. Aimed at the problem of maintaining a multi-UAV formation, reported herein was the design of a collaborative guidance law based on the classic missile-type parallel-approach method. According to (i) the flight mission characteristics of a lead aircraft following a route and wingmen following the lead aircraft and (ii) the general characteristics of the route, a guidance law Within a certain allowable range of control error, this paper reports research on a cooperative formation control method for fixed-wing UAVs. The effect of the actual flight of the formation meets the design requirements and the requirements of formation control for some specific tasks.

Conclusions
The management of UAVs entering and departing from formation is an indispensable part of intelligent and information-based UAV operations. Aimed at the problem of maintaining a multi-UAV formation, reported herein was the design of a collaborative guidance law based on the classic missile-type parallel-approach method. According to (i) the flight mission characteristics of a lead aircraft following a route and wingmen following the lead aircraft and (ii) the general characteristics of the route, a guidance law was designed for three guidance modes: straight line, turning, and circling. The design of the entry and departure management of the formation led to strategies for the access and departure management of multi-machine independent assembly, dissolution, and special cases, thereby allowing effective management and control of the formation. An entry and departure adjustment strategy for special circumstances was given. The hardwarein-the-loop simulation and measured flight verification showed that the guidance law is practical in meeting UAV formation flight-control requirements for specific flight missions. Although it was tested on low-speed UAV platforms, it is still applicable to high-speed UAV platforms. Future work will focus the algorithm on high-speed platform validation. Data Availability Statement: Data available on request due to restrictions eg privacy or ethical. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to commercial use.

Conflicts of Interest:
The authors declare no conflict of interest.