跳轉到

5.2 SPECTRAL理論的逼近觀點#

本章導言#

Fredholm 用行列式、Hilbert 用二次型、Riesz 用緊性——三種語言背後有一個共同的技術事實：連續核可以被退化核一致逼近（1.2 節），於是無窮維問題成為「有限維問題列的極限」。本節把這個逼近觀點精緻化：緊算子=有限秩算子的範數極限（Hilbert 空間），譜的收斂性（ \(\sigma(T_n) \to \sigma(T)\) ？），以及 Galerkin 法的數學骨架。逼近觀點也是數值分析的靈魂：實際計算永遠是有限維的，譜理論必須回答「有限維近似何時收斂到無窮維真相」。

有限秩逼近#

定理 5.2.1： \(H\) 有正交基時， \(T \in \mathcal{B}(H)\) 緊 \(\iff\) \(T\) 是有限秩算子的一致極限。

證明：設 \(\{e_n\}\) 正交基， \(P_n\) = 投影到 \(\operatorname{span}\{e_1, \ldots, e_n\}\) 。 \(T\) 緊時 \(T P_n \to T\) （一致收斂）：否則存在 \(\epsilon > 0\) 與單位 \(x_n \perp e_1, \ldots, e_n\) （即 \(x_n \in \operatorname{ran}(I - P_n)^\perp\) 意義下的單位向量）使 \(\|Tx_n - TP_nx_n\| = \|T(I - P_n)x_n\| \ge \epsilon\) ；但 \((I - P_n)x_n \rightharpoonup 0\) 弱收斂，與緊性（弱收斂列的像強收斂——Hilbert 空間中緊=完全連續）矛盾。 \(TP_n\) 有限秩。∎

這個證明展示了緊性的第二張面孔：緊=壓制「弱逃逸」。無窮維中向量可以「往無窮遠座標跑」（弱收斂到 0），緊算子把這種逃逸壓成強收斂到 0——所以截斷逼近有效。非緊算子（如 \(I\) ）壓不住： \(e_n \rightharpoonup 0\) 但 \(\|Ie_n\| = 1\) 。

譜的收斂性#

逼近觀點的核心問題： \(T_n \to T\) （一致）時， \(\sigma(T_n) \to \sigma(T)\) ？答案是上半連續而非連續：

定理 5.2.2：(1) 若 \(\lambda \in \rho(T)\) 、 \(T_n \to T\) ，則對足夠大的 \(n\) ， \(\lambda \in \rho(T_n)\) 且 \(R^{(n)}_\lambda \to R_\lambda\) （Neumann 級數控制 \(\|(T_n - T)R_\lambda\| < 1\) ）。(2) 但 \(\lambda_n \in \sigma(T_n)\) 可以收斂到 \(\rho(T)\) 中的點——譜可以「消失」。

例（譜消失）： \(H = l^2\) ， \(T_n =\) 左位移 \(L\) 的截斷 \(+\) 一點微調？經典例子：取 \(T_n = J_{n}\) （ \(n \times n\) 右上角 Jordan 塊嵌入 \(l^2\) ，其餘為 \(0\) ）。 \(\sigma(T_n) = \{0\}\) （冪零），但 \(T_n \to S^*\) （左位移）弱？不——換更好的例子： \(T_n = \frac{1}{n}\) (位移)型。簡潔的例子：對角列 \(T_n = \operatorname{diag}(1, \frac12, \ldots, \frac1n, 0, 0, \ldots)\) 收斂（一致）到 \(D = \operatorname{diag}(\frac1k)_k\) 。這裡譜 \(\{1, 1/2, \ldots, 1/n\} \cup \{0\} \to \{0\} \cup \{1/k\}\) ——譜穩定收斂。反例：有限秩投影 \(P_n\) 強收斂到 \(I\) ， \(\sigma(P_n) = \{0,1\}\) 但 \(\sigma(I) = \{1\}\) ——點 \(0\) 消失。一致收斂下 \(0\) 點也會消失（ \(T_n = \frac{1}{n}I\) ： \(\sigma = \{1/n\} \to \emptyset\) ？不， \(\to \{0\}\) = \(\sigma(0)\) ，穩定）。真正的不穩定性出現在非正規算子與強收斂情形；一致收斂+正規算子時譜依 Hausdorff 距離連續（Weyl–Hausdorff 定理）。∎

定理（Weyl–Hausdorff）：正規算子 \(T_n \to T\) 一致收斂，則 \(\sigma(T_n) \to \sigma(T)\) 按 Hausdorff 度量。正規性再次是可計算性的保障。

例題#

例題 5.2.1（Galerkin 法）：用 \(\{e_k\}_{k=1}^n\) 逼近積分方程 \(\varphi = \lambda K\varphi + f\) （ \(K\) 緊自伴、對稱連續核）。設 \(K\) 的特徵對 \((\mu_j, e_j)\) （6.2 節）， \(f = \sum f_j e_j\) 。求 Galerkin 解並證明收斂。

解：Galerkin：求 \(\varphi^{(n)} = \sum_{j\le n} c_j e_j\) 使殘差與子空間正交。展開後方程精確可解：

\[ c_j = \frac{\lambda \mu_j c_j + f_j}{1} \implies c_j = \frac{f_j}{1 - \lambda\mu_j} \quad (\lambda \mu_j \neq 1), \]

即 \(\varphi = \sum_j \frac{f_j}{1 - \lambda\mu_j} e_j\) 的部分和——在特徵基下，積分方程是無窮多個標量方程（這就是 Hilbert 1906 年的洞察！）。收斂性： \(\mu_j \to 0\) （緊性），故 \(1 - \lambda\mu_j \ge 1 - |\lambda|\mu_1 > 0\) （若 \(\lambda \mu_1 \neq 1\) ），係數 \(|c_j| \le \frac{|f_j|}{1 - |\lambda|\mu_1}\) ，而 \(\sum |f_j|^2 < \infty\) （ \(f \in L^2\) ），故級數收斂。∎

此例是 5.3 節擇一定理的譜式證明： \(\lambda\mu_j = 1\) （某個特徵值倒數）時，需要 \(f_j = 0\) （正交條件）——Fredholm 擇一在譜語言中是「分母為零」。

例題 5.2.2：證明緊自伴算子的特徵值可以由 Rayleigh 商的極小極大過程逐步逼近（Courant–Fischer 的無窮維版）： \(\mu_k = \min_{\dim S = k-1}\max_{x \perp S, \|x\| = 1}|\langle Tx, x\rangle|\) 。

解：用 1.3 節 Hilbert 極大化迭代+歸納： \(\mu_1 = \|T\|\) 可達（3.2）；在 \(V_1^\perp\) （ \(V_1\) =第一特徵空間）上重複，得 \(\mu_2 \le \mu_1\) ； \(\mu_k \to 0\) 由緊性（若 \(\mu_k \ge \epsilon\) ，特徵向量正交列的像無收斂子列）。極小極大公式由「特徵空間內的 max + 正交補上的 min」的張力論證。∎

歷史啟示#

  • 逼近觀點是 Hilbert 學派的本能：Hilbert 的整個積分方程理論就是「無窮多變數二次型=有限二次型的極限」。1950 年代 Courant 學派把它變成數值方法（Galerkin 1915 已提出、Ritz 1909 變分版），1990 年代變成有限元工業。
  • 譜收斂的不對稱性：預解集穩定（開+Neumann 控制），譜只能「冒出來」不能「莫名消失」（一致收斂下）；正規性給 Hausdorff 連續。數值線性代數的誤差理論建立在這些不對稱性上。
  • 弱收斂與緊性的配合：有限秩逼近的證明用「弱逃逸被緊性壓制」——第 2 章的弱拓撲（2.3）在第 5 章發力，這是抽象理論回報具體計算的時刻。

與其他章節及領域的關聯#

  • 5.1 節：緊性=有限秩可逼近（本節精確化）。
  • 5.3 節：擇一定理的譜式證明（例題 5.2.1 的推廣）。
  • 6.2 節：Hilbert–Schmidt 定理把例題 5.2.1 的「特徵基下的標量方程」變成一般定理。
  • 7.3 節：Friedrichs 擴張=「二次型逼近」思想在無界算子的應用。
  • 第 9 章：Galerkin 半離散化是發展方程數值解的標準方法。
  • 數值分析：譜收斂理論是特徵值算法（Lanczos 1950、Arnoldi 1951）的誤差分析基礎；偽譜（4.1 場景二）處理非正規情況的病態。

練習#

  1. 證明 Hilbert 空間中「緊⇒弱收斂列的像強收斂」與「有限秩可逼近」等價。
  2. 對例題 5.2.1，取 \(\lambda = 1/\mu_1\) （共振），描述解的集合與正交條件的具體形式。
  3. 構造正規算子列 \(T_n \to T\) （一致）使 \(\sigma(T_n) \not\to \sigma(T)\) ？若做不到，證明 Weyl–Hausdorff 定理的特例（對角算子情形）。
  4. 用 Galerkin 法在 \(\operatorname{span}\{1, x, x^2\}\) 上逼近例題 1.1.1 的 Volterra 方程，比較數值與精確解 \(e^x\) 。
  5. 思考題：為什麼 Lanczos 算法對對稱矩陣快而對非對稱矩陣需要 Arnoldi？從「正規性=譜 Hausdorff 連續」出發討論。

進階研讀：三個深化場景#

場景一：Ritz 法與變分原理。Ritz（1909）用極小化 Rayleigh 商在有限維子空間求特徵值近似；誤差理論（Szegő 1915：Ritz 值從上方逼近特徵值）連接到變分法。Ritz 值與 Lanczos 值在數值分析中統一為 Krylov 子空間方法——1.3 節的極大化迭代百年後成為最快的特徵值算法。

場景二：譜逼近的抽象理論。Haraux（1979）、Kato 的譜逼近理論、以及 Mosco 收斂（1969：二次型的收斂拓撲）處理「變分極限中的譜」。物理應用：週期結構的均勻化（homogenization）——微觀係數極快振盪時，譜的極限行為由 Mosco 極限的二次型決定。

場景三：非正規逼近的陷阱。Greenbaum–Trefethen（1990s）展示：非正規算子的 Ritz 值可以完全誤導（偽譜遠大於譜），「計算出的特徵值」可能只是瞬態增益的影子。這個警告重塑了數值分析的哲學：對非正規系統，譜不是穩定性的完整故事——回頭看 4.3 節 Jordan 塊的 \(\|T^n\| = n\) ，有限維已埋下伏筆。

場景四：譜逼近在偏微分方程特徵值問題中的實作。Laplace 特徵值 \(\lambda_n \to \infty\) （非緊方向！）： \(\Delta\) 無界，但其逆緊——Galerkin 對 \(\Delta\) 本身用「弱形式+協調逼近」，收斂理論（Babuška–Osborn 1991）是 5.2 節譜逼近論的雙曲版本：Ritz 值從上方逼近（因 \(\Delta\) 下有界、極小極大成立）。誤差 \(\lambda_n - \lambda_n^{(h)} \sim h^{2k}\) （ \(k\) =多項式階）——譜收斂速率理論使特徵值計算（模態分析、聲學設計）有嚴格誤差界。「無界算子的譜逼近=其逆的緊譜逼近+譜映射倒轉」：理論分工的教科書案例。

本節方法論清單#

  • 對角化優先：先找特徵基（或譜投影），方程退化為標量級數——一切截斷皆為係數截斷。
  • 截斷即投影：Galerkin=正交投影 \(P_n\) 作用於方程兩側，殘差與子空間正交。
  • 一致有界是收斂的通行證： \(\|(I - \lambda K_n)^{-1}\|\) 的一致控制保證近似解的穩定性（Quasi-optimal）。
  • 譜的單向漂移：Ritz 值從一側逼近真譜（自伴時從上、或依分類），誤差有符號——變分結構的紅利。
  • 弱逃逸的壓制：有限秩逼近成立靠「弱收斂列的像強收斂」——緊性是逼近論的燃料。
  • 非正規的陷阱：Ritz 值可能只是偽譜的影子，殘差檢驗與偽譜圖是必要的驗屍程序。
  • 從解析到數值：Hilbert 1906 的「有限維逼近列」是今日有限元/譜方法/Lanczos 的共同祖先。

本節譜收斂速查#

  • 預解集穩定： \(\lambda \in \rho(T)\) 時 \(T_n \to T\) 給 \(\lambda \in \rho(T_n)\) （ \(n\) 大）且 \(R^{(n)}_\lambda \to R_\lambda\) 。
  • 譜只冒不消失（一致收斂）： \(\lambda_n \in \sigma(T_n)\) 、 \(\lambda_n \to \lambda\) ⇒ \(\lambda \in \sigma(T)\) 。
  • 正規性給 Hausdorff 連續： \(\operatorname{dist}_H(\sigma(T_n), \sigma(T)) \to 0\) 。
  • 強收斂可丟譜： \(P_n \to I\) 強收斂而 \(0 \in \sigma(P_n) \setminus \sigma(I)\) 。
  • 緊算子的截斷： \(\sigma(T_n) \cup \{0\} \to \sigma(T)\) 的 Hausdorff 收斂（ \(T_n = TP_n\) ）。
  • Ritz 值的單側性：自伴情形極小極大給出誤差符號——變分結構的紅利。

名詞速查表#

  • 有限秩逼近：緊算子=有限秩算子的一致極限。
  • 弱逃逸： \(x_n \rightharpoonup 0\) 而 \(\|x_n\| = 1\) ；緊算子壓制之。
  • Weyl–Hausdorff 定理：正規算子譜的一致收斂連續性。
  • Galerkin 法：向有限維子空間的投影近似；特徵基下成為標量方程組。
  • Rayleigh 商： \(R(x) = \langle Tx,x\rangle/\|x\|^2\) ；極小極大刻畫特徵值。
  • Ritz 值：Rayleigh 商在 Krylov 子空間中的極值。

年表速查#

  • 1906：Hilbert 的「有限維逼近無窮維」方法論。
  • 1909：Ritz 變分法；Galerkin（1915）法。
  • 1918：Riesz 緊算子公理化（逼近觀點的抽象化）。
  • 1915–1935：Szegő、Courant 誤差理論。
  • 1950–1951：Lanczos、Arnoldi 迭代特徵值算法。
  • 1969：Mosco 收斂。
  • 1990s：Greenbaum–Trefethen 偽譜警告。

深化問答三則#

問：為什麼「有限秩可逼近」在 Hilbert 空間免費而在一般 Banach 空間是深刻問題？ 答：Hilbert 空間有正交投影族 \(P_n \to I\) （強），且 \(TP_n \to T\) 由緊性保證。一般 Banach 空間可能沒有任何有限秩算子列強收斂到 \(I\) （無逼近性質 AP，Enflo 1973 反例）。這再次顯示內積結構的特權：正交投影是「最好」的有限維近似器，其存在依賴內積。

問：Galerkin 法為什麼對特徵基下的積分方程「精確可解」？ 答：因為特徵基同時對角化方程與投影：方程在特徵基下已是標量方程組，Galerkin 投影只是截斷級數。「對角化→標量化→截斷」是譜方法的一般套路；其成敗繫於特徵基是否已知（解析情形）或可高效近似（Lanczos 情形）。

問：譜的「冒出」與「消失」哪個對數值計算更危險？ 答：消失。冒出的偽特徵值會被殘差檢驗抓住（ \(\|Tx - \lambda x\|\) 大）；消失的真特徵值是無聲失蹤——尤其非正規算子的強收斂逼近（如離散化微分算子）可以丟掉整段連續譜（Krein–Matsaev 效應），數學物理的數值工作中必須用偽譜或解析延拓驗證。

（本節完）