Abstract
In this paper, a switching strategy for recursive fractional variable-order derivative is proposed. This strategy can be interpreted as an explanation of order switching mechanism for this particular type of derivative. Additionally, important properties of variable fractional order derivatives, required for prove the main result, are introduced both in a difference equation and a matrix form. Duality between the recursive and standard variable-order derivative is detailed derived. Based on the switching scheme, an analog realization of the recursive variable-order derivative definition is presented. Experimental results obtained for the analog realization are compared to the numerical results.
Similar content being viewed by others
1 Introduction to Fractional Calculus
Fractional calculus is a generalization of traditional integer order integration and differentiation actions onto non-integer order. The idea of such a generalization has been mentioned in 1695 by Leibniz and L’Hospital. In the end of 19th century, Liouville and Riemann introduced first definition of fractional derivative. However, only just in late 60’ of the 20th century, this idea drew attention of engineers. Fractional calculus was found a very useful tool for modeling the behavior of many materials and systems, especially those based on the diffusion processes. Ones of such devices that can be modeled more efficiently by fractional calculus are ultracapacitors. Models of these electronic storage devices, which capacity can reach even thousands of Farads, based on fractional order models, were presented in [4, 5]. Another system that can be successfully modelled using fractional calculus methods is the heat transfer process in both homogeneous and non-homogeneous media [6, 25]. Numerical realization of fractional dynamic systems is more complicated than for integer order case, that is why analog modeling methods for fractional order systems are widely used, e.g., in fractional order chaotic systems [1–3].
Theoretical background on fractional calculus can be found in [9–11, 15, 19, 22].
Recently, the case, when the order is changing in time, started to be intensively developed. The variable fractional order behavior can be met, for example, in chemistry (when the properties of the system are changing due to chemical reactions), electrochemistry, and other areas. In [20], experimental studies on an electrochemical example of physical fractional variable-order system are presented. In [18], the variable-order equations were used to describe a history of drug expression. Papers [20, 31, 32, 34] present methods of numerical realization of fractional variable-order integrators or differentiators. The fractional variable-order calculus can be also used to obtain variable-order fractional noise [21], and to obtain new control algorithms [12]. Some properties of such systems are presented in [13]. In [29], the variable-order interpretation of the analog realization of fractional orders integrators, realized as domino ladders, was presented. The applications of variable-order derivatives and integrals can be found also in signal processing [22].
The description of the variable-order systems is more complicated than constant-order systems. In the literature [8, 33], three general types of variable-order derivative definitions can be found, however, these definitions were given without interpretation and derivation. In [26], the switching scheme, numerically identical to the 2nd type of fractional derivative definition was introduced. This scheme can be used as an interpretation of this type of definition. In [27, 28], the recursive definition of variable-order derivative was introduced, and numerical results of comparison to other known definitions were presented.
In this work, the switching scheme corresponding to the recursive definition is given and proven. This switching scheme represents an interpretation of the recursive fractional variable-order derivative and exhibits the physical mechanism of order changing characteristic for this type of definition. The matrix approach to the recursive definition is also introduced, and the duality between this new definition and the already well-known 1st type definition of variable-order derivative is established. This property of duality is used in order to prove the main result of this paper, moreover, it constitutes an important result itself. By the duality, we mean the composition property of two types of derivatives, which for the variable-order of the opposite sign yields identity. Equivalence between the reductive-switching strategy and recursive variable-order derivative is tested based on both analog and numerical implementations.
The rest of the paper is organized as follows. Section 2 defines the problem to solve in this work. In Sect. 3, the fractional variable-order Grunwald–Letnikov type derivative definitions are recalled. Section 4 presents some properties of the 1st type of fractional variable-order derivative definition, which will be useful to develop the main result of the paper. In Sect. 5, a recursive (4th type) definition for variable-order differ-integral is given. Section 6 introduces a matrix form of the recursive type definition of variable-order derivative. The duality between the 1st and 4th type of fractional variable-order derivative definitions is established in Sect. 7. In Sect. 8, the main contribution of this paper, i.e., the identity of the recursive definition to the reductive-switching scheme is presented. Numerical verification of the equivalence between switching-order scheme and 4th type derivative is presented in Sect. 9. Finally, in Sect. 10, an analog realization of the fractional variable-order integral, according to the proposed type of derivative, is presented.
2 Problem Statement
We have build the electrical circuit (see Fig. 2) that realizes the reductive-switching scheme (more details in Sect. 10). The reductive-switching order occurs when the initial chain of derivatives is reduced according to changing the variable-order (more details are presented in Sect. 8). This switching scheme assumes that changing order is obtained by reducing chain of integrators from the input side. The simple reductive-switching case (with only one switch, i.e., one change of order) is presented in Fig. 1.
2.1 Observation
Let us assume the following, given by formula (1), description of the fractional variable-order derivative (more details about the origin and properties of this formula are given in Sect. 5)
The experimental results compared to numerical implementation of formula (1) are presented in Fig. 3.
Observation The formula given by (1) is able to describe with very high accuracy the reductive-switching behavior of the electrical circuit.
Problem formulation The problem that will be solved in this work, is to prove that the formula given by (1) is equivalent to the reductive-switching scheme, and to validate these phenomena using experimental electrical circuits.
In order to prove the above stated problem, we present below the required background from fractional calculus and introduce necessary properties of recursive fractional variable-order derivative.
3 Fractional Variable-Order Grunwald–Letnikov Type Derivatives
As a base of generalization of the constant fractional order \(\alpha \in \mathbb {R}\) derivative onto variable-order case, the following definition is taken into consideration:
where \(h>0\) is a step time, and \(n=\lfloor t/h \rfloor \).
The matrix form of the fractional order derivative is given as follows [16, 17]:
where
\(W(\alpha ,k)\in \mathbb {R}^{(k+1)\times (k+1)}\), \(w_{\alpha ,i}=\frac{(-1)^i\left( {\begin{array}{c}\alpha \\ i\end{array}}\right) }{h^\alpha }\), and \(h=t/k\), \(k\) is the number of samples.
For the case of order changing with time (variable-order case), three types of definitions can be found in the literature [8, 33]. The first one is obtained by replacing in (2) a constant-order \(\alpha \) by variable-order \(\alpha (t)\). In that approach, all coefficients for past samples are obtained for present value of the order, and is given as follows:
Definition 1
The 1st type of fractional variable-order derivative is defined as follows:
The second type of definition assumes that coefficients for past samples are obtained for order that was present for these samples. Identity of this definition to the particular (additive) switching scheme was presented in [26]. In this case, the definition has the following form:
Definition 2
The 2nd type of fractional variable-order derivative is defined as follows:
The third definition is less intuitive and assumes that coefficients for the newest samples are obtained, respectively, for the oldest orders. For such a case, the following definition applies:
Definition 3
The 3rd type of fractional variable-order derivative is defined as:
4 Some Properties of the 1st Type of Fractional Variable-Order Derivative
In order to proof the main result of this paper, i.e., the equivalence of the matrix form of the 4th type of variable-order fractional derivative and the reductive-switching scheme, we give the following additional results, which themselves state new contributions in this research area.
4.1 Matrix Approach of the 1st Type of Derivative
The matrix form of the 1st type of fractional-variable order derivative is given by [30]
where
4.2 Proposed Output-Switching Scheme of the 1st Type of Derivative
Let us introduce the following output-switching (o-s) scheme presented in Fig. 4 based on the chain of derivatives blocks related by the following switching rule. The switches \(S_i\), \(i=0,\ldots ,k\), take the following positions
and the output of such a structure we denote \(^\text {o-s}_0\mathrm {D}^{\alpha (t)}_tf(t)\).
Based on Fig. 4 we have the following result.
Lemma 1
The numerical description of the multiple output-switching scheme, when we switch between orders \(\alpha _0,\ldots ,\alpha _{k}\) every \(ih\) time instant, is the following:
where
and
Proof
The proof we will start at analysis signals from the input side. The output signal after the \(\alpha _k\) block has the following form:
The additional block \(\hat{\alpha }_{k-1}\) is connected in whole time despite of time \(k\), that is why this signal is represented by the matrix \(M(\hat{\alpha }_{k-1},k,k-1)\). At the output of this block, we obtain derivative of order \(\alpha _{k-1}\) until time \(k\) when the order is equal to \(\alpha _{k}\).
Repeating analogously, we get the output signal from the block of derivative \(\hat{\alpha }_1\) in the following form
Finally, we obtain the output signal from the block of derivative \(\hat{\alpha }_0\) in the following form:
Combining all this together, we get (7), completing the proof. \(\square \)
Theorem 1
Matrix approach of the 1st type derivative given by (5) is equivalent to the output-switching scheme given by (7), i.e.,
Proof
For simplicity, let us assume the case of one switch between orders, say \(\alpha _1\) and \(\alpha _2\) occurring at time \(T = \tau h\), \(\tau \in \mathbb {N}_+\), we have the following matrix form based on Lemma 1:
where
and
The matrix product
where \(A(\alpha _2)\in \mathbb {R}^{(k-\tau +1)\times \tau }\) and \(B(\alpha _2)\in \mathbb {R}^{(k-\tau +1)\times (k-\tau +1)}\) are suitable sub-matrices of \(W(\alpha _2,k)\), obviously corresponds to \(^{\mathcal {A}}W(\alpha ,k)\) given by (6) for
i.e., in (6) we have \(\alpha {(ih)} = \alpha _1\) for \(i=0,\ldots ,\tau -1\), and \(\alpha {(jh)} = \alpha _2\) for \(j=\tau ,\ldots ,k\). The prove of multiple-switching case can be obtained by simple analogy to the proof of one switching case. \(\square \)
5 Recursive Definition for Variable-Order Differ-Integral
The fractional constant-order difference definition (that can be obtained from Definition 1 for \(\alpha (t)=\mathrm {const}\)) can be rewritten as follows [27]:
where \(\Delta ^{n}_x(z)\) being Z-transform of the signal difference of order \(\alpha \) of variable \(x(t)\). It can be also rewritten as
which gives the following relation
and can be represented in the time domain as
Finally, it can be rewritten in the following form
This type of difference is obtained from all values of previous differences. For variable-order case, we can obtain the following definition:
Definition 4
The 4th type of fractional variable-order difference is defined as follows:
For a continuous time domain case, the 4th difference definition can be rewritten in the following form [27].
Definition 5
The 4th type of fractional variable-order derivative is defined as follows:
Remark 1
For a fractional constant-order \(\alpha =\mathrm {const}\), the fractional derivative given by Definition 5 is numerically identical with constant-order fractional derivative given by (2).
Numerical results of the 4th type derivatives, of variable-order \(\alpha _3(t)\) given by
compared to the other types of derivative definitions are presented in Figs. 5, 6, 7 and 8.
It can be seen from Figs. 5 and 6 that the derivative results in the case of 3rd and 4th type definition, for input signal being Heaviside step function, i.e., \(f(t) = 1(t)\), are very near to each-other. It follows from the fact that for the constant signal the sequence of samples is insignificant. However, for the linear input signal \(f(t) = t\cdot 1(t)\), the difference between these both definitions becomes clearly visible, what is seen in Figs. 7 and 8.
As it can be noticed in Figs. 7 and 8 for the linear input signal the differences between the 4th and others types of derivatives (i.e. 1st, 2nd, and 3rd), are indeed well visible and significant. More plots comparing all the definitions are presented in [27].
6 Matrix Approach for the 4th Type of Fractional Difference Definition
Let us recall the 4th type of variable-order difference definition in the following form:
for \(l=0,1,2,\ldots ,k\). This definition can be rewritten in the matrix form, what is given by the following lemma:
Lemma 2
Fractional difference of 4th type given by Definition 4 can be expressed in the following matrix form:
where
and
and for \(r = 1,\ldots ,k\)
where
and \(v_{-\alpha _r,i} = (-1)^i{-\alpha _r \atopwithdelims ()i}\), for \(i=1,\ldots ,r\), i.e., the \(j\)th element of \(q_r\) is
Proof
It is obtained after consecutive evaluating (10) for each time step \(l=0,1,\ldots ,k\). First, for \(l = 0\), we can write
Next, for \(l=1\):
for \(l=2\):
and, generally, for \(l = r\), we have
Finally, for \(l=k\):
where \(q_{k} = (-v_{-\alpha _k,k}, \ldots ,-v_{-\alpha _k,1})\). Combining all this together, completes the proof. \(\square \)
Remark 2
For step time \(h=1\), we have \(\mathfrak {Q}(\alpha _0,0) = I_{k+1,k+1}\), which implies \(\mathfrak {Q}_0^k = \mathfrak {Q}(\alpha _k,k)\cdots \mathfrak {Q}(\alpha _1,1)\).
Remark 3
From direct multiplication of (13), the following matrix product
is given by
where
Example 1
For \(m=0\) and \(n=k\), from (3), we get the following matrix
where \(\mathfrak {q}_{i,j}\) for \(1\le j < i \le k+1\) are given by (15).
Remark 4
Matrix (3) can be written as
where
where
Example 2
For \(n=m\), from (3), we get
where \(\mathfrak {q}_{m+1,i} = (q_m)_i\) for \(i=1,\ldots ,m\). Then, we obtain
which matches to \(\mathfrak {Q}(\alpha _r,r)\) given by (13) for \(m=r\).
Remark 5
The entries \(\mathfrak {q}_{i,j}\) and \(q_{i,j}\) of (3) and (18), respectively, are equal to each other, for \(i=m+1,\ldots ,n+1\) and \(j = 1,\ldots ,m\), i.e., \(\mathfrak {q}_{i,j}=q_{i,j}\). Moreover, matrix \(Q_m^n\) does not depend on \(h\).
Example 3
For \(k=4\), \(m=1\), and \(n=2\), using (14) or (3), we get
where \((q_i)_j\) stands for \(j\)th entries of vector \(q_i\). More precisely, it is given by
Remark 6
A natural extension of the above considerations (by taking the limit \(h\rightarrow 0\)), is the following form of the 4th type of variable-order derivative definition:
where \(\mathfrak {Q}_0^k\) is given by (12).
7 Duality Between the 1st and 4th Type of Fractional Variable-Order Derivatives
In this section, we present the duality between the 1st and 4th type of fractional variable-order derivatives, which will be used to derive the main result of this paper, i.e., the equivalence of the 4th definition and the multiple reductive-switching case. By the duality, we mean the composition property of two types of derivatives, which for the variable-order of the opposite signs yields identity. The duality property, constitutes also an important value itself in the fractional calculus domain.
Lemma 3
The inverse of \(\mathfrak {Q}_m^n\) is
where
or
where
Proof
We will show that the product of matrices \(\left( \mathfrak {Q}_m^n\right) ^{-1}\) and \(\mathfrak {Q}_m^n\) given by (22) and (17), respectively, is an identity, i.e., \(\left( \mathfrak {Q}_m^n\right) ^{-1}\mathfrak {Q}_m^n = I_{k+1,k+1}\). Let us calculate
where \(\bar{D}\) and \(D\) are the diagonal matrices given in (22) and (17), respectively. Then, we have to show that \(\bar{Q}_m^nQ_m^n=I_{k+1,k+1}\), because, obviously \(\bar{D}D=I_{k+1,k+1}\). Denote by \((\bar{Q}_m^n)^i\) the \(i\)th row of \(\bar{Q}_m^n\), and by \((Q_m^n)_j\) the \(j\)th column of \(Q_m^n\). Then, one has to calculate
for \(i = m+1,\ldots ,n\) and \(j = 1,\ldots ,m\). From (18) and (15), the straightforward calculation gives
Substituting (25) to (24), we get \((\bar{Q}_m^n)^i(Q_m^n)_j = 0\) for \(i = m+1,\ldots ,n\) and \(j = 1,\ldots ,m\).
Similarly, for \(i = m+1,\ldots ,n\) and \(j = m+1,\ldots ,n+1\), and \(i\ge j\), we have
Because from (26), for \(i=j\), we obtain \((\bar{Q}_m^n)^i(Q_m^n)_i = q_{i,i}\), and \(q_{i,i}=1\), then \((\bar{Q}_m^n)^i(Q_m^n)_i = 1\). Since for \(i>j\) from (18) and (15) we have
after substituting (27) into (26), we get \((\bar{Q}_m^n)^i(Q_m^n)_j = 0\). The case, for \(i = m+1,\ldots ,n\) and \(j = m+1,\ldots ,n+1\), where \(i<j\), is obvious, i.e., \((\bar{Q}_m^n)^i(Q_m^n)_j = 0\). The other subblock multiplications are evident, so the proof is finished. \(\square \)
Example 4
Using (21) and (22), the inverse of \(\mathfrak {Q}_0^k\) is
or
where \(\bar{q}_i\), \(0\le i\le k\), are given by (23).
Example 5
For \(n=m\), from (22) we get
which is indeed the inverse of (19).
Remark 7
Because the structures of the matrices \(\mathfrak {Q}_0^k\) and \(\left( \mathfrak {Q}_0^k\right) ^{-1}\) are different, in the sense that in the first case each entry of the same column, say \(i\)th, is multiplied by the same term \(h^{-\alpha _{i-1}}\), \(i=1,\ldots ,k+1\), and in the second case each entries of the same row, say \(j\)th, is multiplied by the same term \(h^{\alpha _{j-1}}\), \(j=1,\ldots ,k+1\), the following, in general, occurs
which implies
where \(\tilde{\alpha }= \{\alpha _0,\ldots ,\alpha _k\}\). It means that for a variable-order difference (derivative) the semigroup property does not hold. However, for a constant-order, i.e., \(\alpha _0=\cdots =\alpha _k\), and, respectively, \(\alpha (t) = \mathrm {const}\), this property holds.
Let us recall the definition of the fractional variable-order difference of the 1st type:
Lemma 4
The inverse matrix \(\left( \mathfrak {Q}_0^k\right) ^{-1}\) with orders \(\{\alpha _0,\ldots ,\alpha _k\}\) defines the 1st type of variable-order difference definition for orders \(\{-\alpha _0,\ldots ,-\alpha _k\}\), i.e.,
where \(y_i = ^\mathcal {D}\Delta ^{\alpha _i} x_i\), for \(i = 0,\ldots ,k\).
Proof
It is sufficient to show that \(\left( \mathfrak {Q}_0^k\right) ^{-1}y\) defines the 1st type of definition given by (29), i.e., it corresponds to the right-hand side of (29). Let us denote by \(\left( \mathfrak {Q}_0^k\right) ^{-1}_i\) the \(i\)th row of matrix \(\left( \mathfrak {Q}_0^k\right) ^{-1}\). Then, from (28), we get
In general, for any \(0\le i\le k\), we get
which shows the desired equivalence, i.e., \(\left( \mathfrak {Q}_0^k\right) ^{-1}_i y=^\mathcal {A}\Delta ^{-\alpha _i} y_i\), for \(0\le i\le k\).
The same can also be proved by comparing matrix \(\left( \mathfrak {Q}_0^k\right) ^{-1}\) with the matrix \(^\mathcal {A}W(\alpha ,k)\) defining the 1st derivative definition (see [30]). \(\square \)
Remark 8
Obviously, since (30) holds, the following is also true
where \(y(ih) = _0^{\mathcal {A}}\mathrm {D}^{\alpha (t)}_{ih} f(t)\), for \(i = 0,\ldots ,k\).
Theorem 2
The following dualities hold
where \(\tilde{\alpha }= \{\alpha _0,\ldots ,\alpha _k\}\).
Proof
In both cases, i.e., in the difference and derivative case, it simply follows from the composition of matrices \(\left( \mathfrak {Q}_0^k\right) ^{-1}\) and \(\mathfrak {Q}_0^k\), which obviously yields identity. \(\square \)
8 Main Result—Equivalence of Reductive-Switching Order Case
In this section, the equivalence between the reductive-switching scheme and the 4th type of fractional variable-order derivative definition will be presented, for clarity and simplicity, in two parts, namely, in the case of simple (one) switching, and later in the general case—multiple switching.
8.1 Simple Reductive-Switching Order Case
The reductive-switching (r-s) order case occurs when the initial chain of derivatives is reduced according to changing the variable-order.
The idea, introduced in [30], is depicted in Fig. 9, where the switches \(S_i\), \(i=1,2\), change their positions depending on an actual value of \(\alpha (t)\). If we want to switch from \(\alpha _1\) to \(\alpha _2\), then, before switching time \(T\), we have: \(S_1 = a\), \(S_2 = b\), and after this time: \(S_1 = b\) and \(S_2 = a\). At the instant time \(T\), the derivative block of complementary order \(\hat{\alpha }_{1}\) is disconnected from the front of the derivative block of order \(\alpha _2\), where
The numerical scheme for reductive-switching case can be obtained similarly to additive-switching scheme (see details in [30]).
Lemma 5
[30] For a reductive-switching order case, when the switch from order \(\alpha _1\) to order \({\alpha _{2}}\) occurs at time \(T = \tau h\), \(\tau \in \mathbb {N}_+\), the numerical scheme has the following form:
where
and
The order \({\hat{\alpha }_{1}}\), appearing above, is given by relation (31).
Theorem 3
In the case of simple switching (between two orders), matrix description (20) corresponding to the 4th type of fractional variable-oreder derivative is equivalent to the matrix form of simple reductive-switching order derivative scheme given in Lemma 5.
Proof
Let us consider the switching order given by (32). In the case of simple reductive-switching scheme based on Lemma 5, we have the following switching description
where \(A(\alpha _2)\in \mathbb {R}^{(k+1-\tau )\times \tau }\), and where we have used the already known property, resulting \(W(\alpha _2,\tau -1)W(\hat{\alpha }_1,\tau -1) = W(\alpha _2+\hat{\alpha }_1,\tau -1) = W(\alpha _1,\tau -1)\) (see, e.g., [16]).
From the other hand, for the 4th type of derivative and the same switching between orders given by (32), we have the following switching matrix form based on (3):
where \(B(\alpha _2)\in \mathbb {R}^{(k+1-\tau )\times \tau }\), \(C(\alpha _2)\in \mathbb {R}^{(k+1-\tau )\times (k+1-\tau )}\), and \(D(\alpha _1)\in \mathbb {R}^{\tau \times \tau }\) are the suitable block matrices. Assuming for a moment that if we would not change the order, i.e., \(\alpha _1 = \alpha _2 =\alpha \), we would obtain \(^\mathcal {D}P = W(\alpha ,k)\) according to Remark 1. From this, we conclude that \(D(\alpha ) = W(\alpha ,\tau -1)\) and \(C(\alpha ) = W(\alpha ,k-\tau )\), and then it gives rise to
Now, we have to show that \(B(\alpha _2)W(\alpha _1,\tau -1) = A(\alpha _2)W(\hat{\alpha }_1,\tau -1)\). Assume again for a moment that we do not switch the order, i.e., we have \(\alpha _1 = \alpha _2 =\alpha \), both in the case of 4th definition and simple-reductive case. Thus, using again Remark 1, we get
because \(\hat{\alpha }_1 = 0\), and then \(W(\hat{\alpha }_1,\tau -1) = I_{\tau ,\tau }\), from which we conclude that
and then, using (33)
Finally, we get \(^\mathcal {D}\!\!P =\, ^\text {r-s}\!P\), which ends the proof. \(\square \)
8.2 Multiple Reductive-Switching Order Case
The reductive-switching order case occurs when the initial chain of derivatives is reduced according to changing the variable-order (Fig. 10).
The switching rule in the multiple-switching case is an analogous extension of the simple-switching case described in Sect. 8.1.
Lemma 6
For a multiple reductive-switching order case, when the switch between orders \(\alpha _0,\ldots ,\alpha _{k}\) occurs every \(ih\) time instant, the numerical scheme has the following form:
where
and
The proof of Lemma 6 is a straightforward extension of the proof of Lemma 5 given in [30].
Theorem 4
Matrix description (20) corresponding to the 4th type of fractional-order derivative definition is equivalent to the matrix form of multiple reductive-switching order scheme given by Lemma 6, i.e.,
Proof
In order to proof this theorem, we will use the introduced in Sect. 7 fact about the duality between 1st and 4th type of fractional variable-order derivatives; it means that we have
where \(^\mathcal {A}W(-\alpha ,k)\) is equal to the matrix product given by (7) expressed for \(-\alpha \). Thus, instead of calculating
we perform the following composition
which yields the identity, and thereby ends the proof. \(\square \)
9 Numerical Verification of Equivalence Between Switching-Order scheme and 4th Type Derivative
Numerical verification of reductive-switching scheme is done in Simulink by comparing it to the numerical implementation of the recursive fractional variable-order derivative. It is depicted in Fig. 11, where: \(\mathrm {D}i\), \(i=1,\ldots ,8\), are fractional constant-order derivative blocks, \(\mathrm {S}i\), \(i=1,\ldots ,7\), are switching blocks, and FVOD is the numerical implementation of the recursive fractional variable-order derivative (for more details see [23]).
The effects of simulation according to Fig. 11 are depicted in Fig. 12, where the equivalence between recursive definition and reductive-switching scheme is also verified.
10 Analog Realization of Switching-Order Scheme
In this section, an analog circuit that corresponds to the proposed numerical definition of 4th type is presented.
10.1 Analog Realization of the Half-Order Integral
In this paper, the following method of half-order integrator implementation, introduced in [24], and meticulously investigated in [14, 29], will be used. The scheme of this method is presented in Fig. 13.
Based on the algorithm described in details in [24], the experimental circuit boards of half-order impedances were constructed. The circuit board that has been used in experimental setup is presented in Fig. 14. The circuit consists of \(200\) discrete elements with the following values: \(R_1=2.4\,\)k\(\mathrm {\Omega }\), \(R_2=8.2\,\)k\(\mathrm {\Omega }\), \(C_1=330\,\)nF, and \(C_2=220\,\)nF.
10.2 Analog Realization of the \(0.25\) Order Impedance
Method mentioned in Sect. 10.1 can be extended to build a fractional order integrator of order \(0.25\). This can be done by replacing the capacitors in the scheme in Fig. 13 by half-order integrators, which can be \(0.5\) order domino ladders. This gives an impedance of order \(\alpha =0.25\), which corresponds to a quarter-order integrator.
In Fig. 15, the scheme of the approximation of a quarter-order integrator is shown; \(Z_{0.5}\) is the impedance of half-order domino ladder. As a quarter-order impedance, the circuit board presented in Fig. 16 has been used. The board contains about 5,000 discrete elements and was designed according to the scheme shown in Fig. 15. The main ladder includes \(50\) branches with the following resistors’ values: \(R_1=2.4\,\)k\(\mathrm {\Omega }\), \(R_2=8.2\,\)k\(\mathrm {\Omega }\). The half-order integrators have been used in the quarter ladders denoted as \(Z_{0.5}\) on the scheme.
10.3 Experimental Setup
Analog model of switching system, used in experimental setup, corresponds to the switching scheme given in Fig. 9. It contains four TL071 operational amplifiers: two (\(A_1\) and \(A_3\)) in fractional integrator configuration, and two (\(A_2\) and \(A_4\)) in inverting amplifier configuration. All of operational amplifiers were supplied with external symmetrical voltage source \(\pm 12\,V\). The fractional integrators were realized similarly to configuration of analog realization of integer order integrator, but in place of capacitor the fractional impedances \(Z_1\), \(Z_2\) were used. The fractional order impedances were the domino ladder approximations presented in Sects. 10.1 and 10.2. In such a configuration, obtained output signal is inverted, that is why another operational amplifier in configuration of inverted amplifier is used. The operational amplifier TL071 was specially used, because it has ability to connect special circuit for voltage offset reduction; the scheme of this circuit is described in amplifier data-sheet [7], Fig. 3]. Reduction of this offset voltage is very important, because of integration character of the circuits, and integrated offset could have significant influence in simulation error. We have equipped all used operational amplifiers in Input Offset-Voltage Null Circuit and we tuned it before obtaining experimental results.
As a realization of switches \(S_1\) and \(S_2\), integrated analog switches DG303 were used. Fractional order impedances with order equal to \(-0.5\) or \(-0.25\), used in experiments, are shown in Figs. 14 and 16, respectively. The experimental circuit (for data acquisition and DG303 switches control) was connected to the dSPACE DS1104 PPC card with a PC. Resistors \(R\) allow to adjust a gain of integrators, and impedances \(Z_1\) and \(Z_2\) are used to obtain desired switching-order configuration (according to the reductive-switching scheme presented in Fig. 9). At the beginning, both integrators are connected in the series (switch \(S_1\) is in position \(a\), and switch \(S_2\) is in position \(b\)), and after order switch in time \(T\) the integrator at front of the chain (based on impedance \(Z_1\)) is disconnected (switch \(S_1\) is in position \(b\), and switch \(S_2\) is in position \(a\)).
10.4 Experimental Results
Example 6
Integrator with switched orders from \(\alpha =-0.5\) to \(\alpha =-0.25\).
In this case, according to the scheme in Fig. 17, the structure has the following parameters: \(Z_1\) and \(Z_2\) are the quarter-order integrators, \(R=100\,\)k\(\mathrm {\Omega }\). The identification results were obtained by numerical minimization of time responses square error for constant orders, time interval \(t\in [0,1.5]\), with sampling time \(T_\mathrm {s}=0.001\) s, and input signal \(u(t)=1(t)\) being the Heaviside step function. After identification process, the following models for orders \(-0.5\) and \(-0.25\), respectively, in time domain, were obtained:
which gives rise to the following variable-order integrator:
or using the duality property:
where (for the switching time \(T=0.7\) s.)
and
The experimental results compared to numerical implementation of the 4th type of variable-order derivative definition are presented in Figs. 18 and 19.
Example 7
Integrator with switched orders from \(\alpha =-1\) to \(\alpha =-0.5\).
In this configuration, according to the scheme in Fig. 17: \(Z_1\) and \(Z_2\) are the half-order integrators, \(R=100\,\)k\(\mathrm {\Omega }\). The identification results were obtained by numerical minimization of time responses square error with sampling time \(T_\mathrm {s}=0.001\) s, and input signal \(u(t)=0.01 \cdot 1(t)\). After identification process, the following models for orders \(-1\) and \(-0.5\), respectively, in time domain, were obtained:
The switching time was equal to \(0.7\) s. The experimental results are presented in Fig. 3.
Example 8
Integrator with switched orders from \(\alpha =-0.5\) to \(\alpha =-0.25\) of ramp function input signal.
In this example, configuration of the experimental setup is similar to the one already used in Example 6.
The identification results were obtained by numerical minimization of time responses square error with sampling time \(T_\mathrm {s}=0.001\) s, and input signal \(u(t)=0.1 \cdot t1(t)\). After identification process, the following models for orders \(-0.5\) and \(-0.25\), respectively, in time domain, were obtained:
The switching time was equal to \(0.7\) s. The experimental results are presented in Figs. 20 and 21.
Example 9
Integrator with switched orders from \(\alpha =-1\) to \(\alpha =-0.5\) of ramp function input signal.
In this example, configuration of the experimental setup is similar to the one already used in Example 7.
The identification results were obtained by numerical minimization of time responses square error with sampling time \(T_\mathrm {s}=0.001\) s, and input signal \(u(t)=0.1 \cdot t1(t)\). After identification process, the following models for orders \(-1\) and \(-0.5\), respectively, in time domain, were obtained:
The switching time was equal to \(0.7\) s. The experimental results are presented in Figs. 22 and 23.
11 Conclusions
The paper presented a recursive definition of variable-order differ-integral stated in the form of difference equation, and also in the derived equivalent matrix form. This new definition is also expressed in the form of fractional variable-order derivative. The main result of the paper is a derivation of the identity between the presented recursive definition and the reductive-switching scheme of variable-order derivative. Moreover, in the paper, the equivalence between the 1st type of variable-order derivative definition and the output-switching scheme is also derived and presented. The obtained results allow us to better understand the behavior and peculiarity of this type of definitions of variable-order derivative. Based on this knowledge, it can give rise to more appropriate choice of definitions type, depending on particular application. The paper presents also an analog circuit that corresponds to the proposed numerical definition. Our obtained experimental results show high accuracy of the proposed method of analog modeling. The obtained results can be used, in the future, to improve existing algorithms in control and signal processing areas.
References
D. Chen, Z. Sun, X. Ma, L. Chen, Circuit implementation and model of a new multi-scroll chaotic system. Int. J. Circuit Theory Appl. 407–424 (2012). doi:10.1002/cta.1860
D. Chen, C. Liu, C. Wu, Y. Liu, X. Ma, Y. You, A new fractional-order chaotic system and its synchronization with circuit simulation. Circuits, Syst. Signal Process. 31(5), 1599–1613 (2012). doi:10.1007/s00034-012-9408-z
D. Chen, C. Wu, H. Iu, X. Ma, Circuit simulation for synchronization of a fractional-order and integer-order chaotic system. Nonlinear Dyn. 73(3), 1671–1686 (2013)
A. Dzielinski, G. Sarwas, D. Sierociuk, Time domain validation of ultracapacitor fractional order model. In 49th IEEE Conference on Decision and Control (CDC), 2010, pp. 3730–3735. doi:10.1109/CDC.2010.5717093.
A. Dzielinski, G. Sarwas, D. Sierociuk, Comparison and validation of integer and fractional order ultracapacitor models. Adv. Differ. Equ. 2011, 11 (2011). doi:10.1186/1687-1847-2011-11
A. Dzieliski, D. Sierociuk, Fractional order model of beam heating process and its experimental verification, in New Trends in Nanotechnology and Fractional Calculus Applications, ed. by D. Baleanu, Z.B. Guvenc, J.A.T. Machado (Springer, Netherlands, 2010), pp. 287–294
Texas Instruments, TL07x Low-Noise JFET-Input Operational Amplifiers (2014). http://www.ti.com/lit/ds/symlink/tl071.pdf. Accessed 30 Sept 2014
C. Lorenzo, T. Hartley, Variable order and distributed order fractional operators. Nonlinear Dyn. 29(1–4), 57–98 (2002). doi:10.1023/A:1016586905654
K. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differenctial Equations (Wiley, New York, 1993)
C.A. Monje, Y. Chen, B.M. Vinagre, D. Xue, V. Feliu, Fractional-order Systems and Controls. Springer, Berlin (2010).
K.B. Oldham, J. Spanier, The Fractional Calculus. Academic Press, Waltham, MA (1974).
P. Ostalczyk, T. Rybicki, Variable-fractional-order dead-beat control of an electromagnetic servo. J. Vib. Control 14(9–10), 1457–1471 (2008)
P. Ostalczyk, Stability analysis of a discrete-time system with a variable-, fractional-order controller. Bull. Polish Acad. Sci. 58(4), 613–619 (2010)
I. Petras, D. Sierociuk, I. Podlubny, Identification of parameters of a half-order system. IEEE Trans. Signal Process. 60(10), 5561–5566 (2012). doi:10.1109/TSP.2012.2205920
I. Podlubny, Fractional Differential Equations. Academic Press, Waltham, MA (1999).
I. Podlubny, Matrix approach to discrete fractional calculus. Frac. Calc. Appl. Anal. 3, 359–386 (2000)
I. Podlubny, A. Chechkin, T. Skovranek, Y. Chen, B.M. Vinagre Jara, Matrix approach to discrete fractional calculus II: partial fractional differential equations. J. Comput. Phys. 228(8), 3137–3153 (2009)
L. Ramirez, C. Coimbra, On the variable order dynamics of the nonlinear wake caused by a sedimenting particle. Physica D-Nonlinear Phenom. 240(13), 1111–1118 (2011)
S. Samko, A. Kilbas, O. Maritchev, Theory and Applications, Fractional integrals and derivative (Gordon and Breach, London, 1987)
H. Sheng, H. Sun, C. Coopmans, Y. Chen, G.W. Bohannan, Physical experimental study of variable-order fractional integrator and differentiator. In Proceedings of the 4th IFAC Workshop Fractional Differentiation and its Applications FDA’10 (2010).
H. Sheng, H. Sun, Y. Chen, T. Qiu, Synthesis of multifractional gaussian noises based on variable-order fractional operators. Signal Process. 91(7), 1645–1650 (2011)
H. Sheng, Y. Chen, T. Qiu, Signal Processing Fractional Processes and Fractional-Order Signal Processing (Springer, London, 2012)
D. Sierociuk, Fractional Variable Order Derivative Simulink Toolkit (2012). http://www.mathworks.com/matlabcentral/fileexchange/38801-fractional-variable-order-derivative-simulink-toolkit. Accessed 30 Sept 2014
D. Sierociuk, A. Dzielinski, New method of fractional order integrator analog modeling for orders 0.5 and 0.25. In Proceedings of the 16th International Conference on Methods and Models in Automation and Robotics (MMAR), 2011, pp. 137–141. Miedzyzdroje, Poland (2011). doi:10.1109/MMAR.2011.6031332.
D. Sierociuk, A. Dzielinski, G. Sarwas, I. Petras, I. Podlubny, T. Skovranek, Modelling heat transfer in heterogeneous media using fractional calculus. Philos. Trans. Royal Soc. A 371(1990), 20120146 (2013).
D. Sierociuk, W. Malesza, M. Macias, Equivalent switching strategy and analog validation of the fractional variable order derivative definition. In Proceedings of European Control Conference ECC’2013, Zurich pp. 3464–3469 (2013).
D. Sierociuk, W. Malesza, M. Macias, On a new definition of fractional variable-order derivative. In Proceedings of the 14th International Carpathian Control Conference (ICCC), (Rytro, Poland 2013), pp. 340–345. doi:10.1109/CarpathianCC.2013.6560566.
D. Sierociuk, W. Malesza, M. Macias, Switching scheme, equivalence, and analog validation of the alternative fractional variable-order derivative definition. In Proceedings of the 52nd IEEE Conference on Decision and Control, Florence, Italy, December 10–13, 2013.
D. Sierociuk, I. Podlubny, I. Petras, Experimental evidence of variable-order behavior of ladders and nested ladders. IEEE Trans. Control Syst. Technol. 21(2), 459–466 (2013). doi:10.1109/TCST.2012.2185932
D. Sierociuk, M. Macias, W. Malesza, Analog modeling of fractional switched order derivative using different switching schemes. IEEE J. Emerg. Sel. Topics Circuits Syst. 3(3), 394–403 (2013)
C.C. Tseng, Design and application of variable fractional order differentiator. In Proceedings of the 2004 IEEE Asia-Pacific Conference on Circuits and Systems, vol. 1, pp. 405–408 (2004).
D. Valerio, J.S. da Costa, Variable-order fractional derivatives and their numerical approximations II-complex orders. In Proceedings of the Symposium on Fractional Signals and Systems, Caparica (2009).
D. Valerio, J.S. da Costa, Variable-order fractional derivatives and their numerical approximations. Signal Process. 91(3, SI), 470–483 (2011). doi: 10.1016/j.sigpro.2010.04.006.
D. Valerio, B. Vinagre, J. Domingues, J.S. da Costa, Variable-order fractional derivatives and their numerical approximations I-real orders. In Proceedings of the Symposium on Fractional Signals and Systems, Caparica, 2009 (2009).
Acknowledgments
This work was supported by the Polish National Science Center with the decision number DEC-2011/03/D/ST7/00260.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
About this article
Cite this article
Sierociuk, D., Malesza, W. & Macias, M. On the Recursive Fractional Variable-Order Derivative: Equivalent Switching Strategy, Duality, and Analog Modeling. Circuits Syst Signal Process 34, 1077–1113 (2015). https://doi.org/10.1007/s00034-014-9895-1
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00034-014-9895-1