Muitas funções têm suas expressões escritas como um "polinômio infinito", isto é, uma soma infinita de termos do tipo
.
Portanto, uma série (de funções) na forma
é chamada série de potência (em torno de
), onde
é chamado coeficiente da série. Ademais, o
desenvolvimento de funções em séries de potências é aplicado para diversas finalidades, sendo uma delas a resolução de equações diferenciais. A seguir é
apresentada uma descrição baseada em Guidorizzi (2001).
Suponha que a função f tem derivadas de todas as ordens e que possa ser representada por uma série de potências, com raio de convergência R, conforme
representado na Eq. abaixo:
%3D%20%5Csum%5E%7B%5Cinfty%7D_%7Bn%3D0%7D%20c_n%20(x-x_0)%5En%20%2C%20%5Cmid%20x-x_0%5Cmid%20%3CR%20%25%5Clabel%7Bserie%7D)
%3D%20c_0%20%2B%20c_1%20(x-x_0)%2B%20c_2%20(x%2Bx_0)%5E2%2B%20%5Ccdots%20%2B%20c_n%20(x-x_0)%5En%20%2B%20%5Ccdots)
Calculando as derivadas da função f(x) e fazendo
em
, tem-se:
%3D%20%5Cfrac%7Bf%5E%7B(0)%7D(x_0)%7D%7B0!%7D%5C%5C%20c_1%20%3D%20f'(x_0)%3D%20%5Cfrac%7Bf%5E%7B(1)%7D(x_0)%7D%7B1!%7D%20%5C%5C%20c_2%20%3D%20f''(x_0)%3D%20%5Cfrac%7Bf%5E%7B(2)%7D(x_0)%7D%7B2!%7D%5C%5C%20%25%5Chspace%7B1cm%7D%20%5Cvdots%20%5Chspace%7B1.5cm%7D%20%5Cvdots%20%20%5C%5C%20%5Chspace%7B-2cm%7D%5Cvdots%20%5Chspace%7B1.cm%7D%5C%5C%20c_n%20%3D%20f%5En%20(x_0)%3D%20%5Cfrac%7Bf%5E%7B(n)%7D(x_0)%7D%7Bn!%7D)
Logo, se f(x) tem expansão em série de potências centradas em
, então ela deve ser da forma:
%3D%20%5Csum%5E%7B%5Cinfty%7D_%7Bn%3D0%7D%20%5Cfrac%7Bf%5E%7B(n)%7D(x_0)%7D%7Bn!%7D(x-x_0)%5En)
A série da equação é chamada série de Taylor da função f(x) centrada em
.
Como consequência, é possível realizar a aproximação das derivadas usando a expansão de Taylor para funções de uma variável.
Essa aproximação pode ser feita de três formas:
- Avançada;
- Atrasada;
- Centrada.
Realizando a expansão de Taylor para f(x) em torno do ponto
, e considerando
, o que implica em
, tem-se:
%3D%20%5Csum%5E%7B%5Cinfty%7D_%7Bn%3D0%7D%20%5Cfrac%7Bf%5E%7B(n)%7D(x%20_0)%20%7D%7Bn!%7D%20(%5Cunderbrace%7Bx-x_0%7D_%7Bh%7D)%5En%3Df(x_0)%2Bf'(x_0)h%2B%5Cfrac%7Bf''(x_0)h%5E2%7D%7B2!%7D%2B%5Cfrac%7Bf'''(x_0)h%5E3%7D%7B3!%7D%2B%5Ccdots)
%3D%20%5Csum%5E%7B%5Cinfty%7D_%7Bn%3D0%7D%20%5Cfrac%7Bf%5E%7B(n)%7D(x%20_0)%20%7D%7Bn!%7D(%5Cunderbrace%7Bx-x_0%7D_%7Bh%7D)%5En%3Df(x_0)-f'(x_0)h%2B%5Cfrac%7Bf''(x_0)h%5E2%7D%7B2!%7D-%5Cfrac%7Bf'''(x_0)h%5E3%7D%7B3!%7D%2B%5Ccdots)
Tomando como base a Eq.
e colocando
em evidência, obtem-se:
%3D%5Cfrac%7Bf(x_0%2Bh)-f(x_0)%7D%7Bh%7D-%5Cfrac%7Bf''(x_0)h%7D%7B2!%7D-%20%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D-%5Ccdots)
Logo,
%5Ccong%5Cfrac%7Bf(x_0%2Bh)-f(x_0)%7D%7Bh%7D)
é uma aproximação para a derivada de primeira ordem de f(x) em
por diferenças avançadas com erro de truncamento de ordem h, dado por:
%3D-%5Cfrac%7Bf''(x_0)h%7D%7B2!%7D-%20%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D-%5Ccdots)
Agora, tomando como base a Eq.
e colocando
em evidência, obtem-se:
%3D%5Cfrac%7Bf(x_0-h)-f(x_0)%7D%7Bh%7D-%5Cfrac%7Bf''(x_0)h%5E2%7D%7B2!h%7D%2B%20%5Cfrac%7Bf'''(x_0)h%5E3%7D%7B3!h%7D-%5Ccdots)
%3D%5Cfrac%7Bf(x_0)-f(x_0-h)%7D%7Bh%7D%2B%5Cfrac%7Bf''(x_0)h%7D%7B2!%7D-%20%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D%2B%5Ccdots)
Logo,
%5Ccong%5Cfrac%7Bf(x_0)%20-%20f(x_0-h)%7D%7Bh%7D)
é uma aproximação para a derivada de primeira ordem de f(x) em
por diferenças atrasadas com erro de truncamento de ordem h, dado por:
%3D%2B%5Cfrac%7Bf''(x_0)h%7D%7B2!%7D-%20%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D%2B%5Ccdots)
Por fim, tomando como base as Eqs.
e
, tem-se:
-f(x_0%20-h)%3D2f'(x_0)h%2B%5Cfrac%7B2f'''(x_0)h%5E3%7D%7B3!%7D%2B%5Ccdots)
Colocando
em evidência, obtem-se:
%3D%5Cfrac%7Bf(x_0%20%2Bh)-f(x_0%20-h)%7D%7B2h%7D-%5Cfrac%7B2f'''(x_0)h%5E3%7D%7B3!2h%7D%2B%5Ccdots)
-f(x_0%20-h)%7D%7B2h%7D-%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D%2B%5Ccdots)
Logo,
%5Ccong%5Cfrac%7Bf(x_0%20%2Bh)-f(x_0%20-h)%7D%7B2h%7D)
é uma aproximação para a derivada de primeira ordem de f(x) em
por diferenças centradas com erro de truncamento de ordem
, dado por:
%5E2%3D-%5Cfrac%7Bf'''(x_0)h%5E2%7D%7B3!%7D%2B%5Ccdots)