Interpolação

PROF. HERON JR.

Objetivo

Interpolar uma função f(x) consiste em aproximar essa função por uma outra função g(x), escolhida entre uma classe de funções definidas (polinômios). g(x) é usada em substituição à função f.

Problemática

Essa necessidade de efetuar esta substituição surge quando: Quando são somente conhecidos os valores

numéricos da função para um conjunto de pontos e é necessário calcular o valor de um ponto no tabelado.

Quando a expressão da função é complicada de mais para ser integrada ou diferenciada.

Em equação

Consideremos n+1 valores distintas: x0, ..., xn (nós da interpolação) e os valores de f(x) nesses pontos: f(x0), ..., f(xn).

Queremos determinar a função g(x) tal que:g(x0)=f(x0)....g(xn)=f(xn)

Graficamente

Classe de funções

Em nosso caso, consideramos a função g(x) com um elemento da classe de funções polinomiais.

Tentaremos aproximar a função f(x) a partir de um conjunto de valores com uma função do tipo: a0+a1x+...+anxn

Interpolação polinomial

Dados os n+1 pontos (x0,f(x0)), ..., (xn,f(xn)), queremos aproximar f(x0) por um polinômio pn(x) de grau menor ou igual a n: f(xk)=pn(xk) ; k=0,1,...n

Interpolação polinomial

Considerando que p o polinômio escreve-se pn(x)= a0+a1x+...+anxn , a condição f(xk)=pn(xk) ; k=0,1,...,n produz o sistema seguinte de n+1 equações , n+1 variáveis:

0 1 0 0 0

0 1 1 1 1

0 1

... ( )

... ( ).........

... ( )

nn

nn

nn n n n

a a x a x f x

a a x a x f x

a a x a x f x

Interpolação polinomial: matriz

A matriz do sistema é:

Essa matriz é uma matriz de Vandermonde, desde que x0, ..., xn são pontos distintos, temos det A¹0. Então o sistema admite uma solução única.

0 0

1 1

1 ...1 ...1 ... ... ...1 ...

n

n

nn n

x xx x

A

x x

Prova

Podemos proceder da forma seguinte: O determinante pode ser considerado como um polinômio em x0:

E um polinômio de grau n com n raízes: x1 a xn, ele pode ser escrito aP(xi-x0); i¹0

0 0

21 10 1 0 2 0 0

1 ...1 ...

...1 ... ... ...1 ...

n

nn

n

nn n

x xx x

x x x

x x

a a a a

Determinante de Vandermonde

O determinante da matriz de Vandermonde pode ser escrito da forma seguinte:

0 0

1 1

0

1 ...1 ...

( )1 ... ... ...1 ...

n

n

j ii j n

nn n

x xx x

x x

x x

Interpolação polinomial: teorema

Em outros termos podemos dizer que:

Existe um único polinômio pn(x) de grau n tal que pn(xk)=f(xk), k=0,1,...,n desde que xi¹xj por j¹k.

Obter pn(x)

Para obter o polinômio pn(x), existem diversos métodos, o mais direto sendo a resolução do sistema linear.

A escolha do método depende de várias condições: a estabilidade do sistema, performance computacional, ...

Resolução do sistema

Vamos encontrar o polinômio de grau 2 que interpola os pontos da tabela:

Considerando p2(x)=a0+a1x+a2x2. Temos o sistema:

x -1 0 2f(x) 4 1 -1

0 1 2

0 0 1 2

0 1 2

a -a +a =47 2a =1 a 1, a , a3 3

a +2a +4a =-1

Condicionamento

A determinação dos coeficientes pela resolução do sistema é um processo simples, mas o sistema pode ser mal condicionado e sua resolução com numeração a ponto flutuante produzir resultados errados.

Existem outros métodos para determinar os polinômios de interpolação. Como existe uma solução única, qualquer método que determina uma solução, determina a solução única.

Forma de Lagrange

Considerando os n+1 pontos (x0,y0=f(x0)), ..., (xn,yn= f(xn)) e o polinômio interpolador pn(x). Lagrange propôs de representar o polinômio pn(x) da forma:pn(x)=y0L0(x)+..+ynLn(x), onde Lk(x) são polinômios de grau n e a condição pn(xi)=yi, i=0,...,n seja satisfeita.

Forma de Lagrange

A melhor forma de ter a condição: pn(xi)=yi

i=0,...,n, é impor:

Por isso, consideramos:

k i

1L (x )=

0se k ise k i

¹

0 1 1 1

0 1 1 1

( )( )...( )( )...( )( )

( )( )...( )( )...( )k k n

kk k k k k k k n

x x x x x x x x x xL x

x x x x x x x x x x

Forma de Lagrange

O numerador de Lk(x) é um produto de n fator em x, então Lk(x) é de grau n.

Podemos verificar também que:

A forma de Lagrange para o polinômio interpolador é:

k i

1L (x )=

0se k ise k i

¹

0,

0

0,

( )( ) ( ), ( )

( )

n

jnj j k

n k k k nk

k jj j k

x xp x y L x L x

x x

¹

¹

Interpolação linear

Interpolação de dois pontos (x0,y0=f(x0)) e (x1,y1=f(x1)).

Usando a forma de Lagrange, temos:0 1 0 0 11

0 10 1 1 0 1 0

( ) ( ) ( )( )( )( ) ( ) ( )n

x x x x y x x yx xp x y yx x x x x x

Exemplo

Seja a tabela:

Temos:

x -1 0 2f(x) 4 1 -1

2 2 222 2 7 2( ) 4 1

3 2 6 3 3nx x x x x xp x x x

21 2

00 1 0 2

20 2

11 0 1 2

20 1

22 0 2 1

( )( ) ( 0)( 2) 2( )( )( ) ( 1 0)( 1 2) 3

( )( ) ( 1)( 2) 2( )( )( ) (0 1)(0 2) 2

( )( ) ( 1)( 0)( )( )( ) (2 1)(2 0) 6

x x x x x x x xL xx x x x

x x x x x x x xL xx x x x

x x x x x x x xL xx x x x

Forma de Newton

Considerando os n+1 pontos (x0,f(x0)), ..., (xn,f(xn)) e o polinômio interpolador pn(x). Newton propôs de representar o polinômio pn(x) da forma:

pn(x)=d0+d1(x-x0)+d2(x-x0)(x-x1)+...+dn(x-x0)...(x-xn-

1)

Os coeficientes dk, k=0,...,n são diferenças divididas de ordem k entre os pontos (xj,f(xj)), j=0,...,k

Operador diferenças divididas

f(x) é uma função tabelada em x0,...,xn. Os operadores de diferenças divididas são

definidos por:0 0

1 00 1

1 0

1 2 0 10 1 2

2 0

1 0 10

0

[ ] ( ) 0[ ] [ ][ , ] 1

[ , ] [ , ][ , , ] 2

[ ,..., ] [ ,..., ][ ,..., ] n n

nn

f x f x ordemf x f xf x x ordemx xf x x f x x

f x x x ordemx x

f x x f x xf x x ordemn

x x

Operador diferenças divididas

x Ordem 0 Ordem 1 Ordem 2 ... Ordem n

x0 f[x0]

f[x0,x1]

x1 f[x1] f[x0,x1,x2]

f[x1,x2]

x2 f[x2] f[x1,x2,x3]

f[x0,...,xn]

f[xn-2, xn-1, xn]

.... f[xn-1, xn]

xn f[xn]

Operador diferenças divididas

Exemplo:x Ordem 0 Ordem 1 Ordem 2 Ordem 3 Ordem 4

-1 1

0

0 1 -1/2

-1 1/6

1 0 0 -1/24

-1 0

2 -1 0

-1

3 -2

x f(x)-1 10 11 02 -13 -2

Forma de Newton

Podemos provar que as diferenças divididas satisfazem a propriedade seguinte:

Onde j0, ..., jk é qualquer permutação de 0, ..., k.00[ ,..., ] [ ,..., ]

kk j jf x x f x x

Forma de Newton

Forma de Newton para o polinômio interpolador: Seja uma função f(x) contínua e com tantas

derivadas contínuas necessárias num intervalo [a,b].

Sejam a=x0<x1<...<xn=b Vamos construir o polinômio pn(x) que interpola

f(x) em x0, ..., xn, construindo sucessivamente os polinômios pk(x), k=0,...,n

Forma de Newton

Para xÎ[a,b], x¹x0 Temos:

Podemos notar que E0(x) é o erro cometido aproximando f(x) por p0(x)

00

0

0 0 0

0 0

( ) ( )[ , ]

( ) ( ) ( ) [ , ]( ) ( )

f x f xf x x

x xf x f x x x f x x

p x E x

Forma de Newton

0 1 00 1 1 0

1

0 0 1 0 0 1 0 1

1 1

[ , ] [ , ][ , , ] [ , , ]

( ) ( ) ( ) [ , ] ( )( ) [ , , ]( ) ( )

f x x f x xf x x x f x x xx x

f x f x x x f x x x x x x f x x xP x E x

Forma de Newton

1 0 2 1 00 1 2 2 1 0

2

0 0 1 0 0 1 0 1

2

0 1 2 0 1 2

2

[ , , ] [ , , ][ , , , ] [ , , , ]

( ) ( ) ( ) [ , ] ( )( ) [ , , ]( )

( )( )( ) [ , , , ]( )

f x x x f x x xf x x x x f x x x xx x

f x f x x x f x x x x x x f x x xP x

x x x x x x f x x x xE x

Forma de Newton

Continuando assim para todos pk(x), temospn(x)=f(x0)+(x-x0)f[x0,x1]+...

+(x-x0)..(x-xn-1) f[x0,...,xn]

O erro é dado por:En(x)=(x-x0)..(x-xn)f[x0,...,xn,x]

Forma de Newton

Considerando a tabela: x -1 0 2f(x) 4 1 -1

x Ord 0 Ord 1 Ord 2-1 4

-30 1 2/3

-12 -1

23

2 7( ) 13 3

p x x x

Estudo do erro

A aproximar a função f(x) por um polinômio, comete-se um erro:

En(x)=f(x)-pn(x)

Estudo do erro

Teorema:Sejam x0<...<xn, seja f(x) com derivadas até

ordem (n+1) para x no intervalo [x0,xn].Em qualquer ponto x do intervalo [x0,xn], o

erro é dado por:( 1)

0 0( )( ) ( )....( ) , [ , ]

( 1)!

nx

n n x nfE x x x x x onde x xn

Î

Estudo do erro

Do teorema precedente, podemos deduzir que:

Dois corolários: Se f(n+1)(x) é contínua em [x0,xn],

Se além disso, x1-x0=x2-x1=...=xn-xn-1=h

( 1)

0 1 0 0( )[ , ,..., , ] , [ , ], [ , ]

( 1)!

nx

n n nff x x x x x x x x x xn

Î Î

( 1)10 1 [ 0, ]

( ) ( )....( ) , ( ( ) )( 1)!

nnn n n x x xn

ME x x x x x M max f x

n

Î

1

1( )4( 1)

n

n nhE x Mn

Estudo do erro

Se a função é dada na forma de uma tabela, só podemos estimar o valor absoluto do erro.

Mas a tabela de diferencias divididas é construída até ordem n+1, podemos usar o maior valor destas diferenças como aproximação para:

Nesse caso, o valor do erro pode ser majorado com:

1

( 1)!nM

n

0 1( ) ( )...( ) max( )n n diferencias divididas de ordem nE x x x x x

Interpolação inversa

Trata-se de, conhecendo um valor y de Î(f(x0),f(xn)), aproximar um valor de x tal que f(x)=y. Uma solução consiste em interpolar f(x) é em seguida

resolver a equação f(x)=y. No caso de graus elevados (>2), a resolução da equação pode ser difícil e não temos avaliação do erro cometido.

Uma outra solução consiste em efetuar uma interpolação inversa, ou seja determinar um polinômio interpolador de f-1(x). Com a interpolação inversa, podemos calcular uma avaliação do erro cometido.

A interpolação inversa só poder ser feita com uma função monótona.

Grau do polinômio

Trata-se de determinar o grau do polinômio para interpolar uma função em um ponto: Deve-se construir a tabela de diferenças

divididas. Se na vizinhança do ponto de interesse, as

diferenças divididas de ordem k são praticamente constantes, podemos concluir que um polinômio de grau k é suficiente.