3、琼斯矩阵与斯托克斯参数：极化状态的数学表征方法
各位同学，大家好。今天我们来聊聊极化状态的数学工具。说实话，刚入行那会儿，我觉得极化就是个“电场怎么抖”的问题，没必要搞那么复杂。直到我第一次设计超表面极化转换器，仿真结果和实测对不上，折腾了两周才发现是极化表征方式没选对。嗯，从那以后，我再也不敢小看这些数学工具了。
3.1 为什么需要两种数学工具？
你想想看，我们描述一个电磁波的极化，本质上是在描述电场矢量在空间中的轨迹。这个轨迹可以是直线、圆，或者椭圆。但问题来了——我们怎么把这个轨迹“翻译”成数学语言，让计算机能算、让工程师能设计？
这里有两个主流方案：
- 琼斯矩阵：用复振幅描述极化，适合处理相干光、单色波
- 斯托克斯参数：用强度描述极化，适合处理非相干光、部分极化波
我个人习惯把琼斯矩阵比作“精确制导导弹”，斯托克斯参数则是“地毯式轰炸”。前者精度高但娇气，后者鲁棒性强但信息量略少。在超材料设计中，我通常先用琼斯矩阵做精确优化，再用斯托克斯参数验证鲁棒性。
3.2 琼斯矩阵：极化转换的“传递函数”
琼斯矩阵的核心思想很简单：把入射波的极化状态看作一个二维复向量，经过超表面后，出射波的极化状态就是另一个二维复向量。这个映射关系用一个2×2的复矩阵来描述，就是琼斯矩阵。
数学上，我们这样写：
E_out = J · E_in
其中：
E_out = [E_x_out, E_y_out]^T
E_in = [E_x_in, E_y_in]^T
J = [[J_xx, J_xy],
[J_yx, J_yy]]
这里每个J分量都是复数，包含了幅度和相位信息。我在项目中遇到过最典型的场景是设计线极化转圆极化的超表面。这时候，琼斯矩阵的理想形式是：
J_ideal = [[1, 0],
[0, i]]
这个矩阵意味着：x极化波通过后幅度不变，y极化波通过后相位滞后90°。两者叠加，就形成了圆极化。
关键点：琼斯矩阵的四个元素不是独立的。对于无源、无增益的超表面，它们必须满足能量守恒和互易性约束。具体来说，|J_xx|² + |J_yx|² ≤ 1，且J_xy = -J_yx（在特定坐标系下）。
3.3 斯托克斯参数：从强度看极化
琼斯矩阵虽然精确，但它有个致命弱点——它只适用于完全极化波。现实世界中，大部分电磁波都是部分极化的，比如太阳光、LED灯光。这时候，斯托克斯参数就派上用场了。
斯托克斯参数由四个实数组成：
S0 = |E_x|² + |E_y|² (总强度)
S1 = |E_x|² - |E_y|² (水平/垂直极化差)
S2 = 2·Re(E_x·E_y*) (45°/135°极化差)
S3 = 2·Im(E_x·E_y*) (左旋/右旋极化差)
这四个参数可以直观地映射到庞加莱球上。S1、S2、S3就是球面上的三维坐标，S0是球的半径。完全极化波落在球面上，部分极化波落在球内部。
我记得有一次调试一个宽带极化转换器，用琼斯矩阵算出来的性能很漂亮，但实测总是差一截。后来改用斯托克斯参数分析，才发现是部分极化分量在作怪——超表面在某些频点产生了退极化效应。这个坑，我踩过。
实用技巧：在超材料设计中，我建议同时计算琼斯矩阵和斯托克斯参数。琼斯矩阵帮你理解“理想情况”，斯托克斯参数帮你评估“实际表现”。两者结合，才能做出靠谱的设计。
3.4 两种工具的转换关系
既然两种工具各有千秋，那它们之间怎么转换？其实很简单：
| 从琼斯向量到斯托克斯参数 | 从斯托克斯参数到琼斯向量 |
|---|---|
| S0 = |E_x|² + |E_y|² | |E_x| = sqrt((S0 + S1)/2) |
| S1 = |E_x|² - |E_y|² | |E_y| = sqrt((S0 - S1)/2) |
| S2 = 2·Re(E_x·E_y*) | Δφ = arctan(S3 / S2) |
| S3 = 2·Im(E_x·E_y*) | 注意：仅适用于完全极化波 |
这里有个坑：从斯托克斯参数反推琼斯向量时，会丢失绝对相位信息。因为斯托克斯参数只关心相对相位差，不关心绝对相位。所以，如果你需要知道波前的绝对相位（比如做干涉设计），还是得用琼斯矩阵。
3.5 知识体系总览
下面这张图是我自己总结的极化表征知识框架，帮你理清思路：
3.6 实战中的选择建议
说了这么多，到底什么时候用琼斯矩阵，什么时候用斯托克斯参数？我根据自己的经验，给你几个建议：
- 设计阶段：用琼斯矩阵。因为你需要精确控制每个单元的幅度和相位响应，琼斯矩阵的复数值正好满足需求。
- 测试阶段：用斯托克斯参数。因为实际测量中，你很难直接测到复振幅，但强度测量很容易。斯托克斯参数可以直接从强度测量结果中提取。
- 宽带设计：两者都要用。琼斯矩阵帮你优化中心频点，斯托克斯参数帮你评估带宽内的稳定性。
避坑指南：我曾经在一个项目中，直接用琼斯矩阵优化了一个超表面，仿真结果完美。但加工回来后，测试结果一塌糊涂。后来发现，问题出在加工误差导致的退极化效应上。如果当时先用斯托克斯参数评估一下鲁棒性，就不会走这个弯路了。
好了，关于极化状态的数学表征，今天就聊到这里。记住一句话：琼斯矩阵是“显微镜”，看得细；斯托克斯参数是“广角镜”，看得全。两者结合，才能把极化转换器设计得既精准又可靠。
公众号:蓝海资料掘金营，微信deep3321