1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
| %%Constantes
Icc=8.362;
Is=1.1e-9;
q=1.6e-19;
Rs=9;
a=1.714;
K=1.38e-23;
T=298;
Rp=690;
eps=0.0001;
%%parametres de l'algorithme
I=0;
V=0;
pas=0.01;
%vecteur tension
for i=1:200; %calclul des valeurs de la tension
f=I;
f(i)=-I+Icc-Is*(exp(q*(V+Rs*I)/a*K*T)-1)-(V+Rs*I/Rp);
fp(i)=-1-Is*(q*Rs/a*K*T)*(exp(q*(V+Rs*I)/a*K*T))-(Rs/Rp);
I=I-f(i)/fp(i);
if (abs(I-f))<eps;
break
end
courant(i)=I;
tension(i)=V;
V=V+pas;
end
figure
plot(tension,courant)
title('Adjudted I-V curve');
xlabel('V(V)');
ylabel('I(A)'); |
Partager