V01a03 Calculo Fracionario Aplicado

Cálculo Fracionário Aplicado ao Problema da
Tautócrona∗
Pedro Felipe Pavanelo Ramos †
Rubens de Figueiredo Camargo †
Resumo
Apresentamos neste trabalho um estudo introdutório sobre integrais e derivadas
de ordens arbitrárias, o assim chamado cálculo de ordem não-inteira, popularizado
com o nome de Cálculo Fracionário. Em particular, discutimos e resolvemos o
clássico problema da tautócrona, isto é, o problema de se determinar a curva na qual
o tempo gasto por um objeto para deslizar, sem fricção, em gravidade uniforme até
seu ponto de mı́nimo é independente de seu ponto de partida, utilizando integrais e
derivadas de ordem não inteira.
Palavras Chave: Cálculo Fracionário, Integral Fracionária, Derivada Fracionária
de Riemann-Liouville.
Introdução
O cálculo fracionário tem sua origem em 30 de Setembro de 1695 em uma carta
escrita por l’Hospital ao seu amigo Leibniz [13, 18], na qual o significado de uma
derivada de ordem meio é proposto e discutido.
A resposta de Leibniz ao seu amigo, somada à contribuição de inúmeros bri-
lhantes matemáticos como Euler, Lagrange, Laplace, Fourier, Abel, Heaviside, Liou-
ville, entre outros, levaram às primeiras definições de derivadas e integrais de ordens
não-inteiras, que no final do século XIX, devido primordialmente as definições pro-
postas por Riemann-Liouville e Grünwald-Letnikov, pareciam estar completas[16].
Até próximo ao final do século passado o desenvolvimento do cálculo fracionário
deu-se estritamente no campo da matemática pura, sem grandes aplicações em ou-
tras áreas. Contudo, em 1969 M. Caputo, em seu livro Elasticità e Dissipazione [13],
resolveu problemas de viscoelasticidade utilizando uma nova definição, proposta por
ele, para a derivada de ordem fracionária. Além disso, o autor utilizou sua definição
para descrever problemas de sismologia. Por outro lado, a assim chamada deri-
vada fracionária de Grünwald-Letnikov, mostrou-se bastante eficiente para resolver
problemas numéricos.
A partir das definições introduzidas por Caputo e Grünwald-Letnikov, que têm
um caráter local, nas últimas décadas diversos autores mostraram que a modelagem
∗
Trabalho realizado como parte do relatório final do projeto de iniciação cientı́fica entitulado: Das
Transformadas Integrais ao Cálculo Fracionário aplicado à Engenharia de Produção, Bolsa Fapesp, pro-
cesso 2011/07241-8. sob a Orientação do Prof. Dr. Rubens de Figueiredo Camargo.
†
Email: pedrof ramos@hotmail.com e rubens@fc.unesp.br, respectivamente, curso de Engenharia de
Produção, UNESP-Bauru e docente do departamento de Matemática da UNESP de Bauru.
________________________________________________________________________
RAMOS, P. F. P.; CAMARGO, R. F. Cálculo fracionário aplicado ao problema da tautócrona. C.Q.D. - Revista Eletrônica Paulista de
Matemática, Bauru, v. 1, p. 15-22, dez. 2012.
DOI: 10.21167/cqdvol1201223169664pfprrfc1522 - Disponível em: http://www2.fc.unesp.br/revistacqd/index.jsp
15
feita a partir do cálculo fracionário oferece uma descrição mais fina de fenômenos
naturais que aquela feita a partir do cálculo usual, tendo em vista que derivadas
fracionárias proporcionam uma excelente descrição para efeitos de memória e pro-
priedades hereditárias [13].
A maneira canônica de se utilizar esta importante ferramenta é substituir a
derivada de ordem inteira de uma equação diferencial parcial, que descreve um
determinado fenômeno, por uma derivada de ordem não-inteira. Vários resultados
importantes e generalizações foram obtidos através desta técnica, em diversas áreas
do conhecimento, tais como: mecânica dos fluidos, fenômenos de transporte, redes
elétricas, probabilidade, biomatemática, dentre outros [1, 5, 8, 9, 10, 11, 12, 13, 14,
17, 18].
No presente trabalho, introduzimos os conceitos de integral e de derivada de
ordens não inteira, segundo Riemann Liouville e utilizamos estes conceitos para
resolver o problema da tautócrona, que consiste em se determinar a curva na qual
o tempo gasto por um objeto para deslizar, sem fricção, em gravidade uniforme até
seu ponto de mı́nimo é independente de seu ponto de partida [15].
1 Integral Fracionária
Vamos introduzir o conceito de integral fracionária como uma generalização da
integral de ordem inteira. Para tanto, vamos mostrar que uma integral de ordem
n, de uma função integrável f (x), pode ser vista como um produto de convolução
de Laplace entre f (x) e a função Gel’fand-Shilov de ordem n. Em seguida, vamos
utilizar a generalização do conceito de fatorial através da função gama para gene-
ralizar o conceito de integral de ordem inteira.
Definição 2: Função de Gel’fand-Shilov
Sejam n ∈ N e ν 6∈ N, definimos a função de Gel’fand-Shilov como
 tn−1
  ν−1
 t
se t ≥ 0 se t ≥ 0
φn (t) := (n − 1)! e φν (t) := Γ(ν)
0 se t < 0 0 se t < 0.
 
Definição 3: Operador Integral de ordem n

Definimos a integral de ordem inteira através do operador I definido como
Z t
If (t) = f (t1 ) dt1 .
0
Desta forma temos

Z tZ t1
I 2 f (t) = I[If (t)] = f (t2 ) dt2 dt1 .
0 0
De maneira análoga para o operador integral de ordem n, isto é,

Z t Z t1 Z t2 Z tn−1
I n f (t) = ··· f (tn ) dtn dtn1 . . . dt2 dt1 .
0 0 0 0
Definição 4: Convlução de Laplace

Definimos a convolução, ou produto de convolução de Laplace, de f (t) e g(t),
denotada por f (t) ∗ g(t), como sendo [7]
________________________________________________________________________
16
Z t Z t
f (t) ∗ g(t) = f (t − τ )g(τ ) dτ = f (τ )g(t − τ ) dτ .
0 0
Teorema 1: Integral de ordem n

Seja f (x) uma função integrável, então
Z t t
(t − τ )n−1
Z
I n f (t) = φn (t) ∗ f (t) = φn (t − τ )f (τ ) dτ = f (τ ) dτ. (1.0.1)
0 0 (n − 1)!
Demonstração: Vamos provar o teorema por indução em n. Para n = 1 temos

que Z t Z t
(t − τ )1−1
If (t) = f (τ ) dτ = f (τ ) dτ = φ1 (t) ∗ f (t).
0 0 (1 − 1)!
Para concluir basta mostrar que se I n f (t) = φn (t) ∗ f (t) então I n+1 f (t) =
I[I n f (t)] = φn+1 (t) ∗ f (t). Por hipótese indutiva temos que
Z t Z tZ u
n+1 (u − τ )n−1
I f (t) = I[φn (t) ∗ f (t)] = φn (u) ∗ f (u)du = f (τ ) dτ du.
0 0 0 (n − 1)!
Pelo teorema de Goursat [13] podemos mudar a ordem de integração, i.e.,

Z t Z t
(u − τ )n−1

n+1
I f (t) = f (τ ) du dτ.
0 τ (n − 1)!
Calculando a integral entre colchetes podemos escrever

Z t
n+1 (t − τ )n
I f (t) = f (τ ) dτ = φn+1 (t) ∗ f (t),
0 n!
que é justamente o resultado desejado, a partir do qual apresentamos a definição
que se segue.
Definição 5: Integral de ordem arbitrária segundo Riemann-Liouville.
Seja f (t) uma função integrável, utilizamos a generalização da função fatorial
pela função gama [4] para definirmos a integral de ordem ν de f (t), denotada por
I ν f (t), por Z t
(t − τ )ν−1
I ν f (t) = φν (t) ∗ f (t) = f (τ ) dτ. (1.0.2)
0 Γ(ν)
2 Derivada Fracionária
Há várias formas de se introduzir a derivada de ordem não-inteira como uma
generalização para a derivada de ordem inteira, dentre elas podemos citar a definição
de Riemann-Liouville, que é a mais conhecida e a de Caputo, que é mais restritiva,
mas parece ser mais adequada para o estudo de problemas fı́sicos [9, 10, 11, 12]. Além
disso, destacamos a definição de Grünwald-Letnikov que é mais apropriada para se
utilizar em problemas numéricos. Estas e outras definições podem ser encontradas
em detalhes no livro de Podlubny [18].
No presente trabalho, estamos interessados na resolução de uma equação inte-
grodiferencial fracionária, relacionada ao problema da tautócrona, por esta razão
apresentamos apenas a derivada fracionária segundo Riemann-Liouville, que como
veremos, baseia-se no fato da derivação ser o operador inverso à esquerda da inte-
gração e na lei dos expoentes de Lagrange [13].
Definição 6: Derivada fracionária segundo Riemann-Liouville.
Sejam f (t) uma função diferenciável, m ∈ IN e α 6∈ IN tais que m − 1 < Re(α) < m.
________________________________________________________________________
17
A derivada de ordem α no sentido de Riemann-Liouville é definida como sendo
a derivada de ordem inteira de uma integral fracionária, de forma que a lei dos
expoentes faça sentido, isto é,
D α f (t) = D m I m−α f (t) = D m [φm−α (t) ∗ f (t)]. (2.0.3)
3 O Problema da Tautócrona
Apresentamos aqui a primeira aplicação do cálculo fracionário presente na li-
teratura, ou seja, a solução proposta por Abel para o problema da Tatócrona ou
curva isocrônica, isto é, o problema de se determinar a curva na qual o tempo gasto
por um objeto para deslizar, sem fricção, em gravidade uniforme até seu ponto de
mı́nimo é independente de seu ponto de partida.
Vamos resolver este problema de duas maneiras distintas a primeira utilizando
o cálculo usual e a segunda utilizando o cálculo fracionário.
A solução proposta por Abel parte do Princı́pio da Conservação de Energia
que estabelece que a quantidade total de energia em um sistema isolado permanece
constante, ou ainda, que a soma entre energia potencial gravitacional e a energia
cinética é constante.
Por hipótese, no problema da tautócrona, a partı́cula se move sem atrito e con-
sequentemente sua energia cinética é exatamente igual à diferença entre a energia
potencial em seu ponto inicial e a energia potencial no ponto em que se encontra.
Sendo m a massa do objeto, v(t) sua velocidade no instante t, y0 a altura em que
foi abandonado e y(t) a altura no instante t, sabemos que as energias cinética e
potencial são dadas, respectivamente, por
1
m v2 e m g y.
2
Como a partı́cula está restrita a mover-se sobre a curva, sua velocidade é simples-
mente v = ds/dt, onde s é a distância medida ao longo da curva. Utilizando o
princı́pio da conservação de energia podemos escrever
2
1 ds
m = mg(y0 − y)
2 dt
ds p
= ± 2g(y0 − y)
dt
ds
dt = ±p
2g(y0 − y)

1 ds
dt = ± √ (y − y0 )−1/2 dy,
2g dy
onde a última passagem é devida à regra da cadeia. Note que a função s(y) descreve
a distância remanescente na curva em termos da altura remanecente y, como a
distância e a altura diminuem a medida que o tempo passa devemos tomar apenas
o sinal negativo na equação anterior, isto é;

1 ds
dt = − √ (y − y0 ) −1/2
dy. (3.0.4)
2g dy
Integrando a equação (3.0.4) em ambos os lados de y0 a zero temos
Z 0 Z 0
1 ds
τ = T (y0 ) = dt = − √ (y − y0 )−1/2
dy,
y0 2g y0 dy
________________________________________________________________________
18
ou seja,
y0
1 ds
Z
τ=√ (y − y0 )
−1/2
dy, (3.0.5)
2g 0 dy
na qual τ é o tempo de descida. Esta equação integrodiferencial pode ser resolvida
de várias formas, a seguir apresentamos a clássica solução que utiliza a transfor-
mada de Laplace e posteriormente a solução proposta por Abel utilizando o cálculo
fracionário.
3.1 Transformada de Laplace

Aplicando a transformada de Laplace em ambos os lados da equação (3.0.5)
podemos escrever
Z y0
1 −1/2 ds
L[τ ] = √ L (y0 − y) dy
2g 0 dy

1 1 ds
= √ L √ L ,
2g y dy
a última igualdade é devida ao teorema da convolução.
Sendo s o parâmetro da transformada de Laplace, temos que L[1] = 1/s e
1 √
L √ = π s−1/2 , logo
y √
ds τ 2g −1/2
L = √ s .
dy π
A fim de explicitar a solução, isto é, ds/dy, aplicamos a transformada de Laplace
inversa de onde obtemos √
ds τ 2g 1
= √ . (3.1.1)
dy π y
3.2 Solução de Abel

A solução proposta por Abel baseia-se na observação de que a integral da equação
√
(3.0.5), exceto pelo fator multiplicativo 1/Γ(1/2) = 1/ π é exatamente a definição
de integração fracionária de ordem 1/2. Desta forma aplicando a derivada segundo
Riemann-Liouville de ordem 1/2 em ambos os lados da equação temos1
√ Z y0
d1/2 π d1/2

1 ds
τ=√ (y − y0 )−1/2
dy . (3.2.1)
dy 1/2 2g dy 1/2 Γ(1/2) 0 dy
Lembrando que a derivada de Riemann-Liouvile de ordem 1/2 de 1 = y 0 vale

√
( π y)−1/2 podemos escrever
√
τ π ds
√ =√ ,
πy 2g dy
ou seja, √
ds τ 2g 1
= √
dy π y
1
Lembre-se que a derivada fracionária de Riemann-Liouville é o operador inverso à esquerda da integral
fracionária.
________________________________________________________________________
19
que é o mesmo resultado obtido na equação (3.1.1). Note que a solução via cálculo
fracionário é bem mais simples e direta que a solução clássica. Essa equação
(variáveis separáveis) nos mostra que
√
τ 2g 1
ds = √ dy,
π y
de onde segue, após a integração, que

√
2 τ 2g 1/2
s(y) = y .
π
3.3 Representação gráfica

A fim de ilustrar a solução do problema da tautócrona apresentamos a seguir
uma representação gráfica da curva correspondente com cinco massas abandonadas
em pontos diferentes, e verificamos que todas atingem o vértice ao mesmo tempo.
S(t)
t
0.5 1.0 1.5
-0.2
-0.4
-0.6 t=0
-0.8
-1.0
S(t)
t
0.5 1.0 1.5
-0.2
-0.4
-0.6 t = 0.6τ
-0.8
-1.0
S(t)
t
0.5 1.0 1.5
-0.2
-0.4
-0.6 t = 0.9τ
-0.8
-1.0
Figura 1: Posição das partı́culas para diferentes valores de t.
________________________________________________________________________
20
S(t)
2.0
1.5
y0 = 0.2
1.0
y0 = 0.4
0.5
y0 = 0.6 y0 = 1
y0 = 0.8
t
0.2 0.4 0.6 0.8 1.0
Figura 2: Caminho percorrido por cada partı́cula.
4 Conclusões
Concluı́mos que além de ser uma ferramente precisa para refinar a descrição de
fenômenos naturais, o cálculo fracionário também pode ser utilizado para obter de
maneira simples e elegante a solução de problemas de cálculo variacional, tais como
o da Tautócrona.
5 Agradecimentos
Gostarı́amos de agradecer à FAPESP por ter financiado este projeto, processo
2011/07241-8. Além disso, agradecemos ao Prof. Dr. Alexys Bruno Alfonso por
sua indispensável ajuda com a parte gráfica e aos pareceristas por suas valiosas
contribuições e correções.
Referências
[1] A. Luisa Soubhia, R. Figueiredo Camargo, E. Capelas de Oliveira, J. Vaz Theo-
rem for Series in Thee-Parameter Mittag-Leffler Function, Fractional Calculus
and Applied Analysis, bf13, (2010).
[2] H. Bernhard The Bernoulli family, In: Wussing, H. Arnold, W. Biographien
bedeutender Mathematiker. Berlin, 1983.
[3] W. E. Boyce and R. C. DiPrima, Equações Equações Diferenciais Elementares
e Problemas de Valores de Contorno, Oitava Edição, LTC, Rio de Janeiro,
(2006).
[4] E. Capelas de Oliveira, Funções Especiais com Aplicações, Editora Livraria da
Fı́sica, São Paulo, (2005).
[5] L. Debnath, Recent Applications of Fractional Calculus to Science and Engi-
neering, Int. J. Math. Math. Sci., 54, 3413-3442, (2003).
[6] Felix Silva Costa ; Capelas de Oliveira, E. ; R. Figueiredo Camargo,; J.
Vaz. Integrais de Mellin-Barnes e a Função de Fox. TEMA. Tendências
em Matemática Aplicada e Computacional, v. 11, p. 01/01, (2011).
[doi:10.5540/tema.2011.012.02.0157]
________________________________________________________________________
21
[7] R. Figueiredo Camargo, Do Teorema de Cauchy ao Método de Cagniard, Dis-
sertação de Mestrado, Universidade Estadual de Campinas, Campinas, (2005).
[8] R. Figueiredo Camargo, E. Capelas de Oliveira and F. A. M. Gomes, The
Replacement of Lotka-Volterra Model by a Formulation Involving Fractional
Derivatives, R. P. 10/07, (2007).
[9] R. Figueiredo Camargo, Ary O. Chiacchio and E. Capelas de Oliveira, Differen-
tiation to Fractional Orders and the Fractional Telegraph Equation, J. Math.
Phys., 49, 033505, (2008) [doi:10.1063/1.2890375].
[10] R. Figueiredo Camargo, R. Charnet and E. Capelas de Oliveira, On
Some Fractional Green’s Functions, J. Math. Phys. 50, 043514 (2009).
[doi:10.1063/1.3119484].
[11] R. Figueiredo Camargo, Ary O. Chiacchio, R. Charnet and E. Capelas de
Oliveira, Solution of the fractional Langevin equation and the Mittag-Leffler
functions, J. Math. Phys. 50, 063507 (2009). [doi:10.1063/1.3152608].
[12] R. Figueiredo Camargo, E. Capelas de Oliveira and J. Vaz Jr. On anoma-
lous diffusion and the fractional generalized Langevin equation for a harmonic
oscillator J. Math. Phys. 50, 123518 (2009). [doi:10.1063/1.3269587].
[13] R. Figueiredo Camargo, Cálculo Fracionário e Aplicações, Tese de Doutorado,
Unicamp, Campinas, (2009).
[14] R. Figueiredo Camargo,; Oliveira, Edmundo Capelas ; Vaz, Jayme. On the
Generalized Mittag-Leffler Function and its Application in a Fractional Tele-
graph Equation. Mathematical Physics, Analysis and Geometry , v. 14, p. 274,
(2011). [doi:10.1007/s11040-011-9100-8]
[15] R. Figueiredo Camargo e E. Capelas de Oliveira. Introdução ao Cálculo de
Ordem não Inteira., livro em fase final de editoração.
[16] K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Frac-
tional Differential Equations, John Wiley & Sons, Wiley-Interscience, (1993).
[17] G. Pedroso Pires Abreu e R. Figueiredo Camargo, Modelagem Matemática: Das
Transformadas Integrais ao Cálculo de Ordem não-Inteira, XXIII Congresso de
Iniciação Cientı́fica da UNESP, Bauru, (2011).
[18] I. Podlubny, Fractional Differential Equations, Mathematics in Science and
Engineering, Vol.198, Academic Press, San Diego, (1999).
________________________________________________________________________
22

V01a03 Calculo Fracionario Aplicado

Enviado por

Direitos autorais:

Formatos disponíveis

V01a03 Calculo Fracionario Aplicado

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

V01a03 Calculo Fracionario Aplicado

Enviado por

Direitos autorais:

Formatos disponíveis

Cálculo Fracionário Aplicado ao Problema da

Definição 2: Função de Gel’fand-Shilov

Sejam n ∈ N e ν 6∈ N, definimos a função de Gel’fand-Shilov como

Definição 3: Operador Integral de ordem n

Desta forma temos

De maneira análoga para o operador integral de ordem n, isto é,

Definição 4: Convlução de Laplace

Teorema 1: Integral de ordem n

Demonstração: Vamos provar o teorema por indução em n. Para n = 1 temos

Pelo teorema de Goursat [13] podemos mudar a ordem de integração, i.e.,

Calculando a integral entre colchetes podemos escrever

D α f (t) = D m I m−α f (t) = D m [φm−α (t) ∗ f (t)]. (2.0.3)

3.1 Transformada de Laplace

3.2 Solução de Abel

Lembrando que a derivada de Riemann-Liouvile de ordem 1/2 de 1 = y 0 vale

de onde segue, após a integração, que

3.3 Representação gráfica

Figura 1: Posição das partı́culas para diferentes valores de t.

Figura 2: Caminho percorrido por cada partı́cula.

Você também pode gostar