function Ia = MSX60(Va,Suns,Tac)
%PVCELL Summary of this function goes here
% Detailed explanation goes here
k=1.38e-23;
q=1.6e-19;
A=1.2;
Vg=1.12;
Ns=36;
T1=273+25;
Voc_T1=21.06/Ns;
Isc_T1=3.8;
T2=273+75;
Voc_T2=17.05/Ns;
Isc_T2=3.92;
Tak=273+Tac;
TrK=273+25;
Iph_T1=Isc_T1*Suns;
a=(Isc_T2-Isc_T1)/Isc_T1*1/(T2-T1);
Iph=Iph_T1*(1+a*(Tak-T1));
Vt_T1=k*T1/q;
Ir_T1=Isc_T1/(exp(Voc_T1/(A*Vt_T1))-1);
Ir_T2=Isc_T2/(exp(Voc_T2/(A*Vt_T1))-1);
b=Vg*q/(A*k);
Ir=Ir_T1*(Tak/T1).^(3/A).*exp(-b.*((1./Tak)-(1./T1)));
X2v=Ir_T1/(A*Vt_T1)*exp(Voc_T1/(A*Vt_T1));
dvdi_Voc=-1.15/Ns/2;
Rs=-dvdi_Voc-1/X2v;
Vt_Ta=A*k*Tak/q;
Vc=Va/Ns;
Ia=zeros(size(Vc));
for j=1:5;
Ia=Ia-(Iph-Ia-Ir.*(exp((Vc+Ia.*Rs)./Vt_Ta)-1))./(-1-(Ir.*(exp((Vc+Ia.*Rs)./Vt_Ta)-1)).*Rs./Vt_Ta);
end
if Ia<=0
Ia=0;
end
评论2
最新资源