Чому матриця проекцій ортогональної проекції симетрична?


16

Я зовсім новачок у цьому, тому сподіваюся, що ви пробачте мене, якщо питання буде наївним. (Контекст: я вивчаю економетрику з книги "Економетрична теорія та методи" Девідсона та МакКіннона , і, схоже, це не пояснюють; я також роздивився оптимізаційну книгу Луенбергера, яка займається прогнозами на трохи більш просунутому рівні, але без удачі).

Припустимо , що у мене є ортогональну проекцію з асоційована проекційної матриці . Мені цікаво проектувати кожен вектор у у якийсь підпростір .PPRnARn

Запитання : чому випливає, що , тобто симетричний? Який підручник я міг би подивитися на цей результат?T PP=PTP


Відповіді:


13

Це фундаментальний результат лінійної алгебри на ортогональних проекціях. Порівняно простий підхід полягає в наступному. Якщо є ортонормальними векторами, що охоплюють -вимірний підпростір , а - матриця із стовпцями , то Це випливає безпосередньо з того факту, що ортогональну проекцію на можна обчислити через ортонормальну основу як З формули вище випливає, що це м U п × р у я Р = U U Т . x A A m i = 1 u i u T i x . Р 2 = Р Р Т = Р .u1,,ummAUn×pui

P=UUT.
xAA
i=1muiuiTx.
P2=P і цеPT=P.

Можна також навести інший аргумент. Якщо - матриця проекцій для ортогональної проекції, то, за визначенням, для всіх Отже, для всіх . Це показує, що , звідки 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

PxyPy.

0=(Px)T(yPy)=xTPT(IP)y=xT(PTPTP)y
x,yRnР = ( Р Т ) Т = ( Р Т Р ) Т = Р Т Р = Р Т .PT=PTP
P=(PT)T=(PTP)T=PTP=PT.

Дякуємо за Ваш проникливий коментар! Якось стаття у Вікіпедії, в якій згадувалося щось про самозалежність оператора проекції, мене відкинула, оскільки ваші докази не такі складні. :) До речі, у вас є улюблений текст лінійної алгебри, який стосується подібних матеріалів?
weez13

Елементарна книга лінійної алгебри, яку я знаю найкраще, не охоплює цього. Найкращі мені відомості - це передові книги з функціонального аналізу. Лінійна алгебра зробила правильну книга виглядає добре, але я не знаю.
NRH

Примітка: відповідь NRH передбачає , що . Тобто єдиний випадок, коли (як стверджується в рівності ), коли оскільки для будь-якої лінійної карти і вектора ,Це насправді не впливає на результат доказування, оскільки це означає, що має значення в обох випадках, але я вважав, що це варто згадати. x=xT(Px)T=xPT(Px)T(yPy)=xPT(IP)yx=xTPx
(Px)T=xTPT.
PTPTP=0
Мілан Моссе

1
@Milan Дякую, що помітили. Зауважте, що для можливе лише тоді, коли , що нецікаво. Просто сталося те, що на передостанній лінії було втрачено декілька транспозитів . Я відновив відсутні місця транспозиції, щоб зробити алгебру правильною. x=xTxRnn=1x
whuber

2

Спроба геометричної інтуїції ... Нагадаємо, що:

  1. Симетрична матриця є самосуміжним.
  2. Скалярний добуток визначається лише компонентами у взаємному лінійному просторі (і незалежними від ортогональних компонентів будь-якого з векторів).

Те, що ви хочете "побачити", - це те, що проекція є самопримикаючою, таким чином, симетричною - наступною (1). Чому це так? Розглянемо скалярний добуток вектора з проекцією другого вектора : . Дотримуючись (2), добуток залежатиме лише від складових в проміжку проекції . Отже, продукт повинен бути таким самим, як , а також після того ж аргументу.xAyx,AyxyAx,AyAx,y

Оскільки є самосуміжним - він симетричний.A


Дуже дякую! Перш ніж прочитати ваш коментар, я сильно розгубився, чому саме тут самосуд є вирішальним. Тепер у мене є якийсь підказки, дякую!
weez13
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.