scikit-rf绘制串联RLC模型S2P文件Z参数幅值图异常但Y参数图正常的原因排查咨询

阿华AIGC实验室

2026-4-7

我理解你遇到的这个问题确实让人困惑——明明Y参数图完全符合预期，Z参数却出现噪声大、数值异常的情况。结合scikit-rf的特性、S参数的物理意义以及Murata这类元件模型的常见细节，咱们来拆解核心原因和常见陷阱：

核心原因：双端口Z/Y参数的物理意义与你的需求不匹配

这是最关键的问题，你混淆了双端口网络参数的定义和你实际想要的串联元件本身的阻抗：

Y参数（Y11）的巧合正确性：双端口Y11的定义是端口2短路时，端口1的输入导纳（Y11 = I₁/V₁，V₂=0）。对于Murata提供的串联模型S2P文件，元件是串联在两个端口之间的，当端口2短路时，端口1的输入导纳恰好等于元件的导纳（Y = 1/Z），所以plot_y_mag(m=0,n=0)绘制的Y11幅值就是1/|Z|，完全符合你对元件导纳的预期。
Z参数（Z11）的物理意义偏差：双端口Z11的定义是端口2开路时，端口1的输入阻抗（Z11 = V₁/I₁，I₂=0）。对于串联元件的双端口模型，当端口2开路时，输入阻抗是无穷大（因为元件串联在开路端口），而低频时电容的阻抗本身很大，会导致Z11的计算出现数值爆炸（S11接近1，Z参数转换公式Z = Z0*(1+S)/(1-S)的分母趋近于0），这就是你看到的“噪声大、数值异常”的直接原因——你绘制的根本不是元件本身的阻抗，而是开路状态下的双端口输入阻抗。

其他常见陷阱与排查点

1. 参考阻抗（Z0）的不匹配或误读

scikit-rf加载S2P文件时会自动读取文件头的Z0字段（比如# Z0 50），但如果文件的实际参考阻抗与读取值不一致，会直接导致Z/Y参数的转换错误：

可以通过print(ntwk.z0)检查加载的参考阻抗是否与文件标注一致（Murata的S参数文件通常默认50Ω）。
如果存在误读，可以手动指定：ntwk.z0 = 50（确保与文件的参考阻抗匹配）。

2. 数值精度问题（S参数接近1时的计算不稳定）

当元件阻抗远大于参考阻抗（比如低频时的电容），S11会趋近于1，此时Z参数转换公式Z = Z0*(1+S11)/(1-S11)的分母趋近于0，会引发数值溢出或噪声，进一步放大Z11的异常。而Y参数的转换公式Y = (1-S)/(Z0*(1+S))此时分子分母同时趋近于0，比值反而稳定（等于1/Z），这也是Y参数图正常的另一个原因。

3. 对“串联模型S2P”的误解

Murata的“串联简单模型”S2P是双端口伪模型，用来模拟元件串联在两个端口之间的行为，并非真正的双端口网络。你需要的是元件本身的阻抗，而非双端口网络的Z参数，因此直接调用plot_z_mag(m=0,n=0)从一开始就偏离了需求。

解决方案：绘制元件本身的阻抗而非双端口Z11

你可以通过以下几种正确的方式获取并绘制串联元件的阻抗：

方法1：从Y参数推导元件阻抗（最直接）

既然Y11的幅值是1/|Z|，直接取倒数即可得到元件阻抗的幅值：

import skrf as rf
import matplotlib.pyplot as plt
import numpy as np

ntwk = rf.Network("GCM31CR71A226KE02_DC0V_25degC_series_simple.s2p")
# 从Y11计算元件阻抗：Z = 1/Y11
z_series = 1 / ntwk.y[:, 0, 0]
z_mag = np.abs(z_series)

# 绘制图形
fig, ax = plt.subplots()
ax.plot(ntwk.frequency.f, z_mag)
ax.set_xscale('log')
ax.set_title('Magnitude of Series Capacitor Impedance')
ax.set_xlabel('Frequency (GHz)')
ax.set_ylabel('Impedance Magnitude (Ohms)')
ax.grid(True, which='both', linestyle=':', linewidth=0.8)
plt.show()

方法2：将双端口转换为单端口（短路端口2）

通过短路端口2，将双端口网络转换为单端口网络，此时单端口的输入阻抗就是元件本身的阻抗：

import skrf as rf
import matplotlib.pyplot as plt

ntwk = rf.Network("GCM31CR71A226KE02_DC0V_25degC_series_simple.s2p")
# 短路端口2，得到单端口网络
ntwk_short = ntwk.short(port=2, mode='z0')
# 绘制单端口输入阻抗的幅值
ntwk_short.plot_z_mag()
ax = plt.gca()
ax.set_xscale('log')
ax.set_title('Magnitude of Series Capacitor Impedance (Short Port 2)')
ax.set_xlabel('Frequency (GHz)')
ax.set_ylabel('Impedance Magnitude (Ohms)')
ax.grid(True, which='both', linestyle=':', linewidth=0.8)
plt.show()

方法3：从S21直接计算元件阻抗

根据串联双端口模型的S参数与元件阻抗的推导关系：Z = 2*Z0*(1/S21 - 1)，可以直接计算：

import skrf as rf
import matplotlib.pyplot as plt
import numpy as np

ntwk = rf.Network("GCM31CR71A226KE02_DC0V_25degC_series_simple.s2p")
z0 = ntwk.z0[0]  # 获取参考阻抗
s21 = ntwk.s[:, 1, 0]  # 提取S21参数
z_series = 2 * z0 * (1/s21 - 1)
z_mag = np.abs(z_series)

# 绘制图形
fig, ax = plt.subplots()
ax.plot(ntwk.frequency.f, z_mag)
ax.set_xscale('log')
ax.set_title('Magnitude of Series Capacitor Impedance (from S21)')
ax.set_xlabel('Frequency (GHz)')
ax.set_ylabel('Impedance Magnitude (Ohms)')
ax.grid(True, which='both', linestyle=':', linewidth=0.8)
plt.show()