| 5 | 1/1 | 返回列表 |
| 查看: 2054 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
hollya至尊木蟲(chóng) (文壇精英)
葛城?{美
|
[交流]
從HPPC 抽出 電池參數(shù)的mat lab code 已有5人參與
|
||
|
前面發(fā)過(guò)一些擬合eis以便得到歐姆電阻,電容,反應(yīng)電阻的code:) 有童鞋說(shuō)沒(méi)有能測(cè)eis的儀器,但是又想知道以上幾個(gè)參數(shù)怎么辦吖?其實(shí)是有替代辦法的。比如,只要跑一下充電或放電的pulse(就是hppc那種),也可以用電池模型(非物理模型)擬合出來(lái)的啦。這里我就寫(xiě)了一個(gè)等效電路一階模型,用于擬合電流pulse得到上面三個(gè)參數(shù)。當(dāng)然,如果是二階的,也并沒(méi)有本質(zhì)上的復(fù)雜,只是多兩個(gè)參數(shù)。 這個(gè)方法比eis麻煩一些,但是對(duì)儀器要求比較低。所要的數(shù)據(jù)是三列矩陣,分別為時(shí)間,電流,電壓。 還有是就是要事先獲得soc-ocv關(guān)系,這個(gè)做個(gè)gitt也就有了。 好啦, good luck。。 ps: 小木蟲(chóng)需要部分重寫(xiě)了,adobe 已經(jīng)不再支持flash player。現(xiàn)在沒(méi)有辦法上圖:( 淚目。。。。 小木蟲(chóng),前向きに行ってみよう! 以下是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 |

新蟲(chóng) (小有名氣)
新蟲(chóng) (知名作家)

木蟲(chóng) (著名寫(xiě)手)

至尊木蟲(chóng) (文壇精英)
葛城?{美

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +8 | 小西笑嘻嘻 2026-03-03 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃?xún)|口芝士 2026-03-02 | 7/350 |
|
|
[考研] 成績(jī)276,專(zhuān)業(yè)代碼0856求調(diào)劑 +7 | 小陳朵 2026-03-03 | 7/350 |
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
[考研] 環(huán)境工程專(zhuān)碩307求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +4 | 詞凝Y 2026-03-02 | 4/200 |
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 材料化工調(diào)劑 +13 | 今夏不夏 2026-03-01 | 16/800 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 302材料工程求調(diào)劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[基金申請(qǐng)] 成果系統(tǒng)訪(fǎng)問(wèn)量大,請(qǐng)一小時(shí)后再?lài)L試。---NSFC啥時(shí)候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來(lái)了 2026-02-28 | 4/200 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-03-01 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-02-28 | 5/250 |
|