这款软件是由Francesco Ciucci课题组开发的,软件下载地址 https://github.com/ciuccislab (链接打不开?建议学会科学上网)
目前提供了MATLAB和Python两个版本,(点击蓝色文字可直接下载)
说明书:User manual_v7
Python版本:Python-Download – DRTTOOLS
matlab版本:matlab-Download – DRTTOOLS
常规的EIS数据解析方法仍然需要先验假设的等效电路的拟合技术,如通过利用Zview软件构建电路基元从而进行拟合,正往往需要一定的经验并需要综合考虑文献中的报道来构建电路。DRT技术作为近两年比较新兴的分析手段,可通过反卷积技术提取电化学系统中的弛豫时间分布。即不通过先验假设而是通过理论计算的方法进行去卷积,从而得到各部分的阻抗。
经典文献
Distribution of Relaxation Times Analysis of High-Temperature PEM Fuel Cell Impedance Spectra
Matlab版本
操作示意图如下:
1、打开matlab,打开找到下载文件中的DRTtools.m,打开后点击运行
2、运行后点击import,弹出窗口中找到示例阻抗数据,import file samples→text files→LIB_data.txt
3、点击Run中的Simple Start,即可出现DRT数据
4、如何绘图如何得到积分阻抗?
没想到这个软件的关注度这么高,很多来私信留言的同学,现在我将方法写在这里,如有错误请批评指正
这是计算的数据,同学们可进行参照
首先我们要明确的是在DRT计算结果中的频率f和τ的转换关系为τ=1/(2πf),不用质疑,没有问题。这样就轻松画出来log10(f)和γ(τ)的关系图了
其次最重要的就是积分问题,怎么得到每个峰对应的阻抗数值呢?观察式子
从小学数学我们能知道,其中γ(lnτ)=τG(τ)。上述这串公式的详细推导过程建议看这篇文献Influence of the Discretization Methods on the Distribution of Relaxation Times Deconvolution_ Implementing Radial Basis Functions with DRTtools.pdf
因此做积分时需要对τ先求ln(τ) (不是log10(τ)!!)再进行积分,如下表所示的C列,对C列和D列画图再进行每个峰的积分即可(注意此处的横坐标为ln(τ) !)
我们也可以求出累计的阻值(注意此处的横坐标仍然是ln(τ) )
我们也可以用横坐标为频率的图,比较一下我们累计计算的阻值和实际Zview分析得到的|Z|(注意得到的|Z|我减去了0.1124),结果如下图所示(注意此处的横坐标是log10(f) )
至于文章里的颜色填充啊什么的,用XPSPEAK41去分个峰,再导入到originpro里画个图就好了
Python版本(推荐)
操作示意图如下:
推荐下载Anaconda软件 (如果选用VS code依赖的库有点多,需要挨个安装库)
下载地址:https://www.anaconda.com/download/success
下载完成后,安装过程到后半段有点慢尤其是到接近95%的时候,耐心等待,估计10分钟吧
安装完成后打开开始菜单
打开Anaconda Prompt后,
输入:pip install cvxpy
安装完成cvxpy库之后,打开上面图片中的Spyder软件
把下载的pyDRTtool.py文件直接拖进来
点击Run→Run file
完活
请教一下博主