24小时热门版块排行榜    

查看: 2135  |  回复: 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
回复此楼
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
5楼: Originally posted by jm0 at 2021-03-13 10:06:01
安卓毫无压力。。。

安卓已经是小木虫最后的阵地

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
6楼2021-03-13 10:38:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

ncomm

新虫 (知名作家)

只为遇见你
2楼2021-03-12 13:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

烟圈o0

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
好久不见(虽然看不懂内容
未经省察的人生是不值一提的
3楼2021-03-12 16:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hollya

至尊木虫 (文坛精英)

葛城?{美

引用回帖:
3楼: Originally posted by 烟圈o0 at 2021-03-12 16:36:29
好久不见(虽然看不懂内容

哎,啥时候能发图,我再回水区

发自小木虫IOS客户端
でも私たちはきっと1000回もメールをやりとりして、たぶん心は1センチくらいしか近づけませんでした。
4楼2021-03-12 20:07:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 13/650 2026-03-15 12:20 by NSFC2026我来了
[考研] 化学调剂0703 +6 啊我我的 2026-03-11 6/300 2026-03-14 22:24 by shangxh
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
请填处理意见