编程与线性代数
先來了解線性代數是什么東東？
在大學數學學科中，線性代數是最為抽象的一門課，從初等數學到線性代數的思維跨度比微積分和概率統計要大得多。很多人學過以后一直停留在知其然不知 其所以然的階段，若干年之后接觸圖形編程或機器學習等領域才發現線性代數的應用無處不在，但又苦于不能很好地理解和掌握。的確，多數人很容易理解初等數學 的各種概念，函數、方程、數列一切都那么的自然，但是一進入線性代數的世界就好像來到了另一個陌生的世界，在各種奇怪的符號和運算里迷失了。
我在初接觸線性代數的時候簡直感覺這是一門天外飛仙的學科，一個疑問在我腦子里浮現出來：線性代數到底是一種客觀的自然規律還是人為的設計？
如果看到這個問題，你的反應是“這還用問，數學當然是客觀的自然規律了”，我一點兒都不覺得奇怪，我自己也曾這樣認為。從中學的初等數學和初等物理 一路走來，很少人去懷疑一門數學學科是不是自然規律，當我學習微積分、概率統計時也從來沒有懷疑過，唯獨線性代數讓我產生了懷疑，因為它的各種符號和運算 規則太抽象太奇怪，完全對應不到生活經驗。所以，我還真要感謝線性代數，它引發了我去思考一門數學學科的本質。其實，不止是學生，包括很多數學老師都不清 楚線性代數到底是什么、有什么用，不僅國內如此，在國外也是這樣，國內的孟巖寫過《理解矩陣》，國外的Sheldon Axler教授寫過《線性代數應該這樣學》，但都還沒有從根本上講清楚線性代數的來龍去脈。對于我自己來講，讀大學的時候沒有學懂線性代數，反而是后來從編程的角度理解了它。很多人說數學好可以幫助編程，我恰好反過來了，對程序的理解幫助了我理解數學。
本文的目標讀者是程序員，下面我就帶各位做一次程序員在線性代數世界的深度歷險！既然是程序員，在進入線性代數的領域之前，我們不妨先從考察一番程序世界，請思考這樣一個問題：
計算機里面有匯編、C/C++、Java、Python等通用語言，還有Makefile、CSS、SQL等DSL，這些語言是一種客觀的自然規律還是人為的設計呢？
為什么要問這樣一個看起來很蠢的問題呢？因為它的答案顯而易見，大家對天天使用的程序語言的認識一定勝過抽象的線性代數，很顯然程序語言雖然包含了 內在的邏輯，但它們本質上都是人為的設計。所有程序語言的共同性在于：建立了一套模型，定義了一套語法，并將每種語法映射到特定的語義。程序員和語言實現 者之間遵守語言契約：程序員保證代碼符合語言的語法，編譯器/解釋器保證代碼執行的結果符合語法相應的語義。比如，C++規定用new A()語法在堆上構造對象A，你這樣寫了C++就必須保證相應的執行效果，在堆上分配內存并調用A的構造函數，否則就是編譯器違背語言契約。
從應用的角度，我們能不能把線性代數視為一門程序語言呢？答案是肯定的，我們可以用語言契約作為標準來試試。假設你有一個圖像，你想把它旋轉60 度，再沿x軸方向拉伸2倍；線性代數告訴你，“行！你按我的語法構造一個矩陣，再按矩陣乘法規則去乘你的圖像，我保證結果就是你想要的”。
實際上，線性代數和SQL這樣的DSL非常相似，下面來作一些類比：
模型和語義：SQL是在低級語言之上建立了關系模型，核心語義是關系和關系運算；線性代數在初等數學之上建立了向量模型，核心語義是向量和線性變換
語法：SQL為每種語義定義了相應的語法，如select, where, join等；線性代數也定義了向量、矩陣、矩陣乘法等語義概念相應的語法
編譯/解釋：SQL可以被編譯/解釋為C語言；線性代數相關概念和運算規則可以由初等數學知識來解釋
實現：我們可以在MySQL、Oracle等關系數據庫上進行SQL編程；我們也可以在MATLAB、Mathematica等數學軟件上進行線性代數編程
所以，從應用的角度看，線性代數是一種人為設計的領域特定語言(DSL)，它建立了一套模型并通過符號系統完成語法和語義的映射。實際上，向量、矩陣、運算規則的語法和語義都是人為的設計，這和一門語言中的各種概念性質相同，它是一種創造，但是前提是必須滿足語言契約。
為什么要有線性代數？
可能有人對把線性代數當成一門DSL不放心，我給你一個矩陣，你就把我的圖形旋轉了60度沿x軸拉伸了2倍，我總感覺不踏實啊，我都不知道你“底 層”是怎么做！其實，這就像有的程序員用高級語言不踏實，覺得底層才是程序的本質，老是想知道這句話編譯成匯編是什么樣？那個操作又分配了多少內存？別人 在Shell里直接敲一個wget命令就能取下一個網頁，他非要用C語言花幾十分鐘來寫一堆代碼才踏實。其實，所謂底層和上層只是一種習慣性的說法，并不 是誰比誰更本質。程序的編譯和解釋本質上是不同模型間的語義映射，通常情況下是高級語言映射為低級語言，但是完全也可以把方向反過來。Fabrice Bellard用JavaScript寫了一個虛擬機，把Linux跑在JavaScript虛擬機上，這就是把機器模型往JavaScript模型上映射。
建立新模型肯定依賴于現有的模型，但這是建模的手段而不是目的，任何一種新模型的目的都為了更簡單地分析和解決某一類問題。線性代數在建立的時候，它的各種概念和運算規則依賴于初等數學的知識，但是一旦建立起來這層抽象模型之后，我們就應該習慣于直接利用高層次的抽象模型去分析和解決問題。
說到線性代數是為了比初等數學更容易地分析和解決問題，下面我們通過一個例子來實際感受一下它的好處：
給定三角形的頂點(x1, y1), (x2, y2), (x3, y3)，求三角形的面積。
初等數學中三角形面積最著名的計算公式是area = 1/2 * base * height ，當三角形有一條邊恰好在坐標軸上時我們就很容易算出它的面積。但是，假如同樣一個三角形我們把坐標軸旋轉一下，讓它的邊不在坐標軸上，怎么辦？我們還能得到它的底和高嗎？答案肯定是可以的，但是就明顯復雜了，而且還要分很多種情況去分別討論。
相反，如果我們用線性代數知識來解決這個問題就非常輕松。在線性代數中兩個向量a，b的叉積(Cross Product)是一個向量，其方向與a，b垂直，其大小等于a，b構成的平行四邊形的面積:
我們可以把三角形的邊視為向量，所以三角形的面積等于兩個邊向量的叉積除以二的絕對值：
area = abs(1/2 * cross_product((x2 - x1, y2 - y1), (x3 - x1, y3 - y1)))
注：abs表示取絕對值，cross_product表示兩個向量的叉積。
這樣一個在初等數學里面有點兒小難的問題在線性代數中瞬間搞定！可能有人會說，你直接基于叉積來做，當然簡單了，但是叉積本身不是也挺復雜的嗎？你把它展開試試看呢？是的，模型的作用就是把一部分復雜性隱藏到模型中，使得模型的使用者可以更加簡單地解決問題。曾經有人質疑C++太復雜，C++之父Bjarne Stroustrup這樣回答：
Complexity will go somewhere: if not the language then the application code.
在特定環境下，問題的復雜性是由其本質決定的，C++把一部分的復雜性納入了語言和標準庫，目的是使得應用程序更為簡單。當然，并非所有場合C++ 都使得問題更加簡單，但是從原理上講，C++的復雜性是有道理的。除了C++，Java、SQL、CSS等各種語言和框架莫不如是，想象一下，如果不使用 數據庫，動不動就自己去做數據存儲和管理是多么復雜啊！這樣我們就不難理解為什么線性代數要定義叉積這樣奇怪的運算了，它和C++把很多常用的算法和容器 納入STL是同一道理。同樣的，甚至你還可以在線性代數中定義自己想要的運算拿來復用。所以，數學一點兒不死板，它和程序一樣是活活潑潑的，你理解了它的 來龍去脈就能駕馭自如。說到這里，我們就順便回答一個很常見的疑惑：
線性代數的點積、叉積還有矩陣運算都很奇怪，為什么要定義這些運算呢？它們的定義又為什么是這個樣子呢？
其實，和程序復用一樣，線性代數定義點積、叉積和矩陣運算是因為它們的應用非常廣，有很大的復用價值，可以作為我們分析和解決問題的基礎。比如，很多問題都涉及到一個向量到另一個向量的投影或是求兩個向量的夾角，那么就會考慮專門定義點積(Dot Product)這個運算：
點積概念的提出屬于設計，有發揮創造的余地；一旦設計定了，具體公式就不能隨意發揮了，必須符合邏輯，保證它映射到初等數學模型的正確性。這就像一門高級語言可以定義很多概念，什么高階函數、閉包等等，但是它必須保證映射到底層實現時在執行產生的效果符合其定義的規范。
線性代數好在哪里？
上面說了，線性代數是一種高層次抽象模型，我們可以采用學習一門程序語言的方法去學習它的語法和語義，但是這一認識不只針對線性代數，它是對每一門數學學科通用的，可能有人會有疑問：
微積分、概率論也是高層次抽象，那么線性代數這種高層次抽象的特點在哪里呢？
這就問到了根本上，線性代數的核心：向量模型。我們在初等數學中學習的坐標系屬于笛卡爾所提出的解析模型，這個 模型很有用，但同時也有很大的缺點。坐標系是人為加上的虛擬參考系，但是我們要解決的問題，比如求面積，圖形旋轉、拉伸等應用都是和坐標系無關的，建立一 個虛擬的坐標系往往無助于解決問題，剛才三角形面積的例子就是這樣。
向量模型很好地克服了解析模型的缺點，如果說解析模型代表了某種“絕對性”的世界觀，那么向量模型就代表了某種“相對性”的世界觀，我推薦把向量模型和解析模型看作對立的兩種模型。
向量模型中定義了向量和標量的概念。向量具有大小和方向，滿足線性組合法則；標量是只有大小沒有方向的量（注：標量的另一種更深刻的定義是在坐標變換中保持不變的量）。向量模型的優點之一是其坐標系無關性， 也就是相對性，它在定義向量和運算規則的時候從一開始就拋開了坐標系的束縛，不管你坐標軸怎么旋轉，我都能適應，向量的線性組合、內積、叉積、線性變換等 等運算全部都是坐標系無關的。注意，所謂坐標系無關性不是說就沒有坐標系了，還是有的，剛才三角形例子的頂點就是用坐標表示的，只是在解決問題的時候不同 的坐標系不會構成影響。用一個比喻，Java號稱平臺無關，不是說Java就是空中樓閣，而是說你用Java編程時底層是Linux還是Windows往 往對你沒有影響。
向量模型有什么好處呢？除了剛才三角形面積問題是一個例子，下面我再舉一個幾何的例子：
給定三維坐標系中的一點(x0, y0, z0)和一個平面a*x + b*y + c*z + d = 0，求點到平面的垂直距離？
這個問題如果是要從解析幾何的角度去解決幾乎復雜到沒法下手，除非是平面恰好是過坐標軸的特殊情況，但是如果從向量模型考慮就很簡單：根據平面方程，平面的法向量(Normal Vector)是v=(a, b, c)，設從平面上任意一點(x, y, z)到(x0, y0, z0)的向量為w，那么通過內積dot_product(w, v)算出w到v的投影向量p，其大小就是(x0, y0, z0)到平面a*x + b*y + c*z + d = 0的垂直距離。這里用到了向量模型的基本概念：法向量，投影向量，內積，整個問題解決過程簡潔明快。
下面再給大家留一道相似的練習題（熟悉機器學習的朋友可能會發現這是線性代數在線性分類中的應用）：
給定n維空間中的兩點(a1, a2, ... an)，(b1, b2, ... bn)和一個超平面c1*x1 + c2*x2 ... + cn*xn + d = 0，請判斷兩點在超平面的同側或異側？
離開向量，下面我們要請出線性代數的另一個主角：矩陣(Matrix)。
線性代數定義了矩陣和向量、矩陣和矩陣的乘法，運算規則很復雜，用來做什么也不清楚，很多初學者都不能很好地理解，可以說矩陣是學好線性代數的攔路 虎。遇到復雜的東西，往往需要先避免一頭陷入細節，先從整體上把握它。其實，從程序的角度看，無論形式多么奇怪，它無非是一種語法，語法必然對應了語義， 所以理解矩陣的重點在于理解其語義。矩陣的語義不止一種，在不同的環境中有不同的語義，在同一環境中也可以有不同的解讀，最常見的包括：1)表示一個線性 變換；2)表示列向量或行向量的集合；3)表示子矩陣的集合。
矩陣作為一個整體對應的是線性變換語義：用矩陣A乘以一個向量v得到w，矩陣A就代表了v到w的線性變換。比如，如果想要把向量v0按逆時針方向旋轉60度得到v'，只需要用旋轉變換矩陣（Rotation Matrix)去乘v0就可以了。
除了旋轉變換，拉伸變換也是一種常見的變換，比如，我們可以通過一個拉伸矩陣把向量沿x軸拉伸2倍（請試著自己給出拉伸矩陣的形式）。更重要的是，矩陣乘法有一個很好的性質：滿足結合率。這就意味著可以對線性變換進行疊加，比如，我們可以把“沿逆時針旋轉60度”的矩陣M和“沿x軸拉伸2倍”的矩陣N相乘，得到一個新矩陣T來代表“沿逆時針旋轉60度并沿x軸拉伸2倍”。這是不是很像我們Shell中把多個命令通過管道進行疊加呢？
上面重點介紹了向量模型的坐標系無關性，除此之外，向量模型的另一優點是：線性性，因而它能用來表示線性關系，下面我們來看一個熟悉的Fibonacci數列的例子：
Fibonacci數列定義為：f(n) = f(n-1) + f(n-2), f(0) = 0, f(1) = 1；問題：輸入n，請給出求f(n)的時間復雜度不超過O(logn)的算法。
首先，我們構造兩個向量v1 = (f(n+1), f(n))和v2 = (f(n+2), f(n+1))，根據Fibonacci數列性質，我們可以得到從v1到v2的遞推變換矩陣：
并進一步得到：
這樣就把線性遞推問題轉化為了矩陣的n次冪經典問題，在O（log n）時間復雜度內解決。除了線性遞推數列，初等數學中著名的n元一次方程組問題也可以轉化為矩陣和向量乘法形式更容易地解決。這個例子是想說明，凡是滿足 線性關系的系統都是向量模型的用武之地，我們往往可以把它轉化為線性代數得到簡潔高效的解決方案。
總之，向量模型是整個線性代數的核心，向量的概念、性質、關系、變換是掌握和運用線性代數的重點，線性代數視為一門特定領域的程序語言，在初等數學基礎上建立了向量模型，定義了一套語法和語義，符合程序語言的語言契約。
END
∑編輯?|?Gemini
來源?|?今日頭條
微信公眾號“算法數學之美”，由算法與數學之美團隊打造的另一個公眾號，歡迎大家掃碼關注！
更多精彩：
?指數函數和自然對數
?科研這條路：一位數學博士給本科生的建議
?華裔天才數學家陶哲軒自述
?暢游一下分析數學
?主宰這個世界的10大算法
?瞎扯現代數學的基礎
?如何向5歲小孩解釋什么是支持向量機（SVM）？
?自然底數e的意義是什么？
?費馬大定理，集驚險與武俠于一體
?簡單的解釋，讓你秒懂“最優化” 問題
?一分鐘看懂一維空間到十維空間
??本科、碩士和博士到底有什么區別？
?小波變換通俗解釋
算法數學之美微信公眾號歡迎賜稿
稿件涉及數學、物理、算法、計算機、編程等相關領域，經采用我們將奉上稿酬。
投稿郵箱：math_alg@163.com
總結
- 上一篇： 众咖齐聚！2018智能汽车跨界融合高峰论
- 下一篇： 推荐几个牛逼技术号