每日一式 · No.015
✦ 第三季 · 倒帶鍵
x = A−1·b
反矩陣與聯立方程 · The Inverse & Linear Systems
施工可以倒帶嗎?—— 解聯立方程,就是問:「哪個輸入,造出了這個輸出?」
📜 公式的身世
被冠錯名的中國算法
「高斯消去法」這名字,其實冤枉了兩邊。《九章算術》第八章「方程」 —— 約成書於西元前後 —— 用「遍乘、直除」一步步消元解聯立方程,與今天的高斯消去法同一套骨架,早了一千八百年;為了消元時的減不過來,書裡還發明「正負術」,是人類最早系統性使用負數的紀錄。高斯本人則是拿這套老方法算小行星軌道與大地測量。名字掛誰,史官說了算;但只要有人開始解「幾個未知數綁在一起」的帳,就一定會再發明它一次。
🌍 它今天在哪裡上工
🍱
團膳的配料聯立 —— 蛋白質、熱量、預算三條限制,幾種食材各配多少?正是解聯立 —— 二戰美軍的「飲食問題」還催生了整門線性規劃。
⚡
電路板的體檢 —— 每個節點列一條基爾霍夫方程,幾百個未知電流一次解 —— 電路模擬軟體天天在做超大型消去法。
📊
諾貝爾獎的投入產出表 —— Leontief 把整個國家經濟寫成矩陣:鋼鐵要多少電、電要多少煤 —— 解反矩陣算政策連鎖效應,拿下 1973 年諾貝爾經濟獎。
012 教矩陣「往前施工」:輸入 x,輸出 Ax。今天反過來問:看到輸出 b,能不能倒推回輸入?能倒帶的變換,它的倒帶鍵就叫 A⁻¹(反矩陣);而國中就認識的「聯立方程」,真面目正是這個反問題。
🍱用你的生意講:加課那張表是「配方 → 成本」的正向試算;今天是稽核日 —— 帳上只剩「這批便當共花 b 元、含 g 克蛋白質」,反推當初用了幾份米、幾支雞。逆向工程 = A⁻¹;會計、偵探、工程師,天天都在解 Ax = b。
四章:先玩倒帶鍵本人(含 2×2 的口訣公式,det 就藏在分母 —— 伏筆);再揭聯立方程的兩種圖像(交點 vs 配方);然後直面 det=0 的倒帶失效世界(無解、或無限多解);最後說實話:電腦根本不按倒帶鍵 —— 高斯消去法,一步步走給你看。
Part 1 · 拉了就懂（白話）
I
第一章 · 倒帶鍵：A 做的，A⁻¹ 拆掉
旋轉的倒帶是反著轉、放大的倒帶是縮小、剪切的倒帶是反向剪 —— 合起來：A⁻¹A = I。
⏪每個變換是一道工序,反矩陣就是它的還原工序:轉 45° 的倒帶是轉 −45°;放大 1.5 倍的倒帶是縮回 1/1.5;往右剪 0.8 的倒帶是往左剪 0.8。連做「施工 → 倒帶」,小魚分毫不差回到原位 —— 數學上寫成一句:A⁻¹·A = I(什麼都沒做的單位矩陣)。
互動 · 播放「施工 → 倒帶」，看小魚完璧歸趙
玩玩看，注意這些
- 三顆晶片各播一次:施工把網格帶走,倒帶把它原路押回 —— 驗收欄永遠是 [1 0 / 0 1]。
- 看倒帶鍵長相:旋轉的 A⁻¹ 是「角度變號」;放大的是「倍率倒數」;剪切的是「剪量變號」—— 倒帶 = 把每個動作的參數取反。
- 順帶一提順序:先 B 後 A 的倒帶,是先拆 A、再拆 B:(AB)⁻¹ = B⁻¹A⁻¹ —— 先穿襪後穿鞋,脫的時候先脫鞋。
✦ 2×2 的倒帶口訣(一輩子夠用):A⁻¹ = (1/det)·[ d −b / −c a ] —— 對角互換、反對角變號、整體除以 det。注意分母:det 住在分母裡。它若是零…… 第三章見。
一句話帶走:A⁻¹ 是 A 的還原工序,A⁻¹A = I;2×2 口訣「換對角、反號、除 det」—— 而 det 在分母,是一顆定時伏筆。
II
第二章 · 聯立方程的兩種圖像：交點，或配方
2x+y=b₁、x+3y=b₂ —— 列圖像是兩條線的交點；欄圖像是把兩支箭頭調成 b 的配方。
🎯國中解聯立,老師只給過一種畫面:每條方程是一條線,解 = 交點(列圖像)。第三季要補上另一半,而且是更有力的一半:把方程組寫成 Ax = b,用加課的欄觀讀 —— b 是「x 份第一欄 + y 份第二欄」調出來的;解 (x, y),就是那份配方!拉動 b 的兩顆滑桿,同一題、兩種眼鏡,交替著看。
互動 · 2x + y = b₁ , x + 3y = b₂ （A = [2 1 / 1 3]）
怎麼看這張圖
- 列圖像:兩條線隨 b 平移,金點 = 交點 = 解。熟悉、直觀 —— 但到了 100 條方程、100 維,線就畫不動了。
- 切到欄圖像:金箭頭 = A 的第一欄 (2,1)、青箭頭 = 第二欄 (1,3),星星 = 目標 b。虛線接力:走 x 份金 + y 份青,不多不少落在星星上 —— 解就是配方。b=(5,10) 時配方 = (1, 3):一份金、三份青。
- 把 b 拉到讓 x 或 y 變負:配方箭頭反著走 —— 「負的份量」在幾何裡毫無懸念。欄圖像的好處:它在任何維度都成立,而且直通「b 到底調不調得出來」這個生死問題 —— 下一章。
一句話帶走:解聯立 = 找「哪個輸入」:列圖像看交點,欄圖像看配方 —— 解 (x,y) 就是把 A 的兩欄調成 b 的份量表。
Part 2 · 帶入深度
III
第三章 · det = 0：倒帶失效的世界
兩欄疊成一條線 —— 調不出線外的 b（無解），線上的 b 又有無限種調法（無限多解）。
🚫012 那句「det=0 = 壓扁、不可逆」,今天兌現成災難現場。換一顆病矩陣 A = [1 2 / 0.5 1]:它的兩欄 (1, 0.5) 與 (2, 1) 指同一個方向 —— 兩種原料,其實是同一味。於是能調出來的 b,只剩一條線(金色):b 落在線外 → 無解,怎麼調都到不了;b 恰好在線上 → 無限多解,因為兩味原料可以互相頂替。det = 1×1 − 2×0.5 = 0,倒帶鍵的分母,死了。
互動 · 病矩陣的世界：能到達的 b，只有那條金線
det = 1×1 − 2×0.50（倒帶鍵報廢）
狀態—
幾何—
玩玩看，注意這些
- 欄圖像:兩支箭頭疊在同一方向,世界塌成一條金線。拉 b 離開線 → 狀態轉紅「無解」;把 b₂ 慢慢調到 b₁ 的一半 → 星星落線,狀態轉綠「無限多解」,並亮出自由方向 (2, −1):沿它滑,配方隨你換、成品不變。
- 切列圖像:兩條方程線永遠平行(斜率同為 −1/2)—— 分開時無交點(無解);b₂ = b₁/2 時整條重合(每一點都是解)。兩種眼鏡,同一場災難。
- 災難的分類完成:det ≠ 0 → 唯一解;det = 0 → 無解或無限多解,再無第四種。這是 2×2 的「秩」的故事,高維同款。
⚠️ 工程師的加映(比 det=0 更常咬人):det 接近 0 的矩陣叫病態 —— 倒帶鍵還在,但抖:輸入 b 差 0.1%,解可以噴掉 100%。量測有雜訊時,病態比奇異更危險。這件事的正式度量叫條件數 —— 名字先掛著,下一課 SVD 給它正身。
一句話帶走:det=0 ⇒ 兩欄共線 ⇒ 可達的 b 塌成一條線:線外無解、線上無限解;而 det≈0 的病態矩陣,是工程日常裡更陰險的敵人。
IV
第四章 · 說實話：電腦不按倒帶鍵 —— 高斯消去法
實務上沒人先算 A⁻¹ 再乘 b。一步步消成三角形，再回頭代 —— 走給你看。
⚙️殘酷的行內真相:x = A⁻¹b 是理論的寫法,不是電腦的做法。先算整顆 A⁻¹ 要多花約三倍力氣,還會放大浮點誤差 —— 所以 numpy 的 solve() 和 inv() 是兩個函數,前者永遠是對的選擇。電腦實際做的是 1810 年代的老手藝:高斯消去 —— 用列運算把矩陣消成三角形(加課列觀:每一步都是「拿一列去修另一列」),再從最後一行回代上來。按「下一步」,親手走完一題 3×3。
互動 · 解 2x+y+z=7 , x+3y+2z=13 , x+z=4
進度0 / 6
這一步起點：增廣矩陣 [A | b]
怎麼看這張圖
- 前三步是消去:金色是主元列(工具),青色是被修的列 —— 目標是把左下角全消成 0,矩陣變上三角。
- 後三步是回代:最後一行只剩一個未知數,z 先落地;帶著 z 往上一行解 y;再帶著 (y,z) 解 x —— (1, 2, 3),收工。
- 成本帳:消去法約 n³/3 次乘法;先算 A⁻¹ 再乘要 ~n³ —— 貴三倍,還更不準。教科書寫 A⁻¹b 是為了「講理」,程式寫 solve(A, b) 是為了「做事」。
📦 LU 分解(工廠版消去):把消去的動作打包存成 A = L·U(下三角×上三角)—— 之後同一顆 A、換一百個 b,每次只要便宜的回代。排程、模擬、電路求解天天重用同一結構,靠的就是這招。
🏗️結構與電路
桁架受力、KCL 電流方程 —— 工程分析的心臟就是一場大型 Ax=b;有限元素法把橋樑切成百萬條方程。
📈最小平方法
方程比未知數多(量測過剩)時無精確解 —— 改解 AᵀAx = Aᵀb 求「最接近」。下一課 SVD 會給它最漂亮的答案。
🧮克拉瑪公式
用行列式直接寫解 —— 理論優雅(2×2 口訣的親戚),但 n 大時慢到天荒地老;考試用,實戰不用。
🤖科學計算的地基
天氣預報、晶片模擬、AI 訓練的每一步底層,都踩在「快速穩定地解 Ax=b」上 —— LAPACK 是現代文明的隱形樑柱。
一句話帶走:消去成三角、回代得解 —— 比倒帶鍵便宜三倍又穩;LU 讓同一顆 A 重用百次。理論寫 A⁻¹b,實戰寫 solve(A,b)。
把整堂課收成一句話 ——
解聯立方程 = 問「哪個輸入,造出了這個輸出」。
倒帶鍵 A⁻¹:換對角、反號、除 det;
兩種圖像:交點,或配方;
det=0:線外無解、線上無限;
而電腦,用兩百年前的消去法,又快又穩。
唯一沒說的是:無解的時候,難道就認了嗎?
—— 下一課,SVD 連這個都有答案。