Exact Solution of Impulse Response to a Class of Fractional Oscillators and Its Stability

Oscillator of single-degree-freedom is a typical model in system analysis. Oscillations resulted from differential equations with fractional order attract the interests of researchers since such a type of oscillations may appear dramatic behaviors in system responses. However, a solution to the impulse response of a class of fractional oscillators studied in this paper remains unknown in the field. In this paper, we propose the solution in the closed form to the impulse response of the class of fractional oscillators. Based on it, we reveal the stability behavior of this class of fractional oscillators as follows. A fractional oscillator in this class may be strictly stable, nonstable, or marginally stable, depending on the ranges of its fractional order.


Introduction
Fractional systems gain increasing attention in applied sciences, ranging from mechanical engineering to electrical engineering, see, for example, 1-11 .Recall that stability is an essential property of systems, see for example 12 , for the stability of conventional systems of integer order and 13-15 for fractional systems.
One of the typical models used in system analysis is the oscillator of single-degreefreedom 16  where ω ω 2 0 − b 2 is angular frequency.Equation 1.3 implies a damped oscillation.In the case of b 0, that is, the zero damping, 1.2 reduces to The impulse response to the above system is which corresponds to a free oscillation.
Recently, research on fractional oscillators has attracted considerable interests, see for example, 17-22 .In the stability analysis of fractional oscillators, the authors in 17-19 studied a class of fractional oscillators expressed by They concluded that the above oscillator may be strictly stable as if it is a damped oscillator.
In this paper, we focus on another class of fractional oscillators that were first introduced by Lim and Muniandy 23 .It satisfies x t e t , β > 0. 1.7 The solutions to 1.7 for e t being a white noise in time domain and frequency domain are obtained in 21 , which are further extended to the oscillator with two fractional indexes in 22, 24, 25 .However, its stability remains an unsolved issue.This paper shows that this type of fractional oscillator is strictly stable when 0 < β < 1, nonstable when β > 1, and marginally stable for β 1.
The remainder of the paper is organized as follows.In Section 2, the impulse response of the fractional system 1.7 in the closed form is proposed.Stability analysis is given in Section 3. Discussions and conclusions are in Section 4.

Impulse Response
For t > 0 and v > 0, denote by D v t the fractional derivative of Caputo type 26 defined by where Γ is the Gamma function.For simplicity, we write 0 D v t by D v below.One can regard D 2 ω 2 0 β as a shifted fractional derivative of D 2ν .By taking the binomial expansion, one gets where β m is the Pochhammer symbol, that is, The fractional oscillator satisfying D 2 ω 2 0 β x t e t , β > 0 has the impulse response function h t , which is the solution to D 2 ω 2 0 β h t δ t .Denote the Laplace transform of h t by H s .Then, we have where L is the operator of the Laplace transform.Expanding the right side of the above using the binomial series yields Therefore, the impulse response to 1.7 is given by where L −1 is the operator of the inverse Laplace transform.
Another form of h t , which may be more convenient for the stability analysis, is expressed below.Considering the Laplace transform pair given in 27 , we have where J β−1/2 ω 0 t is the Bessel function of the first kind of order β − 1/2.Therefore, we have the impulse response given by

Stability Analysis
The above discussion allows us to obtain the following results concerning the stability of the fractional oscillator under consideration.Before we discuss these stability properties, we first recall the criteria of stability based on the principle of bounded-input and bounded-output BIBO .A system is said to be stable if which implies lim t → ∞ h t 0 and poles of L h t are located on the left-hand portion of the s-plane.A system is said to be nonstable if h t is increasing, and accordingly poles of L h t are located on the right-hand portion of the s-plane.One says that a system is neutral if poles of L h t are on the complex jω-axis 28, 29 .
Note that 27

3.2
Thus, according to 3.2 for v 1/2 and considering 2.8 for β 1, h t; β in 2.8 reduces to the simple case with impulse response corresponds to the free oscillator.That is, This leads to the following remark.
Remark 3.1.h t; β reduces to the impulse response to the ordinary oscillator for β 1. Figure 1 indicates the plot of h t; 1 for ω 0 1.The ordinary oscillator is neutral.One notes that h t; β is unbounded if β > 1.As a matter of fact, from 3.2 , we have Since β > 1 implies v > 1/2, one immediately sees that both the right side and the left one on the above expression are, respectively, unbounded as t → ∞. Figure 2 shows the oscillations for various values of β for w 0 1.Thus, we have the following remark.
The other interesting thing is that the h t; β becomes an oscillation with decreasing amplitude if 0 < β < 1.In fact, 3.5 Thus, Consequently, we have which implies that the poles of L h t; β for 0 < β < 1 are located on the left of the s-plane.Therefore, the following remark is an obvious consequence.Figure 3 gives plots of h t; β for several values of β for ω 0 1.
The stability of the fractional oscillators expressed by 1.7 is summarized in Table 1.

Discussions And Conclusions
As previously noted in 25 , a fractal time series can be considered as a solution to a fractional differential equation driven by a white noise.Thus, there may be research niche for other series, for example, those discussed in 30-43 .In this paper, we have presented two forms, that is, 2.6 and 2.8 , of the impulse response to the fractional system expressed by 1.7 .Such a type of fractional oscillators has dramatic performances in its stability.We have revealed that the system with 1.7 contains three subclasses of oscillators.Ordinary free oscillators are a special case of 1.7 for β 1.It corresponds to fractional oscillators of strictly stable for 0 < β < 1 and non stable if β > 1.Note that there is no damping term in the discussed oscillators in form.However, in the case of 0 < β < 1, they are strictly stable, as though they were ordinary oscillators equipped with a certain damping.On the other side, for β > 1, they are non stable.

Figure 1 :
Figure 1: Plot of the impulse response of the ordinary oscillator.
0 is the mass, c the damping constant, k > 0 the stiffness and e t the forcing function.Let 2b c/m ≥ 0 and ω 0 k/m.Then, we rewrite the above by