Cálculo Numérico Prof. Guilherme Amorim 29/10/2013 Aula 3 – Arredondamento e Operações.
Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.
-
Upload
kevin-cornelio -
Category
Documents
-
view
221 -
download
3
Transcript of Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.
![Page 1: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/1.jpg)
Cálculo Numérico
Prof. Guilherme Amorim07/01/2014
Aula 19 – Interpolação – Parte 1
![Page 2: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/2.jpg)
Pergunta Em um dado problema, foram
encontrados os seguintes pontos:
Como poderíamos estimar os valores de f(x) para x=0,5? E para 1,2? E 1,6?
![Page 3: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/3.jpg)
O que é Interpolação? Dado o tabelamento
Como encontrar o valor de f para algum x entre [x0; xn]? Se x pertence à tabela, OK. E se x não pertence à tabela?
![Page 4: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/4.jpg)
Interpolação
Para resolver essa situação, determinaremos uma função
P que passe exatamente nesses pontos tabelados e a
usaremos para calcular o valor aproximado de f.
![Page 5: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/5.jpg)
Interpolação No nosso contexto, vamos aproximar f
por um polinômio. P terá, portanto, a forma:
![Page 6: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/6.jpg)
Exemplo 1 Dados os pontos tabelados abaixo,
descreva um polinômio que passe exatamente pelos pontos tabelados.
xi 1 2f(xi)
3 4
![Page 7: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/7.jpg)
Exemplo 1 Uma reta é um polinômio?
Sim, então podemos aproximar por uma reta...
![Page 8: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/8.jpg)
Exemplo 1 Há outras formas de aproximar os
pontos tabelados por um polinômio?
Entretanto...
![Page 9: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/9.jpg)
Polinômio interpolador (Def. 5.1)
![Page 10: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/10.jpg)
Voltando ao Exemplo 1 Dados os pontos tabelados abaixo,
descreva um polinômio que passe exatamente pelos pontos tabelados.
A reta y=x+2 é o único polinômio (polinômio interpolador) que passa por (1,3) e (2,4).
xi 1 2f(xi)
3 4
![Page 11: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/11.jpg)
Exercício (Exemplo 5.1)
![Page 12: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/12.jpg)
Exercício (Exemplo 5.1)
![Page 13: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/13.jpg)
Comentários “Essa forma de aproximar funções
(interpolação) só será desejável caso tenhamos certeza sobre a corretude dos valores da tabela, pois, de outra forma, não temos uma boa explicação para as perguntas: Por que a preocupação de passar por pontos duvidosos? Não seria melhor um ajustamento?” [1]
![Page 14: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/14.jpg)
Pergunta? Qual a diferença entre Interpolação e
Ajustamento? Ajustamento
Num ajustamento, nós construímos uma curva que se aproxima (se ajusta) dos pontos.
Podemos extrapolar a análise para além dos extremos. Interpolação
Na interpolação, nós construímos uma curva que passa por todos os pontos.
Não podemos extrapolar a análise para além dos extremos. Ajuda a estimar pontos entre [x0; xn]
![Page 15: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/15.jpg)
Qual a diferença entre Interpolação e Ajustamento?
Fonte: [3]
![Page 16: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/16.jpg)
Voltando ao Exemplo 5.1...
Nós utilizamos sistemas lineares para resolver o problema.
Existe outra forma?
![Page 17: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/17.jpg)
Polinômio Interpolador de Lagrange (Teo. 5.1)
![Page 18: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/18.jpg)
Exemplo Qual o polinômio interpolador de
Lagrange para o problema do Exemplo 5.1?
Fica faltando explicitar L...
![Page 19: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/19.jpg)
Exercício no Quadro Detalhar o valor de Li abaixo
Indicar os valores de L0, L1 e L2 para n=2.
![Page 20: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/20.jpg)
Exemplo 5.1 (Lagrange)
![Page 21: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/21.jpg)
Exemplo 5.1 (Lagrange)
![Page 22: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/22.jpg)
Exemplo 5.1 (Lagrange)
![Page 23: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/23.jpg)
Exercício de sala - Exemplo 5.2
![Page 24: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/24.jpg)
Exercício de sala - Exemplo 5.2
![Page 25: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/25.jpg)
Teorema 5.1 (Prova)
Existência
![Page 26: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/26.jpg)
Teorema 5.1 (Prova) Existência
Tomando os valores de Li do Exemplo 5.2:
Notar que Li(xk) só assume os valores 1 ou 0 .
![Page 27: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/27.jpg)
Teorema 5.1 (Prova) Unicidade
(I) D(x) é um polinômio de grau não superior a n. (II) D(xi) = Pa(xi) – Pb(xi) = f(xi) – f(xi) = 0; i=0, 1, ..., n. Logo, D(x) tem n+1 zeros.
Absurdo! De (I) e (II) temos que D(x) = 0 e, logo, Pa(x) = Pb(x) Ou seja, é ÚNICO!
![Page 28: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/28.jpg)
Exercícios
![Page 29: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/29.jpg)
Bibliografia [1] Silva, Zanoni; Santos, José Dias.
Métodos Numéricos, 3ª Edição. Universitária, Recife, 2010.
[2] Ruggiero, Márcia; Lopes, Vera. Cálculo Numérico – Aspectos Teóricos e Computacionais, 2ª Edição. Pearson. São Paulo, 1996.
[3] Kiusalaas, Jaan; Numerical Methods in Engineering with Python. 2ª edição. 2010.
![Page 30: Cálculo Numérico Prof. Guilherme Amorim 07/01/2014 Aula 19 – Interpolação – Parte 1.](https://reader036.fdocument.pub/reader036/viewer/2022062400/570638461a28abb8238f33d1/html5/thumbnails/30.jpg)