S参数解析 — 信号完整性(SI)评估中的散射矩阵理论与实践
理论与物理
什么是S参数
老师，S参数是用来表示高频电路特性的对吧？我经常听到S11、S21这些，但它们具体代表什么含义，说实话我有点模糊……
问得好。简单来说，S参数（Scattering Parameter，散射参数）是用复数表示的当信号输入到高频网络（传输线、连接器、封装等）时，有多少被反射，有多少能通过的量。
用日常例子比喻的话，对着隧道喊话时，“回声返回的声音”相当于S11，“传到另一头的声音”相当于S21。
原来如此，是反射和透射的组合啊。但是为什么不用普通的电阻或阻抗，而特意要用S参数呢？
低频时可以用电压、电流来描述电路，但到了GHz频段，波长变得与布线长度相当，因此需要将信号作为“行波、反射波”来处理，而不是“电压、电流”。S参数正是基于这种行波的描述方法。
具体来说，PCIe Gen5（32 GT/s）的奈奎斯特频率是16 GHz，波长在电路板上约为6 mm。即使是10 cm的布线，也有十几个波长那么长，所以无法用集总参数电路来描述。
散射矩阵的定义
数学上是怎么定义的呢？
$N$端口网络的S参数，是用矩阵来描述各端口的入射波 $a_i$ 与反射波 $b_i$ 之间的关系。对于2端口的情况：
这里各元素的含义如下：
- $S_{11}$ — 端口1的输入反射系数。当端口2用基准阻抗 $Z_0$ 端接时，端口1入射的波有多少被反射
- $S_{21}$ — 端口1→端口2的正向传输系数。信号有多少能通过
- $S_{12}$ — 端口2→端口1的反向传输系数（对于无源元件，通常 $S_{12} = S_{21}$）
- $S_{22}$ — 端口2的输出反射系数
入射波和反射波定义为归一化功率波：
$Z_0$ 是指50 Ω吗？
大多数情况下是50 Ω。因为网络分析仪（VNA）的端口阻抗是50 Ω。但在SI分析中，有时也会用100 Ω差分对特性阻抗进行归一化，所以使用S参数时，务必确认基准阻抗非常重要。
回波损耗与插入损耗
在实际的SI设计中，经常看到“回波损耗”、“插入损耗”这些词。它们和S11、S21有什么关系？
将S参数的绝对值转换为dB，就是“损耗”表示法。
举一些实际工作中的参考标准：
- 回波损耗 > 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等），相同长度下可以抑制到一半左右。
也就是说，插入损耗越接近零越好，回波损耗越大越好，对吧？符号方向相反，有点容易混淆……
没错。这是实际工作中经常混淆的一点。另一点需要注意的是符号的惯例。有些工具中IL = 20 log|S21|（没有负号的负值）来表示。一定要确认自己看到的图表纵轴定义。
无源性与因果律
前辈说“S参数一定要检查无源性和因果律”，这是什么意思？
这是SI分析中最重要的质量标准。有两个物理约束条件：
1. 无源性条件 — 网络不产生能量：
2端口情况：$|S_{11}|^2 + |S_{21}|^2 \leq 1$（对所有频率成立）
2. 因果律条件 — 输出不会出现在输入之前：
- 时域脉冲响应 $s(t)$ 在 $t < 0$ 时为零
- 频域上，实部和虚部满足希尔伯特变换对（Kramers-Kronig关系）
如果这些条件被违反会怎样？
如果直接导入到SPICE或IBIS-AMI中，会导致非物理性的振荡、仿真发散、眼图的异常过冲。这是实际工作中多次遇到的典型问题。
对策包括：
- 应用Ansys SIwave、Keysight ADS、Cadence Sigrity等工具内置的无源性强制算法
- 因果律违反通常是由于测量数据带宽不足或低频外推错误造成的。需要对DC点和高频进行适当的外推
混合模式S参数
我听说在使用差分信号的设计中，有和普通S参数不同的表示方法。
对于差分对，会将4端口的标准S参数（单端）转换为混合模式S参数来使用。分为差模（DD）、共模（CC）和模式转换（DC/CD）：
- $S_{dd21}$ — 差分插入损耗（SI设计者最关注的量）
- $S_{dd11}$ — 差分回波损耗
- $S_{cd21}$ — 差模→共模转换（衡量偏斜或不对称性的指标）
- $S_{cc21}$ — 共模插入损耗（影响EMI特性）
例如，USB4或PCIe Gen6的规范中定义了 $S_{dd21}$ 和 $S_{dd11}$ 的掩模（极限值曲线），不满足此条件的通道即为不符合规范。
散射矩阵的数学性质总结
- 互易性：对于无源线性网络，$\mathbf{S} = \mathbf{S}^T$（等于其转置）。即 $S_{ij} = S_{ji}$。除了铁氧体器件等非互易元件外，SI设计中处理的几乎所有结构都满足此条件。
- 无耗性：$\mathbf{S}^H \mathbf{S} = \mathbf{I}$（酉矩阵）。理想无耗网络的条件，但实际PCB传输线存在介质损耗和导体损耗，因此严格来说不成立。
- 无源性：$\sigma_{\max}(\mathbf{S}) \leq 1$（最大奇异值小于等于1）。等价于所有特征值的绝对值小于等于1。
- 因果律：傅里叶逆变换后的时间响应在零延迟之前为零。实部与虚部满足Kramers-Kronig关系。
基准阻抗与S参数的转换
将基准阻抗为 $Z_{01}$ 的S参数重新归一化到 $Z_{02}$：
$$ \mathbf{S}' = \mathbf{A} (\mathbf{S} - \mathbf{\Gamma}) (\mathbf{I} - \mathbf{\Gamma} \mathbf{S})^{-1} \mathbf{A}^{-1} $$其中 $\Gamma_i = (Z_{02,i} - Z_{01,i}) / (Z_{02,i} + Z_{01,i})$。将50 Ω测量数据转换为100 Ω差分分析时需要用到。
数值解法与实现
S参数的获取方法
S参数在CA分析中是怎么求出来的呢？我很难想象像结构分析那样划分网格求解……
主要有三种方法：
- 3D全波电磁解析（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）对精度至关重要
如果是3D EM解析，具体是如何得到S参数的呢？
以HFSS基于FEM的解法为例：
- 用四面体网格离散化分析区域。电磁场用边单元（Nedelec单元）进行插值
- 求解麦克斯韦方程组的弱形式（变分形式），计算各频率下的E场和H场分布
- 从端口面的模式场与解的内积中，提取入射波振幅 $a_i$ 和反射波振幅 $b_i$
- 在端口 $j$ 以外均匹配端接的条件下，由 $S_{ij}(f) = b_i / a_j$ 构成S参数
简而言之，就是从求解麦克斯韦方程组的结果中提取端口面的波振幅。
频域与时域的转换
S参数是频域数据对吧？像眼图这样的时域评估，是如何转换过去的呢？
通过逆傅里叶变换（IFFT）转换为脉冲响应 $s_{21}(t)$，再与比特模式进行卷积。这是通道仿真的基本流程：
应用IFFT时的实际注意事项：
- 对DC点（0 Hz）的数据进行适当外推。如果缺失，IFFT结果会产生偏移
- 高频侧的外推 — 希望数据的最高频率达到奈奎斯特频率的3~5倍
- 应用窗函数 — 使用高斯窗或凯撒窗来抑制吉布斯现象（振铃）
- 等间隔频率采样 — IFFT需要等间隔数据。非等间隔数据需要先进行插值
去嵌入技术
用VNA测量时，不是连电缆和测试夹具的影响也包含进去了吗？怎么去除这些影响呢？
这就是去嵌入技术。为了只提取DUT（被测器件）的S参数，需要将已知的夹具部分S参数去除。
代表性的方法：
- TRL校准（直通-反射-传输线） — 从三个已知标准件移动校准面。在板级测量中可靠性最高
- AFR（自动夹具移除） — 仅通过直通结构来估计夹具特性的算法
- 2xThru法 — 从两倍长的夹具结构，通过时域选通分离出DUT
在高速SerDes的规范验证（如CEI-112G等）中，测试板的去嵌入精度直接影响合规性判定。如果去嵌入出错，会导致“不符合规范”的误判。
S参数与其他网络参数的关系
除了S参数，还有Z（阻抗）、Y（导纳）、ABCD（传输）参数。它们之间可以相互转换，根据用途区别使用：
- Z参数 — 可以直接读取电路元件的阻抗。便于构建等效电路模型
- ABCD参数 — 可以通过矩阵乘积计算级联网络的整体特性。最适合多级滤波器设计
- S参数 → Z参数转换：$\mathbf{Z} = Z_0 (\mathbf{I} + \mathbf{S})(\mathbf{I} - \mathbf{S})^{-1}$
实践指南
Touchstone文件的运用
S参数的数据文件，我经常看到一种叫“Touchstone”的格式。比如.s2p、.s4p。这是什么机制？
Touchstone（旧称Citifile）是S参数数据交换的行业标准格式。扩展名的数字代表端口数：
- .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”表示幅度-相位格式，“R 50”表示基准阻抗50 Ω。此外也可以使用RI（实部-虚部）或DB（dB-相位）格式。
Touchstone 2.0（2009年制定）增加了以下内容：
- 端口名称的定义（[Port Names]部分）
- 网络数据的结构化
- 噪声参数的整合
- 支持每个端口不同的基准阻抗
实际工作中的注意事项：注意端口编号的顺序。EM解析工具和SPICE工具对端口编号的分配约定可能不同。
なった
詳しく
報告