Método de Euler

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

1

UNIVERSIDADE FEDERAL DE PERNAMBUCO


CENTRO ACADÊMICO DO AGRESTE
NÚCLEO DE TECNOLOGIA ENGENHARIA CIVIL

PEDRO VÍTOR FLORENTINO DE BRITO

MÉTODO DE EULER

Caruaru

2022
2

1 OS ALGORITMOS

1.1 Algoritmo do Método de Euler

Para o Método de Euler, o seguinte algoritmo foi desenvolvido no software


MATLAB:

clear
clc
n = 5;
dt = 0.2;
y = zeros(1,n);
t = zeros(1,n);
y(1) = 2; t(1) = 0;
for i = 1:n-1
y(i+1) = y(i)+dt*(6*y(i)+t(i).^2);
t(i+1) = t(i)+dt;
end

Para critério de testes, foi utilizada a EDO y’=6y+t2, com 5 interações e um


passo de 0.2. Como a solução analítica da EDO adotada não era conhecida,
foram feitas mudanças nos parâmetros do código (passo e interações) com o
intuito de aferir o erro do código.

1.2 Algoritmo do Método de Euler Melhorado

Já para o Método de Euler Melhorado, foi desenvolvido o seguinte código,


também no MATLAB:

clear
clc
x0 = 0; y0 = 2; h = 0.2;
t(1) = x0; y(1) = y0;
for i=1:4
k1 = (y(i));
k2 = (y(i)+h*(k1));
y(i+1) = y(i)+(h/2)*(k1+k2);
t(i+1) = h+t(i);
end
for i=1:5
yp(i) = 2+3*y(i)*t(i)+(t(i)^3)/3;
end
3

Da mesma forma que o método anterior, o método melhorado segue


passo de 0.2 com 5 interações. Semelhantemente, como também não se
conhece a solução analítica da EDO adotada neste caso (y’=3y+t^2), foi feita
uma comparação com o método tendo passo 0.1 e 10 interações, de modo a
permitir aferição do erro.

2 RESULTADOS

2.1 Sobre o erro

Para determinar a precisão do algoritmo do Método de Euler, faz-se


necessário determinar o erro, foi utilizada a seguinte equação:
|𝑆𝑎𝑛 − 𝑆𝑎𝑝|
𝐸𝑟𝑟𝑜𝑟 =
𝑆𝑎𝑛
O resultado foi, então, convertido para percentual. O percentual médio de
erro de cada método foi determinado pela média dos erros de cada passo.
Como as soluções analíticas das EDO utilizadas não eram conhecidas, a
solução encontrada por cada método utilizando passo = 0.1 e 10 interações foi
adotada como solução analítica para viabilizar os cálculos.

2.2 Resultados do Método de Euler

Após a obtenção dos resultados pelo algoritmo desenvolvido, a seguinte


tabela foi elaborada
Passo Aproximação Analítica Erro (percentual)
0.0 2,0000 2,0000 0,000%
0.1 3,2000 -
0.2 4,4000 5,1210 14,079%
0.3 8,1976 -
0.4 9,6880 13,1252 26,188%
0.5 21,0163 -
0.6 21,3456 33,6510 36,568%
0.7 53,8776 -
0.8 47,0323 86,2532 45,472%
0.9 138,0691 -
0.10 -
Erro Médio: 24,461%

O desvio resultante do erro pode ser visualizado no gráfico abaixo:


4

Solução aproximada x Solução Analítica


160,0000
140,0000
120,0000
100,0000
80,0000
60,0000
40,0000
20,0000
0,0000
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.10

Analítica Aproximação

2.3 Resultados do Método de Euler Melhorado

Do mesmo modo que foi feito no método de Euler, aqui também foi
elaborada uma tabela e um gráfico a partir dos resultados obtidos, de modo a
facilitar a visualização e análise dos dados.

Passo Aproximação Analítica Erro (percentual)


0.0 2,0000 2,0000 0,000%
0.1 2,6633 -
0.2 3,4667 3,4679 0,035%
0.3 4,4376 -
0.4 5,5935 5,5995 0,107%
0.5 6,9840 -
0.6 8,6091 8,6255 0,190%
0.7 10,5629 -
0.8 12,8043 12,8401 0,279%
0.9 15,5064 -
0.10 -
Erro Médio: 0,122%
5

Aproximação x Analítica
20,0000

15,0000

10,0000

5,0000

0,0000
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.10

Analítica Aproximação

3 CONCLUSÃO

Como foi possível observar, o Método de Euler Melhorado apresenta


uma precisão muito maior que o método tradicional. Com ele, os resultados
obtidos tiveram divergência mínima, quase desprezível, dos valores de
referência (a solução analítica adotada), atestando a alta confiabilidade do
algoritmo.
O Método de Euler tradicional, entretanto, também tem seu valor. Ainda
que tenha apresentado uma taxa de erro que aumenta em velocidade
relativamente alta, o algoritmo ainda pode ser utilizado para cálculos mais
simples e rápidos, que não demandem tanta precisão.
Além disso, é importante relembrar que sem o método original, os
demais métodos mais precisos jamais teriam sidos desenvolvidos, o que
reforça a construção do conhecimento como algo sempre contínuo.

Você também pode gostar