Next Article in Journal
Bimodal-Distributed Binarized Neural Networks
Next Article in Special Issue
Unilateral Laplace Transforms on Time Scales
Previous Article in Journal
Machine Learning Combined with Restriction Enzyme Mining Assists in the Design of Multi-Point Mutagenic Primers
Previous Article in Special Issue
Reachability and Observability of Positive Linear Electrical Circuits Systems Described by Generalized Fractional Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Systems’ Identification Based on Implicit Modulating Functions

by
Oliver Stark
1,2,3,*,
Marius Eckert
2,3,4,
Albertus Johannes Malan
2 and
Sören Hohmann
2
1
Sartorius Stedim Bitech GmbH, 37099 Göttingen, Germany
2
Faculty of Electrical Engineering and Information Technology, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
3
Faculty of Electrical Engineering and Information Technology, Institute of Control Systems (IRS), 76131 Karlsruhe, Germany
4
Robert Bosch GmbH, 70839 Gerlingen, Germany
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(21), 4106; https://doi.org/10.3390/math10214106
Submission received: 27 September 2022 / Revised: 23 October 2022 / Accepted: 24 October 2022 / Published: 3 November 2022
(This article belongs to the Special Issue Fractional Differential Equations and Control Problems)

Abstract

:
This paper presents a new method for parameter identification based on the modulating function method for commensurable fractional-order models. The novelty of the method lies in the automatic determination of a specific modulating function by controlling a model-based auxiliary system, instead of applying and parameterizing a generic modulating function. The input signal of the model-based auxiliary system used to determine the modulating function is designed such that a separate identification of each individual parameter of the fractional-order model is enabled. This eliminates the shortcomings of the common modulating function method in which a modulating function must be adapted to the investigated system heuristically.

1. Introduction

Fractional calculus is derived from the field of mathematics and the research on this field is still ongoing [1,2]. In recent years, fractional calculus is being used more and more in the engineering field. Complex physical, chemical e.g., Lithium-ion batteries [3], or biological systems e.g., blood alcohol model [4] are increasingly being described by fractional-order models. Fractional differential equations are also used to describe parts of electrical circuits or refine the description of friction equations in mechanical systems. An overview of technical areas is provided by [5]. Due to the non-locality and memory of fractional integration and differentiation, these models approximate partial differential equations more accurately than classic integer-order models [6]. Over the last years, identification methods considering fractional-order models based on the modulating function method (see e.g., [7,8,9,10,11,12,13,14]) have become more prevalent. A benefit of the modulating function method is that a system of algebraic equations may be solved for the parameter identification instead of a set of differential equations and the measured signal need not be differentiated [15]. In current approaches, a generic modulating function is chosen and parameterized heuristically. This procedure constitutes time-consuming educated guessing which does not generalize or provide clues for new applications, especially in the fractional order case. Thus, parameterizing the modulating function is a significant shortcoming of this method. For parameter identification, the algebraic equations are collected in a linear system which has to be solved. Because the modulating function is chosen and parameterized a-priori, linear dependencies between the individual algebraic equations can occur, depending on the measured signal. In this case, the parameter identification cannot be performed, which constitutes is a significant shortcoming of this method.
In this paper, we propose a systematic procedure for parameter identification using an implicit determination of a modulating function considering the current measurements. This procedure allows for direct parameter identification, without requiring the explicit calculation of modulating function. To the best of the authors’ knowledge, it is the first time that a method that automatically determines a modulating function for the parameter identification of a fractional-order model is presented. While this idea has been described for integer-order models [16], this result does not easily generalize to the non-integer order case.
The paper is structured as follows: In Section 2, the basics of fractional calculus are provided. In Section 3, the modulating function method is recapped and an additional property for the modulating function method, enabling the separate identification of the parameters, are given. Afterward, the model-based auxiliary system, which is central in determining a modulating function automatically, is derived in Section 4. The model-based auxiliary system is used to transfer the heuristic determination of a modulating function into a control problem which allows an automatic determination of the parameter-specific modulating function. A solution to the control problem is provided in Section 5, which leads to the parameter identification and represents the main contribution of this paper. Additionally, the error occurs due to superposed noise on the input and output signal is analyzed. A numerical example in Section 6 demonstrating the efficacy of the proposed algorithm concludes the paper.

2. Preliminaries

2.1. Fundamentals of Fractional Calculus

Throughout this paper, consider the smooth functions f : ( , t ] R and h : [ t , ) R with the properties f ( t ) = 0 , t t 0 and h ( t ) = 0 , t t e . Let α R > 0 and t 0 t 1 < t < t 2 t e hold. Furthermore, · describes the floor function and denotes the biggest integer smaller or equal to the argument and is used to state a modified ceiling-function · : = · + 1 . The functions f and h are assumed to be Lebesgue integrable on the integration interval defined by t 0 , t e (see [17]). In addition, the assumption is made that f and h are α -times absolutely continuous on the derivation interval t 0 , t e , where α is the derivation order (see [17]).
Remark 1. 
With regard to the parameter identification, the time variables t 0 , t 1 , t 2 , t e , and t are interpreted as follows. While t 0 describes the point of time before the system is at rest, t e describes the point of time after the system is at rest. The time variable t 1 marks the beginning of the identification, t 2 notes the end of the identification, and t is the independent integration variable.
The uninitialized as well as initialized fractional operators are described in [18] and are given in the following definition. The initialized fractional operators consist of the uninitialized fractional operator and time-variant initialization functions η : R 5 R or ψ : R 5 R , which have been proven to be necessary in order to describe fractional differential equations (FDEs) correctly. Differing from the notation used in [18], in this paper, an additional operator index on the top left of the fractional operator is used (see [19]). This index indicates that the function is integrated or differentiated w.r.t. the named variable. The index on the bottom left represents the lower bound of the fractional integral and the index on the bottom right is the upper bound. The order of the integration or derivative is given by the top right index. The calculus of two-variable functions can be directly extended from the fractional operators for functions depending on one variable (see [20]).
In the following definition, the uninitialized and the initialized fractional integration along with the Riemann-Liouville ( R L ) and the Caputo (C) fractional derivatives according to [21] are summarized.
Definition 1. (Fractional Operators.) 
  • Uninitialized Fractional Integral
    i t α t 0 τ f ( τ ) : = 1 Γ ( α ) t 0 t f ( ν ) ( t ν ) ( 1 α ) d ν
  • Initialized Fractional Integral
    I t α t 1 τ f ( τ ) : = i t α t 1 t f ( t ) + η f , α , t 0 , t 1 , t
  • Uninitialized Riemann-Liouville Fractional Derivative
    d R L t α t 0 τ f ( τ ) : = d d t α i t α α t 0 τ f ( τ )
  • Uninitialized Caputo Fractional Derivative
    d C t α t 0 τ f ( τ ) : = i t α α t 0 τ d d τ α f ( τ )
  • Initialized Fractional Derivative
    D t α t 1 τ f ( τ ) : = d t α t 1 τ f ( τ ) + ψ f , α , t 0 , t 1 , t
where Γ is the Gamma function (see [22] (pp. 1–6)).
Remark 2. 
The definition of the Riemann-Liouville (3) and the Caputo (4) fractional derivative are connected by
d R L t α t 0 τ f ( τ ) = d C t α t 0 τ f ( τ ) + k = 0 α d R L t α t 0 τ f k ( τ ) k ! τ = t 0 t t 0 k
which is given in [17].
Remark 3. 
In (2) and in (5), the application of the initialization function η or ψ necessitates knowledge of the function f in the interval t t 0 , t 1 (see [18]). To ensure that a system is at rest, the initialization function implies that f has to be known from t 0 according to [18]. Therefore, an exact initialization is not possible in practical applications.
Remark 4. 
In (5), the specific operator depends on the utilized uninitialized operator (3) or (4). Thus, D t α t 1 t f t may be described by either D R L t α t 1 t or D C t α t 1 t and must correspond to the uninitialized operator (see [18]). According to the uninitialized operator, the related initialization function has to be employed (see [18]).
In addition to the left-sided definitions, the equivalent right-sided definitions of the integrations (1) and (2) and derivatives (3)–(5) exist.
Definition 2. (Uninitialized Right-Sided Fractional Integral.)
i t e α t τ h ( τ ) : = 1 Γ ( α ) t t e h ( ν ) ( ν t ) ( 1 α ) d ν .
The other definitions follow directly by inserting (7) instead of (1) into the definitions (2)–(5). These are explicitly given in Appendix A. The left- and right-sided fractional operators are connected by a reflection operator Q .
Definition 3. (Reflection Operator.)
The reflection operator Q maps a left-sided function f onto a right-sided one f ˜ (see [23]):
Q f t : = f t 1 + t 2 t = f ˜ t .

2.2. Fractional-Order Models

We define the commensurable fractional-order model using the initialized Caputo fractional derivative (5)
D C t 2 n α t t y ˜ o t + i = 0 n 1 a i D C t 2 i α t t y ˜ o t = j = 0 m b j D C t 2 j α t t u ˜ o t
where a i , b j R are unknown parameters. The number of parameters n , m N 0 can either be determined by white-box modeling or be specified for black-box modeling. We thus assume n and m are known. The fractional order α R > 0 is usually unknown and has to be identified as well.
The input u ˜ o and output signal y ˜ o in (9) are a noisy observation of the undisturbed input u
u ˜ o t = u ˜ t + u ˜ e t
and of the undisturbed output y
y ˜ o t = y ˜ t + y ˜ e t ,
to which the reflection operator (8) has been applied. The measurement disturbances are represented by u ˜ e and y ˜ e .
Remark 5. 
The reflected signals are used instead of the original signals since this allows a left-sided auxiliary system to be derived (see Section 4).
Regarding the input signal u ˜ , we assume that the system is persistently excited, which is a necessary condition for parameter identification [24] (p. 250). In Appendix C, we show that the condition of persistent excitation collapses into a necessary condition on the signal u ˜ (see Lemma A3).
Assumption A1. (Known Fractional Order.)
In this paper, it is assumed that the order is provided by another method like in [25,26]. From this assumption, it directly follows that the fractional orders are known for parameter identification. Furthermore, we consider stable systems and, because of the extended Matignon’s theorem (see [26]), α ( 0 , 2 ) is assumed.
Despite the restrictions placed by Assumption A1, we note that the derivations are not limited solely to this case.

3. Parameter-Specific Identification Using Modulating Function Method

The modulating function method is a well-known method for the parameter identification of fractional-order models (see e.g., [7,8,9,10,11,12]). The modulating function method was derived in [27] for the integer-order case and has been transferred to the fractional-order case (see e.g., [7,8,9,10,11,12]). In Section 3.1, we shortly recap the modulating function approach to introduce the necessary notation. Based on the recap, a parameter-specific modulating function is defined in Section 3.2. The properties of such a parameter-specific modulating function are used to transfer the identification problem into a control problem in Section 5.

3.1. Fractional Modulating Function Method

Suppose, a modulating function φ exists which fulfills Assumption A2.
Assumption A2. (Properties of Modulating Functions.)
Suppose J : = i N > 0 i n and K : = k N > 0 i J : k i α . The modulating function φ : R R fulfills the properties:
( P 1 ) i : i φ t C n α ( [ t 1 , t 2 ] ) , ( P 2 ) i : i d d t k i t i α i α t 1 t φ t t = t 1 = d d t k i t i α i α t 1 t φ t 2 t = t 1 = 0 and ( P 3 ) i : i d R L t i α t 1 t φ ( t ) = 0
where C indicates the class of continuously differentiable functions, i J , and k K .
The property P 3 ensures that the influence of the initialization function ψ in (9) is eliminated. Multiplying (9) with a modulating function φ and integrating the resulting equation by parts results in
i = 0 n a i t 1 t 2 y ˜ o t d R L t i α t 1 t φ t d t = j = 0 m b j t 1 t 2 u ˜ o t d R L t j α t 1 t φ t d t .
without loss of generality we assume a n = 1 . This results in the well-known lemma for the identification of the unknown parameters (see [12]). Note, the mathematical operations changes the Caputo initialized fractional derivative used in (9) into uninitialized Riemann-Liouville fractional derivative (3) in (12).
Lemma 1. (Parameter Identification Applying Modulating Function Method.)
Suppose t , t 1 R , T Δ , T I R > 0 , the modulating-function φhas property P 3 , the system
D C t 2 n α t t y ˜ o t + i = 0 n 1 a i D C t 2 i α t t y ˜ o t = j = 0 m b j D C t 2 j α t t u ˜ o t ,
h N 0 , h n + m + 1 , L : = l N 0 l h , l L , t l = t 1 + l T Δ , and p : = a 0 a n 1 b 0 b m . Furthermore, suppose  M : = m 0 m l m h , where m l is given by
m l : = t l t l + T I y t φ t d t t l t l + T I y t d R L t ( n 1 ) α t 1 t φ t d t t l t l + T I u t φ t d t t l t l + T I u t d R L t m α t 1 t φ t d t .
Is M regular, we get the parameters out of
z = M p
where z : = z 0 z l z h and
z l : = t l t l + T I y t d R L t n α t 1 t φ t d t .
Proof. 
The proof can be found in [12]. □
We note that, depending on the a-priori chosen modulating functions and their parameterization, the equations of system (14) can be linearly dependent, which presents an obstacle for the identification of parameters.

3.2. Parameter-Specific Modulating Function

Sometimes, it is favorable to identify just one parameter at a time. Therefore, we introduce a parameter-specific modulation function approach. Suppose m o t 2 = e s holds true for each parameter p s , where p s represents the s-th element of the parameter vector p : = a 0 a n 1 b 0 b m . Call φ s parameter-specific modulating. In summary, in line with the assumption P 1 P 3 we define:
Definition 4. (Parameter-Specific Modulating Function Set.)
Suppose J : = i N > 0 i n and K : = k N > 0 i J : k i α . The parameter-specific set W s is defined as follows:
W s : = W 1 W 2 W 3 s
W 1 : = φ s : R R φ s t = 0 t t 1
W 2 : = φ s : R R i J k K : d R L t i α + k i α t 1 t φ s ( t ) t = t 2 = 0
W 3 s : = φ s : R R m o t 2 = e s
where the fractional derivative is given in (3),  m o is given by
m o t 2 : = t 1 t 2 y ˜ o t φ s t d t t 1 t 2 y ˜ o t d R L t ( n 1 ) α t 1 t φ s t d t t 1 t 2 u ˜ o t φ s t d t t 1 t 2 u ˜ o t d R L t m α t 1 t φ s t d t .
In (20), e s R ( n + m + 1 ) × 1 consists of the elements
e s j = 1 for j = s 0 e l s e
and · j describes the j-th element of a vector, j S and S : = s N > 0 s n + m + 1 .
Theorem 1. (Parameter-Specific Identification.)
Suppose S : = s N > 0 s n + m + 1 , φ s W s , the fractional derivative (3) and p s : = p s where · s describes the s-th element of a the parameter vector p : = a 0 a n 1 b 0 b m .
The s-th parameter of (12) is calculated using a parameter-specific modulating function φ s as follows:
t 1 t 2 y ˜ o t d R L t n α t 1 t φ s t d t = p s .

4. Model-Based Auxiliary System

Previous approaches using the modulating function method require such a function to be known a-priori. As motivated in the introduction and in Section 3, finding a function that meets P 1 P 3 is often a tedious task, especially for fractional systems.
In the following two sections, we, therefore, present our main result which allows parameter identification without requiring an explicit a-priori modulating function. Instead, we propose in this section the use of an auxiliary system that implicitly contains the requirements of a modulating function. In the sequel, we then demonstrate how applying an appropriate control to this dynamical system can be used to automatically retrieve a valid modulating function and achieve parameter identification at the same time.
In Section 4.1, we introduce notation for the compact representation of an auxiliary system. The auxiliary system used for parameter identification is formally introduced in Section 4.2.

4.1. Notations for the Model-Based Auxiliary System

Before defining the model-based auxiliary system, we here introduce notations that allow for a more compact description of the requirements needed for such a system. Afterward, we introduce additional notations regarding the dimensions of subsystems contained in the model-based auxiliary system and we define normalization parameters.
The main idea behind the model-based auxiliary system is that every expression of the modulating function can be represented as a combination of the n α -th derivative d R L t n α t 1 t φ t and a fractional integration with corresponding order. For example, (19), considering i J , results in
d R L t i α + k i α t 1 t φ s ( t ) = i t λ α + ( n λ ) α k t 1 t d R L t n α t 1 t φ t
where the fractional derivative is given in (3),  λ = 1 , 2 , , n , and k = 0 , 1 , , n λ α . To enable a compact definition of the model-based auxiliary system, a vector collecting all possible derivative orders of (19) is given first.
Definition 5. (Vector of Boundary Term Orders.)
γ : = β 1 β 2 β n
for the derivative orders
β k : = n k α + k α n k α + k α n k α 1 + k α . . . 1 + k α .
Furthermore, suppose the two normalization parameters w : = min k N > 0 α k 1 and r : = min k N > 0 max γ k 1 , where γ is defined by (25). This normalization parameters ensure that all derivative orders of the model-based auxiliary system fulfill the requirement α ( 0 , 1 ] for fractional state spaces (see [28]). Additionally, κ : = n + m + 1 , C : = k = 1 n n k α , and z : = r C are defined for a more compact representation of the subsystem dimensions.

4.2. Model-Based Auxiliary System

Equipped with the notations introduced in the previous section, we now define the model-based auxiliary system. The model-based auxiliary system is constructed by interconnecting the derivatives of the modulating function d R L t i α t 1 t φ t and d R L t j α t 1 t φ t (see (12)), the identification equation of the modulating function method (12), and the resulting boundary terms P 1 P 3 . In this subsection, we start by defining the states associated with each of these parts and proceed with describing their respective state dynamics.
Remark 6. 
We denote the four interconnected parts of the model-based auxiliary system as its constitutive subsystems and assign each a specific identifier. Subsystem describes the equation resulting of the modulating function method (12), subsystem maps the connection between the derivatives of the modulating function in (12). Finally, subsystems Δ and take the derivatives in P 2 into account.
Definition 6. (Fractional State Vectors.)
Suppose x : R R κ × 1 , x : R R n w × 1 , x Δ : R R n α × 1 , the fractional derivative (3), and x : R R z × 1 .
x t = x t x t x Δ t x t
where
x t : = m o t
x t : = x 0 t x n 1 t , x k t : = d R L t k α t 1 t φ t d R L t α w + k α t 1 t φ t d R L t ( w 1 ) α w + k α t 1 t φ t
x Δ t : = d R L t n α n α t 1 t φ t d R L t n α + 1 n α t 1 t φ t d R L t n α 1 t 1 t φ t
x t : = x 1 t x C t , x k t : = x γ k d R L t γ k r t 1 t x γ k d R L t ( r 1 ) γ k r t 1 t x γ k
and
x γ t : = d R L t γ t 1 t φ t .
In Definition 6, the states ∘ represent the derivatives of the modulating function which occur in the basic equation of the modulating function method (12). Hence, the fractional state equations for the subsystem ∘ are as given in the following Lemma 2.
Lemma 2. (Subsystem ∘.)
Suppose A : R R κ × n w and ν = n or ν = m + 1 , depending on whether the output or input signal is considered. The subsystem is described using the state vector of subsystem □:
x ˙ t : = y ˜ t · H n × n w u ˜ t · H ( m + 1 ) × n w x t
= A t x t
where
H : = h 1 h ν , h k l = 1 if l = ( k 1 ) w + 1 0 e l s e .
Proof. 
Using the definition (21), it follows directly that
x ˙ t = y ˜ o t φ s t u ˜ o t d R L t m α t 1 t φ s t .
The expressions related to the modulating function are equivalent to the 1 ( st ) , ( w + 1 ) ( th ) , ⋯, ( ν 1 ) w + 1 ( th ) elements of x t which is selected by (35). □
The connection between the derivatives of the modulating function in (12) are mapped by subsystem □ considering that all derivatives are represented as a combination of the n α -th derivative d R L t n α t 1 t φ t and a fractional integration with corresponding order.
Lemma 3. (Subsystem □.)
Suppose A R n w × n w and b R n w × 1 . Then,
d R L t α w t 1 t x t : = J n w × n w x t + 0 n w 1 × 1 1 u * t
= A x t + b u * t
where J is a Jordan matrix with dimensions n w × n w  which only has eigenvalues of 0.
Proof. 
The connection of the derivatives of the modulating function is a chain of integrators for the input u * t = d R L t n α t 1 t φ t , whereby each integrator is of order α w . This chain can then be written as a system with a Jordan matrix where all eigenvalues are zero. □
While subsystem Δ takes the boundary terms (19) for i = n into account, subsystem represents the boundary terms for i = 0 , 1 , n 1 .
Lemma 4. (Subsystem Δ.)
Suppose A Δ R n α × n α and b Δ R n α × 1 . Then,
x ˙ Δ t : = J n α × n α x Δ t + 0 n α × 1 1 u * t
= A Δ x Δ t + b Δ u * t
where J is a Jordan matrix with dimensions n α × n α and has only eigenvalues of 0.
Proof. 
The proof is analogous to the proof of Lemma 3 with the difference that the order of each integrator is 1. □
Lemma 5. (Subsystem .)
Suppose A k R r × r , b k R r × 1 , A R z × z and b R z × 1 . The subsystem can be subdivided into C subsystems of dimensions r × r :
d R L t γ k r t 1 t x k t : = J r × r x k t + 0 ( r 1 ) × 1 1 u * t
= A k x k t + b k u * t
where J is a Jordan matrix with dimensions r × r  which only has eigenvalues of 0. Connecting all subsystems results in
d R L t γ r t 1 t x t : = A 1 0 r × r 0 r × r A C x t + b 1 b C u * t
= A x t + b u * t
Proof. 
The proof is analogous to the proof of Lemma 3, where the order of each integrator is γ k w . □
Finally, using the defined states and their dynamics comprising the model-based auxiliary system, we can now construct a full system description. We achieve this by simply combining all matrices and vectors of the respective subsystems into the time-variant system matrix
A t : = 0 κ × κ A t 0 κ × n α 0 κ × z 0 n w × κ A 0 n w × n α 0 n w × z 0 n α × κ 0 n α × n w A Δ 0 n α × z 0 z × κ 0 z × n w 0 z × n α A
and the input vector of the model-based auxiliary system
b : = 0 κ × 1 b b Δ b
with the input signal
u * t : = d R L t n α t 1 t φ t .
Remark 7. 
The model-based auxiliary system is derived without any restrictions on the fractional order. Hence, the system must be also valid for the integer-order case. Considering α = 1 , the fractional states of subsystem Δ and subsystem are equivalent to the fractional states of subsystem and the model-based auxiliary system reduces to the subsystem and subsystem □. This result is in line with the results in [16].

5. Implicit Determination of the Modulating Function

In this section, it is shown that the model-based auxiliary system can be used to determine a parameter-specific modulating function φ s . For this, in Section 5.1 it is shown that the model-based auxiliary system represents the parameter-specific modulating function set given by Definition 4. Due to this connection, the determination of a parameter-specific modulating function φ s is transferred into a control problem. In Section 5.2, the identification error resulting from a noisy observation of the input and output signals is analyzed.

5.1. Control-Based Identification

The connection between the parameter-specific modulating function set and the model-based auxiliary system (45)–(47) is given in the following theorem.
Theorem 2. (Representation of the Parameter-Specific Modulating Function Set.)
If a model-based auxiliary system as stated in Section 4.2 exists and if a control input u * can be calculated such that the uninitialized model-based auxiliary system
x t δ ( t t 1 ) ψ x , α w , t 0 , t 1 , t x Δ t δ ( t t 1 ) ψ x , γ r , t 0 , t 1 , t = 0 κ × 1 0 n w × 1 0 n α × 1 0 z × 1
can be transferred into the final state
x t 2 x t 2 x Δ t 2 x t 2 = e s 0 n w × 1 0 n α × 1 0 z × 1 ,
then the final state (49) represents the parameter-specific modulating function set W s .
Proof. 
Due to (48),
x t δ ( t t 1 ) ψ x , α w , t 0 , t 1 , t x Δ t δ ( t t 1 ) ψ x , γ r , t 0 , t 1 , t = 0 κ × 1 0 n w × 1 0 n α × 1 0 z × 1 φ s t ˜ = 0 t ˜ t 1
holds and fulfills (18). The final state (49) fulfills the terminal property (19):
x Δ t 2 = 0 n α × 1 d R L t n α + k n α t 1 t φ s t t = t 2 = 0
where k = 0 , 1 , , n α and
x t 2 = 0 z × 1 d R L t i α + k i α t 1 t φ s t t = t 2 = 0
where k = 0 , 1 , , n α and i = 0 , 1 , , n 1 .
In addition to Definition 4, the final state also ensures that all derivatives of the modulated functions which arise due to the normalization fulfill the terminal property (19)
x t 2 = 0 n w × 1 d R L t k α w t 1 t φ s t t = t 2 = 0
where k = 0 , 1 , , n w 1 . Also, resulting from the final state (49)
x t 2 = e s m o t 2 = e s .
Hence, (20) is fulfilled as well. □
Following Theorem 2, the heuristic determination of a parameter-specific modulating function φ s can be considered as a control problem. It is sufficient to find an input signal u * that steers the auxiliary system from the uninitialized state (48) into the final state (49). Such an input signal
u s * t : = d R L t n α t 1 t φ s t
is called a parameter-specific control input and marked with the lower index s. This enables the replacement of the parameter-specific modulating function d R L t n α t 1 t φ s t by the parameter-specific control input u s * in the parameter-specific identification (23).
Lemma 6. (Control-Based Identification.)
Suppose u s * t = d R L t n α t 1 t φ s t . The parameter-specific identification (23) is equivalent to
t 1 t 2 y ˜ o t u s * t d t = p s .
Proof. 
By inserting (55) into (23), (56) follows directly. □
Remark 8. 
The parameter-specific modulating function can be derived from (55) by n α -times integration:
φ s t = i t n α t 1 t u s * t .
Because the identification problem is transferred into a control problem and to attain a final state (49), the model-based auxiliary system has to be controllable. The controllability is investigated in Appendix C.

5.2. Analysis of the Identification Error

In the following, the identification error, which occurs if the input and the output signal are superposed by noise (10) and (11), is calculated. It is also shown that the energy-optimal control stated in Appendix B minimizes the upper bound of the identification error. For this, it is assumed that a parameter-specific input signal u s * exists and, hence, the identification Equation (56) holds true. The remaining derivatives of the parameter-specific modulating function are provided by the model-based auxiliary system.
Lemma 7. (Identification Error.)
Suppose a parameter-specific control input u s * and the parameter-specific identification Equation (23). Further, the notation p ^ s indicates the parameter estimate for noisy observations of the output signal. The identification error between the estimated parameter p ^ s and the original parameter p s of the fractional-order model is
Δ p s : = p ^ s p s
Δ p s = t 1 t 2 y ˜ e t u s * t d t m e t 2 p
where m e is given by
m e t 2 : = t 1 t 2 y ˜ e t φ s t d t t 1 t 2 y ˜ e t d R L t n 1 α t 1 t φ s t d t t 1 t 2 u ˜ e t φ s t d t t 1 t 2 u ˜ e t d R L t m α t 1 t φ s t d t .
Proof. 
Due to the assumption that a parameter-specific control input u s * exists, the control input u s * may be used instead of the parameter specific modulating function φ s in (23). Assuming the observation of the input and output signals are superposed by additive noise as in (10) and (11), (59) can be separated into a summand which depends on the measured signal and a summand which depends on the disturbance. This leads to
t 1 t 2 y ˜ o t u s * t d t t 1 t 2 y ˜ e t u s * t d t = m o t 2 m e t 2 p
where m o is given by (21). Considering W 3 s for the measured data y ˜ o t and u ˜ o t , (61) results in
p ^ s t 1 t 2 y ˜ e t u s * t d t = p s m e t 2 p .
Rearranging (62) with regard to (58) completes the proof. □
The identification error (59) depends on the input signal u s * of the model-based auxiliary system. Therefore, the choice of the input signal influences the identification quality.
Lemma 8. (Minimal Bound of Identification Error.)
Suppose the identification error Δ p s (59) and a parameter-specific control input u s * . Using the energy-optimal control (A5) to determine the parameter-specific control input u s * leads to a minimal upper bound of the identification error
| Δ p s | P t 1 , t 2 , p , y ˜ e , u ˜ e · 2 Γ ( 2 ς ) i t 2 2 2 ς t 1 t J ( u s * , t 1 , t 2 )
where P t 1 , t 2 , p , y ˜ e , u ˜ e is given by
P t 1 , t 2 , p , y ˜ e , u ˜ e = 1 + i = 1 n ( t 2 t 1 ) 2 i Γ ( i ) ( 2 i 1 ) 2 i | p n i + 1 | t 1 t 2 | y ˜ e t | 2 + j = n m n ( t 2 t 1 ) 2 j Γ ( j ) ( 2 j 1 ) 2 j | p 2 n j + 1 | t 1 t 2 | u ˜ e t | 2
and ς : = min α w , γ r .
Proof. 
For the sake of clarity, the proof of Lemma 8 is given in Appendix B. □
In summary, if a modulating function φ belongs to the parameter-specific set W s , each parameter of a fractional-order model (9) can be identified separately. The automatic determination of such a parameter-specific modulating function φ s is transferred into a control problem by formulating a model-based auxiliary system (45)–(47). Notethat the parameter identification can directly be performed using the control input (55) and that no explicit calculation of the modulating function is necessary. The model-based auxiliary system interconnects the derivatives of the modulating function, the identification equation of the modulating function method (12), and the resulting boundary terms. For a parameter-specific identification to take place, the model-based auxiliary system must be steered to the final state (49) in Theorem 2. The modulating function, and thus the parameters, automatically adjusts to new measurements since the control input continually adjusts to ensure the final state (49) is achieved. Furthermore, applying an energy-optimal control to the model-based auxiliary system reaches the final state while minimizing the upper bound of the parameter identification error in the presence of noisy input and output signals.

6. Numerical Example

The parameter identification based on implicit modulating functions is illustrated in the following. For this purpose, we consider the system
D C t 2 α t t y ˜ o t + a 0 y ˜ o t = b 0 u ˜ o t
where n = 1 , m = 0 , and α = 0.35 are assumed to be known, and a 0 = 2 as well as b 0 = 1 are unknown. The simulation is started at t 0 = 0 s and the identification at t 1 = 5 s . The duration of simulation is 15   s with a sampling time of T s = 0.01   s and the duration of identification is 5   s . A mean-free, pseudo-random binary sequence with an amplitude of 1 is used as an input signal u t [24] (p. 165). The input, output and noisy output ( SNR = 34.7   dB ) signals are shown in Figure 1. A green dotted line marks the starting time of the identification.
To state the model-based auxiliary system, the fractional orders of the boundary terms must first be calculated. Evaluating (25) yields γ = 1.35 and r = 2 . Because n α = 1 , the following fractional pseudo state space yields:
x ˙ t d R L t α t 1 t x t x ˙ Δ t d R L t γ r t 1 t x t = 0 0 y ˜ t 0 0 0 0 0 u ˜ t 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 x t x t x Δ t x t + 0 0 1 1 0 1 u * t
Depending on the requested parameter, the final pseudo state for a 0 is
x ( t 2 ) = 1 0 0 0 0 0
and for b 0
x ( t 2 ) = 0 1 0 0 0 0 .
The parameter is identified by evaluating (23) and (56) respectively. To calculate the parameter-specific control input, (A32) has to be evaluated. Therefore, the matrix approach of [29] is used.
Because one measurement of the input and output signals can be used to identify all parameters, the parameter-specific control input as well as the trajectories of the parameter estimates for a 0 and b 0 are illustrated in Figure 2 for the noise-free case and in Figure 3 for the case with a noisy observation of the output signal. The final estimations a ^ 0 = 2 and b ^ 0 = 1 in the absence of noise and a ^ 0 = 2.005 and b ^ 0 = 0.999 in the case with a noisy observation of the output signal can be read at the final time of the identification t 2 = 10   s .
The results are compared to the parameter identification using a heuristically adapted modulating function. Because the system is not at rest when the identification is started, the spline-type modulating function has to be used (see [12]). For the approach which is described in this paper, only the starting time and the duration of the identification may be chosen freely. In addition to the start and end times, the use of spline-type modulating functions requires that the number of splines and the order of the modulating function is chosen. It should be noted that the maximum order of the modulating function depends on the number of splines and the minimum order depends on the fractional order of the system (see [12]). All the parameters of the spline-type modulating function have to be fixed a-priori. This is not necessary using the newly proposed method based on implicit modulating functions, from which the modulating function can be determined automatically.
In this example, the method described in [12] is applied using a spline-type modulating function with 20 splines and of order 5 for the comparison. The starting time t 1 = 5   s and the identification horizon 5   s are chosen as for the parameter-specific identification method. To derive the independent equations for the linear system (15), the identification horizon has to be shifted. The shifting time is an additional parameter for the identification which is set to 1   s in this example. Using the same input u ˜ o and output signal y ˜ o illustrated in Figure 1, the parameters are identified to a ^ 0 = 1.903 and b ^ 0 = 0.969 . Regarding the number of shifts, each shift leads to a new independent equation. The linear system (15) consits of n + m + 1 equations which means that at least n + m shifts are necessary to state the linear system (15). In this example, one shift is sufficient because of n = 1 and m = 0 (see (65)). The shifting of the identification horizon makes longer measurements necessary. Considering fractional systems with more parameters, the extension of the measurement duration can be significant. Because of the extension of the measurement duration, more data are considered for parameter identification. Nonetheless, the error of the identified parameters is with 4–5% for the method described in [12] significant greater than the error (approx. 0.3%) made with the parameter-specific approach which is described in this paper.

7. Conclusions

The main contribution of this paper covers the automatic generation of modulating functions by transferring the parameter identification into a control problem, as is provided in Definition (6). For this purpose, a model-based auxiliary system, which connects the fractional-order model, the modulating function as well as the boundary terms resulting from the application of the modulating function method, is defined in this paper. Instead of adjusting a generic modulating function for each parameter individually, it is sufficient to calculate an input signal which steers the model-based auxiliary system from the uninitialized state into the final state (49) corresponding to the parameter to be identified. If the input signal for the steering process is designed as an energy-optimal control, the upper bound of the identification error is minimized in the presence of noisy observations of the input and output signal.
Regarding the practical applicability, the controllability of the model-based auxiliary system is also investigated in this paper and depends on the orders of the fractional-order model and the chosen input signal u ˜ e of system (9).

Author Contributions

Conceptualization, O.S.; methodology, M.E. and O.S.; software, O.S.; validation, M.E. and O.S.; investigation, M.E. and O.S.; resources, S.H.; writing—original draft preparation, O.S.; writing—review and editing, A.J.M. and S.H.; visualization, O.S.; supervision, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge support by the KIT-Publication Fund of the Karlsruhe Institute of Technology.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CCaputo
GLGrünwald-Letnikov
RLRiemann-Liouville

Appendix A. Right-Sided Definitions of Fundamentals

In the following, the right-sided definitions of the initialized integration and derivatives of Section 2.1 are given using the notation of [19].
Definition A1. (Right-Sided Fractional Operators.)
  • Initialized Right-Sided Fractional Integral
    I t 2 α t τ h ( τ ) : = i t 2 α t τ h ( τ ) + η rs h , α , t 2 , t e , t
  • Uninitialized Right-Sided Riemann-Liouville Fractional Derivative
    d R L t e α t τ h ( τ ) : = 1 α d d t α i t e α α t t h ( t )
  • Uninitialized Right-Sided Caputo Fractional Derivative
    d C t e α t τ h ( τ ) : = 1 α i t e α α t t d d t α h ( t )
  • Initialized Right-Sided Fractional Derivative
    D t e α t τ h ( τ ) : = d t 2 α t τ h ( τ ) + ψ rs h , α , t 2 , t e , t
where Γ is the Gamma function (see [22] (pp. 1–6)).

Appendix B. Error Minimization with Energy-Optimal Control

In this section, the proof of Lemma 8 is given. Lemma 8 states that a parameter-specific input signal u s * calculated as an energy-optimal control leads to a minimal upper bound of the identification error (59) which occurs due to the presence of a noisy observation for the input and output. With the definition of the R-Matrix (see [28]), the prerequisite is given to define the energy-optimal control, which was originally introduced in [28].
Definition A2. (Energy-Optimal Control.)
Suppose a fractional state space as given in [28], l N > 0 , K : = k N > 0 k l , α R l × 1 , α k ( 0 , 1 ] and ς = min α . The energy-optimal control is defined by
u t = t 2 t 2 2 ς B t R α , rs t , t 2 , A K 1 t 1 , t 2 · x t 2 + t 1 t 2 R α τ , t , A ψ x , α , t 0 , t 1 , τ d τ
where
K t 1 , t 2 = t 1 t 2 t 2 τ 2 2 ς R α τ , t 2 , A B τ · B τ R α , rs τ , t 2 , A d τ ,
k K , R α is the left-sided R-Matrix, and R α , rs the right-sided R-Matrix of the fractional state space (see [28]). Then, (A5) steers the fractional state space from an initialization function vector ψ to any given final state vector x t 2 and minimizes the specific performance index regarding the control energy
J u , t 1 , t 2 = 1 2 t 1 t 2 t 2 τ 2 ς 2 u τ u τ d τ .
Lemma 8 states that an upper bound for the identification error exists and that a parameter-specific control input u s * calculated as an energy-optimal control (A5) minimizes this upper bound. This is shown in the following proof of Lemma 8.
Proof. 
Starting from the identification error (59) and inserting the parameter-specific control input (55) yields
Δ p s = t 1 t 2 y ˜ e t d R L t n α t 1 t φ t d t d t m e t 2 p .
To calculate an upper bound of (A8), the triangle inequality is applied
| Δ p s | | t 1 t 2 y ˜ e t d R L t n α t 1 t φ t d t | + | m e t 2 p | .
Further, applying the Cauchy-Schwarz inequality results in
| Δ p s | i = 1 n t 1 t 2 | y ˜ e t | 2 t 1 t 2 | d R L t ( i 1 ) α t 1 t φ s t | 2 · | p i | + j = 0 m t 1 t 2 | u ˜ e t | 2 t 1 t 2 | d R L t j α t 1 t φ s t | 2 · | p n + j + 1 | + t 1 t 2 | y ˜ e t | 2 t 1 t 2 | u s * t | 2 .
Next, an upper bound is calculated for the integrals which depends on the parameter-specific modulating function in (A10). First, the fractional derivative is separated into a fractional integration of order k α and a fractional derivative of order n α
t 1 t 2 | d R L t ( n k ) α t 1 t φ s t | 2 d t = t 1 t 2 | i t k α t 1 t d R L t n α t 1 t φ s t | 2 d t .
Second, the absolute value is taken under the fractional integral and the order of integration is shifted
t 1 t 2 | d R L t ( n k ) α t 1 t φ s t | 2 d t i t k α t 1 t t 1 t 2 | u s * t | 2 d t .
Third, the integrals are evaluated for α = 1 , considering the mean value theorem [30]
t 1 t 2 | d R L t ( n k ) α t 1 t φ s t | 2 d t ( t 2 t 1 ) 2 k Γ ( k ) ( 2 k 1 ) 2 k t 1 t 2 | u s * t | 2 d t .
Inserting this into (A10) yields
| Δ p s | P t 1 , t 2 , p , y ˜ e , u ˜ e t 1 t 2 | u s * t | 2
where P t 1 , t 2 , p , y ˜ e , u ˜ e is given in (64).
The integral of the parameter-specific control input u s * is rewritten using fractional integrals
t 1 t 2 | u s * t | 2 = 2 Γ ( 2 ς ) i t 2 2 2 ς t 1 t i t 2 2 ς t 1 t Γ ( 2 ς ) 2 | u s * t | 2 ( t 2 t )
= 2 Γ ( 2 ς ) i t 2 2 2 ς t 1 t J ( u s * , t 1 , t )
where ς = min α w , γ r . Because the integrand of J ( u s * , t 1 , t ) is positive on [ t 1 , t 2 ] and, hence, J ( u s * , t 1 , t ) increases with increasing time, the maximum is reached at t = t 2 , i.e.,
t 1 t 2 | u s * t | 2 2 Γ ( 2 ς ) i t 2 2 2 ς t 1 t J ( u s * , t 1 , t 2 )
and, therefore, the identification error is bounded by
| Δ p s | P t 1 , t 2 , p , y ˜ e , u s * · 2 Γ ( 2 ς ) i t 2 2 2 ς t 1 t J ( u s * , t 1 , t 2 ) .
Comparing J u s * , t 1 , t 2 in (A18) and (A7) completes the proof. □

Appendix C. Controllability of the Model-Based Auxiliary System

Because the determination of a parameter-specific modulating function is transferred into a control problem of a fractional state space, the controllability of the model-based auxiliary system has to be analyzed. For this, the Gramian, which in turn is based on the R-Matrix, (see [28]) is needed. Hence, the R-Matrix of the model-based auxiliary system is stated first and the corresponding Gramian is given thereafter. Finally, using the calculated Gramian, the controllability of the model-based auxiliary system is determined. Throughout this section, suppose w : = min k N > 0 α k 1 , C : = k = 1 n n k α , and r : = min k N > 0 max γ k 1 , where γ is defined in (25).
First, the R-Matrix of the model-based auxiliary system is evaluated.
Lemma A1. (R-Matrix of the Model-Based Auxiliary System.)
Suppose κ : = n + m + 1 , z = r · C , ϑ : = κ + n w + n α + z , P y ˜ : R R n × n w , P u ˜ : R R ( m + 1 ) × n w , P : R R n w × n w , P Δ : R R n α × n α , P : R R z × z , and R α : R R ϑ × ϑ . Then the R-Matrix of the model-based auxiliary system with the system matrix (45) results in
R α t , t 2 , A = I n × n 0 n × ( m + 1 ) P y ˜ t 0 n × n α 0 n × z 0 ( m + 1 ) × n I ( m + 1 ) × ( m + 1 ) P u ˜ t 0 ( m + 1 ) × n α 0 ( m + 1 ) × z 0 n w × n 0 n w × ( m + 1 ) P t 0 n w × n α 0 n w × z 0 n α × n 0 n α × ( m + 1 ) 0 n α × n w P Δ t 0 n α × z 0 z × n 0 z × ( m + 1 ) 0 z × n w 0 z × n α P t
where
P y ˜ t = H n × n w k = 0 n w 1 i t 2 k + 1 α w t t y ˜ o t J n w × n w k ,
P u ˜ t = H ( m + 1 ) × n w k = 0 n w 1 i t 2 k + 1 α w t t u ˜ o t J n w × n w k ,
P t = k = 0 n w 1 t 2 t k + 1 α w 1 Γ k + 1 α w J n w × n w k ,
P Δ t = k = 0 n α 1 t 2 t k J n α × n α k ,
P τ = P 1 t 0 r × r 0 r × r P C t , P j t = k = 0 r 1 t 2 t k + 1 γ j r 1 Γ k + 1 γ j r J r × r k
and J represents a Jordan matrix with the corresponding dimensions and which only has eigenvalues of 0.
Proof. 
Because the system matrix (45) is nilpotent, the fractional Peano-Baker series (see [28]) and, hence, the sum in the R-Matrix is finite. The matrices P , P Δ and P can be derived separately because of the structure of (45). For this purpose, each matrix of the subsystems is inserted into the fractional Peano-Baker series and the sum is evaluated. For subsystem □,
F α w t = k = 0 n w 1 t 2 t k α w Γ k α w + 1 J n w × n w k
results. Using an interim result of the composition with the fractional derivatives in [22] (pp. 59–60)
d R L t 1 α t τ τ t β = Γ β + 1 Γ β + α τ t β 1 + α
leads directly to the matrices P , P Δ and P .
To calculate P y ˜ and P u ˜ , the results of subsystem □ (A25) have to be taken into account which leads to (A27)–(A29) for the output signal:
P y ˜ t = H n × n w k = 0 n w 1 t t 2 y ˜ τ τ t k α w Γ k α w + 1 d τ J n w × n w k
= H n × n w k = 0 n w 1 1 Γ k α w + 1 t t 2 y ˜ τ τ t 1 k α w + 1 d τ J n w × n w k
= H n × n w k = 0 n w 1 i t k α w + 1 t t y ˜ t J n w × n w k
where H is given in (35).
Determining P u ˜ is equivalent to the steps before. □
Lemma A2. (Gramian of the Model-Based Auxiliary System.)
The Gramian of the model-based auxiliary system with the input vector (46) and the R-Matrix (A19) results in
G t 1 , t 2 = t 1 t 2 R α τ , t 2 , A b b R α τ , t 2 , A d τ
= t 1 t 2 ϱ τ , t 2 , A ϱ τ , t 2 , A d τ
where ϱ is given by
ϱ t , t 2 , A = [ i t 2 n α t t y ˜ o t , i t 2 ( n 1 ) α t t y ˜ o t , , i t 2 α t t y ˜ o t , i t 2 n α t t u ˜ o t , i t 2 ( n 1 ) α t t u ˜ o t , , i t 2 ( n m ) α t t u ˜ o t , t 2 t n α 1 Γ n α , t 2 t ( n w 1 ) α w 1 Γ ( n w 1 ) α w , , t 2 t α w 1 Γ α w , t 2 t n α 1 Γ n α , t 2 t n α 2 Γ n α 1 , , 1 , t 2 t γ 1 1 Γ γ 1 , t 2 t ( r 1 ) γ 1 r 1 Γ ( r 1 ) γ 1 r , , t 2 t γ 1 r 1 Γ γ 1 r , , t 2 t γ C 1 Γ γ C , t 2 t ( r 1 ) γ C r 1 Γ ( r 1 ) γ C r , , t 2 t γ C r 1 Γ γ C r ]
and
ϱ t , t 2 , A = ϱ t , t 2 , A .
Proof. 
The proof is given in [28]. □
Finally, the controllability of the model-based auxiliary system can be analyzed. Regarding the input signal u ˜ o we assume a non-vanishing input signal.
Assumption A3. (Non-Vanishing Input Signal.) 
Unless every coefficient b j R vanishes, we assume
j = 0 m b j D C t 2 j α t t u ˜ o t 0 .
To ensure that the model-based auxiliary system is fully controllable, two requirements have to be fulfilled. The first one is related to the fractional orders and the second one concerns the input signal.
Definition A3. (Sets of Fractional Orders.)
Suppose J : = i N 0 i n 1 , K 1 : = k 1 N > 0 k 1 n w ,
K 2 : = k 2 N > 0 k 2 n α , K 3 : = k 3 N > 0 k 3 r and
K 4 , i : = k 4 N 0 k 4 k α . The fractional orders of subsystem □ are described by
M 1 : = m 1 R k 1 K 1 : m 1 = k 1 α w 1 ,
of subsystem Δ by
M 2 : = m 2 N 0 k 2 K 2 : m 2 = n α k 2 ,
of subsystem ⋄ by
M 3 : = M 3 , 0 M 3 , n 1 ,
M 3 , i : = m 3 R k 3 K 3 k 4 K 4 , i : m 3 = k 3 i α k 4 + ( n i ) α r 1 .
To analyze if all fractional orders are different and if the input signal excite the system persistently, a special function is defined.
Definition A4. (Q-Function.)
Q t : = k = 1 n w q n w k + 1 t 2 t k α w 1 Γ k α w + k = 0 n α q n w + κ + n α k + 2 t 2 t k + i = 0 n 1 k = 0 j α l = 1 r q i , k , l t 2 t l i α k + ( n i ) α r 1 Γ l i α k + ( n i ) α r .
Lemma A3. (Controllability of the Model-Based Auxiliary System.)
The model-based auxiliary system is completely controllable if the intersection of the sets M 1 , M 2 and M 3 depending on fractional order α and the system order n is empty
M 1 M 2 M 3 =
and if the input signal u ˜ t of the fractional-order model (9) is chosen, such that
d t 2 n w 1 w α + 1 t t Q t + j = 0 m q n w + n + j + 1 d t 2 j α t t u ˜ t 0
is only fulfilled for the trivial solution. In (A41), the fractional derivative operator is equilvalent to the definition of the fractional derivatives of the investigated system (see (5), (9) and Rem. 3).
Proof. 
If the Gramian (A30) of the model-based auxiliary system is regular or the time functions of ϱ t , t 2 , A are linearly independent, the model-based auxiliary system will be complete controllable (see [31]).
Thus, it has to be shown that ϱ t , t 2 , A q 0 holds true only for q = 0 . In the following, q 0 is assumed. Thereby, the entries in (A32) without input or output signals are equivalent to the first term of the sum in (A41). Because they are of polynomial nature, single terms can eliminate each other if and only if the exponents are equal. When considering non-vanishing coefficients, Q t 0 can never occur if M 1 M 2 M 3 = holds true.
Assuming that all elements of q k = 0 for the polynomial type elements of (A32), the second term of the sum in (A41) can be transformed into the form (9) by a right-sided differentiation of order n α
q n w + 1 y ˜ t q n w + n D C t 2 ( n 1 ) α t t y ˜ t + q n w + n + 1 u ˜ t + + q n w + κ D C t 2 m α t t u ˜ t 0
which is a FDE of order ( n 1 ) α . It is assumed that (9) is explicitly given and y ˜ t and u ˜ t can only solve FDEs of order n α . Thus, y ˜ t can only vanish if u ˜ t is a homogeneous FDE
q n w + n + 1 u ˜ t + + q n w + κ D C t 2 m α t t u ˜ t 0
which is excluded by Assumption A3.
Therefore, if the input signal is chosen according to (A41) and the intersection of the sets (A35)–(A37) depending on fractional order α and system order n is empty, the model-based auxiliary system is fully controllable. □
Remark A1. 
For the integer order case with α = 1 , the sets M 2 and M 3 belonging to the subsystem Δ and subsystem are empty, because the model-based auxiliary system only consists of the subsystem and subsystem □. Thus, the controllability of the model-based auxiliary system only depends on the input signal u ˜ t according to (A41). This result is in line with the results in [16].

References

  1. Cesarano, C. Generalized special functions in the description of fractional diffusive equations. Commun. Appl. Ind. Math. 2019, 10, 31–40. [Google Scholar] [CrossRef] [Green Version]
  2. Bin-Mohsin, B.; Rafique, S.; Cesarano, C.; Javed, M.Z.; Awan, M.U.; Kashuri, A.; Noor, M.A. Some General Fractional Integral Inequalities Involving LR-Bi-Convex Fuzzy Interval-Valued Functions. Fractal Fract. 2022, 6, 565. [Google Scholar] [CrossRef]
  3. Ospina Agudelo, B.; Zamboni, W.; Monmasson, E. A Comparison of Time-Domain Implementation Methods for Fractional-Order Battery Impedance Models. Energies 2021, 14, 4415. [Google Scholar] [CrossRef]
  4. Singh, J. Analysis of fractional blood alcohol model with composite fractional derivative. Chaos Solitons Fractals 2020, 140, 110127. [Google Scholar] [CrossRef]
  5. Kanth, A.S.V.R.; Garg, N. Computational Simulations for Solving a Class of Fractional Models via Caputo-Fabrizio Fractional Derivative. Procedia Comput. Sci. 2017, 125, 476–482. [Google Scholar] [CrossRef]
  6. Allafi, W.; Zajic, I.; Burnham, K.J. Identification of Fractional Order Models: Application to 1D Solid Diffusion System Model of Lithium Ion Cell. In Progress in Systems Engineering; Selvaraj, H., Zydek, D., Chmaj, G., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 63–68. [Google Scholar]
  7. Aldoghaither, A.; Liu, D.Y.; Laleg-Kirati, T.M. Modulating Functions Based Algorithm for the Estimation of the Coefficients and Differentiation Order for a Space-Fractional Advection-Dispersion Equation. SIAM J. Sci. Comput. 2015, 37, A2813–A2839. [Google Scholar] [CrossRef] [Green Version]
  8. Dai, Y.; Wei, Y.; Hu, Y.; Wang, Y. Modulating Function-Based Identification for Fractional Order Systems. Neurocomputing 2016, 173, 1959–1966. [Google Scholar] [CrossRef]
  9. Eckert, M.; Kupper, M.; Hohmann, S. Functional Fractional Calculus for System Identification of Battery Cells. At–Automatisierungstechnik 2014, 62, 272–281. [Google Scholar] [CrossRef]
  10. Gao, Z. Modulating Function-Based System Identification for a Fractional-Order System with a Time Delay Involving Measurement Noise Using Least-Squares Method. Int. J. Syst. Sci. 2017, 48, 1460–1471. [Google Scholar] [CrossRef]
  11. Liu, D.Y.; Laleg-Kirati, T.; Gibaru, O.; Perruquetti, W. Identification of Fractional Order Systems Using Modulating Functions Method. In Proceedings of the American Control Conference (ACC), Washington, DC, USA, 7–19 June 2013; pp. 1679–1684. [Google Scholar]
  12. Stark, O.; Kupper, M.; Krebs, S.; Hohmann, S. Online Parameter Identification of a Fractional Order Model. In Proceedings of the 57th IEEE Conference on Decision and Control (CDC), Miami, FL, USA, 17–19 December 2018; pp. 2303–2309. [Google Scholar]
  13. Lu, Y.; Zhang, J.; Tang, Y.G. Parameter Identification of Fractional Order Systems Using a Collocation Method Based on Hybrid Functions. J. Dyn. Syst. Meas. Control 2020, 142, 081007. [Google Scholar] [CrossRef]
  14. Zhang, B.; Tang, Y.; Zhang, X.; Zhang, C. Parameter Identification of Fractional Order Systems Using a Hybrid of Bernoulli Polynomials and Block Pulse Functions. IEEE Access 2021, 9, 40178–40186. [Google Scholar] [CrossRef]
  15. Janiczek, T. Generalization of the Modulating Functions Method into the Fractional Differential Equations. Bull. Pol. Acad. Sci. Tech. Sci. 2010, 58, 593–599. [Google Scholar] [CrossRef]
  16. Schmid, C.; Roppenecker, G. Parameteridentifikation für LTI-Systeme mit Hilfe signalmodellgenerierter Modulationsfunktionen. At–Automatisierungstechnik 2011, 59, 521–528. [Google Scholar] [CrossRef]
  17. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  18. Lorenzo, C.F.; Hartley, T.T. Initialized Fractional Calculus; Technical Report; NASA Glenn Research Center: Cleveland, OH, USA, 2000. [Google Scholar]
  19. Eckert, M.; Kölsch, L.; Hohmann, S. Fractional Algebraic Identification of the Distribution of Relaxation Times of Battery Cells. In Proceedings of the 54th IEEE Conference on Decision and Control (CDC), Osaka, Japan, 15–18 December 2015; pp. 2101–2108. [Google Scholar]
  20. Almeida, R.; Malinowska, A.B.; Torres, D.F.M. A Fractional Calculus of Variations for Multiple Integrals with Application to Vibrating String. J. Math. Phys. 2010, 51, 033503-1–033503-12. [Google Scholar] [CrossRef] [Green Version]
  21. Lorenzo, C.F.; Hartley, T.T. Initialization, Conceptualization, and Application in the Generalized (Fractional) Calculus. Crit. Rev. Biomed. Eng. 1998, 35, 447–553. [Google Scholar] [CrossRef]
  22. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, UK, 1998. [Google Scholar]
  23. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications, 1st ed.; CRC Press: Basel, Switzerland; Philadelphia, PA, USA, 1993. [Google Scholar]
  24. Isermann, R.; Münchhof, M. Identification of Dynamic Systems: An Introduction with Applications; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  25. Stark, O.; Karg, P.; Hohmann, S. Iterative Method for Online Fractional Order and Parameter Identification. In Proceedings of the 2020 59th IEEE Conference on Decision and Control (CDC), Jeju, Korea, 14–18 December 2020; pp. 5159–5166. [Google Scholar]
  26. Victor, S.; Malti, R.; Garnier, H.; Oustaloup, A. Parameter and Differentiation Order Estimation in Fractional Models. Automatica 2013, 49, 926–935. [Google Scholar] [CrossRef]
  27. Shinbrot, M. On the Analysis of Linear and Nonlinear Dynamical Systems from Transient-Response Data; Technical Report; National Advisory Committee for Aeronautics: Moffett Field, CA, USA, 1954. [Google Scholar]
  28. Eckert, M.; Nagatou-Plum, K.; Rey, F.; Stark, O.; Hohmann, S. Controllability and Energy-Optimal Control of Time-Variant Fractional Systems. In Proceedings of the 57th IEEE Conference on Decision and Control (CDC), Miami, FL, USA, 17–19 December 2018; pp. 4607–4612. [Google Scholar]
  29. Podlubny, I.; Chechkin, A.V.; Skovranek, T.; Chen, Y.; Jara, B.M.V. Matrix Approach to Discrete Fractional Calculus II: Partial Fractional Differential Equations. J. Comput. Phys. 2009, 228, 3137–3153. [Google Scholar] [CrossRef]
  30. Bronštejn, I.N.; Semendjaev, K.A.; Musiol, G.; Mühlig, H. Handbook of Mathematics, 6th ed.; Springer: Berlin/Heidelberg, Germany; New York, NY, USA; Dordrecht, The Netherlands; London, UK, 2015. [Google Scholar]
  31. Kailath, T. Linear Systems; Prentice-Hall: Englewood Cliffs, NJ, USA, 1980. [Google Scholar]
Figure 1. Input and output signals for parameter identification.
Figure 1. Input and output signals for parameter identification.
Mathematics 10 04106 g001
Figure 2. Parameter-specific control input and course of unknown parameter in the noise-free case.
Figure 2. Parameter-specific control input and course of unknown parameter in the noise-free case.
Mathematics 10 04106 g002
Figure 3. Parameter-specific control input and course of unknown parameter in the case of a noisy observation of the output signal.
Figure 3. Parameter-specific control input and course of unknown parameter in the case of a noisy observation of the output signal.
Mathematics 10 04106 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

Stark, O.; Eckert, M.; Malan, A.J.; Hohmann, S. Fractional Systems’ Identification Based on Implicit Modulating Functions. Mathematics 2022, 10, 4106. https://doi.org/10.3390/math10214106

AMA Style

Stark O, Eckert M, Malan AJ, Hohmann S. Fractional Systems’ Identification Based on Implicit Modulating Functions. Mathematics. 2022; 10(21):4106. https://doi.org/10.3390/math10214106

Chicago/Turabian Style

Stark, Oliver, Marius Eckert, Albertus Johannes Malan, and Sören Hohmann. 2022. "Fractional Systems’ Identification Based on Implicit Modulating Functions" Mathematics 10, no. 21: 4106. https://doi.org/10.3390/math10214106

APA Style

Stark, O., Eckert, M., Malan, A. J., & Hohmann, S. (2022). Fractional Systems’ Identification Based on Implicit Modulating Functions. Mathematics, 10(21), 4106. https://doi.org/10.3390/math10214106

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