Recurrent Interval Type-2 Fuzzy Wavelet Neural Network with Stable Learning Algorithm: Application to Model-Based Predictive Control

Fuzzy neural networks, with suitable learning strategy, have been demonstrated as an effective tool for online data modeling. However, it is a challenging task to construct a model to ensure its quality and stability for non-stationary dynamic systems with some uncertainties. To solve this problem, this paper presents a novel identification model based on recurrent interval type-2 fuzzy wavelet neural network (RIT2FWNN) with new learning algorithm. The model benefits from both advantages of recurrent and wavelet neural networks such as use of temporal data and fast convergence properties. The proposed antecedent and consequent parameters update rules are derived using sliding-mode-control-theory. To evaluate the proposed fuzzy model, it is utilized to design a nonlinear model-based predictive controller and is applied for the synchronization of fractional-order time-delay chaotic systems. Using Lyapunov stability analysis, it is shown that all update rules of the parameters are uniformly ultimately bounded. The adaptation laws obtained in this method are very simple and have closed forms. Some stability conditions are derived to prove learning dynamics and asymptotic stability of the network by using an appropriate Lyapunov function. The efficacy and performance of the proposed method is verified by simulation examples.


Introduction
In the study of nonlinear dynamical system identification, conventional modeling approaches may not be suitable to be used due to the lack of precise, formal knowledge about the system, strongly nonlinear behavior, a high degree of uncertainty, and time-varying characteristics. Most of problems require methods that can handle qualitative and quantitative information with varying precision and complexity. These considerations impose extra demands on the effectiveness of process modeling techniques [30]. Modelfree methods like fuzzy logic and neural network are universal estimators that can be applied to approximate the behavior of the nonlinear dynamical systems [20]. Fusion of these techniques into a system enhances the capability of them like fuzzy neural network [8,18].
Since the output of dynamic systems is a function of time-delayed input and time-delayed output, recurrent neural network (RNN) is suitable choice for identifying their behavior. These networks have recurrent loops in their topologies that make it possible to temporarily store information, capture dynamic response of the system and enhance the approximation accuracy of the network [15,18]. In addition, combination of RNN and type-2 fuzzy sets is applied to various practical applications [11,17,18] to enhance noise tolerance of the system; different methods for their parameter tuning are discussed. The recurrent fuzzy structure that is used in [17] formed an external loop and internal feedback by feeding the rule firing strength of each rule to other rules and itself. The structure benefits from variable-dimensional Kalman filter and gradient descent (GD) algorithm for tuning its parameters. The recurrent Takagi-Sugeno-Kang (TSK) fuzzy structure in [10] obtains from locally feeding the firing strength of each fuzzy rule back to itself. An iterative linear support vector regression algorithm is applied to tune all parameters. In addition interval type-2 fuzzy system applied to represent the complex nonlinear plant in practical industrial systems [25,26]. For example on fault detection problem, the author in [25] utilized it to design an event-triggered fault detection filter for nonlinear networked systems which the parameter uncertainty is captured by type-2 fuzzy membership functions.
Due to the number of possible nonlinear interactions in nonlinear dynamical system is theoretically infinite, a nonlinear function should be chosen which is rich enough to describe a nonlinear process with good accuracy. Moreover, certain types of functions can efficiently approximate only certain nonlinear relationships [30]. Some basis functions have important general properties to deal with nonlinearity and uncertainties like wavelet representation.
Wavelet function as the activation function can enhance the advantages of neural networks for faster learning ability and wavelet decomposition for identification purposes. In literature, wavelet neural network (WNN) and its integration with fuzzy logic to determine an optimal definition of premise and consequent part of fuzzy rules, are applied in identification and control of nonlinear dynamical systems [5,19,20]. Discrete wavelet transform and WNN which is trained by back propagation and GD algorithms has also been used to improve the pattern recognition effects of sEMG signals [6].
Stability and convergence are fundamental in numerical analysis for online modeling, identification and control tasks. One of the foremost stability analysis techniques is the direct implementation of Lyapunov's stability theory. Another way of designing a robust and stable system is to use the variable structure systems (VSS) theory, which constructs the parameter adaptation mechanism and a rigorous stability analysis [9,29].
Sliding mode control (SMC) approach as a class of VSS has high performance in dealing with uncertainties and imprecision. Since robustness and invariance properties to matched uncertainties are the most significant properties of an SMC, the use of it in artificial neural networks (ANN) or fuzzy neural networks (FNN) can ensure convergence and stability of the learning algorithm [32]. As the result has shown in [29], SMC improves the performance of systems based on soft-computing techniques which utilize the gradient-based training strategies.
In literature, several research works use the diffusion of SMC approach into ANNs. In [12] a sliding mode incremental learning algorithm is used for tuning the parameters of interval type-2 FNN (IT2FNN) where an adaptive learning rate with an adaptation law is derived. An adaptive controller for speed control of induction motor which utilizes IT2FNN and SMC-based learning algorithm is proposed in [21]. In another study, a new learning algorithm for radial basis function neural networks is demonstrated that is based on applying fast terminal sliding mode [13]. In [33], a novel identification and control scheme using multitime scale recurrent high-order neural networks is presented. The scheme uses modified optimal bounded ellipsoid based weight's updating laws to identify a nonlinear systems.
Recently, the synchronization of the fractional-order chaotic systems have been studied for many applications in science and technology. It has been observed that many applications arising in various fields of science and engineering are described more accurately by fractional-order differential equations. Fractional calculus is a mathematical topic which deals with derivatives and integrations with non-integer order and applied to model nonlinear biological systems with complex behavior and long-term memory [2].
Chaos is a nonlinear and deterministic phenomenon which characterized by a subset of features comprising having an unusual sensitivity to initial states, not being periodic, having fractal structures, and being governed by one or more control parameters [7]. On the other hands, the chaotic dynamics of fractional-order differential systems has gained attraction in the investigation of dynamical systems, such as the fractional order of Chua's circuit, Duffing system [1], Lu system, Chen system [27].
In [4], a direct adaptive fuzzy controller is designed to obtain a generalized projective synchronization of two different incommensurate fractional-order chaotic systems in the presence of both uncertain dynamics and external disturbances. Lin and Lee developed an adaptive fuzzy-control scheme incorporating with SMC approach to synchronize two nonlinear fractional-order Duffing-Holmes chaotic systems [16]. They investigated the effect of delay on the chaotic behavior of the fractional-order system for the first time. Mohammadzadeh et.al formulated a robust nonlinear model-based predictive control for synchronization of the fractional-order chaotic systems [23]. By considering a fractional-order Lyapunov function, Mohadeszadeh and Delavari designed an adaptive finite time sliding mode control for chaos synchronization between two identical and nonidentical fractional-order hyper-chaotic systems [22]. The robust adaptive interval type-2 fuzzy-control strategy incorporating Lyapunov stability criterion and H 1 synchronization performance is studied in [14].
The models of real-world systems are not known completely and/or their accuracy are affected by nonlinear and time-varying behavior which can be originated from actual high degrees of uncertainties about the plant or from the plant dynamics, external disturbances and time-varying parameters. Identification performance of such applications is influenced by these factors. This study presents a new robust nonlinear and stable identification model RIT2FWNN to address these issues. In order to show the capabilities of the proposed model in such conditions, it is used to identify a second-order nonlinear time-varying plant and as a model in MPC to estimate the tracking error online. The controller is applied for the synchronization of uncertain fractional-order chaotic systems and the robustness of its synchronization in the presence of external disturbances, and approximation errors are investigated. The model performance in handling the uncertainties and imprecision can be achieved by utilization of recurrent IT2FNN and wavelet activation function. For the parameter adaptation of the RIT2FWNN, sliding mode theory-based supervised online training algorithm is elaborated. This learning algorithm is shown to be stable in the sense of Lyapunov stability theory which derives mathematical stability proofs to ensure finite time stability in RIT2FWNN. Moreover, fast convergence learning speed and the robustness of recurrent proposed structure is guaranteed. The main contributions of this manuscript are as follows: -The proposed structure of RIT2FWNN includes wavelet activation function with its parameters being tuned by a novel sliding-mode-control-theory learning algorithm. -Utilizing the Lyapunov stability theorem which is applied to obtain the parameter tuning algorithms. As it is expected, it is shown that the error value will be smaller in cases where the number of rules and input variables are more. -Since the performance of MPC highly depends on the accuracy of the model, a robust model-based predictive controller scheme is developed with the aid of the proposed model. -The robustness of the developed model-based predictive controller is examined in the case there exist uncertain dynamics and time delays in the dynamic system.
Different symbols defined in this paper are listed in Table 1. The organization of this paper is as follows. In Sect. 2, problem statement and preliminaries are given. The structure of RIT2FWNN and the novel SMC-theory-based learning rules are introduced in Sect. 3. The enhanced model-based predictive controller is presented in Sect. 4. To demonstrate the validation of the proposed method, synchronization of two fractional-order time-delay chaotic systems and identification of a nonlinear dynamic system are investigated in Sect. 5. Finally, the conclusions are given in Sect. 6.

Problem Statement and Preliminaries for Fractional-Order Systems
The commonly used definitions in literature for Fractional-Order operator are Grunwald-Letnikov, Riemann-Liouville, and Caputo definitions. The last one is introduced for engineering applications because its Laplace transform requires integer-order derivatives for the initial conditions. The Caputo's fractional derivative of a function x(t) with respect to time is defined as follows [28]: where m ¼ ½a þ 1, ½a is the integer part of a, D a t is called the a-order Caputo differential operator, and Cð:Þ is the well-known Euler's gamma function: This function can be considered as an extension of the factorial to real number arguments. The mathematical model of the fractional-order chaotic nonlinear systems to be discussed in this paper can be described as: where X ¼ ½x 1 ; x 2 ; . . .; x n T ¼ ½x; x ðaÞ ; x ð2aÞ ; . . .; x ðnÀ1Þa is the system's states vector, y 2 R outputs, u 2 R is the control input with the initial conditions, uð0Þ ¼ 0 and yð0Þ ¼ 0. If a 1 ¼ a 2 ¼ Á Á Á ¼ a n ¼ a the above system is called a commensurate order system. Then equivalent form of the above system is specified as: The unknown function f(X, t) is a bounded smooth nonlinear function which specifies system dynamics and xðtÞ is the external bounded disturbance. The main objective is to force the system output y to follow a given bounded reference signal y d while assuring under certain constraint, all involved signals are bounded. To quantify this objective, the reference signal vector y d and the tracking error vector e are defined as, By substituting (6) into (4), the control system in the state space domain is obtained as follows: where The equivalent control effort can be considered as follows: where u s is predictive control signal in previous sample times which is designed based on model predictive control. K ¼ ðk 1 ; k 2 ; . . .; k n Þ T 2 R is chosen such that the stability condition jargðeigðAÞÞj [ qp=2, 0\q\1 is satisfied, eig(A) is the eigenvalues of the system state matrix given in (7). To design u s , dynamic of the tracking error is modeled by the proposed RIT2FWNN with its premise and consequent part as well adaptive learning rate are updated using the proposed online learning algorithm based on SMC. Fig.1 shows a conceptual diagram of the control scheme. The structure of RIT2FWNN implements a recurrent wavelet fuzzy model is organized into layers whose consequent part is nonlinear function of the input variables. Each rule in RIT2FWNN has the following form: R r : If x 1 isÃ 1j Á Á Á x i isÃ ih and Á Á Á and x I isÃ Il then: where I is the number of input variables, h is the number of membership function and f r is the output of rth rule ðr ¼ f1; 2; . . .; NgÞ.Ã ik is the kth type-2 fuzzy membership function (MF) related to ith input variable. q is the weight coefficient between the input and the hidden layer. Consequent part of above rule involved a wavelet function of input variables. Wavelets are defined by a family of functions a and b ða [ 0; b 2 RÞ as follows: where W r ðxÞ demonstrates the family of wavelets obtained from the single wðxÞ function, called as a mother wavelet and localized in both time and frequency space, by dilation and translation, a r ¼ fa r1 ; a r2 ; . . .; a rI g and b r ¼ fb r1 ; b r2 ; . . .; b rI g, respectively and x r ¼ fx 1 ; x 2 ; . . .; x I g are input variables. In this manuscript, among several families of wavelets, Mexican Hat is considered as mother wavelet function. It is derived from a function that is proportional to the second derivative function of the Gaussian probability density function.
where, z ri ¼ x i Àb ri a ri . RIT2FWNN consists of seven layers. In the first layer, nodes representing input linguistic variables are fed into the network. In the second layer, each node corresponds to one linguistic term. For each entering input variable, type-2 MFs are used which have uncertain standard deviation and fixed center. The membership degree l ik ðx i Þ and l ik ðx i Þ are calculated according to (13): where c ik is center of type-2 MF, r ik and r ik are the upper and lower standard deviation of the kth type-2 MF of ith input. Moreover, n ik is the recurrent parameter defined as (14) which store the past information of the network. It is to be noted that the feedback weights of nodes in layer 2 are interval values.
where h ik and h ik are considered as feedback weights of the nodes. Nodes in the third layer correspond to one fuzzy rule and perform a fuzzy meet operation on inputs from layer 2 to obtain upper and lower firing strength as follows: Layer 4 determines the normalized values of the lower and the upper firing strength corresponding to each node in layer 3: Layer 5 is consequent layer. In this layer, Nodes compute the product of normalized firing strengthw r ,w r and wavelet function of input variables. Layer 6 involves two summation blocks, one of them is for upper and the other is for lower outputs of the previous layer. Layer 7 calculates the output of the network using (17): where q is the design parameter which enables to adjust the contribution of the lower or the upper values of rules depending on identification requirements for the system [3]. In literature, q has either been considered to be a constant or a time-varying parameter [12]. In this paper, the adaptation laws for the parameters and the proof of the stability of the learning process are given using a timevarying q.
Remark 1 Main improvement made in the proposed approach with respect to existing approaches in literature [25,26] are use of wavelet functions which makes it possible to benefit from a local nonlinear general function approximator. This in turn will improve function approximation property of the system and provides means for a tight control system. On the other hand, the existence of recurrent term in type-2 fuzzy membership functions allows the use of a short term memory which captures the dynamics of the system and improves the performance of controller.

Sliding Mode Online Learning Algorithm
Based on the principles of the SMC theory, the zero value of the learning error coordinate can be defined as a timevarying sliding surface, i.e.
which is the condition that guarantees when the system is on the sliding surface S the RIT2FWNN output y N ðtÞ will track the desired output signal y(t) for all time t [ t h where t h is the hitting time of eðtÞ ¼ 0. The online learning algorithm for the adaptation of the parameters of RIT2FWNN should be derived in such a way that the sliding mode condition of the above definition will be enforced.
Theorem 1 If the learning algorithm for the parameters of upper and lower Gaussian MFs and weights of feedback loops are chosen as follows: and the adaptation of the parameters of wavelet functions and coefficient in the consequents part of fuzzy rules are chosen as follows: qw r þ ð1 À qÞw r ðqw r þ ð1 À qÞw r Þ T ðqw r þ ð1 À qÞw r Þ asgnðeÞ qw r þ ð1 À qÞw r ðqw r þ ð1 À qÞw r Þ T ðqw r þ ð1 À qÞw r Þ asgnðeÞ ð25Þ _ q ri ¼ À 1 W ri ðzÞ qw r þ ð1 À qÞw r ðqw r þ ð1 À qÞw r Þ T ðqw r þ ð1 À qÞw r Þ asgnðeÞ ð26Þ and the parameter q is updated as follows: FðW ÀWÞ T asgnðeÞ ð27Þ and the adaptive learning rate ðaÞ changes as follows: The parameter c is considered as the learning rate for the adaptive learning rate a which has a small positive real value. The parameter m should be very small not to interrupt the adaptation mechanism. Then, given an arbitrary initial condition e(0), the learning error e(t) will converge to zero within a finite time t h .
Proof The reader is referred to Appendix. According to the following equation which obtained from the proof, the error value will be smaller in cases where the number of rules and input variables are more.
h In Eq. (28) the parameter c is the learning rate for the adaptive learning rate a. Since the first term of this equation is positive, the second term is considered to avoid bursting in parameter a. The value of the parameter m should be adjusted in such a way that the adaptation mechanism not to be disturbed.
Despite aforementioned advantages of simplicity and robustness of SMC, it suffers from a technical problem, called chattering. Sensitivity to noise is another phenomena caused by the presence of signum function in the structure of SMC. In this manuscript, in order to alleviate the problem associated with the chattering and sensitivity to noise, a continuous approximation method is used which smooths the discontinuity caused by the signum function as follows: where d is a small positive number [32].

Model Predictive Control Design
The target of the following Eq.
where u s is predictive control signal, eðt À k l Þ; l ¼ 1; . . .; r are the tracking errors, Du s ðtÞ ¼ u s ðtÞ À u s ðt À 1Þ, R and Q are positive definite weighting matrices. N c and N p are control and prediction horizon. In order to minimize the cost function (31), the output predictions over the horizon must be computed. RIT2FWNN is used to estimate the tracking error e(t). The consequent parameters of RIT2FWNN are tuned based on the updating rule proposed in previous section. By considering the Diophantine equation corresponding to the prediction forêðt þ k þ 1jtÞ: The last two terms of Eq. (32) depend on past values of the process output and input variables. That is, they correspond to the free response e r of the process considered if the control signals are kept constant and are computed as follows: e r ðt þ kjtÞ ¼ RIT2FWNN½eðt þ k À 1Þ; . . .; eðt þ k À rÞ; u s ðt À 1Þ ð33Þ The first term of Eq. (32) depends only on future values of the control signal and can be interpreted as the forced response which is obtained by (34). That is, the response obtained when the initial conditions are zero eðt À kÞ ¼ 0, Du s ðt À kÞ ¼ 0.
where # i , i ¼ 0; . . .; k À 1 are the step response coefficients that are obtained using unit step on RIT2FWNN model.
Step response can be estimated as follows [24]: gðt À 1Þ ¼ê step ðt þ kjtÞ Àê r ðt þ kjtÞ du s ðtÞ ð35Þ where du s ðtÞ is the step size andê step ðt þ kjtÞ is calculated by RIT2FWNN like as e r . Equation (32) can be rewritten as follows: such that f k ¼ G kp ðz À1 ÞDu s ðt À 1Þ þ F k ðz À1 ÞeðtÞ. Then the step ahead prediction of the system output on data up to time N p in matrix form described as follows: eðt þ 1jtÞ eðt þ 2jtÞ . . .  The predictions can be expressed in condensed form as follows: The optimization problem (31) can be written as follows: By making the gradient of J in (38) equal to zero: Then the optimum is:

Simulation Studies and Discussions
In this section, two illustrative examples are presented to demonstrate the applicability and feasibility of the proposed method and to confirm the theoretical results.
In these examples, the designed controller is applied to synchronize two identical and nonidentical uncertain fractional-order chaotic systems in the presence of external disturbances. Besides, to analyze how the proposed structure measures the mathematical model of a system from measurements of the system inputs and outputs, an example of nonlinear system identification is considered.
Example 1 As the first example, the proposed identification model, is applied to identify a second-order nonlinear time-varying plant [12] described as: where xð1Þ ¼ yðk À 1Þyðk À 2Þyðk À 3Þuðk À 1Þ, xð2Þ ¼ y ðk À 3Þ À bðkÞ, xð3Þ ¼ cðkÞuðkÞ and xð4Þ ¼ aðkÞ þ yðk À 2Þ 2 þ yðk À 3Þ 2 . The time-varying parameters a, b and c are defined as: aðkÞ ¼ 1:2 À 0:2 cosð2pk=TÞ bðkÞ ¼ 1 À 0:4 sinð2pk=TÞ where T ¼ 1000 is the time span of the test. This example has been implemented with three inputs which are fuzzified with three Gaussian type-2 MFs with a fixed center and uncertain standard deviation. The inputs are the delayed signal from the plant outputyðt À T 0 Þ; yðt À 2T 0 Þ and input signal u(t), with the period of discretization selected as T 0 ¼ 1 ms. The incoming network signals have been normalized to be in the range ½À1; 1. The input signal u(k) is as follows: uðkÞ ¼ sinðpk=25Þ; k\250: 1:0; 250 k\500: À1:0; 500 k\750: 0:3 sinðpk=25Þ þ 0:1 sinðpk=32Þ þ0:6 sinðpk=10Þ; 750 k\1000: Figure 2a compares the plant output with the output yielded by RIT2FWNN. Figure 2b shows the root mean square error (RMSE) values versus each epoch number which illustrates the proposed SMC-based learning algorithm is stable. Figure 2c, d show the adaptation of parameter a and q, respectively. Figure 2c shows the learning rate a will be stable at infinity. Table 2 compares the performance of RIT2FWNN with two different approaches and [12]. These results reveal the benefit of using recurrent structure and SMC-based learning algorithm of RIT2FWNN by achieving smaller RMSE values. It is to be noted that by increasing the nonlinearity of the RIT2FWNN by using wavelet functions in the consequents part of fuzzy rules and feedback loops on type-2 MFs, the accuracy of other learning methods such as gradient-based methods are lowered considerably.
Furthermore, the accuracies of extended Kalman filter and the proposed SMC-based learning method are similar. As previously discussed, the adaptation laws obtained in the proposed SMC-based learning method have explicit forms which do not have any derivative computations of the output of RIT2FWNN concerning trainable parameters and manipulation of some high-order matrices. This is the main reason why the computation time for extended Kalman filter and gradient-based methods is higher than SMCbased methods (see Table 2). Additionally, further results are provided to demonstrate the benefits of the proposed model with its type-1 counterpart. As can be seen, the proposed type-2 model achieves smaller test and train RMSE than that by the type-1.

Remark 2
As the simulation results show, the proposed model have better RMSE results than others. However it is more complicated in network structure. The best advantages of the proposed RIT2FWNN with respect to other such as [12], is the rate of convergence. Most controllers in real-world application are based on model and it is needed to be achieved faster. Hence, the degree of freedom in the proposed model is increased to obtain a faster convergence speed. Furthermore, several experiments have been done on the speed of convergence of different models. The average arising times are approximately 29 and 38 s for the proposed and [12] respectively, which confirm the favorable convergence behavior of the proposed model versus time.
Example 2 As the second example, we will verify the effectiveness of the proposed method to synchronize two different uncertain fractional-order Duffing-Holmes chaotic time-delay systems [1]. The master and slave systems dynamics are given as follows, respectively: D a x 1 ðtÞ ¼ 2:5x 2 ðtÞ D a x 2 ðtÞ ¼ À ð 1 2:5 x 1 ðtÞÞ 3 À 1 2:5 x 1 ðtÞ À 0:1x 2 ðtÞ þ 0:01x 1 ðt À 0:001Þ þ 0:01x 1 ðtÞ 2 ðt À 0:001Þ þ 0:01x 2 ðt À 0:001Þ þ 25 cosð1:29tÞ: D a y 1 ðtÞ ¼2:2y 2 ðtÞ D a y 2 ðtÞ ¼ À 1 2:0 y 1 ðtÞ 3 À 1 1:8 y 1 ðtÞ À 0:1y 2 ðtÞ þ0:01y 1 ðt À 0:001Þ þ 0:01y 1 ðtÞ 2 ðt À 0:001Þ þ0:01y 2 ðt À 0:001Þ þ 25 cosð1:29tÞ þ xðtÞ þ uðtÞ: where xðtÞ ¼ 0:7 sinðtÞ is the external bounded disturbance, and u(t) is control input of the slave system. The initial conditions are xð0Þ ¼ ½x 1 ð0Þ; x 2 ð0Þ T ¼ ½0; 0 and yð0Þ ¼ ½y 1 ð0Þ; y 2 ð0Þ T ¼ ½1; À1. The simulation sample time is 0.001 and the fractional derivative order is considered as a ¼ 0:98. The control effort of the slave can be obtained as follows: where K ¼ ½900; 30. The following cost function is considered to achieve the control signal u s : Each three inputs of RIT2FWNN are fuzzified with two Gaussian type-2 MFs with a fixed center and uncertain standard deviation. The premise and consequent parameters of RIT2FWNN are tuned based on the proposed online SMC-based learning algorithm. Figures 3 and 4 show the control trajectory and tracking error, respectively. The trajectories of the states x 1 , y 1 and x 2 , y 2 are depicted in Fig. 5 that demonstrate synchronization is perfectly achieved. The three-dimensional phase portrait, synchronization performance, of the chaotic master and slave systems is shown in Fig. 6. It can be seen that the synchronization performance of the proposed method in the presence of external disturbances and is successfully realized and the control trajectory is also admissible and bounded. The mean square error (MSE) values for the synchronization errors (e 1 , e 2 ) are given in Table 3. It can be seen that the performances of the proposed method in this paper are significantly better than other approaches; (1) sliding mode technique [16], (2) nonlinear model predictive controller (NMPC) [23]. Additionally, tracking error results in the presence of pulse disturbance are depicted in Fig. 7. As can be seen, the controller is robust against external disturbances.

Example 3
In this example, we will apply the proposed method to synchronize two nonidentical fractional-order time-delay chaotic systems. The master system is as follows [31]: D a y 1 ðtÞ ¼ y 3 À 3y 1 ðtÞy 2 ðt À 0:001Þ þ x m 1 ðtÞ D a y 2 ðtÞ ¼ 1 À 0:1y 2 ðtÞ À y 2 1 ðt À 0:001Þ þ x m 2 ðtÞ D a y 3 ðtÞ ¼ Ày 1 ðt À 0:001Þ À y 3 ðtÞ þ x m 3 ðtÞ The initial conditions are yð0Þ ¼ ½y 1 ð0Þ; y 2 ð0Þ; y 3 ð0Þ T ¼ ½0:1; 4; 0:5 and the simulation sample time is 0.001.The slave system is Liu fractional-order time-delay chaotic system which its dynamics are as follows [31]: The trajectories of the states x k , y k where k ¼ 1; 2; 3 and the control signals are shown in Figs. 8 and 9, respectively. As one can see from the figure, the tracking performance is  good even in the presence of disturbances and unknown functions of dynamics of the system. The designed controller can synchronize effectively two nonidentical fractional-order time-delay chaotic systems. From Fig. 10, it can be observed that the tracking errors e 1 ðtÞ, e 2 ðtÞ and e 3 ðtÞ converge both to small intervals around zero. Table 4 indicates that the obtained MSE results are comparable with the results in [31]. Figure 11 shows the tacking error results in the presence of pulse disturbances. It is shown that the developed controller preserves the robustness properties against external disturbances.

Conclusion and Future Work
In this paper, a novel robust model RIT2FWNN with fully SMC-based learning algorithm is proposed to design a model-based predictive controller to deal with synchronization between two different fractional-order time-delay chaotic systems with subject to uncertainty and external additive disturbances. This synchronization is effectively obtained which demonstrate the superior performance of the proposed methodology. The model is tested against identification of a nonlinear and time-varying dynamic system. Since the output of these system depends on different step delayed outputs and inputs, the task of identification of such systems is difficult. The simulation results indicate that the proposed approach is quite useful in modeling unknown function of dynamic systems, having not only favorable tracking performance but also produce strong robustness and faster convergence speed. The Lyapunov stability theorem is applied to derive the parameter update rule of the model to guarantee all parameters uniformly ultimately bounded. The proposed RIT2FWNN owns some advantages:  However, the proposed RIT2FWNN network has disadvantages such as fixed structure which makes it complex spatially for systems with higher numbers of inputs. As future work, there are many challenging issues remaining to investigate, including the choice of higher order sliding mode control in the process of training, employing self-structuring methods to achieve optimized construction of RIT2FWNN network.  The time derivative of (16) is calculated as follows: considering the following equations for _ w r and _ w r : thus the Eq. (51) changes as follows: Furthermore, the time derivative of upper and lower of type-2 MFs are calculated as (54): and the time derivative of the recurrent parameter is as follows: considering _ l ik ðx i Þ ¼ ÀA ik _ A ik where A ik and _ A ik are as follows: By using Maclaurin series expansion: D ik and D ik are limited as follows: thus the J r and J r in (53) are as follows: by substituting (19)-(23) in the above equation: In order to check the stability condition, the following Lyapunov function is considered: The stability condition is met when the time derivative of V is negative definite ( _ V\0). So, by differentiating (17) by substituting (53): by considering (60): since P N r¼1w r ¼ 1 and P N r¼1wr ¼ 1 Eq. (65) can be summarized and then by substituting _ q according to (27), q ri ÀW ri ðzÞ _ x i a ri " " 3 À z 2 ri 1 À z 2 ri z ri qw r þ ð1 À qÞw r where, a Ã !
which is regarded as an unknown parameter and is determined during the tuning of learning rate. Thus, the adaptive learning rate ðaÞ changes as follows: in which m has small real value. The time derivative of Lyapunov function (72) is changed as follow: consequently, In this case in order to have a negative definite time derivative of the Lyapunov function, it is required that the following condition holds.
Hence, it is concluded that error converges to a small neighborhood of zero, in which and stays there. Furthermore, the radius of this neighborhood can be made as small as desired using appropriate values for a Ã and m. This concludes the proof.