UFMT – ICET MATEMÁTICA (2022/2) – Prof. Geraldo L. Diniz, Cáculo numérico – Turmas BF e T.

Respostas da lista II - interpolação

3. Calcular o valor aproximado do número de habitantes de Cuiabá, no ano de 1995, sabendo que em 1970 haviam
150 mil habitantes, em 1980 haviam 250 mil habitantes, em 1990 haviam 500 mil habitantes e em 2000 eram 600
mil habitantes. Use um polinômio de Lagrange, para fazer a estimativa populacional.
habitantes
Solução: Para reduzir a ordem de grandeza dos números, fazendo t = Ano − 1970 e p = , se obtém a
1000
tabela 1 a seguir.
Com base na tabela 1, o polinômio interpolador de Lagrange é dado por:
(t−10)(t−20)(t−30) (t−0)(t−20)(t−30)
PL (t) = 150 (0−10)(0−20)(0−30) + 250 (10−0)(10−20)(10−30) + ···
(t−0)(t−10)(t−30) (t−0)(t−10)(t−20)
Tabela 1: Dados reduzidos do problema. +500 (20−0)(20−10)(20−30) + 600 (30−0)(30−10)(30−20) ⇒
3 2
t 0 10 20 30 PL (t) = −0.05t + 2.25t − 7.5t + 150.
p 150 250 500 600 Daı́, para o ano de 1995 ⇒ t = 25, logo PL (25) = 587.5. Portanto, em
1995 haveriam 587.500 habitantes.
4. A tabela 2 apresenta a demanda máxima diária de energia elétrica numa cidade.

Tabela 2: Demanda de energia elétrica da cidade.

Data 21 de janeiro 31 de janeiro 10 de fevereiro 20 de fevereiro
Pico de demanda (Mw) 10 15 20 13

Faça uma interpolação com esta tabela para as opções a seguir, ache a demanda no dia 12 de fevereiro e no item c)
encontre a data e pico de máxima demanda.
Solução:
Para facilitar o tratamento dos dados, convertendo as datas para o tempo t em dias e considerando a data inicial
21 de janeiro como sendo t = 0, então se obtém a nova tabela 3.

Tabela 3: Demanda de energia elétrica da cidade.

t (dias) 0 10 20 30
Pico de demanda (Mw) 10 15 20 13

(a) polinômio linear; com base na tabela 3, então o polinômio linear que interpola a data indicada (12 de fevereiro
t − 30 t − 20
⇒ t = 22) é dado por P1 (t) = 20 + 13 ⇒ P1 (t) = −0.7t + 34. Logo, P1 (22) = 18.6 Mw.
20 − 30 30 − 20
(b) polinômio quadrático; para o polinômio quadrático, são necessários 3 pares de pontos, no caso seriam os 3
últimos pares de pontos, para a tabela 3 resulta:
(t − 20)(t − 30) (t − 10)(t − 30) (t − 10)(t − 20)
P2 (t) = 15 + 20 + 13 ⇒
(10 − 20)(10 − 30) (20 − 10)(20 − 30) (30 − 10)(30 − 20)
P2 (t) = −0.06t2 + 2.3t − 2 ⇒ P2 (22) = 19.56 Mw.
(c) polinômio de Lagrange; agora usando os 4 pares de pontos da tabela 3 o polinômio de Lagrange é dado por
(t−10)(t−20)(t−30) (t−0)(t−20)(t−30) (t−0)(t−10)(t−30) (t−0)(t−10)(t−20)
PL (t) = 10 (0−10)(0−20)(0−30) + 15 (10−0)(10−20)(10−30) + 20 (20−0)(20−10)(20−30) + 13 (30−0)(30−10)(30−20) ⇒
PL (t) = −0.002t3 + 0.06t2 + 0.1t + 10 ⇒ PL (22) = 19.94 Mw.
Agora para achar a data e o pico de demanda, calculando a derivada primeira de PL (t) se tem
PL′ (t) = −0.006t2 + 0.12t + 0.1 e resolvendo a equação quadrática −0.006t2 + 0.12t + 0.1 = 0 ⇒ t1 = −0.801 ou
t2 = 20.8. Daı́, tomando a derivada segunda calculada nos pontos t1 e t2 , se tem PL′′ (x) = −0.012x + 0.12 ⇒
PL′′ (−0.801) ≈ 0.13 > 0 e PL′′ (20.8) ≈ −0.13 < 0. Portanto, o pico de máxima demanda será em t = 20.8 dias,
ou seja, no dia 11 de fevereiro, cuja demanda será aproximadamente de 20.04 Mw.
2 Soluções da 2ª lista – interpolação

(d) polinômio interpolador de Newton; para obter o polinômio de Newton é necessário obter a tabela de diferenças
dividas, apresentada na tabela 4 a seguir, obtida a partir da tabela 3.

Tabela 4: Tabela de diferenças divididas obtidas da tabela 3.

t D ∆1 ∆2 ∆3
0 10
0.5
10 15 0
0.5 -0.002
c0
20 20 -0.06
c
-0.7 1
30 13

Agora, para obter o polinômio Newton para a data indicada (12 de fevereiro ⇒ t = 22), como 20 < t < 30 ⇒
PN (t) = 20 − 0.7(t − 20) ⇒ PN (22) = 20 − 0.7(22 − 20) = 18.6 Mw (mesmo valor obtido no item a).

5. Obtenha uma estimativa para x = 0.032 usando o polinômio interpolador de Newton para a tabela 5 e calcule o
limitante para o erro cometido.

Tabela 5: Dados coletados para f (x) desconhecida.

x 0.000 0.100 0.300 0.400
y 1.000 0.761 0.067 -0.367

Com base na tabela 5, se obtém a tabela de diferenças divididas (tabela 6).

Tabela 6: Tabela de diferenças divididas obtidas da tabela 5.

x y ∆1 ∆2 ∆3
0.000 1.000 c0
0.032
-2.390 c1
c2
0.100 0.761 -3.600
-3.470 1.750 c3
0.300 0.067 -2.900
-4.340
0.400 -0.367

Como 0 < 0.032 < 0.1 ⇒ PN (x) = c0 + c1 (x − x0 ) + c2 (x − x0 )(x − x1 ), escolhendo polinômio de grau 2, para que
se possa usar ∆3 na estimativa do limitante do erro. Assim, se tem PN (x) = 1 − 2.39(x − 0) − 3.6(x − 0)(x − 0.1) ⇒
∆3 (x − 0)(x − 0.1)(x − 0.3)
PN (0.032) ≈ 0.931 e o limitante do erro ϵ = ⇒ ϵ ≤ 0.00017
3!
6. Uma hidrelétrica tem capacidade máxima de 60 Mw, que é determinada por 3 geradores de, respectivamente, 30
Mw, 15 Mw e 15 Mw. A demanda de energia da cidade atendida por esta hidrelétrica varia num ciclo de 24 h e é
em função dela que o engenheiro operacional distribui as tarefas dos geradores.
Sabe-se que a demanda mı́nima ocorre entre 1 e 5 h da manhã e a demanda máxima entre 13 e 17 horas, encontre
essas demandas máximas e mı́nimas e em que horas ocorrem, a partir da tabela 7, usando interpolação quadrática
e lagrangeana.

Tabela 7: Demandas de geração de energia elétrica.

hora 2:00 3:00 4:00 5:00 13:00 14:00 15:00 16:00 17:00
Demanda (Mw) 16.4 15.2 14.9 16.0 28.0 36.5 43.0 34.9 31.2
Prof. Geraldo L. Diniz 3

Para simplificar os cálculos, principalmente para se obter o polinômio de Lagrange, ao invés de trabalhar com a
tabela toda, com base nas informações do problema, a tabela 7 pode ser subdividida em duas partes, uma com os
dados da manhã e outra com os dados da tarde. Daı́, considerando os dados da parte da manhã se observa que a
demanda mı́nima deve ocorrer entre 3:00 e 5:00.

– o que fornece os 3 pares de pontos para a interpolação quadrática, dada por:

(h − 4)(h − 5) (h − 3)(h − 5) (h − 3)(h − 4)
P2 (h) = 15.2 + 14.9 + 16.0 ⇒ P2 (h) = 24.5 − 5.2h + 0.7h2 .
(3 − 4)(3 − 5) (4 − 3)(4 − 5) (5 − 3)(5 − 4)
Agora, calculando a derivada primeira de P2 (h) e igualando a zero, para obter o ponto de mı́nimo, se chega a
5, 2
equação P2′ (h) = 1, 4h − 5, 2 ⇒ 1, 4h − 5, 2 = 0 ⇒ h = ≈ 3, 7. Logo, usando polinômio de grau 2, a menor
1, 4
demanda seria por volta de 3.7 horas, ou seja, 3:42. Além disso, a demanda nesta hora será P2 (3.7) ≈ 14.8
Mw.
– considerando todos os pontos do perı́odo da manhã, para a interpolação lagrangeana, se tem
(h − 3)(h − 4)(h − 5) (h − 2)(h − 4)(h − 5) (h − 2)(h − 3)(h − 5) (h − 2)(h − 3)(h − 4)
PL (h) = 16.4 +15.2 +14.9 +16.0 ⇒
(2 − 3)(2 − 4)(2 − 5) (3 − 2)(3 − 4)(3 − 5) (4 − 2)(4 − 3)(4 − 5) (5 − 2)(5 − 3)(5 − 4)
2 3 ′ 2
PL (h) = 19.5−1.283h−0.3h +0.083h , cuja derivada primeira é PL (h) = 0.25h −0.6h−1.283 que tem raı́zes h1 = −1.36
e h2 = 3.8.
É fácil verificar que h1 é ponto de máximo enquanto h2 é ponto de mı́nimo (verifique!), através da derivada segunda
calculada nestes dois pontos. Logo, pela interpolação de Lagrange a demanda mı́nima será em h = 3.8 horas, ou seja,
3:48. Nesta hora a demanda, por Lagrange, será PL (3.8) ≈ 14.9 Mw. A figura 1 ilustra as interpolações obtidas com os
dados da tabela 7 para o perı́odo da manhã.

(a) Interpolação quadrática. (b) Interpolação por Lagrange.

Figura 1: Gráficos das interpolações com os dados da tabela 7 perı́odo da manhã.

Para encontrar a hora e pico de demanda máxima, usando os dados do perı́odo da tarde, se tem

– para interpolação quadrática, verificando os dados do perı́odo da tarde, a demanda será máxima entre 14:00 e
16:00, daı́ o polinômio quadrático é dado por
(h − 15)(h − 16) (h − 14)(h − 16) (h − 14)(h − 15)
P2 (h) = 36.5 + 43.0 + 34.9 ⇒
(14 − 15)(14 − 16) (15 − 14)(15 − 16) (16 − 14)(16 − 15)
2
P2 (h) = −1.587.5 + 218.2h − 7.3h , cuja derivada primeira é P2′ (h) = −14.6h + 218.2, então P2′ (h) = 0 ⇒
− 14.6h + 218.2 = 0 ⇒ h = 14.9. Logo, a demanda máxima será as 14.9 horas, ou seja, 14:54, cujo pico neste horário,
pela interpolação quadrática, é dado por P2 (14.9) ≈ 43.02 Mw.
– para a intrepolação lagrangeana, usando todos os valores do perı́odo da tarde, se tem
(h − 14)(h − 15)(h − 16)(h − 17) (h − 13)(h − 15)(h − 16)(h − 17)
PL (h) = 28.0 + 36.5 +
(13 − 14)(13 − 15)(13 − 16)(13 − 17) (14 − 13)(14 − 15)(14 − 16)(14 − 17)
(h − 13)(h − 14)(h − 16)(h − 17) (h − 13)(h − 14)(h − 15)(h − 17)
+ 43.0 + 34.9 +
(15 − 13)(15 − 14)(15 − 16)(15 − 17) (16 − 13)(16 − 14)(16 − 15)(16 − 17)
(h − 13)(h − 14)(h − 15)(h − 16)
+ 31.2 ⇒ PL (h) = 1.32(291.011 − 34.045h + h2 )(164.369 − 25.550h + h2 ),
(17 − 13)(17 − 14)(17 − 15)(17 − 16)
4 Soluções da 2ª lista – interpolação

cuja derivada primeira é PL′ (h) = 5.27(h−16.7)(h−14.9)(h−13.1). Logo, suas raı́zes são h1 = 13.1, h2 = 14.9 e h3 = 16.7.
Para verificar qual será o ponto de máximo, usando a derivada segunda, dada por PL′′ (h) = 15.8(h − 15.9)(h − 13.8),
então PL′′ (13.1) ≈ 35.5 > 0, PL′′ (14.9) ≈ −17.4 < 0 e PL′′ (16.7) ≈ 34.1 > 0.
Logo, h1 e h3 são pontos de mı́nimo, enquanto h2 é ponto de máximo. Portanto, o pico de demanda será as 14.9 horas,
ou seja, as 14:54. Neste momento a demanda máxima será PL (14.9) ≈ 43.1 Mw. A figura 2 ilustra as interpolações
obtidas com os dados da tabela 7 para o perı́odo da tarde.

(a) Interpolação quadrática. (b) Interpolação por Lagrange.

Figura 2: Gráficos das interpolações com os dados da tabela 7 perı́odo da tarde.

(*) A critério de comparação, usando o método de Newton para determinar o horário e o pico de demanda máxima
no periodo da tarde, a resolução é a que se segue.
Com base na tabela 8 a demanda máxima deve ocorrer entre
14 e 17 horas. Assim, o polinômio de Newton para esta
interpolação é dado por:
Tabela 8: Diferenças divididas PN (t) = c0 +c1 (t−14)+c2 (t−14)(t−15)+c3 (t−14)(t−15)(t−16),
t D ∆1 ∆2 ∆3 ∆4 em que c0 = 36.5, c1 = 6.5, c2 = −7.3 e c3 ≈ 3.17. Daı́,
13 28 PN (t) = 3.17t3 − 149.95t2 + 2354.78t − 12238.6
8.5 Agora, para obter os pontos crı́ticos de PN (t), vem
14 36.5 0
c -1 PN′
(t) = 6.34t2 − 299.9t + 2354.78 ⇒
6.5
c 1
-2.1
′
PN (t) = 9.51(t − 14.76)(t − 16.77), logo PN ′
(t) = 0 ⇔ t = 14.76
15 43 -7.3
c 2
1.31667 ou t = 16.77.
c3 Para determinar se é ponto de máximo ou mı́nimo, usando a de-
-8,1 3.1667
rivada segunda vem
16 34.9 2.2 ′′ ′′ ′′
PN (t) = 19.02t−299.9 ⇒ PN (14.76) = −19.16 < 0 e PN (16.77) =
-3.7
19.07 > 0, ou seja, a demanda máxima ocorre para t = 14.76.
17 31.2
Portanto, a demanda máxima será PN (14.76) = 43.49 e ocorrerá
as 14:45.
Além disso, é possı́vel determinar o erro na estimativa do valor desta demanda máxima, que é dado por

M4
ξ≤ (t − 14)(t − 15)(t − 16)(t − 17)
4!
em que M4 é o maior valor das diferenças de ordem 4 (em valor absoluto). Assim,
1.31667
ξ≤ (14.76 − 14)(14.76 − 15)(14.76 − 16)(14.76 − 17) ≈ 0.028
24

Cuiabá, 20 de março de 2023.