Discussion :

[xelux]

Bonjour tout le monde,

J'aimerai tracer la courbe caractéristique (V,I) d'une cellule solaire modélisé par :



régit par l'équation implicite suivante reliant le courant :



Je considère connu tous les paramètres de cette équation et je désire trace la courbe (V,I) ...
J'ai essayer plusieurs methodes (fzeros(I-f(I,V)) ; newton's method ) sans succès , si vous pouvez m'aider.

Cordialement

[lucaswalker]

essaye ezplot
http://www.mathworks.it/it/help/matlab/ref/ezplot.html

[xelux]

ezplot aussi ne fonctionne pas non plus , voilà mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
%% numerical solving of an implicit equation 
 
T=300;
k=1.3806488 * 10^-23;
q=1.60217657 * 10^-19;
n=1;
I_0=3.049804310227623e-62;  
I_L=0.118;
Rs=90;
Rsh=3000;
 
syms y x
ezplot('y-I_L+I_0*exp(q*(x+y*Rs)/(n*k*T)-1)+(x+y*Rs)/Rsh)')

Et je recois le message d'erreur suivant :

cannot be plotted in the xy-plane.

[xelux]

J'ai essaié ce bout de code, mais le resultat est loin d'être satisfaisant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
T=300;
k=1.3806488 * 10^-23;
q=1.60217657 * 10^-19;
n=1;
I_0=3.049804310227623e-62;  
I_L=0.118;
Rs=90;
Rsh=3000;
 
annot=sprintf('y-%e+%e*exp(%e*(x+y*%e)/(%e*%e*%e)-1)+(x+y*%e)/%e',I_L,I_0,q,Rs,n,k,T,Rs,Rsh)
syms y x
ezplot(annot)

[le fab]

salut
ezplot permet de tracer une fonction, pas une chaine de caractères !

je vois que tu a la toolbox symbolic, alors tu peux faire truc comme :

Code :

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

1
2
3
4
5
a = 2;
b = 3;
syms x;
y = a*x+b;
ezplot(y);

sans la symbolic toolbox, cela donnerait :

Code :

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

1
2
3
4
a = 2;
b = 3;
f = @(x) a*x+b;
ezplot(f);

[xelux]

Tout d'abord je te remercie d'avoir essayé de m'aider :

Cependant j'ai deux remarque à faire, la première est que passer en argument une chaîne de caractère ne pose pas de problème, c'est d'ailleurs ce qui est proposé par la doc.


la deuxième est que vu que c'est une fonction implicite, je ne peux pas avoir y=f(x) mais plutot y=f(x,y) ce qui pose bien des soucis à titre d'information pour une résolution analytique et même numérique.

Par contre, si vous pouvez me proposer quelque chose de mieux adapter, je suis preneur !

[lucaswalker]

Bonjour,

Ce n'est pas temps le fait que la fonction soit implicite, mais plutôt qu'elle est difficile a résoudre.
En utilisant l'outil de calcul formel mupad, un résultat sort, et semble plus cohérent que celui de matlab. Sur quels intervalles V et I cherche tu la solution? Si a Mupad la syntaxe est la suivante:

Code :

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

1
2
3
4
5
6
7
8
9
10
T:=300;
k:=1.3806488 * 10^-23;
q:=1.60217657 * 10^-19;
n:=1;
I_0:=3.049804310227623e-62; 
I_L:=0.118;
Rs:=90;
Rsh:=3000;
f:=(x,y)->y-I_L+I_0*exp(q*(x+y*Rs)/(n*k*T)-1)+(x+y*Rs)/Rsh;
plot(f(x,y)=0)

[xelux]

Merci au moins là c'est décroissant mais normalement une pente dècroissante de vrait apparaitre au voisinnage des tensions elevés.
Je dois m'interesser aux valeurs positives du courant et de la tension, typiquement entre 0 et 0.5A pour le courant , et entre 0 et 100V pour la tension.