Next Article in Journal
Semi-Dilute Dumbbells: Solutions of the Fokker–Planck Equation
Previous Article in Journal
Filtering-Based Instance Selection Method for Overlapping Problem in Imbalanced Datasets
Previous Article in Special Issue
Fractional SIR-Model for Estimating Transmission Dynamics of COVID-19 in India
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Stability and Exponential Decay of Processes in Nonlinear Feedback Systems with Different Fractional Orders

Faculty of Electrical Engineering, Bialystok University of Technology, Wiejska 45A, 15-351 Białystok, Poland
*
Author to whom correspondence should be addressed.
Submission received: 10 June 2021 / Revised: 28 June 2021 / Accepted: 2 July 2021 / Published: 12 July 2021
(This article belongs to the Special Issue Global Stability Analysis Non-Linear Systems)

Abstract

:
The global stability of continuous-time multi-input multi-output nonlinear feedback systems with different fractional orders and interval matrices of positive linear parts is investigated. New sufficient conditions for the global stability of this class of positive nonlinear systems are established. Sufficient conditions for the exponential decay of processes in fractional nonlinear systems are given. Procedures for computation of a gain matrix characterizing the class of nonlinear elements are proposed and illustrated by examples.

1. Introduction

In positive systems inputs, state variables and outputs take only nonnegative values for any nonnegative inputs and nonnegative initial conditions [1,2,3]. Examples of positive systems are industrial processes involving chemical reactors, heat exchangers and distillation columns, storage systems, compartmental systems, water and atmospheric pollutions models. A variety of models having positive behavior can be found in engineering, management science, economics, social sciences, biology and medicine, etc. An overview of state of the art positive systems theory is given in the monographs [1,2,3,4,5].
Mathematical fundamentals of the fractional calculus are given in the monographs [4,5,6,7]. The positive fractional linear systems have been investigated in [4,5,8,9,10,11,12,13,14,15,16,17,18,19]. Positive linear systems with different fractional orders have been addressed in [11,12,19]. Descriptor positive systems have been analyzed in [13,20] and their stabilization in [18,19]. Linear positive electrical circuits with state feedbacks have been addressed in [20,21]. The superstabilization of positive linear electrical circuits by state feedbacks have been analyzed in [22] and the stability of nonlinear systems in [21,23]. The global stability of nonlinear systems with negative feedbacks and positive not necessary asymptotically stable linear parts has been investigated in [9,24,25]. The global stability of nonlinear standard and fractional positive feedback systems has been considered in [23].
In this paper the global stability of continuous-time multi-input multi-output nonlinear feedback systems with different fractional orders and interval matrices of positive linear parts will be addressed, and sufficient conditions for exponential decay of processes in fractional nonlinear systems with different orders will be proposed.
The paper is organized as follows. In Section 2 the basic definitions and theorems concerning the positive linear systems with different fractional orders are recalled. The stability of fractional interval positive linear systems is analyzed in Section 3. New sufficient conditions for the global stability of these feedback nonlinear systems with interval matrices of positive linear parts are established in Section 4. In Section 5, a procedure for calculation of a gain matrix characterizing the class of nonlinear elements is presented and illustrated by numerical examples. Sufficient conditions for the exponential decay of processes in fractional nonlinear systems with different orders are proposed in Section 6. Concluding remarks are given in Section 7.
The following notation will be used: —the set of real numbers, n × m —the set of n × m real matrices, + n × m —the set of n × m real matrices with nonnegative entries and + n = + n × 1 , M n —the set of n × n Metzler matrices (real matrices with nonnegative off-diagonal entries), I n — the n × n identity matrix, A ( : , n ) —the sum of all elements of nth column.

2. Positive Different Fractional Orders Linear Systems

Consider the fractional continuous-time linear system
d α x ( t ) d t α = A x ( t ) + B u ( t ) ,   0 < α < 1
y ( t ) = C x ( t ) + D u ( t ) ,
where x ( t ) n , u ( t ) m , y ( t ) p are the state, input and output vectors and A n × n , B n × m , C p × n , D p × m ,
d α x ( t ) d t α = 1 Γ ( 1 α ) 0 t x ˙ ( τ ) ( t τ ) α d τ ,   x ˙ ( τ ) = d x ( τ ) d τ
is the Caputo fractional derivative and
Γ ( z ) = 0 t z 1 e t d t ,   Re ( z ) > 0
is the gamma function [4,5].
Definition 1.
[4,5] The fractional system (1), (2) is called (internally) positive if x ( t ) + n and y ( t ) + p , t 0 for any initial conditions x ( 0 ) + n and all inputs u ( t ) + m , t 0 .
Theorem 1.
[4,5] The fractional system (1), (2) is positive if and only if
A M n ,   B + n × m ,   C + p × n ,   D + p × m .
Definition 2.
[4,5] The positive fractional system (1), (2) (for u ( t ) = 0 ) is called asymptotically stable (the matrix A is Hurwitz) if
lim t x ( t ) = 0   for   any       x ( 0 ) + n .
Theorem 2.
[4,5] The positive linear system (1), (2) is asymptotically stable if and only if one of the following equivalent conditions is satisfied:
 (1) 
All coefficient of the characteristic polynomial
p n ( s ) = det [ I n s A ] = s n + a n 1 s n 1 + + a 1 s + a 0
are positive, i.e.,  a i > 0 for  i = 0 , 1 , , n 1 .
 (2) 
There exists strictly positive vector  λ T = [ λ 1 λ n ] T , λ k > 0 , k = 1 , , n such that
A λ < 0   or   λ T A < 0 .
Theorem 3.
The positive system (1), (2) is asymptotically stable if the sum of entries of each column (row) of the matrix A is negative.
Proof. 
Using (8) we obtain
A λ = a 11 a 1 n a n 1 a n n λ 1 λ n = a 11 a 1 n λ 1 + + a n 1 a n n λ n < 0 0
and the sum of entries of each column of the matrix A is negative since λ k > 0 , k = 1 , , n . The proof for rows is similar. □
Now consider the fractional linear system with two different fractional orders
d α x 1 ( t ) d t α d β x 2 ( t ) d t β = A 11 A 12 A 21 A 22 x 1 ( t ) x 2 ( t ) + B 1 B 2 u ( t ) ,
y ( t ) = [ C 1 C 2 ] x 1 ( t ) x 2 ( t ) ,
where 0 < α , β < 1 , x 1 ( t ) n 1 and x 2 ( t ) n 2 are the state vectors, A i j n i × n j , B i n i × m , C i p × n i ; i, j = 1,2; u ( t ) m is the input vector and y ( t ) p is the output vector. Initial conditions for (10) have the form
x 1 ( 0 ) = x 10 ,   x 2 ( 0 ) = x 20   and   x 0 = x 10 x 20 .
Remark 1.
The state Equation (10) of fractional continuous-time linear systems with two different fractional orders has a similar structure to the 2D Roeesser type models.
Definition 3.
The fractional system (10), (11) is called positive if x 1 ( t ) + n 1 , x 2 ( t ) + n 2 and y ( t ) + p , t 0 for any initial conditions x 10 + n 1 , x 20 + n 2 and all input vectors u + m , t 0 .
Theorem 4.
[11,12] The fractional system (10), (11) for 0 < α < 1 ; 0 < β < 1 is positive if and only if
A ˜ = A 11 A 12 A 21 A 22 M N ,   B ˜ = B 1 B 2 + N × m ,   C ˜ = [ C 1 C 2 ] + p × n ( N = n 1 + n 2 ) .
Theorem 5.
[4] The positive fractional system (10), (11) is asymptotically stable if and only if one of the following equivalent conditions is satisfied:
 (1) 
All coefficients of the characteristic polynomial
det [ I n s A ˜ ] = s n + a ˜ n 1 s n 1 + + a ˜ 1 s + a ˜ 0
are positive, i.e., a ˜ i > 0 for i = 0 , 1 , , n 1
 (2) 
There exists a strictly positive vector λ = [ λ 1 λ n ] , λ k > 0 , k = 1 , , n such that
A ˜ λ < 0 or λ T A ˜ < 0 .
Theorem 6.
The positive system (10), (11) is asymptotically stable if the sum of entries of each column (row) of the matrix  A ˜ is negative.
Proof. 
Proof is similar to the proof of Theorem 3. □
Theorem 7.
The solution of the Equation (10) for 0 < α < 1 ; 0 < β < 1 with initial conditions (12) has the form
x ( t ) = x 1 ( t ) x 2 ( t ) = Φ 0 ( t ) x 0 + 0 t M ( t τ ) u ( τ ) d τ ,
where
M t = Φ 1 t B 10 + Φ 2 t B 01 = Φ 11 1 t Φ 12 1 t Φ 21 1 t Φ 22 1 t B 1 0 + Φ 11 2 t Φ 12 2 t Φ 21 2 t Φ 22 2 t 0 B 2 = Φ 11 1 t B 1 + Φ 12 2 t B 2 Φ 21 1 t B 1 + Φ 22 2 t B 2 = Φ 11 1 t Φ 12 2 t Φ 21 1 t Φ 22 2 t B 1 B 2
and
Φ 0 ( t ) = k = 0 l = 0 T k l t k α + l β Γ ( k α + l β + 1 ) ,
Φ 1 ( t ) = k = 0 l = 0 T k l t ( k + 1 ) α + l β 1 Γ ( k + 1 ) α + l β ,
Φ 2 ( t ) = k = 0 l = 0 T k l t k α + ( l + 1 ) β 1 Γ k α + ( l + 1 ) β ,
T k l = I n for k = l = 0 A 11 A 12 0 0 for k = 1 , l = 0 0 0 A 21 A 22 for k = 0 , l = 1 T 10 T k 1 , l + T 01 T k , l 1 for k + l > 1
Proof. 
Proof is given in [12]. □
Note that if α = β then from (18) we have
Φ 0 α = β ( t ) = k = 0 A k t k α Γ ( k α + 1 ) .

3. Stability of Fractional Interval Positive Linear Systems

Consider the fractional interval positive linear continuous-time system
d α x d t α = A x , 0 < α < 1 ,
where x = x ( t ) n is the state vector and the interval matrix A M n is defined by
A ¯ A A ¯   or   equivalently   A [ A ¯ , A ¯ ] .
Definition 4.
The fractional interval positive system (23) is called asymptotically stable if the system is asymptotically stable for all matrices  A M n satisfying the condition (24).
By condition (8) of Theorem 2, the positive system (23) is asymptotically stable if there exists a strictly positive vector λ > 0 such that the condition (8) is satisfied.
For two fractional positive linear systems
d α x d t α = A ¯ x ,   A ¯ M n
and
d α x d t α = A ¯ x ,   A ¯ M n
there exists a strictly positive vector λ + n such that
A ¯ λ < 0   and   A ¯ λ < 0
if and only if the systems (25), (26) are asymptotically stable.
Theorem 8.
If the matrices A ¯ and A ¯ of fractional positive systems (25), (26) are asymptotically stable then their convex linear combination
A = ( 1 k ) A ¯ + k A ¯   for   0 k 1
is also asymptotically stable.
Proof. 
By condition (8) of Theorem 2, if the fractional positive linear systems (25), (26) are asymptotically stable then there exists a strictly positive vector λ + n such that (27) holds. Using (8) and (27) we obtain
A λ = [ ( 1 k ) A ¯ + k A ¯ ] λ = ( 1 k ) A ¯ λ + k A ¯ λ < 0
for 0 k 1 . Therefore, if the positive linear systems (25), (26) are asymptotically stable then their convex linear combination (28) is also asymptotically stable. □
Theorem 8.
The interval positive systems (23) are asymptotically stable if and only if the positive linear systems (25), (26) are asymptotically stable.
Proof. 
By condition (8) of Theorem 2 if the matrices A ¯ M n , A ¯ M n are asymptotically stable, then there exists a strictly positive vector λ + n such that (8) holds. The convex linear combination (28) satisfies the condition A λ < 0 if and only if (29) holds. Therefore, the interval system (23) is asymptotically stable if and only if the positive linear system is asymptotically stable. □
Example 1.
Consider the fractional interval positive linear continuous-time system (23) with the matrices
A ¯ = 3 2 2 4 , A ¯ = 2 1 1 3 .
Using the condition (8) of Theorem 2, we choose λ = [ 1 1 ] T and we obtain
A ¯ λ = 3 2 2 4 1 1 = 1 2 < 0
and
A ¯ λ = 2 1 1 3 1 1 = 1 2 < 0 .
Therefore, the matrices (30) are Hurwitz.
These considerations can be easily extended to positive different fractional orders linear systems (10), (11).

4. Global Stability of Fractional Nonlinear Positive Feedback Systems

Consider the m-input p-output (MIMO) nonlinear feedback system shown in Figure 1 which consists of the positive fractional linear part, the nonlinear element with the matrix characteristic u = f ( e ) and the feedback with positive gain matrix H. The positive fractional linear part is described by the equations (25), (26) with the interval matrices
A ˜ [ A ¯ , A ¯ ] M n , B ˜ [ B ¯ , B ¯ ] + n × m , C ˜ [ C ¯ , C ¯ ] + p × n .
It is assumed that the interval matrix A ˜ M n is Hurwitz.
The characteristic f ( e ) of the nonlinear element satisfies the condition
0 f ( e ) K e   or   0 u K e ,
where
u = u 1 u m , K = k 11 k 1 p k m 1 k m p , e = e 1 e p
and
f ( 0 ) = 0 ,   u i = f ( e 1 , , e p ) k i 1 e 1 + + k i p e p , i = 1 , , m .
Remark 2.
The matrix K forms m, p-dimensional convex surfaces which constricts possible dynamics of nonlinear elements, i.e., for m = p = 1 the nonlinear characteristicf(e) is presented in Figure 2.
Definition 5.
The fractional nonlinear positive system is called globally stable if it is asymptotically stable for all nonnegative initial conditions  x ( 0 ) + n .
The following theorem gives sufficient conditions for the global stability of the positive nonlinear system.
Theorem 9.
The fractional nonlinear system consisting of the positive asymptotically stable linear part described by (25), (26) with interval matrices (33), the nonlinear element satisfying the condition (34) and the feedback with positive gain matrix  H + m × p is globally stable if the sum of entries of each column (row) of the matrix
( 1 q ) A ¯ + q A ¯ + B ˜ K H C ˜ = A ¯ + B ¯ K ¯ H C ¯ M n for q = 0 A ¯ + B ¯ K ¯ H C ¯ M n for q = 1
is negative.
Proof. 
The proof will be accomplished by the use of the Lyapunov method [26,27]. As the Lyapunov function V ( x ) we choose
V 1 ( x 1 ) + V 2 ( x 2 ) = λ 1 T x 1 + λ 2 T x 2 0   for   x = x 1 x 2 + n , λ = λ 1 λ 2 + n ,
where λ is a strictly positive vector, i.e., λ i j > 0 , i = 1 , 2 ; j = 1 , , n .
Using (37) and (25) we obtain
d α V 1 ( x 1 ) d t α + d β V 2 ( x 2 ) d t β = [ λ 1 T λ 2 T ] d α x 1 d t α d β x 2 d t β = λ T ( A ˜ x + B ˜ u ) = λ T ( A ˜ x + B ˜ f ( e ) ) λ T ( A ˜ + B ˜ K H C ˜ ) x
since u = f ( e ) K e = K H C ˜ x .
From (39) it follows that d α V 1 ( x 1 ) d t α + d β V 2 ( x 2 ) d t β < 0 if the sum of entries of each column (row) of the matrix (37) is negative (Theorem 3) and the nonlinear positive system is globally stable. □

5. Procedure and Example

To find the K + m × p satisfying the condition (37) for the nonlinear positive system, the following procedure can be used.
Procedure 1.
Step 1.
Using the matrices A ¯ ,   B ¯ ,   C ¯ of the positive linear system and the matrix H compute the maximum value of the matrix K ¯ such that the sum of all entries of each column (row) of the matrix
A ¯ ¯ = A ¯ + B ¯ K ¯ H C ¯
is negative. Entries of the matrix K can be computed as the solution of the linear matrix equation
G k = h ,
where the matrix G and the column vector h are defined by the sum of entries of each column (row) of the matrix (40) and vector k contains components of matrix K. If m   p > n or/and rank   G n , then we choose arbitrarily m   p rank   G nonnegative entries of the matrix K.
Step 2.
Using the matrices A ¯ ,   B ¯ ,   C ¯ of the positive linear system and the matrix H compute the maximum value of the matrix K ¯ such that the sum of all entries of each column (row) of the matrix
A ¯ ¯ = A ¯ + B ¯ K ¯ H C ¯
is negative.
Step 3.
Taking into consideration K ¯ computed in Step 1 and K ¯ from Step 2, find the desired K + m × p for which the matrices A ¯ ¯ and A ¯ ¯ are Hurwitz, i.e., the characteristic f(e) satisfy the condition (34).
Remark 3.
The conditions of Theorem 2 can be also used to compute the entries of the matrix K . Usually in this case the computations are more complicated.
Example 2.
Consider the nonlinear feedback system shown in Figure 1 with the interval matrices of the positive linear part
A ¯ = 15.8 0.2 0.9 0.4 14.1 0.5 0.5 0.4 20.1 , A ¯ = 7.9 0.1 0.7 0.2 7.9 0.3 0.1 0.2 12.1 , B ¯ = 0.5 0.4 0.3 0.5 0.8 0.6 , B ¯ = 0.4 0.3 0.2 0.4 0.6 0.5 , C ¯ = 0.4 0.6 0.8 0.6 0.4 0.6 , C ¯ = 0.3 0.5 0.6 0.4 0.3 0.5
and the gain matrix
H = 2 1 1 5 .
From (43) it follows that m = p = 2 and u 1 = f ( e 1 , e 2 ) k 11 e 1 + k 12 e 2 , u 2 = f ( e 1 , e 2 ) k 21 e 1 + k 22 e 2 , then we are looking for
K = k 11 k 12 k 21 k 22 + 2 × 2 ,
for which the nonlinear feedback system is globally stable.
Using Procedure 1 we obtain:
Step 1.
Using (40), (43) and (44) we determined that the sum of columns are the following:
A ¯ ¯ ( : , 1 ) = 2.24 k 11 + 5.44 k 12 + 2.1 k 21 + 5.1 k 22 14.9 , A ¯ ¯ ( : , 2 ) = 2.56 k 11 + 4.16 k 12 + 2.4 k 21 + 3.9 k 22 13.5 , A ¯ ¯ ( : , 3 ) = 3.52 k 11 + 6.08 k 12 + 3.3 k 21 + 5.7 k 22 18.7
and taking into consideration (5.2), we have
G = 2.24 5.44 2.1 5.1 2.56 4.15 2.4 3.9 3.52 6.08 3.3 5.7 , k = k 11 k 12 k 21 k 22 , h = 14.9 13.5 18.7 .
Since rank   G = 2 , then we choose m p rank   G = 2 elements of the vector k (entries of the matrix K) as k 21 = 1 , k 22 = 1 . The solution for (41) with (47) is the following
k 11 k 12 k 21 k 22 = 1.309 0.869 1 1
and for
K ¯ < = 1.309 0.869 1 1
the system with A ¯ ,   B ¯ ,   C ¯ is stable.
Step 2.
Similar to Step 1, using (42), (43) and (44) we can write the linear matrix Equation (41) in the form
1.1 2.53 1.2 2.76 1.43 2.2 1.56 2.4 1.87 3.41 2.04 3.72 k 11 k 12 k 21 k 22 = 7.6 7.6 11.1 .
Since rank   1.1 2.53 1.2 2.76 1.43 2.2 1.56 2.4 1.87 3.41 2.04 3.72 = 2 , then we choose two elements of the vector k as k 11 = 1 , k 12 = 1 . The solution for (50) is
k 11 k 12 k 21 k 22 = 1 1 1.013 1.002 ,
and
K ¯ < = 1 1 1.013 1.002
Therefore, the system with A ¯ ,   B ¯ ,   C ¯ is stable.
Step 3.
Taking into consideration (49) and (52), we have that for
K = 1 1 1 1
the matrices A ¯ ¯ and A ¯ ¯ are Hurwitz since
A ¯ ¯ ( : , 1 ) = 0.01 , A ¯ ¯ ( : , 2 ) = 0.01 , A ¯ ¯ ( : , 3 ) = 0.06 , A ¯ ¯ ( : , 1 ) = 0.02 , A ¯ ¯ ( : , 2 ) = 0.48 , A ¯ ¯ ( : , 3 ) = 0.1 .
Therefore, for nonlinear elements satisfying the condition
0 0 f ( e 1 , e 2 ) 1 1 1 1 e 1 e 2 ,
the nonlinear feedback system with interval matrices (43) of positive linear parts and for the gain matrix (53), the nonlinear system is globally stable.

6. Exponential Decay of Processes in Nonlinear Feedbacks Systems

Consider the nonlinear feedback system shown in Figure 3 which consists of the positive linear part, the nonlinear element with characteristic u = f ( e ) and positive dynamical feedback. The linear part is described by the equations
x ˙ = A x + B u , y = C x ,
with interval matrices
A ¯ A A ¯ , B ¯ B B ¯ , C ¯ C C ¯ ,
where x = x ( t ) + n 1 ,   u = u ( t ) + ,   y = y ( t ) + is the state, input and output vectors of the system (56) and A M n 1 ,   B + n 1 × 1 ,   C + 1 × n 1 . It is assumed that the matrix A of (56) has all eigenvalues s k with real parts smaller than
Re s k < γ ,   i . e . ,   γ > 0 , k = 1 , , n .
The characteristic of the nonlinear element is shown in Figure 2 and it satisfies the condition
0 f ( e ) e k < .
The positive feedback system is described by the equations
z ˙ = F z + G y , e = H z ,
with interval matrices
F ¯ F F ¯ , G ¯ G G ¯ , H ¯ H H ¯ ,
where z = z ( t ) + n 2 ,   e = e ( t ) + are the state vector and output vectors. It is assumed that the matrix F + I n γ M n 2 is also asymptotically stable. From (56) and (60) we have
x ˙ z ˙ = A ^ x z + B ^ u ,
where
A ^ = A 0 G C F M n , B ^ = B 0 + n × 1 , n = n 1 + n 2 .
The following theorem gives sufficient conditions for the exponential decay of transient values in the positive feedback nonlinear system faster than e γ   t .
Theorem 10.
The state variables of the nonlinear system consisting of the positive linear part (56), the nonlinear element satisfying the condition (34), (35), (36) and positive asymptotically stable dynamical feedback system (60) are decaying exponentially faster than  e γ   t if the matrix
A ¯ + I n 1 γ k B ¯ H ¯ G ¯ C ¯ F ¯ + I n 2 γ M n
is asymptotically stable.
Proof. 
The proof will be accomplished by the use of the Lyapunov method [26,27]. As the Lyapunov function V ( x , z ) we choose
V ( x , z ) = λ T x z 0   for   x + n 1 , z + n 2 ,
where λ + n is a strictly positive vector, i.e., λ k > 0 , k = 1 , , n .
It is well-known that if the matrix A M n is asymptotically stable, then the state variables of the system x ˙ = ( A + I n γ ) x are decaying exponentially faster than e γ   t .
Using (56), (60) and (65) we obtain
V ˙ ( x , z ) = λ T x ˙ z ˙ = λ T A 0 G C F x z + B 0 u = λ T A 0 G C F x z + B 0 k H z = λ T A k B H G C F x z 0
since B u = B f ( e ) k B H z by the condition (6.2).
From (66) it follows that V ˙ ( x , z ) 0 if the matrix (64) is asymptotically stable, and therefore the state variables decay exponentially faster than e γ   t . □
Theorem 6.1 can be applied to solve the following two problems.
Problem 1.
Given matrices A, B, C and F, G, H of the positive systems (56), (60) and the nonlinear characteristic u = f ( e ) of the nonlinear element. Knowing the value of k satisfying the condition (59), check if the transient processes in the nonlinear system decay faster than e γ   t .
Problem 2.
Given matrices A, B, C and F, G, H of the positive systems (56), (60) and the nonlinear characteristic u = f ( e ) of the nonlinear element. Find the maximal value of k for which the characteristic u = f ( e ) of the nonlinear element satisfies the condition (59) and the transient values of the nonlinear system decay faster than e γ   t .
The Problem 1 can be solved by the use of the following:
Procedure 2.
Step 1.
Knowing the characteristic u = f ( e ) , find the minimal value of k satisfying the condition (59).
Step 2.
Using Theorem 10, find the sum of entries of each column (row) of the matrix (64). If all these sums are negative, then the transient processes in the nonlinear system decay faster than e γ   t .
The Problem 2 can be solved by the use of the following:
Procedure 2.
Step 1.
Using Theorem 10, find the sum of entries of each column (row) of the matrix (64).
Step 2.
Find the maximal value of k c ( k r ) for which the sums of entries of all columns (rows) of (64) are negative.
Step 3.
Find k max = min ( k c , k r ) .
In this case, the transient process in the nonlinear system decrease faster than e γ   t for all nonlinear characteristics u = f ( e ) satisfying the condition
0 < f ( e ) < k max e .
Remark 4.
The value of k max depends only on the first n 1 rows and of the last n 2 columns of the matrix (64).
Example 3.
Consider the nonlinear system shown in Figure 3 with linear positive parts described by (56) and (60) with
A ¯ = 6.5 1 2 6.2 , A ¯ = 7 1.5 2.3 6.5 , B ¯ = 0.5 0.7 , B ¯ = 0.5 0.7 , C ¯ = [ 0.4 0.5 ] , C ¯ = [ 0.5 0.6 ]
and
F ¯ = 6 2 1.6 7 , F ¯ = 6.2 2.2 1.8 7.3 , G ¯ = 0.8 0.6 , G ¯ = 1 0.8 , ( 6.10 ) H ¯ = [ 0.4 0.2 ] , H ¯ = [ 0.5 0.4 ]
respectively, and the nonlinear element with characteristics satisfying the condition (59).
Case 1.
Using k = 1 , check the global stability of the nonlinear system for γ = 2 . In this case, using (6.5), (6.9) and (6.10) we obtain
A ¯ + I n 1 γ k B ¯ H ¯ G ¯ C ¯ F ¯ + I n 2 γ = 4.5 1 0.3 0.24 2 4.2 0.4 0.32 0.5 0.6 4 2 0.4 0.48 1.6 5
The sums of the entries of columns of the matrix (6.11) are: column 1: = 1.6 , column 2: = 2.12 , column 3: = 1.7 , column 4: = 2.44 . Therefore, by Theorem 6.1, the nonlinear system is globally stable.
Case 2.
Find the maximal value of k max satisfying the condition (4.2) for which the transient process in the nonlinear system decreases faster than e γ   t . Using Procedure 2 we obtain the following:
Step 1.
The sums of entries of each column (row) of the matrix
A ¯ + I n 1 γ k B ¯ H ¯ G ¯ C ¯ F ¯ + I n 2 γ = 4.5 1 0.3 k 0.24 k 2 4.2 0.4 k 0.32 k 0.5 0.6 4 2 0.4 0.48 1.6 5  
are: column 1: = 1.6 , column 2: = 2.12 , column 3: = 0.7 k 2.4 , column 4: = 0.56 k 3 , row 1: = 2.96 , row 2: = 1.48 , row 3: = 0.54 k 3.5 , row 4: = 0.72 k 2.2 .
Step 2.
From Theorem 10 we have: for column 3: k < 3.428 and for column 4: k < 5.357 and for row 1: k < 6.482, row 2: k < 3.0555.
Step 3.
The desired value of k is k max = min ( k c , k r ) = 3.0555 .
Therefore, the transient process in the nonlinear system with characteristics satisfying the condition (34), (35), (36) for k < 3.0555 decreases faster than e γ   t .
Remark 3.
From matrix (6.5) and the computation procedure, it follows that the k depends only on the matrices F, G, H and is independent of the matrices A, C, G.

7. Concluding Remarks

The global stability of continuous-time nonlinear feedback systems with different fractional orders and interval matrices of positive linear parts has been investigated. New sufficient conditions for the global stability of this class of positive nonlinear systems are established (Theorem 9). The procedure for the calculation of a gain matrix characterizing the class of nonlinear elements is presented and illustrated by numerical example. Sufficient conditions for the exponential decay of processes in nonlinear systems have been proposed (Theorem 10) and illustrated by numerical example. The considerations can be extended to the discrete-time nonlinear systems with different fractional orders and interval matrices of positive linear parts. Further investigation could address the extension of the considerations to nonlinear different fractional order systems with e time-varying linear parts.

Author Contributions

Conceptualization, T.K. and Ł.S.; methodology, T.K.; formal analysis, T.K.; writing—original draft preparation, Ł.S.; writing—review and editing, Ł.S.; supervision, T.K.; funding acquisition, T.K. and Ł.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Science Centre in Poland under work No. 2017/27/B/ST7/02443.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

The study did not involve humans.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berman, A.; Plemmons, R.J. Nonnegative Matrices in the Mathematical Sciences; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
  2. Farina, L.; Rinaldi, S. Positive Linear Systems: Theory and Applications; J. Wiley: New York, NY, USA, 2000. [Google Scholar]
  3. Kaczorek, T. Positive 1D and 2D Systems; Springer: London, UK, 2002. [Google Scholar]
  4. Kaczorek, T. Selected Problems of Fractional Systems Theory; Springer: Berlin, Germany, 2011. [Google Scholar]
  5. Kaczorek, T.; Rogowski, K. Fractional Linear Systems and Electrical Circuits; Springer: Cham, Switzerland, 2015. [Google Scholar]
  6. Ostalczyk, P. Discrete Fractional Calculus; World Scientific: River Edgle, NJ, USA, 2016. [Google Scholar]
  7. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  8. Busłowicz, M.; Kaczorek, T. Simple conditions for practical stability of positive fractional discrete-time linear systems. Int. J. Appl. Math. Comput. Sci. 2009, 19, 263–269. [Google Scholar] [CrossRef] [Green Version]
  9. Kaczorek, T. Absolute stability of a class of fractional positive nonlinear systems. Int. J. Appl. Math. Comput. Sci. 2019, 29, 93–98. [Google Scholar] [CrossRef] [Green Version]
  10. Kaczorek, T. Analysis of positivity and stability of fractional discrete-time nonlinear systems. Bull. Pol. Acad. Sci. Techn. 2016, 64, 491–494. [Google Scholar] [CrossRef] [Green Version]
  11. Kaczorek, T. Positive linear systems with different fractional orders. Bull. Pol. Acad. Sci. Techn. 2010, 58, 453–458. [Google Scholar] [CrossRef]
  12. Kaczorek, T. Positive linear systems consisting of n subsystems with different fractional orders. IEEE Trans. Circuits Syst. 2011, 58, 1203–1210. [Google Scholar] [CrossRef]
  13. Kaczorek, T. Positive fractional continuous-time linear systems with singular pencils. Bull. Pol. Acad. Sci. Techn. 2012, 60, 9–12. [Google Scholar] [CrossRef] [Green Version]
  14. Kaczorek, T. Stability of fractional positive nonlinear systems. Arch. Control. Sci. 2015, 25, 491–496. [Google Scholar] [CrossRef]
  15. Mitkowski, W. Dynamical properties of Metzler systems. Bull. Pol. Acad. Sci. Techn. 2008, 56, 309–312. [Google Scholar]
  16. Ruszewski, A. Stability of discrete-time fractional linear systems with delays. Arch. Control Sci. 2019, 29, 549–567. [Google Scholar]
  17. Ruszewski, A. Practical and asymptotic stabilities for a class of delayed fractional discrete-time linear systems. Bull. Pol. Acad. Sci. Techn. 2019, 67, 509–515. [Google Scholar]
  18. Sajewski, Ł. Decentralized stabilization of descriptor fractional positive continuous-time linear systems with delays. In Proceedings of the 22nd International Conference Methods and Models in Automation and Robotics, Międzyzdroje, Poland, 28–31 August 2017; pp. 482–487. [Google Scholar]
  19. Sajewski, Ł. Stabilization of positive descriptor fractional discrete-time linear systems with two different fractional orders by decentralized controller. Bull. Pol. Acad. Sci. Techn. 2017, 65, 709–714. [Google Scholar] [CrossRef] [Green Version]
  20. Borawski, K. Modification of the stability and positivity of standard and descriptor linear electrical circuits by state feedbacks. Electr. Rev. 2017, 93, 176–180. [Google Scholar] [CrossRef] [Green Version]
  21. Kaczorek, T.; Borawski, K. Stability of Positive Nonlinear Systems. In Proceedings of the 22nd International Conference Methods and Models in Automation and Robotics, Międzyzdroje, Poland, 28–31 August 2017. [Google Scholar]
  22. Kaczorek, T. Superstabilization of positive linear electrical circuit by state-feedbacks. Bull. Pol. Acad. Sci. Techn. 2017, 65, 703–708. [Google Scholar] [CrossRef] [Green Version]
  23. Kaczorek, T. Global stability of positive standard and fractional nonlinear feedback systems. Bull. Pol. Acad. Sci. Techn. 2020, in press. [Google Scholar]
  24. Kaczorek, T. Analysis of positivity and stability of discrete-time and continuous-time nonlinear systems. Comput. Probl. Electr. Eng. 2015, 5, 11–16. [Google Scholar]
  25. Kaczorek, T. Global stability of nonlinear feedback systems with positive linear parts. Int. J. Nonlinear Sci. Numer. Simul. 2019, 20, 575–579. [Google Scholar] [CrossRef]
  26. Lyapunov, A.M. Obscaja Zadaca ob Ustoicivosti Dvizenija; Gostechizdat: Moskwa, Russia, 1963. [Google Scholar]
  27. Leipholz, H. Stability Theory; Academic Press: New York, NY, USA, 1970. [Google Scholar]
Figure 1. The nonlinear feedback system.
Figure 1. The nonlinear feedback system.
J 04 00025 g001
Figure 2. Characteristic of single-input single-output nonlinear element.
Figure 2. Characteristic of single-input single-output nonlinear element.
J 04 00025 g002
Figure 3. The nonlinear feedback system.
Figure 3. The nonlinear feedback system.
J 04 00025 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kaczorek, T.; Sajewski, Ł. Global Stability and Exponential Decay of Processes in Nonlinear Feedback Systems with Different Fractional Orders. J 2021, 4, 328-340. https://0-doi-org.brum.beds.ac.uk/10.3390/j4030025

AMA Style

Kaczorek T, Sajewski Ł. Global Stability and Exponential Decay of Processes in Nonlinear Feedback Systems with Different Fractional Orders. J. 2021; 4(3):328-340. https://0-doi-org.brum.beds.ac.uk/10.3390/j4030025

Chicago/Turabian Style

Kaczorek, Tadeusz, and Łukasz Sajewski. 2021. "Global Stability and Exponential Decay of Processes in Nonlinear Feedback Systems with Different Fractional Orders" J 4, no. 3: 328-340. https://0-doi-org.brum.beds.ac.uk/10.3390/j4030025

Article Metrics

Back to TopTop