Stable and quadratic-optimal parallel-distributed-compensation controller design for time-varying Takagi–Sugeno fuzzy model System: A complementary computational approach

A complementary computational approach is proposed for the time-varying Takagi–Sugeno fuzzy model system (TVTSFMS). The proposed approach integrates orthogonal-functional approach (OFA), hybrid Taguchi genetic algorithm (HTGA), and a stabilizability condition (SC) for use in designing stable and quadratic-optimal parallel-distributed-compensation (SQOPDC) controllers for optimal control problems. First, the SC was set according to linear matrix inequalities (LMIs). Next, OFA was used to derive an algorithm that only required algebraic computation to solve the TVTSFMS. Finally, The HTGA could be used to search the SQOPDC controller for the TVTSFMS. The SQOPDC controller obtained by the proposed complementary computational approach was evaluated in a case study of a vibratory pendulum design; the successful design verified the usability of the proposed hybrid intelligent computing method.


Introduction
The plant in an actual control system is mostly nonlinear. Therefore, scholars have developed many different approaches to designing controllers in nonlinear control systems in order to solve controller design problems in actual control systems. 1 In the Takagi-Sugeno fuzzy model, each fuzzy rule for the local linear control system is represented by a linear dynamic equation. [2][3][4][5][6][7] The overall nonlinear control system is then formalized by combining the fuzzy rules. Therefore, linear control theory is generally used in controller design. For nonlinear time-varying systems, it is difficult to design the controllers. Using Takagi-Sugeno fuzzy model (TSFM), the overall nonlinear time-varying system can be formalized by combining the fuzzy rules such that the linear control theory can be used in the controller design facilitating the controller easy implementation in practical situation. However, an important research issue is the design of stable and quadraticoptimal parallel-distributed-compensation (SQOPDC) controllers for a time-varying Takagi-Sugeno fuzzy model system (TVTSFMS) with a minimized performance index. Though the PI/PD/PID controller can be easily implemented in practical situation, it is very difficult to design the PI/PD/PID controller for the time-varying Takagi-Sugeno fuzzy model system (TVTSFMS) with a minimized performance index. In general, linear matrix inequalities (LMIs) are applied to design the SQOPDC controller of time-invariant Takagi-Sugeno fuzzy model systems. However, the LMIs are not directly applicable for solving the problem of designing an SQOPDC controller for a TVTSFMS. [8][9][10][11][12][13][14] Therefore, the objective of this study was to develop a complementary computational approach in which orthogonal-functional approach (OFA) is used to convert the SQOPDC controller design problem of the TVTSFMS into an algebraic computation problem. 15 To simplify the controller design problem, the proposed approach also integrates hybrid Taguchi genetic algorithm (HTGA) [16][17][18][19] and LMIs in the design process to ensure that the TVTSFMS can be stabilized in a closed loop. Here it should be noticed that, for the TVTSFMS, the OFA or the LMI technique cannot be applied alone to find the SQOPDC controllers, and it is also very difficult to apply the genetic algorithms alone to find the SQOPDC controllers. So, in this paper, we complementarily fuse the OFA, the HTGA and the LMI technique to solve the design problem to be studied. The proposed integrative method fusing the OFA, the HTGA and the LMI technique belongs to the hard-computingassisted soft-computing category, where the OFA and the LMI technique belong to the hard computing constituents and the HTGA is one of the soft computing constituents. That is, the main contribution of this paper is to integrate the OFA, the HTGA, and an LMI-based stabilizability condition (SC) for use in designing SQOPDC controllers for optimal control problems. A system block for the closed-loop system is shown in Figure 1.

Problem statement
A TVTSFMS can be expressed as follows: R i : IF z 1 ðtÞ is M i1 and … and z g ðtÞ is M ig with the initial state vector xð0Þ , whereR i ði ¼ 1, 2,…, N Þ denotes the i-th implication, N is the number of fuzzy rules, xðtÞ ¼ ½x 1 ðtÞ, x 2 ðtÞ,…, x n ðtÞ T denotes the n-dimensional state vector, uðtÞ ¼ ½u 1 ðtÞ, u 2 ðtÞ,…, u p ðtÞ T denotes the pdimensional input vector, z i ðtÞ ði ¼ 1, 2,…, gÞ are the premise variables, A i ðtÞ and B i ðtÞ ði ¼ 1, 2,…, N Þ are, respectively, the n × n and n × p consequent time-varying matrices, and M ij ði ¼ 1, 2,…, N and j ¼ 1, 2,…, gÞ are the fuzzy sets. The TVTSFMS in equation (1) has the following two characteristics: (i) time-varying elements a ijk ðtÞ in A i ðtÞ belonging to ½a ijk , a ijk , where a ijk and a ijk are given constants and (ii) time-varying elements b ijk ðtÞ in B i ðtÞ belonging to ½b ijk , b ijk , where b ijk b ijk are given constants.
Equation (1) can be rewritten as (2) in which zðtÞ ¼ ½z 1 ðtÞ, z 2 ðtÞ,…, x g ðtÞ T denotes the g-di-  20 and where F i is the local feedback gain matrices, which are expressed as In equation (3), the stabilizability problem is whether F i of the SQOPDC controller meets the pre-specified SC. Therefore, the following sections propose an SC based on LMI to formalize the instances of equation (3) that can be stabilized.
Theorem: To stabilize SQOPDC controller F i of the TVTSFMS in equation (3), there exists a symmetric positive definite matrix P that enables the controller to satisfy the following LMI where are constant matrices in which E j iαβ ¼ E iαβ denotes an n × n constant matrix with 1 in the αβ-th entry and 0 elsewhere and and V iαβ denote an n × p constant matrix with 1 in the αβ-th entry and 0 elsewhere.
Proof: See Appendix. Therefore, the question considered in this study is how to specify a value for F i of the SQOPDC controller such that (i) equation (3) meets the requirements of the LMIs in equation (4) and (ii) Control performance under the initial state conditions can be optimized by minimizing the following performance index where t f is the shortest time interval. Therefore, the SQOPDC controller design process is specifying F i and then executing the following two steps to optimize stability and control: Step 1: Verify that the SC of the LMI in equation (4) is the constraint condition.
That is, the design problem of the SQOPDC controller for the TVTSFMS is a constrained-optimization problem. In the next section, we will complementarily fuse the OFA, the HTGA and the presented LMI-based SC to solve the SQOPDC controller design problem of the TVTSFMS in equation (1), where the performance index J in equation (6) subject to the constraint of SC in equation (4) is considered to be directly minimized, where if there is no solution for equation (4), a penalty is given in the HTGA.

SQOPDC controller design
and Next, an orthogonal function (OF) can be used to obtain x ðkÞ s P s ðtÞ ¼x ðkÞ PðtÞ , and where PðtÞ ¼ ½P 0 ðtÞ, P 1 ðtÞ,…, P mÀ1 ðtÞ T represents the basis vector of OF, P i ðtÞ represents OF, and x ðkÞ s , A ðkÞ is and B ðkÞ is represent the coefficient matrices.
By substituting uðtÞ ¼ À P N i¼1 h i ðzðtÞÞF i xðtÞ and xðtÞ 1 x ðkÞ PðtÞ into equation (6), performance index J can be rewritten as where h i ðz k Þ ¼ h i ðzðk t f ÞÞ, and where constant matrix 13 From the orthogonal condition and the properties of OF, the product of any two OFs can be acquired and can be expressed by the following formula 15 P a ðtÞP b ðtÞ ¼ X mÀ1 s¼0 ξ abs P s ðtÞ: The consequent output only requires inference in an extremely short time interval. Therefore, integrating equation (2) yields The coefficient matrix of equation (14) can then be expressed as Given a set of local feedback gain matrices fF 1 , F 2 ,…, F N g,x ðkÞ can then be acquired by the algebraic calculation.
The above algebraic calculation clearly indicated that specifying a fF 1 , F 2 ,…, F N g enables the calculation ofx ðkÞ , which in turn facilitates the use of equation (12) to calculate the performance index. Therefore, this control problem can be converted into the following a static-constrained optimization problem If (i) constraints jf ijk j ≤ C ijk , are met and if (ii) the LMI of equation (4) is established (where C ijk is a given positive real number based on actual engineering considerations), the OFA and LMI can be used to obtain a constrained algebraic equation that represents the stability of the TVTSFMS and the SQOPDC control design problem. This representation substantially simplifies the SQOPDC control problem. The HTGA can then be used to determine the optimal solution for equation (16). [16][17][18][19] Remark 1: From the results mentioned above, we can see that, by using the OFA, we can transform the SQOPDC control problem for the TVTSFMS into a static optimization problem represented by the algebraic equations. Then, by incorporating the LMI-based SC for this TVTSFMS, the static optimization problem becomes a static constrainedoptimization problem represented by the algebraic equations with constraint of LMI-based SC. This means that the main characteristic and contribution of this technique, bringing the OFA and the LMI-based SC together, is that it reduces the mixed H 2 /LMI PDC controllers design problem to that of solving a static constrained-optimization problem of algebraic form; thus greatly simplifying the mixed H 2 /LMI PDC controllers design problem and facilitating the work of applying genetic algorithms to solve the SQOPDC controllers design problem. The proposed integrative method, that complementarily fuses the OFA, the HTGA and the LMIbased SC, belongs to the hard-computing-assisted softcomputing category.

Remark 2:
For the problem of designing the SQOPDC controllers of the TVTSFMS under the criterion of minimizing a performance index. But, the LMI-based approaches [8][9][10][11][12][13][14] cannot be applied to find the SQOPDC controllers of the TVTSFMS under the criterion of directly minimizing a performance index. In addition, though the genetic algorithms can solve the complex static-optimization problems that are not easy to analyze mathematically, it is also very difficult to only use the genetic algorithms to solve the dynamic-optimization problem of designing the SQOPDC controllers of the TVTSFMS. Hence, summing up the above statements and reasons, we can see that it is worth while to present a complementarily integrative approach fusing the HTGA, the LMI-based SC to find the SQOPDC controllers of the TVTSFMS under the criterion of minimizing a performance index.

Illustrative examples
Consider a vibratory pendulum system (Figure 2), 21     In equation (17), the proposed complementary computational approach is used to design the SQOPDC controller, to obtain a symmetric positive definite matrix P in equation (4), and to minimize J in equation (16). The OF type considered in this case was the Legendre function. Table 1 lists the average, median, minimum, maximum, and standard deviation of the performance indices obtained in 100 computations performed using the proposed HTGA. For comparison, the table also displays the results for 100 computations performed using a traditional genetic algorithm (TGA). 22 Table 1 and Figure 3 show the results of the performance comparisons of the HTGA and TGA. Compared to the TGA, the HTGA had (i) superior average and median performance indices, (ii) smaller SDs in performance indices, and (iii) superior convergence. Accordingly, HTGA yielded and more stable solutions. Therefore, compared with TGA, the HTGA is more effective for routine use in designing an SQOPDC controller. Figure 4 shows the x 1 ðtÞ, x 2 ðtÞ and uðtÞ responses of the vibratory pendulum system obtained after implementation of the SQOPDC controller designed using the HTGA ðF 1 ¼ ½ 0:2661 1:0889 and

Conclusions
This study developed an algebraic algorithm for using OFA to solve a TVTSFMS controller design problem. An algebraic algorithm was integrated in HTGA to design an SQOPDC controller for a TVTSFMS constrained by LMIs to minimize the performance index. The OFA was also used to convert the performance index to algebraic form. Application of the proposed approach in the design of a vibratory pendulum system verified the effectiveness of the approach for designing the SQOPDC controller for the TVTSFMS.
For the specified feedback gain matrices F i , if there exists a symmetric positive definite matrix P U T ijl P þ PU ijl < 0, then _ V ðxðtÞÞ < 0, and "xðtÞ ≠ 0. (25) Therefore, we can conclude that, to design an SQOPDC controller for the TVTSFMS (equation (3)), a series of F i in a symmetric positive definitive matrix P must be specified that simultaneously satisfies the LMI of equation (4). This completes the verification.