Real-time nonlinear parameter estimation and tracking control of unmanned aerial vehicles in closed-loop

The real-time unknown parameter estimation and adaptive tracking control problems are investigated in this paper for a six degrees of freedom (6-DOF) of under-actuated quadrotor unmanned aerial vehicle (UAV). A virtual proportional derivative (PD) controller is designed to maintain the translational dynamics. Two adaptive schemes are developed to handle the attitude dynamics of the UAV with several unknown parameters. In the beginning, a classical adaptive scheme (CAS) using the certainty equivalence principle is proposed and designed. The idea is to design a controller for an ideal situation by assuming the unknown parameters were known. Then the unknown parameters are replaced by their estimation. A theoretical analysis is provided to ensure the trajectory tracking of the adaptive controller. However, an inherent drawback of this scheme is that there is no guarantee for the estimated parameters to converge to the actual values. To address this issue, a new adaptive scheme (NAS) is developed as the next step by adding a continuously differentiable function to the control structure. The proposed technique guarantees handling of the parametric uncertainties with an appropriate design manifold. A rigorous analytical proof, numerical simulation analyses, and experimental validation are presented to show the effectiveness of the proposed control design.

www.nature.com/scientificreports/ tracking problem, especially for the attitude dynamics of UAVs. One of the common approaches is the feedback linearization method as developed in 6,7 . However, the parameters of the UAV are not always available or will change over time for the feedback control design in various practical situations. These unknown parameters may cause more complicated technical challenges in designing a proper controller. In general, there are two major research lines to tackle this issue. The first research direction is dealing with robust controllers. The idea behind this approach is to propose a feedback controller to dominate the uncertainties in the system dynamics. In this way, the controller guarantees to handle the uncertainties within a particular bound. This disadvantage of this approach is that the bound of the uncertainties should be available a priori to ensure the stability of the closed-loop control system. One of the popular methods in this research line is the sliding mode control (SMC). SMC is widely implemented in various practical situations as it is less sensitive to parametric uncertainties and disturbances as presented in [8][9][10] . However, chattering is a common issue in designing SMC for the autonomous system. Several studies have been investigated to reduce the chattering as well as to compensate the effect of unmodelled dynamics as presented in [11][12][13] .
The second major research line is adaptive control method. This approach is useful to handle the unknown constant parameters in the system dynamics. The traditional approach or CAS was initially proposed in the literature using the certainty equivalence principle. The idea behind this approach is to cancel the nonlinear terms containing unknown constant parameters. This scheme suggests a two-step control design procedure to handle the uncertainties. A controller under an ideal condition is designed in the first step, where all parameters are assumed to be known for feedback controller. In the second step, every unknown constant parameter is replaced by its estimation generated by an adaptive law along the gradient of Lyapunov function. The perfect cancellation of the nonlinear term is deemed to be achieved by driving the estimated parameter to converge to the actual value. This technique has a simple structure for practical implementation; however, it contains an inherent drawback. In this technique, there is no guarantee that the parameter estimation error converges to zero as it fully relies on the system dynamic states. The estimated parameters are updated constantly (even though they reach their actual values) as long as the states of system dynamics don't reach the equilibrium point. Nonetheless, the parameter estimation adaptive law will stop the update as soon as the system dynamics states reach the equilibrium point. For some more information on this technique, an interested reader may refer to in 5,14,15 .
Another adaptive technique to deal with uncertain parameters is the model reference adaptive control (MRAC). This technique also works using the certainty equivalence principle. A reference model is added to the control structure to estimate the unknown constant parameters as presented in 16,17 . However, this method has a major drawback to guarantee the stability of the closed-loop system as investigated in 18 . To rectify this problem, L 1 adaptive control scheme was proposed by adding a linear filter in the whole control structure 19 . This technique was studied further in 20 for the collaborative UAVs. Note that this scheme doesn't guarantee the estimate parameters to converge to the actual value of the unknown parameters as it was tuned for the worstcase condition. Application of this scheme in UAVs can be found in 21 . The mismatch between the estimate parameters and the actual values are used to update the parameters of the system. As a result, this approach can be implemented mainly in limited cases.
Model predictive control (MPC) is one of the common methods proposed in the literature to handle the uncertainties. For example, adaptive MPC with extended state observer was studied in 22 for UAVs in a networked setting. The unknown disturbances and system uncertainties were handled successfully. However, the controller was not designed to estimate the unknown parameters in the class of the nonlinear function. The difficulties in estimating the unknown parameters of a nonlinear function in closed-loop has been avoided by linearizing the system dynamic in this method. It is not surprising as this is an essential assumption in designing a controller using the linear MPC technique. Another result in which MPC is used to handle the disturbances and sensor faults is reported in 23 .
Another common technique used to handle the uncertainties is the intelligent computation approach. For example, neural network (NN) 24,25 was designed for autonomous agents in a networked setting and genetic algorithm (GA) 26 is proposed to address the uncertainties in a robotic manipulator. However, this approach has two major issues. The first issue is that both the tracking error and the parameter estimation error are not converging to zero asymptotically and will have a residual error. This is caused by the mismatch between the approximator and the unknown nonlinear function. The second issue is that the intelligent computation requires a high performance embedded computer in many cases to estimate the uncertainties. As a result, this approach can only be implemented in limited practical scenarios.
Immersion and invariance (I &I) was developed to address the open problem in CAS as studied in 27,28 . This technique also has a two-step control design. In the first step, the system was designed to satisfy the input to state stability (ISS) condition, where the nonlinear term with unknown parameters are assumed to be zero. Then the unknown nonlinear terms are taken into account in the second step. The adaptive term is designed by replacing the unknown parameter with its estimation and an additional appropriate continuous function design to assist the parameter estimation algorithm to converge in the right direction. In this situation, stability can be guaranteed by driving the mismatch estimation error to a specific manifold. Note that the stability condition is derived by looking at the error dynamics for the parameter estimation rather than the closed-loop system dynamics. Hence, it is reasonable to apply the proposed controller under the ISS condition.
Application of I &I adaptive to handle the unknown parameters in UAVs can be found in 21,29 . However, the I &I adaptive law has the major drawback of requiring the mismatch error between the actual value and the estimated parameters. As a result, the technique contains a significant bottleneck to be implemented in practical situations. To solve this issue, we develop two adaptive schemes for a quadrotor UAV with unknown constant parameters in this study. A traditional or CAS is designed in the first scheme to handle the unknown parameters using the certainty equivalence principle. This scheme suggests two steps to estimate the unknown parameters of the system dynamic. First, the controller is designed for the system under the ideal condition by assuming all the www.nature.com/scientificreports/ system parameters to be available for the feedback control design. Then the unknown parameters are replaced by their estimated values updated by adaptation law in the second step. As mentioned above, this classical scheme contains an inherent drawback to estimate the unknown parameters. Inspired by the studies in 27,30,31 , we develop NAS in the second scheme to handle multiple uncertain parameters in the rotational dynamics of the UAV. An appropriate manifold is designed to eliminate the effect of unknown parameters, where it is not necessary to drive the steady-state parametric estimation error to zero, but it follows a pre-defined manifold of UAV states. Here the ISS condition is relaxed compared to the result in 27 . The stability analysis is presented on both the closed-loop trajectory tracking and parameter estimation error dynamics. The mismatch between the actual values and estimate parameters is not required for both feedback control design and adaptation law. In this way, we can remove the inherent drawback in CAS. This extension is a very important extension for the current approach to deal with unknown constant parameters in practical implementation. A virtual PD controller is designed to track the desired position in the outer loop of the proposed nested control structure. The major contributions of this study are summarized as follows: • The asymptotic stablisation and trajectory tracking problem for UAVs as a fully coupled and nonlinear underactuated system with unknown parameters is studied using adaptive control schemes. • A new adaptive control scheme (NAS) for a UAV with unknown parameters is proposed to estimate the unknown parameters of the UAV in closed-loop so that it guarantees the estimated parameters of the UAV converge to a known value from which the value of true parameters can be derived. This is an extension compared to the classical scheme (CAS) in which there is no guarantee that the parameter estimation error converges to zero. • Through extensive numerical simulations the superiority of NAS are evaluated and compared with various existing methods and implemented on a real UAV in an experimental setting. • Some simulations and practical implementations for mini-drone are conducted to verify the performance of our scheme. Especially in the practical test, we demonstrate that our new scheme has a simple structure that makes it possible to be tested in a mini-drone with limited hardware computation.
The remainder of this paper is organized as follows. A dynamic model of UAV is presented in "System dynamics of UAV" section. Following that, the position tracking control for the outer loop and the two proposed adaptive schemes for the attitude dynamics are presented in "Proposed control design" section. Then in "Simulation and experimental results" section, with extensive numerical analyses and simulation results the efficacy of the proposed controllers are demonstrated. Finally, the paper is concluded in "Conclusion and directions for future work" section with a conclusion and some suggestions for future work.

System dynamics of UAV
The 6-DOF dynamic model of the quadrotor UAV is represented by the following states where η 1 = x y z T ∈ R 3×1 is the absolute position of UAV composed of forward (x), lateral (y) and vertical (z) motions; η 2 = φ θ ψ T ∈ R 3×1 is the attitude composed by three Euler angles i.e. roll ( φ ), pitch ( θ ) and yaw ( ψ ) motions; ν 1 = u v w T ∈ R 3×1 is a linear velocity; and ν 2 = p q r T ∈ R 3×1 is an angular velocity.
The dynamic of the UAV can be derived as a result of the coupling between the inertial and body frames. This coupling can be expressed by the transformation matrices J 1 (η 2 ) and J 2 (η 2 ) Both φ and θ are assumed between − π 2 to π 2 . As results, cos φ and cos θ are non-zero, hence J . The dynamic coupling between position and orientation vectors can be represented by the following transformation The use of both body fixed frame and the inertial (earth) coordinates are required to study 6-DOF motions of UAV. Both position and orientation are described with respect to the inertial frame. In another side, both linear and angular velocities are described with respect to the body frame. The illustration of these two coordinate frames can be seen in Fig. 1.
As derived in 21 , the translational dynamic of the quadrotor UAV can be represented as (2) www.nature.com/scientificreports/ where g ∈ R , u ∈ R , m ∈ R , and k t ∈ R and are gravity acceleration, thrust force, mass, and translational drag coefficient, respectively. Vector z e = 0 0 1 T ∈ R 3×1 is the unitary vector in z direction. The attitude dynamics is expressed by where I R ∈ R and k r ∈ R are propeller inertia and rotational drag coefficient, respectively. Vector τ = τ p τ q τ r T ∈ R 3×1 shows the torques acting on the body frame in roll, pitch, and yaw directions. The matrix I M = diag I x I y I z ∈ R 3×3 is a diagonal inertia matrix. The relative angular speed ∈ R generated by the motors is The thrust force u, generated by each propeller can be calculated as where f ι ∈ R is the upward-lifting force generated by every rotor and k ι ∈ R is a positive constant. From (2) and (3), it can be seen that the number of control inputs is less than the number of degrees of freedom, and hence the UAV is a nonlinear under-actuated system. The UAV considered in this paper has a cross configuration. The propeller forces and the input control signals acting on the UAV body are related according to (6) where l ∈ R is the arm length and d ∈ R is the drag factor.
The attitude dynamic in (3) with an additional external disturbance can be rewritten in the linearly parameterized form as where In (7), I R and k r are unknown parameters to be estimated for the feedback control design. Also, hδ is an external disturbance acting on the body frame, in which h is a known vector as a function of time and δ is an unknown scalar. The full dynamic equations of the UAV can be represented by looking at Eqs. (1), (2) and (7) together. As a nonlinear function of the UAV states is present in all equations and the translational and rotational motions are coupled together with these three equations, the nonlinear fully coupled nature of the dynamic equations is exploited for the control design.

Proposed control design
As discussed in the previous section, the quadrotor UAV is an under-actuated system, in which the six system states should be controlled by four control inputs. There is also a strong coupling between the translational and attitude dynamic of UAV. In this section, a nested control strategy for trajectory tracking of the quadrotor www.nature.com/scientificreports/ dynamic is proposed. The position tracking error is controlled by a PD controller in the outer loop. However, two adaptive control schemes are proposed to control the attitude dynamic in the presence of uncertain parameters of the system. As the first step, CAS is designed using the certainty equivalence principle. In this technique, the adaptive law is designed to force the parameter estimation error to converge to zero. Nevertheless, in the second technique, NAS is proposed by adding a certain continuous differentiable function in the structure of the adaptive control design. Contrary to CAS, here the estimation error is converging to a manifold rather than zero.
Translational control design. The tracking controller for the translational dynamic is designed by defining the tracking error of the system as where η 1 and η 1 d are the error vector position and the desired vector position, respectively. The double integrator dynamics of (8) can be written as where K P ∈ R 3×3 and K D ∈ R 3×3 are the diagonal matrices denoting the gains for the proportional and derivative terms, respectively. By selecting K P and K D to be positive definite matrices, the system dynamics (9) satisfies the Routh-Hurwitz stability criterion by having lim t→∞η1 (t) = 0 . The dynamics in (8) can be rewritten as Inspired by 32 ,the virtual control input U =η 1 = U 1 U 2 U 3 T is defined and substituted in (2). This yields or By expanding (12), it would be straightforward to verify that Assuming that cos θ = 0 , the pitch rotation can be generated from (13) as Squaring both sides of (12) yields

As a result
From (14) and (15), we have www.nature.com/scientificreports/ and by substituting (18)- (19), the roll rotation can be derived as From (15), the total thrust can be computed as follows The desired roll φ d and pitch θ d angles can be calculated from (20) and (16) is a similar way Attitude control design. In this section, the main contribution of the paper is explained. The presence of uncertain parameters in the attitude dynamics will make the attitude control of UAV a challenging problem. If all parameters of the system are known for the feedback control design, the nonlinear terms could be cancelled by introducing a simple feedback linearization technique. However, in many practical situations, several parameters of the UAV such as I R , k r and δ are often unknown. Consequently, a proper controller such as adaptive techniques is required to handle the uncertainties. CAS, relying on the certainty equivalence principle, has two major steps in the design procedure. First, a controller is designed for the system without any uncertain parameters. In this way, uncertain parameters are assumed to be readily available for the feedback control design. In the second step, the parameters estimated by the adaptive law are used as the values of the uncertain parameters. Before presenting the main results, the desired trajectory is defined as ν 2 d = p d q d r d T ∈ R 3×1 . Therefore, the trajectory error can be calculated from e = ν 2 − ν 2 d . The reference model is introduced by where α ∈ R is a positive constant and ν 2 is the state of the reference model. This dynamics leads to the desired trajectory ν 2 d . Hence, the tracking error dynamics (7) and (24) can be written as A classical adaptive scheme. In the first scheme, we present CAS to estimate the uncertain parameters in the attitude dynamics. Let Î R (t) ∈ R , k r (t) ∈ R and δ (t) ∈ R be the estimation of unknown parameters I R , k r and δ , respectively. The adaptation law is deemed to be successful if Theorem 1 (CAS) Consider the attitude dynamics (7). The proposed control law is designed as where Î R , k r , δ are updated by the following adaptation laws for some α, γ 1 , γ 2 , γ 3 > 0 and ν 2 = ν 2 −ν 2 . Then the time-derivative of u = m U 1 + k t mẋ (cos φ sin θ cos ψ + sin φ sin ψ) + U 2 + k t mẏ (cos φ sin θ sin ψ − sin φ cos ψ)  (28) and (31), we can see that ν 2 (t) , Ĩ R , k r , and δ are bounded. To show the tracking error ν 2 is driven asymptotically to zero, we calculate the second time-derivative of Lyapunov function V (ν 2 ,Ĩ R ,k r ,δ) as From (31) it can be seen that ν 2 is uniformly bounded, and hence V (ν 2 ,Ĩ R ,k r ,δ) is bounded. This implies that V (ν 2 ,Ĩ R ,k r ,δ) is uniformly continuous. By using Barbalat's Lemma, we have lim t→∞ν2 (t) = 0 , and this implies that lim t→∞ e(t) = 0 . This completes the proof of the theorem.
It can be observed from (32) that the decay rate of the Lyapunov function depends on the model reference tracking error ν 2 (t) and we can guarantee the stability of the close-loop system even without having ν 2 (t) equals zero. On the other hand, looking at the adaptation law (28) along with the Lyapunov function (29) clarifies that since we are minimising all errors together for a non-zero ν 2 (t) , the adaptive rule updates itself constantly even the estimated parameters reach the actual value of the unknown parameter. Nonetheless, when ν 2 (t) converges to zero, we see from (28) that the estimated parameters converge to a constant value and there is no change in the Lyapunov function according to (30) and this may lead to a non-zero Lyapunov function in (29). Consequently, the convergence of the parameter estimation error to zero cannot always be guaranteed using the controller proposed in Theorem 1. Moreover, due to existent of ν 2 for the feedback control design, the controller is not economical for practical implementation. The schematic block diagram of the UAV control system designed according to Theorem 1 is presented in Fig. 2. Although the extra state in (24) can be removed from the control system structure, using CAS 14 , the controller in 14 suffers from the numerical stability due to the presence of the tracking error derivatives.

Theorem 2 (NAS) Consider the attitude dynamic in
(42) www.nature.com/scientificreports/ and 3, we can see that the control structure in NAS is simpler than CAS. This is due to the fact that the reference model is required in CAS to avoid ν 2 d for the feedback control design and generate ν 2 instead in the adaptation law in (28). The presence of ν 2 d in the control input may cause some numerical issues in both simulation and practical results. Nonetheless, the NAS is designed in a way that ν 2 is only required for the stability analysis and it doesn't appear in the adaptation law in (35).

Simulation and experimental results
In this section, the performance of our designed control system is evaluated numerically and demonstrated by practical implementation on a real drone. A virtual PD controller is implemented numerically as the translational controller. One way to design the gains of the PD controller is suggested in "Translational control design" section. The tuned values for these gains are K P = K D = diag(1) . The control gains for the rotational motion are designed according to the Theorems 1 and 2. The gains of (27) and (28) are tuned to be The gains of the controller in (34) and (35) are These gains α , γ 1 , γ 2 and γ 3 are usually tuned on a case-by-case basis and looking at the system response. Generally speaking, the higher the gains, the better the steady-state tracking results will be. However, we should always consider a trade-off to avoid control actuator saturation due to high gain or amplify measurement noises in the system outputs. The external disturbance h = sin(t) 0.8 sin(t) 0.6 sin(t) T with an unknown constant δ = 0.2 is added to the system dynamics. The parameters of the UAV used for both numerical simulation and practical implementation of the controller are listed in Table 1. The hardware platform used in this experiment is the parrot mambo minidrone as illustrated in Fig. 4. The simulation results using the proposed schemes are illustrated in Figs. 5, 6 and 7. The trajectory tracking errors for these results are plotted in Fig. 8. From the figures, we can see that both translational and rotational states of the quadrotor can follow the desired trajectories. This verifies the performance of the controller developed in "Proposed control design" section. The convergence of the parameter estimation errors is illustrated in Fig. 9 . We can see from this figure that NAS is capable in driving the parameter estimation errors to zero. However, CAS is not successful in converging the parameter estimation errors to their actual values. This is an example of the situation where there is no guarantee that the estimated parameters converge to the actual value of unknown parameters or perhaps the algorithm needs more time to converge to zero using CAS. Basically, this is caused by the inherent drawback of CAS, in which the parameter estimation error relies on the model tracking error.   www.nature.com/scientificreports/ In order to evaluate the performance of the proposed algorithm compared to the existing methods, various robust and adaptive control methods in which the system uncertainties are dominated have been reviewed in "Related works and main contributions" section. In these approaches only the tracking controller is designed to maintain the effect of the uncertainties on the system performance in a certain bound, however, in all these scenarios the unknown parameters of the UAV are not estimated. The classical adaptive method, i.e. Theorem 1, is developed using the certainty equivalence principle, and hence the boundedness assumption on the uncertainties is removed. The adaptive controller is not able to guarantee the estimated parameters converge to the actual value of the unknown parameters. This is an inherent drawback of the standard/classical adaptive method. However, in NAS, a new adaptive controller by adding an extra continuous derivative function is proposed to assist the algorithm to find the right direction by estimating the true value of parameters. NAS guarantees handling the parametric uncertainties with an appropriate design manifold. In the literature this problem is dealt with by linearizing the nonlinear term with unknown parameter. Therefore, the difficulty in estimating the unknown parameters of the system in a nonlinear dynamics is avoided. Another method proposed in the literature is by using intelligent computation such as neural networks and genetic algorithms. The tracking control and parameter estimation are achieved but with some level of residual error. A qualitative comparison of these different approaches is summarised in Table 2.
For further assessment of the proposed scheme in front of possible faults in the UAV, the algorithms are numerically evaluated for a scenario in which the UAV actuators experience an abrupt failure. The actuator faults can be represented by a constant parameter, named control effectiveness level of the actuator as explained in 34 . By changing this value between zero and one, it is possible to model different levels of actuator faults. Although the aim of this investigation is not to develop a fault tolerant control system with controller reconfiguration ability, here we evaluate the robustness of the controller for two scenarios for the actuator fault. These scenarios are illustrated in Fig. 10. In the first scenario, it is assumed that the UAV is operating with the nominal performance until 2s and all of the sudden four propellers are working with 10% of the full torque value τ . The parameters of the controllers are tuned to achieve the best nominal performance as shown in Fig. 5 and no changes on the control parameters is applied when an abrupt fault happens. The numerical results in Fig. 11a confirm that the angular rates in NAS is still capable to track the desired trajectories very closely, while CAS fails to maintain a reasonably well tracking performance for the UAV. In the second scenario, we assume multiple abrupt fault happens on τ and the control effectiveness level of this actuator is reduced. As before, the parameters of the controllers are tuned to achieve the best nominal performance according to Fig. 5 and no changes on the control parameters is applied when the abrupt faults happen. The numerical results in Fig. 11b show that the angular rates in NAS is still capable to track the desired trajectories very closely in this degradation scenario, while CAS fails to keep a good tracking performance for the UAV after the second level. For both scenarios the fitness of the tracking error in the steady state behavior of the rotational dynamic is calculated from t = 15s to t = 25s and the results are listed in Tables 3, 4, and 5. The fitness is calculated using the following formula  www.nature.com/scientificreports/ where ν 2 i is the state of rotational dynamics and ν 2 d i is the desired trajectory of ν 2 i . The simulation and fitness results, confirm that there are some major differences between NAS and CAS as listed below Remark 4 As can be seen from Fig. 9, the new scheme (NAS), is able to guarantee the convergence of unknown estimated parameters in the closed loop. There are many approaches proposed to estimate the unknown parameters, but the results for the systems with nonlinear dynamics are still limited in the literature as we stated in "Related works and main contributions" section (Related Works and Main Contributions). As a result, in NAS, we can generate an appropriate control input for the system by converging the adaptive estimations to the true values of the parameters while this is not the case for the CAS algorithm. Tables 4 and 5 and Fig. 11 in the new revision, the NAS shows much more robustness against abrupt actuator faults thanks to the presence of the additional continuous function.

Remark 5 As can be inferred from
In order to validate the performance of the controller experimentally, a Parrot Mambo minidrone was used. The Mambo minidrone platform provides an ideal solution for rapid prototyping of quadrotor control systems. It uses a combination of ultrasonic, air pressure, optical flow, gyroscope, and accelerometer sensors in order to estimate its position and attitude. The states of the Mambo minidrone are estimated using a Kalman filter for sensor fusion. MATLAB and Simulink were used to develop the proposed control system using the Parrot Minidrones Support package. The support package also provides a simulation platform for safe testing of the control    www.nature.com/scientificreports/ systems before deployment to the physical quadrotor. Once the controller was developed and tuned within MATLAB and Simulink, Simulink coder is used to convert the control system into C code, before uploading it wirelessly to the Mambo minidrone via Bluetooth connection. The control loop operates with a sample time of T = 0.005s . Data is then recorded on the Mambo minidrone for the duration of the flight. After the flight is completed, data from the flight can be downloaded to MATLAB via Bluetooth. Figure 4 shows the Mambo minidrone in the testing environment. This was a 3m × 3m × 2.4m gazebo with netting secured around it for safe testing of the Mambo minidrone. In addition to this, further software limitations were placed on the parrot minidrone to ensure safe operation. Firstly, saturation blocks were introduced on motor command signals to ensure that the signals remained within the maximum region for motor control. A crash detection feature was also present within the code that would cause the motors to shut off if the absolute value of x or y position of the minidrone was larger than a selected value. Finally, a crash would also be detected if the rate of change of the quadrotor positions was also above a selected value. Finally, the performance of NAS is also evaluated in a real scenario on an experimental UAV. The experimental results presented in Figs. 12, 13, 14 and 15 confirm that NAS is able to maintain the stable behaviour of the UAV and follow the desired trajectory closely in the experimental setting. More specifically, Figs. 12 and 13 illustrate the attitude rates of the UAV in the local coordinate system and the Euler rotational angles of the UAV in the global coordinate system respectively. Moreover, the profile of the translational motion of the UAV in the three dimensional space is illustrated in Figs. 14 and 15. To get a better view of the performance of the proposed adaptive algorithm in a real-life setting with parametric uncertainties, the fitness values for x, y and z in the steady state condition and from time 23s to 30s are calculated and listed in Table 6. The average fitness     Table 7, where the rise and peak times of x,y and z are calculated from 10s, 20s, and 0s, respectively. Here, the rise time of each position state is the time required to reach 10% and 90% of the desired position.   www.nature.com/scientificreports/

Conclusion and directions for future work
A fully tracking control for 6-DOF of UAV with unknown parameters is developed in this paper. A virtual PD controller is proposed to handle the tracking control of translational motions. Two adaptive schemes are designed and compared for rotational dynamics with the presence of unknown parameters in the nonlinear dynamics. The propeller inertia and translational drag are unknown for feedback control design. Also, an external disturbance containing an unknown constant is added to the dynamical system. CAS using the certainty equivalence principle is designed to handle the uncertain parameters. However, this scheme contains an inherent drawback, where the estimated parameters cannot be guaranteed to converge to the actual value of unknown parameters.
To solve this issue, we develop NAS by adding a continuous function to the control structure. The adaptive law is able to guarantee to handle uncertain parameters with the appropriate design of the manifold. Moreover, this scheme has a simple structure and is more practical to be implemented. Several simulations are presented for a mini-drone to demonstrate the effectiveness of our scheme. We also conduct an experimental test to the real UAV using NAS to validate our design. It will be interesting to extend this scheme for multiple heterogeneous UAV with fully unknown time-varying parameters. Figure 15. The profile of x, y and z in the experimental test using NAS in 3D.