应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 化学化工区 » 电化学 » 化学电源 » 电池或电容器等效电路模拟（1）
101/1返回列表
查看: 3405  |  回复: 19
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

电池或电容器等效电路模拟（1）
untitled.jpg
回复此楼

» 本帖已获得的红花（最新10朵）

烟圈o0 VStart

» 猜你喜欢
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
1楼2020-03-04 03:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
2楼: Originally posted by 糖果里的L at 2020-03-04 07:01:35
赞

谢谢

发自小木虫IOS客户端
回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
3楼2020-03-04 07:11:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
4楼: Originally posted by 烟圈o0 at 2020-03-04 07:56:01
优秀

承蒙照顾吖

发自小木虫IOS客户端
回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
5楼2020-03-04 08:12:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
6楼: Originally posted by 小小淘金者 at 2020-03-04 08:38:21
厉害

术业有专攻吖

发自小木虫IOS客户端
回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
7楼2020-03-04 08:39:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
10楼: Originally posted by VStart at 2020-03-04 17:17:06
1
硬核彗星！真想和你个专业！:

转行，快！

发自小木虫IOS客户端
回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
12楼2020-03-04 19:43:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
8楼: Originally posted by muxinmuxin at 2020-03-04 11:58:14
不会灌水的美女不是好程序员

听着像：不会开车的医生不是好厨子：）
水区人才济济，别的板块也要人才吖

发自小木虫IOS客户端
赞一下 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
13楼2020-03-04 20:01:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
9楼: Originally posted by 糖果里的L at 2020-03-04 14:57:27
另外，可不可以请教如何编写拟合的matlab代码？在你上述等效电路的基础上，通过拟合实测数据，获得电路图里的各个参数，一般都是用最小二乘法，谢谢
...

可以，我会再写个（2）

发自小木虫IOS客户端
赞一下 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
14楼2020-03-04 20:02:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
16楼: Originally posted by 当年刀刀 at 2020-03-05 01:39:20
我虽然学的电气  觉得楼主是学霸

这种级别的电路，在你电气专业里都是幼儿园级别的吧：）
赞一下 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
17楼2020-03-05 03:47:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美
引用回帖:
18楼: Originally posted by 当年刀刀 at 2020-03-06 00:32:47
对楼主的爱慕又多了一分
...

这个版块也太冷清了，要不是水区跑过来的虫友撑场面，都没有人瞥一眼的
赞一下 回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
19楼2020-03-06 00:39:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hollya 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国基在研影响申请结果吗 +6 WOWO159357 2024-06-26 6/300 2024-06-26 17:11 by wspglt
[基金申请] 博后网上办公系统维护 +12 取名字烦人 2024-06-25 18/900 2024-06-26 16:48 by 海阔天空999
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[基金申请] 基金评审时，评审专家还回去看申请人代表作之外的文章吗？ +5 孤单12站 2024-06-26 6/300 2024-06-26 16:05 by 孤单12站
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +9 地理学1995 2024-06-24 13/650 2024-06-26 16:02 by kyukitu
[基金申请] 博后面上和特助今天出吗？ +41 逗您玩 2024-06-21 78/3900 2024-06-26 16:00 by HAPPY_0225
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[有机交流] 做什么表征可以检测塑料中碳氢氧的含量 5+3 pzr的sci之路 2024-06-25 5/250 2024-06-26 10:13 by zqdsb
[有机交流] 过柱子，产品，杂质在是 pe:ea=100:1 也一起出来？ +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[有机交流] 对苯乙烯磺酰氯的合成机理 25+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[基金申请] 能看出是否上会了吗 +10 articlefan 2024-06-23 15/750 2024-06-25 16:05 by 请慎重修改昵称
[基金申请] 焦虑没有毛线用，默默前行是王道！ +4 漠上藜梭 2024-06-24 8/400 2024-06-25 14:32 by 漠上藜梭
[有机交流] 高温酯化反应喷料 20+4 hl24678 2024-06-21 7/350 2024-06-25 09:45 by ghj-hebust
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[有机交流] 三颈瓶加热 5+4 wzjwx888 2024-06-20 4/200 2024-06-24 10:11 by w87437985
[基金申请] 国自然青年基金，1A4B能上会吗？青年和面上的上会标准是一样的吗？ +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
[基金申请] F03青年基金函评结果 +5 暨阳一只柴 2024-06-19 6/300 2024-06-23 14:30 by adsqsj
[基金申请] 面上博士后 +11 jsqy 2024-06-19 13/650 2024-06-21 17:12 by sizhouyi
信息提示
关闭
请填处理意见
关闭