Warum ist eine Projektionsmatrix einer orthogonalen Projektion symmetrisch?

16

Ich bin ziemlich neu darin, also hoffe ich, dass Sie mir vergeben, wenn die Frage naiv ist. (Kontext: Ich lerne Ökonometrie aus Davidson & MacKinnons Buch "Econometric Theory and Methods" und sie scheinen dies nicht zu erklären. Ich habe mir auch Luenbergers Optimierungsbuch angesehen, das sich mit Projektionen auf einem etwas fortgeschritteneren Niveau befasst, aber ohne glück).

Nehmen wir an, dass ich eine orthogonale Projektion zugeordnet ist Projektionsmatrix . Ich bin daran interessiert, jeden Vektor in in einen Unterraum projizieren .PPRnARn

Frage : Warum folgt , ist symmetrisch? Welches Lehrbuch könnte ich mir für dieses Ergebnis ansehen?P=PTP

weez13
quelle

Antworten:

13

Dies ist ein grundlegendes Ergebnis der linearen Algebra bei orthogonalen Projektionen. Ein relativ einfacher Ansatz lautet wie folgt. Wenn orthonormale Vektoren sind, die einen dimensionalen Unterraum überspannen , und die Matrix mit den als Spalten ist, dann ist Dies folgt direkt aus der Tatsache, dass die orthogonale Projektion von auf als orthonormale Basis von als berechnet werden kann Es folgt direkt aus der obigen Formel m A U n × p U i P = U U T . x A A m i = 1 u i u T i x . P 2 = P P T = P .u1,,ummAUn×pui

P=UUT.
xAA
i=1muiuiTx.
P2=P und dasPT=P.

Es ist auch möglich, ein anderes Argument anzugeben. Wenn eine Projektionsmatrix für eine orthogonale Projektion ist, dann per Definition für alle Folglich ist für alle . Dies zeigt, dass , woraus Px,yRn

PxyPy.

0=(Px)T(yPy)=xTPT(IP)y=xT(PTPTP)y
x,yRnPT=PTP
P=(PT)T=(PTP)T=PTP=PT.
NRH
quelle
Vielen Dank für Ihre aufschlussreichen Kommentare! Irgendwie hat mich der Wikipedia-Artikel, in dem etwas über die Selbstadjustiertheit des Projektionsoperators erwähnt wurde, abgelehnt, da Ihre Beweise nicht so schwierig sind. :) Übrigens, hast du einen Lieblingstext für lineare Algebra, der sich mit solchen Dingen befasst?
weez13
Das Buch über die elementare lineare Algebra, von dem ich am besten weiß, behandelt dies nicht. Die besten Referenzen, die ich kenne, sind fortgeschrittene Bücher zur Funktionsanalyse. Das Buch mit der linearen Algebra sieht gut aus, aber ich weiß es nicht.
NRH
Anmerkung: Die Antwort von NRH geht von . Das heißt, der einzige Fall, in dem (wie in der Gleichung ) beansprucht wird, liegt vor, wenn weil für jede lineare Abbildung und Vektor giltDies hat keinen wirklichen Einfluss auf das Ergebnis des Beweises, da die Implikation, dass ist, in beiden Fällen zutrifft, aber ich dachte, es wäre erwähnenswert. x=xT(Px)T=xPT(Px)T(yPy)=xPT(IP)yx=xTPx
(Px)T=xTPT.
PTPTP=0
Milan Mosse
1
@Milan Danke, dass du es bemerkt hast. Beachten Sie, dass für nur möglich ist, wenn , was uninteressant ist. Was einfach passiert ist, dass ein paar Transponierungen auf den in der vorletzten Zeile verloren gegangen sind . Ich habe die fehlenden Transponierten wiederhergestellt, um die Algebra zu korrigieren. x R n n = 1 xx=xTxRnn=1x
whuber
2

Ein Versuch der geometrischen Intuition ... Erinnern Sie sich daran:

  1. Eine symmetrische Matrix ist selbstadjunkt.
  2. Ein Skalarprodukt wird nur durch die Komponenten im gegenseitigen linearen Raum (und unabhängig von den orthogonalen Komponenten eines der Vektoren) bestimmt.

Was Sie "sehen" wollen, ist, dass eine Projektion sich selbst anpasst, also symmetrisch ist - nach (1). Warum ist das so? Betrachten Sie das Skalarprodukt eines Vektors mit der Projektion eines zweiten Vektors : . Nach (2) hängt das Produkt nur von den Komponenten von im Bereich der Projektion von . Das Produkt sollte also dasselbe sein wie und auch nach demselben Argument.A y x , A y x y A x , A y A x , y xAyx,AyxyAx,AyAx,y

Da selbstadjunkt ist, ist es symmetrisch.A

JohnRos
quelle
Danke vielmals! Bevor ich Ihren Kommentar las, war ich ziemlich verwirrt darüber, warum Selbstadjointness hier entscheidend ist. Jetzt habe ich eine Ahnung, danke!
weez13