13.07.2015 Views

Design and Implementation of On-board Electrical Power ... - OUFTI-1

Design and Implementation of On-board Electrical Power ... - OUFTI-1

Design and Implementation of On-board Electrical Power ... - OUFTI-1

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

TaK = 273 + TaC;%Photo-current thermal coefficientK0 = 2.72e-4/Isc_Tref;%Photo-current (G=1 et T=Tref)Iph = Isc_Tref * G * (1 + K0*(TaK - Tref));%Diode saturation current (T=Tref)Id_Tref = Isc_Tref / (exp(q*Voc_Tref/(n*k*Tref))-1);%Diode saturation current (T=Tak)Id = Id_Tref * (TaK/Tref)^(3/n) * exp(-(Vg * q/(n*k))*(1/TaK - 1/Tref));%Calculation <strong>of</strong> serie resistance RsVpmax_Tref = 2.427; %voltage at maximum power for T=TrefIpmax_Tref = 0.529; %current at maximum power for T=TrefRs = (Voc_Tref-Vpmax_Tref)/Ipmax_Tref;%Iterative calculation <strong>of</strong> IscI = zeros(size(V));for i=1:10 \%number <strong>of</strong> iteration = 10for j=1:length(V);I(j) = I(j) - (Iph - I(j) - Id*(exp(q*(V(j)+I(j)*Rs)/(n*k*Tref)) -1))/...(-1 - (Id*(exp(q*(V(j)+I(j)*Rs)/(n*k*Tref)) -1))*q*Rs/(n*k*Tref));endendThe second function gives the voltage at the terminal <strong>of</strong> the used Schottky diode as afuntion <strong>of</strong> the current flowing in the diode.vmax = 3.5;function vv = v_diode(vi)% by Pierre Thirion% Use <strong>of</strong> function : V = v_diode(I)% V = Voltage on diode terminals [V]% V is computed by linear interpolation, on the basis <strong>of</strong> real measurements% <strong>of</strong> the I-V curve <strong>of</strong> the diode.% I = Current in the diode [A]% mesures (except the first)d_v = [0 0.141 0.213 0.227 0.24 0.25 0.26 0.27 0.280.3 0.32 0.34 0.36 0.38 0.39 0.4 0.41 0.42 0.44 0.460.48 0.5 0.52 0.54 0.56 0.58 ] ;d_i = [0 0.00001 0.00002 0.00003 0.0005 0.00007 0.0001 0.00017 0.000230.00047 0.00095 0.00206 0.0036 0.0095 0.015 0.021 0.031 0.04 0.072 0.1220.184 0.27 0.358 0.46 0.605 0.9 ];vv = zeros(size(vi));for k = 1:(length(vi)-1)104

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!