close all; clear all; clc;
% influence de Facteur de qualit� sur Les caract�ristiques I(V)
k=1.38e-23; q=1.60e-19; Vg=1.12; Nc=25; Rsh=25;
E=1000; En=1000; T=273+25; Vth =k*T/q; S1=0.3132; S=S1; Ns=1; Np=1;
T1=273+25; Vco1=21/Nc; Icc1=3.80; T2=273+75; Vco2=17/Nc; Icc2=3.92;
%------------
for A =1.25:0.5:2.25;
Iph=Icc1*(E/En)*(S/S1)*(1+((Icc2-Icc1)/Icc1*1/(T2-T1)*(T-T1)));
Is=(Icc1/(exp(Vco1/(A*Vth*T1/T))-1))*(T/T1).^(3/A).*exp(-(Vg*q/(A*k)).*(1./T-1/T1));
Rs=(1.15/2)/Nc -1/((Icc1/(exp(Vco1/(A*Vth*T1/T))-1))/((A*Vth)*T1/T)*exp(Vco1/(A*Vth*T1/T)));
%...................
% method de Newton:
%...................
% I=Iph-Is*(exp((V+I.*Rs)./(A*Vth))-1)-(V+I*Rs)/Rsh);
% f(I)=Iph-I-Is*(exp((V+I*Rs)/(A*Vth))-1)-(V+I*Rs)/Rsh =0;
% I=I-f(I)/f'(I)
V=0:0.1:25;
I=zeros(size(V/Nc));
if I<0; I=0
end
I=I-(Np*Iph-I-Np*Is*(exp((V/(Nc*Ns)+(Rs/Np)*I)./(A*Vth))-1)-(V/(Nc*Ns)+(Rs/Np)*I)*(Np/Rsh))./(-1-(Np*Is)*(((Rs/Np)/(A*Vth))*exp((V/(Nc*Ns)+(Rs/Np)*I)./(A*Vth))-1)-(Rs/Rsh));
p=I.*V;
max(p)
%........................
title('caracteristique P(V)')
xlabel('tension (V)')
ylabel('puissance (W)')
grid;
axis([0 22 0 70])
%------------
hold on
plot(V,p);
end