Advanced Mathematics and Numerical Modeling of IoT (Internet of Things)
View this Special IssueResearch Article  Open Access
HyeokJune Jeong, Myunggwon Hwang, Hanmin Jung, Youngguk Ha, "Mathematical Modeling of a Multilayered DriftStabilization Method for MicroUAVs Using Inertial Navigation Unit Sensor", Journal of Applied Mathematics, vol. 2014, Article ID 747134, 11 pages, 2014. https://doi.org/10.1155/2014/747134
Mathematical Modeling of a Multilayered DriftStabilization Method for MicroUAVs Using Inertial Navigation Unit Sensor
Abstract
This paper proposes a multilayered quadrotor control method that can move the quadrotor to the desired goal while resisting disturbance. The proposed control system is modular, convenient to design and verify, and easy to extend. It comprises three layers: a physical layer, a displacement control layer, and an attitude control layer. The displacement control layer considers the movement of the vehicle, while the attitude control layer controls its attitude. The physical layer deals with the physical operation of the vehicle. The two control layers use a mathematical method to provide minute stepbystep control. The proposed control system effectively combines the three layers to achieve drift stabilization.
1. Introduction
Unmanned aerial vehicles (UAVs) are expected to become a major part of the aviation industry as they can perform tasks such as traffic control, video recording, reconnaissance, and surveillance. Concomitant with developments in computer science, automatic control, sensors, and communications technologies, the quadrotor, in particular, is being evaluated as a suitable platform for small UAVs. Among its many advantages is the fact that it can move both vertically and horizontally, it can be made very small, and it can carry a variety of electrical devices.
A quadrotor flies by means of four propellers, which are controlled by an automatic system programmed in a microprocessor. This means that the control system has minute control over its flight. Consequently, it is relatively easier to fly a quadrotor than other aircrafts such as helicopters and airplanes.
However, controlling the desired behavior of the quadrotor is never easy. Because it flies in air, the quadrotor has to overcome inertia and undesirable wind disturbance.
Consequently, in order for the quadrotor to be useful, it has to have an effective control system.
Efforts have previously been made to solve this problem; however, this goal has not yet been definitively achieved. The control systems proposed thus far are more complicated than effective. Some researchers have proposed visionbased control systems; however, such systems are inefficient and heavy and therefore not suitable for small UAVs.
In this paper, we propose a quadrotor control system that uses the four motors and nine EOFaxis inertial navigation sensors, which measure the quadrotor’s attitude and movement, to achieve minute control.
The control system operates on the basis of positioning; that is, it can move the quadrotor to any desired point while resisting external forces such as wind. Of course, drift stabilization is possible in any current position. In particular, the proposed control system is implemented in a modular form, so its design efficiency and performance can be easily verified.
The proposed controller significantly increases the quadrotor’s stability and hovering performance and thereby facilitates the use of the quadrotor for various applications.
2. Related Work
Over the past decade, much research has been devoted to quadrotor control systems, a large portion of which has focused on the use of vision. Grzonka et al. [1] pioneered work in this area. Their research focused on indoor quadrotor flight control systems that are able to pilot the quadrotor in indoor spaces via realtime image processing [2]. However, even though their research is excellent, the stability of their control system is in doubt. Bills et al. [3] also proposed a control system that uses vision; however, their proposed system cannot guarantee stable quadrotor movement.
Romero et al. [4] proposed a control system that uses optical signals. The proposed system is stable because it only traces certain points. Gu et al. [5] focused their research on systems flying in formation, with each UAV designed hierarchically according to its role, such as flight leader or wingman. These systems have the advantage of stability but are limited to operating only in particular situations.
Zhang et al. [6] proposed a control system that combines vision with an IMU sensor. The system is stable; however, the efficiency of the control system is poor—it performs well for small vehicles but is too large for microquadrotors.
Consequently, in spite of research efforts expended to date, a suitable implantable control system for microUAVs has not been realized. We believe our proposed multilayered driftstabilization method can provide an effective solution that resolves these problems.
3. Design of the Proposed Control System
In this section, we discuss the design of the proposed control system in terms of its layers and the algorithms used by each layer.
3.1. The Combined Layers
Our proposed multilayered control system has advantages in terms of its design, verifiability, and extendibility. These advantages result from the fact that each control layer is in charge of an assigned function. That is, each layer is considered to be abstracted.
The architecture of our proposed multilayered driftstabilization control system is depicted in Figure 1. It comprises three layers. The first layer consists of physical components such as motor, frame, rotor, and battery. The rotation of the motor is controlled in accordance with a control value. The sensor calculates such parameters as the angular velocity and acceleration and passes the feedback values to the other layers.
The second layer, the “attitude control layer,” is in charge of the attitude angle, “roll, pitch, and yaw.” The attitude control layer receives attitude feedback and an objective angle. However, the layer does not care about how the actual angle is calculated. This layer is simply responsible for ensuring that the quadrotor is at the correct angle received.
Finally, the “displacement control layer” is in charge of the movement of the vehicle. This layer receives displacement feedback and control signals (reference) that it uses to calculate an appropriate angle to pass to the attitude control layer. The displacement control layer does not care how the angle of the vehicle is controlled. This is why we say that the control system design is abstracted.
Consequently, the proposed multilayered control system can be a powerful and convenient system, in which each section only fulfills its specified roles.
3.2. The Physical Layer
A quadrotor has four rotors and four motors on a rigid framework. A sensor is mounted at the center of the frame for accurate operations. Other components, such as circuit and batteries, are mounted in the same position as the sensor to ensure a balanced center of gravity.
Our proposed control system is designed for the shape shown in Figure 2 and takes into consideration the dynamics of this shape.
Despite the shape of the camera, the quadrotor is assumed to be symmetrical in quality and structure. The physical characteristics of the quadrotor are listed in Table 1.

In order to establish the dynamic model of the quadrotor, we can make the following general assumptions.(i)Gravity and resistance of the quadrotor are not affected by flight altitude and other factors.(ii)Thrust in all directions is proportional to the square of the rotor speed.(iii)The quadrotor is a symmetrical rigid body.(iv)The origin of the inertial coordinate system is in the same position as the geometric center and the centroid of the quadrotor.
Two main effects are taken into consideration: generation of the thrust and the drag force. The thrust, , produced by each motor is a force calculated as where is the thrust coefficient, is the air density, is the rotor disk area, and is the blade radius. Further, the drag force is defined as where is the drag force, is the drag force coefficient, and is the speed of the quadrotor.
In the fixed coordinates of the body, the direct inputs are revolutions per minute (RPM) commands for the motors. The resultant outputs are direction thrusts in these coordinates.
However, the outputs under consideration are the attitude and position. To eliminate this gap, four control variables are defined as follows: where , , , and are the respective rotational speed of each rotor, , , , and are the lift forces of the motor’s axis, is the length of the quadrotor, is the total lift force, is the rolling moment, is the pitching moment, and is the pitching moment.
3.3. The Attitude Control Layers
Let us now derive a mathematical model for the quadrotor shown in Figure 2 (see Figure 4). But this layer is not our main interest. For that reason, this layer is designed using existing excellent research. Descriptions, expressions, sentences, and equations also are quoted in those papers (especially reference [7]).
The origin of the inertial coordinate system E is the initial position of the quadrotor. The positive direction of the axis is the designated heading of the quadrotor and is perpendicular to the horizontal plane.
This coordinate system is used to study the relative movement of ground and quadrotor. The axis is perpendicular to the plane. The quadrotor’s spatial coordinates (, , ) can be obtained through the inertial coordinate system.
The origin of quadrotor coordinate system B () is the center of the quadrotor, and is parallel to the center connection of the front and rear rotors and the positive direction points to the front. is parallel to the center connection of the left and right rotors and the positive direction points to the right. The axis is perpendicular to the plane; the positive direction is the direction conforming to the right hand rule.
These two coordinates can be converted to each other through transition matrix .
Euler angles are defined as follows.
Pitch angle : angle between the axis and the projection of in the plane.
Yaw angle : angle between the axis and the projection of in the plane.
Roll angle : angle between the axis and the projection of in the plane.
Consequently, we can obtain the transition matrix , which is from the quadrotor coordinate system to the inertial frame. Consider
With attitude angles defined as in Figure 2, the transformation matrix from the inertial coordinates to the fixed body coordinates is
Then, we can obtain
By Newton’s second law of motion
By Newton’s second law, the dynamic equation of the quadrotor, the line motion equation can be obtained. It is defined as follows: where, , , and is the air resistance.
According to the relationship between Euler angle and angular velocity of the quadrotor, the following result can be obtained:
Expressed with regard to , we obtain
Then, this equation can be expressed as follows:
Finally, the following expression is obtained:
The quadrotor was previously assumed to be symmetrical in quality and structure, so the inertia matrix. So it can be defined as a diagonal matrix: where , , and are the rotary inertia around the , , and axes, respectively.
By calculating the angular momentum, we can obtain the three axial components’ angular motion equations of in the quadrotor coordinate system: , , and . Consider
After simplification, the formula becomes
Equations (3)–(14) are mathematically considered state equations and can be implemented at the attitude control layer.
3.4. The Displacement Control Layer
The displacement control layer controls the movement of the quadrotor to the desired location while resisting disturbances. At this layer, the control system adjusts the attitude of the quadrotor to make it move. That is, this layer receives the input value for position, which it then uses to calculate the proper angle for the quadrotor.
The algorithm is derived from situations such as those shown in Figure 5. The is known, and the force of axis can be obtained by flowing equations. Consider
The quadrotor is tilted to the axis. This angle is defined as theta. The force is generated in the axis by the rotation of the rotor:
Equation (19) is Newton’s second law on the axis. Consider
Equation (21) can be obtained by substituting in (18) and (19). Then, an angle can be obtained for the desired acceleration. The respective domain and range of the sine function must be satisfied (22). The physical meaning of (21) is as follows.(i)In order to obtain a large acceleration, a large angle is necessary.(ii)A large mass requires a large angle.(iii)If the motor output is large, the desired angle is reduced (Figure 3).
Given the angle of the quadrotor, the respective acceleration can be derived easily; however, this system needs the appropriate angle of the quadrotor’s respective displacement. Therefore, we propose another equation derived from (18).
If the system is designed to simply resist the acceleration, it cannot achieve its desired target because of the value of the instantaneous acceleration. Our proposed control system needs to know the tendency of the movement due to acceleration. For this reason, the acceleration in (18) is replaced by a different function:
The conversion of expression (23) does not use the mathematical meaning of equal. The transform function is inferred from (18): where is a desired location in the absolute coordinate system and is the current position in the absolute coordinate system. The function (24) means the distance from the target point, in this sense called “error” in control engineering: where is the previous position in the absolute coordinate system and is the acceleration of the disturbance of the axis. This means that the current position is the sum of the previous position, the current attitude of the aircraft, and the acceleration of the disturbance.
That is, the transform function refers to the error value in the PID feedback control shown in Figure 6, and the feedback control system can minimize the error by adjusting the process control outputs. Consider where is a variable of integration and takes values from time zero to the present time, . , , and signify the gain, a tuning parameter. Consider
The appropriate angle can be calculated using (28). The physical meaning of this expression is as follows (Figure 7, left block).(i)A large mass requires a large angle.(ii)If the motor output is large, the desired angle is reduced.(iii)The larger the control signal, the more the angle tilts.
The displacement control layer enables the quadrotor to move to the desired point while resisting disturbance. Invariably, it also makes it possible for the quadrotor to hover. However, this control layer calculates only the angle for the respective distance error. In other words, this layer is designed to faithfully conduct an indigenous mission.
4. Implementation and Experimentation
We implemented our proposed control system in MATLAB Simulink and conducted several simulations (Figure 7). The experimental values obtained are listed in Tables 1 and 2.

Optimization is not within the scope of this paper. As a result, approximate values were used. To verify our proposed multilayered control system, only layers 1 and 2 were implemented in the first set of experiments.
In Figure 7, the left block is the attitude control layer and the right block is the physical layer. The remaining small blocks are the input and output gates.
Figure 8 shows the results of the first set of experiments. As expected, in situations where there is a continuous disturbance, the quadrotor is able to maintain its attitude. However, it can be seen that the quadrotor is being gradually pushed by the disturbance. The quadrotor control system is intended to give results such as those shown in Figure 8.
(a) Continuous disturbance
(b) Position of quadrotor
(c) Attitude of quadrotor ( , )
In the second set of experiments, we added layer 3 to the system.
The left block is the displacement control layer. It receives position feedback and outputs the appropriate angle (Figure 9).
Figure 10 shows the inside block that implements (26) and (27). With the addition of layer 3 to the system, the results depicted in the figure were obtained.
It can be seen in Figure 11 that the attitude changes to resist the disturbance. This appropriate change in the attitude enables the quadrotor to have drift stabilization.
(a) Step disturbance of the axis
(b) Position of quadrotor
(c) Attitude of quadrotor ( )
The control system knows the acceleration due to the quadrotor control position; thus, control is possible. Figure 12 shows that it is moved to a point that maintains its objective.
(a) Position of quadrotor ( axis)
(b) Position of quadrotor ( axis)
(c) Position of quadrotor ( axis)
Figure 13 shows that it moves to a point to maintain its objective even though the disturbance is very irregular. When there are large angular changes it is designed to move to a certain point, whereas for small angular changes it is designed to resist the disturbance.
(a) Disturbance of axis
(b) Position of quadrotor ( axis)
(c) Attitude of quadrotor (φ)
The final experimental results show that our proposed multilayered driftstabilization control system can be a useful and powerful solution for microquadrotors.
5. Conclusion
The design and implementation of mathematical modeling of our multilayered driftstabilization method were successful. Each of the layers is designed separately and drift stabilization is achieved by cooperation with each layer.
Our proposed system provides powerful solutions for some of the problems that have not yet been solved. The system can be mounted on small vehicles and can make the quadrotor move to its objective even when there are irregular disturbances. The proposed control system has good performance and is easily expandable. Therefore, it is a useful microUAV control system.
However, this control system presently has an integration error. In order to overcome this limitation, we plan to implement a process that corrects it using vision. The system also lacks optimization, particularly in the PID term. Therefore, we hope to conduct further studies from which we can present optimized results that further enhance the usefulness of our proposed control system.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research was supported by the Basic Science Research Program of the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (Grant no. 2012006817).
References
 S. Grzonka, G. Grisetti, and W. Burgard, “A fully autonomous indoor quadrotor,” IEEE Transactions on Robotics, vol. 28, no. 1, pp. 90–100, 2012. View at: Publisher Site  Google Scholar
 S. G. Fowers, D. Lee, B. J. Tippetts, K. D. Lillywhite, A. W. Dennis, and J. K. Archibald, “Vision aided stabilization and the development of a quadrotor micro UAV,” in Proceedings of the IEEE International Symposium on Computational Intelligence in Robotics and Automation (CIRA '07), pp. 143–148, Jacksonville, Fla, USA, June 2007. View at: Publisher Site  Google Scholar
 C. Bills, J. Chen, and A. Saxena, “Autonomous MAV flight in indoor environments using single image perspective cues,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '11), pp. 5776–5783, Shanghai, China, May 2011. View at: Publisher Site  Google Scholar
 H. Romero, S. Salazar, and R. Lozano, “Realtime stabilization of an eightrotor UAV using optical flow,” IEEE Transactions on Robotics, vol. 25, no. 4, pp. 809–817, 2009. View at: Publisher Site  Google Scholar
 D. L. Gu, G. Pei, H. Ly, M. Gerla, and X. Hong, “Hierarchical routing for multilayer adhoc wireless networks with UAVs,” in Proceedings of the 21st Century Military Communications Conference (MILCOM '00), pp. 310–314, October 2000. View at: Google Scholar
 T. Zhang, Y. Kang, M. Achtelik, K. Kühnlenz, and M. Buss, “Autonomous hovering of a vision/IMU guided quadrotor,” in Proceedings of the IEEE International Conference on Mechatronics and Automation (ICMA '09), pp. 2870–2875, Changchun, China, August 2009. View at: Publisher Site  Google Scholar
 J. Li and Y. Li, “Dynamic analysis and PID control for a quadrotor,” in Proceedings of the IEEE International Conference on Mechatronics and Automation (ICMA '11), pp. 573–578, College of Mechanical and Electrical Engineering China Jiliang University, August 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 HyeokJune Jeong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.