WEEK 3 MATHEMATICAL MODELLING OF A BATTERY
T=(0:50:50000);##TIME VARYING FROM ZERO TO 50000 IN STEPS OF 50##
cr=zeros(1,1001);
dod=zeros(1,1001);##DEPTH OF DISCHARGE ##
v=zeros(1,1001);
cs=zeros(1,1001);
i=100;
v0=2.135;
noofcells=5;
capacity=50;
k=1.045;
deltat=50;
te=25;
for n=2:1001
peucap=((capacity/3)^k)*3;
cr(n)=cr(n-1)+((i^k*deltat)/3600);
dod(n)=cr(n)/peucap;
if dod(n)>0.99
dod(n)=1;
end
v(n)= noofcells*(v0-(10^-4)*(273+(te))*dod(n));
if dod(n)>0.99
dod(n)=1;
end
if v(n)>0
cs(n)=cs(n-1)+((i*deltat)/3600);
else
cs(n)=cs(n-1);
end
end
plot(cs,v,'b.')
axis([0.55 3.5 7]);
PLOT :
