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

用戶注銷 (職業(yè)作家)
新蟲 (知名作家)

木蟲 (著名寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程085601調(diào)劑求老師收留 +9 | 強木木木 2026-03-07 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 283求調(diào)劑,材料、化工皆可 +5 | 蘇打水7777 2026-03-11 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | 18880831720 2026-03-11 | 5/250 |
|
|
[考研] 274求調(diào)劑0856材料化工 +9 | z2839474511 2026-03-11 | 10/500 |
|
|
[考研] 一志愿安徽大學材料工程專碩313分,求調(diào)劑的學校 +7 | Yu先生 2026-03-10 | 9/450 |
|
|
[考研] 工科0856專碩化學工程269能調(diào)劑嗎 +9 | 我想讀研11 2026-03-10 | 9/450 |
|
|
[考研] 泣血叩求調(diào)劑恩,愿以丹心報師恩 +5 | Iuruoh 2026-03-11 | 5/250 |
|
|
[考研] 考研材料與化工,求調(diào)劑 +7 | 戲精丹丹丹 2026-03-09 | 7/350 |
|
|
[考研] 420求調(diào)劑 +3 | 莫向外求11 2026-03-10 | 3/150 |
|
|
[考研] 307求調(diào)劑 +6 | 超級伊昂大王 2026-03-10 | 6/300 |
|
|
[考研] 一志愿天津大學,英一數(shù)二305分求調(diào)劑,四六級已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
[考研] 一志愿湖師大化學289求調(diào)劑 +5 | XMCMM3.14159 2026-03-10 | 5/250 |
|
|
[考研] 337一志愿華南理工材料求調(diào)劑 +5 | mysdl 2026-03-07 | 5/250 |
|
|
[考研] 求調(diào)劑 一志愿蘇州大學,0856化工323分 | 本科應化 | 有專利/競賽/科研助手經(jīng)歷 | +7 | 橙子cyx 2026-03-06 | 9/450 |
|
|
[考研] 085600材料與化工,一志愿廣州985,求調(diào)劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
[考研] 301求調(diào)劑 +5 | 一二LV 2026-03-07 | 5/250 |
|
|
[考研] 化學290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
[考研] 070300化學求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|