| 查看: 2124 | 回复: 9 | |||
hollya至尊木虫 (文坛精英)
葛城?{美
|
[交流]
从HPPC 抽出 电池参数的mat lab code 已有5人参与
|
|
前面发过一些拟合eis以便得到欧姆电阻，电容，反应电阻的code：） 有童鞋说没有能测eis的仪器，但是又想知道以上几个参数怎么办吖？其实是有替代办法的。比如，只要跑一下充电或放电的pulse（就是hppc那种），也可以用电池模型（非物理模型）拟合出来的啦。这里我就写了一个等效电路一阶模型，用于拟合电流pulse得到上面三个参数。当然，如果是二阶的，也并没有本质上的复杂，只是多两个参数。 这个方法比eis麻烦一些，但是对仪器要求比较低。所要的数据是三列矩阵，分别为时间，电流，电压。 还有是就是要事先获得soc-ocv关系，这个做个gitt也就有了。 好啦， good luck！！！ ps： 小木虫需要部分重写了，adobe 已经不再支持flash player。现在没有办法上图：（ 泪目！！！！！ 小木虫，前向きに行ってみよう！ 以下是code： function [pbest]=rc1fit1(p0,data) % p0 in the sequence of r0,r1,c1 ydata= data(:,3) % voltage pbest=nlinfit(data,ydata,@rc1_v,p0) % use nlinfit() to fit plot(data(:,1),data(:,3),'o',data(:,1),rc1_v(pbest,data)) % plot xlabel('s') ylabel('v') legend('experiment','fit') end function ycal = rc1_v(p0,data) % the function returns terminal voltage (ycal) at time 't' load('sococveve') % soc_ocv relation soc=sococveve(:,1) % needs to specify for other cells ocv=sococveve(:,2) % needs to specify for other cells dt=diff(data(:,1)) i=data(:,2) capa=5 % need to be modified accordingly r0=p0(1) r1=p0(2) c1=p0(3) % this function assume r1 and c1 and r0 constant ocv(1)=data(1,3) % when t=0 soc(1)=interp1(ocv,soc,ocv(1))% return soc(1) at t=0 ycal(1)=ocv(1) % when t=0 ir1(1)=0 % when t=0, instantaneous current through paralleled r1 is 0 for ii=1:length(data(:,1))-1 ir1(ii+1)=exp(-dt(ii)/(r1*c1))*ir1(ii)+(1-exp(-dt(ii)/(r1*c1)))*i(ii) % dominant function soc(ii+1)=interp1(ocv,soc,ocv(ii))+dt(ii)*i(ii+1)/3600/capa ocv(ii+1)=interp1(soc,ocv,soc(ii+1)) ycal(ii+1)=ocv(ii+1)+ r0*i(ii+1)+ r1*ir1(ii+1) end ycal=ycal' end |
回复此楼» 本帖已获得的红花（最新10朵）
烟圈o0» 猜你喜欢
310 070300化学求调剂
已经有6人回复
-
福建师范大学 院士团队 招2026级博士1名
已经有10人回复
-
分析化学论文润色/翻译怎么收费?
已经有286人回复
-
广西师范大学化学与药学学院杨秀林课题组2026年拟招收博士研究生1名
已经有25人回复
-
26年储能、电池方向博士申请
已经有1人回复
-
求调剂，本科吉大，一志愿华工339
已经有2人回复
-
意大利博洛尼亚大学招聘研究助理2名（2D材料/锂离子与钠离子电池方向）
已经有1人回复
-
关于碳量子点抗菌的研究
已经有14人回复
-
南大化学调剂
已经有11人回复
-
硫酸钠溶液浓度与密度曲线
已经有0人回复
-
0703调剂
已经有2人回复
ncomm
新虫 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 3217.5
- 散金: 755
- 红花: 21
- 沙发: 9
- 帖子: 5142
- 在线: 447.7小时
- 虫号: 5263265
- 注册: 2016-11-21
- 性别: GG
- 专业: 其他无机非金属材料
2楼2021-03-12 13:51:45
烟圈o0
木虫 (著名写手)
- 应助: 178 (高中生)
- 金币: 3253.7
- 散金: 1035
- 红花: 45
- 沙发: 4
- 帖子: 2691
- 在线: 1502.6小时
- 虫号: 986567
- 注册: 2010-03-31
- 性别: GG
- 专业: 电化学
） 
3楼2021-03-12 16:36:29
hollya
至尊木虫 (文坛精英)
葛城?{美
- ECEPI: 1
- 应助: 744 (博后)
- 金币: 44165.3
- 散金: 1143
- 红花: 127
- 沙发: 780
- 帖子: 24910
- 在线: 835小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
4楼2021-03-12 20:07:20
匿名
用户注销 (职业作家)
- 应助: 10 (幼儿园)
- 金币: 1595.5
- 散金: 1198
- 沙发: 271
- 帖子: 3764
- 在线: 1697.9小时
- 虫号: 0
- 注册: 2020-07-12
- 性别: GG
- 专业: 藏学
 本帖仅楼主可见 |
5楼2021-03-13 10:06:01
hollya
至尊木虫 (文坛精英)
葛城?{美
- ECEPI: 1
- 应助: 744 (博后)
- 金币: 44165.3
- 散金: 1143
- 红花: 127
- 沙发: 780
- 帖子: 24910
- 在线: 835小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
6楼2021-03-13 10:38:37
7楼2021-06-18 08:58:05
8楼2021-06-18 08:58:12
9楼2024-11-13 15:45:00
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我司有EIS设备，可以测0-600Ah的电芯以及模组，还有AI大模型ECM拟合软件，以及DRT分析软件，欢迎大家咨询， v：zisheng977897 发自小木虫Android客户端 |
10楼2024-11-13 15:47:22
