近远场变换这个话题，我在做天线测量和电磁散射分析的那几年里反复啃过很多遍，每次卡住的地方几乎都绕不开同一个源头——麦克斯韦方程。很多人一上来就急着推近场到远场的积分公式，结果推到一半发现符号对不上、相位因子多出来一个、或者远场方向图怎么算都不对称，回头一查，问题往往出在对麦克斯韦方程组里那几个旋度方程的理解上。这篇内容就是想把“近远场变换”这条链路里麦克斯韦方程到底扮演什么角色讲透，从方程本身到它如何支撑起整个变换过程，把中间那些容易被跳过的推导细节和实操中的坑都摊开来说。不管你是刚接触天线近场测量的新手，还是已经能跑通仿真但对其所以然说不清楚的同行，应该都能从里面找到对自己有用的部分。核心关键词近远场变换、麦克斯韦方程、旋度方程、远场近似、等效原理会自然穿插在各个环节里，不会刻意堆砌。

1. 为什么近远场变换必须从麦克斯韦方程讲起

1.1 近远场变换到底在解决什么问题

先把场景摆清楚。天线测试里，我们很难直接在远场区测量，因为远场条件要求测试距离满足 r > 2D²/λ，一个口径一米、工作在几个吉赫兹的天线，远场距离轻松就到几十米甚至上百米，暗室根本放不下。所以工程上的常规做法是在近场区（通常是几个波长到十几个波长的距离）用探头扫一个面——平面、柱面或者球面——把近场的幅度和相位采下来，然后通过近远场变换算法，把这组近场数据“推”到远场去，得到远场方向图。这个“推”的过程，本质上就是求解一个电磁边值问题：已知某个闭合面（或开放面）上的切向场，求外部无源区域的辐射场。

而描述电磁场在空间中如何分布、如何传播、如何由源产生的，就是麦克斯韦方程。近远场变换不是某种经验拟合，它是有严格电磁理论支撑的数学变换，这个支撑的根就在麦克斯韦方程组里。你可以把近场测量数据看作“因”，远场方向图看作“果”，连接因果的那条逻辑链，就是麦克斯韦方程加上等效原理。跳过方程直接背变换公式，等于记住了药方却不知道药理，一旦测量面不是标准平面、或者探头补偿没做对，就完全不知道该从哪里排查。

1.2 从方程到变换的逻辑链条

整条逻辑链我习惯这样梳理：麦克斯韦方程组描述场与源的关系；在无源区域（近场测量面之外、远场观察点所在的空间），方程组退化为齐次形式，场满足矢量亥姆霍兹方程；利用等效原理，把测量面上的切向场等效为一层面电流和一层面磁流；这些等效源在外部区域产生的场，可以用矢量位或者直接对等效源做积分求得；最后对积分施加远场近似（r 远大于源尺寸和波长），得到远场方向图与近场分布之间的傅里叶变换关系。

这条链上每一环都依赖麦克斯韦方程：齐次波动方程来自方程组的旋度运算，等效原理来自方程组的边界条件形式，远场近似里场的横向性（远场只有 θ 和 φ 分量，没有 r 分量）同样来自方程组在远区的渐近行为。所以我说，麦克斯韦方程不是近远场变换的“背景知识”，而是它的“骨架”。理解了骨架，公式就是自然长出来的肉；不理解骨架，公式就是一堆需要死记的符号。

1.3 一个常见的认知误区

我见过不少同行，包括早期的我自己，把近远场变换当成纯粹的信号处理问题——近场数据做一次二维 FFT 就出远场，跟电磁理论没关系。这个理解在平面近场、标准采样、探头理想的情况下确实能跑通，因为平面近远场变换在数学上恰好表现为傅里叶变换。但一旦遇到柱面或球面扫描，或者需要做探头方向图补偿，纯 FFT 的思路立刻失效，必须回到等效源和辐射积分的框架里。而那个框架的每一块砖，都是麦克斯韦方程烧出来的。所以我建议，哪怕你日常只用平面近场，也把方程层面的推导过一遍，遇到异常数据时你会感谢自己多花的这几个小时。

2. 麦克斯韦方程组在变换中的核心角色拆解

2.1 四个方程各自管什么

先把方程组摆出来，用频域、时谐约定 e^{jωt} 的形式（这是天线领域最常用的约定，注意和有些教材的 e^{-jωt} 约定不同，符号差异会直接影响后续相位因子的正负，这个坑后面会专门讲）：

• ∇ ×E= -jωμH（法拉第定律）
• ∇ ×H=J+ jωεE（安培定律含位移电流）
• ∇ ·E= ρ/ε（高斯定律）
• ∇ ·H= 0（磁通连续性）

在近远场变换关心的无源区域，J= 0、ρ = 0，前两个方程变成 ∇ ×E= -jωμH和 ∇ ×H= jωεE，后两个变成 ∇ ·E= 0、∇ ·H= 0。这四个齐次方程就是整个变换的理论起点。法拉第定律和安培定律负责“旋度”层面的耦合，说明电场和磁场互相激发、互相缠绕；两个散度方程负责“无源”约束，说明在无源区场线闭合、没有起点终点。近远场变换里所有关于场传播、模式展开、远场横向性的结论，都是从这四个方程推出来的。

2.2 旋度方程如何导出矢量波动方程

把法拉第定律两边取旋度：∇ × ∇ ×E= -jωμ ∇ ×H。右边代入安培定律 ∇ ×H= jωεE，得到 ∇ × ∇ ×E= -jωμ · jωεE= ω²μεE。再用矢量恒等式 ∇ × ∇ ×E= ∇(∇ ·E) - ∇²E，而无源区 ∇ ·E= 0，于是得到 ∇²E+ k²E= 0，其中 k² = ω²με。这就是矢量亥姆霍兹方程。磁场同理。

这一步看着是纯数学操作，但它的物理含义很关键：它告诉我们，无源区里的每一个场分量都满足标量亥姆霍兹方程，而标量亥姆霍兹方程的解是平面波 e^{-jk·r} 的叠加。近远场变换里把近场展开成平面波谱、柱面波谱或球面波谱，理论依据就是这一步。换句话说，你在变换算法里做的那些模式展开，不是数学技巧，而是波动方程通解的必然形式。理解了这一点，你就明白为什么平面近场变换天然对应二维傅里叶变换——因为平面波谱展开本身就是傅里叶变换的形式。

2.3 边界条件与等效原理的方程来源

等效原理常被当成一个独立的定理来讲，其实它直接来自麦克斯韦方程的边界条件。在两种介质的分界面上，切向电场连续、切向磁场连续（无表面电流时），或者切向场的不连续量等于表面电流/磁流密度。近远场变换的做法是：在测量面 S 上，把实际场的切向分量E_t 和H_t 取出来，然后假设 S 面以内场为零、S 面以外场不变，为了维持这个假设，必须在 S 面上放置等效面电流J_s =n×H_t 和等效面磁流M_s = -n×E_t（n为指向外部区域的单位法向）。这组等效源在外部区域产生的场，与实际场完全一致。

这个等效关系的合法性，完全建立在方程组的边界条件形式上。所以当有人问我“等效原理为什么成立”，我的回答永远是：它不是假设，它是麦克斯韦方程在分界面上的必然结论。你在代码里构造等效源的那几行，实际上是在数值上实现边界条件。

2.4 远场近似里方程给出的横向性

远场区 r → ∞ 时，场的行为有一个非常重要的性质：E和H都垂直于传播方向r̂，且H= (1/η)r̂×E，η 为波阻抗。这个横向性不是额外假设，它可以从旋度方程在远区的渐近展开严格推出。直观理解：远场局部看就是平面波，而平面波在无源区满足k·E= 0（来自散度方程）和H= (1/η)k̂×E（来自旋度方程）。近远场变换最后一步把积分化简成只含 θ、φ 分量，靠的就是这个性质。如果你在代码里发现远场结果多出一个径向分量，那基本可以断定是远场近似没施加到位，或者等效源构造时法向取错了。

3. 从方程到可执行变换的完整推导路径

3.1 矢量位方法：把等效源变成辐射积分

直接对等效源积分求场，最常用的中介是矢量位。定义磁矢量位A和电矢量位F：

• A(r) = (μ/4π) ∫J_s(r') e^{-jkR}/R dS'
• F(r) = (ε/4π) ∫M_s(r') e^{-jkR}/R dS'

其中 R = |r-r'|，r是观察点，r' 是源点。这两个积分就是等效源在观察点产生的矢量位，核函数 e^{-jkR}/R 是标量亥姆霍兹方程在无界空间的格林函数——又一次回到波动方程。得到A和F后，场由下式给出：

E= -jωA- j(1/ωμε) ∇(∇ ·A) - (1/ε) ∇ ×F

H= -jωF- j(1/ωμε) ∇(∇ ·F) + (1/μ) ∇ ×A

这两个式子是矢量位到场的标准转换，推导过程同样只用旋度方程和散度方程。实操中我们很少真的去算 ∇(∇ ·A) 这种二阶微分，而是利用远场近似把微分算子简化成代数运算，这就是下一步。

3.2 远场近似：微分算子如何变成乘法

远场条件下 R ≈ r -r̂·r'，1/R ≈ 1/r，格林函数近似为 e^{-jkr}/r · e^{jkr̂·r'}。把这个近似代入矢量位积分，再代入场的转换式，经过一番化简（中间用到 ∇ → -jkr̂的远场替换，即远场区微分算子等效为乘以 -jkr̂），最终得到：

E_far(r) ≈ -jk e^{-jkr}/(4πr) · [ η(J_s 的横向分量) + (M_s 的横向分量) ×r̂] 的积分形式

更常见的写法是引入辐射矢量N和L：

N= ∫J_s e^{jkr̂·r'} dS'，L= ∫M_s e^{jkr̂·r'} dS'

则远场E= -jk e^{-jkr}/(4πr) · (ηN+L×r̂) 的横向部分。这里 e^{jkr̂·r'} 就是平面波谱展开里的相位核，对平面近场来说，r̂的横向分量与扫描面上的波数对应，整个积分恰好是二维傅里叶变换。这就是“平面近场变换等于 FFT”这一说法的严格来源——它不是巧合，是远场近似下辐射积分自然退化成的形式。

3.3 平面、柱面、球面三种扫描的方程视角差异

三种扫描面的变换公式形式不同，但从麦克斯韦方程的角度看，差异只在于用哪种坐标系下的波动方程通解来展开近场。

平面扫描之所以最常用，是因为它的展开最直接、FFT 最快，但对波束很宽的天线，平面扫描截断误差大，因为宽波束的能量在扫描面边缘还没衰减完。柱面和球面扫描用对应的波谱展开，本质上都是把近场投影到该坐标系下波动方程的本征函数基上，再传播到远场。选哪种，取决于天线的波束宽度和你的暗室条件，不是随便挑的。

3.4 相位约定：一个必须提前钉死的细节

前面提过 e^{jωt} 和 e^{-jωt} 两种约定。用 e^{jωt} 时，法拉第定律是 ∇ ×E= -jωμH，平面波是 e^{-jk·r}；用 e^{-jωt} 时，法拉第定律变成 ∇ ×E= jωμH，平面波是 e^{+jk·r}。两种约定下，远场积分里的相位核指数符号相反。如果你在推导时混用了两种约定，最典型的表现是远场方向图的相位出现整体翻转或者镜像，幅度可能看起来正常，但相位一塌糊涂。我的做法是：项目一开始就在文档头写死约定，所有公式、代码注释统一，绝不中途换。这个习惯帮我省过至少两次通宵排查。

4. 实操落地：把方程变成能跑的代码

4.1 数据准备与预处理

近场测量数据通常是扫描面上每个采样点的复数切向场，两个正交极化分量。拿到数据后，第一件事是检查采样间隔是否满足奈奎斯特条件：平面扫描时，采样间隔 Δx 必须小于 λ/2，否则平面波谱会出现混叠，远场方向图在大角度区域会冒出虚假的旁瓣。我实测过，Δx 取 0.4λ 到 0.45λ 比较稳妥，留一点余量应对探头尺寸带来的等效采样展宽。

预处理还包括：去掉直流偏置（如果有）、对齐相位参考、检查扫描面边缘场强是否已衰减到足够低（一般要求边缘场比峰值低 30 dB 以上，否则截断误差会污染远场低旁瓣）。这些检查不做，后面方程推得再对，结果也不可信。

4.2 等效源构造的代码实现

以平面扫描为例，扫描面在 z = 0，法向n=ẑ（指向远场所在半空间）。等效面电流和磁流：

import numpy as np # Ex, Ey: 扫描面上采到的切向电场，形状 (Ny, Nx) # Hx, Hy: 扫描面上采到的切向磁场 # n = z_hat，指向 +z 半空间 Js_x = Hy # (n × H)_x = Hy Js_y = -Hx # (n × H)_y = -Hx Ms_x = -Ey # -(n × E)_x = -Ey Ms_y = Ex # -(n × E)_y = Ex

注意符号：J_s =n×H_t，M_s = -n×E_t，这两个式子的符号错了，远场极化会整体旋转 90 度或者反相。我建议第一次实现时，用一个已知的简单源（比如一个短偶极子的解析近场）做自检，确认等效源构造和后续积分能复现解析远场，再上真实数据。

4.3 平面近场到远场的 FFT 实现

平面扫描下，辐射矢量N和L的积分在离散化后就是二维 FFT：

# kx, ky: 波数域采样，由扫描面尺寸和采样间隔决定 # kx = 2*pi*np.fft.fftfreq(Nx, dx) # ky = 2*pi*np.fft.fftfreq(Ny, dy) Nx_spec = np.fft.fft2(Js_x) * dx * dy Ny_spec = np.fft.fft2(Js_y) * dx * dy Lx_spec = np.fft.fft2(Ms_x) * dx * dy Ly_spec = np.fft.fft2(Ms_y) * dx * dy # 可见区掩膜：kx^2 + ky^2 <= k^2 k = 2*np.pi/lambda_ kz = np.sqrt(k**2 - kx**2 - ky**2 + 0j) # 复数，隐失波取衰减 visible = (kx**2 + ky**2) <= k**2 # 远场分量（省略常数因子，实际需乘 -jk*exp(-jkr)/(4*pi*r)） # 这里演示核心结构，具体分量组合按 3.2 节公式

关键点：kz必须用复数开方，可见区内是实数（传播波），可见区外是纯虚数（隐失波，指数衰减，对远场无贡献，应置零或保留但不参与远场）。很多新手直接用实数开方，可见区外出现 NaN，或者忘了掩膜，隐失波混进远场，导致大角度方向图异常抬升。

4.4 探头补偿的方程依据

实际测量用的是有方向图的探头，不是理想点源。探头在扫描面上每个位置接收到的信号，是探头方向图与当地场分布的卷积。补偿的方法是在波数域除以探头方向图（或者用探头传输函数的逆），这一步的依据是：探头接收过程可以建模为当地场与探头等效接收方向图的乘积在探头口径上的积分，在平面波谱域表现为乘法。所以补偿在谱域做除法最干净。但要注意，探头方向图在可见区边缘可能接近零，直接除会放大噪声，实操中会在低信号区域加正则化或者截断。

4.5 一个完整的自检流程

我习惯用这套流程验证实现是否正确：第一步，构造一个解析已知的源（比如均匀口径分布），用解析公式算近场，喂给变换代码，比对远场主瓣宽度和第一旁瓣电平；第二步，把近场数据沿扫描面平移半个波长，远场相位应该线性变化，幅度不变；第三步，改变扫描面尺寸，远场主瓣应稳定，只有旁瓣细节随截断变化。三步都过，基本可以认为方程到代码的链路是通的。这套自检我每次重构变换模块都会跑一遍，比直接上暗室数据省时间得多。

5. 常见问题与排查实录

5.1 远场方向图不对称

最典型的原因有两个：一是等效源符号错，导致某个极化分量反相，方向图出现镜像不对称；二是扫描面采样中心没对准天线相位中心，引入线性相位斜坡，远场波束偏移。排查方法：先看近场相位分布是否关于中心对称，如果有明显斜坡，先做相位中心校正；再看等效源构造的符号，用 4.2 节的自检源验证。

5.2 大角度区域出现虚假抬升

多半是采样间隔不够（混叠）或者隐失波没滤干净。检查 Δx 是否小于 λ/2，检查 FFT 后可见区掩膜是否施加。还有一种可能是扫描面边缘场强没衰减够，截断导致谱泄漏，表现为大角度区域整体抬升，这种只能靠加大扫描面或者加窗抑制。

5.3 相位整体翻转或镜像

相位约定混用，或者 FFT 的移位操作（fftshift）用错。检查推导和代码是否统一用 e^{jωt} 约定，检查相位核指数符号，检查 fftshift 是否在正确的位置调用。

5.4 常见问题速查表

5.5 几条踩坑心得

第一，永远先用解析源自检，别拿真实数据当第一个测试用例，真实数据出问题你分不清是算法错还是数据错。第二，相位约定写在代码注释第一行，团队协作时这条能救命。第三，可见区掩膜和隐失波处理是平面近场变换里最容易漏的一步，漏了不会报错，只会悄悄污染结果。第四，探头补偿的正则化参数要试，除得太狠噪声放大，除得太松补偿不到位，一般从信噪比最差的频点开始调。

6. 方程视角下的扩展与边界

6.1 为什么球面近场变换更“方程化”

球面扫描的变换用球面波谱展开，本质是把近场投影到球谐函数基上，每个模式对应一个球面波本征解。这个展开直接来自球坐标下矢量亥姆霍兹方程的分离变量解，方程味最浓。它的好处是对任意波束宽度都稳，坏处是模式数随电尺寸平方增长，计算量大。理解球面展开，反过来能加深对平面展开的理解——平面波谱不过是球面波谱在局部近似下的极限形式。

6.2 近场变换与逆散射的方程同源性

近远场变换是“已知面上场求外部场”，逆散射是“已知外部场求内部源或介质分布”，两者都建立在麦克斯韦方程和等效原理上，只是求解方向相反。做过近远场变换的人转去做逆问题，会发现很多工具是共通的：等效源、格林函数、波谱展开。这也是我建议把方程基础打牢的原因——它不只服务一个任务，而是一整类电磁问题的共同底座。

6.3 数值实现里方程被“隐藏”的地方

现代商业软件把方程都封装好了，用户点几下就出远场。这当然高效，但也让人失去对异常的判断力。我的经验是，至少自己手写一遍平面近场变换的核心几十行，把等效源、FFT、掩膜、远场分量组合都亲手敲一遍，跑通一次解析自检。这几十行代码带来的理解，比读十篇综述都实在。方程不会因为你用了软件就消失，它只是藏在了按钮后面，出问题时还得靠它来定位。

最后分享一个我自己的习惯：每次变换结果异常，我先不看代码，先在纸上把等效源和远场积分的关系写一遍，确认符号和相位核，再回去对代码。十次里有七八次，问题在纸上就暴露了。方程这东西，平时觉得抽象，真到排查的时候，它是最靠得住的参照系。