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

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.


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).

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: 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: where fractional parameter ρ > 0, n = α + 1 and γ = x 1−ρ d dx . 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.
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 where L{ f } is the usual Laplace transform of f .
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 + bg has Laplace transform for any constants a and b, and has the following linear property The reference [9] gave some ρ-Laplace transforms of elementary functions as shown below.
is the spectral radius of the matrix A. (2) and (3) have been proved, please refer to [9] for details. Only (4) will be proved below From (2) we can get According to the famous Neumann series, we can know if lim In that case, here Equation (10) can be written in the following form where we can find an s and let lim when lim n→∞ ( A s α ) n = 0, the above formula holds.

Definition 2 ([7]
). Let f and 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 The following lemma gives the commutativity of ρ-convolution of two functions.
The following theorem is the ρ-Laplace transformation of the left generalized fractional derivative starting from 0.

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: 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. (17) is given by: (20) and E α,α (A( t ρ ρ ) α ) is the Mittage-Leffler matrix function, Γ(x) = ∞ 0 e −t t x−1 dt is the gamma function.

Theorem 6. The solution of Equation
Proof. Firstly, ρ-Laplace transform is performed on both sides of Equation (17). Then, when n = 1, through Theorem 5, we can write Let x 0 = ( 0 I 1−α,ρ x)(0) and multiply the identity matrix I on both sides of the above formula, then we have From (4) of Theorems 3 and 4, We can write the above form as follows Therefore,

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.

Proof. From the expansions:
It follows that E α (A( t ρ ρ ) α ) ∈ R n×n 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 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 where e T i e j = 0 and e T i Ψe j = Ψ(i, j). From Ψ(i, j) < 0, we can get a ij < 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.

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). (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
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 ].
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 Proof. The solution of (17) for t ≥ 0 has the form (18), let x(0) = 0, t = t f then we obtain 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 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 × nm dimensional matrix is row full rank.
Proof. Using the well-known Cayley-Hamilton theorem it is possible to write the transition matrix in the form 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 Then where we let [B AB A 2 B · · · A n−1 B] = Q r . If the system is reachable γ 0 ...γ n−1 can be obtained from Equation (40), when rank Q r = n.

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 is a monomial matrix.
Proof. Assuming u(t) = 0 and we have 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 Integrating (43) on the interval [0, t f ], we obtain 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 np × n dimensional matrix is column full rank.
Proof. Let u(t) = 0, the solution of the system described by (17) and (18) is Using the well-known Cayley-Hamilton theorem (37) and substituting (47), we have Let [C CA · · · CA 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.

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

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: The Systems Equations (49)-(51) can be rewritten into the following system where According to (36), From Theorem 11 it follows that the fractional RC circuit is reachable. Also, from Theorem 13 it follows that the fractional RC circuit is observable.

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: and choose The system Equations (58)-(60) can be rewritten into the following systems equations: where the fractional RL circuit is observable.

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.