Interpolação PROF. HERON JR.. Objetivo Interpolar uma função f(x) consiste em aproximar essa...
-
Upload
madalena-palha-wagner -
Category
Documents
-
view
214 -
download
0
Transcript of Interpolação PROF. HERON JR.. Objetivo Interpolar uma função f(x) consiste em aproximar essa...
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.