Abstract
This paper presents a control scheme that allows height position regulation and stabilization for an unmanned planar vertical takeoff and landing aircraft system with an inverted pendular load. The proposed controller consists of nested saturations and a generalized proportional integral (GPI). The GPI controls the aircraft height and the roll attitude; the latter is used as the fictitious input control. Next, the system is reduced through linear transformations, expressing it as an integrator chain with a nonlinear perturbation. Finally, the nested saturation function-based controller stabilizes the aircraft’s horizontal position and the pendulum’s angle. Obtaining the control approach was a challenging task due to the underactuated nature of the aircraft, particularly ensuring the pendulum’s upright position. The stability analysis was based on the second method of Lyapunov using a simple candidate function. The numerical simulation confirmed the control strategy’s effectiveness and performance. Additionally, the numerical simulation included a comparison against a PD controller, where its corresponding performance indexes were estimated, revealing that our controller had a better response in the presence of unknown disturbances.
1. Introduction
The control of underactuated mechanical systems is a widely studied field and continuously increases knowledge, mainly because controlling this kind of system is challenging since it has fewer controllers than degrees of freedom. The inverted pendulum system is a classic example of an underactuated system. It consists of a freely spinning load around an axis and attached to a base that freely moves forward and backward—multiple authors have proposed control laws to stabilize this kind of system in its inverted position. Block and Spong [1] proposed the partial feedback collocated and noncollocated control law for the stabilization around the origin of the acrobot and the pendubot. Fantoni and Lozano [2] proposed a nested saturation control for the wheeled inverted pendulum that enables stabilization around the origin. Ibañez and Frias [3] proposed a nested saturation control for the nonlinear perturbed wheeled inverted pendulum, expressed it as a chain of integrators. The proposed control demonstrates asymptotical stability around the origin through the Lyapunov method when the pendulum angle is in the upper half-plane. The Furuta pendulum is another challenging system that has attracted the attention of several researchers, who have presented many interesting approaches to control this kind of pendulum. For instance, in [4], the authors proposed a Lyapunov-based control method for the stabilization around the origin, Furuta pendulum stabilization around the origin, while in [5], the authors introduced an active disturbance rejection-based control and its stability analysis. Another exciting example of underactuated systems is the unmanned aerial vehicle (UAV) such that this kind of system has been of high interest in the present century because of wide applications for different fields such as farming, photography, exploration, and military [6–8]. A typical example of a UAV is the planar vertical takeoff and landing (PVTOL) aircraft, a simplified model of the actual vertical takeoff and landing aircraft [9], which encompasses almost all the dynamics found in real UAVs. The PVTOL has been used as a suitable benchmark to test new and existing controllers because it behaves like the well-known quadrotor in a two-dimensional plane. There exist many works that tackle the PVTOL stabilization problem. For instance, Aguilar-Ibanez et al. [10] proposed a control scheme using controlled Lagrangian for PVTOL stabilization. Fantoni et al. [11] introduced a control scheme using the PVTOL attitude as fictitious control and nested saturation technique for the stabilization at the origin. In [12], Lozano et al. showed that once controlling the PVTOL aircraft height, the system resembles an inverted pendulum, and using a change of coordinates, it can be seen as an integrator chain with nonlinear disturbance. Also, in this study, the authors present a nested saturation-based control and demonstrated asymptotic stabilization at the origin. In [13], based on the image-based visual servoing method and the backstepping technique, the authors presented a novel control strategy to force a vertical takeoff and landing (VTOL) aircraft. Works [14, 15] can be valuable sources for the readers interested in this problem.
Combining the inverted pendulum and UAV systems adds a degree of freedom, obtaining a system noneasy to stabilize. Hehn et al. [16] proposed this problem in 2011 and named it the flying inverted pendulum, consisting of an inverted pendulum attached to a quadcopter. The control goals of this study are stabilization at the origin and tracking a circular trajectory for the quadcopter, using linear quadratic regulator (LQR) control in both cases. The fact that the pendulum weights less than 5 of the UAV weight allows us to separate its corresponding dynamics. Some works found in the specialized literature deal with the control of the flying inverted pendulum using different control ideas [17–20], and some others consider the control of UAVs carrying loads (commonly known as the flying crane). This problem is closely related and relevant to the central control problem in this study. For instance, Nicotra et al. [21] showed that the linearized model of a quadcopter with a suspended load could be stabilized at the origin using nested saturation functions. Pizetta et al. in [22] proposed a total feedback control with an auxiliary controller to accomplish tracking trajectory; almost all the references therein were developed to test the PVTOL system indoors, mainly to avoid counteracting the undesirable effect of the wind, which is not an easy task. However, techniques for nonlinear systems can be applied to obtain robust controllers, as the controller developed in this study. Of these techniques, perhaps the most used are backstepping control [23, 24], fuzzy control [25–27], active disturbance rejection control approach [28], and others [29]. To the authors’ knowledge, a general solution for stabilizing this type of system has not yet been reported in the literature.
In this context, we propose a control scheme for a PVTOL aircraft system with an inverted pendular load (PVTOL-ASIPL). This scheme mainly consists of a generalized proportional integral (GPI) controller and a nested saturation-based control. The GPI controller accomplishes height and roll attitude control, using the roll attitude angle as fictitious input control. Then, proposing a set of convenient linear transformations and reducing the system can be expressed as a chain of integrators with a nonlinear perturbation. Finally, we design a nested saturation function-based controller to stabilize the horizontal position and the pendulum angle. Applying the second method of Lyapunov assures boundedness of the whole state and asymptotic convergence to the origin.
The main contributions of this study are as follows:(i)An algorithm control that uses a fictitious control, and we propose a combination of a GPI controller and a controller-based saturation for the takeoff and landing maneuvers(ii)A set of convenient transformations, in which a high-order system can be expressed as a various low-order system(iii)A control strategy for the PVTOL aircraft system with an inverted pendular load controls the height, roll attitude, horizontal position, and roll angle simultaneously, even in the presence of exogenous disturbance
The organization of the rest of this study is as follows. Section 2 introduces the model of PVTOL-ASIPL externally perturbed, obtained from the Euler–Lagrange formalism. Section 3 develops the control scheme design, while Section 4 presents the results of the numerical simulations. Finally, Section 5 is devoted to the concluding remarks and future work.
2. Dynamic Model
This section presents the dynamic model of the PVTOL aircraft system with an inverted pendular load (see Figure 1). The dynamic equations were obtained by the Euler–Lagrange formalism as follows:where and are, respectively, the system kinetic and potential energies. Besides, the inverted pendular load base is in the PVTOL aircraft center of mass, , is defined as the angle between the PVTOL aircraft and the horizontal axis, and is the angle between the inverted pendular load and the vertical axis. Finally, the inverted pendular load center of mass is defined as

Therefore, Lagrangian of the system can be expressed aswhere is the PVTOL aircraft mass, is the pendular load mass, is the gravity force, is the PVTOL aircraft inertia, and is the inverted pendular load inertia.
The Euler–Lagrange equations of motion for the PVTOL-ASIPL system are in the form ofwhere is the generalized coordinate vector, is the input control vector, and is the external disturbance vector; without loss of generality, and .
Developing Euler–Lagrange equation (4) leads to
Then, the equations in (5a) represent the PVTOL-ASIPL, and they can be expressed in a compact form aswithwhere is a positive semidefinite matrix, is the Coriolis matrix, is the gravity matrix, is the input control matrix, and is the external disturbance matrix. The terms , , and are assumed to be unmodelled owing to the external disturbances and are defined as follows [30, 31]:
Finally, because is not a singular matrix, it is possible to represent the dynamical model of the PVTOL aircraft system with an inverted pendular load aswhere
2.1. Problem Statement
The control goal consists in proposing an algorithm to accomplish stabilization at the origin of the PVTOL aircraft system with an inverted pendular load. The maneuvers are trajectory tracking tasks involving a step-by-step procedure, consisting of (1) the stabilization of the coordinate ; (2) the stabilization of the coordinate ; and (3) further stabilization of coordinates and , even in the presence of disturbances due to the aerodynamic effects.
We have expressed the system in its minimal representation form, which will allow us to decouple it and simplify it. Instead of working with a complex system, we work with a few simple systems that embody the dynamics of the original one. We are now in a position to design and propose the control scheme in the following section.
3. Control Scheme
This section establishes the framework to solve the main control problem. To this end, please consider that the input control acts over plus disturbance . So, it is possible to design a control law for using a fictitious controller for the PVTOL aircraft with an inverted pendular load [9, 11]. Then, through linear transformations, a GPI law is used for the PVTOL takeoff and landing maneuvers. Once the GPI law stabilized the PVTOL aircraft height, a change of coordinates allows expressing the system as a chain of integrators nonlinearly perturbed, allowing to propose the nested saturation function-based stabilizing controller. Figure 2 presents the schematic diagram of the closed-loop system.

3.1. Controlling the PVTOL Aircraft Attitude
It is clear that the third equation (9c) consists of a double-chain integrator with nonlinear disturbance and control input . Then, a control law for tracking trajectory is searched for .
3.1.1. Control Statement
The dynamical equation for the roll attitude angle can be expressed aswhere is a lumped generalized disturbance input.
Also, according to Lozano Hernández et al. and Fliess et al. [31, 32], to overcome the lack of available measurements of , an integral reconstructor can be proposed, and using the local approximation of the disturbance input, it is possible to propose the control input aswhere is a smooth reference signal for the state , which is twice differentiable. Also, and , and the relation between the actual value of and is expressed by
Substituting (12) into (11) and expressing the resulting dynamics in terms of the tracking error, the following dynamics are obtained:where , with , are selected such that characteristic polynomial is Hurwitz, reducing the undesirable effects of the nonlinear disturbances [33, 34]. So, the dynamic error is exponentially asymptotically stable at the origin. This fact allows using as a fictitious control for subsystems (9a), (9b), and (9d). This proposal was previously used in other works dealing with the stabilization of PVTOL aircrafts [9, 12].
3.2. Simplified PVTOL Aircraft System with an Inverted Pendular Load
After and applying the controller , the following system of equations represents the PVTOL aircraft system with an inverted pendular load:withwhere and are the input controls.
Thus, we propose the following control laws [12]:where are auxiliary control inputs.
To obtain the dynamic model, the following change of coordinates is applied to model (15a)–(15c):
Therefore, models (15a)–(15c) transform into the following system:
Hence, system (19) takes the compact formwith
Then, given system (19) and ignoring the effect of the disturbances, the following partial feedback control can be proposed [1, 9]:where and are new auxiliary control inputs. Thus, the system defined by equation (19), in a closed loop with control model (23), leads us to obtain the following system:
Notice that the above system is the reduced model for the PVTOL aircraft system with an inverted pendular load, with and as the control inputs.
3.3. Stabilization of Height
A GPI controller with a saturation function is applied to obtain the height position control, allowing the tracking control to accomplish the takeoff and landing of the PVTOL aircraft system with an inverted pendular load.
3.3.1. Control Statement
Consider the vertical displacement described by (25), and let be the control input defined as (26):where is a smooth reference signal twice derivable and .
Also, is a saturation function defined as
Thus, the tracking trajectory error is exponentially asymptotically stable.
Proof. The proof is split into two parts. First, it is proven that a saturation function in finite time bounds the error dynamics . Then, we demonstrate exponentially asymptotic stability.
3.3.2. Error Bounded
The tracking trajectory error is defined aswhere is the desired height position.
Let us define the following state variables as
Therefore, the dynamic error is transformed into the following system:where is the control input.
Now, to use the second method of Lyapunov, consider the following candidate function:which is positive definite, with time derivative
Thus, the first and second time derivatives of are known and bounded as . Besides, in a close neighbourhood where and are such that and , the saturation function fulfills . So, parameters were designed as in such that is dominant, and . Finally, the last terms of equation (31) satisfy the following:
Therefore, computing the time derivative of the candidate Lyapunov function leads toand if and are small enough such that , then and are bounded in finite time. Thus, after time , when this condition is satisfied, control law (26) takes the following structure:
From the above, equation (24b) is expressed aswhere is a lumped generalized disturbance input.
Then, to overcome the lack of available measurements of , the following integral reconstructor is introduced [31]:
Using the local approximation of the disturbance input, the relation between the actual value of and is expressed by
Now, differentiating equation (38) and substituting into (36), the control law takes the formwith and .
Substituting (39) into (36) is possible to express the dynamic error as follows:
The coefficients are selected such that the polynomial is Hurwitz, eliminating the nonlinear disturbances [33]. Then, tracking trajectory error is exponentially asymptotically stable to zero.
3.4. Control for (, )
This section proposes a nested saturation-based controller strategy to stabilize the PVTOL-ASIPL horizontal position and roll angle. Notice that the nested saturation-based controller strategy allows stabilizing nonlinear systems that can be approximately expressed as a chain of integrators [35–37]. Thus, our stability problem will be solved as follows. First, we introduce a linear transformation for the stabilizing controller. Then, we demonstrate that the proposed controller guarantees boundedness and converges to zero, in finite time, of the whole state.
3.5. Expressing the PVTOL-ASIPL as a Chain of Integrators
After the application of the controller , system (25) can be reduced to the subsystem (, ) as follows:
Then, introduce the following coordinates’ change:
System (42) takes the following form:where is considered as a nonlinear disturbance and is an unknown constant that depends on variables , and (Note that .). Notice that subsystem (43) is similar to the cart-pole system plus a nonlinear disturbance [38].
3.5.1. Control Statement
Based on work [3], we define the following state variables , , , and . Then, we express the dynamic system as
To express system (44) as a chain of integrators with a nonlinear perturbation and propose a controller for the stabilization of the subsystem (, ), applying the decoupling theorem [3, 39], the following global nonlinear transformation is introduced:
Hence, the transformed system into a chain of integrators is given by
3.6. Nested Saturation Function-Based Controller
Inspired by Teel [40], a linear transformation is proposed to obtain the stabilizing controller for system (46) as follows:
Applying transformation (47) to system (46) leads tofor which the following nested saturation function-based stabilizing controller is proposed:where is defined by equation (26), and with is given by
Remark 1. Note that when and when and . Besides, when and when and . Additionally, when and .
3.7. Whole State Boundedness
Now, we prove that the proposed controller (49) ensures whole state boundedness. Moreover, the bound of each state directly depends on the designed parameters of the controller.
Step 1. A positive definite function is defined asThe time derivative of is expressed bywhere and are selected such that . It is clear that when ; therefore, there exists a finite time such thatThus, when this condition is satisfied, the control law (49) takes the following structure:
Step 2. Behaviour analysis of : let us consider the following positive definite function:Differentiating it with respect to time and after substituting (54) into , the following expression is obtained:To ensure that is achieved, the following conditions must be satisfied:Then, there exists a finite time after whichTherefore, is bounded, and the control law takes the following form:
Step 3. The following positive definite function is introduced:Differentiating and after substituting (60) into the second equation of (48), the following is obtained:where and must satisfy the relation It is easy to observe that if , . Hence, there exists a finite time after whichTherefore, is bounded, and the control is revealed to be
Step 4. Substituting (63) into the first equation of (48), we obtainTo demonstrate that is bounded, a positive definite function is defined as follows:Differentiating along the trajectories of (66) leads towhere must be selected so that and to achieve . Therefore, there exists a finite time such thatConsequently, is also bounded. Finally, the values of parameters , and can be determined by the following inequalities:From the above conditions, the set of control parameters can be selected aswhere is directly related to the magnitude of the system disturbances.
3.8. Whole State Convergence to Zero
Here, we prove that the closed-loop system, provided by (48) and (49) and satisfying (70), is asymptotically stable.
Notice that, after , the control law (49) is no longer saturated; that is,and the closed-loop system turns out to be
To demonstrate whole state convergence to zero, the following Lyapunov function is used:where , and differentiating along the trajectories of equation (72), it is obtained thatandwith being positive definite with .
From Remark 1, it is shown that the following relations are satisfied:
So, fulfills
Hence, the previous inequality is strictly negative definite, and the following constraints are obtained:
Therefore, if the restrictions of (77) are fulfilled, is strictly negative, and the state vector exponentially converges to zero after .
The following proposition summarizes the previous discussion, which is the main result of our study.
Proposition 1. Consider the PVTOL aircraft system with an inverted pendular load as described in (9) and in a closed loop with controllers (12), (17a) and (17b), (23), (25), and (49). Then, the closed-loop system is exponentially asymptotically stable provided that the control parameters , , , and satisfy inequalities (69), and and are selected such that the characteristic polynomial is Hurwitz.
Finally, in Figure 3, the steps’ sequence is shown, obtained from the control laws and .
Having designed the control scheme for the PVTOL aircraft system with an inverted pendular load and carried out its convergence analysis, we test its effectiveness through numerical simulations in the following section. It is worth mentioning that it would be ideal for testing the scheme by conducting actual experiments—unfortunately, the construction of the needed prototype is still in progress.

4. Numerical Simulations
To test the controllers’ performance, we carried out some numerical simulations using the MATLAB-Simulink program, and the results were obtained based on the numerical method of Runge–Kutta of the fourth order with a fixed step of 0.01 s. The physical parameters of the system are kg, kg, m, and . Also, the tuning parameters proposed for each controller are listed in Table 1. Notice that the controller parameter values and were selected such that the error dynamics is equal to the desired closed-loop polynomial . Besides, the parameters of were chosen such that [40]. Finally, the initial conditions were set as follows: m, m, rad, rad, m/s, m/s, rad/s, and rad/s.
In the first experiment, the control strategy simultaneously carries out height position by performing trajectory regulation, stabilizing the horizontal position and roll angle for the PVTOL and the pendulum angle. In this case, the desired trajectory, , was proposed as
To test the effectiveness of the control strategy introduced in this study, we carried out a comparison test against the classical PD controller, with gains tuned as and . Figure 4 shows the outcome of this experiment, where we can observe that both strategies successfully achieve the height position regulation through the trajectory tracking task and the stabilization of the horizontal position, the roll angle, and the pendulum angle. Also, the plot in this figure shows that our controller converges to the desired values faster than the PD controller does. Figure 5 shows the behaviours of the angular velocities, and Figure 6 shows the tracking error between and and the control inputs and .

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)
Using the same setup as before, we run a second experiment, but in this case, the system is affected by external disturbances with parameters fixed as , , and . Figure 7 shows the corresponding plots, where we can see that our controller is capable of accomplishing height tracking and position regulation and, simultaneously, stabilizing the horizontal position, roll angle, and pendulum angle. Also, it can be seen that the PD controller used in the first experiment regulates the system slower than our controller does and exhibits undesirable oscillations. Therefore, our controller has better performance, maneuverability, and whole stabilization, even when the system is affected by external disturbances. Figure 8 shows the system angular velocities, and Figure 9 shows the tracking error behaviour between and and the control inputs and . Please notice that the proposed control has adequate energy management, according to the obtained performance index shown in Figure 10.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(a)

(b)
5. Conclusions
This work presents a nested saturation function-based controller, combined with a GPI controller, to stabilize the PVTOL aircraft system with an inverted pendular load. The model of this system was derived using Euler–Lagrange formalism. The main contribution consists of using a fictitious control, and then a GPI controller is proposed for the aircraft angle . Several linear transformations and coordinate changes were introduced to express the original system to a minimal representation. To accomplish the takeoff and landing maneuvers, we propose a GPI controller to track the desired trajectory. After stabilizing the PVTOL height, the system was represented as a chain of integrators plus nonlinear disturbance, allowing us to use nested saturation functions to design a controller to stabilize the horizontal position and pendulum angle. The stability analysis was carried out using the second method of Lyapunov, using a simple candidate function. Designing the control scheme was not an easy task because the PVTOL system with an inverted pendular load is underactuated, and ensuring the pendulum’s upright position makes this problem even harder to solve. We ran numerical simulations to assess the performance of our control scheme, having obtained convincing results. These simulations included a comparison against the well-established PD control strategy. The performance index of both controllers was computed to compare them, and the outcome revealed that our strategy has a better performance than the PD controller. It is important to note that the performance indexes were estimated in the presence of nonlinear perturbations, which means that the proposed controller behaves well even in this undesirable yet unavoidable realistic scenario. It is important to note that the controller, based on a GPI controller and nested saturation functions, allows us to perform takeoff maneuvers in the presence of exogenous disturbances.
In future work, we will explore a design to estimate the disturbance due to wind or robust techniques for parametric uncertainties of the system. In addition, it is worth mentioning that an experimental platform that allows configuring the PVTOL with an inverted pendular load system has been designed, whose construction is in process, Thus, experimental implementation of the control scheme proposed herein is considered.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This research was funded by Secretaria de Investigación y Posgrado-Instituto Politécnico Nacional grants nos. 20210268, 20210253, and 20211168. C. Alejandro Villaseñor Rios thanks the support from the CONACYT. Cesar Alejandro Villaseñor Rios, Octavio Gutiérrez-Frías, Carlos Aguilar-Ibanez, and Miguel S. Suarez-Castanon are at Instituto Politécnico Nacional and they are the ones to express their gratitude for the support received.