| 查看: 3399 | 回复: 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» 猜你喜欢
硫化物电解质压片
已经有9人回复
-
请问大家锂电三元层间距的计算
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有59人回复
-
BO4的YQ答辩通知发布了吗？
已经有6人回复
-
“催化材料信息学”博士后 松山湖材料实验室联合招聘！待遇从优！
已经有5人回复
-
“催化材料信息学”博士后招聘！松山湖材料实验室联合招聘！待遇从优！
已经有5人回复
-
“催化材料信息学”博士后招聘！松山湖材料实验室联合招聘！待遇从优！
已经有5人回复
-
“催化材料信息学”博士后招聘！松山湖材料实验室联合招聘！待遇从优！
已经有5人回复
-
“催化材料信息学”博士后招聘！松山湖材料实验室联合招聘！待遇从优！
已经有5人回复
-
催化材料信息学博士后招聘 总包40-60w/年 松山湖材料实验室刘利峰研究员、张磊教授
已经有2人回复
-
催化材料信息学博士后招聘 总包40-60w/年 松山湖材料实验室刘利峰研究员、张磊教授
已经有1人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
另外，可不可以请教如何编写拟合的matlab代码？在你上述等效电路的基础上，通过拟合实测数据，获得电路图里的各个参数，一般都是用最小二乘法，谢谢 发自小木虫Android客户端 |
9楼2020-03-04 14:57:27
2楼2020-03-04 07:01:35
hollya
至尊木虫 (文坛精英)
葛城?{美
- ECEPI: 1
- 应助: 744 (博后)
- 金币: 43445.1
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.9小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
3楼2020-03-04 07:11:50
烟圈o0
木虫 (著名写手)
- 应助: 177 (高中生)
- 金币: 3155.8
- 散金: 1035
- 红花: 43
- 沙发: 4
- 帖子: 2643
- 在线: 1468小时
- 虫号: 986567
- 注册: 2010-03-31
- 性别: GG
- 专业: 电化学
4楼2020-03-04 07:56:01
hollya
至尊木虫 (文坛精英)
葛城?{美
- ECEPI: 1
- 应助: 744 (博后)
- 金币: 43445.1
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.9小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
5楼2020-03-04 08:12:31
小小淘金者
至尊木虫 (知名作家)
- 应助: 4 (幼儿园)
- 金币: 22984.3
- 散金: 1558
- 红花: 16
- 沙发: 66
- 帖子: 8359
- 在线: 1211.6小时
- 虫号: 3942665
- 注册: 2015-06-27
- 性别: GG
- 专业: 基础物理学
6楼2020-03-04 08:38:21
hollya
至尊木虫 (文坛精英)
葛城?{美
- ECEPI: 1
- 应助: 744 (博后)
- 金币: 43445.1
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.9小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
7楼2020-03-04 08:39:43
muxinmuxin
新虫 (知名作家)
- 应助: 0 (幼儿园)
- 金币: 4187.9
- 散金: 6524
- 红花: 57
- 沙发: 67
- 帖子: 9841
- 在线: 1498.9小时
- 虫号: 7526450
- 注册: 2017-11-26
- 性别: GG
- 专业: 工程热物理与能源利用
8楼2020-03-04 11:58:14
10楼2020-03-04 17:17:06
20