3、电磁仿真基础：S参数、Y参数、Z参数的定义与转换，以及电感提取的常用仿真设置
在RFIC（射频集成电路）设计中，片上电感（Spiral Inductor）的性能直接决定了VCO（压控振荡器）、LNA（低噪声放大器）和PA（功率放大器）等关键模块的指标。然而，片上电感的寄生效应（如衬底损耗、金属趋肤效应、邻近效应）无法通过简单的解析公式精确建模，必须依赖电磁（EM）仿真工具进行参数提取。本章将深入剖析S参数、Y参数、Z参数的本质定义与相互转换关系，并给出针对片上电感提取的仿真设置流程与工程实践要点。理解这些参数是后续进行等效电路建模、去嵌入（De-embedding）和设计优化的重要基础。
上图展示了从端口激励到参数后处理的完整链路。在RFIC设计中，我们通常将片上电感视为一个二端口网络（或单端口网络，取决于接地方式）。电磁仿真工具首先在指定频段内求解Maxwell方程组，得到端口间的S参数（散射参数）。随后，通过数学变换，我们可以获得Y参数（导纳参数）和Z参数（阻抗参数）。对于电感提取，最核心的步骤是将仿真得到的S参数转换为Z参数或Y参数，进而从虚部中解析出电感值L和品质因数Q。下面我们将逐一展开这些参数的定义、转换公式及仿真设置细节。
3.1 S参数、Y参数、Z参数的核心定义与物理意义
（1）S参数（散射参数）：S参数是基于入射波和反射波的概念定义的，非常适合高频测量。对于二端口网络，S11表示端口1的反射系数（输入回波损耗），S21表示从端口1到端口2的正向传输系数（插入损耗/增益）。S参数可以直接由矢量网络分析仪（VNA）测量得到，也是EM仿真工具的标准输出格式（Touchstone文件，即.s2p文件）。S参数的优点是物理意义直观，且在高频下易于测量；缺点是它依赖于参考阻抗（通常为50Ω），且不便于直接用于电路仿真中的串并联元件计算。
（2）Y参数（导纳参数）：Y参数定义了端口电压与端口电流之间的关系，即 [I] = [Y]·[V]。Y11是端口1的输入导纳（端口2短路），Y21是正向转移导纳。对于片上电感，Y参数特别有用，因为电感的一端通常接地（衬底），此时Y11直接反映了电感与地之间的等效导纳。通过Y11的虚部，我们可以直接提取电感值：L = 1 / (2πf · Im(Y11))。Y参数的优势在于，并联元件的计算非常方便（并联导纳直接相加）。
（3）Z参数（阻抗参数）：Z参数定义了端口电压与端口电流之间的关系，即 [V] = [Z]·[I]。Z11是端口1的输入阻抗（端口2开路），Z21是正向转移阻抗。Z参数在分析串联元件时非常直观。对于电感提取，Z11的虚部即为感抗：L = Im(Z11) / (2πf)。Z参数适合分析电感在低频下的行为，但在高频下，由于开路条件难以实现，Z参数通常由S参数或Y参数转换而来。
知识结构对比表：三种参数的特性与适用场景
| 参数类型 | 定义基础 | 测量/仿真便捷性 | 电感提取公式 | 主要应用场景 |
|---|---|---|---|---|
| S参数 | 入射波/反射波 | 最易测量（VNA） | 需先转换为Y或Z | 高频测试、系统级联、阻抗匹配 |
| Y参数 | 短路导纳 | 适合并联结构分析 | L = 1/(2πf·Im(Y11)) | 电感接地模型、π型等效电路 |
| Z参数 | 开路阻抗 | 适合串联结构分析 | L = Im(Z11)/(2πf) | T型等效电路、低频建模 |
3.2 参数转换的数学方法与工程实现
在电磁仿真后处理中，我们通常直接获得S参数。为了提取电感值，需要将其转换为Y参数或Z参数。转换公式基于归一化阻抗和单位矩阵。设参考阻抗为Z0（通常为50Ω），S参数矩阵为[S]，则Y参数矩阵可通过以下公式计算：
Y = Y0^(1/2) · (I - S) · (I + S)^(-1) · Y0^(1/2)其中，Y0 = 1/Z0 为参考导纳，I为单位矩阵。对于二端口网络，该公式可以展开为显式表达式。同理，Z参数矩阵为Y参数的逆矩阵：Z = Y^(-1)。在Python的scikit-rf库或MATLAB的RF工具箱中，这些转换均有现成函数。但在工程实践中，我们更关心的是如何从转换后的参数中准确提取电感的L和Q值。
操作流程：从S参数到电感值的标准步骤
- 数据准备：从EM仿真工具（如HFSS、Cadence EMX）导出.s2p文件，确保频点范围覆盖目标工作频率（例如，对于5GHz的RFIC，建议扫描1GHz~15GHz）。
- 参数转换：使用脚本或工具将S参数转换为Z参数（或Y参数）。推荐使用scikit-rf库，代码示例如下：
import skrf as rf
# 读取Touchstone文件
net = rf.Network('inductor.s2p')
# 转换为Z参数
z_net = net.z
# 提取Z11的虚部
freq = net.frequency.f
Z11_imag = z_net[:, 0, 0].imag
# 计算电感值
L = Z11_imag / (2 * 3.14159 * freq)
# 计算Q值（单端激励，假设端口2接地）
Q = z_net[:, 0, 0].imag / z_net[:, 0, 0].real- 去嵌入处理：如果仿真结构中包含了测试焊盘（Pad）和互连线，必须通过Open/Short/Through结构进行去嵌入，以消除寄生影响。这一步至关重要，否则提取的L值会偏大，Q值会偏小。
- 数据验证：将提取的L和Q值绘制成曲线，观察其随频率的变化趋势。理想情况下，L值在低频段应保持平稳，Q值应呈现先升后降的峰值特性。
经验提醒：在参数转换时，务必注意参考阻抗的一致性。如果S参数是在非50Ω系统下仿真的，转换前需重新归一化。此外，当电感工作频率接近自谐振频率（SRF）时，Z11的虚部会从正变为负，此时电感表现为电容特性，提取公式不再适用。
3.3 电感提取的常用仿真设置与实战案例
在RFIC设计中，片上电感通常采用顶层厚金属（如Al或Cu）制作，以降低串联电阻。仿真设置的合理性直接影响提取精度。以下是针对片上螺旋电感的EM仿真关键设置清单：
| 设置项 | 推荐值/方法 | 说明 |
|---|---|---|
| 求解器类型 | 全波（FEM）或矩量法（MoM） | 对于5GHz以上频段，建议使用全波求解器以考虑辐射效应 |
| 端口类型 | Wave Port（波端口）或Lumped Port（集总端口） | 对于片上电感，推荐使用Lumped Port并设置端口阻抗为50Ω |
| 参考地 | 设置理想导体边界（PEC）或衬底网格 | 需包含衬底损耗（电阻率约10 Ω·cm） |
| 频点设置 | 线性扫描，步长0.1GHz或0.5GHz | 确保覆盖自谐振频率（SRF） |
| 网格剖分 | 趋肤深度自适应剖分 | 在金属厚度方向至少剖分2~3层网格 |
实战案例：5GHz VCO用片上电感提取
假设我们需要设计一个中心频率为5GHz的VCO，其谐振回路需要一颗感值约为2nH、Q值大于15的电感。我们采用HFSS进行仿真，结构为3.5圈八边形电感，外径120μm，金属宽度10μm，间距2μm，顶层金属厚度3μm。
- 初始现象：初次仿真后，提取的L值在5GHz处为2.3nH，但Q值仅为10，低于设计指标。
- 分析判断：Q值偏低通常由金属损耗和衬底损耗共同导致。检查仿真设置发现，衬底电阻率设置过高（误设为100 Ω·cm），且未开启金属粗糙度模型。
- 实际操作：将衬底电阻率修正为10 Ω·cm，并添加表面粗糙度（RMS=0.3μm）模型。重新仿真后，5GHz处Q值提升至16.5，L值微调至2.1nH，满足设计要求。
- 结果对比：修正前后，L值变化不大（2.3nH→2.1nH），但Q值显著改善（10→16.5）。这说明衬底损耗对Q值影响巨大，而对L值影响较小。
- 风险与代价：增加粗糙度模型会显著增加仿真时间（约30%），但为了精度，这是必要的代价。
- 复盘经验：在电感EM仿真中，衬底参数和金属表面模型必须与实际工艺吻合。建议在工艺设计套件（PDK）中直接调用经过校准的衬底文件。
风险提示：切勿直接使用EM仿真工具默认的“真空”环境。必须手动添加衬底（Silicon）和介质层（Oxide）的介电常数与损耗角正切。否则，提取的Q值会严重偏大，导致流片后性能不达标。
3.4 常见误区与总结复盘
在电感参数提取过程中，工程师常犯以下错误：
- 误区一：忽略去嵌入。直接提取包含Pad和引线在内的S参数，导致L值偏大。解决方法是制作Dummy结构（Open/Short）进行两步去嵌入。
- 误区二：端口设置不当。使用Wave Port时，端口尺寸设置过大或过小都会引入寄生模式。对于小尺寸片上电感，Lumped Port更合适。
- 误区三：频点范围过窄。如果扫描频率未超过自谐振频率，就无法观察到L值的下降趋势，从而误判电感带宽。
- 误区四：混淆Y参数与Z参数提取公式。Y参数提取公式为L = 1/(2πf·Im(Y11))，而Z参数提取公式为L = Im(Z11)/(2πf)。两者在数值上等价，但若网络中存在较大并联寄生电容，使用Y参数提取的L值更接近物理定义。