Discussion :

[powersolda]

Salut tout le monde.

S'il vous plait, aidez moi pour tracer cette figure :
Soit trois vecteurs :
i : vecteur de courant de GPV.
v: vecteur de tension de GPV.
E: vecteur de variation d'éctairment.

j'ai tracé cette figure en 3D mais ce n'est pas comme je veux (avec mesh comme illustré sur l'image attachée)

le lien de l'image
http://www.4shared.com/photo/auX2yzPm/Sans_titre_2.html
merci a l'avance.

[Jean Dumoncel]

Bonjour,

il n'y a pas d'image attaché à ton message.

Peux-tu nous montrer le code utilisé?

[powersolda]

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
for E=100:100:1000
 
    k=1;i=3.8;Pmax=[];pmax=0;ii=3.8;
    while ii>=0
 
        y=fsolve('tens',y0,[],Isc,Io,Vth,Rs,i);
        v=real(y(1));
        DT=T-Tr;
        DI=alpha*i*(E/Er)*DT+((E/Er)-1)*Isc;
        DV=-beta*v*DT-Rs*DI;
        ii=i+DI;
        ipv(k)=ii;
        vv=v+DV;
        vpv(k)=vv;
        y0=[ii];
        p(k)=vv*ii;
        if p(k)>pmax
            vmpp=[vpv(k)];
            impp=[ipv(k)];
            pmpp=[p(k)];
            pmax=p(k);
            Pmax=[pmax p];
        end
 
        TT(k) = E;
        k=k+1;
        i=i-0.05;
 
    end
    a = meshgrid(vpv);
    b= meshgrid(ipv);
    c = meshgrid(TT);
    figure(3)
    mesh(c,a,b)
    axis([0 1000 0 25 0 4]);
    figure(1)
    plot3(vpv,TT,ipv);axis([0 25 0 1000 0 4]);
    hold on
    title('caractéristque I-V dun GPV')
    xlabel('La tension (en V)')
    ylabel('Eclairement(en W/m2)')
    zlabel('Le courant (en A)')
    plot3(vmpp,E,impp,'r*');
    grid on
    figure(2)
    plot3(vpv,TT,p);axis([0 25 0 1000 0 62]);
    hold on
    title('caractéristque P-V dun GPV')
    xlabel('La tension (en V)')
    ylabel('Eclairement(en W/m2)')
    zlabel('La puissance (en W)')
 
    plot3(vmpp,E,pmpp,'r*');
    grid on
end

[Jerome Briot]

Je ne suis pas sûr mais essaie en replaçant :

Code :

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

1
2
3
    a = meshgrid(vpv);
    b= meshgrid(ipv);
    c = meshgrid(TT);

par

Code :

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

    [a,b,c] = meshgrid(vpv, ipv, TT);

De plus, certains crochets sont inutiles dans ton code :

Code :

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

1
2
3
4
5
6
y0=[ii];
        p(k)=vv*ii;
        if p(k)>pmax
            vmpp=[vpv(k)];
            impp=[ipv(k)];
            pmpp=[p(k)];

Et pour finir, seuls les graphiques construits lors de la dernière itération sur E seront affichés.

[powersolda]

salut merci pour votre aide
mais il reste le mème problème aussi avec votre aidé