Erros e incertezas experimentais

Não consegue ver as equações correctamente? Mude https para http no endereço desta página e recarregue.

Grandezas físicas

Um dos principais objectivos da Física Experimental consiste na medição quantitativa de grandezas físicas. A palavra medição designa o acto de medir, do qual resulta uma medida (ou mais do que uma, no caso de se repetir o processo). É fundamental classificar os principais tipos de grandezas encontradas:

Tipos de grandezas físicas

Directas	São aquelas cujo valor se obtém com uma medição, não sendo necessário envolver os valores de outras grandezas físicas. Exemplos: comprimento $L$ tempo $t$ temperatura $T$ massa $m$
Indirectas	São aquelas que envolvem a medição de duas ou mais grandezas, que por sua vez podem ser directas (D) ou indirectas (I). Exemplos: velocidade (escalar) $v$ - envolve comprimento (D) e tempo (D) área $A$ - envolve comprimento (D) e largura (D) densidade $\rho$ - envolve massa (D) e volume (I)
Escalares	São caracterizadas por um número e pela unidade de medida que a define. Exemplos: temperatura $T$ (K) tempo $t$ (s) massa $m$ (kg)
Vectoriais	Além de um número e uma unidade de medida, é necessário também saber a direção e o sentido destas grandezas. Exemplos: posição espacial $\vec{r}$ (m) velocidade $\vec{v}$ (m/s) força $\vec{F}$ (N)

Definições fundamentais

Incerteza

Em física experimental, um dos conceitos mais importantes é a incerteza de uma medição. A incerteza é uma expressão quantitativa da dúvida que existe na medição, reflectindo aspectos como os limites dos instrumentos e dos métodos experimentais. Em qualquer resultado experimental é indispensável indicar a incerteza associada, uma vez que esta mede a fiabilidade dos resultados e permite comparar medições. Assim, contabilizar correctamente as incertezas é essencial para garantir a validade das conclusões tiradas a partir dos dados experimentais.

Precisão e exactidão

Na linguagem coloquial os termos precisão e exactidão usam-se como sinónimos, mas no método científico experimental traduzem conceitos muito diferentes. Pode existir uma medida exacta e não precisa, ou outra precisa mas não exacta (ver ilustração). O grande mérito de um experimentalista será obter simultaneamente a melhor precisão e a melhor exactidão possíveis.

A precisão de uma medição é o grau da concordância entre determinações repetidas: está relacionado com a variação das sucessivas medidas A exactidão é tanto maior quanto menor for a distância entre a medida (ou a média de determinações repetidas) e um valor “verdadeiro”, “nominal”, “tomado como referência” ou “aceite”. A esta distância damos o nome de erro experimental.

Numa actividade experimental, em regra geral, o valor verdadeiro das grandezas físicas não é conhecido a priori, pelo que naturalmente também não é possível calcular o valor do erro experimental. Nas actividades laboratoriais de LIFE existem algumas excepções em que este valor “verdadeiro/referência” é conhecido com grande precisão/exactidão (e.g Exp. Thomson, Exp. Millikan, Velocidade da Luz, etc). Outras há em que não se conhece o valor verdadeiro (e.g. carga de uma gota de óleo electrizada, temperatura da sala, índice de refração de um material transparente, etc).

Erros sistemáticos e aleatórios

As fontes para a incerteza experimental podem ser muito variadas, mas podem ser classificadas dois tipos principais: os erros de natureza sistemática e os de natureza aleatória. Note-se que em física experimental um erro não significa um engano ou uma falsidade, mas sim a diferença entre um valor medido para uma grandeza e o seu valor "verdadeiro".

Erros sistemáticos

Erros aleatórios

Conduzem em geral a valores sistematicamente desviados (para valores superiores ou inferiores) do valor da grandeza a medir, contribuindo para uma menor exactidão. Podem ser originados por más condições de calibração dos instrumentos de medida o uso destes instrumentos em condições diferentes das que são recomendadas leituras sistematicamente incorrectas do observador (e.g. paralaxe) utilização de um método físico que não é adequado à descrição da experiência. Devem ser corrigidos e minimizados sempre que possível. Só a comparação dos resultados obtidos com outros instrumentos de referência (calibração) pode elucidar se esses erros foram suficientemente reduzidos.

Resultam das flutuações aleatórias que se observam nos resultados obtidos para diferentes leituras, contribuindo para uma menor precisão. Podem ser originados por falta de sensibilidade dos instrumentos e do observador leituras incorrectas (mas não sistemáticas) ruído (vibrações mecânicas ou eléctricas) processos estatísticos intrínsecos ao fenómeno observado (por exemplo declínio radioactivo). A análise estatística das flutuações está fora do âmbito da LIFE, mas importa referir que, habitualmente, considera-se o valor médio dos erros aleatórios como zero. Isto é importante pois, ao repetirem-se as medições e fazendo a média aos $N$ resultados, os erros aleatórios compensam-se, reduzindo-se assim a contribuição aleatória. Os erros aleatórios podem e devem ser sempre caracterizados, mas dado o seu carácter estocástico, não podem ser eliminados totalmente, pelo que qualquer medição tem associada uma incerteza experimental.

Em conclusão, podemos resumir todos estes conceitos nestes pontos:

1. Toda a medição experimental é sujeita a um erro experimental. Só por grande coincidência o valor numérico obtido pela medição é igual ao valor verdadeiro da grandeza.
2. Antes da experiência, devemos identificar e corrigir os erros sistemáticos de todas as grandezas directas e das constantes utilizadas, de modo a minimizar os erros sistemáticos e aumentar a exactidão. No final, a comparação do valor médio obtido com o valor da mesma grandeza tabelado, nas mesmas condições físicas (ou proveniente de outras experiências), permite estimar o desvio à exactidão do valor obtido, que pode ser estimado em percentagem como

$$\mathrm{desvio}(\%)=\left|\frac{\text { valor }_{\text {conhecido }}-\text { valor }_{\text {medido }}}{\text { valor }_{\text {conhecido }}}\right| \cdot 100$$

3. Porque existem sempre erros aleatórios, toda a medição é afectada de uma incerteza, que indica o grau de precisão. Obrigatoriamente em todos os resultados tem de se apresentar sempre o valor mais provável da grandeza, mais a respectiva estimativa numérica da incerteza. Exemplo: $v_{\mathrm{som(ar)}} =$ 343.5 ± 0.6 m/s
4. Quando se calculam grandezas indirectas a partir das medições directas, utilizando as equações físicas, as incertezas propagam-se, gerando uma incerteza do resultado final.

Veremos nas próximas secções como se pode, de uma forma simplificada, calcular e representar os valores mais prováveis para as grandezas directas e indirectas e as respectivas incertezas.

Resolução e sensibilidade

A resolução de um instrumento de medição é o menor intervalo mensurável com esse instrumento. É uma característica do seu desempenho em termos de qual o menor detalhe ou mudança que o instrumento consegue detectar. Por exemplo, na linguagem comum emprega-se o termo "resolução de um écran" para designar o nível de detalhe com que um aparelho reproduz imagens digitais. Em física experimental a resolução é estimada tendo em conta a menor escala ou algarismo exibido pelo instrumento e varia consoante este seja analógico ou digital:

Instrumentos analógicos	Considera-se que a resolução é metade da menor escala do instrumento. Exemplo: um voltímetro com uma escala graduada com divisões de 1 V tem uma resolução de 0,5 V, uma vez que a olho nu é possível perceber se uma dada medida está mais próxima de um traço (por exemplo, 10,0 V) ou do ponto médio entre dois traços (por exemplo 9,5 V).	A resolução desta escala analógica é 0,5 V.
Instrumentos digitais	Considera-se que a resolução corresponde à útima casa decimal exibida pelo instrumento, uma vez que essa é a incerteza sobre qual o arrendondamento que foi feito. Exemplo: Um voltímetro digital que mostre uma leitura de 12,6 V pode corresponder a um valor real contido entre 12,55 V e 12,65 V, pelo que a resolução é 0,1 V	A resolução desta escala digital é 0,1 V.

Por regra, a incerteza de uma medida (única) realizada com um instrumento é igual à sua resolução. No entanto, se a leitura do instrumento não permanecer constante – por exemplo, se a agulha de um voltímetro digital oscilar ou se os dígitos de um voltímetro digital variarem – a regra já não é válida e a incerteza deve ser estimada, usando bom senso, a partir do intervalo de variação.

A sensibilidade, por outro lado, é uma indicação do mínimo sinal detectável pelo instrumento, isto é, qual o valor mínimo que é necessário atingir para que uma leitura seja registada. Por exemplo, qual a menor massa que é necessário colocar no prato de uma balança para que esta registe o seu peso? Esse valor é a sua sensibilidade. Qual a menor divisão da escala da balança? Esse valor é a sua resolução.

Valor médio e incerteza nas medições experimentais

Normalmente, numa medição não se adquire apenas uma única medida de uma dada grandeza, mas sim um dado número $N$ que pode ser pequeno ou grande, consoante a importância de se conhecer o valor da grandeza com boa precisão e/ou exactidão. Tomando o valor médio de um conjunto de medidas, o efeito dos erros aleatórios pode ser atenuado, uma vez que os desvios de sinal oposto irão cancelar-se. No entanto, o efeito dos erros sistemáticos não é afectado pelo número de medições, permanencendo constantes. Para corrigir estes erros é preciso investigar as suas causas e corrigi-las.

Valor médio

A repetição de uma medição da variável $x$ nas mesmas condições experimentais conduz a uma distribuição aleatória de resultados em torno de um valor médio $\bar{x}$ (média aritmética), que pode ser considerado como o melhor valor obtido nesta medida. Por exemplo, para $N$ medidas $x_1,x_2,...$ da grandeza $x$ temos

$$\bar{x}=\frac{x_1+x_2+...+x_N}{N}=\frac{\sum_i x_i}{N}$$	$\quad\quad$ (Valor médio)

Num grande número de situações, esta repetição realizada $N$ vezes nas mesmas condições experimentais conduz a um valor médio que se aproxima do “verdadeiro” valor da grandeza à medida que $N$ aumenta. Para o cálculo da incerteza associada a esse valor médio devemos distinguir se se trata de uma grandeza directa ou indirecta.

Grandezas directas: determinação da incerteza

Devemos distinguir duas situações, dependendo do valor de $N$:

$N$ pequeno $(1<N<10)$

Para um número reduzido de medições, a incerteza deve ser estimada usando um majorante $\Delta x$, que será o maior desvio em relação ao valor médio. Define-se o desvio de cada medida individual $x_i$ como a diferença absoluta $\Delta x_i=|\bar{x}-x_i|$, pelo que a incerteza da medição é $$\Delta x=\max|\bar{x}-x_i|$$ Majorante dos desvios O resultado final neste caso pode apresentar-se numa das seguintes formas: $$\bar{x}\pm\Delta x$$ Resultado final, incerteza absoluta $$\bar{x}\pm\Delta x/\bar{x}$$ Resultado final, incerteza relativa No caso da incerteza relativa, o resultado é expresso em percentagem. Importante: se a incerteza calculada por este método for menor do que a incerteza intrínseca do instrumento (e.g. resolução da escala), a estimativa deve ser substituída por esta última.

$N$ grande $(N\gg 10)$

No caso de se dispôr de um número elevado de medidas é mais adequado empregar métodos estatísticos. Pode calcular-se o desvio padrão $s$, que exprime a dispersão dos resultados: $$s=\sqrt{\frac{\sum_i\left(x_i-\bar{x}\right)^2}{N-1}}$$ $\quad\quad$ (Desvio padrão) O melhor valor para a incerteza do valor médio $u$, é dado pelo desvio padrão da média, $u=s\sqrt{N}$, também chamado erro padrão ou erro padrão da média: $$u=\sqrt{\frac{\sum_i\left(x_i-\bar{x}\right)^2}{N(N-1)}}$$ $\quad\quad$ (Incerteza do valor médio) O resultado final neste caso (para um número elevado de determinações nas mesmas condições experimentais) deve apresentar-se como: $$\bar{x}\pm u$$ $\quad\quad$ (Apresentação de resultado final)

Exemplo. Nos trabalhos experimentais de LIFE tipicamente lida-se com um número pequeno de medições $(1<N<10)$, efectuadas manualmente. Considere-se o seguinte conjunto de cinco medidas de uma dada grandeza: # $t/\mathrm{s}\,(±0,1\,\mathrm{s})$ $\bar{t}/\mathrm{s}$ $\Delta t_i/\mathrm{s}$ 1 $31,0$ $31,4$ $|31,4-31,0|=0,4$ 2 $31,8$ $|31,4-31,8|=0,4$ 3 $30,6$ $|31,4-30,6|=0,8$ 4 $32,2$ $|31,4-32,2|=0,8$ 5 $31,4$ $|31,4-31,4|=0,0$ O maior dos desvios é 0,8 s, pelo que o resultado deve ser apresentado na forma $$t=31,4\pm 0,8\,\mathrm{s}$$ (incerteza absoluta) $$t=31,4\,\mathrm{s}\pm 2,6\%$$ (incerteza relativa)|}

Grandezas indirectas: determinação da incerteza

Para uma grandeza indirecta $F(X,Y,Z,…)$ sendo $X,Y,Z,…$ grandezas medidas directas, com incertezas que foram estimadas pelas equações acima como sendo $u_X , u_Y, u_Z$ pode estimar-se a incerteza $u_F$ da grandeza $F$ a partir das respectivas derivadas parciais:

$$u_F=\sqrt{\left(\frac{\partial F}{\partial X} u_X\right)^2+\left(\frac{\partial F}{\partial Y} u_Y\right)^2+\left(\frac{\partial F}{\partial Z} u_X\right)^2 \cdots}$$

Quando não é possível fazer uma análise estatística $(1<N<4)$, um majorante do erro da grandeza indirecta $\Delta F$ é calculável a partir de

$$\Delta F=\left|\frac{\partial F}{\partial X}\right| \Delta X+\left|\frac{\partial F}{\partial Y}\right| \Delta Y+\left|\frac{\partial F}{\partial Z}\right| \Delta Z$$

onde $\Delta X,\Delta Y,\Delta Z$, são as incertezas estimadas através dos majorantes dos erros das variáveis correspondentes. Caso estas incertezas sejam relevantes, as derivadas deverão ser calculadas por majoração.

Caso particular: para uma função racional (por ex. $F(X,Y,Z)=cte∙X^a Y^b Z^c$, com $a,b,c$ inteiros) o majorante do erro relativo pode ser dado simplesmente pela soma dos majorantes dos erros relativos das variáveis multiplicados pelos expoentes em valor absoluto:

$$\Delta F/F=|a|\cdot\frac{\Delta X}{X}+|b|\cdot\frac{\Delta Y}{Y}+|c|\cdot\frac{\Delta Z}{Z}$$

Exemplo. Consideremos a velocidade escalar $v=x/t$. É uma grandeza indirecta cujo medição envolve a medição das grandezas directas comprimento $\bar{x}\pm\Delta x$ e tempo $\bar{t}\pm\Delta t$. Para calcular a incerteza associada à velocidade, calculamos as respectivas derivadas parciais: $$\left|\frac{\partial v}{\partial x}\right|=\frac{1}{t}\quad\quad\left|\frac{\partial v}{\partial t}\right|=\frac{x}{t^2}$$ A majoração das derivadas faz-se calculando os seus valores na "pior" (maior valor numérico) situação, ou seja, maximizando os numeradores e minimizando os denominadores: $$\left|\frac{\partial v}{\partial x}\right|_\mathrm{maj}=\frac{1}{t-\Delta t}\quad\quad\left|\frac{\partial v}{\partial t}\right|_\mathrm{maj}=\frac{x+\Delta x}{(t-\Delta t)^2}$$ Para um número pequeno de medições obtemos a expressão para a incerteza do valor médio, $$\Delta v=\frac{1}{t-\Delta t} \Delta x+\frac{x+\Delta x}{(t-\Delta t)^2} \Delta t$$ Usando o "método expresso" do caso particular, uma vez que $v=x^1t^{-1}$ podemos escrever $$\frac{\Delta v}{v}=\frac{\Delta x}{x}+\frac{\Delta t}{t}\rightarrow\Delta v=\frac{\Delta x}{t}+\frac{x\Delta t}{t^2}$$ Majorando os quocientes, voltamos a obter a expressão calculada explicitamente pelas derivadas parciais. Assim, este método é muito mais prático e rápido. Aplicando a valores concretos, suponhamos que a tabela usada no exemplo anterior lista os tempos medidos para a duração de um percurso $\bar{x}=10,0±0,1\,\mathrm{mm}$, com duração média $\bar{t}=31,4±0,8\,\mathrm{s}$. A tabela abaixo indica as velocidades e respectivas incertezas, calculadas usando o método de propagação de incertezas. # $t/\mathrm{s}\,(±0,1\,\mathrm{s})$ $v/(\mathrm{mm/s})$ $\bar{v}/(\mathrm{mm/s})$ $\Delta\bar{v}/(\mathrm{mm/s})$ 1 $31,0$ $0,318$ 0,318 $$\frac{1}{31,4-0,8}\times 0,1+\frac{10+0,1}{(31,4-0,8)^2}\times 0,8$$ 2 $31,8$ $0,314$ 3 $30,6$ $0,327$ 4 $32,2$ $0,311$ 5 $31,4$ $0,311$

Combinação de resultados

Em muitas ocasiões pretende-se determinar o valor de uma dada grandeza física que é medida usando diferentes parâmetros, como forma de aumantar a gama de observações e minimizar as incertezas. Por exemplo, constantes físicas (velocidade da luz no vácuo, carga do electrão, constante de Planck, etc) ou propriedades materiais (índice de refracção de um vidro, período de oscilação de um pêndulo, etc) podem ser medidas usando diferentes parâmetros experimentais. O último passo consiste em combinar os resultados obtidos no conjunto de medições na forma de um "valor final" (e respectiva incerteza) para a experiência. Este processo pode ser feito através da média simples ou da média ponderada.

Média simples

Tendo os resultados de $N$ medições na forma $x_1±\Delta x_1,x_2±\Delta x_2,…x_N±\Delta x_N$, no caso de as incertezas $\Delta x_1, … \Delta x_N$ serem comparáveis podemos usar a média simples,

$$\bar{x}=\frac{x_1+x_2+...+x_N}{N}=\frac{\sum_i x_i}{N}$$	$\quad\quad$ (Média simples)

A incerteza do valor médio é dada pela regra de propagação de incertezas para medições independentes:

$$\Delta\bar{x}=\frac{\sqrt{\Delta x_1^2+\Delta x_2^2+...+\Delta x_N^2}}{N}=\frac{\sum_i \Delta x_i}{N}$$	$\quad\quad$ (Incerteza do valor médio)

Exemplo. Considere-se o seguinte conjunto de medições: $$x_1=1,0±0,1; x_2=1,1±0,2; x_3=1,2±0,2$$ . Como as incertezas são comparáveis, podemos aplicar as expressões acima e obter o valor da média simples e a sua incerteza: $$\bar{x}=\frac{1,0+1,1+1,2}{3}=1,1$$ $$\Delta\bar{x}=\frac{\sqrt{0,1^2+0,2^2+0,2^2}}{3}=0,1$$ O resultado final é apresentado na forma $x=1,1±0,1$

Média ponderada

No caso de as incertezas $\Delta x_1, … \Delta x_N$ serem significativamente diferentes deverá ter-se esse facto em consideração no cálculo do valor final; isto é, uma medição com uma incerteza pequena deverá ter mais peso no resultado final do que uma medição com uma incerteza grande. Nestes casos usamos a média ponderada, em que o peso de cada contribuição é dado por $w_i=1/\Delta x_i^2$:

$$\bar{x}=\frac{x_1 w_1^2+x_2 w_2^2+...+x_N w_N^2}{w_1^2+w_2^2+…+w_N^2}=\frac{\sum_i x_iw_i}{\sum_i w_i}$$	$\quad\quad$ (Média ponderada)

A incerteza do valor médio é dada pelo inverso da soma dos pesos:

$$\Delta\bar{x}=\sqrt{\frac{1}{w_1+w_2+…+w_N}}=\sqrt{\frac{1}{\sum_i w_i}}$$	$\quad\quad$ (Incerteza do valor médio)

Exemplo. Considere-se o mesmo conjunto de medições do exemplo acima. Comecemos por calcular os pesos: $$w_1=\frac{1}{0,1^2}=100\quad\quad w_2=\frac{1}{0,2^2}=25\quad\quad w_3=\frac{1}{0,3^2}=25$$ Assim, o resultado cujo incerteza é metade das outras tem um peso quatro vezes superior. Aplicando as expressões acima para a média ponderada obtemos: $$\bar{x}=\frac{1,0\times 100+1,1\times 25+1,2\times 25}{100+25+25}=1,05$$ $$\Delta\bar{x}=\sqrt{\frac{1}{100+25+25}}≈0,08$$ O resultado final é apresentado na forma $x=1,05±0,08$

Representação de resultados da medição de grandezas

Regras fundamentais:

• O resultado de qualquer medição deve ser apresentado na seguinte forma:

valor mais provável $$\pm$$ incerteza unidades físicas

• Como normalmente o valor da incerteza é determinado aproximadamente, em regra deverá ser indicado apenas com um ou dois algarismos significativos. Exemplo: 0,2 ou 0,21 é correcto, mas 0,213 não é.
• Por sua vez, o valor mais provável deve usar as mesmas casas decimais, arredondando-se o algarismo mais à direita. Exemplo: $x=2,25±0,15\, \mathrm{m}$ é uma representação correcta, mas $x=2,255±0,15\, \mathrm{m}$ não é.
• O uso de notação científica facilita o seguimento das regras acima e evita ambiguidades. Exemplo: em vez de apresentar o resultado na forma $x=2346±14\,\mathrm{m}$, deve-se apresentar na forma $x=(2,35\pm0,01)\times 10^2\,\mathrm{m}$
• Para as unidades físicas deverá usar-se o Sistema Internacional de Unidades. As unidades são apresentadas em tipo de letra romano (isto é, nem itálico, nem negrito) e separadas dos valores numéricos por um espaço. Exemplo: $x=(2,35\pm0,01)\times 10^2m$ apresenta os dois tipos de erros a evitar.
• Para o separador decimal, apesar de em Portugal vigorar o uso da vírgula (,) também é aceitável utilizar o ponto (.).

Bons exemplos	$R = 0.185\pm 0.030\, \mathrm{m}$ $Temp = 297.0\pm 0.5\,\mathrm{K}$ $v = 344.3\pm 0.4 \,\mathrm{m\cdot s}^{-1}$ $B = (5.92\pm 0.08)\times10^{-4}\,\mathrm{ T}$ $q/m = (1.77\pm 0.07)\times10^{11}\,\mathrm{ C\cdot kg}^{-1}$ $e = 0.050\pm 0.001 \,\mathrm{mm}$ ou $e=50\pm 1 \,\mathrm{m}$
Maus exemplos	$B = (5.9297887668888668898\pm 0.08) 10^{-4} \,\mathrm{ T}$ $Temp = 297\pm 0.0005$ $q/m = (1.8\pm 0.07789) 10^{11}\,\mathrm{ C\cdot kg}^{-1}$

Algarismos significativos

Com excepção do caso em que todos os números envolvidos são inteiros, não é possível representar o valor de uma grandeza com exactidão ilimitada. Diz-se que uma representação de um número tem $n$ algarismos significativos quando se admite um erro na casa decimal seguinte. Por exemplo:

• 1/7 = 0,14 tem dois algarismos significativos
• 1/30 = 0,0333 tem três algarismos significativos

Note-se que a posição da vírgula não afecta o número de a.s.

Regras

• Algarismos zero à esquerda não contam para o total de a.s. – exemplo: 0,00044 ( 2 a.s.)
• Algarismos zero à direita contam para o total de a.s. – exemplo: 12,00 (4 a.s.)
• Algarismos 1–9 e zeros entre eles são sempre a.s. – exemplo: 1203,4 (5 a.s.)
• Potências de dez são ambíguas, e devem ser representadas usando notação decimal – exemplo: 800 é ambíguo, 8,00$\times$10$^2$ é correcto (3 a.s.)
• As constantes têm um número arbitrário de a.s.

Soma e subtracção

O resultado deve manter o número de casas decimais do operando com o menor número de casas decimais

Exemplo: 105,4+0,2869+34,27 = 139,9569 = 140,0=1,400$\times$10$^2$

Multiplicação e divisão

O resultado deve manter o mesmo número de algarismos significativos do operando com o menor número de algarismos significativos.

Exemplo: 7,325$\times$8,14 = 59,6255 = 59,6

Raízes quadradas

O número de a.s. é igual ao de partida.

Exemplo: √92 = 9,59166 = 9,6

Incertezas nas representações gráficas: ajuste linear

A análise de resultados é frequentemente facilitada se se usarem representações gráficas das leis matemáticas que supostamente descrevem os fenómenos físicos em observação. O ajuste é particularmente simples se se tratar de uma lei linear, onde se pode fazer um ajuste visual.

De uma forma mais sistemática deve usar-se o método dos mínimos quadrados, que consiste na determinação analítica de qual a recta $y=a+b\cdot x_i$ que se desvia o menos possível do conjunto de pontos experimentais. Na sua forma mais simples, sendo $(x_i,y_i)$ as coordenadas dos $N$ pontos pretende-se determinar $(a,b)$ tal que

$$\chi^2=\sum_{i=0}^N\left(y_i-y\right)^2=\sum_{i=0}^N\left(y_i-a-b \cdot x_i\right)^2$$

seja mínimo. As condições de estacionariedade desta função $χ^2=F(a,b)$, dependente dos dois parâmetros $(a,b)$, podem ser descritas como $\partial(χ^2)/\partial a=0,\partial(χ^2)/\partial b=0$ e $\partial(χ^2)/\partial b^2=0$. As duas primeiras equações resultam em

$$\begin{aligned} & \sum_{i=0}^N\left(y_i-a-b x_i\right)=0 \Leftrightarrow \sum y_i-N a-b \sum x_i=0 \\ & \sum_{i=0}^N x_i\left(y_i-a-b x_i\right)=0 \Leftrightarrow \sum x_i y_i-a \sum x_i-b \sum x_i^2=0 \end{aligned}$$

A resolução deste sistema de duas equações permite obter os valores de $a$ e $b$:

$$a=\frac{\left(\sum x_i\right)^2 \sum y_i-\sum x_i \sum x_i y_i}{N \sum x_i^2-\left(\sum x_i\right)^2} \quad b=\frac{N \sum x_i y_i-\sum x_i \sum y_i}{N \sum x_i^2-\left(\sum x_i\right)^2}$$

A grande maioria dos programas de cálculo e as calculadoras científicas incorporam estas expressões para calcular os parâmetros de ajuste $a$ e $b$. Apenas os programas mais avançados para gráficos e análise de dados científicos (e.g Origin , Fitteia ou Qtiplot ) permitem também calcular as estimativas das incertezas $u_a$ e $u_b$.

Ajuste linear manual

É possível também obter um ajuste linear aproximado fazendo um traçado manual, com o rigor possível. Podemos usar como ponto de partida o ponto médio por onde passa a recta. Consideremos o sistema de duas equações acima; tomando a primeira e dividindo por $N$,

$$\frac{\sum y_i}{N}-a-b \frac{\sum x_i}{N}=0 \Leftrightarrow \bar{y}=a+b \bar{x}$$

em que $\bar{y}$ e $\bar{x}$ são respectivamente as médias de cada um dos conjuntos de valores. Daqui conclui-se que a recta que corresponde ao melhor ajuste passa pelo ponto médio $(\bar{x},\bar{y})$.

O passo seguinte consiste em traçar as rectas de maior $(y=a_1+b_1 x)$ e menor $(y=a_2+b_2 x)$ inclinação que, passando por este ponto, melhor se ajustam aos pontos medidos e suas incertezas. Por fim, a recta do melhor ajuste e o respectivo erro é obtida pela média desta duas rectas, de acordo com

$$a=\frac{a_1+a_2}{2} \quad \varepsilon_a=\frac{\left|a_1-a_2\right|}{2} \quad b=\frac{b_1+b_2}{2} \quad \varepsilon_b=\frac{\left|b_1-b_2\right|}{2}$$

Exemplo: Considere-se o conjunto de pontos da tabela abaixo e a sua representação no gráfico.

$x$ $y$ $\epsilon_y$ 1 2,4 1,1 2 5 0,5 3 6 0,7 4 7,7 0,8 5 10,4 0,9 6 10,5 1,5 Ponto médio: $\bar{x}=3,5; \bar{y}=7,0$

Traçamos duas rectas que passem pelo ponto médio (3,5; 7,0) e que correspondam visualmente (e com bom senso) ao maior e menor declive que contenham os pontos da medição.

Recta 1 (azul): $y=1,4x+2,05$ Recta 2 (verm.): $y=2,0x-0,05$

Por fim, calculam-se os coeficientes da recta que bissecta estas duas, dados pelas expressões acima.

Melhor ajuste: $$\begin{align} a &=\frac{2,05+(-0,05)}{2}=1,0\\ ϵ_a &=\frac{|2,05-(-0,05)|}{2}=1,05\\ b &=\frac{1,4+2}{2}=1,7\\ ϵ_b &=\frac{|1,4-2}{2}=0,35 \end{align}$$ Recta final (verde): $y=1,7x+1,0$

Como comparação, os valores calculados pelo método dos mínimos quadrados dão para a mesma recta o resultado $y=1,67x+1,16$.

Bibliografia

• John R. Taylor, An Introduction to Error Analysis: The Study of Uncertainties in Physical Measurements, University Science Books; 2nd edition (August 1, 1996)
• V. Thomsen, Precision and The Terminology of Measurement. The Physics Teacher Vol. 35, pp.15-17, Jan. 1997.
• Ifan Hughes and Thomas Hase, Measurements and their Uncertainties: A practical guide to modern error analysis, Oxford University Press (July 1, 2010)