1.2 物理学基础与数学工具
1.2 物理学基础与数学工具
本节摘要：SOURCE 1.2 将 Maxwell 方程、傅里叶变换与近远场判据翻译为 PCB 可操作的铜箔间距、叠层与带宽估算。
故障场景：100 ps 边沿下「集总参数」仿真与实测不符
工程师用 SPICE 集总模型预测 3.3 V CMOS 输出，实测却出现 GHz 辐射峰。SOURCE 解释：t_r=100\,\text{ps} 时有效带宽 f_{eff}\approx 0.35/t_r\approx 3.5\,\text{GHz}，位移电流项不可忽略，须用分布参数/传输线模型。
同样的错误还常发生在去耦电容建模上：把 0402 电容当成纯容性元件，忽略其 ESL≈0.4 nH 的安装电感。在 3.5 GHz 处，单颗 100 nF 电容的阻抗已由 1/(\omega C)=0.45\,\Omega 反转为 \omega L\approx 8.8\,\Omega——高于目标阻抗，滤波失效。本节的所有工具，都是为了避免「量级错判」这类错误。
一、Maxwell 在 PCB 上的降维
| 方程项 | PCB 含义 | 设计映射 |
|---|---|---|
| \nabla\cdot\mathbf{D}=\rho_v | 电荷守恒 | PDN 瞬时供电荷能力 |
| \nabla\cdot\mathbf{B}=0 | 无磁单极 | 地平面完整性 |
| \nabla\times\mathbf{E}=-\partial\mathbf{B}/\partial t | 变化磁场生电场 | 串扰、EMI 环路面积 |
| \nabla\times\mathbf{H}=\mathbf{J}+\partial\mathbf{D}/\partial t | 电流生磁场 | 高频走线、过孔电感 |
50 Ω 微带线的 Z_0=\sqrt{L'/C'} 正是 Maxwell 在铜箔-FR4-参考平面边界下的本征解。工程上最常用的是 IPC-2141 微带线近似公式：
其中 H 为介质厚度（mil）、W 为线宽、T 为铜厚、\varepsilon_r 为介质介电常数。常用 FR4 的 \varepsilon_r 取 4.2–4.6，损耗正切 D_f 约 0.015–0.020；高频板材如 Megtron 6 的 D_f 可低至 0.002，适合 25 Gbps 以上通道。
二、傅里叶：从边沿到频谱风暴
真实数字信号的有效谐波带宽：
PCIe 6.0 若 t_r<15\,\text{ps}，频谱延伸至 30 GHz 以上，与 CISPR/FCC 辐射频段高度重叠。阻抗不连续引起的振铃会在频域引入间隔 1/(2\tau) 的离散谱线——同一芯片在不同走线长度下可能在完全不同频率超标。
典型的带宽估算可以直接用脚本完成，便于在设计评审时快速给出「这条信号会辐射到多高频率」的结论：
# 信号带宽与谐波估算（阶跃边沿模型） def bandwidth(tr_ps: float) -> float: """由 10%-90% 上升时间估算有效带宽，单位 GHz""" return 0.35 / (tr_ps * 1e-3) # tr_ps 单位 ps，换算为 ns # 常见接口的边沿与带宽 interfaces = [ ("PCIe Gen3", 35, 10.0), ("PCIe Gen4", 25, 14.0), ("DDR4-3200", 120, 2.9), ("USB 3.2", 100, 3.5), ("1000BASE-T", 300, 1.17), ("SATA 6G", 90, 3.9), ] print(f"{'接口':<12}{'tr/ps':<8}{'BW/GHz':<8}") for name, tr, _ in interfaces: bw = bandwidth(tr) print(f"{name:<12}{tr:<8}{bw:<8.2f}") # 三次谐波是否进入 FM 广播频段（88-108 MHz）？ # 例：48 MHz 时钟，tr=1 ns → BW=0.35 GHz，三次谐波 144 MHz 在 FM 之上 # 例：12 MHz 时钟，tr=2 ns → 三次谐波 36 MHz，接近 30 MHz 起始辐射限值频段 print("48 MHz 时钟三次谐波:", 48 * 3, "MHz")
这段脚本的价值不在于计算本身，而在于把「边沿—带宽—谐波—标准频段」的链条显式化。任何一个 48 MHz 晶振配 1 ns 边沿的 IO，三次谐波 144 MHz 都会落在 CISPR 32 B 类 30 MHz–1 GHz 考核频段内，必须由布局质量而不是运气来保证达标。
三、近场 vs 远场
近场以感应为主，远场以辐射为主；判据与波长 \lambda 相关。场边界经验判据为距离 d=\lambda/(2\pi)：d 之内是近场（电偶极近场 E\propto 1/r^3，磁偶极近场 H\propto 1/r^3），之外是远场（E,H\propto 1/r，波阻抗恒为 377 Ω）。
| 频点 | 波长 λ | 近/远场边界 λ/2π | 诊断要点 |
|---|---|---|---|
| 30 MHz | 10 m | 1.59 m | 大多数 PCB 尺寸均在近场 |
| 144 MHz | 2.08 m | 33 cm | 近场探头可精确定位 |
| 1 GHz | 30 cm | 4.8 cm | 近场区很窄，须精细扫描 |
| 5 GHz | 6 cm | 9.5 mm | 只能逐段扫，注意探头姿态 |
PCB 调试多用近场探头定位 dipole-like / loop-like / slot-like 辐射源类型：
- dipole-like（电场探头测到强 E 场）：源头是共模电压，常见于跨分割走线与长线缆；
- loop-like（磁场探头测到强 H 场）：源头是环路电流，常见于高速信号—回流环路；
- slot-like（沿缝隙两侧扫描出现峰）：源头是开槽/缺口，常见于地平面槽与屏蔽罩缝隙。
四、微带线阻抗计算示例
微带线是 PCB 上最常见的高速传输结构，其特性阻抗由线宽、介质厚度与介电常数共同决定。设计阶段的第一步是「由目标阻抗反推线宽」，这一步可以用近似公式快速完成，再由板厂按精确场解微调：
import math def microstrip_z0(w_mil: float, h_mil: float, t_mil: float, er: float) -> float: """IPC-2141 微带线阻抗近似，单位 mil，返回 Ω""" return (87 / math.sqrt(er + 1.41)) * math.log(5.98 * h_mil / (0.8 * w_mil + t_mil)) # 6 层板典型参数：H=4.0 mil（半固化片）、T=1.2 mil（0.5 oz 铜厚） # 目标阻抗 50 Ω，FR4 εr≈4.3 h, t, er = 4.0, 1.2, 4.3 for w in [5.0, 6.0, 7.0, 8.0, 9.0]: z = microstrip_z0(w, h, t, er) print(f"线宽 {w:.1f} mil → Z0 = {z:6.1f} Ω") # 结论示例：H=4.0 mil 时约 7.5 mil 线宽可得 50 Ω # 若换 1.6 oz 铜厚（T=2.4 mil），同样线宽阻抗略降，须重新核算
这个计算说明了一个关键点：阻抗不是单独决定的，而是「线宽—介质厚度—铜厚」三元平衡。改叠层厚度而不改线宽，目标阻抗就会漂移；同样，同一线宽在两层不同介质厚度的参考平面之间会得到不同 Z0。这也是为什么层叠结构一旦确定，就应冻结信号层与参考层的间距参数。
⚠️ 常见坑：低频用集总、高频仍用同一模型——边沿变快时必须切换分布参数思维。
💡 关键直觉：所有「噪声」都是 Maxwell 在非理想边界下的强迫响应。
一节小结
- f_{eff}\approx 0.35/t_r 决定与 EMC 标准的频谱重叠；
- 地平面分割 = 强迫回流绕行 → 共模辐射；
- 工具链：Maxwell 定边界，FFT 定频段，近场定位置；
- 微带线 Z0 公式把板材参数（H、W、T、εr）直接映射为线宽决策；
- 近/远场边界 λ/2π 决定探针选型与扫描策略。