Enunciado

Considere o problema de valor inicial

  y '   =   – 10 y   +   2,5 t 2   +   0,5 t ,   y ( 0 ) =   4 .  

(a) Encontre a solução y   =   ϕ ( t ) .

(b) A análise de estabilidade no texto sugere que, para esse problema, o método de Euler só é estável para h   <   0,2 . Confirme se isto é verdade, aplicando o método de Euler a esse problema para 0   ≤   t   ≤   5 com tamanhos de passos próximos de 0,2 .

(c) Aplique o método de Runge-Kutta a esse problema para 0   ≤   t   ≤   5 com diversos tamanhos de passos. O que você pode concluir sobre a estabilidade desse método?

(d) Aplique o método de Euler inverso a esse problema para 0   ≤   t   ≤   5   com diversos tamanhos de passos. Qual o tamanho de passo necessário para que o erro em t   =   5 seja menor do que 0,01 ?

Passo 1

(a) Precisamos resolver um problema y ' = - 10 y   + 2,5 t ^ 2 + 0,5 t   com a condição inicial y ( 0 ) =   4 . Primeiro, resolva a equação homogênea correspondente y   ' =   -   10 y , teremos que a solução é:

y h = c e - 10 t

Agora encontre a solução específica. É da forma y p = A + B t + C t 2 , uma vez que deve ser uma solução do problema inicial, substituímos pelo problema para encontrar A e B .

y p ' = 10 y p + 2,5 t 2 + 0,05 t

B + 2 C t = - 10 A - 10 B t - 10 C t 2 + 2,5 t 2 + 0,5 t

Se:

A = B = 0

C = 1 4

A solução do ODE inicial é dada como y = y h + y p , portanto, a solução geral é y ( t ) = c - 10 t   +   0,25 t 2 . Mas temos uma condição inicial y ( 0 ) = 4 e podemos determinar c :

y 0 = c = 4

c = 4

A solução será:

y t = 4 e - 10 t + 0,25 t 2

Passo 2

Para a resolução do passo b usaremos as seguintes fórmulas:

x n + 1 = x n + h f n     ; ( 4 )

x n + 1 y n + 1 = x n y n + h f ( t n , x n , y n ) g ( t n , x n , y n )     ; ( 5 )

(b) O método Euler afirma que t n = t 0 = n h e y n   +   1 - y n   +   h . f   ( t n ,   y n ) , em que d y d t =   f   ( t ,   y ) e a condição inicial y ( t 0 ) = y 0 . O primeiro passo pode ser dado por:

t 1 = t 0 + h = 0,1

y 1 = y 0 + 0,1 - 10 y 0 + 2,5 . t 0 2 + 0,5 t 0 = 0

Nós calculamos os passoas até y ( 0,4 ) . Os valores obtidos estão na tabela a seguir.

Passo 3

(c) Para este problema, usaremos:

y n + 1 = y n + h k n 1 + 2 h k n 2 + k n 3 + h k n 4 6  

Podemos ter então :

t 1 = t 0 + h = 0 + 0,1

k n 1 = f t n , y n = h f 0,1 = - 4

k n 2 = f t n + 0,5 h ; y n + 0,5 h k n 1 = h f 0,05 ; 2 = - 2

k n 3 = f ( t n + 0,5 h ; y n + 0,5 h k n 2 = h f 0,05 ; 3 = - 3

k n 4 = f ( t n + h ; y n + h k n 3 = h f 0,1 ; 1 = - 1

y 1 = y 0 + h k n 1 + 2 k n 2 + k n 3 + k n 4 6 = 1,5

Repetimos os passos até encontrar y ( 5 ) :

Resposta

  1. y t = 4 e - 10 t + 0,25 t 2
  2. 6,2500
  3. 6,2501, estável com h = 0,1 .
Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

3 - x^' = - t x - y - 1 , y ' = x ; x ( 0 ) = 1 , y ( 0 ) = 1
9 - Considere o problema de valor inicial x ' = f ( t ,x ,y ) e y ' = g ( t ,x ,y ) , com x t 0 = x 0 e y ( t 0 ) = y 0 .A generalização do método de previsão e correção de Adam-Moulton da Seção 8.4 é
8- Considere o problema de valor inicial x ' ' + t 2 x ' + 3 x = t , x 0 = 1 , x ' 0 = 2 Transforme esse problema em um sistema de duas equações de primeira ordem e determinevalores aproximados da sol
5 - x^' = x ( 1 - 0,5 x - 0,5 y ) , y^' = y ( - 0,25 + 0,5 x ) ; x ( 0 ) = 4 , y ( 0 ) = 1
4 - x^' = x - y + x y , y^' = 3 x - 2 y - x y ; x ( 0 ) = 0 , y ( 0 ) = 1