方法说明 / 2026.07.24

为什么 EIS 频率和 DRT 弛豫时间不能直接对照？

DRT 横轴不是把 EIS 中的每个频率点逐个换算得到的。要看懂它，需要分清物理上的特征时间和程序实际使用的数值网格。

2026 年 7 月 24 日Oparic 团队约 3 分钟阅读

先说结论

如果 DRT 图上的弛豫时间范围和你按 EIS 频率算出的范围对不上，先不用判断为异常。理想 RC 公式告诉我们一个峰大约对应什么频率；DRT 横轴则是程序为了求解整条分布建立的数值网格。两者用途不同，端点自然不必相同。

对理想的单一 RC 弛豫过程，响应峰附近满足：

\omega\tau=1,\qquad \omega=2\pi f

所以，峰的特征时间可以近似写成：

\tau\approx\frac{1}{2\pi f}

它表达的是“这个弛豫过程在哪个频率附近响应最明显”：频率越高，对应的过程越快， $\tau$ 越小；频率越低，对应的过程越慢， $\tau$ 越大。

但这不是把每个 EIS 频率点逐个改写成一个 $\tau$ 点。DRT 会同时使用整段阻抗谱，反推出一条连续的弛豫时间分布。真实过程还可能彼此重叠、展宽，或者并不符合理想单一 RC 模型，因此这个公式适合估算峰位，不适合直接规定 DRT 横轴的端点。

Oparic 当前先按每个测量频率建立数值求解节点：

\tau_i=\frac{1}{f_i}

这里的 $1/f$ 是求解器使用的数值坐标，不是把理想 RC 的特征频率改成了 $f=1/\tau$ 。实际计算阻抗响应时仍然使用 $\omega=2\pi f$ ；解释某个峰对应的物理特征频率时，仍应使用 $f\approx1/(2\pi\tau)$ 。

如果测得的最低和最高频率分别记作 $f_{\min}$ 和 $f_{\max}$ ，Oparic 输出的时间网格范围是：

\tau_{\mathrm{grid,min}}=\frac{10^{-1/2}}{f_{\max}},\qquad \tau_{\mathrm{grid,max}}=\frac{10^{1/2}}{f_{\min}}

换句话说，程序先取倒数频率范围，再在对数坐标两端各扩展半个数量级，并在更密的网格上画出拟合结果。这会把分布的边缘也显示出来，但不会凭空增加实验信息。绘图软件还可能在坐标轴两端留白，因此最外侧刻度也不一定就是网格的精确端点。

如果你在 drt_data.csv 中看到 freq_hz，当前这一列记录的是时间网格的倒数 $1/\tau$ 。它不是理想 RC 峰的特征频率；后者应按 $1/(2\pi\tau)$ 换算。

DRT 是从有限且带噪声的阻抗数据中反推连续分布。网格可以向外延伸，数据提供的约束却不会跟着变多，所以接近两端的峰通常比中间区域更需要验证。

先检查高频端的电缆和夹具寄生电感、相位误差与仪器噪声。即使程序同时拟合了串联电感，也不能保证所有高频误差都已消除。

先检查低频响应是否测完整，以及长时间测试中温度、荷电状态或体系本身是否发生漂移。非稳态数据可能形成看似合理的慢过程。

都应回看 Nyquist/Bode 原始数据、阻抗重建和残差，并比较重复测试。条件允许时扩展频率范围：真实且稳定的峰不应只是跟着新的横轴边界移动。

要点：横轴范围是求解器的工作范围，不是可信度保证。峰是否值得解释，取决于它是否被数据充分约束、能否重复，以及是否与原始 EIS 和实验机理相互印证。