Algebra lineare · Proiezione ortogonale

Proiezione ortogonale

La proiezione ortogonale di un vettore b su un sottospazio è il punto p del sottospazio più vicino a b. L’errore e = b − p è perpendicolare al sottospazio e ogni altro punto q è più lontano: |b − q|² = |e|² + |p − q|².

La proiezione ortogonale in questa visualizzazione

La proiezione ortogonale trova il punto di un sottospazio più vicino a un vettore dato. Qui il sottospazio SS è un piano o una retta passante per l’origine O e il vettore è b⃗\vec b. La sua proiezione p⃗\vec p appartiene a SS, e l’errore e⃗=b⃗−p⃗\vec e = \vec b - \vec p è perpendicolare a SS. Puoi muovere b⃗\vec b e modificare SS, confrontare p⃗\vec p con qualsiasi altro punto di SS e vedere quando le ombre di b⃗\vec b su due assi sommate danno p⃗\vec p.

Perché l’angolo retto dà il punto più vicino? Prendi un altro punto qq in SS. Il segmento da qq a p⃗\vec p appartiene a SS, quindi è perpendicolare a e⃗\vec e. I punti b⃗\vec b, p⃗\vec p e qq formano un triangolo rettangolo con l’angolo retto in p⃗\vec p, e Pitagora dà

Ogni altro punto di S è più lontano
∣b⃗−q∣2=∣e⃗∣2+∣p⃗−q∣2|\vec b - q|^2 = |\vec e|^2 + |\vec p - q|^2

Quindi b⃗\vec b è più lontano da qq che da p⃗\vec p, a meno che qq non sia proprio p⃗\vec p. Vale anche il contrario: se e⃗\vec e fosse inclinato verso una direzione di SS, spostare leggermente p⃗\vec p in quella direzione lo avvicinerebbe a b⃗\vec b. In un piano, il punto più vicino è dunque quello in cui e⃗\vec e è perpendicolare sia a u⃗\vec u sia a w⃗\vec w. Su una retta, basta la perpendicolarità a u⃗\vec u.

Trovare p⃗\vec p si riduce a queste condizioni di angolo retto. Su una retta, p⃗\vec p è un multiplo di u⃗\vec u, e la condizione e⃗⋅u⃗=0\vec e\cdot\vec u = 0 determina il multiplo:

Proiezione su una retta
p⃗=u⃗⋅b⃗u⃗⋅u⃗ u⃗\vec p = \frac{\vec u\cdot\vec b}{\vec u\cdot\vec u}\,\vec u

In un piano, scriviamo p⃗=c1u⃗+c2w⃗\vec p = c_1\vec u + c_2\vec w. Le due condizioni e⃗⋅u⃗=0\vec e\cdot\vec u = 0 e e⃗⋅w⃗=0\vec e\cdot\vec w = 0 diventano due equazioni lineari per c1c_1 e c2c_2, dette equazioni normali:

Le equazioni normali
(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

La lettura mostra p⃗\vec p sia in coordinate sia come combinazione di u⃗\vec u e w⃗\vec w. Trascina u⃗\vec u e w⃗\vec w all’interno dello stesso piano: i coefficienti cambiano, ma p⃗\vec p no. La proiezione dipende solo dal sottospazio, non dai vettori usati per descriverlo.

Quando u⃗\vec u e w⃗\vec w sono perpendicolari, le due equazioni si disaccoppiano. Ogni coefficiente dipende da un solo asse e p⃗\vec p è semplicemente la somma delle ombre di b⃗\vec b sui due assi:

Solo con assi perpendicolari
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

Con assi obliqui, questa scorciatoia di solito non funziona. Ogni ombra è calcolata come se l’altro asse non esistesse, quindi la loro somma ss di solito non coincide con p⃗\vec p. In tal caso, la lettura mostra che b⃗−s\vec b - s non soddisfa almeno una delle due condizioni di perpendicolarità. Rendi w perpendicolare a u sostituisce w⃗\vec w con la sua componente perpendicolare a u⃗\vec u, un passo del procedimento di Gram–Schmidt. Il piano resta lo stesso, ma ora i suoi assi sono perpendicolari e le ombre sommate danno p⃗\vec p per ogni b⃗\vec b.

La diagonale passante per (1,1,1)(1, 1, 1) offre un esempio familiare. Proiettare b⃗\vec b su questa retta sostituisce ogni componente di b⃗\vec b con la media di tutte le componenti, e ∣e⃗∣2|\vec e|^2 è la somma dei quadrati degli scarti da tale media. La regressione ai minimi quadrati usa lo stesso tipo di proiezione.

Comandi

  • Il punto su b⃗\vec b lo muove parallelamente a SS, e p⃗\vec p segue. Il piccolo anello oltre il punto muove b⃗\vec b perpendicolarmente a SS, e p⃗\vec p resta fermo.
  • Con Piano selezionato, u⃗\vec u e w⃗\vec w si possono trascinare solo all’interno del loro piano. Per cambiare piano, inserisci nuove coordinate o usa i tasti dei valori.

Vedi anche