3、变换光学基础：坐标变换与材料参数，雅可比矩阵，等效介质理论
好，咱们今天聊点硬核的。隐身斗篷听起来像科幻，但它的数学根基其实非常扎实。说白了，就是一套坐标变换的把戏。我当年第一次接触这个理论时，也觉得玄乎，后来亲手算过一个案例才明白——嗯，其实就是给光“指路”，让它绕着物体走。
3.1 坐标变换的核心思想
想象一下，你手里有一张橡胶膜。上面画了一个点。你用手指把这个点戳到旁边去，橡胶膜就变形了。光在空间中的传播，本质上也可以这样“变形”。
变换光学的思路很简单：我们想要光在物理空间中走一条弯曲的路径，等价于在虚拟空间中走直线。你想想看，如果我能把空间“扭曲”掉，让光以为前方是空的，那它不就绕过物体了吗？
具体怎么做？我们需要一个映射关系：
原始空间 (x, y, z) → 变换后空间 (x', y', z')这个映射不是随便写的。它必须满足麦克斯韦方程的形式不变性。我在项目中遇到过一个问题：随便选了一个映射，结果算出来的材料参数是负值，根本造不出来。所以，映射的选择非常关键。
3.2 雅可比矩阵——变形程度的量化
有了映射，怎么描述“变形”的程度？雅可比矩阵登场了。
雅可比矩阵 J 的定义是：
J = ∂(x', y', z') / ∂(x, y, z)
即：
J_ij = ∂x'_i / ∂x_j这个矩阵的每个元素，表示变换后坐标对原始坐标的偏导数。说白了，就是告诉你：原始空间里一个小立方体，经过变换后变成了什么形状——是拉长了、压扁了，还是扭曲了。
我个人习惯把雅可比矩阵看作“变形指示器”。它的行列式 det(J) 表示体积变化率。如果 det(J) = 1，体积不变；如果 det(J) > 1，体积膨胀；如果 det(J) < 1，体积收缩。
关键公式：变换后的介电常数 ε' 和磁导率 μ' 由下式给出：
ε' = J · ε · J^T / det(J)
μ' = J · μ · J^T / det(J)其中 ε 和 μ 是原始空间的材料参数（通常是真空的 ε₀ 和 μ₀）。
你看，材料参数完全由雅可比矩阵决定。这就是变换光学的核心——设计映射 → 计算雅可比 → 得到材料参数。
3.3 等效介质理论——把复杂结构简化
好，现在问题来了：算出来的材料参数往往是各向异性的，而且数值很极端。比如，某个方向的介电常数是 2.5，另一个方向是 0.8。这种材料自然界不存在，怎么办？
等效介质理论就是来解决这个问题的。它的思想是：用亚波长结构的人工材料，来等效实现所需的电磁参数。
举个例子。我需要一个介电常数为 2.0 的材料。但我手头只有空气（ε=1）和陶瓷（ε=10）。怎么办？我把陶瓷做成小颗粒，均匀分散在空气中。只要颗粒尺寸远小于波长，这个混合物就表现出一个等效介电常数，介于 1 和 10 之间。
常用的等效介质公式有：
| 模型 | 公式 | 适用场景 |
|---|---|---|
| Maxwell-Garnett | ε_eff = ε_h + f·ε_i·(ε_i - ε_h)/(ε_h + L·(ε_i - ε_h)) | 低填充率，球形夹杂 |
| Bruggeman | f·(ε_i - ε_eff)/(ε_i + 2ε_eff) + (1-f)·(ε_h - ε_eff)/(ε_h + 2ε_eff) = 0 | 高填充率，双连续结构 |
| Lichtenecker | ln(ε_eff) = f·ln(ε_i) + (1-f)·ln(ε_h) | 随机分布，对数混合 |
我曾经踩过一个坑：用 Maxwell-Garnett 公式算了一个高填充率的案例，结果实验测出来偏差很大。后来才发现，填充率超过 30% 时，这个模型就不准了。换成 Bruggeman 模型才搞定。
避坑指南：等效介质理论的前提是结构尺寸远小于波长（通常要求 < λ/10）。如果结构太大，会出现衍射效应，等效参数就失效了。我建议你在仿真时先用本征模分析验证一下等效参数的有效性。
3.4 一个完整的计算流程
咱们把上面这些串起来，走一遍完整的流程。假设我们要设计一个二维的隐身斗篷，把半径为 R 的圆柱区域隐藏起来。
- 定义映射：将圆柱区域 (r < R) 映射到环形区域 (R₁ < r' < R₂)。常用的线性映射是 r' = R₁ + (R₂ - R₁)·r/R。
- 计算雅可比：在柱坐标下，雅可比矩阵是对角化的。算出来 J_rr = ∂r'/∂r = (R₂ - R₁)/R，J_θθ = r'/r，J_zz = 1。
- 得到材料参数：代入公式，得到 ε'_r = r/r'，ε'_θ = r'/r，ε'_z = (R₂ - R₁)²/R² · r/r'。μ 同理。
- 简化实现：这些参数是径向变化的，很难直接制造。可以用多层同心环结构来近似，每层用等效介质理论设计。
你看，步骤很清晰。但实际做的时候，每一步都有坑。比如，映射的边界条件要连续，否则会出现反射。我刚开始做时忽略了这一点，仿真结果里斗篷内部亮得跟灯泡似的——光全反射回来了。
3.5 知识体系总览
下面这张图，是我梳理的本章知识结构。你可以把它当作一个导航图，随时回来对照。
3.6 仿真中的实操要点
最后，分享几个我在仿真中积累的经验：
- 网格划分要匹配雅可比：变换后的区域，材料参数变化剧烈。我建议在径向方向加密网格，至少 10 层以上。
- 边界条件用 PML：斗篷外面一定要加完美匹配层，否则边界反射会污染结果。我吃过这个亏，仿真出来的隐身效果被反射波完全破坏了。
- 先做 2D 验证：别一上来就搞 3D。2D 圆柱斗篷的计算量小，调试方便。等 2D 跑通了，再扩展到 3D。
- 材料参数的截断处理：算出来的参数有时会趋于无穷大（比如 r' → 0 时）。实际实现时，需要做截断。我一般把最大值限制在 10 以内，效果损失不大。
注意：变换光学给出的材料参数，在理想情况下是完美的。但实际制造时，总会有误差。特别是等效介质结构的周期性和填充率，必须严格控制。偏差超过 5%，隐身效果就会明显下降。
好了，这一章的内容就到这儿。变换光学是隐身斗篷的理论基石，理解透了，后面的仿真和实验才能有的放矢。你想想看，我们其实是在用数学“欺骗”光——让它以为前方空无一物，实际上物体就在那里。这种思维，真的很妙。
公众号:蓝海资料掘金营，微信deep3321