Discussion :

[nzouclau]

Bonjour à tous.

svp pourriez-vous m'aider à corriger mon programme.

ci-joint code et algorithme.

Cordialement


le code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
% Doctorant : Claude, 
% Directeur de Thèse : Pr Martin & Pr Patrice 
 
% Cet algorithme calcul les paramètres PV d'une cellule solaire
 
% INITIALISATION DES DONNEES EXPERIMENTALES MESUREES
clear;
clc;
 
% Chargement des données expérimentales mesures 
% Tempereature de la surface du module PV
% Irradiance solaire 
donnees=load('extract_T');
% Courant expérimentales mesures 
I_mes=0.1:0.5:7;
% Tension expérimentales mesures 
V_mes=6:1:20;
 
% Parametres du fabriquant sur le module PV
% Nombre de cellule en serie
Nsc=36;
% Constante de bolztman
k=1.38*10^-23; 
% Irradiance aux conditions STC
Gref=1000; 
% Temperature aux onditions STC
Tref=25+237.15;
% Charge de l'électron
q=1.60*10^-19; 
% Coefficient de tension
Kv=-0.0865; 
% Courant de circuit au condition STC
Iref=4.01; 
% Nombre de mesures ou donnees expérimentales enregistrees
Nbre_mesure=length(donnees.T);
 
% On parcourt toutes les mesures enregistrées (mes ou theo) grace a j
 
for j=0:Nbre_mesure
 
    % Initialisation des parametres de depart
 
    % Valeur minimale admise
    % Courant de saturation de depart ou valeur minimale admise
    Is(j)=10^-10;
    % Facteur de qualité de depart  ou valeur minimale admise
    n(j)=1.0; 
    % Resistance serie de depart ou valeur minimale admise
    Rs(j)=0.01;
    % Resistance parallele de depart  ou valeur minimale admise
    Rp(j)=0.01; 
 
    % Pas d'encremantation
    % Courant de saturation à chaque incrementation
    delta_Is=0.001;
    % Facteur de qualité à chaque incrementation
    delta_n=0.001;
    % Résistance serie à chaque incrementation
    delta_Rs=10^-9;
    % Résistance parallele à chaque incrementation
    delta_Rp=0.01;
 
    % Valeur maximale admise
    % Courant de saturation max
    Is_max=7;
    % Facteur de qualité max
    n_max=2; 
    % Resistance serie max
    Rs_max=2;
    % Resistance parallele max
    Rp_max=200;
 
 
    % Pour chaque enregistrement calculons Is, n, Rs et Rp en utilisant N_it
    N_it=1;
 
    % On initialise epsilon qui est la difference entre la valeur mesuree et celle theorique obtenue grace a lambert
    epsilon(j)=5;
 
    % Tension thermique de la jonction pn a la temperature T
    Vt(j)=(k*donnees.T(j))/q;
 
    % Calcul du courant photogenere a la temperature T sous une irradiance G
    Iph(j)=((donnees.G(j))/(Gref))*(Iref+Kv*(donnees.T(j)-Tref))
 
 
 
 
 
    while(N_it<=200 && Rs(j)<Rs_max && Is(j)<Is_max && n(j)<n_max && Rp(j)<Rp_max && epsilon(j)>10^-2)  
 
        %expression du courant theorique Ipv_theo a comparer au courant mésuré I_mes
 
 
        Ipv_theo(j)=(((Is(j)+Iph(j))*Rp(j))-V_mes(j))/(Rp(j)+Rs(j))-(((n(j)*Nsc*Vt(j))/Rs(j))*lambertw((Rs(j)*Rp(j)*Is(j))/(n(j)*Nsc*Vt(j)*(Rs(j)+Rp(j)))*exp(((Rs(j)*Rp(j))*((V_mes(j)/Rs(j))+Iph(j)+Is(j)))/((n(j)*Nsc*Vt(j))*(Rs(j)+Rp(j)))))); 
        % On calcul le courant théorique à partir des précedentes  valeurs de Rs, Is et n calculés plus haut
 
        %Ipv_theo(j)=Is(j)+Iph(j)-((n(j)*Nsc*Vt(j))/(Rs(j)))*lambertw(((Rs(j)*Is(j))/(n(j)*Nsc*Vt(j)))*exp((V_mes(j)+Rs(j)*Iph(j)+Rs(j)*Is(j))/(n(j)*Nsc*Vt(j))));  
        % On calcul le courant théorique à partir des précedentes  valeurs de Rs, Is et n calculés plus haut 
 
        %on calcul l'erreur entre I du module mesure experimentalement au I %theorique
        epsilon(j)=abs(I_mes(j)-Ipv_theo(j));  
 
 
        % Incrementation des different parametres du module pour rapprocher %I mesure a I theo
        % on augmente la précedente résistance série de delta_Rs=0.001  
        Rs(j)=Rs(j)+delta_Rs;
        % on augmente le précedente courant de saturation de delta_Rs=10^-9
        Is(j)=Is(j)+delta_Is;
        % on augmente le précedente facteur de qualité delta_Rs=0.001
        n(j)=n(j)+delta_n;  
        % on augmente le précedente facteur de qualité delta_Rs=10
        Rp(j)=Rp(j)+delta_Rp;
        % On incremente le parametre d'une valeur delta constante
        N_it=N_it+1;
    end
 
    %recuperation de la valeur de Iph
    Iph
    %recuperation de la valeur de Is
    Is
    %recuperation de la valeur de n
    n
    %recuperation de la valeur de Rs
    Rs
    %recuperation de la valeur de Rp
    Rp
end