Álgebra linear · Projeção ortogonal

Projeção ortogonal

A projeção ortogonal de um vetor b num subespaço é o ponto p do subespaço mais próximo de b. O erro e = b − p é perpendicular ao subespaço, e qualquer outro ponto q está mais longe: |b − q|² = |e|² + |p − q|².

A projeção ortogonal nesta visualização

A projeção ortogonal encontra o ponto de um subespaço mais próximo de um vetor dado. Aqui, o subespaço SS é um plano ou uma reta que passa pela origem O, e o vetor é b⃗\vec b. A sua projeção p⃗\vec p pertence a SS, e o erro e⃗=b⃗−p⃗\vec e = \vec b - \vec p é perpendicular a SS. Pode mover b⃗\vec b e alterar SS, comparar p⃗\vec p com qualquer outro ponto de SS e ver quando as sombras de b⃗\vec b sobre dois eixos somam p⃗\vec p.

Porque é que o ângulo reto dá o ponto mais próximo? Tome outro ponto qq em SS. O segmento de qq a p⃗\vec p está em SS, pelo que é perpendicular a e⃗\vec e. Os pontos b⃗\vec b, p⃗\vec p e qq formam um triângulo retângulo com o ângulo reto em p⃗\vec p, e Pitágoras dá

Qualquer outro ponto de S está mais longe
∣b⃗−q∣2=∣e⃗∣2+∣p⃗−q∣2|\vec b - q|^2 = |\vec e|^2 + |\vec p - q|^2

Assim, b⃗\vec b está mais longe de qq do que de p⃗\vec p, a menos que qq seja o próprio p⃗\vec p. O recíproco também é verdadeiro: se e⃗\vec e se inclinasse numa direção de SS, deslocar p⃗\vec p um pouco nessa direção aproximá-lo-ia de b⃗\vec b. Num plano, o ponto mais próximo é, portanto, aquele em que e⃗\vec e é perpendicular tanto a u⃗\vec u como a w⃗\vec w. Numa reta, basta ser perpendicular a u⃗\vec u.

Encontrar p⃗\vec p resume-se a estas condições de ângulo reto. Numa reta, p⃗\vec p é um múltiplo de u⃗\vec u, e a condição e⃗⋅u⃗=0\vec e\cdot\vec u = 0 determina esse múltiplo:

Projeção numa reta
p⃗=u⃗⋅b⃗u⃗⋅u⃗ u⃗\vec p = \frac{\vec u\cdot\vec b}{\vec u\cdot\vec u}\,\vec u

Num plano, escrevemos p⃗=c1u⃗+c2w⃗\vec p = c_1\vec u + c_2\vec w. As duas condições e⃗⋅u⃗=0\vec e\cdot\vec u = 0 e e⃗⋅w⃗=0\vec e\cdot\vec w = 0 tornam-se duas equações lineares para c1c_1 e c2c_2, chamadas equações normais:

As equações normais
(u⃗⋅u⃗) c1+(u⃗⋅w⃗) c2=u⃗⋅b⃗,(w⃗⋅u⃗) c1+(w⃗⋅w⃗) c2=w⃗⋅b⃗(\vec u\cdot\vec u)\,c_1 + (\vec u\cdot\vec w)\,c_2 = \vec u\cdot\vec b,\qquad (\vec w\cdot\vec u)\,c_1 + (\vec w\cdot\vec w)\,c_2 = \vec w\cdot\vec b

A leitura apresenta p⃗\vec p em coordenadas e como combinação de u⃗\vec u e w⃗\vec w. Arraste u⃗\vec u e w⃗\vec w dentro do mesmo plano: os coeficientes mudam, mas p⃗\vec p não. A projeção depende apenas do subespaço, não dos vetores usados para o descrever.

Quando u⃗\vec u e w⃗\vec w são perpendiculares, as duas equações deixam de depender uma da outra. Cada coeficiente depende apenas de um eixo, e p⃗\vec p é simplesmente a soma das sombras de b⃗\vec b sobre os dois eixos:

Só com eixos perpendiculares
p⃗=b⃗⋅u⃗u⃗⋅u⃗ u⃗+b⃗⋅w⃗w⃗⋅w⃗ w⃗\vec p = \frac{\vec b\cdot\vec u}{\vec u\cdot\vec u}\,\vec u + \frac{\vec b\cdot\vec w}{\vec w\cdot\vec w}\,\vec w

Com eixos oblíquos, este atalho geralmente falha. Cada sombra é calculada como se o outro eixo não existisse, pelo que a soma ss costuma ficar afastada de p⃗\vec p. Nesse caso, a leitura mostra que b⃗−s\vec b - s falha pelo menos uma das duas condições de perpendicularidade. Tornar w perpendicular a u substitui w⃗\vec w pela sua componente perpendicular a u⃗\vec u, um passo do processo de Gram–Schmidt. O plano continua a ser o mesmo, mas agora os seus eixos são perpendiculares, e as sombras somam p⃗\vec p para qualquer b⃗\vec b.

A diagonal que passa por (1,1,1)(1, 1, 1) dá um exemplo familiar. Projetar b⃗\vec b nessa reta substitui cada componente de b⃗\vec b pela média de todas elas, e ∣e⃗∣2|\vec e|^2 é a soma dos quadrados dos desvios em relação a essa média. A regressão pelo método dos mínimos quadrados usa o mesmo tipo de projeção.

Controlos

  • O ponto de b⃗\vec b move-o paralelamente a SS, e p⃗\vec p acompanha-o. O pequeno anel para lá do ponto move b⃗\vec b perpendicularmente a SS, e p⃗\vec p fica imóvel.
  • Com Plano selecionado, u⃗\vec u e w⃗\vec w só podem ser arrastados dentro do seu plano. Para mudar o plano, introduza novas coordenadas ou use as teclas de valores.

Ver também