Next Article in Journal
SEIR Mathematical Model of Convalescent Plasma Transfusion to Reduce COVID-19 Disease Transmission
Next Article in Special Issue
Fractional Systems’ Identification Based on Implicit Modulating Functions
Previous Article in Journal
Comparison of Molecular Geometry Optimization Methods Based on Molecular Descriptors
Previous Article in Special Issue
On the Initial Value Problems for Caputo-Type Generalized Proportional Vector-Order Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives

1
College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
2
Faculty of Economics and Business Administration, Sofia University “St. Kl. Ohridski”, 125 Tzarigradsko Chaussee Blvd., Bl. 3, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(22), 2856; https://doi.org/10.3390/math9222856
Submission received: 12 October 2021 / Revised: 6 November 2021 / Accepted: 9 November 2021 / Published: 10 November 2021
(This article belongs to the Special Issue Fractional Differential Equations and Control Problems)

Abstract

:
Positive linear electrical circuits systems described by generalized fractional derivatives are studied in this paper. We mainly focus on the reachability and observability of linear electrical circuits systems. Firstly, generalized fractional derivatives and ρ -Laplace transform of f is presented and some preliminary results are provided. Secondly, the positivity of linear electrical circuits systems described by generalized fractional derivatives is investigated and conditions for checking positivity of the systems are derived. Thirdly, reachability and observability of the generalized fractional derivatives systems are studied, in which the ρ -Laplace transform of a Mittag-Leffler function plays an important role. At the end of the paper, illustrative electrical circuits systems are presented, and conclusions of the paper are presented.

1. Introduction

Fractional differential equations play an important role in the analysis and modeling of various processes. Many classical methods, such as perturbation method [1], Fourier transform and Merlin transform, have important applications in fractional calculus. Fractional calculus has many new diffusion processes in physics [2,3,4]. There are many types of fractional derivatives in fractional calculus, see [5,6,7,8,9,10,11] for details. For the basic principles of fractional calculus and its most interesting applications, see [12,13,14,15,16,17,18]. In recent years, fractional calculus has become a useful and promising tool in modeling and analyzing different dynamic behavior processes. Fractional calculus system has attracted more and more scholars’ attention since its wide application value in the field of science and engineering.
Fractional derivative has more advantages than classical derivative. The first advantage is that the fractional derivative takes into account memory. Memory effect is a basic property of differential equations. This explains the application of fractional derivative in differential equation modeling [14,16]. Another advantage is that the fractional derivative produces many diffusion processes [4]. Due to these advantages, fractional derivative also has many applications in electronics [15,19,20].
The new mathematical model appears with the emergence of fractional derivative. In recent years, researchers began to apply fractional derivatives to describe electrical circuit systems. Many types of fractional electrical circuits have been introduced recently in the literature, see [14,15,19,20] and fractional RL (Resistor-Inductor) and RC (Resistor-Capacitance) circuit modeling in [11,18,21,22,23]. Many studies on the properties of fractional electrical circuit system will use numerical solutions and analytical solutions., because it is an important way to study properties of circuits systems. Reference [18] studies the analytical and numerical solutions of fractional RL and RC circuits using Atangana-Baleanu derivative and bi-order derivative. Reference [21] studies the solutions of fractional electrical RL, LC (Inductor-Capacitance) and RC electrical circuit systems described by Mittag-Leffler fractional derivatives. Fractional RL and LC electrical circuits systems are studied in reference [22]. They also compared the fractional electrical circuit systems with the traditional electrical circuit systems.
In the literature of fractional electrical circuits systems, many studies mainly focus on performance properties of systems, such as stability [11,23,24,25,26,27,28,29,30,31,32]. In [27,28], the asymptotic stability of integer and fractional positive continuous-time linear system with delays is studied, respectively. The necessary and sufficient conditions for the asymptotic stability of integer and fractional positive linear system with delays are given, and it is proved that the asymptotic stability integer and fractional of the systems is independent of delays. By comparison, it can be found that the asymptotic stability conditions and checking methods of integer and fractional positive linear systems with delays are similar. In recent years, generalized fractional derivative has been widely used in stability analysis. For more details, see [6,33,34]. In [29], several stabilities of fractional differential equations described by Caputo type generalized fractional derivatives are studied. In this paper, a new concept of stability, fractional input stability is introduced. Thus, it provides a new idea and method for the stability analysis of fractional differential equations described by generalized fractional derivatives. In [32], the stability of RLC (Resistor-Inductor-Capacitance) electrical circuits described by Caputo-Liouville generalized fractional derivatives is studied. The local asymptotic stability and global asymptotic stability of trivial equilibrium are analyzed. The results of this paper are very valuable for the study of circuit stability. Of course, in addition to stability, electrical circuits systems have other performance properties [35,36,37], such as reachability and observability. In [37], fractional positive discrete-time linear systems are studied in the literature. The necessary and sufficient conditions of positivity, reachability and controllability are given. In [38], the reachability and observability of integer and fractional positive linear electrical circuits are studied.And it is found that the reachability and observability of integer and fractional positive linear electrical circuits are invariant. At this time, reference [39] further extends the result of [38] to positive linear electrical circuits with delays. In [39], the electrical circuit systems are extended to positive linear electrical circuits with delays, and it is found that the reachability and observability of integer and fractional positive linear electrical circuits with delays are similar. Many scholars have studied the properties of electrical circuits described by generalized fractional derivatives and positive linear electrical circuit systems. In this paper, generalized fractional derivatives are applied to positive linear electrical circuit systems.
The generalized fractional derivative units the Riemann-Liouville fractional derivative and Hadamard fractional derivative into a unified form in that it is mediated by an additional fractional parameter, which is more general than the ordinary classical fractional derivative. By observing, we can find that when α is fixed, the smaller the parameter ρ is, the greater the initial slope of the described circuit trajectory is, and the more realistic the described trajectory is. Contributions of this paper are listed below:
(1) We present the positivity of linear electrical circuits systems described by generalized fractional derivatives and obtained conditions for checking positivity of the systems. Also we investigate the checking method of reachability and observability of the generalized fractional derivatives studied.
(2) In references [23,29], the ρ -Laplace transform was performed on the Mittag-Leffler function with a constant λ , but we are extended to the Mittag-Leffler function with matrix A, and the corresponding form of ρ -Laplace transform is obtained.
(3) We investigate the effect of the parameter ρ on the electrical circuits systems. Since the generalized fractional derivative involves parameters, which is more general than the classical fractional derivative. Therefore, we studied the electrical circuits systems by illustrative examples when the fixed order α is unchanged, We let the parameter ρ take different values to observe how the parameter ρ affects the change of state trajectories.
The remainder of the paper is organized as follows. In Section 2, generalized fractional integrals, derivatives and ρ -Laplace transform are reviewed and some new results are given. Positive linear electrical circuits systems described by generalized fractional derivatives are presented in Section 3. In Section 4, reachability of positive linear electrical circuits systems described by generalized fractional derivatives are investigated. Observability of positive linear electrical circuits systems described by generalized fractional derivatives are investigated in Section 5. The illustrative electrical circuits are presented in Section 6. Concluding remarks are given in Section 7.
The following notation will be used in this paper. R is a set of real numbers, R n × m is a set of n × m dimensional real matrix, R + n × m is a set of n × m dimensional real matrix with nonnegative entries and R + n = R + n × 1 , M n is a set of n × n Metzler matrix (real matrix with nonnegative off-diagonal entries).

2. Generalized Fractional Integrals, Derivatives and ρ -Laplace Transform

For the knowledge of classical fractional calculus, we refer to references [7,13], in which there are many types of fractional integrals and derivatives. In this paper, we consider the following forms of generalized fractional integrals:
( a I α , ρ f ) ( x ) = 1 Γ ( α ) a x ( x ρ t ρ ρ ) α 1 f ( t ) d t t 1 ρ
and
( I b α , ρ f ) ( x ) = 1 Γ ( α ) x b ( t ρ x ρ ρ ) α 1 f ( t ) d t t 1 ρ .
where α is the order, ρ is an additional fractional parameter, and integrals (1) and (2) are left generalized fractional integral and right generalized fractional integral, respectively. We can find that when ρ = 1 , the integrals (1) and (2) become Riemann-Liouville fractional integrals. See [13] for details. When the limit is ρ 0 , the integrals (1) and (2) will become Hadamard fractional integrals defined in [7].
The generalized fractional derivatives defined by [6] for order α > 0 are as follows:
( a D α , ρ f ) ( x ) = γ n ( a I n α , ρ f ) ( x ) = γ n Γ ( n α ) a x ( x ρ t ρ ρ ) n α 1 f ( t ) d t t 1 ρ
and
( D b α , ρ f ) ( x ) = ( γ n ) ( a I n α , ρ f ) ( x ) = γ n Γ ( α ) x b ( t ρ x ρ ρ ) n α 1 f ( t ) d t t 1 ρ
where fractional parameter ρ > 0 , n = α + 1 and γ = x 1 ρ d d x .
Derivatives (3) and (4) are left and right generalized fractional derivatives, respectively. We can find that when ρ = 1 , the derivatives (3) and (4) become Riemann-Liouville fractional derivatives. See [13] for details. When the limit is ρ 0 , the derivatives (3) and (4) will become Hadamard fractional derivatives defined in [7].
The reason we study the generalized fractional derivative is that it generalizes the Riemann-Liouville fractional derivative and Hadamard fractional derivative into one form. When the parameters are fixed at different values or take limits, the above derivative is generated as a special case [6]. The generalized fractional derivative is more general than the ordinary classical fractional derivative.
Let’s recall some details of the ρ -Laplace transform and its basic concepts, which plays an important role in the next research.
Definition 1
([9]). The ρ-Laplace transform of a continuous function f : [ 0 , ) R is defined by
L ρ { f ( x ) } ( s ) = 0 e s x ρ ρ f ( x ) d x x 1 ρ , ρ > 0 ,
The integral is valid for all values of s.
Theorem 1
([7]). If the ρ-Laplace transform of a continuous function f : [ 0 , ) R exists, the following relationship holds
L ρ { f ( x ) } ( s ) = L { f ( ( ρ x ) 1 ρ ) } ( s ) ,
where L { f } is the usual Laplace transform of f.
Theorem 2
([7]). If continuous function f : [ 0 , ) R has ρ-Laplace transform for s > c 1 and continuous function g : [ 0 , ) R has ρ-Laplace transform for s > c 2 . Then, continuous function a f + b g has Laplace transform for any constants a and b, and has the following linear property
L ρ { a f ( x ) + b g ( x ) } ( s ) = a L ρ { f ( x ) } ( s ) + b L ρ { g ( x ) } ( s ) , s > m a x { c 1 , c 2 } .
The reference [9] gave some ρ -Laplace transforms of elementary functions as shown below.
Theorem 3.
Suppose L ρ { · } is the ρ-Laplace transform
( 1 ) L ρ { 1 } ( s ) = 1 s , s > 0 . ( 2 ) L ρ { t p } ( s ) = ρ p ρ Γ ( 1 + p ρ ) s 1 + p ρ , p R , s > 0 . ( 3 ) L ρ { e λ t ρ ρ } ( s ) = 1 s λ , s > λ . ( 4 ) L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) } ( s ) = ( s α I A ) 1 , A R n × n , ρ ( A ) < s α and lim n ( A s α ) n = 0 ,
where ρ(A) is the spectral radius of the matrix A.
Proof. 
(1), (2) and (3) have been proved, please refer to [9] for details. Only (4) will be proved below
L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) } ( s ) = k = 0 A k Γ ( k α + α ) ρ k α + α 1 L ρ { t k α ρ + α ρ ρ }
From (2) we can get
k = 0 A k Γ ( k α + α ) ρ k α + α 1 ρ k α + α 1 Γ ( k α + α ) s k α + α = 1 s α k = 0 ( A s α ) k = 1 s α ( I + A s α + A 2 s 2 α + ) .
According to the famous Neumann series, we can know if lim n A n = 0 , then I A is nonsingular and ( I A ) 1 = I + A + A 2 + = k = 0 A k . In that case, here Equation (10) can be written in the following form
1 s α ( I + A s α + A 2 s 2 α + ) = 1 s α ( I A s α ) 1 = ( s α I A ) 1
where we can find an s and let lim n ( A s α ) n = 0 , then I A s α is nonsingular. Then
L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) } ( s ) = 1 s α ( I A s α ) 1 = ( s α I A ) 1
when lim n ( A s α ) n = 0 , the above formula holds. □
Definition 2
([7]). Let fand g be two piecewise continuous functions on each interval [ 0 , T ] and of exponential order. We define the ρ-convolution formula of functions f and g is
( f ρ g ) ( t ) = 0 t f ( ( t ρ τ ρ ) 1 ρ ) g ( τ ) d τ τ 1 ρ
The following lemma gives the commutativity of ρ -convolution of two functions.
Lemma 1
([7]). Let fand g be two piecewise continuous functions on each interval [ 0 , T ] and of exponential order. Then
f ρ g = g ρ f
Theorem 4
([7]). Let fand g be two piecewise continuous functions on each interval [ 0 , T ] and of exponential order e c t ρ ρ . Then
L ρ { f ρ g } = L ρ { f } L ρ { g } s > c .
The following theorem is the ρ -Laplace transformation of the left generalized fractional derivative starting from 0.
Theorem 5
([7]). Let α > 0 and f A C γ n [ 0 , a ] for any a > 0 and 0 I n k α , ρ f , k = 0 , 1 , , n 1 be of ρ-exponential order e c t ρ ρ . Then
L ρ { ( 0 D α , ρ f ) ( t ) } ( s ) = s α L ρ { f ( t ) } k = 0 n 1 s n k 1 ( 0 I n k α , ρ f ) ( 0 ) , s > c .

3. Positive Linear Electrical Circuits Described by Generalized Fractional Derivatives

Consider the linear electrical circuit systems described by the left generalized fractional derivatives as the following system equation:
( 0 D α , ρ x ) ( t ) = A x ( t ) + B u ( t ) , 0 < α 1
y ( t ) = C x ( t ) + D u ( t ) ,
where x ( t ) R n is the state vector, u ( t ) R m is the input vector and y ( t ) R p is the output vector, and A R n × n , B R n × m , C R p × n , D R p × m . The initial conditions for (17)and (18) have the form x 0 = ( 0 I 1 α , ρ x ) ( 0 ) , where x 0 is a given vector function.
Theorem 6.
The solution of Equation (17) is given by:
x ( t ) = ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x ( 0 ) + 0 t ( t ρ τ ρ ρ ) α 1 E α , α ( A ( t ρ τ ρ ρ ) α ) B u ( τ ) d τ τ 1 ρ
where
E α , α ( A ( t ρ ρ ) α ) = k = 0 A k t α ρ k Γ ( α k + α ) ρ α k
and E α , α ( A ( t ρ ρ ) α ) is the Mittage-Leffler matrix function, Γ ( x ) = 0 e t t x 1 d t is the gamma function.
Proof. 
Firstly, ρ -Laplace transform is performed on both sides of Equation (17). Then, when n = 1 , through Theorem 5, we can write
L ρ { ( 0 D α , ρ x ) ( t ) } = A L ρ { x ( t ) } + L ρ { B u ( t ) } , s α L ρ { x ( t ) } ( 0 I 1 α , ρ x ) ( 0 ) = A L ρ { x ( t ) } + L ρ { B u ( t ) } ,
Let x 0 = ( 0 I 1 α , ρ x ) ( 0 ) and multiply the identity matrix I on both sides of the above formula, then we have
s α I L ρ { x ( t ) } I x 0 = A L ρ { x ( t ) } + I L ρ { B u ( t ) } , ( s α I A ) L ρ { x ( t ) } = I x 0 + I L ρ { B u ( t ) } , L ρ { x ( t ) } = ( s α I A ) 1 x 0 + ( s α I A ) 1 L ρ { B u ( t ) } ,
From (4) of Theorems 3 and 4, We can write the above form as follows
L ρ { x ( t ) } = ( s α I A ) 1 x 0 + ( s α I A ) 1 L ρ { B u ( t ) } = L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) } x 0 + L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) } L ρ { B u ( t ) } = L ρ { ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x 0 + ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) ρ B u ( t ) } .
Therefore,
x ( t ) = ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x 0 + ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) ρ B u ( t ) = ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x 0 + 0 t ( t ρ τ ρ ρ ) α 1 E α , α ( A ( t ρ τ ρ ρ ) α ) B u ( τ ) d τ τ 1 ρ .
Definition 3
([12]). The linear electrical circuit described by (17) and (18) is called (internally) positive if x ( t ) R + n and y ( t ) R p for t 0 for any initial conditions x 0 R + n for t 0 and u ( t ) R + m , t 0 .
Theorem 7
([12]). The linear electrical circuit system described by (17) and (18) is (internally) positive if and only if
A M n , B R + n × m , C R + p × n , D R + p × m .
Theorem 8.
Let A R n × n and 0 < α 1 . Then, for t 0 ,
E α ( A ( t ρ ρ ) α ) = k = 0 A k t α ρ k Γ ( α k + 1 ) ρ α k R + n × n
and
E α , α ( A ( t ρ ρ ) α ) = k = 0 A k t α ρ k Γ ( α k + α ) ρ α k R + n × n
if and only if A is a Metzler matrix.
Proof. 
From the expansions:
E α ( A ( t ρ ρ ) α ) = I + A t α ρ Γ ( α + 1 ) ρ α + A 2 t 2 α ρ Γ ( 2 α + 1 ) ρ 2 α + , E α , α ( A ( t ρ ρ ) α ) = I Γ ( α ) + A t α ρ Γ ( 2 α ) ρ α + A 2 t 2 α ρ Γ ( 3 α ) ρ 2 α + .
It follows that E α ( A ( t ρ ρ ) α ) R + n × n and E α , α ( A ( t ρ ρ ) α ) R + n × n for small t > 0 only if A is a Metzler matrix.
The sufficiency will be proved by counter evidence. Suppose Ψ = E α , α ( A ( t ρ ρ ) α ) is not positive, that is, there are i and j such that Ψ ( i , j ) < 0 , i j . (27) can be written in the following form
E α , α ( A ( t ρ ρ ) α ) I Γ ( α ) t α ρ Γ ( 2 α ) ρ α = A [ I + A t α ρ Γ ( 2 α ) Γ ( 3 α ) ρ α + ]
Let e i = ( 0 0 1 0 ) T denote a vector whose i-th component is 1 and the rest is 0. Calculate the limit of t 0 + for the above formula, we can get
lim t 0 + Ψ ( i , j ) t α ρ Γ ( 2 α ) ρ α = lim t 0 + e i T E α , α ( A ( t ρ ρ ) α ) I Γ ( α ) t α ρ Γ ( 2 α ) ρ α e j = e i T A e j = a i j
where e i T e j = 0 and e i T Ψ e j = Ψ ( i , j ) .
From Ψ ( i , j ) < 0 , we can get a i j < 0 , i j . Therefore, matrix A is not a Metzler matrix. Thus, A is Metzler matrix, which means E α , α ( A ( t ρ ρ ) α ) > 0 for t 0 . The proof for Equation (26) is similar. □

4. Reachability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives

In this part, since it is independent of the output term, we only need to consider the fractional electrical circuit systems (17).
Definition 4
([12]). The fractional linear electrical circuit described by Equation (17) is called (internally) positive if x ( t ) R + n and all u ( t ) R + m , t 0 .
Theorem 9
([12]). The fractional linear electrical circuit (17) is (internally) positive if and only if
A M n , B R + n × m .
Definition 5
([29]). If there exists the input u ( t ) R + m for t [ 0 , t f ] , t f > 0 , which steers the state of electrical circuit from x ( 0 ) = 0 to the given final state x f R + n , i.e., x ( t f ) = x f , we will call this fractional positive electrical circuit (17), reachable in time [ 0 , t f ] .
Theorem 10.
The fractional linear positive electrical circuit described by Equation (17) is reachable in the time [ 0 , t f ] if and only if the reachability matrix
R α ( t f ) = 0 t f ( t f ρ τ ρ ρ ) α 1 E α , α ( A ( t f ρ τ ρ ρ ) α ) B B T E α , α T ( A ( t f ρ τ ρ ρ ) α ) d τ τ 1 ρ R + n × n
is a monomial matrix.
The input u ( t ) R + m , t [ 0 , t f ] which steers the state of system from x ( 0 ) = 0 to the given final state x f R + n , is given by
u ( τ ) = B T E α , α T ( A ( t f ρ τ ρ ρ ) α ) R α 1 ( t f ) x f R + m , τ [ 0 , t f ] .
Proof. 
The solution of (17) for t 0 has the form (18), let x ( 0 ) = 0 , t = t f then we obtain
x ( t f ) = 0 t f ( t f ρ τ ρ ρ ) α 1 E α , α ( A ( t f ρ τ ρ ρ ) α ) B u ( τ ) d τ τ 1 ρ
It is well-known that R α 1 ( t f ) R + n × n if and only if the matrix (32) is monomial [12]. Substituting (33) into (34) we obtain
x ( t f ) = 0 t f ( t f ρ τ ρ ρ ) α 1 E α , α ( A ( t f ρ τ ρ ρ ) α ) B B T E α , α T ( A ( t f ρ τ ρ ρ ) α ) R α 1 ( t f ) x f d τ τ 1 ρ = 0 t f ( t f ρ τ ρ ρ ) α 1 E α , α ( A ( t f ρ τ ρ ρ ) α ) B B T E α , α T ( A ( t f ρ τ ρ ρ ) α ) d τ τ 1 ρ R α 1 ( t f ) x f = x f
Therefore, the input (33) steers the state of the electrical circuit from x ( 0 ) = 0 to x ( t f ) = x f . □
Theorem 11.
The fractional linear positive electrical circuit (17) is reachable if and only if the following n × n m dimensional matrix is row full rank.
Q r = [ B A B A 2 B A n 1 B ] , r a n k Q r = n .
Proof. 
Using the well-known Cayley-Hamilton theorem it is possible to write the transition matrix in the form
E α , α ( A ( t ρ ρ ) α ) = i = 0 n 1 a i ( t ρ ) A i
where the coefficients a i ( t ρ ) depend on t, and a i ( t ) 0 , i = 1 , 2 , , n 1 .
Substitute the above formula into (34), we can get
x ( t f ) = i = 0 n 1 A i B 0 t f ( t f ρ τ ρ ρ ) α 1 a i ( t f ρ τ ρ ) u ( τ ) d τ τ 1 ρ
where
0 t f ( t f ρ τ ρ ρ ) α 1 a i ( t f ρ τ ρ ) u ( τ ) d τ τ 1 ρ γ i 1 γ i 2 . . . γ i m = γ i , ( i = 0 , 1 , , n 1 ) .
Then
x ( t f ) = i = 0 n 1 A i B γ i = [ B A B A 2 B . . . A n 1 B ] γ 0 γ 1 . . . γ n 1 ,
where we let [ B A B A 2 B A n 1 B ] = Q r . If the system is reachable γ 0 γ n 1 can be obtained from Equation (40), when r a n k Q r = n . □

5. Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives

This section will study the observability of fractional linear electrical circuits described by differential Equations (17) and (18). The positivity of electrical circuits described by (17) and (18) have been explained in definition 3 and theorem 8. The observability of the circuit will be studied below.
Definition 6
([29]). If by knowing the input u(t) and output y(t) for [ 0 , t f ] it is possible to find the unique x ( 0 ) R + n of the electrical circuit, we will call the fractional positive electrical circuit described by (17) and (18) is (strongly) observable in the interval of [ 0 , t f ] .
Theorem 12.
The fractional positive electrical circuit described by (17) and (18) is observable in the interval [ 0 , t f ] if and only if the matrix
W α = 0 t f ( t ρ ρ ) α 1 E α , α T ( A ( t ρ ρ ) α ) C T C E α , α ( A ( t ρ ρ ) α ) d t R + n × n
is a monomial matrix.
Proof. 
Assuming u ( t ) = 0 and we have
y ( t ) = C ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x ( 0 ) .
Using the value of y ( t ) in [ 0 , t f ] , by weighting, i.e., multiply E α , α T ( A ( t ρ ρ ) α ) C T left on both sides of (42), then
( t ρ ρ ) α 1 E α , α T ( A ( t ρ ρ ) α ) C T C E α , α ( A ( t ρ ρ ) α ) x ( 0 ) = E α , α T ( A ( t ρ ρ ) α ) C T y ( t ) .
Integrating (43) on the interval [ 0 , t f ] , we obtain
W α x ( 0 ) = 0 t f E α , α T ( A ( t ρ ρ ) α ) C T y ( t ) d t
and
x ( 0 ) = W α 1 0 t f E α , α T ( A ( t ρ ρ ) α ) C T y ( t ) d t R + n
if and only if the matrix (41) is monomial [29]. □
Theorem 13.
The fractional linear positive electrical circuit system described by (17) and (18) is observable if and only if the following n p × n dimensional matrix is column full rank.
Q o = C C A . . . C A n 1 , r a n k Q o = n
Proof. 
Let u ( t ) = 0 , the solution of the system described by (17) and (18) is
x ( t ) = ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x ( 0 ) ,
y ( t ) = C ( t ρ ρ ) α 1 E α , α ( A ( t ρ ρ ) α ) x ( 0 ) ,
Using the well-known Cayley-Hamilton theorem (37) and substituting (47), we have
y ( t ) = ( t ρ ρ ) α 1 i = 0 n 1 a i ( t ρ ) C A i x ( 0 ) = ( t ρ ρ ) α 1 [ a 0 ( t ρ ) a 1 ( t ρ ) a n 1 ( t ρ ) ] C C A . . . C A n 1 x ( 0 ) .
Let [ C C A C A n 1 ] T = Q o , since a i ( t ρ ) is a known function, the initial state x ( 0 ) can be uniquely determined according to y ( t ) in finite time [ 0 , t f ] if and only if the matrix Q o is column full rank. □

6. Illustrative Examples

In this section, we will present some generalized fractional derivatives circuit systems examples.

6.1. RC Circuit Systems

The circuit shown in Figure 1 is a fractional RC circuit. It includes source e, resistances R 1 ,   R 2 and R 3 and fractional element Capacitor C 1 and C 2 ; both of them are of α . Denote v ( t ) as voltage, let x 1 ( t ) = v C 1 , x 2 ( t ) = v C 2 and u ( t ) = e is the source voltage; u ( t ) and y ( t ) are input and output vectors, respectively. Set its parameters as follows: C 1 = 3 × 10 4 µF, C 2 = 3 × 10 4 µF, R 1 = 10 Ω , R 2 = 20 Ω , R 3 = 30 Ω . Using Kirchhoff’s law, we can write the following systems:
( 0 D α , ρ x 1 ) ( t ) = R 2 + R 3 C 1 R x 1 ( t ) + R 3 C 1 R x 2 ( t ) + R 2 C 1 R u ( t )
( 0 D α , ρ x 2 ) ( t ) = R 3 C 2 R x 1 ( t ) R 1 + R 3 C 2 R x 2 ( t ) + R 1 C 2 R u ( t )
where R = [ R 1 ( R 2 + R 3 ) + R 2 R 3 ] and choose
y ( t ) = x 1 ( t ) + x 2 ( t )
The Systems Equations (49)–(51) can be rewritten into the following system
( 0 D α , ρ x ) ( t ) = A x ( t ) + B u ( t )
y ( t ) = C x ( t )
where
A = R 2 + R 3 C 1 R R 3 C 1 R R 3 C 2 R R 1 + R 3 C 2 R , B = R 2 C 1 R R 1 C 2 R , C = 1 1 .
From (54), we have
A = 1.515 0.909 0.682 0.909 , B = 0.606 0.227 , C = 1 1 .
According to (36),
Q r = B A B = 0.606 0.71175 0.227 0.20695 , r a n k Q r = 2 .
From Theorem 11 it follows that the fractional RC circuit is reachable.
Also, from Theorem 13 it follows that
Q o = C C A = 1 1 R 3 C 2 R R 2 + R 3 C 1 R R 3 C 1 R R 1 + R 3 C 2 R = 1 1 0.833 0 , r a n k Q o = 2
the fractional RC circuit is observable.

6.2. RL Circuit Systems

The circuit shown in Figure 2 is a fractional RL circuit. It includes sources e 1 , e 2 , resistances R 1 ,   R 2 and R 3 , and fractional element Capacitor L 1 and L 2 , both of them are of α . Denote i ( t ) as current, let x 1 ( t ) = i L 1 , x 2 ( t ) = i L 2 and u ( t ) is the source voltage; u ( t ) and y ( t ) are input and output vectors, respectively. Set its parameters as follows: L 1 = 0.6 H , L 2 = 0.8 H , R 1 = 300 Ω , R 2 = 400 Ω , R 3 = 500 Ω .
Using Kirchhoff’s law, we can obtain the following systems equations:
( 0 D α , ρ x 1 ) ( t ) = R 1 + R 3 L 1 x 1 ( t ) + R 3 L 1 x 2 ( t ) + 1 L 1 e 1
( 0 D α , ρ x 2 ) ( t ) = R 3 L 2 x 1 ( t ) R 2 + R 3 L 2 x 2 ( t ) + 1 L 2 e 2
and choose
y ( t ) = R 1 x 1 ( t ) + R 2 x 2 ( t )
The system Equations (58)–(60) can be rewritten into the following systems equations:
( 0 D α , ρ x ) ( t ) = A x ( t ) + B u ( t )
y ( t ) = C x ( t )
where
A = R 1 + R 3 L 1 R 3 L 1 R 3 L 2 R 2 + R 3 L 2 , B = 1 L 1 0 0 1 L 2 , C = R 1 R 2 , u ( t ) = e 1 e 2 .
From (63), we have
A = 1.3 0.8 0.6 1.1 , B = 1.7 0 0 1.25 , C = 300 400 .
According to (36),
Q r = B A B = 1.7 0 2.2 1 0 1.25 1 1.4 , r a n k Q r = 2 .
From Theorem 11 it follows that the fractional RL circuit is reachable.
Also, from Theorem 13 it follows that
Q o = C C A = 300 400 150 200 , r a n k Q o = 2
the fractional RL circuit is observable.

6.3. RC, RL Circuit Systems with Different Parameters ρ

Consider the RC electrical circuit (49–51) described by the left generalized fractional derivative with C 1 = 3 × 10 4 µF, C 2 = 3 × 10 4 µF, R 1 = 10 Ω , R 2 = 20 Ω , R 3 = 30 Ω and e = 6 . The numerical simulations for the circuit with α = 0.85 and ρ = 0.8 are depicted in Figure 3. In the RC circuit described by (49–51) with α = 0.85 , the relationship between the states and different parameters ρ is shown in Figure 4.
Consider the RL electrical circuit (58–60) described by the left generalized fractional derivative with L 1 = 0.6 H , L 2 = 0.8 H , R 1 = 300 Ω , R 2 = 400 Ω , R 3 = 500 Ω , e 1 = 23.56 and e 2 = 5.89 . The numerical simulations for the circuit with α = 0.85 and ρ = 0.8 are depicted in Figure 5. And the relationship between states and different parameters in RL circuit described by (58–60) with α = 0.85 are depicted in Figure 6.

7. Conclusions

In this work, we use the differential equations based on the generalized fractional derivative to describe linear electrical circuits systems, and discuss the positivity of the electrical circuits systems, methods of checking the positivity of the systems are also presented. The main work is to deduce the checking method of reachability and observability of fractional order positive linear electrical circuit systems described by generalized fractional derivative. The ρ -Laplace transform of a Mittage-Leffler function with matrix form is proved, and the conditions of this transformation are also given. At the end of the paper, illustrative examples of RLC, RC and RL electrical circuits systems are presented, the influence of parameter ρ on the state trajectory is observed, and it is found that parameter ρ has an influence on the slope of the state trajectory. The state trajectory of the circuit system described by the generalized fractional derivative is more realistic.

Author Contributions

Conceptualization, H.Y.; Formal analysis, T.Y.; Investigation, T.Y.; Supervision, H.Y. and I.G.I.; Writing—original draft, T.Y.; Writing—review and editing, H.Y. and I.G.I. All authors have read and agreed to the published version of the manuscript.

Funding

This paper is partly financial support by The National Statistical science research project granted 2021LY083 and National Natural Science Foundation of China granted 62176140.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors are thankful to the respected reviewers for their valuable comments and constructive suggestions towards the improvement of the original paper.

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this paper.

References

  1. Sene, N.; Fall, A. Homotopy perturbation ρ-Laplace transform method and its application to the fractional diffusion equation and the fractional diffusion-reaction equation. Fractal Fract. 2019, 3, 14. [Google Scholar] [CrossRef] [Green Version]
  2. Dos Santos, M.; Gomez, I. A fractional Fokker-Planck equation for non-singular kernel operators. J. Stat. Mech. Theory Exp. 2018, 2018, 123205. [Google Scholar] [CrossRef] [Green Version]
  3. Dos Santos, M. Fractional Prabhakar derivative in diffusion equation with non-static stochastic resetting. Physics 2019, 1, 40–58. [Google Scholar] [CrossRef] [Green Version]
  4. Sene, N.; Abdelmalek, K. Analysis of the fractional diffusion equations described by Atangana-Baleanu-Caputo fractional derivative. Chaos Soli. Fract. 2019, 127, 158–164. [Google Scholar] [CrossRef]
  5. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  6. Katugampola, U. A new approach to generalized fractional derivatives. Bul. Math. Anal.Appl. 2014, 6, 1–15. [Google Scholar]
  7. Fahd, J.; Abdeljawad, T. A modified Laplace transform for certain generalized fractional operators. Results Nonlinear Anal. 2018, 1, 88–98. [Google Scholar]
  8. Hristov, J. On the Atangana-Baleanu derivative and its relation to the fading memory concept: The diffusion equation formulation. In Trends in Theory and Applications of Fractional Derivatives with Mittag-Leffler Kernel; Gomez, J., Torres, L., Escobar, R., Eds.; Springer: Cham, Switzerland, 2019; pp. 175–193. [Google Scholar]
  9. Abdeljawad, T. On conformable fractional calculus. J. Comput. Appl. Math. 2013, 279, 57–66. [Google Scholar] [CrossRef]
  10. Khan, H.; Jarad, F.; Abdeljawad, T.; Khan, A. A singular ABC-fractional differential equation with ρ-Laplacian operator. Chaos Solitons Fractals 2019, 129, 56–61. [Google Scholar] [CrossRef]
  11. Sene, N. Fractional input stability for electrical circuits described by the Riemann-Liouville and the Caputo fractional derivatives. AIMS Math. 2019, 4, 147–165. [Google Scholar] [CrossRef]
  12. Kaczorek, T. Selected Problems of Fractional Systems Theory; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  13. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Application of Fractional Differential Equations, Volume 13; North-Holland Mathematics Studies; Elsevier Science: Amsterdam, The Netherlands, 2006; Volume 204, 523p. [Google Scholar]
  14. Kaczorek, T.; Rogowski, K. Fractional Linear Systems and Electrical Circuits; Springer: Cham, Switzerland, 2015. [Google Scholar]
  15. Aguilar, J. Behavior characteristics of a cap-resistor, memcapacitor, and a memristor from the response obtained of RC and RL electrical circuits described by fractional differential equations. Turk. J. Electr. Eng. Comput. Sci. 2016, 24, 1421–1433. [Google Scholar] [CrossRef] [Green Version]
  16. Gomez-Aguilar, J.; Atangana, A. New insight in fractional differentiation: Power, exponential decay and Mittag-Leffler laws and applications. Eur. Phys. J. Plus 2017, 132, 13. [Google Scholar] [CrossRef]
  17. Aguilar, J.; Baleanu, D. Baleanu, Fractional Transmission Line with Losses. Z. Naturforschung A 2014, 69, 539–546. [Google Scholar] [CrossRef] [Green Version]
  18. Gomez-Aguilar, J.; Escobar-Jimenez, R.; OlivaresPeregrino, V.; Taneco-Hernandez, M.; Guerrero-Ramirez, G. Electrical circuits RC and RL involving fractional operators with bi-order. Adv. Mech. Eng. 2017, 9, 1–10. [Google Scholar] [CrossRef]
  19. Morales-Delgado, V.; Gomez-Aguilar, J.; TanecoHernandez, M.; Escobar-Jimenez, R. Fractional operator without singular kernel: Applications to linear electrical circuits. Int. J. Circuit Theory Appl. 2018, 46, 2394–2419. [Google Scholar] [CrossRef]
  20. Morales-Delgado, V.; Gomez-Aguilar, J.; TanecoHernandez, M. Analytical solutions of electrical circuits described by fractional conformable derivatives in Liouville-Caputo sense. AEU-Int. J. Electron. Commun. 2018, 85, 108–117. [Google Scholar] [CrossRef]
  21. Gomez-Aguilar, J.; Morales-Delgado, V. Electrical circuits RC, LC, and RL described by AtanganaBaleanu fractional derivatives. Int. J. Circuit Theory Appl. 2017, 45, 1514–1533. [Google Scholar] [CrossRef]
  22. Radwan, A.; Salama, K. Fractional-order RC and RL circuits. Circuits Syst. Signal Process. 2012, 31, 1901–1915. [Google Scholar] [CrossRef]
  23. Sene, N.; Gomez-Aguilar, J. Analytical solutions of electrical circuits considering certain generalized fractional derivatives. Eur. Phys. J. Plus 2019, 134, 260. [Google Scholar] [CrossRef]
  24. Si, X.; Yang, H. A new method for judgment and computation of stability and stabilization of fractional order positive systems with constraints. J. Shandong Univ. Sci. Technol. (Nat. Sci.) 2021, 40, 12–20. [Google Scholar]
  25. Si, X.; Yang, H. Optimization approach to the constrained regulation problem for linear continuous-time fractional-order systems. Int. J. Nonlinear Sci. Numer. Simul. 2021, 22, 1–16. [Google Scholar] [CrossRef]
  26. Ji, S.; Li, G. Existence of solution to nonlocal fractional differential inclusions via resolvent operators. J. Shandong Univ. Sci. Technol. (Nat. Sci.) 2021, 40, 103–108. [Google Scholar]
  27. Kaczorek, T. Stability of positive continuous-time linear systems with delays. Bull. Pol. Acad. Sci. Tech. Sci. 2009, 57, 395–398. [Google Scholar]
  28. Kaczorek, T. Stability Tests of Positive Fractional Continuous-time Linear Systems with Delays. TransNav Int. J. Mar. Navig. Saf. Sea Transp. 2013, 7, 211–215. [Google Scholar] [CrossRef]
  29. Sene, N. Stability analysis of the generalized fractional differential equations with and without exogenous inputs. J. Nonlinear Sci. Appl. 2019, 12, 562–572. [Google Scholar] [CrossRef]
  30. Yu, Z.; Sun, Y.; Dai, X. Stability and Stabilization of the Fractional-order Power System with Time Delay. IEEE Trans. Circuits Syst. 2021, 68, 3446–3450. [Google Scholar] [CrossRef]
  31. Chao, C.; Chen, D.; Chiou, J. Stability Analysis and Robust Stabilization of Uncertain Fuzzy Time-Delay Systems. Mathematics 2021, 9, 2441. [Google Scholar] [CrossRef]
  32. Sene, N. Stability analysis of electrical RLC circuit described by the Caputo-Liouville generalized fractional derivative. Alex. Eng. J. 2020, 59, 2083–2090. [Google Scholar] [CrossRef]
  33. Erdelyi, A.; Kober, H. Some remarks on hankel transforms. Q. J. Math. 1940, 11, 212–221. [Google Scholar] [CrossRef]
  34. Herrmann, R. Towards a geometric interpretation of generalized fractional integrals—Erdelyi-Kober type integrals on RN as an example. Fract. Calc. Appl. Anal. 2014, 17, 361–370. [Google Scholar] [CrossRef] [Green Version]
  35. Jerzy, K. Controllability of Fractional Linear Systems with Delays. In Proceedings of the 25th International Conference on Methods and Models in Automation and Robotics (MMAR), Międzyzdroje, Poland, 23–26 August 2021; Institute of Electrical and Electronics Engineers: New York, NY, USA, 2021. [Google Scholar]
  36. Gu, P.; Chen, Y.; Tian, S. Learnability of Linear Fractional-Order ILC Systems. IEEE Trans. Circuits Syst. 2021, 68, 963–967. [Google Scholar] [CrossRef]
  37. Kaczorek, T. Reachability and controllability to zero of positive fractional discrete-time systems. In Proceedings of the European Control Conference (ECC), Kos, Greece, 2–5 July 2007; Institute of Electrical and Electronics Engineers: New York, NY, USA, 2007. [Google Scholar]
  38. Kaczorek, T. Invariant properties of positive linear electrical circuit. Arch. Elektrotech. 2019, 68, 875–890. [Google Scholar]
  39. Yuan, T.; Yang, H. Invariance of reachability and observability for fractional positive linear electrical circuit with delays. Arch. Elektrotech. 2021, 70, 513–530. [Google Scholar]
Figure 1. The fractional RC circuit.
Figure 1. The fractional RC circuit.
Mathematics 09 02856 g001
Figure 2. The fractional RL circuit.
Figure 2. The fractional RL circuit.
Mathematics 09 02856 g002
Figure 3. Voltage trajectories across the capacitors described by the fractional RC circuit (49–51).
Figure 3. Voltage trajectories across the capacitors described by the fractional RC circuit (49–51).
Mathematics 09 02856 g003
Figure 4. States with different parameters in RC circuit (49–51): (a) Trajectory of x 1 with different parameters; (b) Trajectory of x 2 with different parameters.
Figure 4. States with different parameters in RC circuit (49–51): (a) Trajectory of x 1 with different parameters; (b) Trajectory of x 2 with different parameters.
Mathematics 09 02856 g004
Figure 5. Current trajectories across the inductors described by the fractional RL circuit (58–60).
Figure 5. Current trajectories across the inductors described by the fractional RL circuit (58–60).
Mathematics 09 02856 g005
Figure 6. Relationship between states and different parameters in RL circuit (58–60): (a) Trajectory of x 1 with different parameters; (b) Trajectory of x 2 with different parameters.
Figure 6. Relationship between states and different parameters in RL circuit (58–60): (a) Trajectory of x 1 with different parameters; (b) Trajectory of x 2 with different parameters.
Mathematics 09 02856 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuan, T.; Yang, H.; Ivanov, I.G. Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives. Mathematics 2021, 9, 2856. https://doi.org/10.3390/math9222856

AMA Style

Yuan T, Yang H, Ivanov IG. Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives. Mathematics. 2021; 9(22):2856. https://doi.org/10.3390/math9222856

Chicago/Turabian Style

Yuan, Tong, Hongli Yang, and Ivan Ganchev Ivanov. 2021. "Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives" Mathematics 9, no. 22: 2856. https://doi.org/10.3390/math9222856

APA Style

Yuan, T., Yang, H., & Ivanov, I. G. (2021). Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives. Mathematics, 9(22), 2856. https://doi.org/10.3390/math9222856

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