Q13: Em um determinado fenômeno físico, uma grandeza Z depende de duas variáveis de controle x e y. O modelo matemático para esta dependência é Z=a x^{2}+b y+c, onde os parâmetros a, b e c dependem do arranjo experimental. Para se determinar estes parâmetros, foram realizadas 7 medições de Z para diferentes configurações de x e y e os resultados obtidos estão na tabela abaixo:
x
y
Z
-3
2
47.6
-2
2
23.3
-1
5
14.5
0
-5
-18.5
1
1
3.5
2
4
29.2
3
4
54.4
Com base nos dados experimentais acima, ajuste os parâmetros a, b e c de forma a minimizar (no sentido dos mínimos quadrados) o desvio entre o modelo e os dados. Feito isto, calcule o valor de Z previsto pelo modelo ajustado para (x, y)=(8,6). Assinale a melhor aproximação para esse valor, dentre as opções abaixo.
Sugestão: Use algum software para resolver este problema. Pode ser uma planilha ou alguma linguagem de programação (como Python ou Octave, por exemplo). Vários destes ambientes estão disponíveis online.
(A) 342.9
(B) 343.9
(C) 344.9
(D) 345.9
(E) 346.9
(F) 347.9
(G) 348.9
(H) 349.9
Passo 1
O enunciado pede que se ajuste por mínimos quadrados os parâmetros a,b,c do modelo
Z = a x^{2} + b y + c
usando as 7 medições fornecidas na tabela, e depois use o modelo ajustado para prever Z em (x,y)=(8,6).
Peguei os dados da tabela (cada linha é um ponto (x,y,Z)):
Determinar a,b,c que minimizam a soma dos quadrados dos resíduos entre os Z medidos e os Z preditos pelo modelo.
Calcular o valor previsto de Z para x=8 e y=6.
Comparar com as alternativas e escolher a melhor aproximação.
Passo 2
Breve explicação do método (muito prática): para um modelo linear nos parâmetros (aqui linear em a,b,c), montamos a matriz de desenho X cujas colunas correspondem às funções multiplicadas por cada parâmetro. Aqui a coluna 1 será x^{2}, a coluna 2 será y e a coluna 3 será 1 (para c). Escrevendo em forma matricial:
A solução de mínimos quadrados (normal equations) é obtida por
X^{T}X\,\boldsymbol{\beta} = X^{T}\mathbf{Z}.
Aqui \mathbf{Z} é o vetor coluna com as 7 medidas de Z. Vou aplicar isso diretamente com os dados: é um procedimento padrão — soma das colunas, produtos cruzados, resolver o sistema linear 3×3 para \boldsymbol{\beta}.]
Passo 3
Calcular as quantidades necessárias (peguei os x,y,Z da tabela e calculei x^{2} para cada linha).
Os valores de x^{2} por linha são:
9,\,4,\,1,\,0,\,1,\,4,\,9.
Agora as somas e produtos que entram em X^{T}X e em X^{T}\mathbf{Z}: