Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Calculo Numerico

Fazer download em pdf ou txt
Fazer download em pdf ou txt
Você está na página 1de 21

MS211 - Cálculo Numérico

Aula 22 – Métodos Numéricos para Sistemas de Equações


Diferenciais e Equações de Ordem Superior.

Marcos Eduardo Valle

Marcos Eduardo Valle MS211 - Cálculo Numérico 1 / 20


Na aula anterior, apresentamos os métodos de Runge-Kutta para
resolução de um PVI da forma:
(
y 0 = f (x, y ),
y (x0 ) = y0 .

Em termos gerais, um método de Runge-Kutta define

yk +1 = yk + hφ(xk , yk ), ∀k = 0, 1, . . . ,

em que φ é uma função de x e y que depende indiretamente de f e


do tamanho do passo h.

Na aula de hoje, veremos como os métodos de Runge-Kutta podem


ser aplicados também para a resolução numérica de sistemas de
equações diferenciais e equações de ordem superior, ambos com
valores iniciais.
Marcos Eduardo Valle MS211 - Cálculo Numérico 2 / 20
Sistemas de Equações Diferenciais
Um sistema de equações diferenciais com valor inicial


 y10 = f1 (x, y1 , y2 , . . . , yn ), y1 (x0 ) = y1,0 ,

y 0 = f2 (x, y1 , y2 , . . . , yn ),

y2 (x0 ) = y2,0 ,
2
.
.. ..


 .
yn (x0 ) = yn,0 ,

 0
yn = f1 (x, y1 , y2 , . . . , yn ),
pode ser escrito de forma mais compacta como
y0 = f(x, y), y(x0 ) = y0 ,
em que x ∈ R e y(x), f(x, y) e y0 são vetores em Rn escritos como
     
y1 (x) f1 (x, y) y1,0
y2 (x) f2 (x, y) y2,0 
y(x) =  .  , f(x, y) =  .  e y0 =  .  .
     
 ..   ..   .. 
yn (x) fn (x, y) yn,0
Marcos Eduardo Valle MS211 - Cálculo Numérico 3 / 20
Métodos de Euler

O método de Euler para sistemas de equações diferenciais é

yk +1 = yk + hf(xk , yk ), ∀k = 0, 1, . . . ,

em que y1 , y2 , . . . são aproximações para y(x1 ), y(x2 ), . . ..

Exemplo 1
Use o método de Euler para obter uma aproximação numérica da
solução do sistema de equações diferenciais
(
y 0 = z, y (0) = 1
0
z =y +e x z(0) = 0

para x ∈ [0, 0.2] usando h = 0.1.

Marcos Eduardo Valle MS211 - Cálculo Numérico 4 / 20


Resolução: A fórmula do método de Euler fornece
     
yk +1 yk zk
= +h .
zk +1 zk yk + exk

Para k = 0, temos
           
y1 y0 z0 1 0 1
= +h x = + 0.1 = .
z1 z0 y0 + e 0 0 1+1 0.2

Para k = 1, temos
           
y2 y1 z1 1 0.2 1.02
= +h = + 0.1 = .
z2 z1 y1 + ex1 0.2 1 + e0.2 0.4105

Portanto, encontramos as aproximações

y (0.2) ≈ 1.02 e z(0.2) ≈ 0.4105.

Marcos Eduardo Valle MS211 - Cálculo Numérico 5 / 20


Método de Heun
De um modo similar, o método de Heun para sistemas de equações
diferenciais com valor inicial é
h
yk +1 = yk + (k1 + k2 ), ∀k = 0, 1, . . . ,
2
em que
k1 = f(xk , yk ) e k2 = f(xk + h, yk + hk1 ).

Exemplo 2
Use o método de Heun para obter uma aproximação numérica da
solução do sistema de equações diferenciais
(
y 0 = z, y (0) = 1
0
z =y +e x z(0) = 0

para x ∈ [0, 0.2] usando h = 0.1.


Marcos Eduardo Valle MS211 - Cálculo Numérico 6 / 20
Resolução: Para k = 0, o método de Heun fornece:
        
y0 z0 0 0
k1 = f x0 , = = =
z0 y0 + ex0 1 + e0 2
e
        
y0 1 0
k2 = f x0 + h, + hk1 = f 0 + h, +h
z0 0 2
      
1 0.2 0.2
= f 0.1, = = .
0.2 1 + e0.1 2.1052

Logo,
         
y1 y h 1 0.1 0 0.2
= 0 + (k1 + k2 ) = + +
z1 z0 2 0 2 2 2.1052
 
1.01000
=
0.20526

Marcos Eduardo Valle MS211 - Cálculo Numérico 7 / 20


Para k = 1, o método de Heun fornece:
        
y1 z1 0.20526 0.20526
k1 = f x1 , = = =
z1 y1 + ex1 1.01 + e0.1 2.1152
e
   
y1
k2 = f x1 + h, + hk1
z1
    
1.01000 0.20526
= f 0.1 + h, +h
0.20526 2.1152
      
1.03053 0.41678 0.41678
= f 0.2, = = .
0.41678 1.03053 + e0.2 2.2519
Logo,
   
y2 y h
= 1 + (k1 + k2 )
z2 z1 2
       
1.01000 0.1 0.20526 0.41678 1.04110
= + + =
0.20526 2 2.1152 2.2519 0.42362

Marcos Eduardo Valle MS211 - Cálculo Numérico 8 / 20


Método de Runge-Kutta de Ordem 4

Finalmente, o método de Runge-Kutta de ordens 4 para sistemas de


equações diferenciais é descrito por:

h
yk +1 = yk + (k1 + 2k2 + 2k3 + k4 ) ,
6
em que

k1 = f(xk , yk ),
k2 = f(xk + h/2, yk + k1 h/2),
k3 = f(xk + h/2, yk + k2 h/2),
k4 = f(xk + h, yk + k3 h).

Marcos Eduardo Valle MS211 - Cálculo Numérico 9 / 20


Exemplo 3
Use o método de Runge-Kutta de ordem 4 para obter uma
aproximação numérica de
(
y 0 = z, y (0) = 1
0
z =y +e x z(0) = 0

para x ∈ [0, 0.2] usando h = 0.1.

Marcos Eduardo Valle MS211 - Cálculo Numérico 10 / 20


Resolução: Para k = 0, temos
    
y 0
k1 = f x0 , 0 = ,
z0 2
     
h y0 h 0.1000
k2 = f x0 + , + k1 =
2 z0 2 2.0513
     
h y0 h 0.1026
k3 = f x0 + , + k2 = ,
2 z0 2 2.0563
     
y0 0.2056
k4 = f x0 + h, + k3 h = .
z0 2.1154

Logo,  
h 1.0102
y1 = y0 + (k1 + 2k2 + 2k3 + k4 ) = .
6 0.2055

Marcos Eduardo Valle MS211 - Cálculo Numérico 11 / 20


Para k = 1, temos
    
y 0.2055
k1 = f x1 , 1 = ,
z1 2.1154
     
h y1 h 0.3113
k2 = f x1 + , + k1 =
2 z1 2 2.1823
     
h y1 h 0.3146
k3 = f x1 + , + k2 = ,
2 z1 2 2.1876
     
y1 0.4243
k4 = f x1 + h, + k3 h = .
z1 2.2630

Logo,  
h 1.0415
y2 = y1 + (k1 + 2k2 + 2k3 + k4 ) = .
6 0.4241

Marcos Eduardo Valle MS211 - Cálculo Numérico 12 / 20


Modelo Presa-Predador

Considere duas espécies que interagem como presa-predador. Por


exemplo, lobos e coelhos.

Denote por p(t) e q(t) a quantidade (e.g. densidade populacional)


de presas e predadores no instante t.

A dinâmica populacional das duas espécies pode ser descrita pelo


sistema de equações diferenciais, chamado equações de
Lotka-Volterra: (
p0 = ap − bpq,
q 0 = −cq + dpq,
em que a, b, c e d são constantes positivas.

Marcos Eduardo Valle MS211 - Cálculo Numérico 13 / 20


Considere

a = 0.25, b = 0.01, c = 1.00 e d = 0.01,

e a condição inicial

p0 = 80 e q0 = 30.

Equivalentemente, temos
(
p0 = 0.25p − 0.01pq, p(t0 ) = 80,
q 0 = −1.00q + 0.01pq, q(t0 ) = 30.

Nesse caso, encontramos as seguintes aproximações usando os


métodos de Euler, Heun e Runge-Kutta de ordem 4.

Marcos Eduardo Valle MS211 - Cálculo Numérico 14 / 20


Método de Euler

140
h=1.00
h=0.50
130 h=0.25
Exata
120

110
Presa

100

90

80

70

60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 15 / 20
Método de Heun

140
h=1.00
h=0.50
130 h=0.25
Exata
120

110
Presa

100

90

80

70

60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 16 / 20
Método de Runge-Kutta ordem 4

140
h=1.00
h=0.50
130 h=0.25
Exata
120

110
Presa

100

90

80

70

60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 17 / 20
Equações Diferenciais de Ordem Superior
Uma equação diferencial de ordem m
u (m) = g(x, u, u 0 , u 00 , . . . , u (m−1) ),
pode ser escrita como um sistema com m equações diferenciais de
ordem 1:
y0 = f (x, y),
tomando
y1 (x) = u(x), y2 (x) = u 0 (x), ..., ym (x) = u (m−1) .
Especificamente, temos


y10 = y2 ,

y 0

= y3 ,
2
..


 .

 0
ym = g(x, y1 , y2 , . . . , ym ).
Marcos Eduardo Valle MS211 - Cálculo Numérico 18 / 20
Exemplo 4
Escreva o problema de valor inicial
(
u 00 − u = ex ,
u(0) = 1 e u 0 (0) = 0,

como um sistema de equações diferenciais com valor de inicial.

Marcos Eduardo Valle MS211 - Cálculo Numérico 19 / 20


Exemplo 4
Escreva o problema de valor inicial
(
u 00 − u = ex ,
u(0) = 1 e u 0 (0) = 0,

como um sistema de equações diferenciais com valor de inicial.

Resposta: Tomando y = u e z = u 0 , podemos escrever o


problema de valor inicial como
(
y 0 = z, y (0) = 1
0
z =y +e , x z(0) = 0,

que é o mesmo sistema encontrado no Exemplo 1.

Marcos Eduardo Valle MS211 - Cálculo Numérico 19 / 20


Considerações Finais
Na aula de hoje, vimos que os métodos de Runge-Kutta podem ser
aplicados para aproximar a solução de um sistema de equações
diferenciais.

Vimos também como escrever uma equação diferencial de ordem


superior como um sistema de equações diferenciais.

Durante toda abordagem, consideramos h fixo. Porém, a maioria


dos softwares de computação científica usam métodos de passo
variado. Os métodos de passo variado fazem uma estimativa do
erro local a cada iteração. O método de Runge-Kutta-Fehlberg é um
exemplo de método de passo variado.

Muito grato pela atenção!

Marcos Eduardo Valle MS211 - Cálculo Numérico 20 / 20

Você também pode gostar