简介：这份资源围绕无穷小偶极子天线这一经典电磁学模型，提供基于Matlab的完整仿真实现，面向电子信息、通信工程等专业的本科生、研究生及教研人员，用于理解天线近场与远场的辐射特性、场强分布规律及数值计算方法。压缩包共6个文件，约1.08MB，以m脚本文件为核心，配合png结果图、说明文档与嵌套压缩包，便于直接运行、对照结果并快速上手。资源已在平台获得24人学习下载，可作为课程设计、实验报告或科研入门的参考素材。读者可借助main.m脚本复现偶极子天线近区与远区的电磁场计算流程，结合生成的图像直观观察场强随距离与角度的变化趋势，并参考说明文档理解参数设置与仿真思路，从而掌握从理论公式到Matlab编程实现的完整链路，为后续天线阵列、电磁兼容等方向的仿真打下基础。
1. 无穷小偶极子天线的近场远场：为什么它是天线仿真的第一块试金石
很多人第一次在 Matlab 里画偶极子天线方向图，都会下意识把 r 设成一个很大的值，然后心安理得地认为"远场嘛，肯定准"。结果把 r 从 100λ 缩到 0.1λ，方向图几乎没变，于是得出结论：近场远场差不多。这个结论是错的，错在把"无穷小偶极子"当成了普通半波振子。无穷小偶极子（infinitesimal dipole，也叫赫兹偶极子）长度 dl 远小于波长，电流沿长度均匀分布，它是唯一能写出严格闭式解的辐射源，也是所有天线教材推导近场、远场分界的起点。Matlab 在这里的价值不是画一张好看的 3D 图，而是让你亲手把 1/r³、1/r²、1/r 三项拆开，看清哪一项在哪个距离上主导。这套东西适合三类人：正在学天线理论但公式记不住的学生、需要验证 HFSS 或 CST 远场导出数据是否合理的一线工程师、以及要给近场探头或 NFC 线圈做理论基准的人。下面我按"公式怎么落到代码、参数怎么设、哪里会翻车"的顺序讲一遍。
2. 从赫兹偶极子闭式解到 Matlab 可执行公式：先立住理论再动手
2.1 三个场分量和三个距离区间
无穷小偶极子沿 z 轴放置，长度 dl，电流 I₀，位于原点。在球坐标 (r, θ, φ) 下，磁场只有 φ 分量，电场有 r 和 θ 两个分量。标准闭式解（Balanis《Antenna Theory》第 4 章）写成：
- H_φ = j·(I₀·dl·sinθ / 4π) · (k/r) · (1 + 1/(jkr)) · e^(-jkr)
- E_r = η·(I₀·dl·cosθ / 2π) · (1/r²) · (1 + 1/(jkr)) · e^(-jkr)
- E_θ = j·η·(I₀·dl·sinθ / 4π) · (k/r) · (1 + 1/(jkr) - 1/(kr)²) · e^(-jkr)
其中 k = 2π/λ 是波数，η ≈ 377 Ω 是自由空间波阻抗。关键在括号里那几项：E_r 只有 1/r² 和 1/r³ 量级，没有辐射项；E_θ 同时含 1/r（辐射）、1/r²（感应）、1/r³（静场）三项；H_φ 只有 1/r 和 1/r²。这就是近场远场分界的物理来源——不是人为规定，而是三项随距离衰减速度不同。
工程上习惯用 kr 划界：kr ≪ 1 是近场（静场和感应场主导），kr ≫ 1 是远场（辐射场主导），kr ≈ 1 即 r = λ/(2π) ≈ 0.159λ 是过渡区。注意这个 0.159λ 和常说的"远场边界 2D²/λ"不是一回事，后者是针对有限尺寸天线的口径判据，对无穷小偶极子这种点源，2D²/λ 会退化成 0，没有意义。这是第一个容易混的点。
2.2 把公式写成可复用的 Matlab 函数
我一般不会把公式直接塞进脚本，而是写成一个返回三个分量的函数，方便后面扫距离、扫角度。下面这段在 R2016b 以后的版本都能跑，注意用 element-wise 运算：
function [Er, Etheta, Hphi] = infDipoleField(r, theta, I0, dl, lambda) % 无穷小偶极子近场/远场闭式解 % r : 观察点距离 (m)，可为向量 % theta : 极角 (rad)，可为向量 % I0 : 电流幅值 (A) % dl : 偶极子长度 (m)，需满足 dl << lambda % lambda : 工作波长 (m) k = 2*pi/lambda; eta = 377; % 自由空间波阻抗 kr = k .* r; % 防止 r=0 导致除零，实际使用时 r 从 0.001*lambda 起 Er = eta * I0 * dl * cos(theta) ./ (2*pi) .* ... (1./r.^2) .* (1 + 1./(1j*kr)) .* exp(-1j*kr); Etheta = 1j * eta * I0 * dl * sin(theta) ./ (4*pi) .* ... (k./r) .* (1 + 1./(1j*kr) - 1./(kr).^2) .* exp(-1j*kr); Hphi = 1j * I0 * dl * sin(theta) ./ (4*pi) .* ... (k./r) .* (1 + 1./(1j*kr)) .* exp(-1j*kr); end逻辑说明：kr 单独算出来是为了让三个括号项写起来干净，也方便你后面单独把 1/(kr)² 那一项拎出来看。参数上，dl 必须远小于 lambda，我一般取 dl = lambda/1000 到 lambda/50，超过 lambda/10 闭式解就开始失真，因为电流均匀分布的假设崩了。I0 取 1 A 只是为了归一化，实际算方向图时幅值会被约掉。r 不能取 0，物理上原点处场发散，代码里从 0.001λ 起步。
2.3 用一段脚本验证三项主导关系
光看公式没感觉，跑一段扫距离的脚本，把 |E_θ| 的三个分量分别画出来，你会立刻明白近场远场是怎么回事：
lambda = 1; I0 = 1; dl = lambda/100; r = logspace(-3, 2, 500) * lambda; % 0.001λ 到 100λ theta = pi/2; % 赤道面，sinθ=1，E_θ 最大 k = 2*pi/lambda; eta = 377; kr = k*r; % 三项分别计算（去掉公共系数） C = eta * I0 * dl / (4*pi) * k; rad = C ./ r; % 1/r 辐射项 ind = C ./ (r .* kr); % 1/r² 感应项 stat = C ./ (r .* kr.^2); % 1/r³ 静场项 loglog(r/lambda, abs(rad), 'r', r/lambda, abs(ind), 'g', ... r/lambda, abs(stat), 'b', 'LineWidth', 1.5); xlabel('r / \lambda'); ylabel('|E_\theta| 分量幅值'); legend('1/r 辐射','1/r^2 感应','1/r^3 静场','Location','southwest'); grid on; xline(1/(2*pi), 'k--', 'kr=1');跑完你会看到三条斜率不同的直线，在 r ≈ 0.159λ 处感应项和辐射项交叉，静场项在更近处就掉下去了。这就是"近场远场"最直观的证据。参数说明：logspace 保证近场区域采样够密，否则 0.001λ 到 0.01λ 这段会被挤成一条线；theta 取 pi/2 是因为 E_r 在赤道面为零，只剩 E_θ，分析最干净。
3. 近场和远场方向图到底差在哪：用 Matlab 把差异画出来
3.1 远场方向图：标准的 sinθ 甜甜圈
远场（kr → ∞）时，括号里只剩 1，E_θ 退化成：
E_θ ≈ j·η·(I₀·dl·k / 4π) · (sinθ / r) · e^(-jkr)
归一化方向图就是 |sinθ|，与 φ 无关，3D 形状是个甜甜圈。Matlab 里用 meshgrid 扫 θ 和 φ 就能画：
theta = linspace(0, pi, 181); phi = linspace(0, 2*pi, 361); [TH, PH] = meshgrid(theta, phi); % 远场归一化方向图 F = abs(sin(TH)); % 转直角坐标画 3D X = F .* sin(TH) .* cos(PH); Y = F .* sin(TH) .* sin(PH); Z = F .* cos(TH); surf(X, Y, Z, F, 'EdgeColor', 'none'); axis equal; colormap jet; colorbar; title('无穷小偶极子远场方向图');这段没什么坑，唯一注意 theta 从 0 到 pi 而不是 0 到 2pi，否则 z 轴会重复画一遍。远场方向图的最大值在 θ = 90°，零点在 θ = 0° 和 180°，前后比无穷大（理想点源没有后瓣）。
3.2 近场方向图：形状随距离变形
近场里 E_θ 的括号项 (1 + 1/(jkr) - 1/(kr)²) 是个复数，模和相位都随 kr 变。这意味着近场方向图不只是幅值缩放，形状本身会变。我一般固定 r，扫 θ，把 |E_θ| 归一化后和远场叠在一起看：
lambda = 1; I0 = 1; dl = lambda/100; eta = 377; k = 2*pi/lambda; theta = linspace(0.01, pi-0.01, 300); r_list = [0.05, 0.159, 0.5, 5] * lambda; % 近场到远场 figure; hold on; for r = r_list kr = k*r; E = abs(1j * eta * I0 * dl * sin(theta) / (4*pi) * (k/r) * ... (1 + 1/(1j*kr) - 1/(kr)^2) * exp(-1j*kr)); plot(theta*180/pi, E/max(E), 'LineWidth', 1.5); end legend('r=0.05\lambda','r=0.159\lambda','r=0.5\lambda','r=5\lambda'); xlabel('\theta (deg)'); ylabel('归一化 |E_\theta|'); grid on;跑出来你会发现 r = 0.05λ 时方向图明显比 sinθ 窄，r = 0.159λ 时接近但还没完全重合，r = 5λ 时和 sinθ 几乎分不出差别。这就是近场方向图"变形"的量化表现。参数上 r_list 特意选了 0.159λ 这个 kr=1 的点，它是过渡区的标志。注意归一化用的是每个 r 各自的最大值，如果你想看绝对幅值差异，去掉 /max(E) 即可，但那样远场曲线会被压到看不见。
3.3 波阻抗：近场不是 377 Ω
远场里 E_θ/H_φ = η = 377 Ω，这是大家背熟的。但近场里这个比值是复数，且随 kr 变。算一下：
kr = logspace(-2, 2, 200); Zw = 377 * (1 + 1./(1j*kr) - 1./(kr).^2) ./ (1 + 1./(1j*kr)); semilogx(kr, abs(Zw), 'b', kr, real(Zw), 'r--', kr, imag(Zw), 'g-.'); xlabel('kr'); ylabel('波阻抗 (\Omega)'); legend('|Z_w|','Re','Im'); grid on; xline(1,'k--');kr ≪ 1 时 |Z_w| 远大于 377，且虚部主导，说明近场以储能为主，不是辐射；kr ≫ 1 时收敛到 377 实数。这个图是判断"我到底在近场还是远场"最硬的依据，比记 0.159λ 更靠谱，因为不同天线结构过渡点不一样。
4. 避坑与排查：无穷小偶极子仿真里最容易翻车的 5 个地方
4.1 现象：方向图在近场出现"多瓣"，以为算错了
原因：把 E_r 和 E_θ 的模直接相加，或者用 abs(Er) + abs(Etheta) 当总场。E_r 和 E_θ 是正交分量，总场模应该是 sqrt(|Er|² + |Eθ|²)，而且它们相位不同，直接标量相加会造出假瓣。
解决：算总场用 sqrt(abs(Er).^2 + abs(Etheta).^2)，或者干脆只看 E_θ（赤道面 E_r=0）。我一般分析方向图时固定 θ 扫 φ 或固定 φ 扫 θ，避免分量混叠。
4.2 现象：r 取 0.001λ 时结果爆炸或 NaN
原因：公式里有 1/r³ 项，r 太小时浮点数溢出，或者 kr 小到 1/(kr)² 超过 double 精度。
解决：r 下限设成 0.001λ 到 0.005λ，再小物理上也没意义（偶极子本身有尺寸 dl，r < dl 时点源近似失效）。代码里加一句 r = max(r, 0.001*lambda) 兜底。
4.3 现象：dl 取 lambda/4 结果和教材对不上
原因：无穷小偶极子的前提是 dl ≪ λ 且电流均匀。dl 到 lambda/4 时电流分布已经接近正弦，均匀假设崩了，闭式解不再成立。
解决：dl 控制在 lambda/50 以内。如果你要算半波振子，得用正弦电流分布积分，不能套这个公式。这是标题里"无穷小"三个字的硬约束。
4.4 现象：用 2D²/λ 算远场边界，得到 0 或负数
原因：无穷小偶极子 D → 0，2D²/λ 判据失效。这个判据是给口径天线用的。
解决：点源类天线用 kr ≥ 10（即 r ≥ 1.6λ）作为远场工程判据，或者直接看波阻抗是否收敛到 377。别硬套 2D²/λ。
4.5 现象：Matlab 2023 中文注释乱码
原因：R2023 之前默认编码是 GBK，脚本存成 UTF-8 后中文注释变问号。热词里"matlab 2023 的中文注释乱码"就是这个。
解决：在 R2023b 及以后，用 matlab.editor 设置里把编码改成 UTF-8；老版本在脚本开头加 %% 后手动另存为 GBK，或者干脆注释写英文。我现在的习惯是函数头注释用英文，正文说明写在外面，省得跨版本折腾。
5. 进阶：把近场数据导出给 HFSS 对比，以及一个验证小技巧
实际项目里，Matlab 算出来的近场经常要拿去和 HFSS 或 CST 的场监视器结果对比。热词里"hfss导出远场数据"就是这个场景。HFSS 导出的远场数据一般是 .csv 或 .fld，列是 theta、phi、|E|、相位。我的做法是：Matlab 这边在同样的 r、theta、phi 网格上算一遍，两边都归一化到各自最大值，然后算相对误差。如果远场区域误差在 5% 以内，说明模型对得上；如果近场误差大，先检查 HFSS 里场监视器的 r 是不是设在了辐射边界内部——很多人把监视器放在 PML 里面，数据本身就是吸收边界反射后的，当然对不上。
一个我常用的验证技巧：不画方向图，直接画 |E_θ| 随 r 的衰减斜率。远场区域 log-log 斜率应该是 -1（1/r），近场静场区是 -3（1/r³），过渡区在 -1 和 -3 之间。用 polyfit 对 log(r) 和 log(|E|) 做一次线性拟合，斜率对不上就说明你的 r 区间选错了，或者 dl 太大破坏了点源假设。这个比看方向图灵敏得多，方向图在过渡区变化很缓，肉眼容易骗过去。
% 衰减斜率验证 r = logspace(-2, 2, 200) * lambda; kr = 2*pi/lambda * r; E = abs(1./(r) .* (1 + 1./(1j*kr) - 1./(kr).^2)); % 去掉常数系数 p = polyfit(log10(r), log10(E), 1); fprintf('拟合斜率 = %.3f (远场应趋近 -1，近场趋近 -3)\n', p(1));跑出来你会看到整段拟合斜率大概在 -1.5 左右，因为跨越了过渡区。正确做法是分段拟合：r < 0.05λ 取一段，r > 2λ 取一段，分别看斜率。我踩过的坑就是一开始拿整段拟合，斜率 -1.5 看着"挺合理"，实际上把两个区域平均掉了，什么问题都发现不了。分段之后近场段 -2.9、远场段 -1.02，才算真正验证通过。
这套东西值不值得做？如果你只是要一张方向图交作业，直接画 sinθ 就行，五分钟的事。但如果你要建立对近场远场的物理直觉、要验证商业软件的场数据、要给近场耦合问题打基准，那把这三个分量拆开跑一遍是绕不过去的。我自己的习惯是每换一个天线结构，先用无穷小偶极子的闭式解跑一遍衰减斜率和波阻抗，确认自己的观察距离落在哪个区，再去跑全波仿真。这个基准花不了半小时，但能省掉后面几天的"为什么仿真结果和实测对不上"的排查。希望帮到你。
本文还有配套的精品资源，点击获取
