Int. J. of Computers, Communications & Control, ISSN 1841-9836, E-ISSN 1841-9844 Vol. VI (2011), No. 3 (September), pp. 418-427 Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control T.C. Lin, C.H. Kuo, V.E. Balas Tsung-Chih Lin Feng-Chia University, 40724, Taichung, Taiwan E-mail: tclin@fcu.edu.tw Chia-Hao Kuo Ph.D Program in Electrical and Communications Engineering Feng-Chia University, Taichung, Taiwan E-mail: peterqo022@hotmail.com Valentina E. Balas Aurel Vlaicu University of Arad, Romania B-dul Revolutiei 77, 310130 Arad, Romania E-mail: balas@drbalas.ro Abstract: In this paper, in order to achieve tracking performance of un- certain fractional order chaotic systems an adaptive hybrid fuzzy controller is proposed. During the design procedure, a hybrid learning algorithm combining sliding mode control and Lyapunov stability criterion is adopted to tune the free parameters on line by output feedback control law and adaptive law. A weighting factor, which can be adjusted by the trade-off between plant knowl- edge and control knowledge, is adopted to sum together the control efforts from indirect adaptive fuzzy controller and direct adaptive fuzzy controller. To confirm effectiveness of the proposed control scheme, the fractional order chaotic response system is fully illustrated to track the trajectory generated from the fractional order chaotic drive system. The numerical results show that tracking error and control effort can be made smaller and the proposed hybrid intelligent control structure is more flexible during the design process. Keywords: Fractional order chaotic systems; fuzzy logic control, adaptive hybrid control. 1 Introduction in domain Due mainly to its demonstrated applications in numerous seemingly diverse and widespread fields of science and engineering, fractional calculus has gained considerable popularity and im- portance during past three decades [1]- [2]. In control system, due to the fact that the theoretical aspects are well established, fractional order controllers are successfully used to enhance the per- formance of the feedback control loop. It is observed that the description of some systems is more accurate when the fractional derivative is used. Nowadays, many fractional-order differ- ential systems behave chaotically, such as the fractional-order Chua’s system [3], the fractional- order Duffing system [4], the fractional-order system, the fractional-order Chen’s system [5], the fractional-order cellular neural network [6], the fractional-order neural network [7]. The tracking problem of fractional order chaotic systems is first investigated by Deng and Li [21] who carried out tracking in case of the two fractional Lü systems. Afterwards, they studied chaos tracking of the Chen system with a fractional order in a different manner [22]-[24]. Based on the universal approximation theorem, [9]- [20] (fuzzy logic controllers are general enough to perform any nonlinear control actions) there is rapidly growing interest in systematic Copyright c⃝ 2006-2011 by CCC Publications Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control 419 design methodologies for a class of nonlinear systems using fuzzy adaptive control schemes. Like the conventional adaptive control, the adaptive fuzzy control is classified into direct and indirect fuzzy adaptive control categories [9], [17]- [19]. A direct adaptive fuzzy controller uses fuzzy logic systems as controller in which linguistic fuzzy control rules can be directly incorporated into the controller. On the other hand, an indirect adaptive fuzzy controller uses fuzzy descriptions to model the plant in which fuzzy IF-THEN rules describing the plant can be directly incorporated into the indirect fuzzy controller. Moreover, a hybrid adaptive fuzzy controller can be constructed using a weighting factor to sum together the control efforts from indirect adaptive fuzzy controller and direct adaptive fuzzy controller. Although the concept of sliding mode control (SMC) and the theory of fractional order system are well known, their integration, fractional sliding mode control, is an interesting filed of research dwelt on this paper with some applications [8]. The motivation of this paper stands on two driving forces: One, most systems in the reality display behavior characterized best in time domain of fractional operators, the other, the uncertainties on the process dynamics can appropriately be alleviated by utilizing SMC technique. In this paper, by combining the approximate mathematical model, linguistic model descrip- tion and linguistic control rules into a single adaptive fuzzy controller, an adaptive hybrid fuzzy controller is proposed to achieve prescribed tracking performance of fractional order chaotic sys- tems. A new adaptive hybrid fuzzy SMC algorithm incorporated Lyapunov stability criterion is proposed so that not only the stability of adaptive fuzzy control system is guaranteed but also the influence of the approximation error and external disturbance on the tracking error can be attenuated to an arbitrarily prescribed level. This paper is organized as follows: In section2, an introduction to fractional derivative and its relation to the approximation solution will be addressed. Section 3 generally proposes adaptive hybrid fuzzy SMC of uncertain fractional order systems in presence of uncertainty and its stability analysis. In Section 4, application of the proposed method on fractional order expression chaotic system is investigated. Finally, the simulation results and conclusion will be presented in Section 5. 2 Basic definition and preliminaries for fractional order systems The concept of fractional calculus is popularly believed to have steamed from a question raised in the year 1695 by Marquis de L’Hoptial to Gottfried Wilhelm Leibniz. It is a generalization of integration and differentiation to non-integer order fundamental operator, denoted by aD q t , where a and t are the limits of the operator. This operator is a notation for taking both the fractional integral and functional derivative in a single expression defined as [1] aD q t = dq dtq , q > 0 1 q = 0∫ a t (dτ)−q, q < 0 (1) There are some basic definitions for the general fractional and the commonly used defini- tions are Grunwald-Letnikov and Riemann-Liouville [1]. The Grunwald-Letnikov definition is expressed as aD q t f(t) = lim h→0 [t−ah ]∑ j=0 (−1)j ( a b ) f(t − jh) (2) 420 T.C. Lin, C.H. Kuo, V.E. Balas where [.] is the integer part. The simplest and easiest definition is Riemann-Liouville defini- tion given as aD q t f(t) = 1 Γ(n − q) dn dtn ∫ t 0 f(τ) (t − τ)q−n+1 dτ (3) where n is the first integer which is not less q, i.e., n − 1 < q < n, and Γ is the Gamma function. The numerical simulation of a fractional differential equation is not simple as that of an ordinary differential equation. In this paper, the algorithm which is an improved version of Adams-Bashforth-Moulton algorithm to find an approximation for fractional order systems based on predictor-correctors is given. Consider the following differential equation aD q t y(t) = r(y(t), t), 0 ≤ t ≤ T and y (k)(0) = y(k)o ,k = 0,1,2, ...,m − 1 (4) where aD q t y(t) = 1 Γ(m−q) ∫ t 0 f(m)(τ) (t−τ)q−m+1 dτ, m − 1 < q < m dm dtm y(t), q = m (5) and m is the first integer larger the q. The solution of the equation (4) is equivalent to Volterra integral equation [1] described as y(t) = [q]−1∑ k=0 y (k) 0 tk k! + 1 Γ(q) ∫ t 0 (t − λ)q−1r(y(λ),λ)dλ (6) Let h=T/N, tn = nh, n=0,1,2,· · · N . Then (6) can be discretized as follows. yh(tn+1) = [q]−1∑ k=0 y (k) 0 tkn+1 k! + hq Γ(q + 2) r(y p h(tn+1), tn+1) + hq Γ(q + 2) n∑ j=0 aj,n+1r(yh(tj), tj) (7) where predict value yph(tn+1) is determined by y p h(tn+1) = [q]−1∑ k=0 y (k) 0 tkn+1 k! + hq Γ(q) n∑ j=0 bj,n+1r(yh(tj), tj) (8) and aj,n+1 = nq+1 − (n − q)(n + 1)q, j = 0 (n − j + 2)q+1 + (n − j)q+1 − 2(n − j + 1)q+1 1 ≤ j ≤ n 1 j = n + 1 (9) bj,n+1 = hq q ((n + 1 − j)q − (n − j)q) (10) The approximation error is given as max j=0,1,2,···N |y(tj) − yh(tj)| = O(hp) (11) where p=min(2,1+q). Therefore, the numerical solution of a fraction order chaotic system discussed in this paper can be obtained by applying the above mentioned algorithm. Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control 421 3 Adaptive hybrid fuzzy sliding mode control of uncertain frac- tional order chaotic systems In this section, we study adaptive hybrid fuzzy tracking control of uncertain fractional order chaotic systems, i.e., to force output trajectory which is obtained by the algorithm mentioned in section 2 of the response system to track output trajectory of the drive system. Consider a fractional order chaotic dynamic system x(nq) = f(x,t) + g(x,t)u + d(t), y = x1 (12) where x = [x1,x2, ...,xn]T = [x,x(q),x(2q), ...,x((n−1)q)]T is the state vector, f(x,t) and g(x,t) are unknown but bounded nonlinear functions which express system dynamics, d(t) is the exter- nal bounded disturbance, |d(t)| ≤ D, and u(t) is the control input. The control objective is to force the system output y to follow a bounded reference signal yd which is the output trajectory of a drive system, under the constraint that all signals involved must be bounded. To begin with, the reference signal vector y d and the tracking error vector e will be defined as y d = [ yd,y (q) d , ...,y ((n−1)q) d ]T ∈ Rn, e = y d − x = [ e,e(q),e(2q), ...,e((n−1)q) ]T ∈ Rn,e(iq) = y(iq)d − x (iq) = y (iq) d − y (iq) In general, in the space of the error state a sliding surface is defined by s(x,t) = −(ke) = − ( k1e + k2e (q) + ... + kn−1e (n−2)q + e(n−1)q ) (13) where k = [k1,k2, ...,kn−1,1] in which the ki’s are all real and are chosen such that h(r) =∑n i=1 kir (i−1)q,kn = 1 is a Hurwitz polynomial where r is a Laplace operator. The tracking problem will be considered as the state error vector e remaining on the sliding surface s(x,t) = 0 for all t ≥ 0. The sliding mode control process can be classified into two phases, the approaching phase with s(x,t) ̸= 0 and the sliding phase with s(x,t)= 0 for initial error e(0) = 0. In order to guarantee that the trajectory of the state error vector e will translate from the approaching phase to the sliding phase, the sufficient condition s(x,t)ṡ(x,t) ≤ −η > 0 (14) must be satisfied. Two type of control law must be derived separately for those two phases described above. In the sliding phase, it implies s(x,t) = 0 and s(q)(x,t) = 0. In order to force the system dynamics to stay on the sliding surface, the equivalent control u can be derived as follows: If f(x,t) and g(x,t) are known and free of external disturbance, i.e., d(t)=0, taking the derivative of the sliding surface with respective to time, we get s(q) = − ( n−1∑ i=1 cie (iq) + e(nq) ) = − ( n−1∑ i=1 kie (iq) + y (nq)−y(nq) d ) = − ( n−1∑ i=1 kie (iq)−f(x,t)−g(x,t)ueq ) − y(n)d = − n−1∑ i=1 kie (i) + f(x) + b(x)u(t) − x(n)d = 0 (15) Therefore, the equivalent control can be obtained as 422 T.C. Lin, C.H. Kuo, V.E. Balas u = 1 g(x,t) ( n−1∑ i=1 kie (iq) − f(x,t) + y(nq)d ) (16) On the contrary, in the approaching phase, s(x,t) ̸= 0, an approaching-type control uap must be added in order satisfy the sufficient condition (4) and the complete sliding mode control will be expressed as u = u − uap, uap = ψhsgn(s) (17) where ψh ≥ η > 0. To obtain the sliding mode control (17), the system functions f(x,t),g(x,t) and switching parameter ψh must be known in advance. However, f(x,t) and g(x,t) are unknown and external disturbance, d(t) ̸= 0, the ideal control effort (16) cannot be implemented. We replace f(x,t), g(x,t) and uap by the fuzzy logic system f(x|θf), g(y|θg) and h(s|θh) in specified form as [9], [17]- [19], i.e., f(x|θf) = ξ T (x)θf,g(x|θg) = ξ T (x)θg, h(s|θh) = ∅ T (s)θh (18) let |h(s|θh)| = D + ψh + ωmax when s(x,t) is outside the boundary layer. Here the fuzzy basis functions ξ(x) and ∅(s) depend on the fuzzy membership functions and is supposed to be fixed, while θf,θg and θh are adjusted by adaptive laws based on Lyapunov stability criterion. Therefore, depending on plant knowledge and control knowledge, a hybrid adaptive fuzzy con- troller can be constructed by incorporating both fuzzy description and fuzzy control rules using a weighting factor α to combine the indirect adaptive fuzzy controller and the direct adaptive fuzzy controller. Based on the trade-off between plant knowledge and control knowledge, the weighting factor α ∈ [1,1] can be adjusted. Therefore, the total control effort can be expressed as uc = αui + (1 − α)ud (19) where the direct adaptive fuzzy controller ud and the indirect adaptive fuzzy controller ui are given as follows: ud(x) = uD(x|θD) − h(s|θh) g(x,t) and ui(x) = 1 g(x|θg) [ n−1∑ i=1 kie (iq) + y (nq) d − f(x, |θf) − h(s|θh) ] (20) where uD(x|θ) is obtained by fuzzy logic system specified as uD(x|θD) = ξ T (x)θD (21) The optimal parameter estimations θ∗f,θ ∗ g, θ ∗ h and θ ∗ D are defined as θ∗f = arg minθf ∈Ωf [ sup x∈Ωx |f(x|θf) − f(x,t)| ] , θ∗g = arg minθg∈Ωg [ sup x∈Ωx |g(x|θg) − g(x,t)| ] θ∗D = arg minθD∈ΩD [ sup x∈Ωx |uD(x|θg) − u| ] , θ∗h = arg minθh∈Ωh [ sup x∈Ωx |h(s|θh) − uap| ] Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control 423 where Ωf,Ωg,ΩD and Ωx are constraint sets of suitable bounds on θf,θg,θ ∗ h,θD and x respec- tively and they are defined as Ωf = {θf|θf| ≤ Mf}, Ωg = {θg|θg| ≤ Mg}, ΩD = {θD|θD| ≤ MD}, Ωh = {θh|θh| ≤ Mh} and Ωx = {x||x| ≤ Mx}, where Mf,Mg,MD,Mh and are positive constants. By using (20), (21), sliding surface equation (15) can be rewritten as s(q) = ω + α [ f(x|θ∗f) − f(x|θf) ] + α [ g(x|θ∗g) − g(x|θg) ] ui − (1 − α)h(s|θ∗h) −αh(s|θh) − (1 − α)g(x) [uD(x|θ ∗ D) − uD(x|θD)] + αh(s|θ ∗ h) − αh(s|θ ∗ h) (22) +(1 − α)h(s|θ∗h) − (1 − α)h(s|θ ∗ h) + d(t) where the minimum approximation errors is defined as ω = α [ f(x) − f(x|θ∗f) ] + α [ g(x) − g(x|θ∗g) ] ui + (1 − α) [uD(x|θ∗) − uD] (23) If θ̃f = θf − θ ∗ f , θ̃g = θg − θ ∗ g and, θ̃D = θD − θ ∗ D, we have s(q) = −(1 − α)h(s|θ∗h) + ω − αθ̃ T h ∅(s) − αθ̃ T f ξ(x) − αθ̃ T g ξ(x)ui +(1 − α)g(x)θ̃ T Dξ(x) − αh(s|θ ∗ h) − (1 − α)θ̃ T h ∅ + d(t) (24) Following the proceeding consideration, the following theorem can be obtained. Theorem: Consider the fractional order SISO nonlinear chaotic system (12) with control input (19), if the fuzzy-based adaptive laws are chosen as θ (q) f = r1sξ(x), θ (q) g = r2sξ(x)ui, θ (q) D = r3s∅(s) and θ (q) h = −r4sg(x)ξ(x) (25) where ri > 0, i = 1 ∼ 4. Then, the overall adaptive scheme guarantees the global stability of the resulting closed-loop system in the sense that all signals involved are uniformly bounded and the tracking error will converge to zero asymptotically. Proof: In order to analyze the closed-loop stability, the Lyapunov function candidate is chosen as V = 1 2 s2 + α 2r1 θ̃ T f θ̃f + α 2r2 θ̃ T g θ̃g + α 2r4 θ̃ T h θ̃h + (1 − α) 2r3 θ̃ T Dθ̃D (1 − α) 2r4 θ̃ T h θ̃h (26) Taking the derivative of the (26) with respect to time, we get V (q) = ss(q) + α r1 θ̃ T f θ̃ (q) f + α r2 θ̃ T g θ̃ (q) g + α r4 θ̃ T h θ̃ (q) h + (1 − α) r3 θ̃ T θ̃ q + (1 − α) r4 θ̃ T h θ̃ (q) h = −(1 − α)sh(s|θ∗h) + sω − αsθ̃ T h ∅(s) − αsθ̃ T h ξ(x) − αsθ̃ T g ξ(x)ui + (1 − α)sg(x)θ̃ T ξ(x) −ash(s|θ∗h) − (1 − α)sθ̃ T h ∅ + sd(t) + α r1 θ̃ T f θ̃ (q) f + α r2 θ̃ T g θ̃ (q) g + α r4 θ̃ T h θ̃ (q) h + (1 − α) r3 θ̃ T θ̃ q + (1 − α) r4 θ̃ T h θ̃ (q) h ≤ α r1 θ̃ T f ( θ̃ (q) f − r1sξ(x) ) + α r2 θ̃ T g ( θ̃ (q) g − r2sξ(x)ui ) + 1 r4 θ̃ T h ( θ̃ (q) h − r4s∅(s) ) − αs(D + η)sgn(s) + 1 − α r3 θ̃ T ( θ̃ (q) + r3sg(x)ξ(x) ) − (1 − α)s(D + ψh)sgn(s) + sd(t) + sω (27) From the robust compensator ua and the fuzzy-based adaptive laws are given (25), after simple manipulation, we have V (q) ≤ sω − sψhsgn(s) = sω − |s|ψh (28) Using the corollary of Barbalat’s Lemma [16]-[19], we have limt→∞ |s(x,t)| = 0. Therefore, limt→∞ |e(t)| = 0. The proof is completed. 424 T.C. Lin, C.H. Kuo, V.E. Balas 4 Simulation example In this section, we will apply our adaptive hybrid fuzzy sliding mode controller to force the fractional order chaotic gyro response system to track the trajectory of the fractional order chaotic gyro drive system. Example: The fractional order chaotic gyro drive and response systems are given as follows: Drive System: { y (q) 1 = y2 y (q) 2 = −100 ( y1 4 ) + y31 12 − 0.5y2 − 0.05y32 + sin(y1) + 35.5 sin(2t)y1 − x31 6 + d(t) Response System: { x (q) 1 = x2 x (q) 2 = −100 ( x1 4 ) + x31 12 − 0.7x2 − 0.08x32 + sin(x1) + 33 sin(2t)x1 − x31 6 + ∆f(x1,x2) + d(t) + uc(t) where structured uncertainty ∆f(x1,x2) = −0.1sin(x1) and external disturbance d(t) = 0.2cos(πt). The main objective is to control the trajectories of the response system to track the reference trajectories obtained from the drive system. The initial conditions of drive and response systems are chosen as [y1(0),y2(0)]T = [1,−1]T and [x1(0),x2(0)]T = [1.6,0.8]T , respectively. For q=0.95, α = 0.7 and all design constants are specified as k1 = k2 = 1,r1 = 150,r2 = 20,r3 = 1,r4 = 1 and step size h = 0.01. The phase portrait of the drive and response systems for free of control input is given in Figure 1. It is obvious that the tracking performance is bad without control effort supplied to response system. Figure 1: Phase portrait of chaotic drive and response systems The membership functions for xi i=1,2 are selected as follows: µF i1 (xi) = exp [ −0.5 ( xi−4 2 )2] , µF i2(xi) = exp [ −0.5 ( xi−2.7 2 )2] , µF i3(xi) = exp [ −0.5 ( xi−1.2 2 )2] , µF i4 (xi) = exp [ −0.5 ( xi 2 )2], µF i5(xi) = exp[−0.5(xi+1.22 )2], µF i6(xi) = exp[−0.5(xi+2.72 )2], µF i7 (xi) = exp [ −0.5 ( xi+4 2 )2] , From the adaptive laws (25)-(28), the control effort of the response system can be obtained as uc = αui + (1 − α)ud (29) Figure 2 shows the trajectories of the states xi,yi and x2,y2, respectively. Control effort trajectory is given in Figure 3 and phase portrait, tracking performance, of the drive and response Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control 425 Figure 2: The trajectories of the states xi,yi and x2,y2 Figure 3: Trajectory of the control effort Figure 4: Phase portrait, tracking perfor- mance, of the drive and response systems systems is shown in Figure 4. Trajectory of the sliding surface is given is Figure 5. The maximum value of V (q)(t) is -1.711e-4 which is always negative defined and consequently is stable. In order to show the robustness of the proposed adaptive hybrid fuzzy sliding mode control, the control effort is activated at 5 second. The phase portrait, tracking performance, of the drive and response systems is given in Figure 6. Figure 7 shows the trajectories of the states xi,yi and x2,y2 respectively. We can see that a fast tracking of drive and response is achieved as the control effort is activated. Control effort trajectory is given in Figure 8. Trajectory of the sliding surface is given is Figure 9. The maximum value of V (q)(t) is -1.732e-4 which is always negative defined and consequently is stable. Figure 5: Trajectory of the sliding surface Figure 6: Phase portrait, tracking perfor- mance, of the drive and response systems 426 T.C. Lin, C.H. Kuo, V.E. Balas Figure 7: The trajectories of the states xi,yi and x2,y2 Figure 8: Trajectory of the control effort Figure 9: Trajectory of the sliding surface 5 Conclusions A novel adaptive hybrid fuzzy sliding mode controller is proposed to achieve tracking perfor- mance of fractional order chaotic systems in this paper. It is a flexible design methodology by the trade-off between plant knowledge and control knowledge using a weighting factor ? adopted to sum together the control effort from indirect adaptive fuzzy controller and direct adaptive fuzzy controller. Based on the Lyapunov synthesis approach, free parameters of the adaptive fuzzy controller can be tuned on line by output feedback control law and adaptive laws. The simulation example, the output trajectory of the fractional order chaotic response system to tracking the trajectory of the fractional order chaotic drive system, is given to demonstrate the effectiveness of the proposed methodology. Bibliography [1] M. S. Tavazoei, M. Haeri, Synchronization of chaotic fractional-order systems via active sliding mode controller, Physica A, Vol. 387, pp. 57-70, 2008. [2] A.A. Kilbas, H.M. Srivastava and J.J. Trujillo, Theory and applications of fractional differ- ential equations, North-Holland Math. Studies 204, Elsevier, Amsterdam, 2006. [3] I. Petras, A note on the fractional-order Chua’s system, Chaos, Solitons & Fractals, 38 (1), pp. 140-14,7 2008. [4] X. Gao and J. Yu, Chaos in the fractional order periodically forced complex Duffing’s oscil- lators, Chaos, Solitons & Fractals 26, pp. 1125-1133, 2005. [5] J.G. Lu and G. Chen, A note on the fractional-order Chen system, Chaos, Solitons & Fractals 27, pp. 685-688, 2006. Uncertain Fractional Order Chaotic Systems Tracking Design via Adaptive Hybrid Fuzzy Sliding Mode Control 427 [6] P. Arena and R. Caponetto, Bifurcation and chaos in non-integer order cellular neural net- works, Int J Bifurcat Chaos 8 (7), pp. 1527-1539, 1998. [7] Petras I., A note on the fractional-order cellular neural networks. In: Proceedings of the IEEE world congress on computational intelligence, international joint conference on neural networks, Vancouver, Canada; pp.16-21, 2006. [8] S. H. Hosseninnia, R. Ghaderi, A. Ranjbar N., M. Mahmoudian, S. Momani, Sliding mode synchronization of an uncertain fractional order chaotic system, Journal Computers & Math- ematics with Applications, Vol. 59, pp. 1637-1643, 2010. [9] L.A. Zadeh, Fuzzy logic, neural networks and soft computing. Commun. ACM 37 3, pp. 77-84, 1994. [10] X. Z. Zhang; Y.N. Wang; X. F. Yuan, H? Robust T-S Fuzzy Design for Uncertain Non- linear Systems with State Delays Based on Sliding Mode Control, International Journal of Computers Communications & Control, 5(4):592-602, 2010. [11] Balas, M.M., Balas, V.E., World Knowledge for Control Applications by Fuzzy-Interpolative Systems, International Journal of Computers Communications and Control, Vol. 3, Supple- ment: Suppl. S, pp. 28-32, 2008. [12] L. X. Wang, and J. M. Mendel, Fuzzy basis function, universal approximation, and orthog- onal least square learning, IEEE Trans. Neural Networks, vol. 3, no. 5, pp. 807-814, 1992. [13] L. X. Wang, Adaptive Fuzzy Systems and Control: Design and Stability Analysis. Engle- wood Cliffs, NJ: Prentice-Hall, 1994. [14] K. Diethelm, An algorithm for the numerical solution of differential equations of fractional order, Elec. Trans. Numer. Anal. 5, pp. 1-6, 1997. [15] J. L. Castro, "Fuzzy logical controllers are universal approximators," IEEE Trans. Syst., Man, Cybern., 25, pp. 629-635, 1995. [16] S. S. Ge, T. H. Lee, C. J. Harris, Adaptive Neural Network Control of Robotic Manipulators, World Scientific Publishing Co., Singapore, 1998. [17] C. H. Wang, T. C. Lin, T. T.. Lee and H. L. Liu, "Adaptive Hybrid Intelligent Control for Unknown Nonlinear Dynamical Systems", IEEE Transaction on Systems, Man, and Cyber- netics Part B, 32 (5), October, pp. 583-597, 2002. [18] C. H. Wang, H. L. Liu and T. C. Lin, "Direct Adaptive Fuzzy-Neural Control with Observer and Supervisory Control for Unknown Nonlinear Systems", IEEE Transaction on Fuzzy Sys- tems, 10(1) , pp. 39-49, 2002. [19] T. C. Lin, C. H. Wang and H. L. Liu, "Observer-based Indirect Adaptive Fuzzy-Neural Tracking Control for Nonlinear SISO Systems Using VSS and ", Fuzzy Sets and Systems, 143, pp.211-232, 2004. [20] L. A. Zadeh, Knowledge Representation in Fuzzy Logic, IEEE Trans. Knowledge and Data Engineering, 1 (1), pp. 89-100, 1989.