線性代數(shù)學(xué)習(xí)總結(jié)-正交

概念

如果兩個(gè)向量的點(diǎn)積等于0呵曹,那么則說(shuō)這兩個(gè)向量是正交的脏款。(其實(shí)也就是垂直啦)

v\cdot w=0\quad 或者\(yùn)quad v^{T}w=0

正交矩陣的性質(zhì)有著很多實(shí)際的應(yīng)用。首先從4個(gè)子空間之間的正交性開(kāi)始吧退盯。

子空間正交性

這里是接著子空間更進(jìn)一步引述出的結(jié)論彼乌。

  • 行空間垂直與零空間垂直

很明顯,對(duì)于Ax=0的情況下渊迁,矩陣A的每一行都與x垂直

  • 列空間垂直與Left Null SpaceN(A^{T})垂直

與上面一致慰照,對(duì)于A^{T}y=0的情況下,矩陣A的每一列都與y垂直

Orthogonal Complements

對(duì)于一個(gè)子空間V而言琉朽,它的正交補(bǔ)碼包含每個(gè)與子空間V垂直的向量毒租,記作V^{\perp}

投影

這里投影就是指一個(gè)向量在另一個(gè)向量、平面或其他什么東西上的映射漓骚,跟字面意思一樣蝌衔。
首先考慮最簡(jiǎn)單的情況,投影到一條直線跟投影到一個(gè)平面上蝌蹂,為了節(jié)約紙張噩斟,我畫(huà)在一幅圖里了


projection on line/plane

如圖,向量b在向量a或者平面a上的投影為P孤个,如何求出P呢剃允?

假設(shè)P = x'a\quad與P垂直的向量為e
e = b - P=b - x'a
因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=e" alt="e" mathimg="1">與a垂直,有a\cdot e=a\cdot (b-x'a)=0
可以得到x'=\frac{a^{T}b}{a^{T}a}\qquad P=x'a=\frac{a^{T}b}{a^{T}a}a
對(duì)于平面而言齐鲤,求解方法一致斥废,得到x'=\frac{A^{T}b}{A^{T}A}\qquad P=Ax'=A(A^{T}A)^{-1}A^{T}b
得解

注意不要隨便亂消元哈,這不是自然數(shù)给郊。

接著我們考慮一種更復(fù)雜的情況牡肉,映射到子空間。
假設(shè)給出n個(gè)R^{m}空間內(nèi)的線性無(wú)關(guān)的向量組a_1, a_2,...a_n统锤,找出組合p=x'_1a_1+x'_2a_2+...x'_na_n,使得它與向量b的差值最小炭庙。
所謂差值最小饲窿,其實(shí)就是求b先這個(gè)子空間的投影了。那么我們只需要找出投影矩陣就可以了

對(duì)于p=x'_1a_1+x'_2a_2+...x'_na_n而言焕蹄,因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=b" alt="b" mathimg="1">的差值與a_i是垂直的逾雄,有
a_i^{T}(b-Ax')=0\quad即\quad\begin{bmatrix}a_1^{T}\\a_2^{T}\\...\\a_n^{T}\end{bmatrix}\begin{bmatrix}b-Ax'\\b-Ax'\\...\\b-Ax'\end{bmatrix} = 0
得到A^{T}(b-Ax')=0 \quad \rightarrow \quad A^{T}Ax'=A^{T}b
p=Ax'=A(A^{T}A)^{-1}A^{T}b\quad 因此,投影矩陣為\quad P=A(A^{T}A)^{-1}A^{T}

線性擬合

接下來(lái)的內(nèi)容比較偏應(yīng)用一點(diǎn)。對(duì)于方程Ax=b無(wú)解是一種比較常見(jiàn)的情況鸦泳,例如出現(xiàn)線性依賴(lài)等情況银锻。在這種情況下,我們需要求出一個(gè)最優(yōu)解做鹰,該怎么做呢徒仓?
所謂的最優(yōu)解,就是該解與b的差值的平方最小
是的誊垢,求投影就行了。如果差值的平方為0的話(huà)症见,則b就等于Ax喂走。
如下圖

Fitting a straight line

對(duì)于Ax=b\quad 可以表達(dá)為\quad C+D_{ti}=b_i\quad即\quad A=\begin{bmatrix}1&t_1\\1&t_2\\...&...\\1&t_n\end{bmatrix}\begin{bmatrix}C\\D\end{bmatrix}=b
我們需要求出C和D,根據(jù)公式A^{T}Ax'=A^{T}b谋作,有A^{T}A=\begin{bmatrix}1&1&...&1\\t_1&t_2&...&t_m\end{bmatrix}\begin{bmatrix}1&t_1\\1&t_2\\...&...\\1&t_m\end{bmatrix}=\begin{bmatrix}m&\sum_{}t_i\\\sum_{}t_i&\sum_{}t_i^{2}\end{bmatrix}
A^{T}b=\begin{bmatrix}1&1&...&1\\t_1&t_2&...&t_m\end{bmatrix}\begin{bmatrix}b_1\\b_2\\...\\b_m\end{bmatrix}=\begin{bmatrix}\sum_{}b_i\\\sum_{}t_ib_i\end{bmatrix}
那么當(dāng)A^{T}Ax'=A^{T}b的時(shí)候芋肠,\left||Ax-b\right||^2=e_1^2+e_2^2+...+e_m^2最小\begin{bmatrix}m&\sum_{}t_i\\\sum_{}t_i&\sum_{}t_i^{2}\end{bmatrix}\begin{bmatrix}C\\D\end{bmatrix}=\begin{bmatrix}\sum_{}b_i\\\sum_{}t_ib_i\end{bmatrix}求出C、D遵蚜,回代入公式帖池,得解。

線性擬合在應(yīng)用中有較多的使用吭净。不過(guò)容易出現(xiàn)的問(wèn)題就是干擾點(diǎn)睡汹,這些點(diǎn)離中心數(shù)據(jù)太遠(yuǎn),極大的影響了結(jié)果寂殉。通常在應(yīng)用中會(huì)選擇剔除這些點(diǎn)囚巴。

正交基和Gram-Schmidt

上面的過(guò)程可以看到,在求解的時(shí)候是非常麻煩的友扰,計(jì)算量很大彤叉。那么這里我們主要通過(guò)構(gòu)造正交矩陣來(lái)減少運(yùn)算量。
對(duì)于向量v_1,v_2,...v_n而言村怪,如果他們的點(diǎn)積d_i\cdot d_j=0(i\neq j)秽浇,我們則說(shuō)他們是正交的。如果向量的長(zhǎng)度為1甚负,那么我們則稱(chēng)之為正交單位向量柬焕,也叫標(biāo)準(zhǔn)正交。

q_i^Tq_j=\begin{Bmatrix}0\quad when\quad i\neq j\quad (orthogonal\quad vectors)\\1\quad when\quad i= j\quad (unit\quad vectors:\left||q_i\right||=1)\end{Bmatrix}

這樣我們可以很輕松的求解腊敲,因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=Q%5E%7BT%7DQ%3DI" alt="Q^{T}Q=I" mathimg="1">击喂,那么對(duì)于原來(lái)的公式A(A^{T}A)^{-1}A^{T}而言,可以化簡(jiǎn)為AA^{T}碰辅,極大的方便了運(yùn)算懂昂。

那么如何將原來(lái)的矩陣A變成正交單位矩陣Q呢?
如圖

transform

簡(jiǎn)單起見(jiàn)没宾,從低維開(kāi)始考慮凌彬,假設(shè)三維空間中沸柔,我們有向量組成的矩陣,如圖所示铲敛,要將它變?yōu)檎痪仃嚻鋵?shí)就很簡(jiǎn)單了褐澎。

  • 首先我們求出e,用B表示為B=b-\frac{A^{T}b}{A^{T}A}A
  • 其次伐蒋,我們用c減去其投影在B,a上的分量工三,既得到垂直于B,a平面的量了,如下圖
    c
  • 得到C=c-\frac{A^{T}c}{A^{T}A}A-\frac{B^{T}c}{B^{T}B}B

a,B,C組成正交單位矩陣先鱼。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末俭正,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子焙畔,更是在濱河造成了極大的恐慌掸读,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,496評(píng)論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件宏多,死亡現(xiàn)場(chǎng)離奇詭異儿惫,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)伸但,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,407評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門(mén)肾请,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人砌烁,你說(shuō)我怎么就攤上這事筐喳。” “怎么了函喉?”我有些...
    開(kāi)封第一講書(shū)人閱讀 162,632評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵避归,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我管呵,道長(zhǎng)梳毙,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,180評(píng)論 1 292
  • 正文 為了忘掉前任捐下,我火速辦了婚禮账锹,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘坷襟。我一直安慰自己奸柬,他們只是感情好蜂大,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,198評(píng)論 6 388
  • 文/花漫 我一把揭開(kāi)白布奥洼。 她就那樣靜靜地躺著,像睡著了一般宇整。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上桌粉,一...
    開(kāi)封第一講書(shū)人閱讀 51,165評(píng)論 1 299
  • 那天蒸绩,我揣著相機(jī)與錄音,去河邊找鬼铃肯。 笑死患亿,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的押逼。 我是一名探鬼主播步藕,決...
    沈念sama閱讀 40,052評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼挑格!你這毒婦竟也來(lái)了漱抓?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 38,910評(píng)論 0 274
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤恕齐,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后瞬逊,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體显歧,經(jīng)...
    沈念sama閱讀 45,324評(píng)論 1 310
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,542評(píng)論 2 332
  • 正文 我和宋清朗相戀三年确镊,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了士骤。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,711評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡蕾域,死狀恐怖拷肌,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情旨巷,我是刑警寧澤巨缘,帶...
    沈念sama閱讀 35,424評(píng)論 5 343
  • 正文 年R本政府宣布,位于F島的核電站采呐,受9級(jí)特大地震影響若锁,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜斧吐,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,017評(píng)論 3 326
  • 文/蒙蒙 一又固、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧煤率,春花似錦仰冠、人聲如沸。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,668評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至,卻和暖如春木张,著一層夾襖步出監(jiān)牢的瞬間众辨,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 32,823評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工舷礼, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留鹃彻,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,722評(píng)論 2 368
  • 正文 我出身青樓妻献,卻偏偏與公主長(zhǎng)得像蛛株,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子育拨,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,611評(píng)論 2 353

推薦閱讀更多精彩內(nèi)容