| 查看: 3404 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hollya至尊木虫 (文坛精英)
葛城?{美
|
[交流]
电池或电容器等效电路模拟（1）已有9人参与
|
||
|
有很多虫友一直对EIS图的模拟或者拟合有疑问，一些软件又要钱，比如zview. 我觉得这个内容其实不是很复杂，就自己写模型模拟吧。 这里我写了一个简单的matlab代码可以模拟一阶等效电路模型。不知道应该发在计算模拟区还是电化学区，最后觉得可能电化学方面用的更多一些，就发到这里吧：） 这个模型的等效电路为R0-(RL||L)-(R1||C1)-W. R0是系统固有的电阻，RL||L, 是电感-电阻耦，R1||C1是电容电阻耦，W-是Warburg element，代表扩散阻抗。 如果赋值模型：R0=1欧姆; R1=1欧姆; C1=0.1法拉第; L=0.0001; RL=1欧姆; A=1; f=logspace(-2,3)（频率为10mHz-1000Hz) 模拟结果如图：） 版主给加个精华吧，哈哈：） ============================================ 下面是代码（函数形式）：） function [z]=ecm1wr(R0,R1,C1,RL,L,A,f) % This function calcultes the impedance of first order of equivalent % circuit model include Warburg diffusion element. % Model: R0-(RL||L)-(R||C)-W % R0: material, connection, electrolyte resistance % R: charge transfer resistance % C: equalent capacitance % L: inductance % f: experimental frequency % example: R0=1; R1=1; C1=0.1; L=0.0001; RL=1; A=1; f=logspace(-2,3) w=freq(f) zr=res(w,R1) zrl=res(w,RL) zc=cap(w,C1) zl=ind(w,L) zw=war(w,A) z=R0+1./(1./zrl+1./zl)+1./(1./zr+1./zc)+zw plot(real(z),-imag(z),'x') xlabel('Re_Z(Ohm)') ylabel('Im_z(Ohm)') title('Nyquist plot') end function w=freq(f) % transform f to w w=2*pi*f end function zc=cap(w,C) % capacitance zc=1./(j*w*C) end function zl=ind(w,L) % inductance zl=j*w*L end function zw=war(w,A) % Warburg element % A: Warburg coefficient zw=A*(1./w.^(1/2)+1./(j*w.^(1/2))) end function zr=res(w,R) % Resistance zr=R*ones(size(w)) end  untitled.jpg |
回复此楼» 本帖已获得的红花（最新10朵）
烟圈o0» 猜你喜欢
请问有用活性炭或者碳布做负极的吗
已经有0人回复
-
2 5-二氨基苯磺酸的氨基保护
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有170人回复
-
清华大学未来实验室 实习生招募
已经有2人回复
-
络合
已经有2人回复
-
锂离子电池电解液电化学窗口怎么测定？
已经有0人回复
-
氨基酸会被紫外杀菌灯破坏结构吗？具体是甘氨酸和蛋氨酸
已经有0人回复
-
上海交大工业催化方向高薪招聘博士后研究员
已经有18人回复
-
上海交大储能方向高薪招聘博士后研究员
已经有17人回复
-
祈福博后基金
已经有59人回复
-
年龄42岁，想读博怎么办？大家有什么好的建议？是不是很难了
已经有9人回复
小小淘金者
至尊木虫 (知名作家)
- 应助: 4 (幼儿园)
- 金币: 23010.1
- 散金: 1558
- 红花: 16
- 沙发: 66
- 帖子: 8365
- 在线: 1211.6小时
- 虫号: 3942665
- 注册: 2015-06-27
- 性别: GG
- 专业: 基础物理学
6楼2020-03-04 08:38:21
5