ExRes - Interp e MMQ
ExRes - Interp e MMQ
ExRes - Interp e MMQ
2. Use o Teorema do Erro, e determine uma cota superior do erro, para as aprox-
imações calculadas no exercı́cio 1
3. Sabendo-se que f (0, 81) = 16, 94410 f (0, 83) = 17, 56492 f (0, 86) = 18, 50515
e f (0, 87) = 18, 82091, calcule um valor aproximado de f (0, 84), usando:
Pergunta-se:
11. O número de bactérias, por unidade de volume, existente em uma cultura após
x horas é dado na tabela abaixo:
número de horas 0 1 2 3 4 5 6
número de bactérias 32 47 65 92 132 190 275
(a) Ajuste os dados acima a curva y = aebx pelo método dos mı́nimos quadra-
dos.
(b) Quantas horas seriam necessárias para que o número de bactérias por
unidade de volume ultrapasse 2000?
12. Dada a tabela abaixo, faça o gráfico de dispersão dos dados e ajuste uma curva
da melhor maneira possı́vel.
x 0, 5 0, 75 1 1, 5 2, 0 2, 5 3, 0
y −2, 8 −0, 6 1 3, 2 4, 8 6, 0 7, 0
Exer cício 1:
(a) f(x) = cos(x)
P2 ( x ) = L 0 ( x ) f ( x 0 ) + L1 ( x ) f ( x 1 ) + L 2 ( x ) f ( x 2 )
onde :
Portanto,
Erro :
f(x) - P2 (x) = cos(0,45) − P2 (0,45) ≈ 2,822 ⋅ 10 −3
P2 ( x ) = d 0 + d1 ( x − x 0 ) + d 2 ( x −x 0 )( x − x 1 )
onde :
d0 = f [ x 0 ]
f [ x1 ] − f [ x 0 ]
d1 = f [ x 0 , x 1 ] =
x1 − x 0
f [ x1 , x 2 ] − f [ x 0 , x1 ]
d 2 = f [ x 0 , x1 , x 2 ] =
x2 − x0
Vamos montar a seguinte tabela:
Exer cício 2:
M n +1
E n (x) = f(x) − Pn (x) ≤ (x − x 0 )(x − x 1 ) ⋅ ⋅ ⋅ (x − x n )
(n + 1) !
onde :
(n +1)
M n +1 = máx f (x) para x ∈ [x 0 , x n ]
Então ,
f ′′′( x ) máx
E 2 ( x ) = f ( x ) − P2 ( x ) ≤ ( x − x 0 )( x − x 1 )( x − x 2 )
3!
f ′′(0)
′ = 0 ; f ′′(0,6)
′
≈ 0,565 ; f ′′(0,9)
′ ≈ 0,7833
Máximo
0,7833
E 2 (0,45) ≤ (0,45 − 0)(0,45 − 0,6)(0,45 − 0,9) ≈ 3,965 ⋅ 10 −3
3!
Exer cício 3:
(a) Devemos neste item construir por Lagrange P1(x), P2(x), P3(x) tais que:
com x 0 = 0,83 , x 1= 0,86 e x 2 = 0,87 (lembrando que escolhemos para x 0 o valor mais próximo de x)
P1 ( x ) = d 0 + d1 ( x − x 0 )
P2 ( x ) = d 0 + d1 ( x − x 0 ) + d 2 ( x − x 0 )( x − x 1 )
P3 ( x ) = d 0 + d1 ( x − x 0 ) + d 2 ( x − x 0 )( x − x 1 ) + d 3 ( x − x 0 )( x − x 1 )( x − x 2 )
(c) Se a função f(x) é dada na forma de tabela, o valor absoluto do erro |En(x)| só pode ser estimado.
Isto porque, neste caso, não é possível calcular Mn+1; mas, se construirmos a tabela de diferenças
divididas até ordem n+1, podemos usar o maior valor (em módulo) destas diferenças como uma
M n +1
aproximação para no intervalo [x0 , xn].
( n + 1)!
Assim,
Não é possível determinar |E3(x)| porque não temos as diferenças divididas de ordem 4.
Exer cício 4:
Neste exercício, temos pontos xi igualmente espaçados. Sendo h o passo, temos:
x 1 − x 0 = x 2 − x 1 = ⋅ ⋅ ⋅ = x n − x n −1 = h
f ′′( x ) máx
E 1 ( x ) = f ( x ) − P1 ( x ) ≤ ( x − x 0 )( x − x 1 )
2!
Para achar ( x − x i )( x − x i +1 ) máx , basta verificarmos que como se trata de uma parábola, a coordenada
x + x i +1
que contém o valor máximo para w(x) = ( x − x i )( x − x i +1 ) é (x vértice , y vértice ) = i , y vértice
2
Então ,
h h
M x i + x i +1 x i + x i +1 M x i +1 − x i x i − x i +1 M x i +1 − x i x i − x i +1 Mh 2
E1 ( x ) ≤ − xi − x i +1 = = =
2! 2 2 2 2 2 2 2 2 8
Exer cício 5:
h n +1 M n +1
E n (x) = f(x) − Pn (x) <
4(n + 1)
Exer cício 6:
(a) Vamos ordenar a tabela por peso:
Usando P2 ( x ) = d 0 + d1 ( x −x 0 ) + d 2 ( x −x 0 )( x − x1 ) , temos:
5 1
P2 (70) = 173 + (70 − 69) − (70 − 69)(70 − 73) = 174,375 cm
4 24
Ù [ = Â sen(x) + Ã cos(x)
Vamos ajustá-la aos dados da tabela através do Método dos Mínimos Quadrados, fazendo:
4
S(Â , Ã = ∑ [f(x i ) − Â sen(x i ) − Ã cos(x i ) ] 2
i =0
onde :
4
∂S
∂Â
=0 ⇒ 2
i =0
∑
[f(x i ) −Â sen(x i ) − Ã cos(x i ) ] ⋅ [-sen(x i ) ] = 0 (1)
4
∂S
∂Ã
=0 ⇒ 2
i =0
∑
[f(x i ) −Â sen(x i ) − Ã cos(x i ) ] ⋅ [-cos(x i ) ] = 0 (2)
4 4 4
2α ∑ sen
i =0
2
( xi ) + β ∑ sen ( 2x ) = 2∑ f ( x )sen ( x )
i =0
i
i =0
i i (3)
4 4 4
α ∑
i =0
sen ( 2x i ) + 2β ∑
i =0
cos 2 ( x i ) = 2∑ f (x
i =0
i ) cos( x i ) (4)
α ≈ 190,717
β ≈ -15,187
Agora, vamos usar essa equação para achar a altura aproximada de uma pessoa de 70 Kg :
Dados :
( x i , f ( x i )) , i = 0,1,..., m (Tabela de f)
ϕ 0 ( x ), ϕ1 ( x ), ..., ϕ n ( x ) (Funções quaisquer contínuas)
onde c i ∈ R , i = 0,1,..., n
Teríamos então:
c 0 ϕ 0 ( x 0 ) + c 1ϕ 1 ( x 0 ) + ... + c n ϕ n ( x 0 ) = f ( x0 )
c 0 ϕ 0 ( x 1 ) + c 1ϕ 1 ( x 1 ) + ... + c n ϕ n ( x 1 ) = f ( x1 )
. . . .
. . . .
. . . .
c 0 ϕ 0 ( x m ) + c 1ϕ1 ( x m ) + ... + c n ϕ n ( x m ) = f ( x m )
É possível obter um mesmo polinômio que interpola e faz o ajuste de curvas pelo MMQ se o modelo ajustar
m
exatamente os dados. Dessa forma, o mínimo de S(c 0 , c 1 ,..., c n ) = ∑[ f (x
k =1
k ) − g ( x k ) ] 2 será zero e, portanto,
Peso(Kg) 63 69 73 79 82
Altura(cm) 163 173 178 183 188
Velocidade(km/h) 14 16 15 15 14
P2 (x, y) = f(x 0 , y 0 )L 0 (x)L 0 (y) + f(x 0 , y1 )L 0 (x)L1 (y) + f(x 0 , y 2 )L 0 (x)L 2 (y) +
f(x 1 , y 0 )L1 (x)L 0 (y) + f(x 1 , y1 )L1 (x)L1 (y) + f(x 1 , y 2 )L1 (x)L 2 (y) +
f(x 2 , y0 )L 2 (x)L 0 (y) + f(x 2 , y1 )L 2 (x)L1 (y) + f(x 2 , y 2 )L 2 (x)L 2 (y)
Observação: na hora de calcular f(xi , yj), colocamos xi como ponto fixo (que não varia).
Depois, verificamos o valor de f(xi , yj) , a velocidade representada neste exercício, no ponto
yj .
Portanto,