8.EQUAÇÕES DIFERENCIAIS ORDINÁRIAS Parte 1 8.1 – INTRODUÇÃO – PVI’s 8.2 – MÉTODOS DE PASSO SIMPLES 8.2.1 – MÉTODO DE EULER hoje 8.2.2 – MÉTODOS DE TAYLOR 8.2.3 – MÉTODOS DE RUNGE-KUTTA 8.3 – MÉTODOS DE PASSO MÚLTIPLO 8.4 – MÉTODOS PREVISOR-CORRETOR 8.5 – EDO’s DE ORDEM SUPERIOR E SISTEMAS DE EDO’s 8.6 - PVC’s E O MÉTODO DAS DIFERENÇAS FINITAS 8. EDO’s 8.1 INTRODUÇÃO Problemas de Valores Iniciais (PVI’s) Se dada uma EDO de ordem n, a função, assim como suas derivadas até ordem n-1, são especificadas em um único ponto, então temos um problema a valores iniciais. Exemplo: y x 1 y cos x y x com 2 1 y x 2 y y ( 0 ) 1 , y ( 0 ) 2 , y ( 0 ) 3 2 sen x y 8. EDO’s 8.1 INTRODUÇÃO Problemas de Valores no Contorno (PVC’s) Se para uma dada EDO de ordem n, as n condições forem dadas em diferentes pontos, então temos um problema a valores no contorno. Ao contrário dos PVI’s, os PVC’s podem não apresentar unicidade de solução. 8. EDO’s 8.1 INTRODUÇÃO Exemplo de PVC’s. 1- Seja um barra de comprimento L sujeita a uma carga uniforme q. Se em x=0 ela está fixada e em x=L ela está apoiada, então temos o problema y 4 k yq y ( 0 ) y ( 0 ) 0 com y ( L ) y ( L ) 0 8. EDO’s 8.1 INTRODUÇÃO PVC’s sem unicidade na solução. 2- O problema y 0 y ( 1) 0 com y (1) 2 y (1) 0 tem como solução y ( x ) 1 x para todo 8. EDO’s 8.1 INTRODUÇÃO Nesta primeira parte do estudo de EDO’s abordaremos métodos para resolução de PVI’s de primeira ordem. Dado o PVI y f x , y com y x 0 y 0 construiremos x 1 , x 2 , x 3 , ..... , x n , para simplificar igualmente espaçados, ou seja, x i 1 x i h para i 0 ,1, 2 ,....., n 1 e calculamos as aproximações y i y x i neste pontos. 8. EDO’s 8.1 INTRODUÇÃO Se para calcular y i y x i , usamos apenas y i 1 y x i 1 , então dizemos que o Método é de Passo Um ou de Passo Simples. Porém se usarmos mais valores teremos um Método de Passo Múltiplo. Para PVI’s de primeira ordem temos que y x 0 y 0 é uma aproximação inicial para a solução. Problema auto-iniciante. Para Métodos de Passos Múltiplos devemos ter estratégias para as aprox. iniciais. 8. EDO’s 8.1 INTRODUÇÃO Os Métodos de Passos Simples têm as seguintes características: 1) Deve-se calcular os valores de f ( x , y ) e de suas derivadas em muitos pontos. Fator negativo. 2) A estimativa dos erros não é trivial. 8. EDO’s 8.2.1 Método de Euler Considere o PVI y d dt y f t , y com y t 0 y 0 (1) Suponha que exista uma única solução t n do problema no intervalo de interesse. Reescrevendo (1) no ponto t t n d dt t n f t n , t n (2) 8. EDO’s 8.2.1 Método de Euler Direto Aproximando a derivada em (2) pelo quociente de diferenças para frente (ou direto), obtemos t n 1 t n t n 1 t n f t n , t n (3) Substituindo t n 1 y n 1 e t n y n por seus valores aproximados, t n 1 t n h , temos a fórmula de Euler: y n 1 y n f t n , y n t n 1 t n (4) 8. EDO’s 8.2.1 Método de Euler Direto Outra maneira para obter a fórmula de Euler é escrever o problema (1) como uma equação integral. Integrando (1) de t n até t n 1 t n 1 tn ( t ) dt t n 1 f t , ( t ) dt tn ( t n 1 ) ( t n ) t n 1 tn f t , ( t ) dt (5) 8. EDO’s 8.2.1 Método de Euler Direto Se aproximarmos a integral substituindo f t , ( t ) por f t n , ( t n ) E como t n 1 t n h , então t n 1 t n f t n , t n t n 1 t n t n 1 t n h f t n , t n y n 1 y n f t n , y n h 8. EDO’s 8.2.1 Método de Euler Direto A integral em (5) é a área abaixo da curva vermelha. No Método de Euler Direto é a área lilas. y f t , t f t n 1 , t n 1 f t n , t n tn t n 1 8. EDO’s 8.2.1 Método de Euler Direto Note que quanto menor forem as partições, melhor será a convergência do Método de Euler. O Erro da fórmula de Euler pode ser majorado através da fórmula de Taylor. Seja y t em torno de t t n t n h t n t n h t n t n h t n f t n , t n com t t n , t n 1 h 2 ... 2! h t h 2 2! ERRO 8. EDO’s 8.2.1 Método de Euler Direto Note que sendo M max t com t t n , t n 1 onde (t) f ( t , ( t )) e (t) f t ( t , ( t )) f y ( t , ( t )) (t) f t f y f então o erro devido ao truncamento de Euler é majorado por 2 en M h 2 8. EDO’s 8.2.1 Método de Euler Direto Exemplo 1: Considere o problema de valor inicial y 1 t 4 y com y ( 0 ) 1 A solução exata é dada por y ( y ) 3 16 1 4 t 19 e 4t 16 Utilizando a fórmula de Euler (direta) e passos h 0 .05 , h 0 .025 , h 0 .010 e h 0 .001 determine a solução do problema no intervalo 0 t 2 8. EDO’s 8.2.1 Método de Euler Direto Solução por Euler direta de y 1 t 4 y com y (0) 1 t h=0.05 h=0.025 h=0.01 h=0.001 Exata 0.0 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 0.1 1.5475000 1.5761188 1.5952901 1.6076289 1.6090418 0.2 2.3249000 2.4080117 2.4644587 2.5011159 2.5053299 0.3 3.4333560 3.6143837 3.7390345 3.8207130 3.8301388 0.4 5.0185326 5.3690304 5.6137120 5.7754845 5.7942260 0.5 7.2901870 7.9264062 8.3766865 8.6770692 8.7120041 1.0 45.588400 53.807866 60.037126 64.382558 64.897803 1.5 282.07187 361.75945 426.40818 473.55979 479.25919 2.0 1745.6662 2432.7878 3029.3279 3484.1608 3540.2001 8. EDO’s 8.2.1 Método de Euler Direto Note que os erros gerados em t=2.0 são grandes! Para h=0.001, ou seja, 2000 subintervalos, temos um erro acumulado de 1.6% Como 2 h e n 1 ( t ) 2 3 4 t 19 e (t ) 19 e 4t (t ) 16 4t 8. EDO’s 8.2.1 Método de Euler Direto O erro devido ao truncamento local é 2 h 4t e n 1 19 e n para t n t n t n h 2 Para ir de t=1.95 a t=2.0, quando h=0.05, 19 e 0 .0025 e 7 .8 57 . 96 16 19 e 0 .0025 70 .80 8 40 16 Para obter um erro local de truncamento de 0.01 neste problema necessitamos de h=0.0006 em torno de t=2 e h=0.03 em torno de t=0. Tais métodos com erros constante são chamados ADAPTATIVOS. 8. EDO’s 8.2.1 Método de Euler Inverso Uma variante do método de Euler, chamado Método de Euler Inverso, consiste em aproximar a derivada em y d dt y f t , y com y t 0 y 0 (1) pelo quociente de diferenças para trás (ou inverso) t n t n 1 t n t n 1 f t n , t n (6) 8. EDO’s 8.2.1 Método de Euler Inverso Substituindo t n y n e t n 1 y n 1 por seus valores aproximados, e fazendo n n 1 temos a fórmula de Euler inversa y n 1 y n f t n 1 , y n 1 h (7) Note que a fórmula de Euler inversa fornece o valor de y n 1 de forma implícita. 8. EDO’s 8.2.1 Método de Euler Inverso A integral em (5) é a área abaixo da curva vermelha. No Método de Euler Inverso é a área verde. y f t , t f t n 1 , t n 1 f t n , t n tn t n 1 8. EDO’s 8.2.1 Método de Euler Inverso Exemplo 2: Considere o problema de valor inicial y 1 t 4 y com y ( 0 ) 1 A solução exata é dada por y ( y ) 3 16 1 4 t 19 e 4t 16 Utilizando a fórmula de Euler (inversa) e passos h 0 .05 , h 0 .025 , h 0 .010 e h 0 .001 determine a solução do problema no intervalo 0 t 2 8. EDO’s 8.2.1 Método de Euler Inverso Solução por Euler inversa de y 1 t 4 y com y (0) 1 é dada pela fórmula de Euler inversa y n 1 y n h 1 t n 1 4 y n 1 O primeiro passo gera: y 1 y 0 h 1 t 1 4 y 1 1 0.05 1 0 . 05 4 y 1 y 1 1 . 309375 y 2 y 1 h 1 t 2 4 y 2 1 . 309375 0.05 1 0 . 1 4 y 2 y 2 1 . 6929688 Continuando, temos a tabela: 8. EDO’s 8.2.1 Método de Euler Inverso Solução por Euler inversa de y 1 t 4 y com y (0) 1 t h=0.05 h=0.025 h=0.01 h=0.001 Exata 0.0 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 0.1 1.6929688 1.6474375 1.6236638 1.6104634 1.6090418 0.2 2.7616699 2.6211306 2.5491368 2.5095731 2.5053299 0.3 4.4174530 4.0920886 3.9285724 3.8396379 3.8301388 0.4 6.9905516 6.3209569 5.9908303 5.8131282 5.7942260 0.5 10.996956 9.7050002 9.0801473 8.7472667 8.7120041 1.0 103.06171 80.402761 70.452395 65.419964 64.897803 1.5 959.44236 661.00731 542.12432 485.0825 479.25919 2.0 8934.0696 5435.7294 4172.7228 3597.4478 3540.2001 8. EDO’s 8.2.1 Método de Euler Aprimorado Note que tanto o método de Euler direto quanto o inverso geram erros acumulativos quando t cresce. No exemplo o erro foi da ordem de 1.2%. Os Métodos adaptativos de Euler são uma solução, contudo teremos uma sub-rotina para calcular o tamanho do passo para cada n. Fórmula de Euler Aprimorada ou centrada aproxima a função f na integral por uma média. 8. EDO’s 8.2.1 Método de Euler Aprimorado A fórmula de Euler Aprimorada escrevese como: y n 1 y n f t n 1 , y n 1 f t n , y n h (8) 2 Os erros são menores e convergência é mais rápida neste caso. 8. EDO’s 8.2.1 Método de Euler Aprimorado A integral em (5) é a área abaixo da curva vermelha. No Método de Euler Aprimorado é a área amarela. y f t , t f t n 1 , t n 1 f médio f t n , t n tn t n 1 8. EDO’s 8.2.1 Método de Euler Aprimorado A solução por Euler Aprimorado de y 1 t 4 y com y (0) 1 é dada pela fórmula y n 1 y n h 2 1 t n 1 4 y n 1 1 t n 4 y n também conhecida como fórmula de Heun. Calculando temos a tabela: 8. EDO’s 8.2.1 Método de Euler Aprimorado Solução por Euler aprimorado de y 1 t 4 y com y (0) 1 t h=0.05 h=0.025 h=0.01 h=0.001 Exata 0.0 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 0.1 1.5952901 1.6076289 1.6079462 1.6088585 1.6090418 0.2 2.4644587 2.5011159 2.5020618 2.5047827 2.5053299 0.3 3.7390345 3.8207130 3.8228282 3.8289146 3.8301388 0.4 5.6137120 5.7754845 5.7796888 5.7917911 5.7942260 0.5 8.3766865 8.6770692 8.6849039 8.7074637 8.7120041 1.0 60.037126 64.382558 64.497931 64.830722 64.897803 1.5 426.40818 473.55979 474.83402 478.51588 479.25919 2.0 3029.3279 3484.1608 3496.6702 3532.8789 3540.2001 8. EDO’s 8.2.1 Método de Euler Aprimorado O Método de Euler Aprimorado fornece resultados muito melhores do que aqueles de Euler Direto e Inverso. O Método de Euler Aprimorado Adaptativo fornece melhores resultados através da variação no tamanho dos passos. Neste procedimento, variando o tamanho dos passos, mantemos constante o erro de truncamento local da aproximação