Perché una matrice di proiezione di una proiezione ortogonale è simmetrica?


16

Sono abbastanza nuovo per questo, quindi spero che tu mi perdoni se la domanda è ingenua. (Contesto: sto imparando l'econometria dal libro di Davidson e MacKinnon "Teoria e metodi econometrici" , e non sembrano spiegarlo; ho anche guardato il libro di ottimizzazione di Luenberger che tratta le proiezioni a un livello un po 'più avanzato, ma senza fortuna).

Supponiamo che ho una proiezione ortogonale con è associato proiezione matrice . Sono interessato a proiettare ogni vettore in in un sottospazio . P R n A R nPPRnARn

Domanda : perché segue che , cioè è simmetrico? Quale libro di testo posso consultare per questo risultato?T PP=PTP


Risposte:


13

Questo è un risultato fondamentale dell'algebra lineare sulle proiezioni ortogonali. Un approccio relativamente semplice è il seguente. Se sono vettori ortonormali che abbracciano un sottospazio m -dimensionale A , e U è la matrice con le come colonne, quindi Ciò deriva direttamente dal fatto che la proiezione ortogonale di su può essere calcolata in termini di base ortonormale di come u1,...,ummUNUu i P = U U T . x A A m i = 1 u i u T i x . P 2 = P P T = P .n×puio

P=UUT.
XUNUN
Σio=1muiouioTX.
Segue direttamente dalla formula sopra che e cheP2=PPT=P.

È anche possibile fornire un argomento diverso. Se è una matrice di proiezione per una proiezione ortogonale, allora, per definizione, per tutti Di conseguenza, per tutti . Ciò mostra che , da cui x , y R n P x y - P y . 0 = ( P x ) T ( y - P y ) = x T P T ( I - P ) y = x T ( P T - P T P ) y x , y R n P T = PPX,yRn

PXy-Py.

0=(PX)T(y-Py)=XTPT(io-P)y=XT(PT-PTP)y
X,yRnPT=PTP
P=(PT)T=(PTP)T=PTP=PT.

Grazie per i tuoi commenti perspicaci! In qualche modo l'articolo di Wikipedia, che menzionava qualcosa sull'auto-aggiunta dell'operatore della proiezione, mi ha sconcertato, poiché le tue prove non sono così difficili. :) A proposito, hai un testo di algebra lineare preferito che si occupa di questo tipo di cose?
weez13,

Il libro elementare di algebra lineare che conosco il meglio non copre questo. I migliori riferimenti che conosco sono libri avanzati sull'analisi funzionale. L' algebra lineare fatta nel libro giusto sembra buona, ma non lo so.
NRH

X=XT(PX)T=XPT(PX)T(y-Py)=XPT(io-P)yX=XTPX
(PX)T=XTPT.
PT-PTP=0

1
x=xTxRnn=1x

2

Un tentativo di intuizione geometrica ... Ricordiamo che:

  1. Una matrice simmetrica è autoaggiunta.
  2. Un prodotto scalare è determinato solo dai componenti nello spazio lineare reciproco (e indipendente dai componenti ortogonali di uno qualsiasi dei vettori).

XUNyX,UNyXyUNX,UNyUNX,y

UN


Molte grazie! Prima di leggere il tuo commento, ero piuttosto confuso sul perché l'auto-aggiunta è cruciale qui. Ora ho qualche indizio, grazie!
weez13,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.