Research Article | Open Access

Naibao He, Changsheng Jiang, Bin Jiang, Qian Gao, "Terminal Sliding Mode Control with Unidirectional Auxiliary Surfaces for Hypersonic Vehicles Based on Adaptive Disturbance Observer", *Mathematical Problems in Engineering*, vol. 2015, Article ID 285716, 11 pages, 2015. https://doi.org/10.1155/2015/285716

# Terminal Sliding Mode Control with Unidirectional Auxiliary Surfaces for Hypersonic Vehicles Based on Adaptive Disturbance Observer

**Academic Editor:**Carlo Cosentino

#### Abstract

A novel flight control scheme is proposed using the terminal sliding mode technique, unidirectional auxiliary surfaces and the disturbance observer model. These proposed dynamic attitude control systems can improve control performance of hypersonic vehicles despite uncertainties and external disturbances. The terminal attractor is employed to improve the convergence rate associated with the critical damping characteristics problem noted in short-period motions of hypersonic vehicles. The proposed robust attitude control scheme uses a dynamic terminal sliding mode with unidirectional auxiliary surfaces. The nonlinear disturbance observer is designed to estimate system uncertainties and external disturbances. The output of the disturbance observer aids the robust adaptive control scheme and improves robust attitude control performance. Finally, simulation results are presented to illustrate the effectiveness of the proposed terminal sliding mode with unidirectional auxiliary surfaces.

#### 1. Introduction

It is well known that the near space hypersonic vehicle (NHV) as a novel aerospace aircraft plays an important role in military and civil applications [1]. Compared with traditional flight vehicles, the NHV has significant difference in flight airspace, working states, and mission patterns. In addition, it possesses many advantages over traditional spacecraft, such as lower launch cost, maintainability, speed of redeployment, and reusability. Consequently, the NHV concept has garnered considerable interest over the last few years. On the other hand, the NHV dynamics is severely nonlinear, time-varying, highly uncertain, and strongly coupled. It also suffers different external disturbances with changes in the flight environment, propulsion system perturbation, atmospheric disturbances, and other factors. Therefore, robust flight control of the NHV is a challenging problem. Many advanced control schemes have been presented in the literatures for control of NHVs. For example, Chen and Jiang [2] developed robust attitude control for NHVs with time-varying disturbances. Zhao et al. [3, 4] proposed a fault-tolerant control scheme based on a second-order dynamic terminal sliding mode. Their model calls for a near space vehicle (NSV) attitude dynamic system in the presence of model parameter uncertainties and external disturbances. He et al. [5] designed an adaptive output feedback fault-tolerant controller for the longitudinal dynamics of a generic hypersonic flight vehicle. The robust sliding mode control strategy was applied to the NSV attitude dynamics. In [6], an adaptive nonsingular terminal sliding mode approach was presented for the attitude control of NHV in the presence of parameter uncertainties and external disturbances. In [7], a second-order dynamic terminal sliding mode control was proposed for the NSV with high speed, violent changes of aerodynamic parameters, and large external disturbance. The adaptive attitude control was studied for the NSV with dynamic uncertainties using a functional link network (FLN) in [8]. In [9, 10], an adaptive tracking control was proposed for the NSV using Takagi-Sugeno fuzzy models. In [11, 12], a robust back-stepping control was studied for a flexible air-breathing hypersonic vehicle in the presence of input constraint and aerodynamic uncertainty. The disturbance observer-based control method was proposed for uncertain flight dynamics system in [13–16]. In the above-mentioned control approaches, sliding mode control (SMC) has been extensively studied for handling systems with large uncertainties, nonlinearities, and bounded external disturbances. Due to remarkable robust control performance, the SMC has been widely applied in many practical systems [17–23].

In order to alleviate chattering, a new approach called sliding mode control with unidirectional auxiliary surfaces (UAS-SMC) has been proposed to constrain the system states [24–27]. The unidirectional auxiliary surfaces in UAS-SMC are used to construct a positive invariant set. It implies that once the system states come into this positive invariant set, they would never come out again. However, the UAS-SMC scheme has a specific disadvantage. Its drawback is that the convergence speed of UAS-SMC is slower than that of the conventional SMC. This drawback prevents the UAS-SMC strategy from being applied to some high-frequency systems such as the angular rate loop of flight control systems.

Based on the above analysis, the key breakthroughs can be concluded as follows: The aim of this paper is to first propose UAS-TSMC with suitably high convergence speed. The linear switching surface in UAS-SMC strategy is replaced by a terminal switching surface in UAS-TSMC to improve the convergence speed. And the nonlinear unidirectional auxiliary surfaces in UAS-TSMC can be proven to stay in a positive invariant set. Secondly, by combining the nonlinear disturbance observer with the adaptive terminal sliding control via unidirectional auxiliary surfaces (UAS-TSMC-NDO), a composite UAS-TSMC-NDO method is constructed to handle the external disturbances, and it can make the system more stable, alleviate chattering, and strengthen robustness.

This work is motivated by the above discussion. A novel UAS-TSMC-NDO control scheme is proposed for the NHV, which has been applied to the tracking control of attitude angles. The organization of the paper is as follows: Section 2 details the problem statement, and Section 3 gives the design of terminal sliding mode control with unidirectional auxiliary surfaces based on the disturbance observer. Simulation results are presented to demonstrate the effectiveness of our proposed approaches given in Section 4, followed by the conclusion remarks in Section 5.

#### 2. Problem Statement

The attitude mathematic model for hypersonic vehicles developed at NASA Langley Research Center is given by [28]where is the angle of attack, is the sideslip angle, is the bank angle, is the roll rate, is the pitch rate, and is the yaw rate.

According to the singular perturbation theory, the six equations can be divided into the fast loop and the slow loop, respectively. In order to study robust attitude control for the NHV, the kinematic and dynamic attitudes of the hypersonic vehicle are described as the following multi-input and multioutput (MIMO) affine nonlinear system [29]:where is the slow loop state, the attitude angle vector. is the fast loop state, the body-axis angular rate vector. is the compounded external disturbance. is the system control input and is the system output. The detailed expressions of , , , and can be found in [29].

In this paper, the attitude control objective is that the robust control schemes are designed to track a given desired output of the attitude motion. In the attitude control design, we definewhere is the attitude angle tracking error vector, is the attitude angular rate tracking error vector, and is a command signal which will be designed.

In order to proceed with the design of robust attitude control scheme for attitude motions (2a), (2b), and (2c) of the NHV, the following assumptions are required.

*Assumption 1. *For the compounded disturbance , there exists an unknown positive constant such that .

*Assumption 2. *For all , there exist such that , .

#### 3. Design of Terminal Sliding Mode Control with Unidirectional Auxiliary Surfaces

In this section, a terminal sliding mode control with unidirectional auxiliary is developed. First, a terminal sliding mode control algorithm is designedto stabilize the slow loop system and ensure that the unidirectional auxiliary sliding mode converges to the equilibrium surfaces asymptotically. Then, the robust adaptive control will be designed based on disturbance observer that guarantees the fast loop system will remain stable in the presence of uncertainties and external disturbances effecting the NHV.

##### 3.1. Slow Loop Terminal Sliding Mode Controller Design

Considering the slow loop dynamic (3) of the NHV and differentiating yieldedThe detailed design of UAS-TSMC is given as follows.

*Step 1. *The terminal switching surfaces for system (5) are given bywhere , are the designed constants and and are positive odd integers, which satisfy the following condition:Defining , , and , the terminal switching surfaces in (6) can be rewritten as

*Step 2. *Based on the terminal switching surface , , number subspace can be defined as in Figure 1.

As shown in Figure 1, the appropriate points , , , and should be selected on terminal switching surfaces , , where points , are located in the fourth quadrant and points , are located in the second quadrant. The lines , , , and defined in Figure 2 are UAS , , , and , respectively. For the simplicity of discussion, we define that the UAS , , is designed by points , where , . The expression of UAS is given as follows:where , , , .

*Remark 3. *These UAS , , , should constitute a convex quadrilateral which is given byHence, the system states in set imply that the system states satisfy the state constraints. This convex quadrilateral can be proved to be a positively invariant set [30] with control equation (9).

*Step 3. *The UAS in Figure 1 would be used to design control input when the states are moving in number subspace. The compact form of current UAS is given by where

*Step 4. *The UAS-TSMC control law for system (5) is designed aswhere is the approaching law.

*Remark 4. * is a nonsingular sliding mode surface, which does not produce singularities. But it may cause singularity when we compute the time derivative of . Therefore, in the design of the controller, we add the saturation function in order to overcome the singularity problem [31].

*Definition 5 (see [30]). *The set is said to be positively invariant for system (5) with control law in (13). If for all the solution for . If implies for all , we say that is invariant.

*Remark 6. *As shown in Figure 2, the set is a convex set. For the point , if the point locates in number subspace where , there exists .

Lemma 7. *For the point in number subspace where , the UAS in No. subspace is designed by points , where and . The expression of UAS is written aswhere . Then, for point , there exists the fact that and when .*

*Proof. *Define the coordinates of points and as follows:As shown in Figure 2, it is noted that, for the point in number subspace, there exists , , . Then, the coordinate of point can be expressed as follows:In (14) with (16), one hasFrom the condition , , we haveSubstituting (18) into (17), we obtainHence, the first conclusion is as follows:.

When , there exists from (19). Since , and , we can conclude , . Then, we can obtain the second conclusion as follows:This completes the proof.

Theorem 8. *Considering system (5) with the control law in (13), system (5) is asymptotically stable. And the set can be guaranteed as a positively invariant set [30].*

*Proof. *See the Appendix.

##### 3.2. Fast Loop Terminal Sliding Mode Disturbance Observer Design

In this section, a robust flight control approach will be developed based on the combination of nonlinear disturbance observer and terminal sliding mode control with unidirectional auxiliary surfaces. For control of external disturbances and model uncertainties in the fast loop, they cannot be directly used to design the robust sliding mode controller. To efficiently handle them, the nonlinear disturbance observer is proposed to estimate them.

In order to estimate the external disturbance of the near space hypersonic vehicle system, the general nonlinear near space hypersonic vehicle fast loop dynamic error model in (4) will be used. Considering (4) and differentiating yielded

The disturbance observer is proposed as follows [32]:where is an auxiliary variable of the proposed disturbance observer, is the output of the nonlinear disturbance observer, , and is an appropriate dimension gain matrix which is given byDefine the estimate error . Since there is no prior information about the derivative of the disturbance , it is reasonable to suppose that , which implies that the disturbance varies slowly relative to the observer dynamics.

Differentiating , we have

Considering (18) and differentiating with respect to time yieldedThus, we have . It is apparent that the choice of must make the disturbance observer error converge to zero.

According to the design of terminal sliding mode with unidirectional auxiliary surfaces, the robust adaptive control of the hypersonic vehicle based on estimated disturbance observation can be surmised in the following theorem.

Theorem 9. *For NHV dynamic error model (4), the disturbance observer is designed as (22), and there exists a terminal sliding mode with unidirectional auxiliary surfaces (26), under the two terminal switching surfaces (28), (29), the adaptive law (27) and the NDO-based robust flight control law (30). Then the robust stability of the control system in the presence of the parameter uncertainties and external disturbances is guaranteed:where, , , , are designed constants and and are positive odd integers and satisfy .*

*Proof. *Consider the following Lyapunov function:Differentiating with respect to time yieldsThus, the proof is completed.

#### 4. Simulation Study

To illustrate the effectiveness of the proposed robust control scheme, the simulations are given for the near space hypersonic vehicles. The initial conditions are chosen as follows: kg, m/s, m, m, km, , , the attitude angle deg, deg, deg, the attitude angle rate deg/s, deg/s, and deg/s. In order to guarantee the stability of the near space hypersonic vehicles attitude, the reference signals of the attitude angle are chosen as deg, deg, and deg. The disturbances are considered as follows:

For the purpose of comparison for the effectiveness of TSMC between the proposed approach presented in this paper and an existing method, the corresponding simulation results are depicted in Figures 3–5. From Figure 3, it is shown that control torques responses using TSMC method cannot meet system performance requirements. From Figure 4, under the designed UAS-TSMC attitude control scheme, attitude angle rate and attitude angle become unstable quickly. The result shows that there exists seriously error during most of the simulation time. Furthermore, in Figure 5, the UAS-TSMC with NDO algorithm can make the system stable and indicates that all of the responses of attitude angle and angular velocity can track the desired objectives well. Thus, the robust flight control approach based on the UAS-TSMC with NDO can satisfy precision and robustness requirements in the presence of uncertainties and large external disturbances.

#### 5. Conclusion

In this paper, a novel robust sliding mode control with unidirectional auxiliary surface approach has been developed to solve the tracking control problem for the NHV with functional uncertainty and external disturbance. UAS-TSMC is employed to increase the robustness of the control system and improves the control performance. Meanwhile, the nonlinear disturbance observer has been developed to estimate uncertainties and external disturbances. Finally, the proposed method has been applied to near space hypersonic vehicles attitude dynamics, and simulation results have been given to illustrate the effectiveness of the proposed control approach. Studying attitude tracking control and accommodation approaches for NHVs with functional uncertainty and external disturbance will be investigated in our future work.

#### Appendix

*Proof of Theorem 8. *(1) The proof of the asymptotical stability for system (5) is given in this part.

For system (5), we consider the following Lyapunov function candidate:Firstly, it would be proved that and there exists when . From the definition of , it follows that and there exists when . Based on the discussion given in Lemma 7, it is noted that when is located in number subspace where . Then, the condition implies . Hence, there exists when .

Secondly, the conclusion when is given as follows.

It follows from (11) and (13) that the time derivative of is given asSince , the derivative of should be a positive value.

From the Lyapunov function in (A.1), the derivative of is expressed asBased on Lemma 7, the conclusion is obtained that when is located in number subspace and . As when , it implies that when . Then, there exists the fact that when . From and the positive derivative of , the derivative of should be a negative value when .

Thirdly, the proof for the continuity of Lyapunov function in (A.1) is given in this part.

As shown in Figure 6, the point on switching surface is selected to show the continuity of Lyapunov function in (A.1) when the point is switched between number and , , , subspace.

The point is one of the points which is used to design the UAS and . Define the coordinates of point as . From (18), the expressions of UAS and can be given as follows:Since the point is located on the UAS and , (A.5) can be written asFor the point in number subspace, there exists , , as given in Figure 6. By substituting the coordinate of point into the UAS , we haveClearly, there exists by considering (A.5) and (A.6).

Then, the Lyapunov function in number subspace can be rewritten as follows:For the point in number subspace, there exists , , as given in Figure 6. By substituting the coordinate of point into the UAS , it is obtained thatClearly, there exists by considering (A.5) and (A.8).

Then the Lyapunov function in number subspace can be rewritten as follows:Hence, when the point is switched between number and subspaces, there exists . Then, the continuity of Lyapunov function in (A.1) is proved.

From the above demonstration, it is noted that the Lyapunov function in (A.1) is continuous; and there exists when ; when . Hence, system (5) with the control law in (13) is asymptotically stable and converges to . And the state is defined asIt is clear that when . Thus, we obtain that system (5) converges to .

Demonstration for the positively invariant set is given in this part. This conclusion is supported with the reduction to absurdity.

Considering the initial point in number subspace, we suppose that there exists in number subspace, where . Since point in number subspace and , it is noted that the Lyapunov function in (A.1) can be expressed as follows:Since point in number subspace, there exists from Remark 6. Clearly, the conclusion in (A.2) is obtained asIt follows from (A.11) and (A.12) that the continuous Lyapunov function has the result , . Hence, there should exist a point , , where . However, from the above discussion in Part (1), it is noted that . Then, we find the contradiction since the derivative of Lyapunov function cannot be positive and negative at the same time on the point . Hence, the assumption that there exists is incorrect, and implies for all . From Definition 5, the set is a positively invariant set. This completes the proof.

#### Conflict of Interests

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

#### Acknowledgments

The authors wish to thank the editor and anonymous reviewers for useful comments and suggestions. This work is partially supported by JSUT Research Funding (Grant nos. KYY13001, KYY13017) and Innovation Team Funding (Grant no. TDZD13003).

#### References

- M. Young, S. Keith, and A. Pancotti, “An overview of advanced concepts for near-space systems,” in
*Proceedings of the 45th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit*, pp. 1–8, August 2009. View at: Google Scholar - M. Chen and B. Jiang, “Robust attitude control of near space vehicles with time-varying disturbances,”
*International Journal of Control, Automation and Systems*, vol. 11, no. 1, pp. 182–187, 2013. View at: Publisher Site | Google Scholar - J. Zhao, B. Jiang, P. Shi, Z. Gao, and D. Xu, “Fault-tolerant control design for near-space vehicles based on a dynamic terminal sliding mode technique,”
*Proceedings of the Institution of Mechanical Engineers. Part I: Journal of Systems and Control Engineering*, vol. 226, no. 6, pp. 787–794, 2012. View at: Publisher Site | Google Scholar - J. Zhao, B. Jiang, P. Shi, and Z. He, “Fault tolerant control for damaged aircraft based on sliding mode control scheme,”
*International Journal of Innovative Computing, Information and Control*, vol. 10, no. 1, pp. 293–302, 2014. View at: Google Scholar - J. He, R. Qi, B. Jiang, and J. Qian, “Adaptive output feedback fault-tolerant control design for hypersonic flight vehicles,”
*Journal of the Franklin Institute. Engineering and Applied Mathematics*, vol. 352, no. 5, pp. 1811–1835, 2015. View at: Publisher Site | Google Scholar | MathSciNet - R. Zhang, L. Dong, and C. Sun, “Adaptive nonsingular terminal sliding mode control design for near space hypersonic vehicles,”
*IEEE/CAA Journal of Automatica Sinica*, vol. 1, no. 2, pp. 155–161, 2014. View at: Publisher Site | Google Scholar - M. Pu, Q. X. Wu, C. S. Jiang, and J. Zhang, “Near space vehicle control based on second-order dynamic terminal sliding mode,”
*The Journal of Astronautics*, vol. 31, no. 4, pp. 1056–1062, 2010. View at: Publisher Site | Google Scholar - Y. L. Du, Q. X. Wu, C. S. Jiang, and J. Wen, “Adaptive functional link network control of near-space vehicles with dynamical uncertainties,”
*Journal of Systems Engineering and Electronics*, vol. 21, no. 5, pp. 868–876, 2010. View at: Publisher Site | Google Scholar - D. X. Gao and Z. Q. Sun, “Fuzzy tracking control design for hypersonic vehicles via T-S model,”
*Science China Information Sciences*, vol. 54, no. 3, pp. 521–528, 2011. View at: Publisher Site | Google Scholar | MathSciNet - B. Jiang, Z. Gao, P. Shi, and Y. Xu, “Adaptive fault-tolerant tracking control of near-space vehicle using Takagi-Sugeno fuzzy models,”
*IEEE Transactions on Fuzzy Systems*, vol. 18, no. 5, pp. 1000–1007, 2010. View at: Publisher Site | Google Scholar - F. Wang, C. Hua, and Q. Zong, “Attitude control of reusable launch vehicle in reentry phase with input constraint via robust adaptive backstepping control,”
*International Journal of Adaptive Control and Signal Processing*, 2015. View at: Publisher Site | Google Scholar - Q. Zong, F. Wang, B. Tian, and R. Su, “Robust adaptive approximate backstepping control of a flexible air-breathing hypersonic vehicle with input constraint and uncertainty,”
*Proceedings of the Institution of Mechanical Engineers. Part I: Journal of Systems and Control Engineering*, vol. 228, no. 7, pp. 521–539, 2014. View at: Publisher Site | Google Scholar - M. Chen and B. Jiang, “Robust bounded control for uncertain flight dynamics using disturbance observer,”
*Journal of Systems Engineering and Electronics*, vol. 25, no. 4, Article ID 6905958, pp. 640–647, 2014. View at: Publisher Site | Google Scholar - X. Wei and L. Guo, “Composite disturbance-observer-based control and ${H}_{\infty}$ control for complex continuous models,”
*International Journal of Robust and Nonlinear Control*, vol. 20, no. 1, pp. 106–118, 2010. View at: Google Scholar - W. H. Chen, “Nonlinear disturbance observer-enhanced dynamic inversion control of missiles,”
*Journal of Guidance, Control, and Dynamics*, vol. 26, no. 1, pp. 161–166, 2003. View at: Publisher Site | Google Scholar - F. Chen, F. Lu, B. Jiang, and G. Tao, “Adaptive compensation control of the quadrotor helicopter using quantum information technology and disturbance observer,”
*Journal of the Franklin Institute*, vol. 351, no. 1, pp. 442–455, 2014. View at: Publisher Site | Google Scholar - L. Wu, P. Shi, and X. Su,
*Sliding Mode Control of Uncertain Parameter-Switching Hybrid Systems*, John Wiley & Sons, Chichester, UK, 2014. View at: Publisher Site - Y. Wang, Q. Zhang, and L. Zhu, “Active control of hypersonic airfoil flutter via adaptive fuzzy sliding mode method,”
*Journal of Vibration and Control*, vol. 21, no. 1, pp. 134–141, 2015. View at: Publisher Site | Google Scholar - L. Yan and B. S. Ju, “PMSM speed controller using switching algorithm of PD and sliding mode control,” in
*Proceedings of the ICROS-SICE International Joint Conference*, pp. 1260–1266, Fukuoka International Congress Center, Fukuoka, Japan, August 2009. View at: Google Scholar - M. S. Chen, Y. R. Hwang, and M. Tomizuka, “A state-dependent boundary layer design for sliding mode control,”
*IEEE Transactions on Automatic Control*, vol. 47, no. 10, pp. 1677–1681, 2002. View at: Publisher Site | Google Scholar | MathSciNet - V. S. Bandal and P. N. Vernekar, “Design of a discrete-time sliding mode controller for a magnetic levitation system using multirate output feedback,” in
*Proceedings of the American Control Conference (ACC '10)*, pp. 4289–4294, Baltimore, Md, USA, July 2010. View at: Google Scholar - R. Zhang, C. Sun, J. Zhang, and Y. Zhou, “Second-order terminal sliding mode control for hypersonic vehicle in cruising flight with sliding mode disturbance observer,”
*Journal of Control Theory and Applications*, vol. 11, no. 2, pp. 299–305, 2013. View at: Publisher Site | Google Scholar - J. Yang, S. Li, C. Sun, and L. Guo, “Nonlinear-disturbance-observer-based robust flight control for airbreathing hypersonic vehicles,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 49, no. 2, pp. 1263–1275, 2013. View at: Publisher Site | Google Scholar - J. Fu, Q. X. Wu, C. S. Jiang, and Y. Wang, “Robust sliding mode positively invariant set control for a class of nonlinear continuous systems,”
*Acta Automatica Sinica*, vol. 37, no. 11, pp. 1395–1401, 2011. View at: Google Scholar | MathSciNet - J. Fu, Q. X. Wu, C. S. Jiang, and L. Cheng, “Robust sliding mode control with unidirectional auxiliary surfaces for a nonlinear system with state constraints,”
*Control and Decision*, vol. 26, no. 9, pp. 1288–1294, 2011. View at: Google Scholar | MathSciNet - W. I. Ai, Y. H. Wang, L. Zhu, and C. Xu, “Design of sliding mode controller with unidirectional auxiliary surfaces based on improved exponent trending law for post-stall maneuver,”
*Journal of Jilin University (Information Science Edition)*, vol. 33, no. 2, pp. 153–160, 2015. View at: Google Scholar - T. W. Hu, Q. X. Wu, C. S. Jiang et al., “Exponential approach law for sliding mode control with unidirectional auxiliary surfaces,”
*Electronics Optics & Control*, vol. 20, no. 7, pp. 36–41, 2013. View at: Google Scholar - S. A. Snell,
*Nonlinear dynamic-inversion flight control of supermaneuverable aircraft [Ph.D. thesis]*, University of Minnesota, 1991. - M. Chen, B. Jiang, Q. X. Wu, and C. Jiang, “Robust control of near-space vehicles with input backlash-like hysteresis,”
*Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering*, vol. 227, no. 8, pp. 635–644, 2013. View at: Publisher Site | Google Scholar - F. Blanchini, “Set invariance in control,”
*Automatica*, vol. 35, no. 11, pp. 1747–1767, 1999. View at: Google Scholar - Y. Feng, X. Yu, and F. Han, “On nonsingular terminal sliding-mode control of nonlinear systems,”
*Automatica*, vol. 49, no. 6, pp. 1715–1722, 2013. View at: Publisher Site | Google Scholar | MathSciNet - W. H. Chen, “Disturbance observer based control for nonlinear systems,”
*IEEE/ASME Transactions on Mechatronics*, vol. 9, no. 4, pp. 706–710, 2004. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2015 Naibao He 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.