Code Matlab Cellule Photovoltaïque

Salut,

je suis en plein projet fin d'étude et je me trouve face a un problème de programmation matlab pour obtenir les caractéristiques I-V et P-V par la méthode de newton-raphson. j'ai beau essayé plusieurs code mais je n'y arrive pas. je pense que j'ai des problèmes dans le vecteur tension.

Je vous pris de bien vouloir m'aider s'il vous plait dans les délai les plus proches voici ma tentative du code matlab pour Ns = 36 cellule en series.

Va =12;
B = 0.7;
k = 1.38e-23; % Boltzman?s const
q = 1.60e-19; % charge of an electron
A = 1.2; % "diode quality" factor,
Eg = 1.12; % band gap voltage,
Ns = 36; % number of series connected cells (diodes)
T1 = 273 + 30;
Voc = 21.06/Ns; % open cct voltage per cell at temperature T1
Iscr = 3.80; % short cct current per cell at temp T1
Tr = 273 + 25; % reference temp
% when Va = 0, light generated current Iph_T1 = array short cct current
% constant "a" can be determined from Isc vs T
Iphr = Iscr;
a = 0.00065;
Iph = B*(Iphr + a*(Tr - T1));
Vt = k * T1 / q; % = A * kT/q
Vtr= k*Tr/q;
I0r = Iscr / (exp(Voc/(A*Vtr))-1);
b = Eg * q/(A*k);
Ir = I0r * (T1/Tr).^(3) .* exp(-b.*(1/Tr - 1/T1));
Rs = 0.3; % series resistance per cell
Vc = Va/Ns;
% Ia = 0:0.01:Iph;
for j=1:35
Ia(1) = Iscr;
Ia(j+1) = Ia(j) - (Iph - Ia(j) - Ir.*( exp((Vc+Ia(j).*Rs)./Vt) -1))./ (-1 - (Ir.*( exp((Vc+Ia(j).*Rs)./Vt) -1)).*Rs./Vt)

end
for Ia=0
Vc=Voc
end
for Ia=Iscr
Vc=0
end
figure, hold on
plot(Vc,Ia)
for y=2:34
Vc(y+1)=Vc(y);
end
title('Photovoltaic Module I-V curve')
xlabel('Module Voltage (V)')
ylabel('Module current (I)')
axis([0 50 0 5])
figure, hold on
for Tr= 1:10:81
Pa=Vx.*Ia;
plot(Vx,Pa)
end
title('Photovoltaic Module P-V Curve')
xlabel('Module Voltage')
ylabel('Module Power')
axis([0 90 0 90])

S'il vous plait c'est urgent
Merci d'avance