Baixe no formato PDF, TXT ou leia online no Scribd
Fazer download em pdf ou txt
Você está na página 1/ 9
Cálculo Numérico IV
Lista 5 — 2020.2 Ajuste de Curvas
Prof. Americo Cunha
Profa. Nancy Baygorrea
Use o GNU Octave para lhe auxiliar com as contas complicadas. Algumas questões necessitam de justificativa, i.e., não basta fazer contas, o raciocı́nio deve ser explicado através de um texto coerente.
Exercı́cio 1 Considere os pares ordenados apresentados na tabela a seguir.
z 0 1 2 3 4 p(z) 0 1 1 4 4
Encontre a reta que melhor de ajusta a esses pontos no sentido dos mı́nimos quadrados.
Exercı́cio 2
A tabela abaixo apresenta o resultado de um teste realizado com 6 digitadores de idades diferentes. Nesse teste avaliou-se quantas palavras por minuto cada um é capaz de digitar.
idade do digitador palavras por minuto
25 78 30 70 35 65 40 58 45 48 50 42
1. Utilizando o método dos mı́nimos quadrados, encontre a reta y = c0 + c1 x que melhor ajusta os pontos da tabela, sendo x a idade do digitador e y o número de palavras por minuto que ele é capaz de digitar. 2. Use essa reta para estimar quantas palavras por minuto um digitador de 55 anos é capaz de digitar. Sua estimativa parece raozável?
1 Exercı́cio 3 É um fato conhecido de engenharia que a velocidade v de um carro de corrida depende da pressão p no aerofólio dianteiro. Durante uma bateria de testes com um certo carro, engenheiros fizeram medições de como v se relaciona com p. Essas medidas podem ser vistas na tabela a seguir.
v (m/s) 7, 87 11, 50 16, 40 22, 60 32, 80
p (kP a) 14 28 56 112 225
Esses engenheiros desejam encontrar uma função que lhes permita, mesmo sem conhecer por com- pleto a fı́sica do fenômenos observado, calcular a pressão no aerofólio correspondente a uma dada velocidade do carro. Obtenha a função desejada pelos engenheiros através do método dos mı́nimos quadrados, ajustando os pontos da tabela uma curva da forma p(v) = c0 + c1 v 2 .
Exercı́cio 4 A tabela abaixo apresenta o resultado de um experimento feito em laboratório, onde se mediu a posição de uma partı́cula em 6 instantes diferentes.
tempo (s) posição (m)
1,2 1,1 2,3 2,1 3,0 3,1 3,8 4,0 4,7 4,9 5,9 5,9
1. Utilizando o método dos mı́nimos quadrados, encontre a reta r(t) = c0 + c1 t que melhor ajusta os pontos da tabela no sentido dos mı́nimos quadrados, sendo r(t) a posição da partı́cula e t o tempo. 2. Use essa reta para estimar a posição da partı́cula no instante t = 7 s.
Exercı́cio 5 A tabela que se segue apresenta os valores experimentais da pressão p de uma dada massa de gás, que ocupa um certo volume v.
v (litro) 54, 3 64, 8 72, 4 88, 7 118, 6 194, 0
p (atm) 61, 2 49, 5 37, 6 28, 4 19, 2 10, 1
Um gás ideal, que é submetido a um processo adiabátivo reversı́vel, tem essas duas grandezas relacionadas por uma equação da forma p v γ = C, onde γ é o chamado ı́ndice adiabático e C é uma constante do processo. Determine parâmetros C e γ de forma a ajustar a curva p v γ = C aos pontos da tabela. Utilize a curva obtida para estimar a pressão correspondente a v = 100 litros. (Dica: ln p v γ = ln C ⇐⇒ ln p = −γ ln v + ln C)
2 Exercı́cio 6 A tabela abaixo apresenta o resultado de um experimento feito em laboratório, onde se mediu a posição de uma massa, presa a uma mola, em 6 instantes diferentes.
1. Utilizando o método dos mı́nimos quadrados, encontre a curva x(t) = c0 cos (t) + c1 sin (t) que melhor ajusta os pontos da tabela, sendo x a posição da massa e t o tempo. 2. Use a curva ajustada para estimar a posição da partı́cula no instante t = 20 s.
Exercı́cio 7 A tabela abaixo mostra a evolução de uma certa população, ao longo de um perı́odo de 20 anos, que apresenta crescimento exponencial. Utilize o método dos mı́nimos quadrados para ajustar uma curva da forma N = N0 eαt aos pontos dessa tabela.
Considere os pares ordenados apresentados na tabela a seguir.
x 0 1 2 y 1.1 3 8.12
Encontre, utilizando o método dos mı́nimos quadrados, a curva da forma y = a ex + b x que melhor ajusta esses pontos.
3 Questões de múltipla escolha
Exercı́cio 9 Dada a tabela a seguir, obtenha a melhor aproximação para x∗ tal que f (x∗ ) = 0, 6018, utilizando um polinômio do primeiro grau.
x f(x) 30 0,500 35 0,5736 45 0,7071
a) 37, 00 b) 37, 11 c) 37, 98 d) 39, 82 e) Não sei
Exercı́cio 10 A tabela abaixo apresenta o resultado de um teste realizado com 3 digitadores de idades diferentes. Nesse teste avaliou-se quantas palavras por minuto cada um é capaz de digitar.
Idade do Digitador Palavras por Minuto
30 70 35 65 40 58
Utilizando o método dos mı́nimos quadrados, encontre a reta y = c0 + c1 x que melhor ajusta os pontos da tabela, sendo x a idade do digitador e y o número de palavras por minuto que ele é capaz de digitar.
a) y = 100 − x b) y = 11 − 1, 4x c) y = 106, 67 − 2, 97x d) y = 106, 33 − 1, 2x e) Não sei
Exercı́cio 11 Utilizando o método dos mı́nimos quadrados, ajuste uma função do tipo g(x) = a + bx2 aos pontos da tabela que segue e determine g(1, 5).
x f(x) 0,0 1,1 1,0 0,1 2,0 -3,1
a) −0, 1470
4 b) −2.4825 c) −1.2425 d) −0, 1500 e) Não sei
Exercı́cio 12 11 Se a reta y = 10 x é a que melhor ajusta os pontos (1, 1), (2, 2), (3, a) e (4, 4) no sentido dos mı́nimos quadrados, então
a) a = 4.
b) a = 3. c) a = 5. d) a = 0. e) Não sei.
Exercı́cio 13 Qual é a curva da forma y = c0 sin x + c1 cos x que melhor ajusta, no sentido dos mı́nimos quadrados, os pontos da tabela a seguir?
x 0 π/6 √π/4 √π/3 π/2
y 0 1/2 2/2 3/2 1
a) y = cos x b) y = sin x + cos x c) y = sin x d) y = sin x − cos x
e) Não sei.
5 Gabarito da Lista 5
As repostas e soluções apresentadas a seguir são para auxiliar na resolução desta lista, mas não estão isentas de possı́veis erros de digitação ou mesmo de desenvolvimento. Use o gabarito com cautela, exercitando sempre o seu senso crı́tico. Se encontrar algum erro, por favor, reporte ao professor.
Resposta do Exercı́cio 1 A reta que melhor se ajusta ao pontos da tabela é 11 1 p(z) = z− . 10 5
Resposta do Exercı́cio 2 1. Para encontrar os parâmetros de ajuste (c0 e c1 ) é necessário resolver a equação normal associada ao problema de mı́nimos quadrados AT Ac = AT b, onde 1 25 78 1 30 70 1 35 , c = c0 , b = 65 A= . 1 40 c1 58 1 45 48 1 50 42 Desenvolvendo os produtos matriciais temos que
6 225 c0 361 = , 225 8875 c1 12905 donde
c0 1 8875 −225 361 2402/21 = = . c1 2625 −225 6 12905 −253/175 253 Portanto, reta que ajusta os pontos da tabela é y = − 175 x + 2402 21 . 2. Substituindo x = 55 nessa reta obtemos y = − 175 × 55 + 2402 253 21 ≈ 34, 87. Logo, essa reta estima que um digitador de 55 anos é capaz de digitar aproximadamente 35 palavras por minuto.
Resposta do Exercı́cio 3 A função que ajusta os pontos da tabela é dada por
p(v) = 1, 203 + 0, 209 v 2 ,
onde p é medido em kP a e v em m/s.
6 Resposta do Exercı́cio 4 1. Para encontrar os parâmetros de ajuste (c0 e c1 ) é necessário resolver a equação normal associada ao problema de mı́nimos quadrados
c0 1 87, 07 −20, 9 21, 1 −0, 1432542927 = = . c1 85, 61 −20, 9 6 88, 49 1, 050695012 Logo, temos c0 ≈ −0, 1 e c1 ≈ 1, 1. Portanto, reta que ajusta os pontos da tabela é r(t) = 1, 1 t − 0, 1. 2. Substituindo t = 7 s nessa reta obtemos r = 1, 1 × 7 − 0, 1 = 7, 6 m. Logo, essa reta estima que em t = 7 s a partı́cula estará na posição r = 7, 6 m
Resposta do Exercı́cio 5 Os parâmetros de ajuste são C = 16019, 17 atm × litros1,40487 e γ = 1, 40487, de forma que a curva de ajuste é p v 1,40487 = 16019, 17. Essa equação pode ser reescrita como p(v) = 16019, 17 v −1,40487 e pode ser utilizada para estimar a pressão correspondente a v = 100 litros. Assim sendo, temos p(100 litros) ≈ 24, 8 atm.
Resposta do Exercı́cio 6 1. Para encontrar os parâmetros de ajuste (c0 e c1 ) é necessário resolver a equação normal associada ao problema de mı́nimos quadrados AT Ac = AT b, onde cos (3, 8) sin (3, 8) −1, 3 cos (7, 2) sin (7, 2) 1, 2 cos (9, 4) sin (9, 4) c 0 −1, 1 A= , c = , b= . cos (11, 9) sin (11, 9) c 1 0, 5 cos (14, 8) sin (14, 8) 0, 3 cos (18, 5) sin (18, 5) 0, 5 Desenvolvendo os produtos matriciais e arredondando com duas casas decimais temos que
c0 1 2, 13 0, 35 −0, 12 −0, 05 = = . c1 8, 11 0, 35 3, 87 −0, 37 −0, 18 Logo, temos c0 ≈ −0, 05 e c1 ≈ −0, 18. Portanto, a curva que ajusta os pontos da tabela é x = −0, 05 cos (t) − 0, 18 sin (t). 2. Substituindo t = 20 s na equação x = −0, 05 cos (t) − 0, 18 sin (t) obtemos x = −0, 05 cos (20) − 0, 18 sin (20) ≈ −0, 1192015684 m. Logo, essa curva estima que em t = 20 s a partı́cula estará na posição x = −0, 12 m.
Resposta do Exercı́cio 7 Aplicando o logaritmo natural aos dois lado dessa equação obtemos ln N = ln N0 eαt ⇐⇒ ln N = ln N0 + α t.
Para encontrar os parâmetro α e ln N0 é necessário resolver a equação normal associada ao problema
de mı́nimos quadrados AT Ax = AT b, onde 1 0 ln 100 1 5 ln 165 , x = ln N0 , b = A= 1 10 ln 272 . α 1 15 ln 448 1 20 ln 738 Desenvolvendo os produtos matriciais, e arredondando, temos que
5 50 ln N0 28, 03 = , 50 750 α 305, 24 donde
ln N0 4, 61 = . α 0, 10 Assim, N0 ≈ 100, 48. Portanto, a curva de ajuste obtida é N = 100, 48 e0,10 t .
Resposta do Exercı́cio 8 Para encontrar os parâmetros de ajuste (a e b) é necessário resolver a equação normal associada ao problema de mı́nimos quadrados AT Ac = AT b, onde
e0 0 1.1 a A = e1 1 , c = , b= 3 . b e2 2 8.12 Desenvolvendo os produtos matriciais (com duas casas decimais) temos que
62.99 17.50 a 69.25 = , 17.50 5 b 19.24 donde
a 1 62.99 17.50 69.25 1.09 = = . b 8.81 17.50 5 19.24 0.02 Portanto, reta que ajusta os pontos da tabela é y = 1.09 ex + 0.02 x.
