Next Article in Journal
Bimetallic Mo–Fe Co-Catalyst-Based Nano-Carbon Impregnated on PAC for Optimum Super-Hydrophobicity
Next Article in Special Issue
Dynamic Response of the Newton Voigt–Kelvin Modelled Linear Viscoelastic Systems at Harmonic Actions
Previous Article in Journal
Detection of Epileptic Seizures Using Wavelet Transform, Peak Extraction and PSR from EEG Signals
Previous Article in Special Issue
Rational Transfer Function Model for a Double-Pipe Parallel-Flow Heat Exchanger
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Problem of Invariance in Nonlinear Discrete-Time Dynamic Systems

Department of Automation and Control, Far Eastern Federal University, 690095 Vladivostok, Russia
Submission received: 2 May 2020 / Revised: 27 June 2020 / Accepted: 15 July 2020 / Published: 28 July 2020
(This article belongs to the Special Issue Symmetry in Dynamic Systems)

Abstract

:
The paper considers the problem of invariance with respect to the unknown input for discrete-time nonlinear dynamic systems. To solve the problem, the algebraic approaches, called algebra of functions and logic–dynamic approach, are used. Such approaches assume that description of the system may contain non-differentiable functions. Necessary and sufficient conditions of solvability the problem are obtained. Moreover, procedures which find the appropriate functions and matrices are developed. Some applications of such invariance in fault detection and isolation, disturbance decoupling problem, and fault-tolerant control are considered.

1. Introduction and Problem Statement

The problem of invariance with respect (IWR) to the unknown input in nonlinear dynamic systems is of wide theoretical and practical applications, in particular, in fault diagnosis [1,2,3], fault-tolerant control [4,5] and disturbance decoupling [6,7,8,9]. The problem of invariance in these branches is solved mainly for continuous-time systems based on the methods of linear algebra [10] and differential geometry developed in [11]. Such tools are of limited field of application in practice since actual systems may contain different non-differentiable function such as Coulomb friction, saturation, hysteresis and dead zone.
It is well-known from literature that the extensions of the differential geometric tools are not well developed for discrete-time systems in comparison with the continuous time case, see [12,13,14] (disturbance decoupling problem solution), [15] (transitivity and accessibility problems), [16] (non-interactive control). Such extensions are rather complicated and assume that systems under consideration are described by analytical functions. For this reason, such extensions are of limited field of application as well.
In this paper, we investigate the problem of IWR to the unknown input in discrete-time nonlinear dynamic systems. To overcome the difficulty with the above-mentioned practical applications, it is suggested to solve the problem of invariance based on algebra of functions and logic–dynamic approach (LDA). The advantage of these approaches, if compared to well-known linear algebraic [10] and differential geometric methods [10,17], is that they are applicable to non-differentiable systems.
The algebra of functions was developed in [18,19] on the basis of the pair algebra of partitions [20] elaborated for finite automata study. The algebra of functions is intended mainly for discrete-time nonlinear dynamic systems with non-differentiable nonlinearities. Since the algebra of functions is not well-known and demands complex analytical calculations, the webMathematica based software was developed [21].
The LDA was suggested in [22] for solving different problems in dynamic system theory. The advantage of the LDA is that under some limitations on the original system and a class of possible solutions it uses only methods of linear algebra. Moreover, the LDA can be applied to both continuous-time and discrete-time systems with non-differentiable nonlinearities. Shortcoming of the LDA is that functions transforming the system under consideration are required to be linear, see the relation φ ( x ( t ) ) = Φ x ( t ) in Section 4. Such a requirement imposes definite limitations on application of the LDA.
The contributions of this paper can be summarized as follows: (1) some relations describing the problem of invariance with respect to the unknown input in terms of the algebra of functions and the LDA are obtained; (2) links between the algebra of functions relations and the LDA relations are established since such links are interesting and useful for theory and practice, but not well known.
Consider a discrete-time nonlinear dynamic system described by the equations
x ( t + 1 ) = f ( x ( t ) , u ( t ) , w ( t ) ) ,     y ( t ) = h ( x ( t ) ) .
Here x X R n , u U R m and y Y R l are vectors of state, input and output; w ( t ) W R p is the unknown input; f and h are nonlinear functions. Note that f may be non-differentiable function.
The problem of invariance with respect to the unknown input is stated as follows: find a vector function x 0 = φ ( x ) , transforming system (1) into the system
x 0 ( t + 1 ) = f 0 ( x 0 ( t ) , y ( t ) , u ( t ) ) ,     y 0 ( t ) = h 0 ( x 0 ( t ) )
which does not depend on the unknown input. Here x 0 R n 0 , n 0 < n , is the state vector, f 0 and h 0 are some functions to be determined. To solve this problem for system (1) and to develop a procedure of transformation, the algebra of functions will be used.

2. Algebra of Functions

The main definitions and concepts used in this paper are as follows [19,23,24]. Let S be a set of vector functions with the domain X. The elements of algebra of functions are vector functions on S with the following relations, operations and operators: (1) relation of partial preorder ≤, (2) two binary operations × and , (3) binary relation Δ, (4) two operators m and M.
Given α , β S , one says that α β if a function γ exists such that β ( x ) = γ ( α ( x ) ) for all x X . When α β and β α , the functions α and β are called equivalent, denoted as α β .
The relation is reflexive, symmetric and transitive; therefore, this relation divides the set S into the classes of equivalence. Denote by S \ the set of all classes of equivalence; then the relation ≤ is partial order on this set. It can be shown that S \ is a lattice where every two elements α and β have a unique supremum (least upper bound) sup ( α , β ) and a unique infimum (greatest lower bound) inf ( α , β ) . As is customary, we will operate not with sup ( α , β ) and inf ( α , β ) but with two binary operations and × respectively.
In the simple cases, the definition α β = sup ( α , β ) is used to compute α β . The rule for operation × is simple:
( α × β ) ( x ) = ( α ( x )   β ( x ) ) T
Note that two special vector functions 0 and 1 exist; they correspond to the identity and constant functions, respectively, in the sense that for every vector function α S , 0 α 1 .
Consider illustrative example. Let α ( x ) = ( x 1   x 3 + x 4 ) T and β ( x ) = ( x 2   x 3   x 4 ) T , then
( α × β ) ( x ) = ( x 1   x 2   x 3   x 4 ) T = 0   , ( α β ) ( x ) = x 3 + x 4 .
Clearly, ( α × β ) ( x ) α ( x ) ( α β ) ( x ) and ( α × β ) ( x ) β ( x ) ( α β ) ( x ) ; the functions α ( x ) and β ( x ) are incomparable.
Given α , β S , then ( α , β ) Δ if a function f 0 exists such that for all x , u , w X × U × W β ( f ( x , u , w ) ) = f 0 ( α ( x ) , u , w ) . When ( α , β ) Δ , one says that α and β form an ordered pair. Binary relation Δ is used to define the operators m and M.
Operator m ( α ) is a function in S, satisfying the following conditions: (i) ( α , m ( α ) ) Δ , (ii) if ( α , β ) Δ , then m ( α ) β .
Operator M ( β ) is a function in S, satisfying the following conditions: (i) ( M ( β ) , β ) Δ Δ, (ii) if ( α , β ) Δ , then α M ( β ) .
It follows from the last definitions that given α , m ( α ) is the minimal function, that forms a pair with α and given β , M ( β ) is the maximal function, that forms a pair with β .
Lemma 1
[19,23].Let α and β be some functions from S. Then
(i). α β α α × β β α β ,
(ii). α β M ( α ) M ( β ) m ( α ) m ( β ) ,
(iii). M ( α × β ) M ( α ) × M ( β ) , m ( α β ) m ( α ) m ( β ) .
Computation of the operators m and M. It is known from [19,23] that there exists the function γ satisfying the condition ( α × u ( t ) ) f γ ( f ) ; define m ( α ) γ . When β ( f ( x ( t ) , u ( t ) ) ) can be transformed into the form
β ( f ( x ( t ) , u ( t ) ) ) = i = 1 d a i ( x ( t ) ) b i ( u ( t ) ) ,
where a 1 ( x ( t ) ) , a 2 ( x ( t ) ) , …, a d ( x ( t ) ) are arbitrary vector functions and b 1 ( u ( t ) ) , b 2 ( u ( t ) ) , …, b d ( u ( t ) ) are linearly independent functions, then M ( β ) = a 1 × a 2 × × a d .

3. Problem Solution

To solve the problem of IWR to the unknown input for system (1), we at first find a vector function φ 0 with maximal number of functionally independent components such that the function φ 0 ( f ( x , u , w ) ) is independent of the unknown function w ( t ) . The function φ 0 can be obtained by heuristic methods. Actually, this function provides some combination of the function f ( x , u , w ) to achieve independence of the unknown input.
One says that the function φ is ( h , f ) -invariant if φ ( f ( x , u , w ) ) = f 0 ( φ ( x ) , h ( x ) , u , w ) for some vector function f 0 . It is known [19,23] that φ is ( h , f ) -invariant if and only if φ × h M ( φ ) or m ( φ × h ) φ .
Theorem 1.
System (2) is IWR to the unknown input w ( t ) if and only if ( h , f ) -invariant function φ exist such that
φ 0 φ .
Proof. 
Let φ be ( h , f ) -invariant function satisfying (3). By definition of ( h , f ) -invariance, φ ( f ( x , u , w ) ) = f 0 ( φ ( x ) , h ( x ) , u , w ) . Since φ 0 φ , then γ ( φ 0 ) = φ for some function γ , therefore φ ( f ( x , u , w ) ) = γ ( φ 0 ( f ( x , u , w ) ) ) . By definition, φ 0 ( f ( x , u , w ) ) is independent of the unknown function w ( t ) ; as a result, the function f 0 ( φ ( x ) , h ( x ) , u , w ) = γ ( φ 0 ( f ( x , u , w ) ) ) is independent of w ( t ) as well. On the other hand, let φ ( f ( x , u , w ) ) is independent of w ( t ) . Since the function φ 0 has the same property and has maximal number of functionally independent components, then γ ( φ 0 ) = φ for some function γ or φ 0 φ . Since φ is ( h , f ) -invariant, then φ ( f ( x , u , w ) ) = f 0 ( φ ( x ) , h ( x ) , u , w ) for some function f 0 . □
To construct system (2) of maximal dimension, the function φ should be minimal in terms of the partial preorder relation ≤. Such a function can be obtained as follows.
Theorem 2
[19,23].Given φ 0 , compute recursively for i 0 , based on the formula
φ i + 1 = φ i m ( φ i × h ) ,
the sequence of vector functions φ 0 φ 1 The sequence converges in a finite number of steps, since if φ i φ i 1 , the number of components of the function φ i is less than that of the function φ i 1 , i = 1 , 2 , This means that there exists a finite k such that φ k + 1 φ k . The function φ : = φ k is minimal satisfying the condition φ 0 φ .
Hence, to solve the problem of IWR to the unknown input, one has to find the function φ 0 and then to use the recursive (4).
When the problem of fault detection and isolation is considered, one needs to generate co-called residual r ( t ) as a mismatch between system (1) behavior and the model (2) behavior, based on the outputs y ( t ) = h ( x ( t ) ) and y 0 ( t ) = h 0 ( x 0 ( t ) ) , respectively. Such a mismatch is presented in the form r ( t ) = ρ ( y ( t ) ) y 0 ( t ) for some function ρ . When faults are absent, r ( t ) = 0 or ρ ( y ( t ) ) = y 0 ( t ) . It can be shown that the last equality is equivalent to the functional relation
ρ ( h ) = h 0 ( φ ) 1 .
Theorem 3.
If h φ 1 , the relation (5) is true for some nontrivial functions ρ and h 0 .
Proof. 
Let h φ 1 , then h φ h and h φ φ by the definition of operation . By definition of the partial order relation ≤, the nontrivial functions ρ and h 0 exist such that ρ ( h ) = h φ = h 0 ( φ ) 1 . □
Given ρ , to construct system (2) of minimal dimension, satisfying the conditions (3) and (5), ( h , f ) -invariant function φ should be maximal in terms of the partial preorder relation ≤ Such a function can be obtained based on the following Algorithm.
Algorithm
Step 1. Set β 0 : = ρ ( h ) and i : = 0 .
Step 2. Compute the function γ i = M ( β i ) .
Step 3. If the components of the vector function γ i can be expressed in terms of the function h × β 0 × × β i , then go to Step 5. Otherwise, go to Step 4.
Step 4. Find the vector function β i + 1 with minimal number of components such that h × β 0 × × β i + 1 γ i , set i : = i + 1 and go to Step 2.
Step 5. Define φ : = β 0 × × β i .
Solution of the problem of IWR to the unknown input may be simplified significantly when the function φ is sought in a class of linear functions. Actually, this restricts a set of possible solutions, but allows to solve the problem for systems with non-differential nonlinearities by methods of linear algebra. Note that if the problem of invariance for system (1) has a solution with linear function φ , such a solution can be found by the logic–dynamic approach described below [22].

4. Logic–Dynamic Approach

To use the LDA, one has to present system (1) in the form
x ( t + 1 ) = F x ( t ) + G u ( t ) + Ψ ( x ( t ) , u ( t ) ) + L w ( t ) ,     y ( t ) = H x ( t ) ,
where
Ψ ( x ( t ) , u ( t ) ) = C ( φ 1 ( A 1 x ( t ) , u ( t ) ) φ q ( A q x ( t ) , u ( t ) ) ) ;
matrices F and G are used to describe the linear dynamic part of the system (6); H, H * , C and L are constant matrices, the functions φ 1 , …, φ q may be non-differentiable, A 1 , …, A q are constant matrices. The model (6) can be obtained from the original system (1) by some simple transformations [22]. Specifically, the linear part with the matrices F and G, is separated from the nonlinear part (7) containing non-differentiable functions φ 1 , …, φ q and matrices C, A 1 , …, A q .
By analogy with (2), a system IWR to the unknown input is described by
x 0 ( t + 1 ) = F 0 x 0 ( t ) + G 0 u ( t ) + J y ( t ) + C 0 ( φ 1 ( A 01 z 0 ( t ) , u ( t ) ) φ q ( A 0 q z 0 ( t ) , u ( t ) ) ) ,     y ( t ) = H 0 x 0 ( t ) ,
where x 0 R n 0 , z 0 = ( x 0 T y T ) T , n 0 n , F 0 , G 0 , J , C 0 , A 01 , …, A 0 q are matrices to be determined.
Assuming initially q = 1 , we construct system (8). The LDA, which is used for solving this problem, contains three main steps [22].
Step 1. Remove the term Ψ ( x ( t ) , u ( t ) ) from the original system (6).
Step 2. Solve the problem of IWS to the unknown input for the linear part under some linear restriction. This restriction is necessary to find out whether or not the nonlinear term can be designed based on the linear solution.
Step 3. Supplement the solution, obtained at Step 2, by the transformed nonlinear term.
It should be noted that the LDA can be applied to the continuous-time systems as well as to the discrete-time ones. This is possible due to a linear nature of the solution at Step 1 and a linear nature of the restriction at Step 2. Moreover, note that a transformation at Step 3 does not transform the nonlinear functions φ 1 , …, φ q themselves, but transforms their arguments based on the relation (9) and the matrix C into C 0 .
We will assume that x 0 ( t ) = φ ( x ( t ) ) = Φ x ( t ) for some matrix Φ of maximal rank satisfying the following conditions [22]:
Φ F = F 0 Φ + J H ,    G 0 = Φ G ,    Φ L = 0 .
One can show that the relations C 0 = Φ C and
A = A 0 ( Φ H ) ,
corresponding to the term Ψ ( x ( t ) , u ( t ) ) , are true [22]. Clearly, the last relation is equivalent to
rank ( Φ H ) = rank ( Φ H A ) .
When q > 1 , the matrix A in (9) and (10) is replaces by A i , i = 1 , , q . Note that (10) is precisely restriction which is checked at Step 2.

5. Solvability Conditions of Invariance

Before constructing system (8), it is worth to check whether or not such a system exists. To make the appropriate checking, note that analog of the function φ 0 is the matrix L 0 of maximal rank satisfying the condition L 0 L = 0 . Analog of the condition φ 0 φ is the relation Φ L = 0 which is equivalent to Φ = N L 0 for some matrix N .
Replace the matrix Φ in the equation Φ F = F 0 Φ + J 0 H by Φ = N L 0 and transform the result: N L 0 F F 0 N L 0 J 0 H = 0 . The obtained equation has a nontrivial solution with some matrices F 0 and J 0 when the rows of the matrices L 0 F and ( L 0 T   H T ) T are linearly dependent that is equivalent to the rank condition
rank ( L 0 F L 0 H ) < rank ( L 0 F ) + ran k ( L 0 H ) .
Considering analogously the Equation (8), we get the necessary condition
rank ( L 0 H ) = rank ( L 0 H A ) .
Note that the relations (11) and (12) are necessary solvability conditions, i.e., if one of the conditions (11) and (12) is not true, the system IWR to the unknown input does not exist.

6. Problem Solution

To design system (8), we assume that the matrices F 0 and H 0 are sought in the form
F 0 = ( 0 1 0 0 0 0 1 0 0 0 0 0 ) ,   H 0 = ( 1 0 0 0 ) .
As a result, the equation Φ F = F 0 Φ + J H is replaced by k equations:
Φ i F = Φ i + 1 + J 0 i H ,   i = 1   , ,   n 0 1 ,   Φ n 0 F = J n 0 H
where Φ i and J i denote the i-th rows of the matrices Φ and J , respectively, i = 1   , ,   n 0 , n 0 is the number of the matrix Φ rows.
It was shown in [22] that (13) and the condition Φ L = 0 can be transformed into the single equation
( Φ 1 J 1 J 2 J k ) ( W ( k )   L ( k ) ) = 0 ,   k = 1 , 2 ,
where
W ( k ) = ( F k H F k 1 H ) ,   L ( k ) = ( L F L F k 1 L 0 H L H F k 1 L 0 0 0 ) .
To construct the system of maximal dimension, set k : = n p and check the condition
rank ( W ( k )   L ( k ) ) < l k + n .
When (15) is true, then there exists the row ( Φ 1 J 1 J k ) such that (14) has a solution. Then one calculates the matrix Φ based on (13) and checks the condition (10). If it is satisfied, find the matrix A 0 from (9), set n 0 : = k , G 0 : = Φ G , and C 0 : = Φ C . Thus, system (8) IWR to the unknown input w ( k ) has been constructed.
If (15) is not true, set k : = k 1 and continue checking (15) and (10). If (15) and (10) are not true for all k, then the system IWR to the unknown input does not exist and the problem has no solution. Since the dimension n 0 is maximal, the matrix Φ is an analog of the function φ from Theorem 2.
By analogy with general case, in some applications it is necessary to take into account the output function y 0 ( t ) = H 0 x 0 ( t ) and the requirement y 0 ( t ) = R y ( t ) for all t 0 and some matrix R . It can be shown that this is equivalent to the relation
R H = H 0 Φ ,
which is analog of (5). Replace the matrix Φ in (16) by Φ = N L 0 and obtain the equation R H H 0 N L 0 = 0 which has a nontrivial solution when
rank ( H L 0 ) < rank ( H ) + rank ( L 0 ) .
Note that (17) is necessary solvability condition additional to (11) and (12). When the matrix Φ has been obtained, the matrices R and H 0 can be found from (16).
To construct system (8) of minimal dimension, one takes k : = 1 and checks the condition (15). When (15) is satisfied, one calculates the matrix Φ and checks the condition (10). If it is satisfied, find the matrix A 0 from (9), set n 0 : = k , G 0 : = Φ G and C 0 : = Φ C . If (15) is not true, set k : = k + 1 and continue checking (15) and (10).

7. Applications

The property “invariance with respect to the unknown inputs” has many different practical applications. Consider some of them related to the fault diagnosis, disturbance decoupling and fault-tolerant control. Here, the unknown inputs are interpreted as the disturbances and faults.
In the fault diagnosis process, the residual is generated as a result of mismatch between the original system behavior and the reference model behavior. Then a decision is made by evaluation of this residual. System (2) or (8) in the fault diagnosis process is used as a reference model. The residual r ( t ) is generated in the form
r ( t ) = ρ ( y ( t ) ) y 0 ( t )   or   r ( t ) = R y ( t ) y 0 ( t ) .
Different tools for fault diagnosis have been developed: diagnostic observers, parity relations and identification [1,2,3].
The main goal of fault detection process is to construct system (2) or (8) of minimal dimension IWR to the disturbances such that the relation (5) or (16) holds. Consider the simple practical example of an electric servo–actuator described by the equations
x 1 ( t + 1 ) = x 1 ( t ) + k 1 x 2 ( t ) , x 2 ( t + 1 ) = x 2 ( t ) + k 2 x 3 ( t ) + k 3 s i g n ( x 2 ( t ) ) + w ( t ) , x 3 ( t + 1 ) = k 4 x 2 ( t ) + k 5 x 3 ( t ) + k 6 u ( t ) , y 1 ( t ) = x 1 ( t ) ,    y 2 ( t ) = x 3 ( t ) .
Using the LDA model (6), we obtain
F = ( 1 k 1 0 0 1 k 2 0 k 4 k 5 ) ,   G = ( 0 0 k 6 ) ,   H = ( 1 0 0 0 0 1 ) ,   L = ( 0 1 0 ) ,   C = ( 0 k 3 0 ) ,   A = ( 0   1   0 ) ,   φ = s i g n ( A x ) .
It can be shown that Φ = ( k 4   0   k 1 ) , R = ( k 4 k 1 ) , J = ( k 4 k 1 k 5 ) , G 0 = k 1 k 6 , C 0 = 0 . The diagnostic observer description is given by
x * ( t + 1 ) = k 4 y 1 ( t ) k 1 k 5 y 2 ( t ) k 1 k 6 u ( t ) ,     y * ( t ) = x * ( t ) ,
the residual is generated as follows: r ( t ) = k 4 y 1 ( t ) k 1 y 2 ( t ) y * ( t ) . The observer is invariant with respect to the disturbances w ( t ) and allows to detect the faults in sensors and deviation of the coefficients k 1 ,   k 4 ,   k 5 ,   k 6 from their nominal values.
When the problem of fault isolation is solved, a bank of such systems is constructed where some faults are considered as the unknown inputs for each system which is constructed to be IWR to the disturbances and faults, considered as the unknown inputs and sensitive to other faults. As a result, we obtain selective sensitivity to different faults that allows to develop fault isolation process based on so-called matrix of syndromes [2].
Sliding mode observers are often used to solve the fault identification problem [3,25,26]. In [3,25], such observers are constructed on the basis of the original system and then IWR to the disturbances is achieved. In contrast to this approach, we suggest at first to construct system (8) of minimal dimension IWR to the disturbances and then design sliding mode observer based on this system [26]. This allows to reduce sliding mode observer complexity and relaxed the limitation imposed on the original system.
The disturbance decoupling problem can be stated as follows. The purpose is to find a dynamic measurement feedback in such a way that the output-to-be-controlled y * ( t ) = h * ( x ( t ) ) , for t 0 , of the closed-loop system does not depend on the disturbances (unknown inputs) w ( t ) . This problem for nonlinear control systems has been studied in [6,7,8,9,23]. Except [23] the papers [6,7,8,9] study the continuous-time case, the solvability conditions are provided in the papers [6,7,8] on the basis of differential geometric tools.
To solve the disturbance decoupling problem for the initial system (1) or (6), system (2) or (8) of maximal dimension, IWR to the unknown inputs, are constructed at first under some restriction imposed by the function h * ( x ) . Then this system is transformed into special compensator which generates the control for the initial system. The details for system (1) can be found in [23], for system (6) in [27,28].
Fault-tolerant control allows to meet the design purposes when the faults occur or if impossible, to redefine the attainable design purposes [4,5]. Active approaches in fault-tolerant control are fault accommodation and plant reconfiguration. The purpose of fault accommodation is to find a new control law which can attain the predefined control purposes. In system reconfiguration, either the controller or the faulty plant is reconfigured when the faults occur. Both approaches are based on system (2) or (8) of maximal dimension without the relations (5) and (16). Such a system is considered as a dynamic part of the compensator and then it is supplemented by static part, generating a new control. In [24], the problem of faulty plant reconfiguration has been solved based on the disturbance decoupling problem solution.

8. Conclusions

The paper deals with the problem of IWR to the unknown input in discrete-time nonlinear dynamic systems. So-called algebra of functions and logic–dynamic approach are used to solve the problem. The algebra of functions produces a solution in general form but demands analytical calculations. The advantage of the LDA is that only methods of linear algebra are used to solve the problems and the considered system may contain non-differential nonlinearities such as Coulomb friction, backlash and hysteresis. Moreover, the LDA methods can be applied both to the discrete-time and the continuous-time systems.

Funding

This research was funded by Russian Scientific Foundation (project 16-19-00046-P).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ding, S. Data-driven Design of Fault Diagnosis and Fault-tolerant Control Systems; Springer: London, UK, 2014. [Google Scholar]
  2. Blanke, M.; Kinnaert, M.; Lunze, J.; Staroswiecki, M. Diagnosis and Fault-Tolerant Control; Springer: Berlin, Germany, 2006. [Google Scholar]
  3. Alwi, H.; Edwards, C.; Tan, C. Fault Detection and Fault-Tolerant Control Using Sliding Modes; Springer: London, UK, 2011. [Google Scholar]
  4. Tabatabaeipour, S.; Stroustrup, M.; Bak, T. Fault-tolerant control of discrete-time LPV systems using virtual actuators. Int. J. Robust Nonlinear Control 2015, 25, 707–734. [Google Scholar]
  5. Puig, V. Fault diagnosis and fault tolerant control using set-membership approaches: Application to real case studies. Int. J. Appl. Math. Comput. Sci. 2010, 20, 619–635. [Google Scholar] [CrossRef] [Green Version]
  6. Andiarti, R.; Moog, C.H. Output feedback disturbance decoupling in nonlinear systems. IEEE Trans. Autom. Control 1996, 41, 1683–1689. [Google Scholar] [CrossRef]
  7. Battilotti, S. A sufficient condition for nonlinear disturbance decoupling with stability via measurement feedback. In Proceedings of the 36th Conference on Decision & Control, San Diago, CA, USA, 10–12 December 1997; pp. 3509–3514. [Google Scholar]
  8. Isidori, A.; Krener, A.J.; Gori-Giorgi, C.; Monaco, S. Nonlinear decoupling via feedback: A differential gemetric approach. IEEE Trans. Autom. Control 1981, 26, 331–345. [Google Scholar] [CrossRef]
  9. Xia, X.; Moog, C.H. Disturbance decoupling by measurement feedback for SISO nonlinear systems. IEEE Trans. Autom. Control 1999, 44, 1425–1429. [Google Scholar] [CrossRef]
  10. Wonham, W. Linear Multivariable Control: A Geometric Approach; Springer: Berlin, Germany, 1979. [Google Scholar]
  11. Isidori, A. Nonlinear Control Systems; Springer: London, UK, 1995. [Google Scholar]
  12. Grizzle, J. Controlled invariance for discrete-time nonlinear systems with an application to the disturbance decoupling problem. IEEE Trans. Autom. Control 1985, 30, 868–873. [Google Scholar] [CrossRef]
  13. Aranda-Bricaire, E.; Kotta, U. A geometric solution to the dynamic disturbance decoupling for discrete-time nonlinear systems. Kybernetika 2004, 40, 197–206. [Google Scholar]
  14. Kaldmae, A.; Kotta, U. Disturbance decoupling for discrete-time nonlinear systems by static measurement feedback. In Proceedings of the 18th International Conference on Process Control, Tatranska Lomnica, Slovakia, 14–17 June 2011; pp. 135–140. [Google Scholar]
  15. Albertini, F.; Sontag, E. Discrete-time transitivity and accessibility: Analytical systems. SIAM J. Control Optim. 1993, 31, 1599–1622. [Google Scholar] [CrossRef] [Green Version]
  16. Califano, C.; Monaco, S.; Normand-Cyrot, D. Nonlinear non-interactive control with stability in discrete-time framework. Int. J. Control 2002, 75, 11–22. [Google Scholar] [CrossRef]
  17. Conte, G.; Moog, C.H.; Perdon, A.M. Algebraic Methods for Nonlinear Control Systems. Theory and Applications; Springer: Berlin, Germany, 2007. [Google Scholar]
  18. Zhirabok, A.; Shumsky, A. The Algebraic Methods for Analysis of Nonlinear Dynamic Systems; Dalnauka: Vladivostok, Russia, 2008. (In Russian) [Google Scholar]
  19. Shumsky, A.; Zhirabok, A. Unified approach to the problem of full decoupling via output feedback. Eur. J. Control 2010, 16, 313–325. [Google Scholar] [CrossRef]
  20. Hartmanis, J.; Stearns, R. The Algebraic Structure Theory of Sequential Machines; Prentice-Hall: New York, NY, USA, 1966. [Google Scholar]
  21. Available online: http://webmathematica.cc.ioc.ee/webmathematica/NLControl/main/index.html (accessed on 1 July 2014).
  22. Zhirabok, A.; Shumsky, A.; Solyanik, S.; Suvorov, A. Fault detection in nonlinear systems via linear methods. Int. J. Appl. Math. Comput. Sci. 2017, 27, 261–272. [Google Scholar] [CrossRef] [Green Version]
  23. Kaldmae, A.; Kotta, U.; Shumsky, A.; Zhirabok, A. Measurement feedback disturbance decoupling in discrete-time nonlinear systems. Automatica 2013, 49, 2887–2891. [Google Scholar] [CrossRef]
  24. Kaldmae, A.; Kotta, U.; Jiang, B.; Shumsky, A.; Zhirabok, A. Faulty plant reconfiguration based on disturbance decoupling methods. Asian J. Control 2016, 8, 858–867. [Google Scholar] [CrossRef]
  25. Yan, Z.; Edwards, C. Nonlinear robust fault reconstruction and estimation using a sliding modes observer. Automatica 2007, 43, 1605–1614. [Google Scholar] [CrossRef]
  26. Zhirabok, A.; Zuev, A.; Shumsky, A. Diagnosis of linear dynamic systems: An approach based on sliding mode observers. Autom. Remote Control 2020, 81, 211–225. [Google Scholar] [CrossRef]
  27. Zhirabok, A. Disturbance decoupling problem: Logic-dynamic approach-based solution. Symmetry 2019, 11, 555. [Google Scholar] [CrossRef] [Green Version]
  28. Bobko, E.; Zhirabok, A.; Shumsky, A. Method of fault accommodation in technical systems. J. Comput. Syst. Sci. Int. 2016, 55, 735–749. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Zhirabok, A. The Problem of Invariance in Nonlinear Discrete-Time Dynamic Systems. Symmetry 2020, 12, 1241. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12081241

AMA Style

Zhirabok A. The Problem of Invariance in Nonlinear Discrete-Time Dynamic Systems. Symmetry. 2020; 12(8):1241. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12081241

Chicago/Turabian Style

Zhirabok, Alexey. 2020. "The Problem of Invariance in Nonlinear Discrete-Time Dynamic Systems" Symmetry 12, no. 8: 1241. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12081241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop