# Characteristic Roots of a Class of Fractional Oscillators

Authors: Ming Li; S. C. Lim; Carlo Cattani; Massimo Scalia
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

The fundamental theorem of algebra determines the number of characteristic roots of an ordinary differential equation of integer order. This may cease to be true for a differential equation of fractional order. The results given in this paper suggest that the number of the characteristic roots of a class of oscillators of fractional order may in general be infinitely great. Further, we infer that it may also be the case for the characteristic roots of a differential equation of fractional order greater than 1. The relationship between the range of the fractional order and the locations of characteristic roots of oscillators in the complex plane is considered.


---

## Body

## 1. Introduction

Oscillators are an essential component in devices in electron positron collider systems (see, e.g., Zhao et al. [1], Ma et al. [2], Zang et al. [3], Ding et al. [4], Marder et al. [5], Barroso [6], Miller et al. [7], and Lemke [8], just citing a few). As a matter of fact, oscillations are phenomena widely observed in sciences and engineering relating to high energy physics (see, e.g., Akhmediev et al. [9], Bachas [10], Winter et al. [11], Dodonov [12], Tan [13], Diamandis et al. [14], Greenwald et al. [15], Mathews et al. [16], Faiman [17], Cocho et al. [18], Baldiotti et al. [19], Kyu Shin [20], Kirson [21], Clement [22], Sikström et al. [23], Asghari et al. [24], Um et al. [25], Bahar and Yasuk [26], Hassanabadi et al. [27], Bhattacharya and Roy [28], and Saad et al. [29], simply mentioning a few).There are various structures of oscillators, such as Mathieu oscillator (Floris [30]), Liénard type oscillator (Yaşar [31]), relativistic oscillator (Osborne [32]), Schrödinger equation type oscillator (Cornwall and Tiktopoulos [33]), and Duffing oscillator (Baltanás et al. [34] and Erturk and Inman [35]). In fact, oscillators play a role in various fields, ranging from experimental physics to electronics engineering (see, e.g., Riley et al. [36], Soong and Grigoriu [37], Harris [38], Papoulis [39], Bendat and Piersol [40], Devasahayam [41], Karrenberg [42], Edson [43], and Balaban et al. [44]).This research is in the domain of fractional oscillators that attract increasing interests of physicists and engineers. More specifically, we aim at revealing specific properties of characteristic roots of a class of fractional oscillators. In doing so, we first consider an ordinary differential equation of ordern given by
(1)bndny(t)dtn+bn-1dn-1y(t)dtn-1+⋯+b0=x(t),
where n is a natural number and bn is any complex number. We always assume that at least one of the higher coefficients bn≠0 for n>1. The characteristic equation of (1) is given by
(2)B(α)=bnαn+bn-1αn-1+⋯+b0=0.The fundamental theorem of algebra says that the number of roots of (2) is n (G. A. Korn and T. M. Korn [45]). This theorem is stated in the domain of complex variables (Krantz [46]).Suppose that then roots of B(α) are α1,α2,…,αn. For each root r of multiplicity of m, either real or complex, we always consider r the m roots in what follows unless otherwise stated. Using the partial fraction expansion, B(α) can be expressed by
(3)B(α)=bn(α-α1)(α-α2)⋯(α-αn)=bn∏i=1n(α-αi).Now, we rewrite (3) by the following expression:
(4)B(α)={K∏j=1n/2(mjαj2+cjαj+kj),niseven,K(α-αn)∏j=1(n-1)/2(mjαj2+cjαj+kj),nisodd,
where mj, cj, kj, and K are constants. Without loss of generality, we can suppose that the only simple zero of B(α) is αn if n is odd.The factor(mjαj2+cjαj+kj) in (4) corresponds to the oscillator equation in the form
(5)mjd2y(t)dt2+cjdy(t)dt+kjy(t)=x(t).The characteristic equation of (5) is in the form
(6)(mjαj2+cjαj+kj)=0.There are two classes of fractional oscillators. One is in the form (Ryabov and Puzenko [47, Eq. (1)], Ahmad et al. [48], Radwan et al. [49], Drozdov [50, Eq. (9)], Tofighi and Pour [51], Tofighi [52, Eq. (2)], Blaszczyk et al. [53, Eq. (10)], and Narahari Achar et al. [54, 55])
(7)mjd2±εy(t)dt2±ε+cjd±εy(t)dt±ε+kjy(t)=x(t)for0≤ε<1.The other is expressed with the form (Lim et al. [56–58], Muniandy and Lim [59], Eab and Lim [60], and Li et al. [61])
(8)(mjD2y(t)+cjDy(t)+kjy(t))β=x(t)forβ>0.
This research utilizes the form of (8).According to the fundamental theorem of algebra, there are only two characteristic roots with respect to the oscillator equation (5). They are
(9)αj,12=-cj±cj2-4mjkj2mj.One might be carelessly misled to consider that there exist only two characteristic roots regarding the fractional oscillator equation (8) because
(10)(mjαj,122+cjαj,12+kj)β=0.However, we shall show that the number of the roots in the above expression dramatically differs from what in the following expression:(11)(mjαj,122+cjαj,12+kj)=0.The contributions of this paper are twofold. One is to exhibit that the number of the characteristic roots of (8) is in general infinitely great. The other is to reveal the relationship between the range of β and the locations of the characteristic roots of (8) in a complex plane. In addition, if all αj (j=1,…,n) are simple complex pair of roots, the ordinary differential equation of order n (1) and its generalization given by
(12)(bnDny(t)+bn-1Dn-1y(t)+⋯+b0)β=x(t),β>0
may be taken as the product of oscillators of integer order and fractional order in series in the wide sense for n being even, respectively.The rest of the paper is organized as follows. We shall give the results in Section2, including the proof that there are infinite characteristic roots regarding (8), and the explanation that (1) and (12) may be taken as oscillators in series in the wide sense. Discussions are given in Section 3, which is followed by Conclusions.


## 2. Results

### 2.1.  Result 1

The number of the characteristic roots of (8) may be infinitely great.Denote byC the set of complex numbers. Let z∈C. Suppose that a power function is given by
(13)w=zb=ebln(z).Then, the number of different values ofw relies on the value of b for a given z. More precisely, we express that by the following lemmas, which can be found in the literature, such as [45] or Yu [62].Lemma 1.
Ifb is a rational number expressed by the irreducible fraction l/m, where m≥1, the number of values of zb is m.Lemma 2.
Ifb is an irrational number or imaginary number, the number of values of zb is infinitely great.The general expression ofw is in the form
(14)w=zb=ebln(z)=eb[ln|z|+i(argz+2mπ)],i=-1.
Therefore, from Lemma 2, we have the theorem below.Theorem 3.
The number of the characteristic roots of the fractional oscillator (8) is infinitely great if β≠1.Proof.
Letcj = 0. Then, the characteristic roots αj,12 in (9) become the imaginary numbers expressed by
(15)αj,12=±imjkjmj.According to Lemma2, the number of the roots of either (α-αj,1)β or (α-αj,2)β is infinitely great. Thus, Theorem 3 results.


### 2.2. Result 2

Equations (1) and (12) may be taken as oscillators in series.Denote(mjαj2+cjαj+kj) in (4) by Bj(α):
(16)Bj(α)=(mjαj2+cjαj+kj).Without loss of generality,n is assumed to be even. In addition, we suppose that all αj (j=1,…,n) are simple complex pair of roots. Then, we have the theorem below.Theorem 4.
The ordinary differential equation (1) may be taken as an oscillator (i.e., product of oscillators) in the wide sense if n is even and all αj (j=1,…,n) are simple complex pair of roots. By wide sense, one means that it is a system consisting of the product of a series of conventional 2-order oscillators.Proof.
On the one hand,Bj(α) stands for the characteristic equation of the jth oscillator of order 2 since n is even and all αj (j=1,…,n) are simple complex pair of roots. On the other hand, the characteristic equation of (1) can be expressed by
(17)B(α)={K∏j=1n/2Bj(α),niseven,K(α-αn)∏j=1(n-1)/2Bj(α),nisodd.Based on the theory of filter design (Mitra and Kaiser [63]), the system (1) in the case of n being even may be expressed by Figure 1.Figure 1
Oscillators of 2-order in concatenation.Therefore, the system (1) may be expressed by the product of a series of 2-order oscillators.Denote byBβ(α) the characteristic equation of (12). Then,
(18)Bβ(α)={K∏j=1n/2Bjβ(α),niseven,K(α-αn)∏j=1(n-1)/2Bjβ(α),nisodd,
where
(19)Bjβ(α)=(mjαj2+cjαj+kj)β.Thus, the system of fractional order expressed by (12) may be the product of a series of fractional oscillators of (8).


## 2.1.  Result 1

The number of the characteristic roots of (8) may be infinitely great.Denote byC the set of complex numbers. Let z∈C. Suppose that a power function is given by
(13)w=zb=ebln(z).Then, the number of different values ofw relies on the value of b for a given z. More precisely, we express that by the following lemmas, which can be found in the literature, such as [45] or Yu [62].Lemma 1.
Ifb is a rational number expressed by the irreducible fraction l/m, where m≥1, the number of values of zb is m.Lemma 2.
Ifb is an irrational number or imaginary number, the number of values of zb is infinitely great.The general expression ofw is in the form
(14)w=zb=ebln(z)=eb[ln|z|+i(argz+2mπ)],i=-1.
Therefore, from Lemma 2, we have the theorem below.Theorem 3.
The number of the characteristic roots of the fractional oscillator (8) is infinitely great if β≠1.Proof.
Letcj = 0. Then, the characteristic roots αj,12 in (9) become the imaginary numbers expressed by
(15)αj,12=±imjkjmj.According to Lemma2, the number of the roots of either (α-αj,1)β or (α-αj,2)β is infinitely great. Thus, Theorem 3 results.


## 2.2. Result 2

Equations (1) and (12) may be taken as oscillators in series.Denote(mjαj2+cjαj+kj) in (4) by Bj(α):
(16)Bj(α)=(mjαj2+cjαj+kj).Without loss of generality,n is assumed to be even. In addition, we suppose that all αj (j=1,…,n) are simple complex pair of roots. Then, we have the theorem below.Theorem 4.
The ordinary differential equation (1) may be taken as an oscillator (i.e., product of oscillators) in the wide sense if n is even and all αj (j=1,…,n) are simple complex pair of roots. By wide sense, one means that it is a system consisting of the product of a series of conventional 2-order oscillators.Proof.
On the one hand,Bj(α) stands for the characteristic equation of the jth oscillator of order 2 since n is even and all αj (j=1,…,n) are simple complex pair of roots. On the other hand, the characteristic equation of (1) can be expressed by
(17)B(α)={K∏j=1n/2Bj(α),niseven,K(α-αn)∏j=1(n-1)/2Bj(α),nisodd.Based on the theory of filter design (Mitra and Kaiser [63]), the system (1) in the case of n being even may be expressed by Figure 1.Figure 1
Oscillators of 2-order in concatenation.Therefore, the system (1) may be expressed by the product of a series of 2-order oscillators.Denote byBβ(α) the characteristic equation of (12). Then,
(18)Bβ(α)={K∏j=1n/2Bjβ(α),niseven,K(α-αn)∏j=1(n-1)/2Bjβ(α),nisodd,
where
(19)Bjβ(α)=(mjαj2+cjαj+kj)β.Thus, the system of fractional order expressed by (12) may be the product of a series of fractional oscillators of (8).


## 3. Discussions

The previous section says that there are infinite roots inBjβ(α) if β≠1. In the case of β=1, Bjβ(α) reduces to the characteristic equation of the conventional oscillator (5) with two roots only. Thus, the fraction β≠1 dramatically alters the behavior of characteristic roots of oscillators. For facilitating our discussions, we omit the subscript j in what follows if not confused. More precisely, we specifically consider the fractional oscillator in the form
(20)(mD2y(t)+cDy(t)+ky(t))β=x(t)forβ>0.Figure2 shows an RLC resonance circuit in series, where R, L, and C represent resistor, inductor, and capacity, respectively. In Figure 2, I(t) is the electronic current and v(t) the power source. According to the Kirchhoff voltage law, one has
(21)D2I(t)+RLDI(t)+1LCI(t)=1Ldv(t)dt.Figure 2
Illustration of RLC resonance circuit in series.Letω=1/LC and R/L=2b. Denote (1/L)(dv(t)/dt) by e(t). Then, (21) becomes the form
(22)D2I(t)+2bDI(t)+ω2I(t)=e(t).Generalizing (22) to the fractional order β yields
(23)(D2I(t)+2bDI(t)+ω2I(t))β=e(t).Below, we specifically study the circuit in Figure2 with R=0, as indicated in Figure 3.Figure 3
Illustration of LC resonance circuit in series with dampingR=0.In the case of Figure3, (23) becomes the form
(24)(D2I(t)+ω2I(t))β=e(t).Denote byh(t) the impulse response function of (24). Then, using the techniques in fractional calculus and differential equations [64–81], we have (see [61] for details)
(25)h(t)=πΓ(β)(2ω)β-1/2tβ-1/2Jβ-1/2(ωt),β>0,t≥0,
where Jβ-1/2(ωt) is the Bessel function of the first kind of order β-1/2.The following theorems reflect the particularity of roots ofBβ(α).Theorem 5.
If0<β<1, all roots of Bβ(α) are located in the left side of the complex plane.Proof.
Note that(26)Jv(t)=∑m=0∞(-1)mm!Γ(m+v+1)(t2)2m+v.From the above, we have the asymptotic expression in the form(27)Jv(t)~1tfort⟶∞.Applying (27) to (25) produces
(28)limt→∞h(t)=limt→∞πΓ(β)(2ω)β-1/2tβ-1/2Jβ-1/2(ωt)=0,0<β<1.Denote byH(s) the Laplace transform of h(t). Then, according to the final-value theorem, we have
(29)lims→0sH(s)=0,0<β<1.The above implies that all poles ofH(s) except the origin are strictly in the left side of s plane. In the right of the s plane, H(s) is analytic. This completes the proof.Theorem 6.
Ifβ>1, at least, parts of roots of Bβ(α) are located in the right side of the complex plane.Proof.
Note that(30)Jv(t)=(t/2)vΓ(v+1/2)Γ(1/2)×∫-11(1-u2)v-1/2cos(tu)du,Rev>-12.
From the above, we have the following:
(31)-(t/2)vΓ(v+1/2)Γ(1/2)∫-11(1-u2)v-1/2du≤Jv(t)≤(t/2)vΓ(v+1/2)Γ(1/2)∫-11(1-u2)v-1/2du.Sinceβ>1 implies v>1/2, we immediately see that both the right side and the left one on the above expression are respectively unbounded when t→∞. Thus, for β>1, the fractional oscillator (24) is nonstable according to the theory of systems (Gabel and Roberts [77], Dorf and Bishop [78]). Consequently, at least, some of poles of H(s) are in the right of the s plane. Therefore, at least, parts of roots of Bβ(α) are located in the right side of the complex plane.Most of previous discussions take oscillators of fractional order (24) as a specific object. Note that the number of the characteristic roots of differential equation in general in the form of (12) may also be infinitely great. Hence, comes the following theorem in passing.Theorem 7.
Fractional-order differential equation expressed by (12) has infinite characteristic roots if n>1 and if there is at least a pair of roots that are simple complex.Proof.
The characteristic equation of (12) may be decomposed in the form of (18) due to n>1. Because there is at least a pair of roots that are simple complex, the number of the characteristic roots of (19) is infinitely great. Thus, the number of characteristic roots of (12) is infinitely great. This completes the proof.The previous discussions exhibit interesting phenomena of the characteristic roots of the oscillators of the fractional type of (24). In the future, we will work on exploring the answers of the questions described below.(i)
Are all poles ofH(s) with respect to (24) in the right of the s plane when β>1?(ii)
Might there be interesting oscillation behavior of (12) if all cj=0 in (18) and if n is even?


## 4. Conclusions

We have explained that the number of the characteristic roots of fractional-order oscillators of (24) is usually infinitely great. This conclusion has been further inferred to the case of fractional-order differential equation of (12). We have exhibited that all characteristic roots of (24) are strictly located in the left side of the complex plane if 0<β<1 and at least some of characteristic roots of (24) are in the right side of the complex plane if β>1. In the case of β=1, (24) reduce to an ordinary damping-free oscillator.



---

*Source: 853925-2013-10-08.xml*