Use o botão abaixo para reportar erros ou dar sugestões.
Cálculo Numérico - Versão Scilab
Nas seções 10.2 e 10.3, construimos dois métodos numéricos para resolver problemas de valor inicial. No Exercício Resolvido 10.3.1, vimos que o erro do método de Euler de do método de Euler melhorado caem quando se reduz o passo , ademais, o erro do método de Euler melhorado cai conforme o quadrado de , enquando o do método de Euler cai conforme . Este fenômeno motiva a definção de ordem de precisão.
Definição 10.6.1. O erro de truncamento local é definido como o erro introduzido em cada passo pelo truncamento da equação diferencial supondo conhecida a solução exata no início do intervalo. Um método numérico é dito ter ordem de precisão se o erro de truncamento local for da ordem de .
Exemplo 10.6.1. O método de Euler tem erro de truncamento local de ordem . Para obter este resultado, observamos via expansão de Taylor que:
| (10.123) |
Se escolhermos nesta expressão e, portanto, , temos:
| (10.124) |
Agora notamos que o termo principal do erro é dado por , como a derivada segunda da solução não depende de , o erro local de truncamento decresce conforme e assim a ordem de precisão do método é 1.
Definição 10.6.2. O erro de truncamento global é definido como erro acumulado ao longo de todos os passo de resolução, supondo a condição inicial exata.
A relação entre o erro de truncamento global e o erro de truncamento local depende da função envolvida. Diante de suficiente regularidade, o erro acumulado é da mesma ordem de grande do erro de truncamento local acumulado ao longo do processo, isto é, pode ser estimado multiplicando o erro local pelo número de passos. Como o número de passos necessários para calcular a solução de um problema de valor inicial no ponto é dado por , temos que a erro de truncamento global é uma ordem inferior ao erro de truncamento local e equivale à ordem de precisão do método.
Usamos também a notação para o erro de truncamento local e para o erro de truncamento global. De forma que, para o método de Euler, temos:
| (10.125) |
Exemplo 10.6.2. Vamos obter o erro de truncamento local do método de Euler melhorado. Partimos da construção do esquema iterativo de Euler melhorado:
Neste ponto, usamos o erro de truncamento do método de trapézios para aproximar a integral envolvida:
Assim, temos que o erro de truncamento local do método de Euler melhorado é e, portanto, um método de ordem 2.
Esta seção carece de exercícios resolvidos. Clique em e inicie a editá-la agora mesmo. Veja outras formas de participar clicando aqui.
Esta seção carece de exercícios. Clique em e inicie a editá-la agora mesmo. Veja outras formas de participar clicando aqui.
E 10.6.1. Aplique o método de Euler e o método de Euler melhorado para resolver o problema de valor inicial dado por
com passo , , , e para obter aproximações para . Compare com a solução exata dada do problema dada por através do erro relativo e observe a ordem de precisão do método.
| h | |||||||
| Euler | 0,4495791 | 0,4660297 | 0,4675999 | 0,4677562 | 0,4677718 | ||
| 9,1e-03 | 8,9e-04 | 8,9e-05 | 8,9e-06 | 8,9e-07 | |||
| Euler mod, | 0,4686037 | 0,4677811 | 0,4677736 | 0,4677735 | 0,4677735 | ||
| 1,8e-03 | 1,6e-05 | 1,6e-07 | 1,6e-09 | 1,6e-11 | |||
A solução exata vale .
E 10.6.2. Resolva o problema de valor inicial dado por
com passo , , , e para obter aproximações para
| h | |||||||
| Euler | 1,1617930 | 1,1395726 | 1,1374484 | 1,1372369 | 1,1372157 | ||
| Euler mod | 1,1365230 | 1,1372075 | 1,1372133 | 1,1372134 | 1,1372134 | ||
Este texto é disponibilizado nos termos da licença Creative Commons Atribuição-CompartilhaIgual 3.0 Não Adaptada (CC-BY-SA 3.0). Página gerada em 19/8/2020 às 17:36:32.