Interpolação polinomial por partes
A maneira de interpolar muitos pontos, mas usando polinômios de grau baixo
Vamos lá?
O fenômeno de Runge foi um balde de água fria na nossa intenção de construir melhores interpolações através de considerar cada vez mais pontos de interpolação. Nesta aula vamos ver como contornar essa dificuldade.
- 1
- 2
Se temos 21 pontos de interpolação, e consideramos 10 subintervalos, podemos dizer que,..
- A. em cada subintervalo temos 3 pontos que são interpolados por um polinômio de grau 2.
- B. em cada subintervalo temos 2 pontos que são interpolados por um polinômio de grau 2.
- C. em cada subintervalo temos 3 pontos que são interpolados por um polinômio de grau 3.
- D. nem todos os intervalos terão a mesma quantidade de pontos de interpolação.
Só há duas maneiras de evitar o fenômeno de Runge, ou usamos os pontos de Chebyshev como nós de interpolação (veja a aula passada) ou nos restringirmos a polinômio de grau baixo. Como usar os pontos de Chebyshev muitas vezes não é possível, temos que ficar com a segunda opção. Mas como interpolar uma tabela com muitos pontos e ainda assim considerar polinômios de grau baixo? A solução é dividir o problema de interpolação completo em muitos pequenos problemas de interpolação menores.
Na figura abaixo vemos o resultado de interpolar os sete pontos, considerando-os dois a dois, três a três ou quatro a quatro (mantenha o mouse sobre a figura para pausar a animação).
Na figura acima, quando consideramos pontos dois a dois, foram realizadas 6 interpolações com polinômios de grau 1. No segundo caso, realizamos três interpolações, nos subintervalos [0,2], [2,4] e [4,6], em cada uma deles, usando polinômios de grau 2. No terceiro caso, foram apenas duas interpolações, nos subintervalos [0,3] e [3,6], em cada um deles usando polinômios de grau 3.
Como a função interpoladora (linha contínua em vermelho, no gráfico) é construída pela colagem de polinômio, ela não é em si um polinômio. Essa função é dita um polinômio por partes. Em particular, enquanto polinômios são funções infinitamente diferenciáveis, polinômios por partes são funções apenas contínuas, por conta dessa colagem nos extremos dos subintervalos. Claro que nos pontos interiores aos subintervalos a função segue sendo infinitamente diferenciável.
Ainda sobre o exemplo acima, fica evidente que a medida que o grau do polinômio aumentou, a aproximação ficou melhor. Porém, não podemos nos esquecer do fenômeno de Runge. Logo, o usual é não ir muito além de grau 3.
Vamos considerar o caso particular dos polinômios por partes de grau 1, também conhecidos por splines lineares. Suponha que uma função
Da fórmula acima vemos que a medida que
É bem verdade que para que o erro reduza a um nível aceitável, a quantidade de pontos de interpolação pode ser grande. No exemplo do vídeo, a função
Talvez você esteja com a impressão de que splines lineares sejam só uma maneira pedante de dizer ligue os pontos. Talvez seja mesmo. Mas nem por isso devemos menosprezá-las. De fato, todos os gráficos que você já viu neste curso foram construídos com splines lineares! Por exemplo, volte ao topo desta página e veja novamente o primeiro gráfico. O gráfico da função suave (tracejado em preto) não é realmente suave. Isto é apenas uma impressão visual. Se você reparar bem, cada tracinho preto é um segmento de reta, mas por serem pequenos você ficou com a impressão de ver o gráfico de uma função suave.
Considere os comandos abaixo, interpretados no Octave.
n = 101; x = linspace(0,3,n); y = sin(x); plot(x,y)
Ao executá-los, você terá a impressão de ver o gráfico da função seno, mas na verdade estará vendo o gráfico da spline linear que interpola a função seno em
Será que a ilusão de suavidade é sempre suficiente? Nem sempre. Há situações em que precisamos de uma função interpoladora que realmente seja suave e capaz de lidar com muitos pontos, sem correr o risco de observar o fenômeno de Runge. Neste caso, outra abordagem é necessária. Este será tema outra outra aula.
1. Se
2. Em quantos pontos a função cosseno deve ser amostrada para que sua spline linear interpolante tenham erro máximo de
Lembre que a função cosseno é periódica, então basta representá-la dentro de um intervalo finito. Qual intervalo é suficiente?
Se
3. Considere o conjunto
- Mostre que
é um espaço vetorial. - Seja
definido por se e para Mostre que é uma base para o espaço