Métodos numéricos/Equações diferenciais ordinárias: diferenças entre revisões
| [edição não verificada] | [edição não verificada] |
Conteúdo apagado Conteúdo adicionado
Linha 5:
== Problema de condições iniciais ==
Um problema de valor inicial (PVI) nada mais é que uma EDO acompanhada de uma condição inicial. Condição inicial é um ponto x (numa função f(x)) que pertence ao domínio da função f(x) que é solução do PVI. <br>
Por exemplo: se conhecemos a função velocidade <math> v(t)</math> podemos perguntar qual será a posição <math> x(t) </math> de uma partícula no instante <math> t </math>. A equação diferencial a resolver será:<br>
<math> dx/dt = v(t) </math>, <br>
cuja solução é <math> x(t) = x(0)+<int><sub>0</sub><sup>t</sup>v(t)dt </math>.
A posição inicial é necessãria para definir unicamente a solução do problema.
=== Método de Euler ===
A fórmula que rege esse método é a seguinte<br> <math> y_{n+1} = y_n + h*f(y_n,t_n)</math> (1.1),<br> onde y é igual a f(x) e h é um valor denominado passo; e f(y<sub>n</sub>, t<sub>n</sub>) = y'.<br>
Torne a EDO do PVI em função de y<sub>n</sub> (observe que este é um termo da fórmula acima).Determine o número de passos pela fórmula 1.2. Feito isto podemos iniciar os cálculos para determinar o valor do PVI para determinado valor de x.. A fórmula do método será repetida por várias vezes, portanto por ser um método massivo, seria interessante a utilização de planilhas de cálculo. Para utilizarmos 1.1 devemos primeiramente calcular os valores de yn e y´n. Para a primeira linha do cálculo temos que yn = (y da condição inicial) e y´n é calculado a partir da EDO do PVI para os valores x e y da condição inicial.
Calculados estes dois termos, devemos substituí-los em 1.1, tendo assim efetuado a primeira linha do cálculo, o próximo passo é o cálculo da primeira iteração. Para a primeira iteração temos que<math> y_n=(y_n + 1) </math> calculado na linha anterior e y´n é calculado a partir da EDO do PVI para os valores x + h e (yn+1) da linha de comando anterior; calculados estes dois termos, novamente devemos substituídos em 1.1 e assim terminamos a primeira interação. A segunda iteração tem uma sistemática idêntica a da primeira iteração. Assim iterações sucessivas são executadas até que se atinja a iteração limite calculada anteriormente pela fórmula 1.2; alcançado este valor limite de iteração, o y desejado será o valor da fórmula 1.1 da linha anterior.
| |||