透鏡成像、矩陣光學與像差

從成像公式到光線的線性代數

一束來自物點的光，經過透鏡後重新聚於一點——這是「成像」最樸素的直覺。但要把整個光學系統（多片透鏡、不同曲率、不同間距）當成一個可計算、可串接的整體來處理，僅靠初等的成像公式 $\frac{1}{v}-\frac{1}{u}=\frac{1}{f}$ 是不夠的。本篇從近軸近似出發，建立矩陣光學（ABCD 矩陣）的形式體系，再說明為何「完美成像」在真實系統中必然破裂，導入賽德爾（Seidel）像差的物理機制。

在近軸（paraxial）區域，我們只保留光線角度與高度的一階項，光的傳播因此成為一個線性映射。一條光線在任一參考面上由二維狀態向量描述：高度 $y$ 與斜率 $\theta$（嚴格而言用光學方向餘弦 $n\theta$，$n$ 為折射率）。任何光學元件都對應一個 $2\times 2$ 矩陣，作用在 $(y,\ n\theta)^{\mathsf T}$ 上。

近軸近似與折射的線性化

考慮一條光線以小角度 $\theta$ 入射到曲率半徑為 $R$ 的球面，兩側折射率為 $n_1$、$n_2$。Snell 定律 $n_1\sin\theta_1 = n_2\sin\theta_2$ 在近軸下展開 $\sin\theta\approx\theta$，得 $n_1\theta_1 = n_2\theta_2$。利用球面法線與入射高度的幾何關係，可推得單一折射面對「縮減角度」$u\equiv n\theta$ 的變換：

$$ n_2\theta_2 = n_1\theta_1 - \frac{n_2-n_1}{R}\,y, \qquad y_2 = y_1. $$

定義該面的屈光力（power） $P=\dfrac{n_2-n_1}{R}$，折射矩陣為

$$ \mathcal{R}=\begin{pmatrix} 1 & 0 \\ -P & 1 \end{pmatrix}. $$

自由空間中行進距離 $d$（折射率 $n$）時，高度隨斜率累積、角度不變：

$$ \mathcal{T}=\begin{pmatrix} 1 & d/n \\ 0 & 1 \end{pmatrix}. $$

一個系統的整體矩陣即為各元件矩陣由右至左相乘。對薄透鏡（兩折射面夾零距離，置於空氣 $n=1$），把兩個 $\mathcal R$ 相乘並用磨鏡者公式 $\frac{1}{f}=(n-1)\!\left(\frac{1}{R_1}-\frac{1}{R_2}\right)$，得到簡潔的薄透鏡矩陣：

$$ M=\begin{pmatrix} A & B \\ C & D \end{pmatrix}=\begin{pmatrix} 1 & 0 \\ -1/f & 1 \end{pmatrix}. $$

一個關鍵性質：所有這些矩陣的行列式 $AD-BC=n_1/n_2$，在前後同介質時 $\det M=1$。這正是相空間中Liouville 定理（面積守恆，即光學的 Lagrange 不變量守恆）的代數表現，後文會再回到它的力學意涵。

成像條件與系統矩陣

把物面到像面的整段系統寫成 $M_{\text{tot}}=\mathcal{T}_2\,M_{\text{lens}}\,\mathcal{T}_1$。成像的定義是：像點高度只依賴物點高度、與出射角度無關——也就是說，從同一物點發出、角度各異的光線都匯聚於同一像高。這要求總矩陣的 $B$ 元素為零：

$$ B_{\text{tot}}=0. $$

此時 $A_{\text{tot}}$ 就是橫向放大率 $m$，而 $C_{\text{tot}}$ 與整體屈光力相關。對單一薄透鏡代入 $\mathcal T(v)\,M\,\mathcal T(u)$ 展開 $B=0$，立刻還原出高斯成像公式 $\frac1v-\frac1u=\frac1f$。矩陣法的價值在於：任意複雜系統（多片組、厚透鏡、空氣間隔）都用同一規則機械化處理，且能自然定義主平面、節點、焦距等基本參數——焦距由 $f=-1/C$ 讀出，主平面位置由 $A$、$D$ 與焦距決定。

定量小範例

兩片薄透鏡，$f_1=+10\,\text{cm}$、$f_2=+15\,\text{cm}$，間距 $d=5\,\text{cm}$（空氣中）。系統矩陣：

$$ M=\begin{pmatrix} 1 & 0 \\ -1/15 & 1 \end{pmatrix}\!\begin{pmatrix} 1 & 5 \\ 0 & 1 \end{pmatrix}\!\begin{pmatrix} 1 & 0 \\ -1/10 & 1 \end{pmatrix}. $$

逐步計算。先算右兩矩陣（傳播後接第一透鏡，由右至左）：

$$ \begin{pmatrix} 1 & 5 \\ 0 & 1 \end{pmatrix}\!\begin{pmatrix} 1 & 0 \\ -0.1 & 1 \end{pmatrix}=\begin{pmatrix} 0.5 & 5 \\ -0.1 & 1 \end{pmatrix}. $$

再左乘第二透鏡：

$$ \begin{pmatrix} 1 & 0 \\ -1/15 & 1 \end{pmatrix}\!\begin{pmatrix} 0.5 & 5 \\ -0.1 & 1 \end{pmatrix}=\begin{pmatrix} 0.5 & 5 \\ -0.1\overline{3} & 0.6\overline{6} \end{pmatrix}. $$

其中 $C=-\tfrac{1}{15}(0.5)-0.1=-0.1333\,\text{cm}^{-1}$。等效焦距 $f_{\text{eff}}=-1/C=7.5\,\text{cm}$。驗證：合成屈光力 $\frac1f=\frac1{10}+\frac1{15}-\frac{5}{10\cdot15}=0.1333\,\text{cm}^{-1}$，完全一致。

為何完美成像會破裂：賽德爾像差

矩陣光學的所有美好都建立在 $\sin\theta\approx\theta$ 的一階截斷上。把展開保留到三階，$\sin\theta=\theta-\frac{\theta^3}{6}+\cdots$，三次項就破壞了線性映射，使「同物點不同光線」不再嚴格共點。這就是幾何像差的根源。

以波前像差（實際波前與理想球面波的光程差）$W(\rho,h,\phi)$ 展開——$\rho$ 為光瞳半徑、$h$ 為物高、$\phi$ 為方位角——最低階（四次）共有五個獨立係數，即五種賽德爾單色像差：

$$ W = \underbrace{a\,\rho^4}_{\text{球差}} + \underbrace{b\,\rho^3 h\cos\phi}_{\text{彗差}} + \underbrace{c\,\rho^2 h^2\cos^2\phi}_{\text{像散}} + \underbrace{d\,\rho^2 h^2}_{\text{場曲}} + \underbrace{e\,\rho\,h^3\cos\phi}_{\text{畸變}}. $$

橫向光線像差正比於波前對光瞳座標的梯度 $\nabla_\rho W$，因此球差使邊緣光線聚焦較近、產生隨 $\rho^3$ 增長的模糊斑；彗差讓離軸點呈彗星狀拖尾；像散使子午與弧矢焦線分離。此外，介質折射率隨波長變化 $n(\lambda)$，導致色像差——焦距 $f$ 因 $\lambda$ 而異，這是與單色賽德爾像差獨立的另一族缺陷，需以不同色散材料（冕牌＋火石玻璃）的消色差雙合透鏡（achromatic doublet）補償。透鏡設計的核心工作，正是在多自由度（曲率、間距、玻璃）下使這些係數彼此抵消。

深入探討（研究所視角）

矩陣光學的 $\det=1$ 並非偶然，而是幾何光學具有辛結構（symplectic structure）的影子。把光線傳播類比於力學：費馬原理 $\delta\!\int n\,\mathrm ds=0$ 與哈密頓原理 $\delta\!\int L\,\mathrm dt=0$ 在形式上同構，光程扮演作用量的角色。選定光軸 $z$ 為「時間」，橫向座標 $\mathbf q=(x,y)$ 為廣義座標，其共軛動量為 $\mathbf p=n\,\mathrm d\mathbf q/\mathrm ds$（縮減方向餘弦）。光學拉格朗日量

$$ L(\mathbf q,\dot{\mathbf q},z)=n(\mathbf q,z)\sqrt{1+\dot x^2+\dot y^2} $$

經 Legendre 變換得光學哈密頓量 $H=-\sqrt{n^2-p_x^2-p_y^2}$，光線軌跡即 $\mathrm d\mathbf q/\mathrm dz=\partial H/\partial \mathbf p$、$\mathrm d\mathbf p/\mathrm dz=-\partial H/\partial \mathbf q$。近軸極限下 $H$ 展開到二次，運動方程線性化，其解算子恰為 $Sp(2)$ 群中的 $2\times 2$ 辛矩陣——這就是 ABCD 矩陣的群論身分。賽德爾像差則對應 $H$ 的高階項所生成的非線性（保辛但非線性）流，即正則變換的高階生成函數，可用 Lie 變換（Lie algebraic optics，Dragt 等人發展）系統展開到任意階。Lagrange 不變量 $\mathcal H=n(y\bar\theta-\bar y\theta)$ 守恆，正是辛形式 $\mathrm dp\wedge\mathrm dq$ 的不變性。

從更基礎的層面看，幾何光學是波動光學的短波長極限。將純量波 $u=A\,e^{ik_0 S}$ 代入亥姆霍茲方程，令 $k_0\to\infty$，最低階得程函方程（eikonal equation） $|\nabla S|^2=n^2$，其特徵線即光線；次階給出強度沿光線的輸運方程（能流守恆）。這套 WKB 式展開與量子力學中 $\hbar\to 0$ 的半古典近似在數學上完全平行：程函 $S$ 對應作用量、光線對應古典軌道、波前對應等作用量面。因此「成像」在波動語言下是 Fourier 光學問題——理想成像系統是空間頻率的低通濾波器，其點擴散函數（PSF）由光瞳函數的繞射決定，而像差正是光瞳函數中附加的相位因子 $\exp(ik_0 W)$。

近代發展把這套框架推得更遠：自適應光學以變形鏡即時補償波前像差 $W$，源於天文與視光學需求；自由曲面（freeform）與超穎透鏡（metalens）以次波長結構直接設計相位分布，繞過傳統球面在賽德爾係數上的耦合限制；而電子顯微鏡中的球差校正器（aberration corrector）以多極場抵消電子光學的 $C_s$，使解析度突破原子尺度。這些前沿都共享同一核心：把成像理解為相空間中的辛映射加上其上的相位修正，並在波動與幾何兩套語言間自由切換。