Fichas de Trabalho EDO Solucoes
Fichas de Trabalho EDO Solucoes
Fichas de Trabalho EDO Solucoes
Exercício 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Considere a equação diferencial y 0 = y(x2 1) com y(0) = 1 e x 2 [0, 1].
Verifique se a função f (x, y) verifica a condição de Lipschitz.
Calcule as soluções aproximadas usando os métodos de Euler progressivo e regressivo com passo
h = 0.1.
x3 3x
Determine o erro das aproximações, sabendo y = e 3 é a solução da equação diferencial.
Solução:
Pág. 1 de 9
Matemática
Pág. 2 de 9
Matemática
Exercício 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Resolva a equação diferencial do exercício anterior usando o método de Taylor de ordem 2.
Solução:
o problema, a função f (x, y(x)) = y · (x2 1), a fórmula iterativa é dada pela seguinte
expressão:
h2
yi+1 = yi + h · yi · (x2i 1) + · yi · (2xi + (x2i 1)2 )
2
pois, para f (x, y) = y(x2 1), pela regra da derivada da função composta, f 0 (x, y(x)) =
fx0 (x, y) + fy0 (x, t)y 0 (x), obtendo-se: f 0 (x, y) = 2xy + y(x2 1)2 .
Considerando o passo h = 0.1 e a condição inicial y(0) = 1, obtém-se os seguintes valores:
xn 0.00000 0.10000 0.20000 0.30000 0.40000 0.50000 0.60000 0.70000 0.80000 0.90000 1.00000
yn 1.00000 0.90500 0.82074 0.74738 0.68470 0.63234 0.58986 0.55685 0.53308 0.51849 0.51340
y(xn ) 1.00000 0.90514 0.82092 0.74752 0.68477 0.63234 0.58978 0.55673 0.53295 0.51840 0.51342
" 0.00000 0.00014 0.00017 0.00014 0.00007 0.00001 0.00007 0.00012 0.00013 0.00009 0.00001
Comparando os resultados com os obtidos pelos métodos de Euler progressivo e regressivo
(Exercício 1), o valor do erro é significativamente menor (o número de casas decimais exatas
aumentou de 1 para 3) e a distribuição ao longo da solução não apresenta acumulação de
erros.
Pág. 3 de 9
Matemática
Exercício 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Suponha que se pretende resolver o seguinte problema de condições iniciais pelo método de
Euler. (
y0 = y
y(0) = 1, x 2 [0, 1]
Prove que se usar como passo h = 0.002, as aproximações yn de y(xn ) são obtidas com 3 casas
decimais significativas qualquer que seja o xn 2 [0, 1].
Solução:
|Th |
ry n = < |Th |, porque a solução do problema é a função y(x) = ex e para x 2 [0, 1],
|y(xn )|
o valor de y(x) 2 [1, e].
h
Dado que |Th | < maxx2[0,1] |f 0 (x, y)| e que f 0 (x, y) = y 00 = ex , maxx2[0,1] |f 0 (x, y)| =
2
h 0.002
maxx2[0,1] e = e = e. Pelo que |Th | < e =
x 1 e ⇡ 2.7 ⇥ 10 3 < 5 ⇥ 10 3 , para
2 2
qualquer xn 2 [0, 1].
Pág. 4 de 9
Matemática
Exercício 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dado y 0 = x y e y(0) = 2.0, determine y(0.2) e y(0.4).
Solução:
Pág. 5 de 9
Matemática
Exercício 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Utilizando um método de Taylor, determine y(1) com um erro máximo de 5 ⇥ 10 2 sabendo
que max |y (4) (x)| < 11.1,
x2[0,1]
(
y 0 y = sin(x) + cos(2x)
y(0) = 1
Solução:
h2 00 h3
y(xi+1 ) = y(xi ) + h · y 0 (xi ) + y (xi ) + y (3) (xi )
2! 3!
h4 h4 0.05 · 24
Pelo que, |Th | maxx2[0,1] |y (4) | · 11.1 < 5 ⇥ 10 2, portanto, h4 < e
4! 24 11.1
h < 0.573.
Considerando que y 0 = y + sin(x) + cos(2x), y 00 = y 0 + cos(x) 2 sin(2x),
y (3) = y 00 sin(x) 4 cos(2x)
O método iterativo pode ser implementado da seguinte forma:
8 2 3
>
>
> yi+1 = yi + h · C1 + h2 · C2 + h6 · C3
>
<C = y + sin(x ) + cos(2x )
1 i i i
>
> C 2 = C 1 + cos(x i ) 2 sin(2x i)
>
>
:C = C sin(xi ) 4 cos(2xi )
3 2
Pág. 6 de 9
Matemática
Exercício 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dado o PVI, determine y(1) e y(2) utilizando o método de Runge-Kutta de quarta ordem com
h = 0.5. 8
<y 0 = 2y + (x + 1)3
x+1
:y(0) = 3
Solução:
Pág. 7 de 9
Matemática
Exercício 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dado o problema de valores iniciais seguinte calcule y(0.50) e y(0.50), para h = 0.25:
8
> 00 0
<y = 4y 3y x
4
y(0) = 9
>
: 0
y (0) = 73
Solução:
O primeiro passo consiste em reescrever a equação de segunda ordem num sistema de duas
equações de primeira ordem:
8
>
> y1 = y
>
>
<y 0 = y
1 2
> y 0 = f (x, y , y ) = 4y 3y1 x
>
> 2 1 2 2
>
:y (0) = 4 , y (0) = 7
1 9 2 3
Usando o método de Euler para resolução de cada uma das equações do sistema, obtém-se
o seguinte esquema iterativo:
8 4
>
> y1,0 =
>
> 9
>
< 7
y2,0 =
> 3
>
> y = y1,i + h · y2,i
>
> 1,i+1
:
y2,i+1 = y2,i + h(4y2,i 3y1,1 xi )
Pág. 8 de 9