International Journal of Aerospace Engineering
Volume 2018 (2018), Article ID 5306125, 10 pages
https://doi.org/10.1155/2018/5306125

Research Article

Design and Simulation of Eight-Rotor Unmanned Aerial Vehicle Based on Hybrid Control System

Xueqiang Shen, Jiwei Fan, and Haiqing Wang

Northeast Electric Power University, Jilin 132000, China

Correspondence should be addressed to Jiwei Fan; 752341471@qq.com

Received 10 June 2018; Accepted 3 October 2018; Published 5 December 2018

Academic Editor: Jacopo Serafini

Copyright © 2018 Xueqiang Shen 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.

Abstract.

In order to control the position and attitude of unmanned aerial vehicle (UAV) better in different environments, this study proposed a hybrid control system with backstepping and PID method for eight-rotor UAV in different flight conditions and designed a switching method based on altitude and attitude angle of UAV. The switched process of hybrid controller while UAV taking off, landing, and disturbance under the gust is verified in MATLAB/Simulink. A set of appropriate controllers always matches to the flight of UAV in different circumstances, which can speed up the system response and reduce the steady-state error to improve stability. The simulation results show that the hybrid control system can suppress the drift efficiently under gusts, enhance the dynamic performance and stability of the system, and meet the position and attitude of flight control requirements.

1. Introduction

In recent years, multirotor UAVs have attracted increasing attention. As early as World War II, many countries, such as the United Kingdom and the United States, have investigated in the development of UAVs. With its flexibility, simple structure, and low cost, multirotor UAV has broad prospects for development in military and scientific research and many other fields [1, 2]. Compared with fixed wing UAVs, the eight-rotor UAVs can hover, vertically takeoff and landing, and perform indoor flight. The control of the eight-rotor UAV is mainly divided into attitude control and position control. The change of the attitude of the eight rotors will cause a corresponding position change. Thus, the flight attitude control of the eight rotors is critical to the control of the technology as a whole. The common control algorithms are classical PID control [3], fuzzy PID attitude control [4], neural network control [5], sliding mode control [6], backstepping control [7], dynamic surface control [8], visual feedback control, and adaptive control [9]. Those methods are based on the mathematical model of the normal flight of an eight-rotor UAV. However, it is difficult to establish an accurate mathematical model in practice. This situation leads to the discrepancy between the actual control effect of the eight-rotor UAV and the simulation results. In the literature [10], the problem of stabilization of the three-rotor UAV is studied using the method of fuzzy control and pole assignment. However, fuzzy control of the actual system will lead to the decrease of control precision and the variation of the dynamic quality. The improvement of the control accuracy will inevitably increase the fuzzy rules of the quantization series, which will lead to the expansion of the scope of the rule search and the reduction of the decision-making speed, leading to poor real-time performance of the system. In [11], the hybrid high-order sliding mode control method is addressed to control the eight-rotor UAV, but the effect of suppressing the jitter caused by sliding mode control is not very good, and the control precision is insufficient. In [12], the vertical takeoff and control of the four-rotor UAV are studied based on the state feedback and the neural network adaptive hybrid control method. The approach achieves the height control of the four-rotor UAV; however, the attitude angle control is neglected. Therefore, the important problem in the study of the UAV control system is to design a UAV with stronger ability to control the flight attitude and the adaptability to the environment based on the hybrid control theory.

In view of the above shortcomings and disadvantages, this paper presents an eight-rotor UAV controller design based on a hybrid control method. Adopting the idea of hybrid control, the backstepping controller, PID controller, and hybrid switcher are combined to realize the switching of the controller in different flight modes. It not only inherits the reliability, low static error, and interference ability of a traditional PID algorithm but also maintains the advantages of backstepping control method. The two methods complement each other to achieve good control of the UAV in a complex environment and to improve the control accuracy.

2. The Building of the System Model

The eight-rotor UAV is an underdamped system with four inputs and six outputs, according to the distribution of the rotor. Given its overall layout, it belongs to the noncoaxial multirotor aircraft [13]. The eight-rotor UAV works on the principle of changing the eight motors’ velocity to drive the propeller and the paddle. The change of motor speeds results in the alteration of lift and spin and affects the vertical, front, rear, left, and right planes of movement. Through these basic combination of actions, the UAV is able to complete the diverse range of missions. This paper mainly focuses on the study of the uniformly distributed eight-rotor UAV.

In order to study the dynamic performance of the eight-rotor UAV better, a relatively accurate mathematical model of the eight-rotor UAV must be established. Without losing the generality, a sufficient condition for the model is established.

(a)The body is a uniform symmetrical steel structure(b)The inertial coordinate system is the ground coordinate system(c)Air resistance and gravity suffered during flight are not changed by flight altitude(d)Body coordinate system B origin and the body of the center of mass and geometric center are in the same position(e)The lift generated by the propeller is proportional to the square of the motor speed

Based on the above assumptions, the motion force structure and dynamic analysis of the eight-rotor unmanned aerial vehicle are implemented, and the Euler angles are defined and two basic coordinate systems are established: the inertial coordinate system E () and the body coordinate system B (), as shown in Figures 1 and 2.

(a)

(b)

(c)

Figure 1: Sketch of roll angle, pitch angle, and yaw angle.

Figure 2: Force analysis of eight-rotor UAV.

From Figures 1 and 2, the transformation rotation matrix between the body coordinate system and the inertial coordinate system is as follows:

The mathematical model of the line motion of the eight-rotor UAV given by the kinetic analysis of Newton’s second law is as follows [14]:

where is defined as the quality of UAV, as the gravity of body, as the flight speed, as the displacement, as the role of the UAV on the external force, as a single rotor rotation generated by the lift, with as the lift coefficient, and as the motor speed.

From the angular momentum theorem of centroid motion, the angular motion model of the eight-rotor UAV is written as follows:

According to the relationship between the angular velocity of the eight-rotor unmanned aerial vehicle and the Euler angle, (4) can be reduced to

where is the distance from the center of the propeller to the origin of the coordinate system of the body; , , and are the moments of inertia of the , , and shafts, respectively; , , and are the angular velocity components of the -axis, -axis, and -axis, respectively.

The control input for the eight independent control channels of the eight-rotor UAV is determined by the rotor speed as follows:

where is defined as the vertical direction control input quantity, as the control direction for the rolling direction input, as the input for the pitch direction, as the control input for the yaw direction, and as the air resistance coefficient.

In summary, the mathematical model of the eight-rotor UAV system is as follows:

3. The Design of Hybrid Control System

The purpose of the design is to provide the UAV with a set of appropriate control strategies in different states or flight environments, in order to change the flight control curve appropriately [15]. The flight control of the eight-rotor UAV studied in this paper adopts the method of backstepping control and PID control and through the hybrid control method attempts to control the eight-rotor UAV in different environments and improve the flight stability and antijamming capability of the UAV. When the UAV is in normal flight, it can be approximated that the smaller the value of roll angle and pitch angle, the smoother the flight state. However, in the actual situation, when the UAV flies in the takeoff and landing state, subjected to gust and other interference factors, it will often cause the UAV to lose stability or even crash out of control. Therefore, the flight height , roll angle , and pitch angle are taken as the switched variables between the controllers, that is, when the height, roll angle, and pitch angle exceed the set threshold, the backstepping controller will be switched to the PID controller, to achieve enhanced anti-interference ability. When the UAV is faulty or the roll angle and pitch angle exceed the threshold, the backstepping controller or PID controller can be switched to the runaway protection controller for motor braking and landing.

From the above mathematical model of the eight-rotor UAV system, we can see that the angular velocity and attitude angle do not depend on linear motion. We also see that the spatial position of the line motion depends on the angular motion, and there is a semicoupling relationship between the line motion and the angular motion. The control system of the eight-rotor UAV can be divided into the inner ring (attitude loop) and outer ring (position loop) control [16]. Therefore, the eight-rotor unmanned aerial vehicle control structure based on hybrid control system is shown in Figure 3, where are the desired position input and are the desired gesture input.

Figure 3: Eight-rotor UAV control system structure.

According to the flight control requirements, the eight-rotor UAV control system for the discrete state of the collection is , where is the takeoff state, is the normal flight state, is the stability of the flight state, and is out-of-control state. Let the discrete event corresponding to the set of discrete states be . In the formula, is the condition of the switching event that the system changes from continuous state to continuous state ; is expressed as the switching event condition of the system from continuous state to continuous state .The discrete state transition of the hybrid control system is shown in Figure 4.

Figure 4: Discrete state transition diagram of hybrid control system.

When the UAV is in the takeoff or landing state, there is a need for the controller to have a fast response and small overshoot features. When the UAV is in the gust, the controller is needed to have good robustness and a soft control system. When the tilt angle of the UAV is greater than the out-of-control set threshold, the controller switches to the runaway state to prevent a crash occurring. Compared with a single control system, a hybrid control system can provide improved flight status, enhanced the UAV flight stability, and suited to the current flight environment better.

3.1. Backstepping Controller Design

Since the eight-rotor unmanned aerial vehicle (UAV) is an underactuated system that controls the six degrees of freedom of space by adjusting the rotational speed of the eight rotors, there is a nonlinear relationship between the degrees of freedom of the direct excitation portion and the degrees of freedom of the underactuated portion. The backstepping control method is one of the commonly used methods for designing nonlinear system controllers. It is a linear regression design method combining Lyapunov function and controller design. The main design idea is to decompose the complicated system into several subsystems which do not exceed the order of the system. The Lyapunov function and the intermediate virtual control quantity are designed for each subsystem, respectively, until the system stability state is reached, thus completing the entire system controller design [17]. The underactuated system controller designed by backstepping method has the advantages of a short adjustment time and a small overshoot and has been widely used in aircraft, missile, and robot system [18].

The system state of UAV is defined as follows:

The system variables are divided into six control channels: roll channel, pitch channel, yaw channel, -axis channel, -axis channel, and -axis channel, which are defined by the literature [19]. The variables , , , , , and , respectively, are as follows:

The backstepping controller of the rolling channel is

The pitch channel backstepping controller is

The yaw channel backstepping controller is

The height channel backstepping controller is

The -axis channel backstepping controller is

The -axis channel backstepping controller is

where , , , , , , and to are more than 0 constant.

3.2. PID Controller Design

PID control is one of the most widely used control laws in engineering practice. It is a control rate that includes proportional, integral, and differential control. Its advantage is that only three parameter values need to be set and properly debugged to achieve good control effect, especially in the absence of accurate mathematical model and model error, PID control more practical. The control rate for each channel of the eight-rotor UAV is designed as follows:

The rolling channel PID controller is

The pitch channel PID controller is

The yaw channel PID controller is

The height channel PID controller is

The -axis channel PID controller is

The -axis channel PID controller is

3.3. System Stability Analysis

Taking the rolling channel as an example, in the controller stability analysis, the state variable is set as , assuming that the expected roll angle is bounded and changes in the range of ; thus, the roll angle signal tracking error must be bounded. The boundedness of can obtain that in the backstepping control method where the virtual control signal is bounded, the roll angular velocity error is bounded, and the roll control channel input is bounded. Therefore, all signals within the roll channel are bounded.

By defining the amount of error change as , we can derive

Among which, is the system error, and is the amount of virtual control to be determined.

The Lyapunov function is chosen as follows: , where we can derive by taking it into (20) and get

Suppose is the system parameter; in order to meet the stability requirements of the closed loop system, eliminating the coupling term , the closed loop of (20) can be obtained as follows:

The remaining variables are as follows:

For the guidance, is presented as follows:

where is expressed as a systematic error and as the dummy variable to be determined.

The Lyapunov function is chosen as , where we can derive by taking it into (24) and get

Equations (21) and (25) are obtained by the virtual control formula (20).

Tracking error index is gradually stabilized. Substituting the resulting control amount into (25) yields

The error index is asymptotically stable and the tracking error is different from that of the differential homeomorphism . The positive and negative definite of the function are given. According to the Lyapunov stability principle, the eight-rotor UAV backstepping control method is progressively stable. Similarly, other control channels can be seen as gradually stable; we can see that the various subsystems in the hybrid control system are stable.

Consider the model of the general switched linear system as , where , , and is switching signal, and represents the number of subsystems: is switching gain and .

Introduce the following variables:

satisfies the conditions:

If the system switches to subsystem, then . If the system switches to subsystem, then ; it means that there is only one controller works at any time. Then, and .

Hence, satisfying , where is the identity matrix with suitable dimension, the switched system can be equivalent to , adopting the Common Lyapunov function to study the system stability. Some lemmas are required to prove the main results.

Lemma 1 (see [20]). , , and are real matrices with appropriate dimension; if they satisfy , then there always exists .

Lemma 2 (see [19]). (Schur’s lemma) For the given symmetric matrix,

where and are square matrices, and then there exists the following equivalent conditions:

Theorem 1. The switched system is said to be asymptotically stable if there exists positive definite matrix to make

Proof 1. Choosing the Lyapunov function , then we have

by taking a derivative with it.
According to , applying Lemma 1, we have
Then,

It is seen that the system is asymptotically stable, if it satisfies . According to Lemma 2, it is equivalent to

This completes the proof.

Taking the pitching motion of eight-rotor aircraft, we construct its linearized mode as , where and , according to the small disturbance linearization theory in literature [21]. Taking the related parameters of the aircraft, we have

Assume the switching gain as and ; by calculating the linear matrix inequalities, we get

is a positive definite matrix; satisfying the theorem condition, it is obtained that the switched hybrid system is asymptotically stable.

The reason for the instability of the system is that the energy generated by the system is not absorbed during switching, so that the “slow switch” is carried out on the stable subsystem, that is, the average stay of the switching is sufficiently large so that the energy generated at the time of switching is sufficiently absorbed, which can guarantee the stability of the system [22].

4. System Simulation and Analysis

Based on the model of the eight-rotor UAV, the related literature and the relevant parameters of the eight-rotor UAV are shown in Table 1.

Table 1: Main parameters of UAV model.

Parameter
Value
Parameter
Value

Quality ()
1.9
Rotor’s length ()
0.32

Lift coefficient
3.13
Resistance coefficient
7.5

-axis rotary inertia
0.018
-axis rotary inertia
0.018

-axis rotary inertia
0.036
Gravitational acceleration ()
9.8

Based on the above design process and parameters, the MATLAB/Simulink simulation model is established. Taking the pitch channel as an example, the simulation model of the hybrid control system is shown in Figure 5. The simulation model is mainly composed of UAV dynamics model, switcher, controller, and state flow module. Among them, the UAV kinematics module mainly describes the aerodynamic characteristics of the UAV. The switch makes the corresponding controller switch based on the output signal of the state flow. The controller is mainly composed of three parts, such as the PID controller, backstepping controller, and runaway protection module for controlling the attitude and position of the unmanned aerial vehicle so that the UAV can reach the desired state from the initial state. The function of the state flow module is to calculate the number of the selected controller according to the input signal. The PID controller parameters in the figure are , , and ; the backstepping controller parameters are , , , , and .

Figure 5: Pitching channel simulation model.

Hybrid controller performance will be verified by three experiments: fixed point hovering, expected tracking, and antijamming capability test. Due to the complexity of aerodynamic effects and large influence of external conditions, precise aerodynamic analysis of UAV is difficult to model accurately, so the UAV parameters are uncertain. Based on the idea of hybrid control system, the backstepping control method is used to complete the trajectory tracking control of the UAV. The backstepping control method adopts the reverse design method. The Lyapunov function is used to derive the corresponding control rate, which can control the nonlinear system of the -order to make the design process more structured and systematic. However, the method is more dependent on the model aerodynamic parameters, and the change of factors such as external disturbance is slower. So, the backstepping control method is taken as the UAV hybrid control system with the normal flight state controller. In the case of system perturbation, the PID control method is used to suppress the influence of disturbance and improve the control precision. Therefore, the PID control method is used as the control system of the hybrid control system and the uncontrolled state controller to achieve the precise control of the purpose. In the analysis of antijamming capability, random white noise is used as the simulated gust interference source, and the square wave signal is used to simulate the obstacle information. When the system has reached the steady state, it is applied to the corresponding system feedback state variable, and then the effect of the interference information for UVA stability is analyzed. The UAV spatial displacement trajectory simulation results are shown in Figure 6, pitch angle control simulation results are shown in Figure 7, obstacle avoidance control simulation results are shown in Figure 8, and the height control simulation results are shown in Figure 9.

Figure 6: Eight-rotor UAV space displacement trajectory map.

Figure 7: Simulation results of pitch angle control.

Figure 8: Simulation results of obstacle avoidance control.

Figure 9: High control simulation results.

In Figure 6, the initial position coordinates of the spatial displacement trajectory of the eight-rotor UAV are set as . The initial attitude angle is , the expected position coordinate is , the expected attitude angle is , the simulation time is 10 s, and the calculation method is ode4. It can be seen from the figure that the overshoot of the system is small, although there are some errors due to the influence of interference; the trajectory of the response characteristic curve is smooth, the control output distortion is small, and the control of the eight-rotor unmanned aerial vehicle has better control effect in line with the actual flight needs. In Figure 7, the system reaches the expected value of the pitch angle at 0.7 s, the system response time is faster, the rising trend is smooth, and there is no overshoot. White noise interference signals are added between 2 s and 4 s, and the system is switched to PID controller to reach the desired target stability value again at 4.7 s. Besides, the range of interference effect changed slightly, and the steady-state error is almost 0.

In Figure 8, the detection result signal of the UAV external sensor obstacle is simulated by the square wave signal, and the detection result is used as the switching signal of the hybrid control system. Taking the height obstacle as an example, the simulation results show that the square wave interference signal appears at 3.5 s to 4.5 s and 5.5 s to 6.5 s. Based on the preset value of the obstacle sensor detection result, the system will make the switching action in advance to realize the smooth crossing of the obstacle. The simulation results show that the state transition when the unmanned aerial vehicle adopts the obstacle avoidance behavior shows good continuity and stability, the flight trajectory is stable, and the stability control of the eight-rotor UAV is realized.

As can be seen from Figure 9, the system achieves a high expected value at 0.6 s, and the system response time is fast and no overshoot. Adding a white noise interference signal in 2 s to 4 s time, the system switches to the PID controller and reaches a stable state within 5.2 s. The degree that the system affected by the interference amplitude is small and the interference output response curve is smooth. System is affected by the interference signal in 7 s, the rotor motor began to slow down, and the UAV began landing, at 7.9 s. Due to ground gravity and reaction force, there is a slight “reverse” overshoot phenomenon. But the curve of UAV Landing characteristics is smooth, with slightly longer landing time. The advantage is not to cause injuries and other incidents of UAV.

In order to verify the actual performance of hybrid controller, a physical test platform is set up as shown in Figure 10. The platform takes STM32F103 as the main controller chip and MPU6050 chip as the inertial measurement unit to real-time measurement of eight-rotor UAV attitude data. The main controller chip calculates the acceleration and angular velocity data collected by the inertial measurement into the attitude angle and then passes through the main controller to generate the corresponding duty cycle PWM wave to control the speed of the relevant brushless DC motor to realize the correlation of the UAV action. The experimental results show that the hybrid controller can effectively and stably control the roll, pitch, and yaw of the eight-rotor UAV, and it is reliable and satisfies the requirement of attitude control.

Figure 10: Eight-rotor UAV physical map.

5. Conclusion

The continuous development of space technology and the maturation of eight-rotor UAV have led to the availability of low-cost, high-performance, uniquely structured system that is widely used in military and civilian areas. In this paper, the eight-rotor UAV is used as the research object, and the dynamic modeling of the system is expounded in detail. The hybrid control system is based on the backstepping control method and PID control of the eight-rotor UAV in different environments or states. The simulation model of the hybrid control system is built in a MATLAB/Simulink environment. The spatial displacement trajectory of the UAV and the switching process of the pitch angle and the height control channel of the UAV under the interference signal are simulated. The simulation results show that the hybrid control system can choose the appropriate control method according to the change of the environment and better adapt to the current flight environment. Finally, by establishing the physical test platform of the eight-rotor UAV, the reliability, stability, and anti-interference ability of the hybrid control system are verified, and the control requirements of flight attitude and position are satisfied. The stability control of the drone in a complex environment is realized, and the flight stability is improved.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

References

H. Bolandi, M. Rezaei, R. Mohsenipour, H. Nemati, and S. M. Smailzadeh, “Attitude control of a quad-rotor with optimized PID controller,” Intelligent Control and Automation, vol. 4, no. 3, pp. 335–342, 2013.

Q. Hu, Q. Fei, and Q. H. Wu, “Research and application of nonlinear control techniques for quad-rotor UAV,” Journal of University of Science and Technology of China, vol. 42, no. 8, pp. 656–663, 2012.

A. MILHIM and Y. M. ZHANG, “Gain scheduling based PID controller for fault tolerant control of quad-rotor UAV,” in AIAA Infotech@Aerospace 2010, pp. 1–13, Atlanta, GA, USA, April 2010.

L. Zhang and H. Li, “Attitude control of four-rotor aircraft via fuzzy PID,” Computer Simulation, vol. 31, no. 8, pp. 73–77, 2014.

T. Dierks and S. Jagannathan, “Neural network output feedback control of a quadrotor UAV,” in 2008 47th IEEE Conference on Decision and Control, pp. 3633–3639, Cancun, Mexico, December 2008.

B. Sumantri, N. Uchiyama, S. Sano, and Y. Kawabata, “Robust tracking control of a quad-rotor helicopter utilizing sliding mode control with a nonlinear sliding surface,” Journal of System Design and Dynamics, vol. 7, no. 2, pp. 226–241, 2013.

L. I. U. Huanye, L. I. Jian, Y. A. O. Jianguo et al., “Backstepping based adaptive control for a mini rotorcraft with four rotors,” in 2010 Second International Conference on Computer Modeling and Simulation, pp. 472–476, Sanya, Hainan, China, January 2010.

Y. Wang, Q. Wu, and Y. Wang, “Distributed cooperative control for multiple quad-rotor systems via dynamic surface control,” Nonlinear Dynamics, vol. 75, no. 3, pp. 513–527, 2014.

Y. N. Li, T. I. Li, Y. Jiang, and J. I. Fan, “Adaptive PID control of quadrotor based on RBF neural network,” Control Engineering of China, vol. 23, no. 3, pp. 378–382, 2016.

Z. Ali, D. Wang, and M. Aamir, “Fuzzy-based hybrid control algorithm for the stabilization of a trirotor UAV,” Sensors, vol. 16, no. 5, p. 652, 2016.

T. P. Hoang, T. D. Chi, B. P. Thanh, and V. T. Nguyen, “Hybrid terminal sliding mode control and quadrotor’s vision based ground object tracking,” in 2013 International Conference on Control, Automation and Information Sciences (ICCAIS), pp. 334–339, Nha Trang, Vietnam, November 2013.

G. Xia, Y. Liao, and L. Wang, “Research on hybrid control system of quadrotor UAV,” Journal of Computer Applications, vol. 33, no. 3, pp. 858–861, 2013.

D. Mellinger, “Recent advances in quadrotor capabilities,” in 2011 IEEE International Conference on Robotics and Automation, pp. 2964-2965, Shanghai, China, May 2011.

L. D. Minh, “Modeling and control of quadrotor MAV using vision based measurement,” IEEE Transactions on Circuits and Systems, vol. 33, no. 4, pp. 70–77, 2010.

B. Lennartson, M. Tittus, B. Egardt, and S. Pettersson, “Hybrid systems in process control,” IEEE Control Systems Magazine, vol. 16, no. 5, pp. 45–56, 1996.

Q. Geng, H. Shuai, and Q. Hu, “Obstacle avoidance approaches for quad-rotor UAV based on backstepping technique,” in 2013 25th Chinese Control and Decision Conference (CCDC), pp. 3613–3617, Guiyang, China, May 2013.

C. Ha, Z. Zuo, F. B. Choi, and D. Lee, “Passivity-based adaptive backstepping control of quadrotor-type UAVs,” Robotics and Autonomous Systems, vol. 62, no. 9, pp. 1305–1315, 2014.

H. T. Zhen, X. H. Qi, and M. Q. Xia, “Block backstepping attitude controller design of quad-rotor UAV,” Electronics Optics & Control, vol. 20, no. 10, pp. 87–91, 2013.

A. DAS, F. LEWIS, and K. SUBBARAO, “Backstepping approach for controlling a quad-rotor using Lagrange form dynamics,” Journal of Intelligent and Robotic Systems, vol. 56, no. 1-2, pp. 127–151, 2009.

L. Yu, Robust Control-Linear Matrix Inequality Processing Method, Tsinghua University Press, Beijing, China, 2002.

M. Q. Zhang, Flight Control System, Aviation Industry Press, Beijing, China, 1994.

G. Zhai, B. Hu, K. Yasuda, and A. N. Michel, “Disturbance attenuation properties of time-controlled switched systems,” Journal of the Franklin Institute, vol. 338, no. 7, pp. 765–779, 2001.