📚 课程目录第2章 Jones矩阵入门

2. Jones矩阵入门：定义、常见偏振器件的矩阵表示

好，咱们正式开始啃Jones矩阵这块硬骨头。

说实话，我第一次接触Jones矩阵时，觉得这玩意儿不就是个2x2的矩阵嘛，有啥了不起的？后来真做项目才发现——嗯，真香。它把偏振光从“画箭头、看椭圆”的几何模式，直接拉进了“矩阵乘法、线性代数”的代数模式。你想想看，一个2x1的列向量描述光，一个2x2的矩阵描述器件，光通过器件就是矩阵乘向量。就这么简单。

2.1 为什么需要Jones矩阵？

在上一章我们聊了偏振态的几何表示——Jones向量。但光在光学系统中传播，会经过各种器件：线偏振片、波片、旋光器、反射镜……每个器件都会改变偏振态。

如果每个器件都用几何法去画椭圆、算参数，那系统一复杂，人就疯了。

Jones矩阵就是干这个的：它用一个2x2的复数矩阵，完整描述一个偏振器件对入射光Jones向量的线性变换。

核心思想：

入射光 Jones 向量 E_in，经过器件 Jones 矩阵 J，出射光 Jones 向量 E_out = J · E_in

就这么一个乘法，搞定一切。

说白了，Jones矩阵就是偏振光学里的“传递函数”。你只要知道器件的J，输入任意偏振光，输出一算便知。

2.2 Jones矩阵的数学定义

一个Jones矩阵 J 是2x2的复数矩阵：

J = [[Jxx, Jxy],
     [Jyx, Jyy]]

其中每个元素都是复数，表示x和y分量之间的耦合关系：

Jxx：x分量到x分量的增益/相位延迟
Jxy：y分量到x分量的耦合（交叉项）
Jyx：x分量到y分量的耦合（交叉项）
Jyy：y分量到y分量的增益/相位延迟

如果器件是非耦合型（比如理想线偏振片、波片），Jxy和Jyx为0，矩阵退化为对角阵。如果器件有旋光性或者扭曲，那非对角项就出来了。

我的经验：实际项目中，90%的常见器件都是对角阵或可通过对角化处理的。遇到非对角阵时，先别慌——看看能不能通过旋转坐标系变成对角阵。我在做液晶器件仿真时就经常用这招。

2.3 常见偏振器件的Jones矩阵

下面这几个矩阵，我建议你背下来。不是死记硬背，而是理解每个元素为什么是那个值。我在项目中反复用过它们，闭着眼都能写出来。

2.3.1 线偏振片

线偏振片只允许某一方向的电场分量通过。假设透光轴沿x轴方向：

J_pol_x = [[1, 0],
           [0, 0]]

如果透光轴沿y轴：

J_pol_y = [[0, 0],
           [0, 1]]

如果透光轴与x轴成θ角呢？这时候需要旋转矩阵帮忙：

J_pol(θ) = R(-θ) · J_pol_x · R(θ)

其中旋转矩阵：

R(θ) = [[cosθ, -sinθ],
        [sinθ,  cosθ]]

展开后得到：

J_pol(θ) = [[cos²θ,   sinθ·cosθ],
            [sinθ·cosθ, sin²θ   ]]

注意：线偏振片的Jones矩阵是奇异矩阵（行列式为0）。这意味着信息丢失了——你不可能从出射光反推出完整的入射光偏振态。我曾经在测量系统里踩过这个坑，以为加个偏振片能提高信噪比，结果把偏振信息给丢了。

2.3.2 波片（相位延迟器）

波片在两个正交方向之间引入固定的相位差。快轴沿x轴、慢轴沿y轴的波片：

J_waveplate = [[e^(i·δ/2), 0],
               [0,          e^(-i·δ/2)]]

其中δ是相位延迟量：

δ = π/2 → 四分之一波片（QWP）
δ = π → 二分之一波片（HWP）
δ = 2π → 全波片（基本没啥用，除了校准）

具体数值代入：

# 四分之一波片 (QWP)
J_QWP = [[e^(i·π/4), 0],
         [0,          e^(-i·π/4)]]

# 二分之一波片 (HWP)
J_HWP = [[i, 0],
         [0, -i]]   # 或者写成 [[e^(i·π/2), 0], [0, e^(-i·π/2)]]

关键理解：波片不改变光的强度（矩阵是酉矩阵），只改变两个分量之间的相位关系。说白了，它只转偏振态，不衰减光强。

2.3.3 常见器件对比表

器件	Jones矩阵（标准形式）	特性
线偏振片（x轴）	[[1,0],[0,0]]	奇异矩阵，强度衰减
线偏振片（y轴）	[[0,0],[0,1]]	奇异矩阵，强度衰减
四分之一波片（快轴x）	[[e^(iπ/4),0],[0,e^(-iπ/4)]]	酉矩阵，δ=π/2
二分之一波片（快轴x）	[[i,0],[0,-i]]	酉矩阵，δ=π
旋光器（θ角）	[[cosθ, -sinθ],[sinθ, cosθ]]	旋转偏振方向

2.4 知识结构图：Jones矩阵体系

下面这张图帮你理清Jones矩阵在整个偏振光学里的位置。我习惯把这类图贴在工位旁边，写代码时扫一眼就知道自己在哪一层。

2.5 系统级联：多个器件串联

实际光学系统很少只有一个器件。多个偏振器件串联时，总Jones矩阵等于各器件矩阵的逆序乘积：

J_total = J_n · J_(n-1) · ... · J_2 · J_1

注意顺序！光先经过J₁，再经过J₂，最后经过Jₙ。矩阵乘法不满足交换律，顺序错了结果全错。

我踩过的坑：有一次仿真一个偏振干涉仪，死活算不出正确结果。查了两天，发现是把波片和偏振片的顺序写反了。矩阵乘法里，J₂·J₁ ≠ J₁·J₂，这个教训我记到现在。

2.6 用Python验证Jones矩阵

光说不练假把式。咱们写几行Python，验证一下Jones矩阵的效果。这里用NumPy做矩阵运算，非常直观。

import numpy as np

# 定义Jones向量：45°线偏振光
E_in = np.array([[1], [1]]) / np.sqrt(2)

# 定义线偏振片（x轴透光）
J_pol_x = np.array([[1, 0],
                    [0, 0]])

# 定义四分之一波片（快轴x）
J_QWP = np.array([[np.exp(1j * np.pi/4), 0],
                  [0, np.exp(-1j * np.pi/4)]])

# 计算：光先经过QWP，再经过线偏振片
E_mid = J_QWP @ E_in
E_out = J_pol_x @ E_mid

print("入射光：", E_in)
print("经过QWP后：", E_mid)
print("最终出射光：", E_out)

# 验证：出射光应该是x方向线偏振
# 因为QWP把45°线偏转成圆偏，再被x偏振片滤出x分量

运行这段代码，你会看到E_out的y分量为0，只剩x分量。这就是Jones矩阵的威力——几行代码，替代了半页纸的几何推导。

小结一下：

Jones矩阵是2x2复数矩阵，描述偏振器件
常见器件：线偏振片（奇异阵）、波片（酉矩阵）、旋光器（旋转阵）
系统级联：逆序相乘，顺序不能错
Python+NumPy可以快速验证和仿真

掌握了这些，你已经可以处理大部分偏振光学问题了。下一块内容，咱们会深入Jones矩阵的本征分析——为什么有些偏振态通过器件后不变？那些就是器件的“本征偏振态”。不过那是后面的事了，先把今天的矩阵背熟再说。