Teoria

Introdução

Fala meus consagrados, vamos começar a entrar em terras turbulentas e desvendar mistérios que nos fazem até tremer as pernas hehehe já sabe do que eu estou falando?

Já usamos métodos pra achar zero reais de funções, aí passamos pra métodos que transformamos dados em polinômios e até mesmo facilitando a resolução de Integrais.

Mas agora... vamos começar a usar algumas técnicas para resolver Equações Diferenciais Ordinárias... isso mesmo, as temidas EDO’s.

Veja bem, teremos métodos e esquemas numérico pra resolver Problemas de Valor Inicial (PVI) pra encontrar uma função y ( x ) , de modo que teremos algo assim:

2 y x - 2 = x ∙ y ' ( x ) y 1 = 3

Aí, vamos aplicar aquela série famosinha, a tal da série de Taylor para y ( x ) no ponto x k .

Mas calma, então preciso saber essa tal de série de Taylor primeiro né? Correto!! Por isso, nesse tópico, vamos aprender como usar a série de Taylor, não teremos resoluções de EDO ainda. A série vai ser muuiito importante na frente, temos que entender ela primeiro.

Então, vem tranquilo. O que é série de Taylor?

Cara, basicamente, a série de Taylor é usada para aproximar funções, tipo essa aqui f x = e - 2 x . Vamos ver como ela faz isso já já.

Talvez você tenha visto essa matéria em outra disciplina, e pode estar fresquinho aí na sua cabeça, pode passar isso aqui:

Se você não viu ouuuuuu se já viu, mas não lembra muito bem é uma ótima chance para revisar haha.

O que acontece é que em muitas faculdades Taylor é de fato uma matéria de numérico (pelo menos uma introdução a Taylor) então esse tópico é pra vocês!

Série de Taylor

Eu aposto que você já ouviu falar alguma coisa sobre essa série, seja em física, EDO, cálculo 3... A série de Taylor é uma série bem famosa, e está por trás de alguns métodos que já vimos em Cálculo Numérico.

Basicamente ela serve para igualar uma função qualquer a uma soma infinita dessa forma:

f x = ∑ n = 0 ∞ f n a x - a n n !

Onde f n ( a ) é a enésima derivada da função no ponto a e ( n + 1 ) é a quantidade de termos da nossa série.

Nesse caso dizemos que essa é a Série de Taylor de f x em torno de x = a . Além disso, toda Série de Taylor tem um “raio de convergência” R em torno do ponto onde ela é centrada a - R , a + R .

Dentro dessa região a série converge para o valor da função. Fora dessa região, dá tudo errado haha. Em alguns casos, no entanto, R = ∞ e a série aproxima perfeitamente a função em qualquer ponto!

Polinômio de Taylor

Show de bola, mas e daí?

Que raios isso tem a ver com Cálculo Numérico? Hmmm... pensa um pouquinho, aqui em numérico nós estamos seeeeempre tentando fazer algumas aproximações aqui e ali (claro, sempre prezando pra não ter erro né!)

Agora, ao invés de pegarmos a série infinita vamos truncar ela em algum ponto obtendo o polinômio de Taylor de grau n .

Dessa maneira, nós podemos aproximar qualquer função por um polinômio! Essa é a graça do Cálculo Numérico, simplificar nossa vida!

É claro que, associado a esse polinômio nós teremos um erro relativo ao truncamento que fizemos. Basicamente, o polinômio de Taylor de grau n vai ficar assim:

f x = f x 0 + f ' x 0 x - x 0 + f ' ' x 0 2 ! x - x 0 2 + … + f n x 0 n ! x - x 0 n = ∑ k = 0 n f k x 0 k ! x - x 0 k  

E o erro de truncamento vai ser dado por:

E x = f n + 1 ξ x n + 1 ! x - x 0 n + 1

Naquele mesmo esquema da maioria dos métodos que vimos até agora, onde a raiz ξ x é um valor que depende de x e se encontra no intervalo entre x e x 0 . Dificilmente nós iremos calcular esse número, o importante aqui é saber que ele se encontra entre x e x 0 .

De forma geral o que fazemos é encontrar o máximo de f n + 1 ξ x nesse intervalo chegando assim a um limitante do erro cometido:

E x ≤ max k ∈ x 0 , x ⁡ f n + 1 k n + 1 ! x - x 0 n + 1

Obs.: Para obtermos um polinômio de Taylor de grau n de alguma função f x essa função deve ter derivadas até a ordem n contínuas.

Exemplo

Então vamos lá! Pra resolver um exercício vamos basicamente calcular alguns polinômios de Taylor .

Bora lá? Vamos calcular o polinômio de Taylor de graus 1 , 2 , 3 e 4 de f x = e - 2 x com x 0 = 0 .

Como e - 2 x tem todas as derivadas contínuas não precisamos nos preocupar com isso. Calculando as derivadas de primeira até quarta, temos que:

f x = e - 2 x

f ' x = - 2 e - 2 x

f ' ' x = 4 e - 2 x

f 3 x = - 8 e - 2 x

f 4 x = 16 e - 2 x

Aplicando x 0 = 0 na função e em todas as derivadas, ficamos com:

f 0 = 1

f ' 0 = - 2

f ' ' 0 = 4

f 3 0 = - 8

f 4 0 = 16

Beleza, até aí tranquilo. Agora vamos relembrar aquela fórmula que vimos lá encima:

f x = f x 0 + f ' x 0 x - x 0 + f ' ' x 0 2 ! x - x 0 2 + … + f n x 0 n ! x - x 0 n = ∑ k = 0 n f k x 0 k ! x - x 0 k  

Com isso vamos calcular os polinômios com n = 1 ; 2 ; 3   e   4 .

Para n = 1 :

f x = f 0 + f ' 0 x - 0 = 1 - 2 x

n = 2 :

f x = f 0 + f ' 0 x - 0 + f ' ' 0 2 x - 0 2 = 1 - 2 x + 2 x 2

n = 3 :

f x = f 0 + f ' 0 x - 0 + f ' ' 0 2 x - 0 2 + f ( 3 ) 0 6 x - 0 3 = 1 - 2 x + 2 x 2 - 4 3 x 3

n = 4 :

f x = f 0 + f ' 0 x - 0 + f ' ' 0 2 x - 0 2 + f ( 3 ) 0 6 x - 0 3 + f ( 4 ) 0 24 x - 0 4 = 1 - 2 x + 2 x 2 - 4 3 x 3 + 2 3 x 4

Tudo beleza? Olha essa imagem pra ver como à medida que aumentamos o n nossa aproximação vai melhorando:

Olhando o gráfico, percebemos que todas as aproximações tem valor exato em 0, mas isso acontece devido ao fato de que nosso polinômio de Taylor foi construído em torno de x 0 = 0 .

Vamos ver o erro de cada uma? O erro é dado por:

E x ≤ max k ∈ x 0 , x ⁡ f n + 1 k n + 1 ! x - x 0 n + 1

Para n = 1 :

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ' ' 1 + 1 ! x = max k ∈ 0 , x ⁡ 4 e - 2 x 2 ! x

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ' ' 1 + 1 ! = 2 x

Para n = 2 :

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ' ' ' 2 + 1 ! x = max k ∈ 0 , x ⁡ - 8 e - 2 x 3 ! x

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ' ' ' 2 + 1 ! = 4 3 x

Para n = 3 :

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ( 4 ) 3 + 1 ! x = max k ∈ 0 , x ⁡ 16 e - 2 x 4 ! x

E x ≤ max k ∈ 0 , x ⁡ e - 2 x 4 3 + 1 ! = 2 x 3

Para n = 4 :

E x ≤ max k ∈ 0 , x ⁡ e - 2 x ( 5 ) 4 + 1 ! x = max k ∈ 0 , x ⁡ - 32 e - 2 x 5 ! x

E x ≤ max k ∈ 0 , x ⁡ e - 2 x 5 4 + 1 ! = 4 15 x

Com isso, podemos ver em todos os limitantes que, conforme x vai se afastando de 0 , o erro vai aumentando. Isso porque os erros são dados em função de x ! Mas vimos agora que quanto maior o n melhor são as aproximações em torno de x 0 = 0 .

Para valores distantes de x 0 = 0 os maiores n não serão necessariamente melhores!

E é só isso que eu tinha pra falar mesmo galera! Vamos fazer uns exercícios que no próximo tópico veremos uma aplicação desse polinômio de Taylor.

Série de Taylor

Polinômio de Taylor

Exemplo