用MATLAB计算矩形微带天线尺寸：从公式到仿真实战-育师

简介：天线设计入门常从微带天线开始，其尺寸计算是射频工程设计的基础环节。矩形微带贴片天线的物理尺寸依赖中心频率、基板介电常数与厚度三要素，需通过有效介电常数和边缘场修正准确求得。利用MATLAB脚本可实现宽度、长度及馈点位置的快速估算，Hammerstad经验公式是其中的关键补偿模型。这一方法适用于2.4GHz、5.8GHz等频段的PCB天线初步设计，为后续全波仿真与阻抗匹配优化提供可靠初值。以完整可运行的MATLAB代码为例，拆解尺寸推导顺序、单位坑与参数边界，可帮助工程师缩短天线研发迭代周期。

1. 矩形微带天线的尺寸计算：一段 MATLAB 代码背后的完整设计逻辑

矩形微带天线的尺寸计算，看起来只是把频率、介电常数、板厚三四个数字塞进公式，按回车就能得到贴片的 W 和 L。实际上这套几行 MATLAB 代码的工程价值在于：它用哈默斯塔德（Hammerstad）经验公式补偿了边缘场效应，把宽度、有效介电常数、长度修正量和馈点位置按正确顺序一次算全。顺序颠倒、单位混用、忘了回代，最后算出来的尺寸只是“看起来正确”，进仿真和实测立刻露馅。这份资源适合刚接触微带天线、手头有明确板材和频段要求的设计者，也适合要快速估算贴片尺寸再进全波仿真微调的工程师。注意，代码里 f0 单位是 Hz、h 单位是 cm、c 单位是 m/s，单位并不统一，这是最容易翻车的起点。

2. 从三个基础参数推出六个关键尺寸：公式拆解与选用理由

2.1 最少输入参数：f0、er、h，为什么说它们是“设计三要素”

矩形微带天线的几何尺寸主要由三个参数决定：中心频率 f0、基板相对介电常数 er、基板厚度 h。这三个量也是选板材时最先确定的三个数。f0 由系统指标决定，比如 RFID 读写器常用的 2.45GHz、5.8GHz，或者 Wi-Fi 的 2.4GHz/5.2GHz；er 由板材决定，常见的有 FR4 的 4.4、罗杰斯类板材的 3.38 或 2.2；h 由板材厚度决定，常见的是 0.8mm、1.0mm、1.6mm。三者缺一不可，其中 er 和 h 直接进入有效介电常数与边缘场补偿项，对最终长度 L 的影响非常明显。

用这段 MATLAB 代码时，用户只需要改三个赋值语句：f0、er、h。其他变量全部由公式自动推导。这个设计思路是“最小输入”原则——需要用户决定的物理量尽量少，避免在中间步骤引入人为误差。光速 c 可视为常量，但代码里保留为变量，方便想用英寸或毫米制的人自行换算。

2.2 先宽后长：W、eeff、al、L 四步推导的物理依据

公式顺序不是随手排的，它遵循微带贴片天线的传输线模型推导逻辑。第一步先算宽度 W：

W = c / (2 * f0 * sqrt((er + 1) / 2));

这里 (er+1)/2 是“平均介电常数”的近似，适用于宽度远大于厚度、且基板不太厚的情况。W 的物理意义是：贴片等效为一条宽 W 的微带传输线，宽度决定了辐射边的长度，也直接影响辐射电阻和阻抗带宽。W 先算出来，是因为后面的有效介电常数和长度修正都依赖 W/h 这个比值。

第二步算有效介电常数 eeff：

eeff = ((er + 1) / 2) + (er - 1) / (2 * sqrt(1 + 12 * (h / W)));

电场线在贴片边缘会有一部分穿过空气，所以实际的“平均介电常数”介于 1 和 er 之间，并且随 W/h 变化。这个公式是微带线理论中的常用近似式，适用于 W/h 大于 1 的常规薄基板。它比直接用 er 更接近电磁波在贴片下方真实传播时的相速度。代码里 h 和 W 在同一公式中比较，这里就埋了单位坑：W 由 c 算出来单位是米，如果 h 用厘米填，比值直接错掉，后面的修正量全部跟着偏。

第三步算边缘场延长量 al：

al = (0.412 * h * (eeff + 0.3) * ((W / h) + 0.264)) / ((eeff - 0.258) * ((W / h) + 0.8));

这是哈默斯塔德给出的经验公式，用于计算贴片两个辐射边因为边缘场“看起来多出来”的长度。电磁波在贴片开口端不会突然截止，场会向外扩展一小段距离，等效成电长度变长。如果不补偿这个量，天线的谐振频率会比设计值偏低，仿真和实测都会出现中心频率偏移。

第四步算贴片长度 L：

L = (c / (2 * f0 * sqrt(eeff))) - 2 * al;

贴片的电长度近似等于介质中半个波长，几何长度要减去两个边缘延长量（左右各一个）。先用 eeff 算出理想半波长长度，再减掉 2 倍的 al，得到最终的物理长度。这里最容易犯的错是只减一个 al，导致 L 偏长，谐振频率偏低。

2.3 馈点位置 Xf：为什么用 tre 而不是直接用 eeff

同轴线馈电时，馈点位置决定了输入阻抗。贴片边沿阻抗高，中心阻抗低，从边沿往中心移动的过程中存在一个输入阻抗等于 50Ω 的位置。代码里的 Xf 就是沿 L 方向从边沿往中心算出的馈点距离：

tre = ((er + 1) / 2) + (er - 1) / (2 * sqrt(1 + (12 * (h / L)))); Xf = L / (2 * sqrt(tre));

注意 tre 的公式和 eeff 很像，但分母里的特征尺寸从 W 换成了 L。这是因为 L 方向的电场分布与 W 方向不同，计算 Xf 时应该使用 L 方向对应的有效介电常数。文档里把这个变量命名为 tre，实际它就是 L 方向上的有效介电常数近似。用 eeff 代替 tre 算出的 Xf 偏差可能在 0.1mm 量级，虽然不大，但对于追求低驻波比的馈点设计来说，0.1mm 的偏差足以让 S11 最深点移动十几兆赫兹。

3. 坐到 MATLAB 前复现：完整脚本、运行结果与逐行参数修改

3.1 一份可以直接粘贴运行的改进脚本

原文档的写法是“打开命令窗口，把代码全部复制粘贴进去”，对于一次性计算够用，但不方便改参数和排查单位问题。我一般会把它整理成带注释的脚本文件，并把所有单位统一到米制。下面这份脚本保留了原文档的全部公式，只把 h 的赋值和注释做了修正：

% 矩形微带天线尺寸估算脚本 % 所有长度单位统一为 米（m） % 修改下面三个参数即可，其他行勿动 f0 = 2.45e9; % 中心频率，单位 Hz er = 3.38; % 基板相对介电常数 h = 0.0016; % 基板厚度，单位 m（1.6mm 板材） c = 3e8; % 光速，单位 m/s % Step 1: 贴片宽度 W = c / (2 * f0 * sqrt((er + 1) / 2)); % Step 2: 有效介电常数（W 方向） eeff = ((er + 1) / 2) + (er - 1) / (2 * sqrt(1 + 12 * (h / W))); % Step 3: 边缘场延长量（Hammerstad 经验公式） al = (0.412 * h * (eeff + 0.3) * ((W / h) + 0.264)) / ((eeff - 0.258) * ((W / h) + 0.8)); % Step 4: 贴片长度 L = (c / (2 * f0 * sqrt(eeff))) - 2 * al; % Step 5: L 方向有效介电常数 tre = ((er + 1) / 2) + (er - 1) / (2 * sqrt(1 + (12 * (h / L)))); % Step 6: 馈点距辐射边沿的距离 Xf = L / (2 * sqrt(tre)); % 显示结果，单位均为米 disp(['W = ', num2str(W), ' m']); disp(['L = ', num2str(L), ' m']); disp(['Xf = ', num2str(Xf), ' m']);

这段脚本把原文档的“手输命令”变成了可重复执行的文件。跟原版相比，改动点只有一个：h 从 0.005 改成 0.0016，并且在注释里写明是 1.6mm 板材。为什么这么改？因为原文档说明 h 的单位是 cm，那么 0.005cm 是 0.05mm，这相当于超薄柔性基板，和大多数刚性 PCB 板材的厚度不在一个量级。如果你手头是 1.6mm 的板子，h 应该填 0.0016。如果手头是 0.8mm 板，h 填 0.0008。别直接抄原文档的 0.005。

3.2 运行结果与数值核对方法

以 f0=2.45GHz、er=3.38、h=1.6mm 为例，脚本输出大约为：

W = 0.0414 m L = 0.0328 m Xf = 0.0093 m

也就是 W 约 41.4mm，L 约 32.8mm，Xf 约 9.3mm。这几个数值是有物理意义的：2.45GHz 自由空间波长约 122mm，贴片宽度 41mm 接近半波长但不到，长度 33mm 因为介电常数而被压缩，这些都是正常量级。核对时可以记住一个经验规则：矩形微带贴片的两个边长通常在 0.25 到 0.5 个自由空间波长之间。如果算出来的尺寸明显超出这个范围，先回头查单位，尤其是 h 和 W 的比值是否合理。

3.3 参数修改的三种典型情况

把这个脚本当作模板，最常见的三种改法是换频段、换板材、换板厚。

修改场景	f0	er	h	预期效果
2.4GHz 换 2.45GHz	2.4e9 → 2.45e9	不变	不变	W、L 略变小
FR4 板（er=4.4）	不变	3.38 → 4.4	不变	W、L 明显变小
1.6mm 换 0.8mm 板	不变	不变	0.0016 → 0.0008	L 变化不大，Xf 略变

换频段时，f0 直接出现在 W 和 L 的分母上，频率越高尺寸越小，这是线性关系。换介电常数时，er 变大导致相速度变慢、波长变短，W 和 L 都变小。换板厚时，h 主要影响 eeff 和 al，对 L 的影响在毫米级，但对馈点位置和阻抗带宽影响不可忽略。改完参数后，建议先把输出结果按毫米读一遍，再决定是直接打样还是先进仿真。

4. 把公式用到自己的频段和板材上：参数边界与走向判断

4.1 趋势对照：f0、er、h 对尺寸的影响

为了在改参数时不盲目，我通常会在算完一组值之后，再扫两个相邻频点或换一种板材做对照，确认尺寸变化方向符合预期。下面是我用同一脚本对几种常见组合做的估算，单位统一为毫米。

参数组合	W（mm）	L（mm）	Xf（mm）
2.45GHz，er=3.38，h=1.6mm	41.4	32.8	9.3
2.40GHz，er=3.38，h=1.6mm	42.2	33.6	9.4
2.50GHz，er=3.38，h=1.6mm	40.6	32.1	9.1
2.45GHz，er=4.4，h=1.6mm	37.3	28.8	8.1

从这个表能看出两个规律。第一，频率从 2.4GHz 升到 2.5GHz，W 和 L 都收缩了 1mm 左右，方向一致。第二，er 从 3.38 换成 4.4，L 一下子收缩了 4mm，比单纯调频率的影响大得多。这说明板材选型比频率微调更敏感。如果仿真发现谐振点偏低，常见的修正手段是微调 L，而如果尺寸严重不对，先怀疑 er 填错了。

4.2 公式的适用边界：薄基板假设与 W/h 比值

这套公式不是万能的，它基于薄基板近似，使用前必须检查两个边界条件。第一是基板电厚度：一般要求 h 远小于工作波长，工程上通常按 h 小于 0.02λ0 到 0.05λ0 来把握。2.45GHz 时空气波长约 122mm，1.6mm 约为 0.013λ0，满足条件；但如果做到 5.8GHz 且基板还是 1.6mm，h 达到 0.031λ0，公式误差会明显增大，这时需要全波仿真修正。第二是 W/h 比值：哈默斯塔德公式隐含 W/h 大于 1 的适用前提，本例 W/h 约 26，远离边界。如果设计的是窄贴片（W/h 接近甚至小于 1），有效介电常数和边缘延长量的精度都会下降，建议换用更严格的计算模型。

另外要注意，这套公式是一阶估算，没有计入地面尺寸、馈电探针引入的电感、贴片拐角的高次模效应。它给出的尺寸适合作为仿真建模的起点，不承诺直接加工就能得到完美驻波。我一般把算出来的尺寸当作“初值”，紧接着做两件事：一是进全波仿真验证谐振频率，二是对 Xf 做参数扫描找最优匹配点。

4.3 从算尺寸到能加工：圆整、公差和留余量

算出来的 W、L、Xf 带着一串小数，不能直接发给板厂。我的习惯是统一圆整到 0.1mm。频率越高，尺寸公差越敏感，2.45GHz 下 0.3mm 的长度误差可能让谐振点移动几兆赫兹到十几兆赫兹，所以不要做 0.01mm 级别的圆整——那是骗自己，板材蚀刻精度根本到不了。圆整后还要给馈点位置留调整空间。常见做法是：先按 Xf 打一个过孔位置，同时预留 1mm 左右的余量，方便后续用仿真扫参决定是否移动。如果条件允许，可以在贴片 L 方向两端各留一点铜皮，实测频率偏高（天线偏短）再往里补，频率偏低再往外锉，这是后文会提到的“后悔药”技巧。

5. 避坑手册：从 MATLAB 报错到中心频率偏移的五个实战记录

5.1 现象：同一块 1.6mm 板，算出的 L 从 30mm 到 33mm 来回跳

原因：h 的单位没统一。原文档写明 h 单位为 cm，示例值 0.005 如果按 cm 理解是 0.05mm，但实际常用板材是 1.6mm，即 0.16cm。如果把 h 直接填成 0.005，等于把 1.6mm 的板子当作 0.05mm 超薄基板算，al 和 eeff 全部失真，最终 L 偏差可达 2mm 以上。如果再把 h 按 m 填成 0.0016 但把 W 按 cm 读，又是一套错数值。

解决：所有长度统一用米。h 直接按米填：1.6mm 填 0.0016，0.8mm 填 0.0008。W、L、Xf 的输出也全部按米读，最后再转毫米。不要在一个公式里混用 cm 和 m，这是最稳妥的。

5.2 现象：Xf 算出来和仿真最优位置差了将近 1mm

原因：计算 Xf 时用了 eeff 而不是 tre。eeff 是用 W 算出来的，对应的场分布方向是宽度方向；Xf 是在长度方向上移动馈点，应该用 tre，也就是用 L 回代的有效介电常数。两者差异虽然只有零点几，但叠加在谐振长度上，馈点阻抗匹配点就会漂移。

解决：严格按第 3.1 节的顺序跑，先算 L，再算 tre，最后算 Xf。不要图省事直接复制 eeff 进 Xf 的公式。

5.3 现象：尺寸都对，同轴线馈电仿真 S11 还是不到 -10dB

原因：Xf 只给出了 L 方向的馈点位置，矩形贴片还有 W 方向。实际设计中馈点可能不在贴片中心线上，而是偏移到某个位置以兼顾两个方向的阻抗。文档代码没有给出 W 方向馈点位置的计算，这部分要用同样的思路补算，或者直接用仿真扫参。

解决：建模时把同轴探针放在 Xf 位置先行验证，如果阻抗偏大或偏小，沿 W 方向微调探针坐标，每次移动 0.5mm 扫一遍 S11，找到最佳点。不要指望公式给的位置一次到位。

5.4 现象：把代码粘进命令窗口直接报错 Undefined function or variable

原因：最常见的是括号或引号被自动替换成全角字符，尤其是从网页或文档复制时。其次是命令窗口里残留了上一次的变量，导致中间变量被旧值覆盖。MATLAB 对全角符号基本不识别，报错信息会指向某一行，但根源往往是肉眼看不到的字符差异。

解决：把代码存成 .m 文件再运行，不要直接粘贴到命令窗口。如果一定要粘贴，先执行 clear all，然后确认所有括号、引号都是半角字符。脚本文件方式也方便下次改参数。

5.5 现象：按算好的尺寸打样，实测中心频率偏了 80MHz

原因：公式是一阶估算，介电常数在不同频率下有色散，板材批次也有偏差，蚀刻精度和过孔工艺都会改变等效电气长度。频率偏低说明天线电气长度偏长，频率偏高说明电气长度偏短。

解决：先仿真微调再打样。实测后用“微调 L”的方式补救：频率偏低，把贴片沿长度方向两端各锉短一点，每次 0.2mm 量级，配合网络分析仪看 S11 最深点移动方向。频率偏高则只能重新做板或在两端补铜箔，这也是为什么我建议每一版都留出 L 方向的修整余量。

6. 算完尺寸先别打样：一轮阻抗扫描把“图纸尺寸”变成“可交付尺寸”

拿到 W、L、Xf 初值后，我给自己定了一条规矩：至少做一轮仿真扫参再下单，这个习惯帮我省掉了不止一次返工。

第一步，按毫米值建仿真模型。W 取 41.4mm、L 取 32.8mm，地面比贴片每边外扩至少 5 到 10mm，同轴探针半径按 0.65mm 左右建模，介质基板高度 1.6mm。第二步，扫 Xf。固定 W 和 L，把探针从 8mm 扫到 10.5mm，步进 0.5mm，记录 S11 的 -10dB 带宽和阻抗圆图轨迹。