Módulo 8 - Quadrados Mínimos

Como já vimos no Módulo 7, uma forma de se trabalhar com uma função definida por uma
tabela de valores é a interpolação. Contudo, a interpolação pode não ser aconselhável
quando:

- É preciso obter um valor aproximado da função em algum ponto fora do intervalo de

tabelamento (extrapolação).

- Os valores tabelados são resultado de experimentos físicos, pois estes valores poderão conter
erros inerentes que, em geral, não são previsíveis.

Surge então a necessidade de ajustar estas funções tabeladas em uma outra função que seja
uma “boa aproximação” e que nos permita “extrapolar” com certa margem de segurança.

Assim, o objetivo deste processo é aproximar uma função f por outra função g, escolhida de
uma família de funções em duas situações distintas:

Domínio discreto: quando a função f é dada por uma tabela de valores

Domínio contínuo: quando a função f é dada por sua forma analítica.

Vamos trabalhar no domínio discreto onde queremos encontrar um ajuste de curvas quando
temos uma tabela de pontos. Sendo x1 , x2 , x3 , ... , xm 𝜖 [ a ,b ], considere a tabela:

xi x1 x2 x3 ... xm
f(xi) f(x1) f(x2) f(x3) ... f(xm)

O problema do ajuste de curvas consiste em: “escolhidas” n funções contínuas g1(x), g2(x),
g3(x), ... , gn (x), contínuas em [a , b], obter n constantes 𝛼1 , 𝛼2 , 𝛼3 ,..., 𝛼𝑛 tais que a função
g(x) se aproxime ao máximo de f(x), sendo g(x) dada por:

𝑔(𝑥) = 𝛼1 . 𝑔1 (𝑥) + 𝛼2 . 𝑔2 (𝑥) + 𝛼3 . 𝑔3 (𝑥) + ⋯ + 𝛼𝑛 . 𝑔𝑛 (𝑥)

Surge então a primeira pergunta: Como escolher as funções contínuas g1(x), g2(x), g3 (x),...,
gn (x)?

Esta escolha pode ser feita observando o gráfico dos pontos tabelados ou baseando-se em
fundamentos teóricos do experimento que forneceu a tabela.

No Módulo 6 aplicamos o método dos mínimos quadrados escolhendo g(x) retas como as
funções de aproximação de f(x). Hoje escolheremos polinômios.

Lembrando: seja o desvio xk dado por dk:

𝑑𝑘 = 𝑓(𝑥𝑘 ) − 𝑔(𝑥𝑘 )

O método dos mínimos quadrados consiste em escolher os coeficientes 𝛼1 , 𝛼2 , 𝛼3 ,..., 𝛼𝑛 de

tal forma que a soma dos quadrados dos desvios seja mínima, isto é:

∑ 𝑑𝑘 2 = ∑(𝑓(𝑥𝑘 ) − 𝑔(𝑥𝑘 ))2

Vejamos então como aplicar esse método e obter uma função de aproximação utilizando um
polinômio como ajuste.
Quadrados Mínimos - Ajuste a um polinômio

O objetivo, agora, é mostrar como ajustar os pontos de uma tabela com n pontos a uma
função polinomial de grau m

𝑃(𝑥) = 𝑎0 . 𝑥 0 + 𝑎1 . 𝑥 1 + 𝑎2 . 𝑥 2 + 𝑎3 . 𝑥 3 + ⋯ + 𝑎𝑚 . 𝑥 𝑚
Onde:
𝑚≤𝑛−1

Neste caso, a soma dos quadrados das distâncias de 𝑦𝑖 à 𝑃(𝑥𝑖 )é dada por:

𝑞 = ∑(𝑦𝑖 − 𝑃(𝑥𝑖 ))2

E depende de m + 1 parâmetros 𝑎0 , 𝑎1 , 𝑎2 , … , 𝑎𝑚

Para minimizar essa função temos que satisfazer às m + 1 condições a seguir:

∆𝑞
= 0 ∀𝑖 = 0,1,2,3, … 𝑚
∆𝑎𝑖

A qual fornece um sistema de m + 1 equações normais.

No caso de a função polinomial ser quadrática, isto é:

𝑃(𝑥) = 𝑎0 . 𝑥 0 + 𝑎1 . 𝑥 1 + 𝑎2 . 𝑥 2

As equações normais são:

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 0 . 𝑦𝑖
0 1 2

𝑖=1 𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 1 . 𝑦𝑖
1 2 3

𝑖=1 𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 2 . 𝑦𝑖
2 3 4

𝑖=1 𝑖=1 𝑖=1 𝑖=1

Podemos simplificar essas equações facilmente:

𝑛 𝑛 𝑛

𝑛. 𝑎0 + (∑ 𝑥𝑖 1 ) . 𝑎1 + (∑ 𝑥𝑖 2 ) . 𝑎2 = ∑ 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛
1 2 3
(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 . 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 2 . 𝑦𝑖
2 3 4

𝑖=1 𝑖=1 𝑖=1 𝑖=1

Observe que este sistema é simétrico. Para resolvê-lo, isto é, para encontrar as incógnitas
𝑎0 , 𝑎1 , 𝑎2 , … , 𝑎𝑚 podemos aplicar qualquer um dos métodos numéricos apresentados em
Tópicos de Matemática Aplicada.

Vejamos a aplicação do Método dos Mínimos Quadrados em exercícios:

Exercício 1: considerando os dados apresentados na tabela a seguir. Para encontrar o polinômio

de mínimos quadrados que aproxima esses dados, polinômio de grau dois do tipo:

𝑷(𝒙) = 𝒂𝟎 . 𝒙𝟎 + 𝒂𝟏 . 𝒙𝟏 + 𝒂𝟐 . 𝒙𝟐

Para isso é necessário montar um sistema. Encontre o sistema que melhor representa as
equações necessárias para encontra o polinômio y(x).

x 1 2 3 4
y -1 4 11 20

Solução:

n=4
𝑛 𝑛 𝑛
1 2
𝑛. 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛
1 2 3
(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 . 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 2 . 𝑦𝑖
2 3 4

𝑖=1 𝑖=1 𝑖=1 𝑖=1

Precisamos então calcular os somatórios. Faremos isso construindo uma tabela com os
exponenciais necessários:

i 𝒙𝒊 𝒙𝒊 𝟐 𝒙𝒊 𝟑 𝒙𝒊 𝟒 𝒚𝒊 𝒙𝒊 . 𝒚𝒊 𝒙𝒊 𝟐 . 𝒚𝒊
2 3 4 1. −1 = −1 1. −1 = −1
1 1 1 =1 1 =1 1 =1 -1
2 2 22 = 4 23 = 8 24 = 16 4 2.4 = 8 4.4 = 16
3 3 32 = 9 33 = 27 34 = 81 11 3.11 = 33 9.11 = 99
4 4 42 = 16 43 = 64 44 = 256 20 4.20 = 80 16.20 = 320
Soma 10 30 100 354 34 120 434

Montando as equações, temos:

4. 𝑎0 + 10. 𝑎1 + 30. 𝑎2 = 34

10. 𝑎0 + 30. 𝑎1 + 100. 𝑎2 = 120

30. 𝑎0 + 100. 𝑎1 + 354. 𝑎2 = 434

Exercício 2: considerando os dados apresentados na tabela a seguir. Para encontrar o polinômio
de mínimos quadrados que aproxima esses dados, polinômio de grau dois do tipo:

𝑷(𝒙) = 𝒂𝟎 . 𝒙𝟎 + 𝒂𝟏 . 𝒙𝟏 + 𝒂𝟐 . 𝒙𝟐

Para isso é necessário montar um sistema. Encontre o sistema que melhor representa as
equações necessárias para encontra o polinômio y(x).

x -1 0 1 2
y 1 1 3 7

Solução:

n=4
𝑛 𝑛 𝑛
1 2
𝑛. 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛
1 2 3
(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 . 𝑦𝑖
𝑖=1 𝑖=1 𝑖=1 𝑖=1

𝑛 𝑛 𝑛 𝑛

(∑ 𝑥𝑖 ) . 𝑎0 + (∑ 𝑥𝑖 ) . 𝑎1 + (∑ 𝑥𝑖 ) . 𝑎2 = ∑ 𝑥𝑖 2 . 𝑦𝑖
2 3 4

𝑖=1 𝑖=1 𝑖=1 𝑖=1

Precisamos então calcular os somatórios. Faremos isso construindo uma tabela com os
exponenciais necessários:

i 𝒙𝒊 𝒙𝒊 𝟐 𝒙𝒊 𝟑 𝒙𝒊 𝟒 𝒚𝒊 𝒙𝒊 . 𝒚𝒊 𝒙𝒊 𝟐 . 𝒚𝒊
1 -1 −12 = 1 3
−1 = −1 −14 = 1 1 −1.1 = −1 1.1 = 1
2 0 02 = 0 03 = 0 04 = 0 1 0.1 = 0 0.1 = 0
3 1 12 = 1 13 = 1 14 = 1 3 1.3 = 3 1.3 = 3
4 2 22 = 4 23 = 8 24 = 16 7 2.7 = 14 4.7 = 28
Soma 2 6 8 18 12 16 32

Montando as equações, temos:

4. 𝑎0 + 2. 𝑎1 + 6. 𝑎2 = 12

2. 𝑎0 + 6. 𝑎1 + 8. 𝑎2 = 16

6. 𝑎0 + 8. 𝑎1 + 18. 𝑎2 = 32
Exercício 3: considerando os dados apresentados na tabela a seguir. Para encontrar o polinômio
de mínimos quadrados que aproxima esses dados, polinômio de grau dois do tipo:

𝑷(𝒙) = 𝒂𝟎 . 𝒙𝟎 + 𝒂𝟏 . 𝒙𝟏 + 𝒂𝟐 . 𝒙𝟐

Para isso é necessário montar um sistema. Encontre o sistema que melhor representa as
equações necessárias para encontra o polinômio y(x).

x -1 0 1
y 5 -1 3

Resposta:
3. 𝑎0 + 0. 𝑎1 + 2. 𝑎2 = 7

0. 𝑎0 + 2. 𝑎1 + 0. 𝑎2 = −2

2. 𝑎0 + 0. 𝑎1 + 2. 𝑎2 = 8