S参数分析 — 信号品质(SI)评估的散射矩阵理论与实践

S参数的理论基础

什么是S参数

比喻一下，向隧道喊话时"反响回来的声音"就是S11，"传到对面的声音"就是S21。

具体来说，PCIe Gen5（32 GT/s）的奈奎斯特频率是16 GHz，在基板上的波长约6 mm。10 cm的配线就包含十多个波长，不能用集中参数电路来描述。

散射矩阵的定义

• $S_{11}$ — 端口1的输入反射系数。端口2由基准阻抗 $Z_0$ 终接时，入射到端口1的波中有多少被反射
• $S_{21}$ — 端口1→端口2的正向透射系数。信号有多少被传输
• $S_{12}$ — 端口2→端口1的反向透射系数（无源元件中 $S_{12} = S_{21}$）
• $S_{22}$ — 端口2的输出反射系数

入射波和反射波定义为归一化功率波：

回波损耗和插入损耗

• 回波损耗 > 10 dB（|S11| < 0.316）→ 反射功率在10%以下。最低要求
• 回波损耗 > 15 dB（|S11| < 0.178）→ 反射功率在3%以下。高速传输路的一般目标
• 回波损耗 > 20 dB（|S11| < 0.1）→ 优秀的匹配。连接器单体的目标

插入损耗随频率增大而增加。例如10 cm的FR-4微带线，在1 GHz约0.5 dB，在10 GHz约3 dB，在20 GHz约7 dB。低损耗基板（如Megtron6）同样长度可降至一半。

无源性和因果律

1. 无源性（接收性）条件 — 网络不产生能量：

• 时间域脉冲响应 $s(t)$ 在 $t < 0$ 时为零
• 频域中，实部和虚部满足Hilbert变换对（Kramers-Kronig关系）

对策是：

• 使用内置无源性强制（Passivity Enforcement）算法的工具（Ansys SIwave、Keysight ADS、Cadence Sigrity等）
• 因果律破坏多数是由低频数据缺失或外推不当引起的。正确进行DC点和高频外推

混合模式S参数

• $S_{dd21}$ — 差动插入损耗（SI设计者最关注的量）
• $S_{dd11}$ — 差动回波损耗
• $S_{cd21}$ — 差动→共模转换（歪斜和非对称性的指标）
• $S_{cc21}$ — 共模插入损耗（影响EMI特性）

例如USB4和PCIe Gen6规范为 $S_{dd21}$ 和 $S_{dd11}$ 定义了掩码（限界值曲线），不满足这些掩码的通道就不符合规范。

S参数的数值计算手法

S参数的获取手法

• 3D全波EM分析（FEM / FDTD / MoM） — 计算连接器、过孔、封装等3D结构的S参数。Ansys HFSS（FEM）、CST（FDTD/FIT）代表。精度最高但计算成本大
• 2.5D平面分析 — 将PCB配线图形作为2.5D结构快速求解。Cadence Clarity、Ansys SIwave、Keysight Momentum等。最适合多层基板S参数提取
• VNA（矢量网络分析仪）实测 — 直接测量实物S参数。校准（SOLT/TRL）精度决定了结果质量

1. 用四面体网格离散化求解域。电磁场用边界单元（Nedelec单元）插值
2. 求解Maxwell方程的弱形式（变分形式），计算各频率下的E场和H场分布
3. 从端口面的模式场与求解结果的内积中，提取入射波振幅 $a_i$ 和反射波振幅 $b_i$
4. 由 $S_{ij}(f) = b_i / a_j$ 构成S参数（端口j除外其他端口由匹配终端）

简单说，就是求解Maxwell方程得到的结果，再从端口面提取波动振幅。

频域与时域的转换

• DC点（0 Hz）数据要正确外推。缺失会导致IFFT结果有直流偏移
• 高频侧外推 — 最好数据的最高频率是奈奎斯特频率的3~5倍
• 窗函数应用 — 用Gauss或Kaiser窗抑制Gibbs现象（振铃）
• 等间隔频率采样 — IFFT需要等间隔数据。非等间隔时要插值

De-embedding技术

代表性的方法：

• TRL校准（Thru-Reflect-Line） — 从3个标准结构进行校准，移动校准面。是PCB级测量的最可靠方法
• AFR（Automatic Fixture Removal） — 仅从Thru结构推估固定装置特性的算法
• 2xThru法 — 从两倍长的固定装置结构，用时域门选分离DUT

高速SerDes的规范验证（CEI-112G等）中，测试板的De-embedding精度直接影响是否符合规范。De-embedding错了就会造成"规范不符"的误判。

S参数的实务应用

Touchstone文件的运用

• .s1p — 1端口（例：天线的输入反射特性）
• .s2p — 2端口（例：单端传输路）
• .s4p — 4端口（例：一对差动信号）
• .s8p — 8端口（例：一对差动对加上串扰）
• .s12p以上 — 多通道SerDes（连接器全端口等）

# GHz S MA R 50
! 频率[GHz]  |S11| ∠S11  |S21| ∠S21  |S12| ∠S12  |S22| ∠S22
1.000  0.05  -12.3  0.98  -5.6  0.98  -5.6  0.04  -15.1
2.000  0.08  -25.1  0.95  -11.2  0.95  -11.2  0.07  -28.3
...

其中"S"表示S参数，"MA"表示幅度-相位（Magnitude-Angle）形式，"R 50"表示基准阻抗50Ω。还可以用RI（实部-虚部）或DB（dB-相位）形式。

• 端口名称定义（[Port Names]段）
• 网络数据的结构化
• 噪声参数的集成
• 每个端口不同基准阻抗的支持

实务上要注意的是端口号顺序。EM分析工具和SPICE工具的端口号分配规约可能不同。差动对的正负搞反的话，$S_{dd11}$ 的结果会完全不同。

从S参数到眼图从S参数到眼图

1. 获取通道S参数 — 通过EM分析或VNA测量获得.s4p
2. 质量检查 — 验证无源性、因果律、相反性。有违反就修正
3. 准备IBIS-AMI模型 — 定义Tx的均衡（FFE/De-emphasis）和Rx的均衡（CTLE/DFE）为AMI模型
4. 执行通道仿真 — 整合S参数+AMI模型，生成统计眼图（BER=$10^{-12}$等）
5. 裕度评估 — 判断眼开口量是否满足规范掩码

这整个流程就是"S参数基础的通道仿真"，是现代高速数字设计的核心。

常见失败及对策

S参数的软件对比

商用工具对比

选择指南

• 分析对象的结构 — 连接器、过孔等3D结构用HFSS/CST/Clarity。PCB配线图形为主用SIwave/Sigrity
• 与现有EDA流程的整合 — Allegro环境用Sigrity/Clarity最自然。Ansys EM环境用HFSS+SIwave。VNA用户用ADS
• 与通道仿真的统一 — 是否能从S参数到IBIS-AMI基础通道仿真一站式完成。Ansys Channel Design/Cadence SystemSI等

开源方面，Python的scikit-rf在S参数读取、转换、绘图、校准方面很全面，Touchstone文件的前处理和自定义分析脚本很有用。

S参数的前沿研究

112 Gbps及以上时代

• 介质损耗（Df）精度 — 毫米波段基板材料的Df成为支配因素。需要用宽频Debye/Djordjevic-Sarkar模型
• 导体表面粗糙度模型 — 导体损耗的30~50%来自表面粗糙。Huray/Hammerstad-Jensen模型的选择很重要
• EM分析网格爆炸 — 按λ/20标准，56 GHz要求单元尺寸<0.05 mm。不用HPC并行无法在实用时间内完成
• PAM4调制对应 — 从NRZ二值转到PAM4四值。SNR裕度减少9.5 dB，通道损耗预算更严格

机器学习辅助S参数补全

• 代理模型 — 参数化的S参数预测。以配线长、介电常数、层结构为参数，用神经网络快速推估S参数。设计空间搜索加速数千倍
• PINN（物理信息神经网络） — 把无源性和因果律融入网络损失函数，生成物理上一致的S参数
• 贝叶斯优化 — 在有限的3D EM分析结果基础上，寻找最优过孔结构和间距
• GAN进行数据扩充 — 当测量数据不足时，生成物理上合理的S参数数据

但这些多数还在研究阶段，设计验证还是要依靠传统的物理基础EM分析。AI适合"搜索"，但"验证"（目前）还需要信物理。

S参数故障排查

S参数品质诊断

1. 无源性检查 — 全频段特征值 ≤ 1 吗？有违反说明测量或分析质量有问题
2. 因果律检查 — IFFT后脉冲响应在 $t < 0$ 为零吗？群延迟不出现负值吗？
3. 相反性检查 — $|S_{12} - S_{21}|$ 足够小吗？很大的话端口定义或De-embedding有误
4. DC连续性 — $S_{21}(0)$ 物理上合理吗？DC耦合应|S21(0)| ≈ 1，AC耦合应|S21(0)| = 0
5. 高频外推 — 最高频率至少是比特速率的3倍吗？
6. 频率点密度 — 采样足够密集以捕捉共振峰吗？

常见错误及对策