Observer Design for Takagi–Sugeno Fuzzy Systems with Unmeasurable Premise Variables Based on Differential Mean Value Theorem †

: In this work, we present the design of an observer for Takagi–Sugeno fuzzy systems with unmeasurable premise variables. Moving away from Lipschitz-based and L 2 attenuation-based methods—which fall short in eliminating the mismatching terms in the estimation error dynamics— we leverage the differential mean value theorem. This approach not only removes these terms but also streamlines the factorization of the estimation error dynamics, making it directly proportional to the estimation error. To ensure the asymptotic convergence of the estimation error, we apply the second Lyapunov theorem, which provides sufficient stability conditions described as linear matrix inequalities. A numerical example applied on a three-tank hydraulic system is presented to demonstrate the observer’s effectiveness.


Introduction
In the industrial sector, cost-effectiveness is paramount.A key strategy to achieve this is using observers to reduce the need for expensive sensors.The Luenberger observer [1] has paved the way for numerous advancements.Takagi-Sugeno (TS) fuzzy systems, which represent nonlinear systems as a weighted sum of linear ones [2], have provided valuable tools for understanding complex dynamics.Using the sector nonlinearity approach, these systems can be accurately described [3].Based on variables in their weight functions, they're classified into measurable premise variables (MPV) and unmeasurable premise variables (UPVs), with the latter being a primary research focus because it represents the largest category of systems.
Designing observers for systems equipped with UPV tends to be more intricate than their measurable counterparts.These complexities stem predominantly from the mismatching terms in the error dynamics.To address such hurdles, the scientific community has forwarded various techniques.Initially, the Lipschitz-based method emerges as a straightforward solution [4,5], yet stumbles when the nonlinear system's Lipschitz constant exceeds an admissible value, thereby introducing pronounced conservatism in Linear Matrix Inequality (LMI) constraints.An alternative, the L 2 -attenuation-based approach [6,7], focuses on minimizing the aforementioned mismatches.Though typically less conservative than the Lipschitz method, there remain instances where the attenuation level is not minimal enough to be accepted, even if the simulation works.This paves the way for the mean value theorem (MVT) method [8,9], allowing the factorization of the estimation error dynamics, which leads to making it proportional completely to the estimation error.Hence, the mismatching terms disappeared from the estimation error dynamics.
The observer based on the MVT, leading to an exact transformation of error dynamics into an LPV system, was introduced in [10].Subsequent works, such as [11,12], represent dynamic error using TS representation.This theorem's applications span various studies: real-time motor control in [13], robust H ∞ control for motors in [14], sensorless control for a PMSM in [15], and automotive slip angle estimation in [12].A notable advancement is in [16], where line integral Lyapunov functions reduce conservatism in MVT observers.
In this paper, we introduce an observer design based on the mean value theorem and validate its efficacy via a simulation on a three-tank hydraulic system.The paper's structure is as follows: Section 2 delves into the TS fuzzy representation and the mean value theorem.Section 3 details the observer design, Section 4 showcases the simulation results, and Section 5 concludes with potential directions for future research.

Preliminaries 2.1. Takagi-Sugeno Fuzzy Representation
Let us consider the following nonlinear system: where x(t) ∈ R n is the state vector, u(t) ∈ R n u is the input vector, and y(t) ∈ R n y represents the output vector.
The TS model reformulates the nonlinear system (1) using a convex combination of linear sub-models.Each i th sub-model follows the given fuzzy rule: where ξ is premise variable, l is the number of premise variables, and M ij is the membership function of the i th fuzzy rule corresponding to the j th premise variable.A i ∈ R n×n , B i ∈ R n×n u , and C i ∈ R n y ×n are known matrices.The global TS fuzzy representation of the nonlinear system (1) is described as follows: where µ i (x(t)) is the weighting functions described by µ i (x(t)) = ∏ l j=1 M ij ξ j and verifies the convex sum property

Differential Mean Value Theorem [10]
Let Φ(x) : R n => R n be a differentiable vector function described as follows: where the set E n is the canonical basis of the vectorial space R n for all n ≥ 1 given by: Let a, b ∈ R n .Then, there are constant vectors z 1 , ..., z n ∈ (a, b), z i ̸ = a, z i ̸ = b for i = 1, ..., n such that the mean value theorem ensures the following relation: Applying sector nonlinearity allows for rewriting the above equation using TS transformation: where H i represents the sub-models of the nonlinear term ∑ n i=1 ∑ n j=1 e n (i)e n (j) T ∂Φ i (zj) , and ) is its weighting function.

Mean Value Theorem Observer-Based Design
Let us consider the following observer for fuzzy system (3): The following theorem provides sufficient conditions described as LMI to ensure the asymptotic convergence of the error dynamic: Theorem 1.The estimation error converges asymptotically toward zero with decay rate α if there exist matrices P = P T ∈ R n x ×n x > 0 and M ∈ R n x ×n y such that the following LMI holds ∀i = 1, ..., q: The observer gain is given by L = P −1 M.
Proof.The dynamics of the estimation error e(t) = x(t) − x(t) are given as follows: Using the mean value theorem on the terms Φ 1 and Φ 2 , the error dynamics become the following: To study the stability of the error dynamics, the quadratic Lyapunov function is used: The derivative of V(t) with respect to t is as follows: To improve the performance of the estimation, the following decay rate is used: By substituting ( 12) and ( 13) in ( 14), the following inequality is obtained: The inequality (15) is not linear due to the product of the variables P and L. However, applying the change in variable M = PL provides a solution for achieving the linear stability conditions outlined in Theorem 1.
In order to solve the inequalities outlined in Theorem 1, we will utilize the Yalmip Toolbox, a well-regarded modeling language in MATLAB designed for formulating opti-mization problems.This toolbox seamlessly integrates with a variety of solvers, including "Mosek", "SDPT3", and "LMILAB", all recognized for their efficiency in handling LMIs.
In Figure 1, the entire procedure is depicted through a graphical representation:

Initialize system parameters
Resolve linear matrix inequalities in Theorem 1

Dynamic Model of the System
The three-tank hydraulic system illustrated in Figure 2, based on Guzman's design [17], features three tanks with equal cross-sectional areas S, connected by pipes with areas S p1,2,3 .Water from a reservoir fills the first and second tanks via pumps P 1 and P 2 , with flow rates u 1 and u 2 .Valves in each tank manage water release, and the system ensures water levels in the order x 1 > x 3 > x 2 .Let us define x(t) = x 1 x 2 x 3 T , u(t) = u 1 u 2 T and the premise variables T .Using these definitions, the system can be represented in the following state-space form: where With a gravitational pull of g = 9.8 m/s 2 , the system has discharge coefficients η 13 = η 32 = 0.456 and η 20 = 0.652.The tubes' cross-sectional areas are S p1 = S p3 = 0.5 × 10 −4 m 2 and S p2 = 0.8 × 10 −4 m 2 , with all tanks having S = 154 × 10 −4 m 2 .Given the constraint x 1 > x 3 > x 2 and these parameters, we derive C 1 = 1.0094 × 10 −4 , C 2 = 2.3092 × 10 −4 , and C 3 = 1.0094 × 10 −4 .

Observer Design for Three-Tank Hydraulic System
In order to apply Theorem 1, the matrices A i and C i have to be determined.According to the mean value theorem, Φ 1 (ξ) and its Jacobian ∂Φ 1 (ξ) ∂x can be defined as follows: where the new premise variables ε i (x) are given by ε 1 (x) = and the limits are L ε1 (j) = 2 12 , L ε2 (k) = 2 4 and L ε3 (d) = 2 14 for j, k, d = 1 : 2.
By replacing every premise variable by its respective limits in a loop, the matrices A i can be obtained as follows for (i = 1, . . ., 8): According to the output equation, which is linear, Φ 2 (ξ) = Cx(t); hence,

Conclusions
This paper introduces a design for the Takagi-Sugeno observer using the mean value theorem, bypassing the commonly used Lipschitz assumption and the L 2 attenuationbased method.While our method adeptly tracks the system's states, it uniquely applies the MVT to systems with nonlinear outputs, contrasting prior works such as [8,9,13,14] that targeted only linear output systems.However, the proposed approach does have limitations, particularly for systems with unknown inputs, directing our future research ambitions.We aim to explore observer designs accommodating unknown inputs and to investigate emerging control system methodologies, like the adaptive fuzzy control for pneumatic active suspensions, as seen in [18].This paves the way for the Takagi-Sugeno observer's broader applications in complex systems.

Figure 1 .
Figure 1.Overall schematic diagram of observer design and implementation.