17.06.2014 Views

View - Romanian Reports in Physics

View - Romanian Reports in Physics

View - Romanian Reports in Physics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

908<br />

Dumitru Baleanu, Jihad H. Asad, Ivo Petras<br />

2<br />

been used widely to study many phenomena <strong>in</strong> eng<strong>in</strong>eer<strong>in</strong>g, physics, as well as <strong>in</strong><br />

other sciences [4–13]. The formulation of the fractional Euler-Lagrange problem<br />

has been recently drew the attention of many authors <strong>in</strong> their works [14–19].<br />

Particularly, f<strong>in</strong>d<strong>in</strong>g the numerical solutions of the equations <strong>in</strong>volv<strong>in</strong>g the left and<br />

right derivatives is still an open problem <strong>in</strong> the field of the fractional dynamics [see<br />

for Ref. 19 and the references there<strong>in</strong>].<br />

Numerical analysis of fractional differential equations has been used by many<br />

authors to solve wide ranges of differential equations [20–25]. In his recent work,<br />

Podlubny [26, 27] shows how we can numerically <strong>in</strong>vestigate differential equations<br />

us<strong>in</strong>g the so-called matrix approach method.<br />

Hav<strong>in</strong>g these th<strong>in</strong>gs <strong>in</strong> m<strong>in</strong>d, <strong>in</strong> this paper, we pay attention to study<br />

numerically the fractional Euler-Lagrange equation of a mechanical physical<br />

system called two electric pendulum.<br />

This work is organized as follows:<br />

In Sect. 2 we discussed briefly the basic def<strong>in</strong>itions of the fractional<br />

derivatives. In Sect. 3 we study the fractional two electric pendulum model. In<br />

Sect. 4 numerical analysis of the obta<strong>in</strong>ed Euler-Lagrange equation of our model is<br />

carried out. F<strong>in</strong>ally, we closed our paper with conclud<strong>in</strong>g remarks.<br />

2. BASIC DEFINITIONS<br />

In this section we discussed the def<strong>in</strong>itions of the fractional derivatives (left<br />

and right Riemann-Liouville fractional derivatives). These def<strong>in</strong>itions are used <strong>in</strong><br />

the Hamiltonian formulation and the solution of examples lead<strong>in</strong>g to the equations<br />

of motion of the fractional order. The left Riemann-Liouville fractional <strong>in</strong>tegral<br />

(LRLFI) is def<strong>in</strong>ed as follows [8, 9]:<br />

t<br />

α<br />

1<br />

α−1<br />

a t<br />

= −<br />

(α)<br />

a<br />

I x() t ( t τ) x(τ)dτ<br />

Γ ∫<br />

(1)<br />

The right Riemann-Liouville fractional <strong>in</strong>tegral (RRLFI) has the form<br />

b<br />

α<br />

1<br />

α−1<br />

t b<br />

= −<br />

(α)<br />

t<br />

I x() t (τ t) x(τ)dτ<br />

Γ ∫<br />

. (2)<br />

The left Riemann-Liouville fractional derivative (LRLFD) reads<br />

n x<br />

α<br />

1 ⎛ d ⎞ f (τ)<br />

aDx f x =<br />

α− n+<br />

1<br />

Γ( n−α) ⎜<br />

d x<br />

⎟<br />

( x−τ)<br />

a<br />

( ) ∫ dτ. (3)<br />

⎝ ⎠<br />

The right Riemann-Liouville fractional derivative (RRLFD) reads

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!