# Decentralized Polynomial Observer Design for Discrete-Time Large-Scale Polynomial T-S Fuzzy System

Authors: Seung-Hun Han; Van-Phong Vu; Manh-Son Tran
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In this paper, a new approach for synthesizing a decentralized observer is proposed to estimate the unmeasurable states of a discrete-time large-scale nonlinear system. The large-scale nonlinear system in this study is modeled in terms of a set of discrete-time polynomial Takagi-Sugeno (T-S) fuzzy subsystems and interconnection terms. This modeling method will assist to reduce significantly the number of fuzzy rules. The interconnection terms are considered as the unknown inputs; then, the unknown input method is employed to design observer for this system. It should be emphasized that the interconnection parts in this paper are arbitrary and their effects are eliminated completely. On the basis of the Lyapunov methodology and SOS (Sum of Square) technique, the conditions for observer design expressed under the framework of SOS are derived in the main theorems. Finally, an illustrative example is presented to show the effectiveness and merit of the proposed method.


---

## Body

## 1. Introduction

Nowadays, due to the development of technology and human’s demand, the complexity and large-scale of the system are significantly increasing. These complicated systems are called large-scale systems which consist of a set of interconnection subsystems. The large-scale system has appeared in many real applications such as communication systems, power systems, ecosystems, and transportation systems. There are plenty of studies focusing on stability analysis, controller design, and observer design for the large-scale nonlinear systems in recent years [1–4]. For example, the methods to design the decentralized controllers for the large-scale nonlinear systems were investigated in [1, 2], where the interconnection terms of these systems were unknown and estimated by using neural networks. In [4], an optimal controller design method based on the H∞ technique was proposed for a large-scale nonlinear system to stabilize the system and deal with the impacts of disturbance and saturated actuator.Recently, a method for modeling nonlinear systems called Takagi-Sugeno fuzzy system was introduced in [5] in which a nonlinear system was represented in terms of a set of IF-THEN rules and linear subsystems. This method came out in 1984 and, until now, the T-S fuzzy model increasingly plays a very important role in modeling the nonlinear systems. Numerous studies focusing on the T-S fuzzy system could be found in books [6, 7]. Moreover, the T-S fuzzy model is recently employed to represent the large-scale nonlinear system which is called “large-scale T-S fuzzy system”. In the past few years, there have been a lot of papers concentrating on solving the issues of the large-scale T-S fuzzy system such as stability analysis, controller design, and observer synthesis [8–15]. In [8], an approach for designing controller was investigated for large-scale nonlinear systems in which the interconnection parts were with the nonlinear form. Moreover, a method for stability analysis and controller design based on H∞ was proposed for the discrete-time large-scale T-S fuzzy system [9] where the piecewise Lyapunov function was employed instead of the common Lyapunov function. In [11], a method that relied on H∞ technique for designing controller was studied for large-scale T-S fuzzy systems with the existence of both disturbance and uncertainties. However, the drawback of the modeling method based on the T-S fuzzy model is that when the number of subsystems or number of linearization terms of each subsystem increases, the number of fuzzy rules will exponentially grow up. This results in the increase of the computational load.Recently, an extended form of the T-S fuzzy model called “polynomial T-S fuzzy model” was introduced in [16]. The major difference between the T-S fuzzy system and the polynomial T-S fuzzy system was that the elements of the system, input, and output matrices of the polynomial T-S fuzzy model were in polynomial form instead of the constant of the T-S fuzzy model. By putting the variables inside the system matrices, it will help to reduce the number of linearization terms and lead to significant decrease in the modeling error and number of fuzzy rules with respect to the T-S fuzzy model. With the aid of the Sum-of-Square (SOS) technique [17], the polynomial T-S fuzzy system becomes an increasingly popular tool for modeling nonlinear systems. In recent years, there have been a lot of researchers focusing on stability analysis, controller design, and observer synthesis for polynomial T-S fuzzy systems [18–20]. For example, in [18], observer and observer-based controller were synthesized for the polynomial T-S fuzzy system with three classifications. Moreover, a method to design a polynomial observer-based controller for the polynomial T-S fuzzy system was introduced in [19] in which the premise variables were unmeasurable and they need to be estimated as well. However, to the best of our knowledge, there is no previous study applying the polynomial T-S fuzzy framework for modeling discrete-time large-scale nonlinear system.In practice, it is very difficult to measure some physic parameters of many practical systems. Moreover, in some cases, using the sensor to obtain the information of the variables is costly or encounters technical problems such as the impact of noise. However, these variables are necessary for monitoring system and designing controller. Therefore, replacing sensors by the observer is a pressing issue and plays a substantial role in the control field. Regarding observer synthesis for large-scale systems, many papers have been published in the past few years [21–29] focusing on solving this issue. For example, the decentralized H∞ filters were synthesized for the nonlinear large-scale system and the nonlinear large-scale system with unknown interconnection terms in [21, 22], respectively. However, the interconnection term in [21] must be modeled in linear form and the interconnection part in [22] must satisfy the bounded constraints. Moreover, an approach to design a robust observer for a multimachine power system expressed in terms of nonlinear large-scale system was concerned in [23] in which the sliding mode technique was employed to eliminate the effect of the uncertainties and estimate the unknown states. In [25], an observer was synthesized to detect faults of the affine large-scale T-S fuzzy system. It should be noted that the interconnection terms were unknown, but the interconnection terms in [25] have to fulfill the norm bounded constraints. Additionally, the filter and observer were synthesized for discrete-time large-scale T-S fuzzy system in [26–28]. However, the interconnection terms in these papers must be expressed under the linear form of the T-S fuzzy model. Otherwise, the proposed methods in [26–29] failed to design the filter for these systems. In other words, there are two forms of the interconnection terms in the previous studies [21, 22, 25–29]: one is known and other is unknown. The advantage of the unknown interconnection terms is that we do not need to know exactly the form of interconnection terms and these terms do not need to linearize to the linear form or T-S fuzzy model in order to reduce the modeling error and complexity. However, the unknown interconnection terms must satisfy the specific constraints such as the bounded constrains in [22, 25]. Regarding knowing interconnection terms, their advantage is that they are unnecessary to fulfill any constraint. However, they have to be able to model under the linear form or T-S fuzzy model [21, 26–29]; it leads to increase in the modeling error, complexity, and computational load during designing observer/controller for the large-scale systems.From above discussions, it is easily seen that employing the T-S fuzzy model [8–15] for representing large-scale nonlinear system will cause to significantly increase the number of fuzzy rules and modeling errors and make the computational load heavier. In addition, the interconnection terms of the large-scale system [22, 25] must be bounded and in [21, 26–29] have to model under the T-S fuzzy model form. However, in practice, in some cases, the upper bound of the interconnection terms is not easy to determine and the interconnection terms of several systems are unable/difficult to be represented in the T-S fuzzy model; thus, the proposed methods in the previous papers are not applicable to these systems. Additionally, to the full extent of our knowledge, the polynomial T-S fuzzy model has not been applied to represent the discrete-time nonlinear system in any previous study. Due to these reasons, we are inspired to study further a new method for synthesizing a decentralized observer to estimate the unknown states of the discrete-time large-scale polynomial T-S fuzzy system with arbitrary interconnection terms. The main contributions of this paper are emphasized in the following aspects:(1)
This is the first time the polynomial T-S fuzzy system used to model the discrete-time large-scale nonlinear systems. This modeling method will help to reduce the number of fuzzy rules, modeling errors, and computational load with respect to large-scale T-S fuzzy systems in [8–15].(2)
In this paper, the interconnection terms are arbitrary and are treated as unknown inputs. The proposed approach allows us to apply for the large-scale systems in which the interconnection terms are not known exactly. In addition, the interconnection term in this paper is unnecessary to satisfy the bounded constraints which are mandatory in [22, 25]. Moreover, with the proposed method, the interconnection terms do not need to linearize to transform them under the T-S fuzzy model form [21, 26–29], which assists to avoid increasing of the fuzzy rules and complexity of the modeled system. The unknown input method [30–34] is employed to design observer for the large-scale polynomial T-S fuzzy system in which the effects of the interconnection terms are eliminated completely and unmeasurable states are estimated asymptotically.(3)
The fuzzy Lyapunov function is applied in this paper instead of the common Lyapunov function to make the conditions for observer design to be more relaxed.The rest of this paper is organized as follows. The considered system and the solved problem in this paper are presented in Section “System Model and Problem Description.” The main theorems for designing a decentralized polynomial observer are shown in Section “Decentralized Observer Design”. A numerical example to prove the effectiveness of the proposed method is mentioned in Section “Illustrative Example”. Finally, several conclusions are stated in Section “Conclusion”.Notations.In this paper, Λ>0(Λ<0) indicates the positive (negative) definite matrix. ΛT stands for the transpose of matrix Λ. Λ-1 represent the inverse of matrix Λ and I indicates the identity matrix. Λ+ stands for the Moore-Penrose pseudoinverse of Λ where Λ+=(ΛTΛ)-1ΛT. (∗) stands for the symmetric term of a matrix.


## 2. System Model and Problem Description

Consider the discrete-time large-scale nonlinear system withN subsystems as follows.Subsystem  l (1) x l k + 1 = h l x l k , u l k + g l x y l k = C l x l k l = 1,2 , … Nwhere xlk and ulk are the vectors of the state variables and inputs of thel-thsubsystem. hlxlk,ulk is the nonlinear function and glx is the interconnection term which expresses the interaction betweenl-th subsystem with others. x is a vector of states which consist of all states xlk.Assumption 1.
Assume that the interconnection part of system (1) satisfies the matching condition glx=Tlϑlx. The matrices Tl and Cl are full column and row rank, respectively. The rank of (ClTl) is equal to the rank of (Tl).Remark 2.
The rank conditions of matricesTl and Cl in Assumption 1 are necessary to guarantee the existence of the general solution of matrix equation in the consequent section. This is a common condition that can be found in many previous papers studying about unknown input method such as [30–34].On the basis of Assumption1 and applying the sector nonlinear or linearization method in [6, 7], system (1) is transformed into the discrete-time large-scale polynomial T-S fuzzy system as follows.Subsystem l-thRule iIfθ1lk is Mi1l, …, and θslk is Misl, then(2a)xlk+1=Ailζlkxlk+Bilζlkulk+Tlϑlx(2b)ylk=Clxlkin which θlk=θ1lk,θ2lk,…,θslk is a vector of the premise variables which is known in advance, Mijl (i=1,2,…,rl; j=1,2,…,s) is the fuzzy set, and rl and s are the number of fuzzy rules and premise variables, respectively.Overall large-scale system (2) is represented.S u b s y s t e m   l (3a) x l k + 1 = ∑ i = 1 r l β i l θ l k A i l ζ l k x l k + B i l ζ l k u l k + T l ϑ l x (3b) y l k = C l x l kwhere Ailζlk,Bilζlk are the time-varying system and input polynomial matrices with time-varying entry ζlk. Tl and Cl are the interconnection and output matrices. ylk is the vector of the outputs. ζlk is the time-varying variables which could be output, external signal or time. βilθlk is the normal weight function and fulfills the following conditions:(4)wilθlk=∏j=1sMijlθlk,wilθlk≥0,βilθlk=wilθlk∑i=1rlwilθlk0≤βilθlk≤1and∑i=1rlβilθlk=1.Suppose that some state variables of system (1) are unmeasurable and the interconnection term glx is arbitrary and the upper bound constraints of glx are unknown as well. However, these state variables are necessary to monitor the system. It should be noted that the interconnection term glx will make the observer design much more difficult. In addition, the interconnection terms glx are arbitrary and the upper bounds of glx are not known; therefore, the methods in [21, 22, 25–29] are unable to employ to design an observer for this system. Therefore, the objective of this study is to propose a new approach to synthesize a decentralized observer for estimating the unknown states and completely eliminate the influences of the interconnection terms.Useful LemmasLemma 3.
Considering two matricesΨ and Λ with appropriate dimensions, the following condition holds for any positive scalar ϵ:(5)ΨTΛ+ΛTΨ≤ϵΨTΨ+ϵ-1ΛTΛLemma 4 (see [35]).
Consider two matricesX and Φ with the compatible dimensions. The following two statements are equivalent.
( 1 ) There exists a positive symmetric matrix P such that (6)XTPX-Φ<0
( 2 ) There exists a positive symmetric matrix P and matrix Z such that(7)-Φ∗ZX-Z-ZT+P<0


## 3. Decentralized Observer Design

In this section, a polynomial observer is synthesized for system (3a) and (3b) to estimate the unknown states and eliminate the effects of the interconnection terms.Let us consider the observer for system (3a) and (3b) as follows:(8a)zlk+1=∑i=1rlβilθlkNilζlkzlk+Gilζlkulk+Lilζlkylk(8b)xl^k=zlk-Jlylkwhere zlk is the vector of the state variables of the observer (8a) and (8b). The polynomial matrices Nilζlk,Gilζlk,Lilζlk and constant matrix Jl are the observer gains which will be determined in the next steps. It should be noted that, in this study, premise variables θl(k) are measured; hence, we can use these premise variables for the observer (8a) and (8b) and do not need to estimate them.Let us define estimation errors of system (3a) and (3b) with observer (8a) and (8b) (9)elk=x^lk-xlkSubstituting (3b) and (8b) into (9) yields(10)elk=zlk-Jlylk-xlk=zlk-I+JlClxlkDenote(11)Hl=I+JlClthen (9) is rewritten as follows:(12)elk=zlk-HlxlkThe dynamic equation of estimation error is obtained(13)elk+1=zlk+1-Hlxlk+1Combining (3a), (8a), and (13) yields(14)elk+1=∑i=1rlβilθlkNilζlkxlk+Gilζlkulk+Lilζlkylk-Hl∑i=1rlβilθlkAilζlkxlk+Bilζlkulk+Tlϑlx=∑i=1rlβilθlkNilζlkelk+NilζlkHl-HlAilζlk+LilζlkClxlk+Gilζlk-HlBilζlkulk+HlTlϑlxTheorem 5.
Based on Assumption1, the estimation error (9) of system (3a) and (3b) with observer (8a) and (8b) converges to zero asymptotically if there exist the polynomial matrices Nilζlk, Gilζlk, Lilζlk, constant matrix Jl, and symmetric matrices Pil such that the following conditions hold:(15)NilζlkHl-HlAilζlk+LilζlkCl=0(16)Gilζlk-HlBilζlk=0(17)HlTl=0(18)v1TPil-ε1ilIv1isSOS(19)-v2TΞilζlk+ε2ilζlkIv2isSOSwhere(20)Ξilζlk=NilζlkTPmlNilζlk-PilHl=I+JlClv1 and v2 are the vectors which are independent on ζlk, ε1il is a positive scalar, ε2ilζlk>0 with ζlk≠0.Proof.
If the conditions (15)-(17) hold, then (14) becomes(21)elk+1=∑i=1rlβilθlkNilζlkelkSelect fuzzy Lyapunov function(22)Vk=∑l=1NVlk=∑l=1N∑i=1rlβilθlkelTkPilelk.From condition (18) of Theorem 5, it implies that Pil is a positive matrix. Thus, Vk>0. From (22), we have(23)∆Vk=Vk+1-VkCombining (22) and (23) obtains(24)∆Vk=∑l=1N∑i=1rlβilθlk+1elTk+1Pilelk+1-∑i=1rlβilθlkelTkPilelkSubstituting (21) into (24) yields (25)∆Vk=∑l=1N∑i=1rl∑j=1rl∑m=1rlβilθlkβjlθlkβmlθlk+1elTkNjlζlkTPmlNilζlk-PilelkDefine(26)Φmlζlk=∑i=1rl∑j=1rlβilθlkβjlθlkelTkNjlζlkTPmlNilζlk-PilelkFrom (25) and (26), one obtains(27)∆Vk=∑l=1N∑m=1rlβmlθlk+1ΦmlζlkFrom (26), it infers that (28)Φmlζlk=∑i=1rlβilθlk2·elTkNilζlkTPmlNilζlk-Pilelk+∑i=1rl∑j=1i<jrlβilθlkβjlθlkelTk·NilζlkTPmlNjlζlk-Pil+NjlζlkTPmlNilζlk-PjlelkApplying Lemma 3 for (28) obtains(29)Φmlζlk≤∑i=1rlβilθlk2·elTkNilζlkTPmlNilζlk-Pilelk+∑i=1rl∑j=1i<jrlβilθlkβjlθlkelTk·NilζlkTPmlNilζlk-Pil+NjlζlkTPmlNjlζlk-PjlelkIt is easily seen that(30)∑i=1rlβilθlk2elTkNilζlkTPmlNilζlk-Pilelk+∑i=1rl∑j=1i<jrlβilθlkβjlθlkelTk·NilζlkTPmlNilζlk-Pil+NjlζlkTPmlNjlζlk-Pjlelk=∑i=1rlβilθlkelTkNilζlkTPmlNilζlk-Pilelktherefore,(31)Φmlζlk≤∑i=1rlβilθlkelTkNilζlkTPmlNilζlk-PilelkFrom (31), it leads to (32)∆Vk=∑l=1N∑i=1rlβilθlk+1Φmlζlk≤∑i=1rlβilθlkelTkNilζlkTPmlNilζlk-PilelkIt is obvious that if condition (19) holds, then ∆Vk<0. It means that the estimation errors (9) approach zero asymptotically. The proof is completed.However, condition (19) is a polynomial Bilinear Matrix Inequality which makes it very difficult to obtain the results in Matlab by using SOS tool. Thus, a transformation to convert (19) into the Polynomial Linear Matrix Inequality is needed.Theorem 6.
Based on Assumption1, the estimated state variables of system (3a) and (3b) with the observer (8a) and (8b) approach the real states asymptotically if there exist the polynomial matrices Nilζlk, Gilζlk, Lilζlk, Zlζlk,Jl, and symmetric matrices Pil such that the following conditions satisfy (33)v1TPil-ε1ilIv1isSOS(34)-v2T-Pml∗Ωilζlk-Zlζlk-ZlζlkT+Pml+ε2ilζlkIv2isSOSin which(35)Ωilζlk=ZlζlkI+UlClAilζlk+QlζlkVlClAilζlk-MilζlkCl(36)Qlζlk=ZlζlkYl(37)Milζlk=ZlζlkKilζlk(38)Ul=-TlClTl+(39)Vl=I-ClTlClTl+(40)Hl=I+JlClObserver gains are obtained as follows:(41)Jl=Ul+YlVl(42)Gilζlk=HlBilζlk(43)Lilζlk=KilζlkI+ClJl-HlAilζlkJl(44)Nilζlk=HlAilζlk-KilζlkCl(ClTl)+ is the Moore-Penrose pseudoinverse of (ClTl), i=1,2,…,rl, l=1,2,…,N.Proof.
From condition (19) of Theorem 5, it infers that (45)Ξilζlk=NilζlkTPmlNilζlk-Pil<0Employing Lemma 4 for condition (45) of Theorem 5 with a slack variable Zlζlk yields(46)-Pml∗ZlζlkNilζlk-Zlζlk-ZlζlkT+Pml<0From condition (17), we have(47)JlClTl=-TlThe general solution of (47) is (48)Jl=Ul+YlVlwhere Ul=-TlClTl+Vl=I-ClTlClTl+(ClTl)+=((ClTl)TClTl)-1(ClTl)T is the Moore-Penrose pseudoinverse of (ClTl)Yl is an arbitrary matrix with the compatible dimension.
Define(49)Kilζlk=Lilζlk+NilζlkJlCombining (15) and (49) yields(50)Nilζlk=HlAilζlk-KilζlkCl(51)Lilζlk=KilζlkI+ClJl-HlAilζlkJlFrom (48) and (50), one obtains(52)Nilζlk=I+Ul+YlVlClAilζlk-KilζlkClSubstituting (52) into (46) yields(53)-Pil∗ZlζlkI+Ul+YlVlClAilζlk-KilζlkCl-Zlζlk-ZlζlkT+Pil<0Let us denote(54)Qlζlk=ZlζlkYl(55)Milζlk=ZlζlkKilζlkFrom (53), (54), and (55), we obtain(56)-Pil∗Ωilζlk-Zlζlk-ZlζlkT+Pil<0where(57)Ωilζlk=ZlζlkI+UlClAilζlk+QlζlkVlClAilζlk-MilζlkClFrom (55), it is seen that condition (19) of Theorem 5 has been transformed into PLMI (34) of Theorem 6 which can be solved easily in Matlab. Hence, the proof is completed.For the sake of understanding easily, the procedure for observer design is briefed as follows.Step 1.
ComputeUl and Vl from (38) and (39).Step 2.
SubstituteAilζlk, Ul, Vl, and Cl into (34), then solving this SOS by using SOS tool of Matlab obtains Pml, Zlζlk, Qlζlk, and Milζlk.Step 3.
CalculateYlandKilζlk from (36) and (37), respectively. After that, Hl is determined by (40).Step 4.
The observer gainsJl,Gilζlk, Lilζlk,andNilζlk, are computed from (41)-(44).Step 5.
With observer gains in Step4, observer (8a) and (8b) is constructed.Remark 7.
In [22, 25], the interconnection terms are required to satisfy the bounded constraint; however, the considered interconnection parts in this work are arbitrary and do not need to satisfy any bounded constraint; therefore, the methods in [22, 25] failed to design the observer for these systems while our method is still applicable. In addition, the method in [21, 26–29] is applicable when the interconnection terms can be linearized under the linear framework or T-S fuzzy model. However, in practice, in some systems, it is very difficult to transform the interconnection terms to the linear or T-S fuzzy model forms; it cause the method in [21, 26–29] to be unable to apply for these systems. Even if these interconnection parts can be linearized under the linear form or T-S fuzzy model as in [21, 26–29], it also has the following drawbacks: (1) the modeling error between the original large-scale nonlinear system and the modeled system will increase and (2) the number of the fuzzy rules, the complexity, and computational load significantly grows up.Remark 8.
Regarding the equality conditions (15)-(17) of Theorem 5, they have been integrated to transform the PBMI (Polynomial Bilinear Matrix Inequality) (19) into the PLMI (Polynomial Linear Matrix Inequality) (34) of Theorem 6. Then by using SOS tool of Matlab, it is easy to find the feasible solution of PLMI (34) and these obtained results are guaranteed satisfying the equality conditions (15)-(17) as well.


## 4. Illustrative Example

Take into account a discrete-time nonlinear large-scale system with three interconnection subsystems as follows:S u b s y s t e m   1 (58a) x 11 k + 1 = - 0.1 x 11 k + 0.1 x 11 k x 21 k cos ⁡ x 11 k - 0.2 x 21 k cos ⁡ x 22 k - 0.2 x 13 k sin ⁡ x 23 k + u 1 k x 21 k + 1 = 0.1 x 11 k - 0.1 x 11 2 k x 21 k - x 21 k cos ⁡ x 22 k - x 13 k sin ⁡ x 23 k + u 1 k y 1 k = x 11 kS u b s y s t e m   2 (58b) x 12 k + 1 = - 0.1 x 12 k - 0.1 x 12 2 k x 22 k + 0.5 x 11 k × sin ⁡ x 12 k + 0.5 x 23 k cos ⁡ x 13 k + u 2 k x 22 k + 1 = - 0.1 x 12 k cos ⁡ x 12 k - 0.5 x 22 k 2 x 11 k sin ⁡ x 12 k + 2 x 23 k c o s x 13 k + u 2 k y 2 k = x 12 kS u b s y s t e m   3 (58c) x 13 k + 1 = - 0.2 x 13 k + 0.1 x 13 k + 1 x 23 k cos ⁡ x 13 k + 0.2 x 11 k sin ⁡ x 21 k - 0.2 x 22 k cos ⁡ x 12 k + u 3 k x 23 k + 1 = - 0.5 x 13 k cos ⁡ x 13 k - 0.1 x 13 2 k x 23 k + x 11 k sin ⁡ x 21 k - x 22 k cos ⁡ x 12 k + u 3 k y 3 k = x 13 kThe interconnection terms:(59)g1x=-0.2x21kcos⁡x22k-0.2x13ksin⁡x23k-x21kcos⁡x22k-x13ksin⁡x23kg2x=0.5x11ksin⁡x12k+0.5x23kcos⁡x13k2x11ksin⁡x12k+2x23kcos⁡x13kg3x=0.2x11ksin⁡x21k-0.2x22kcos⁡x12kx11ksin⁡x21k-x22kcos⁡x12kIn this example, rewriting system (58a), (58b), and (58c) under the framework of system (3a) and (3b) with ζlk=ylk yields(60)xlk+1=∑i=1rlβilθlkAilylxlk+Bilylulk+Tlϑlxylk=Clxlk,l=1,2,3where (61)x1k=x11kx21k;x2k=x12kx22k;x3k=x13kx23k;A11y1=-0.10.1y10.1-0.1y12,B11y1=B21y1=11A21y1=-0.1-0.1y10.1-0.1y12,C1=10,T1=0.21.ϑ1x=-x21kcos⁡x22k-x13ksin⁡x23kA12y2=-0.1-0.1y22-0.1-0.5,A22y2=-0.1-0.1y220.1-0.5,B12y2=B22y2=11C2=10,T2=0.52.ϑ2x=x11ksin⁡x12k+x23kcos⁡x13kA13y3=-0.20.1y3+0.1-0.5-0.1y32,A23y3=-0.2-0.1y3-0.10.5-0.1y32,B13y3=B23y3=11C3=10,T3=0.21.ϑ3x=x11ksin⁡x21k-x22kcos⁡x12kApplying the method in Section 3 and solving the conditions of Theorems 5 and 6 yields(62)N11y1=-0.229∗10-15y12+0.298∗10-12y1-0.815∗10-700.004y12+95.38y1+2.86-0.1y12-0.5y1N21y1=-0.225∗10-15y12+0.308∗10-12y1-0.815∗10-700.004y12+95.36y1+2.59-0.1y12+0.5y1N12y2=0.10∗10-16y22+0.55∗10-7y2+0.154∗10-300.769∗10-10y22-0.55∗10-7y2+165.310.4y22-0.5N22y2=0.30∗10-17y22+0.105∗10-15y2+0.133∗10-30-0.934∗10-9y22-0.124∗10-6y2+213.00.4y22-0.5N13y2=-0.91∗10-14y32+0.619∗10-8y3+0.127∗10-200.0025y32-43.27y3-302.950.2y32-0.5y3-0.5N23y2=-0.406∗10-13y32+0.904∗10-8y3+0.263∗10-200.007y32+31.17y3-455.25-0.2y32+0.5y3+0.5J1=-1-5,J2=-1-4,J3=-1-5,G11=0-4,G21=0-4,G12=0-3,G22=0-3,G13=0-4,G23=0-4L11y1=0-0.5y12-2.5y1+0.6,L21y1=0-0.5y12+2.5y1+0.6L12y2=01.6y22-0.17,L22y2=01.6y22-1.5L13y3=0-0.5y32-2.5y3-2,L23y3=0-0.5y32+2.5y3+4Simulating the system in Matlab, we obtain the results in Figures1–6.Figure 1
Real statex11(k), estimated state x^11(k), and estimation error e11(k).Figure 2
Real statex21(k), estimated state x^21(k), and estimation error e21(k).Figure 3
Real statex12(k), estimated state x^12(k), and estimation error e12(k).Figure 4
Real statex22(k), estimated state x^22(k), and estimation error e22(k).Figure 5
Real statex13(k), estimated state x^13(k), and estimation error e13(k).Figure 6
Real statex23(k), estimated state x^23(k), and estimation error e23(k).The real statesx11(k), x21(k), estimated states x^11(k), x^21(k), and estimation errors e11k,e21(k) of the subsystem 1 are shown in Figures 1 and 2; Figures 3 and 4 illustrated the real states x12(k), x22(k), estimated states x^12(k), x^22(k), and estimation errors e12k,e22(k) of the subsystem 2. Finally, the real states x13(k), x23(k), estimated states x^13(k), x^23(k), and estimation errors e13k,e23(k) are demonstrated in Figures 5 and 6, respectively. From these simulation results, it is obvious that all estimated states approach the real states and the estimation errors converge to zero asymptotically. Hence, we can conclude that the proposed method in this paper is successful to design an observer for the large-scale system (58a), (58b), and (58c).Remark 9.
It should be noted that if we employ the T-S fuzzy model to represent system (58a), (58b), and (58c), then the number of fuzzy rules in each subsystem will double making the number of the conditions for observer design exponentially grows up. Moreover, it is seen that we do not know the upper bounds of the interconnection terms g1x, g2x,g3x in this example; thus, the methods in [22, 25] are unable to apply for this system. In addition, because the interconnection terms glx,l=1,2,3 of system (58a), (58b), and (58c) in this example contain the nonlinear function of states from other interconnection systems, it is not possible to model system (58a), (58b), and (58c) under the framework of the T-S fuzzy model as [8–15, 21, 26–29]. Hence, the methods in [8–15, 21, 26–29] failed to apply to design observer for system (58a), (58b), and (58c).


## 5. Conclusions

A new approach based on the unknown input method to synthesize an observer for the discrete-time large-scale polynomial nonlinear system is studied in this paper. The discrete-time large-scale nonlinear system is modeled in terms of a group of discrete-time polynomial T-S fuzzy subsystems to decrease the number of fuzzy rules. The interconnection parts in this system are arbitrary and they do not require the bounded constraints. By employing the unknown input method, the influences of the interconnection terms are eliminated completely and the unmeasurable state variables of the system are estimated asymptotically. The conditions represented under the form of SOS for synthesizing observer are derived in the two theorems. An illustrative example is presented to demonstrate the success and effectiveness of the proposed method.



---

*Source: 2192768-2019-04-22.xml*