1) Introdução ao conceito de interpolação polinomial (Geral), apresentação
superficial da formula da Interpolação polinomial de Lagrange
- Lucas Rocha
Inicialmente, o foco é obter uma função polinomial cuja curva do gráfico passe por
todos os pontos de um determinado conjunto iniciado em (x 0, y0), . . ., (xn, yn). Ou
seja, encontrar um polinômio P(x) =a0+a1x+. . .+anxn tal que P(x0) =y0, P(x1) =y1, . . .,
P(xn) =yn. O Teorema a seguir comprova a sua existência.
Teorema: Para um conjunto de n pontos de abscissas diferentes, há uma única
função polinomial de grau até n−1 cujo gráfico passa por todos os pontos.
Esse processo de obtenção do polinômio chama-se Interpolação Polinomial. Há
algumas formas de se obter esse polinômio. O primeiro a ser levantado aqui é a
Interpolação Polinomial de Lagrange.
Considere os pontos (x0, y0), . . . ,(xn, yn) e os polinômios a seguir:
L0(x) = (x−x1)(x−x2)(x−x3) . . . (x−xn) / (x0−x1)(x0−x2)(x0−x3) . . . (x0−xn)
L1(x) = (x−x0)(x−x2)(x−x3) . . . (x−xn) / (x1−x0)(x1−x2)(x1−x3) . . . (x1−xn)
L2(x) = (x−x0)(x−x1)(x−x3) . . . (x−xn) / (x2−x0)(x2−x1)(x2−x3) . . . (x2−xn)
...
Ln(x) = (x−x0)(x−x1)(x−x2). . .(x−xn−1) / (xn−x0)(xn−x1)(xn−x2). . .(xn−xn−1)
É possível perceber que Lk(xj) ={1, k = j | 0, k ≠ j}.
Temos que o polinômio interpolador de Lagrange é dado pela seguinte fórmula:
P(x) = a0L0(x) + a1L1(x) +. . .+anLn(x).
Basta encontrar os valores de a0, a1, . . ., an.
Uma vez que o polinômio é aplicado em cada valor xi, o resultado é o seguinte:
P(x0) =a0L0(x0) +a1L1(x0) +. . .+anLn(x0) =y0
P(x1) =a0L0(x1) +a1L1(x1) +. . .+anLn(x1) =y1
P(x2) =a0L0(x2) +a1L1(x2) +. . .+anLn(x2) =y2
P(xn) =a0L0(xn) +a1L1(xn) +. . .+anLn(xn) =yn
A partir das igualdades, encontramos: a0=y0, a1=y1, a2=y2, . . ., an=yn.
Portanto, o polinômio de Lagrange pode ser escrito como:
P(x) =y0L0(x) + y1L1(x) +. . .+ ynLn(x)
Substituindo L0,L1,L2,. . . , na igualdade acima:
P(x) =y0 [(x−x1)(x−x2)(x−x3). . .(x−xn) / (x0−x1)(x0−x2)(x0−x3). . .(x0−xn)] + y1 [(x−x0)(x−x2)
(x−x3). . .(x−xn) / (x1−x0)(x1−x2)(x1−x3). . .(x1−xn)] +. . .+ yn [(x−x0)(x−x1)(x−x2). . .(x−xn−1) /
(xn−x0)(xn−x1)(xn−x2). . .(xn−xn−1)]
2) Demonstração da aplicação da Interpolação polinomial de Lagrange
- Gabriel
Interpolação Polinomial de Newton?
Um outro modo para encontrar o mesmo polinômio interpolador de Lagrange que
passa por todos os pontos de um conjunto é através da Interpolação Polinomial de
Newton. O resultado é o mesmo do polinômio de Lagrange, porém os cálculos são
otimizados em comparação ao primeiro método, principalmente para um elevado
número de pontos.
O primeiro passo é construir uma tabela a partir dos pontos (x1, y1), (x2, y2), (x3, y3), . .
.,(xn, yn). No esquema a seguir são considerados quatro pontos, mas é possível
construir para um número maior de pontos seguindo o mesmo raciocínio.
- Forma de Newton do Polinômio de Interpolação
Utilizando os valores obtidos no topo da tabela, constrói-se o polinômio da seguinte
forma: P(x) =y1+ ∆01(x−x1) + ∆11(x−x1) (x−x2) + ∆21(x−x1) (x−x2) (x−x3) +. . ..
Exemplo:
Encontre o polinômio interpolador de Newton para os pontos (−1,3), (0,5), (2,15). O
primeiro passo é montar a tabela de diferenças divididas, colocando na primeira
coluna os valores das abscissas e na segunda coluna os valores correspondentes
desordenados. As demais colunas seguem a construção anterior.
O próximo passo é aplicando a fórmula, tomando os valores que estão no topo da
tabela. O resultado obtido é mostrado a seguir:
Interpolação Polinomial de Gregory-Newton
Um caso particular do polinômio interpolador de Newton ocorre quando a distância
entre as abscissas dos pontos consecutivos é a mesma. Podemos dizer em outras
palavras que quando os pontos estão igualmente espaçados sobre o eixo x. Neste
caso, podemos utilizar uma tabela como por exemplo:
- Tabelas de diferenças finitas:
Considere os pontos (x1, y1), (x2, y2), (x3, y3), . . .,(xn, yn)tais que
x2−x1=x3−x2=x4−x3=. . .=h. No esquema a seguir são considerados quatro
pontos, entretanto é possível construir para um número maior de pontos seguindo o
mesmo raciocínio.
Polinômio Interpolador de Gregory-Newton
Utilizando os valores obtidos no topo da tabela, constrói-se o polinômio da seguinte
forma:
Exemplo: Encontre o polinômio interpolador de Gregory-Newton para os pontos
(0,2), (2,−2) e (4,2).
A primeira iniciativa é montar a tabela de diferenças finitas, colocando na primeira
coluna os valores das abscissas e na segunda coluna os valores correspondentes
desordenados. As demais colunas seguem a construção anterior.
O próximo passo é aplicando a fórmula, tomando os valores que estão no topo da
tabela. O resultado obtido é mostrado a seguir:
Links:
- Material utilizado na aula virtual;
- [Link]
- [Link]
2_Forma_de_NewtonGregory.pdf
-
- [Link]
- [Link]
[Link]
5) Demonstração da aplicação da Interpolação polinomial de Gregory-
Newton
-Daniel Tawi
Aplicação do Polinômio de Gregori Newton em valores igualmente espaçados:
Determinar o Polinômio interpolador com os seguintes dados:
Utilizamos a formula geral:
P2 (x)= f(x0) + (x – x0) Δf(x0)/h + (x – x0)(x – x1)Δ²f(x0)/2h²
-> Substituindo:
P2 (x)= 0 + (x – 0) 0,5/1 + (x – 0)(x – 1)(-0,25)/2*(1²)
P2 (x)= 0,5x + [x² - x]*(-0,125)
P2 (x)= 0,5x -0,125x² + 0,25x
P2 (x)= -0,125x² + 0,625x