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
-4
51.9
-2
4
-5.2
-1
-1
13.4
0
2
-7.5
1
2
-2.0
2
-2
25.0
3
4
12.0
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,8). 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) 156.9
(B) 157.9
(C) 158.9
(D) 159.9
(E) 160.9
(F) 161.9
(G) 162.9
(H) 163.9
Passo 1
O enunciado dá que a grandeza \,Z\, segue o modelo \,Z=a x^{2}+b y+c\, e pede que ajustemos os parâmetros \,a,b,c\, por mínimos quadrados usando as 7 medições da tabela e, em seguida, que prevejamos o valor de \,Z\, no ponto (x,y)=(8,8). Peguei os valores diretamente da tabela fornecida; os pontos medidos são:
(x,y,Z)=(-3,-4,51{,}9)
(x,y,Z)=(-2,4,-5{,}2)
(x,y,Z)=(-1,-1,13{,}4)
(x,y,Z)=(0,2,-7{,}5)
(x,y,Z)=(1,2,-2{,}0)
(x,y,Z)=(2,-2,25{,}0)
(x,y,Z)=(3,4,12{,}0)
O que se pede:
Achar \,a,b,c\, que minimizam a soma dos quadrados das diferenças entre as medições e o modelo.
Calcular \,Z(8,8)=a\cdot 8^{2}+b\cdot 8 + c\, e escolher a melhor aproximação entre as opções dadas.
Passo 2
Ideia do método (rápida e intuitiva). Para um modelo linear nos parâmetros (aqui linear porque a,b,c aparecem linearmente se tomarmos x^{2} como variável), a forma padrão é usar o método dos mínimos quadrados lineares. Monta-se a matriz de projeto X com colunas correspondentes às funções multiplicadoras dos parâmetros, aqui \,x^{2},\,y,\,1\,. Escrevendo o vetor de parâmetros \theta=\begin{pmatrix}a\\[4pt] b\\[4pt] c\end{pmatrix} e o vetor de observações z, os mínimos quadrados são obtidos resolvendo as equações normais
X^{T}X\,\theta = X^{T} z.
Vou explicar cada termo e calcular tudo passo a passo para chegar nos valores de \,a,b,c\,.]
Passo 3
Construção das colunas de X e cálculo das somas necessárias (peguei cada x,y,Z da tabela).
Para cada medição calculei \,x^{2},\,x^{4},\,x^{2}y,\,x^{2}Z,\,y^{2},\,yZ\, porque aparecem nas somas de X^{T}X e X^{T}z. A tabela auxiliar dos cálculos (cada número está entre parênteses na unidade usada na tabela):
Para (x,y,Z)=(-3,-4,51{,}9): x^{2}=9,\;x^{4}=81,\;x^{2}y=9\cdot(-4)=-36,\;x^{2}Z=9\cdot51{,}9=467{,}1,\;y^{2}=16,\;yZ=-4\cdot51{,}9=-207{,}6.
Para (x,y,Z)=(-2,4,-5{,}2): x^{2}=4,\;x^{4}=16,\;x^{2}y=16,\;x^{2}Z=4\cdot(-5{,}2)=-20{,}8,\;y^{2}=16,\;yZ=4\cdot(-5{,}2)=-20{,}8.
Para (x,y,Z)=(-1,-1,13{,}4): x^{2}=1,\;x^{4}=1,\;x^{2}y=-1,\;x^{2}Z=13{,}4,\;y^{2}=1,\;yZ=-13{,}4.
Para (x,y,Z)=(0,2,-7{,}5): x^{2}=0,\;x^{4}=0,\;x^{2}y=0,\;x^{2}Z=0,\;y^{2}=4,\;yZ=2\cdot(-7{,}5)=-15{,}0.
Para (x,y,Z)=(1,2,-2{,}0): x^{2}=1,\;x^{4}=1,\;x^{2}y=2,\;x^{2}Z=-2{,}0,\;y^{2}=4,\;yZ=-4{,}0.
Para (x,y,Z)=(2,-2,25{,}0): x^{2}=4,\;x^{4}=16,\;x^{2}y=-8,\;x^{2}Z=4\cdot25{,}0=100{,}0,\;y^{2}=4,\;yZ=-50{,}0.
Para (x,y,Z)=(3,4,12{,}0): x^{2}=9,\;x^{4}=81,\;x^{2}y=36,\;x^{2}Z=108{,}0,\;y^{2}=16,\;yZ=48{,}0.
Agora somo coluna a coluna (escrevo os valores como números):