简介本资源是面向本科及硕士阶段电磁场与天线课程教学与科研实践的Matlab仿真案例聚焦无穷小偶极子天线近场与远场辐射特性的数值建模与可视化分析。通过理论公式推导与Matlab编程实现直观呈现电场/磁场空间分布、方向图、场强衰减规律等核心物理特性助力理解天线辐射机理与近远场划分标准。压缩包共6个文件1.08MB含主程序main.m、3幅关键结果图像png、说明文档txt及推广引导图结构简洁、开箱即用所有代码兼容Matlab 2014a/2019a附运行结果截图降低调试门槛。目前已有23人学习下载适合电磁场仿真入门者快速掌握天线场计算流程并可延伸至阵列设计、方向图综合等进阶应用。 做天线仿真这个事很多人第一次上手就卡在“理论公式我全记住了可一写代码就懵”。大一学电磁场时老师让我们用 Matlab 画一个无穷小偶极子天线的近场和远场我当时觉得这有什么难的直接套公式不就行了。结果真动手才发现坐标变换、复数运算、近场远场切换条件、三维可视化每一个环节都有坑。等我全部调通把那张漂亮的 3D 方向图截下来交作业时我才算真正把“辐射”“波阻抗”“归一化方向函数”这些概念连成了一条线。这个项目看起来只是一个小小的仿真作业但它实际上是把《电磁场与电磁波》和《天线原理》两门课里最重要的概念——偶极子辐射、近场感应、远场辐射、方向图、功率密度——全部串起来的一次实操。适合三类人看正在学电磁场与天线课程的学生、刚入门的射频工程师以及想把 Matlab 数值计算能力用在电磁仿真上的硬件爱好者。这篇文章就把我当时的完整思路、代码、踩坑记录全部摊开讲清楚。1. 为什么无穷小偶极子天线是天线仿真的必修课1.1 一个被低估的“标准模型”很多初学者不理解为什么要在“无穷小偶极子”上下功夫——现实中哪有无穷小的天线我又不能拿它去发 WiFi 信号学了有什么用这里我得说个实在话无穷小偶极子Hertzian Dipole不是让你直接用的天线而是从数学上逼出来的最小辐射单元。它的定义很简单一段长度 \( l \) 远小于工作波长 \( \lambda \) 的直导线上面流过一个均匀的时谐电流 \( I_0 \)电流幅度和相位处处相同。因为 \( l \ll \lambda \)所以可以把整个单元看成“一个点”这就是“无穷小”这个说法的来源。这个模型最大的价值在于它是唯一一个能用初等函数写出严格解析解的辐射系统。随便翻一本天线教材第一个推导的就是赫兹偶极子。你把它的场解搞明白了后面所有线天线、阵列天线、面天线本质上都是把这个“最小单元”按照不同方式叠加起来。你可以把无穷小偶极子理解成天线世界的“积木块”单块可能长得不起眼但它是整个摩天大楼的基础。1.2 从麦克斯韦方程组到工程建模的桥梁我当时做这个题目时有个很深的感受课本上麦克斯韦方程组写得漂漂亮亮真让你拿来算天线辐射你根本不知道从哪下手。无穷小偶极子恰好提供了这条路径先引入磁矢位 \( \mathbf{A} \)通过亥姆霍兹方程解出 \( A_z \)再借助 \( \mathbf{E} -j\omega\mathbf{A} - \nabla\Phi \) 和 \( \mathbf{H} \frac{1}{\mu}\nabla \times \mathbf{A} \) 这两个关系式一步步把电场和磁场分量抠出来。正是这套推导让我理解了为什么工程上要用“矢量位”而不是直接解电场。矢量位在球坐标系下只有一个 \( z \) 方向分量方程形式极其简单而直接去解 \( \mathbf{E} \)要面对的是三个耦合分量手推一遍就知道有多痛苦。仿真也是一样你选对了中间变量代码就会好写很多数值稳定性和可读性也会完全不同。另一个重要收获是“模型简化”思维。无穷小偶极子的 \( l \ll \lambda \) 假设把辐射特性收敛到“点源辐射”这一极限情形没有方向性副瓣、没有互耦、没有阻抗带宽这些工程因素干扰。这样你研究近场和远场的物理本质时看到的就是最纯粹的电磁波传播行为不会被次要因素带偏。2. 远场与近场的数学表达及关键参数拆解2.1 完整电场与磁场表达式做这个项目你必须在代码里同时处理近场和远场所以我先把理论公式完整写出来。设偶极子沿 \( z \) 轴放置中心位于原点电流 \( I_0 \) 方向为正 \( z \)长度为 \( l \)则球坐标系中的非零场分量为\( E_r \frac{\eta I_0 l \cos\theta}{2\pi r^2}\left(1 \frac{1}{jkr}\right)e^{-jkr} \)\( E_\theta \frac{j\eta k I_0 l \sin\theta}{4\pi r}\left(1 \frac{1}{jkr} - \frac{1}{(kr)^2}\right)e^{-jkr} \)\( H_\phi \frac{jk I_0 l \sin\theta}{4\pi r}\left(1 \frac{1}{jkr}\right)e^{-jkr} \)其中 \( \eta \sqrt{\mu_0/\varepsilon_0} \approx 120\pi \approx 377\, \Omega \) 是自由空间波阻抗\( k 2\pi/\lambda \) 是波数\( \omega \) 是角频率\( E\)、\(H\) 之间满足 \(\mu H_\phi\) 和 \(E_\theta\) 的量纲匹配。我工作后回看这几条公式发现它们其实在讲三件事括号里的1对应远场辐射项随 \( 1/r \) 衰减第二项\( 1/(jkr) \)是感应近场项随 \( 1/r^2 \) 衰减第三项\( -1/(kr)^2 \)是静电场项只出现在 \( E_\theta \) 和 \( E_r \) 中随 \( 1/r^3 \) 衰减。离天线很近时第三项数值大得吓人但它只在天线周围转圈不会把能量真正送走离天线足够远后第一项占据绝对主导这一项才代表“真正的辐射”。我在仿真时就是通过分别绘制每一项的贡献直观看到了这三者随距离换班的画面。2.2 近场区和远场区的判据既然公式里有 \( kr \) 这个无量纲数那判据就从它来这一点非常关键。当 \( kr \ll 1 \)也就是 \( r \ll \lambda/(2\pi) \) 时场主要被高次项控制称为近场区无功近场区。这个区域里的电场和磁场相位差接近 \( 90°\)所以平均坡印廷矢量的实部很小能量大多储存在天线附近并不辐射出去。当 \( kr \gg 1 \)也就是 \( r \gg \lambda/(2\pi) \) 时1/r 项占主导\( E_\theta \) 和 \( H_\phi \) 同相比值趋近 \( \eta \)这就是远场区夫琅禾费区能量真正向外传播。中间那段过渡区域通常叫做辐射近场区菲涅尔区它没有严格的解析边界工程上一般按场景取经验值。我之前看到有同学直接把“近场”等同于“ \( r\lambda \) ”这其实是粗略做法正确的判据要看 \( r \ll \lambda/(2\pi) \)。 \( \lambda/(2\pi) \approx 0.159\lambda \)也就是大约六分之一波长。我一般仿真时把观察半径取到 \( 3\lambda \) 以上就能同时覆盖近场区和远场区在图上可以清晰看到场分量角色切换。下面这张表是我常用的快速参考区域名称判据主导衰减项能量性质无功近场区\( r \ll \lambda/(2\pi) \)\( 1/r^3 \)、\( 1/r^2 \) 项储能为主几乎不辐射辐射近场区\( r \approx \lambda/(2\pi) \) 到 \( r2D^2/\lambda \)\( 1/r^2 \) 项过渡到 \( 1/r \) 项储能与辐射并存远场区\( r \gg 2D^2/\lambda \)\( 1/r \) 项辐射为主波阻抗趋于 \( \eta \)对无穷小偶极子等效口径 \( D \) 很小\( 2D^2/\lambda \) 几乎可以忽略所以只要 \( r \gg \lambda/(2\pi) \) 就能当远场处理。我做这个题目时把观察点从 \( 0.1\lambda \) 扫到 \( 10\lambda \)画出的电场幅值曲线在约 \( 0.5\lambda \) 之后斜率就基本稳定到 \( -20dB/十倍距 \)这就是远场衰减特征。2.3 为什么方向图函数和辐射电阻也重要做偶极子仿真不能只盯着电场幅值方向图和辐射电阻是另外两个绕不开的结果。无穷小偶极子的归一化方向图函数是\( F(\theta) |\sin\theta| \)也就是说它在 \( \theta90°\)赤道面辐射最强沿着偶极子轴线\(\theta0°\) 和 \(\theta180°\)方向完全没有辐射。这个“面包圈”形状在三维图里特别明显是识别偶极子辐射最直观的特征。辐射电阻则反映天线把多少输入功率转化成了辐射功率。对无穷小偶极子\( R_r 80\pi^2 \left(\frac{l}{\lambda}\right)^2 \)如果 \( l \lambda/50 \)算出来只有 \( 80\pi^2 \times (1/50)^2 \approx 0.316\, \Omega \)。这个值非常小实际中很难匹配所以在工程上不会直接用这么短的偶极子发射功率这正是这个模型的“教学属性”所在。但模拟它时你需要把 \( R_r \) 也放到代码里作为检验仿真结果的一个重要坐标——如果仿真算出的辐射增益和理论值对不上就能很快发现问题。3. Matlab实现路线与核心代码拆解3.1 三条可走的路解析公式、矩量法、联合仿真拿到“Matlab模拟无穷小偶极子天线近场和远场”这个任务大多数人第一反应是直接套公式画图这也确实是最快能出结果的路径。但作为过来人我建议你至少先了解三条路第一完全用解析公式计算场分布。把上一节列出的 \( E_r \)、\( E_\theta \)、\( H_\phi \) 公式原样翻译成 Matlab 代码用ndgrid或meshgrid生成空间网格逐点计算复数场值再取模、相位、功率密度。优点是代码短、无需额外工具箱、逻辑透明适合验证理论缺点是你无法处理任意复杂结构只能算点源、偶极子这类理想模型。第二用 Matlab 自己实现矩量法MoM代码。把偶极子沿长度切分成若干小段每段假设为一个小偶极子建立阻抗矩阵并求解电流分布。这个过程会用到矩阵求逆复杂度高于解析解但能处理有限长度天线的端部效应和阻抗变化。如果不追求工程级精度完全可以在纯 Matlab 中跑通。第三用 Matlab 控制外部电磁仿真软件。比如通过 Matlab 脚本驱动 HFSS 或 CST 建模并取回结果再用 Matlab 做数据后处理和绘图。这条路适合你自己会的求解器已有限而还想复用 Matlab 生态的情况缺点是需要额外授权且学习曲线较陡。我在这次项目里选的是第一条路因为任务本身就是“模拟无穷小偶极子的近场和远场”解析公式就是最匹配的工具能直接在代码里体现出每一项的物理意义。后面如果想扩展到半波偶极子、阵列天线再考虑 MoM 也不迟。3.2 核心代码三维方向图与近场/远场可视化解析法实现的核心是一个双层循环或者网格化运算。我用的代码如下工程上叫“向量化”避免for循环逐点算clear; close all; clc; % 输入参数 lambda 1; % 工作波长单位m k 2*pi/lambda; % 波数 eta 120*pi; % 自由空间波阻抗 I0 1; % 电流幅度单位A l lambda/50; % 偶极子长度远小于波长 % 观察网格 N 80; % 网格密度 r linspace(0.1*lambda, 3*lambda, N); % 径向距离避开原点 theta linspace(0, pi, N); % 俯仰角0~pi phi linspace(0, 2*pi, N); % 方位角0~2pi [R, Theta] ndgrid(r, theta); % 距离-俯仰角网格 % 解析场计算 E_r eta*I0*l*cos(Theta)./(2*pi*R.^2) ... .* (1 1./(1j*k*R)) .* exp(-1j*k*R); E_theta 1j*eta*k*I0*l*sin(Theta)./(4*pi*R) ... .* (1 1./(1j*k*R) - 1./(k*R).^2) .* exp(-1j*k*R); H_phi 1j*k*I0*l*sin(Theta)./(4*pi*R) ... .* (1 1./(1j*k*R)) .* exp(-1j*k*R); % 取模用于绘图 E_r_mag abs(E_r); E_theta_mag abs(E_theta); H_phi_mag abs(H_phi);这段代码的核心是把课本公式一行行翻译过来没有用到任何信号处理或通信工具箱纯基础 Matlab 就能运行。注意1j表示虚数单位避免用i作为变量名导致复数单位被意外覆盖这是初学者常吃的暗亏。画三维方向图时我推荐在球坐标下进行因为方向图本身就是 \( F(\theta)\sin\theta \)天然适合球面展开% 三维辐射方向图 figure(Color, w, Position, [100 100 640 500]); [Phi_mesh, Theta_mesh] ndgrid(phi, theta); R_dir abs(sin(Theta_mesh)); % 归一化方向函数 % Matlabsph2cart坐标系中第二个参数是仰角需把theta换算成elev [X_dir, Y_dir, Z_dir] sph2cart(Phi_mesh, pi/2 - Theta_mesh, R_dir); surf(X_dir, Y_dir, Z_dir, R_dir, EdgeColor, none, FaceAlpha, 0.85); xlabel(x / \lambda); ylabel(y / \lambda); zlabel(z / \lambda); title(无穷小偶极子天线三维方向图 F(\theta)sin\theta); axis equal; colormap(jet); colorbar; view(135, 30);运行这段代码你会看到一个“面包圈”形状——赤道方向凸出来两极方向凹陷到零。这和我们前面说的 \( \theta90°\) 辐射最强完全吻合。近场和远场的可视化建议分两个图分别画不要挤在一起% E_theta近场幅值分布 figure(Color, w); [X_cart, Y_cart, Z_cart] sph2cart(Theta*0, pi/2-Theta, R); % 这里要先转二维 % 直接用二维切面展示更直观 [Phi_cut, R_cut] ndgrid([0, pi/2, pi, 3*pi/2], r); [X_cut, Y_cut, Z_cut] sph2cart(Phi_cut, 0, R_cut); % 计算E_theta模值在该切面上的分布 E_cut abs(1j*eta*k*I0*l*sin(pi/2)./(4*pi*R_cut) ... .* (1 1./(1j*k*R_cut) - 1./(k*R_cut).^2) .* exp(-1j*k*R_cut)); plot(R_cut, E_cut, LineWidth, 1.6); xlabel(距离 r / \lambda); ylabel(|E_\theta| (V/m)); title(E_\theta 幅值随距离变化赤道面); grid on; set(gca, YScale, log);这张图能很清楚地看到在 \( r0.5\lambda \) 时场幅值快速起伏并衰减过了远场区后则平滑地按 \( 1/r \) 规律下降。我通常还会叠加一条理论参考曲线 \( C/r \) 做对比方便判断是否进入了正确的远场区间。3.3 四种可视化方法的取舍近场和远场数据是复数怎么展示也是可以琢磨一下的可视化方式代码难度物理信息适用场景三维方向图曲面简单辐射强度空间分布远场方向图展示二维切面幅值曲线简单幅值随距离变化趋势近场/远场切换验证球坐标场强热力图中等近场幅度分布云图观察近场能量聚集相位等值面动画较复杂波前传播过程理解等相位面运动我做作业时最常用的组合是“三维方向图 切面曲线 热力图”因为三维图直观切面曲线能定量校验衰减趋势热力图能一眼看出近场的能量集中位置。相位动画虽然好看但计算量大而且对初学者来说不如幅值图容易解读。4. 实操过程与可复现的完整仿真流程4.1 参数初始化和网格设计任何仿真项目第一步都是把物理参数想清楚否则后面的代码就是空中楼阁。我把关键参数列成一张表写代码前先填好参数符号取值说明工作波长\( \lambda \)1 m归一化处理便于缩放波数\( k \)\( 2\pi \) rad/m由波长决定偶极子长度\( l \)\( \lambda/50 \)满足远小于波长假设电流幅值\( I_0 \)1 A归一化电流观察半径范围\( r \)0.1\(\lambda\) ~ 3\(\lambda\)覆盖近场和远场网格点数N80平衡精度和计算速度选 \( \lambda1 \) 是为了归一化后面的距离、场强全都可缩放这个技巧在做理论验证时非常实用。网格点数 80 在我的旧笔记本上大概零点几秒跑完如果调到 300 也能跑但三维图会明显变慢内存占用也高。网格设计的核心原则是在“关键变化区域”加密在平缓区域稀疏。近场区场变化剧烈我会在 \( 0.1\lambda \) 到 \( 1\lambda \) 之间用对数或分段线性方式加密远场区变化平缓均匀网格就够。4.2 分段实现代码完整项目建议拆成四个功能模块这样方便调试和复用。我当时的代码目录结构大概是dipole_simulation/ ├─ main.m % 主脚本 ├─ init_params.m % 参数初始化 ├─ compute_fields.m % 计算近场/远场电场磁场 ├─ plot_3d_pattern.m % 画三维方向图 └─ plot_fields_2d.m % 画切面幅值曲线和热力图main.m只负责调用各模块并控制整体流程这样每个函数都能单独测试。写功能模块时注意尽可能向量化Matlab 在数组运算上的性能远高于循环。如果要用循环也建议预分配数组例如E_field zeros(Nr, Ntheta);否则循环中动态增长数组会导致内存反复分配速度慢一个量级。关键模块compute_fields.m的核心代码在前面已经给出了这里补充一个判断近场/远场的辅助函数function region classify_region(r, lambda) Kr 2*pi*r/lambda; if Kr 1 region near; % 无功近场区 elseif Kr 3 region transition; % 过渡区 else region far; % 远场区 end end这里把 \( kr 1 \) 当作近场判据\( kr 3 \) 当作远场判据是有依据的\( kr1 \) 时 \( 1/(kr)^2 \) 项与 1 项等大此后高次项迅速衰减\( kr3 \) 时高次项已经小于 \( 11\% \)可以安全忽略。实际打印各区域的场值进行对比能帮你更直观理解“近场和远场的角色转换”。4.3 结果解读从等相位面到功率密度代码跑出来后不能只截个图就交差。我建议至少解读三个层面第一看幅值变化。在赤道面\( \theta90°\)取一条径向线画 \( |E_\theta| \) 随 \( r \) 的变化曲线对数坐标下观察斜率。在近场区曲线起伏明显因为 \( 1/r^2 \) 和 \( 1/r^3 \) 项剧烈相互影响进入远场区后曲线变成一条斜率为 \( -20dB/十倍距离 \) 的直线这就是 \( 1/r \) 衰减的典型特征。看到这条直线才说明仿真范围确实覆盖到了远场区。第二看相位。对 \( E_\theta \) 的相位做空间分布远场区的等相位面是同心球面相位值随 \( r \) 线性增加近场区的等相位面则明显偏离球面因为感应场分量的相位和辐射场不同。你可以通过angle(E_theta)提取相位并做三维切片那画面很有冲击力。第三看功率密度。算实功率密度P_avg 0.5 * real(E_theta .* conj(H_phi)); % 时间平均坡印廷矢量在近场区\( P_{avg} \) 会有明显振荡甚至负值但那不意味着能量回流只是因为在无功近场区储能项占主导局部坡印廷矢量方向会出现来回摆动到了远场区\( P_{avg} \) 变成平滑的正值且与 \( 1/r^2 \) 成正比此时能量稳定外流。我习惯把功率密度和幅值画在同一张图里两相对照物理过程非常清楚。5. 常见问题与排查技巧实录5.1 角坐标映射错位我在这个项目里踩过的第一个大坑是画出的三维方向图方向不对总觉得“面包圈”的轴线偏了。后来才发现是sph2cart的坐标系约定和物理学中常用的球坐标不一致。Matlab 的sph2cart(az, elev, r)中第二个参数是仰角范围在 \([-\pi/2, \pi/2]\)而天线理论里用的 \(\theta\) 是俯仰角范围在 \([0, \pi]\)两者之间差一个 \(\pi/2\)。转换时必须写成[X, Y, Z] sph2cart(phi, pi/2 - theta, r);否则 \( \theta0 \) 的位置会指向赤道方向而不是北极方向整个方向图就旋转了 90 度。这类坐标约定问题在所有三维绘图软件里都存在遇到图形方向不对时永远先检查坐标变换。5.2 原点奇点和数值溢出公式里带 \( 1/r^2 \) 和 \( 1/r^3 \) 项如果不做处理\( r0 \) 时运算直接产生Inf或NaN绘图就会缺一块。我的做法是让径向网格从 \( r0.1\lambda \) 开始从物理角度这也合理——偶极子本身有长度原点附近也不是理想点源。如果你需要观察更近的区域也要加一个很小的eps偏移比如 \( r 1e-6\lambda \)而不是直接取 0。另一个数值问题是复数运算的浮点误差。当 \( kr \) 很大时exp(-1j*k*r)的相位快速旋转数组里的实部和虚部可能出现大量相互抵消导致的精度损失。这个在偶极子模型里不太明显但如果计算阵列方向图时相位差很小就会暴露出来。建议用vpa或者其他高精度手段做局部验证至少心里有数。5.3 网格分辨率与内存平衡网格数从 80 提升到 500肉眼几乎看不出方向图差异但计算时间从零秒变成了几十秒三维surf也会卡顿。我后来摸索出的经验是方向图画 60~100 个点就够了曲线图可以细分到 500 个点不要全局加密。对想要观察近场细节的区域可以单独做一个局部高分辨率网格对远场均匀稀疏网格即可。如果机器内存有限尽量避免在同一个数组里保存全空间三维复数场可以分区域计算、分区域绘图。5.4 复数直接surf导致的空白图surf函数默认不接受复数输入如果你直接把E_theta传给surf它要么报错要么只画出实部虚部被静默忽略。更麻烦的是plot对复数会默认按实部为横轴、虚部为纵轴绘制出来的图看起来像一条螺旋线和你想的完全不是一回事。所以一切绘图前先取abs()、angle()或real()不要直接拿复数矩阵去画。这是一个特别低级的错误但我在给学弟学妹看代码时经常看到类似写法。5.5 近场为什么画不出“波”的形状有同学跟我说他想画近场的“波状传播”但画出来的图一点波动都没有。问题在于他把 \( kr0.5 \) 到 \( kr500 \) 的全部范围放在线性坐标里导致近场区域的快速变化被远场的缓慢变化压扁根本看不出波纹。解决方法是改用对数坐标显示或者分段裁剪先单独画 \( r \in [0.1\lambda, 2\lambda] \)再画 \( r \in [2\lambda, 100\lambda] \)。另外如果你想看到清晰的驻波或行波特征需要同时画电场实部和虚部只看幅值是看不出相位传播的。我把上面这些问题整理成一个排查表方便你按图索骥症状可能原因检查方法方向图轴线偏转sph2cart仰角/俯仰角混用检查pi/2 - theta原点附近有空白/Inf\( r0 \) 未做偏移将 \( r \) 起点设为 0.1\(\lambda\)曲线忽大忽小难以观察线性坐标压缩近场细节尝试YScale,logsurf报错或画不出复数直接传入先abs()再绘图远场衰减斜率不对观察范围未覆盖远场增大 \( r \) 至 \( 3\lambda\) 以上相位动画不连贯相位跳变 2\(\pi\)使用unwrap处理6. 从偶极子到实用天线下一步该怎么走6.1 半波偶极子和无方向性天线对比跑通了无穷小偶极子你手里其实已经握住了一套可以向上扩展的方法论。最自然的下一步是半波偶极子长度 \( l\lambda/2 \)它的电流不再均匀分布而是近似余弦分布方向图函数变成\( F(\theta) \frac{\cos(\frac{\pi}{2}\cos\theta)}{\sin\theta} \)和 \( \sin\theta \) 相比半波偶极子的方向图在赤道面更扁平方向性系数从 1.5 提升到 1.64辐射电阻也从不到 1 欧姆跳到约 73 欧姆更接近实际工程中的 50 欧姆或 75 欧姆系统。这个对比非常直观地说明了“为什么工程上不用无穷小偶极子而要用半波偶极子”。代码层面你只需要改动方向图函数和场表达式中的电流分布项其余网格、绘图、后处理逻辑全部复用。这也是为什么我强调要把项目拆成模块扩展性会好很多。6.2 Matlab与全波仿真软件的协作Matlab 在理论验证和算法原型上很强但处理复杂结构时工业界更依赖 HFSS、CST 这类全波仿真工具。那么两者怎么协作我见过一种非常高效的用法用 HFSS 建模复杂天线结构导出远场方向图数据和近场场分布数据再用 Matlab 做深度后处理——画方向图、算波束形成权值、做优化算法迭代。HFSS 自带脚本功能可以批量导出结果到文本或.mat文件Matlab 直接读取后处理。对正在找工作或者搞科研的人掌握这条数据链路是很加分的。如果你不想装全波仿真软件也可以直接在 Matlab 里写矩量法代码。对偶极子天线矩量法实现不复杂把偶极子分 N 段每段上的电流用基函数展开构建阻抗矩阵 \( Z_{mn} \)求解电流系数再叠加所有小段的辐射场。这一步做完你对“数值解”和“解析解”的关系理解会再上一个台阶。6.3 我的工具链建议和心得做这类仿真项目我的个人建议是先把解析模型玩透再碰通用电磁仿真软件。理由很简单你只有知道了正确答案长什么样才能判断软件给出的结果到底可不可信。我用 HFSS、CST 都建过很多天线模型经常遇到网格剖分不合理、端口激励设置错误等问题导致结果严重偏离理论值如果心里没有解析解做锚点你根本无法分辨这是真实物理现象还是仿真错误。还有一条是关于版本选择的Matlab 自身版本对这类小项目几乎没有影响R2020a 和 R2024a 跑同样代码都一样快。真正影响效率的是你写代码的方式比如用向量化而不是循环、用模块化而不是单文件堆到底、用clear和clc保持工作区干净。这些习惯比版本本身值钱得多。结尾回到这个项目本身我认为它最大的价值不是画出一张漂亮图而是逼着你去理解“近场和远场到底是怎么分界的”。我在做这个仿真之前对 \( kr \ll 1 \) 和 \( kr \gg 1 \) 只有一个模糊的印象真到代码里设定判据、观察曲线斜率变化时才意识到这个划分背后是电磁场能量从“储能”到“辐射”的物理切换。所以如果你也在做类似的验证千万别只抄代码动手改一改参数比如把 \( \lambda \) 从 1 改成 10把网格密度调大调小把观察距离从 1\(\lambda\) 拉到 10\(\lambda\)你会明显感觉到自己对场分布的理解会发生变化。最后再分享一个小技巧保存仿真图时建议把三维方向图和近场/远场切面图分开保存成矢量图格式这样论文里放大后不会糊。Matlab 里exportgraphics(gca, dipole_pattern.pdf, ContentType, vector)一行就能搞定。这个细节在写期末报告或发小论文时很加分我当年就是吃了截图模糊的亏才养成了这个习惯。本文还有配套的精品资源点击获取