Discussion :

[Esar001]

Bonjour à tous,
je voudrais utliser la fonction [fzero] pour trouver la racine d'une certaine fonction, voici 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
syms Ts lambda 
% Constantes
c = 3e8;
h = 6.6262e-34;
k = 1.38e-23;
% 
L = 1/(lambda^5*(exp(h*c/(lambda*k*Ts))-1));
%
dL = @(Ts) diff(L,'lambda');
% 
lambda = 600e-9;
% 
Ts =  fzero(dL, 5000);

Dans un premier temps je veux calculer la dérivée partielle de L (fonction de lambda et de Ts) par rapport à la variable lambda. J'obtiens donc dL et je veux ensuite grâce à [fzero] calculer, pour une valeur de lambda de 600e-9 fixée, la racine de ma fonction en partant d'une valeur arbitraire 5000. Ma variable est donc Ts ici.A l'éxécution, l'erreur suivante apparait :

??? Undefined function or method 'isfinite' for input arguments of type 'sym'.

Error in ==> fzero at 317
    elseif ~isfinite(fx) || ~isreal(fx)

Error in ==> Untitled4 at 13
Ts =  fzero(dL, 5000);

J'ai eu beau chercher pendant de longues heures, et essayer d'utiliser des fonctions comme [eval] ou [subs] je ne vois pas où est mon erreur.

Merci d'avance pour votre réponse

[Invité]

Bonjour,

Ton problème vient du fait que la fonction FSOLVE n'est pas faite pour des expressions symboliques mais pour des fonctions "normales".
Si tu souhaites garder des expressions symboliques, tu peux regarder la fonction SOLVE, sinon tu peux créer (dans le même style que tu fait déjà) une fonction anonyme du style

Code :

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

dL = @(Ts,lambda) ...

en évitant d'utiliser des variables symboliques.

[Esar001]

Bonsoir Winjerome,
je te remercie déjà pour m'avoir orienté dans ma recherche, néanmoins que ce soit par la première ou la deuxième méthode que tu me proposes j'obtiens toujours des erreurs, voici mon code :

pour la transformation en fonction symbolique:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
syms Ts lambda 
% Constantes
c = 3e8;
h = 6.6262e-34;
k = 1.38e-23;
% 
L = 1/(lambda^5*(exp(h*c/(lambda*k*Ts))-1));
%
dL =@(Ts,lambda)diff(L,'lambda');
% 
lambda = 600e-9;
%
Ts= fzero(dL,5000);

j'obtiens toujours la même erreur

??? Undefined function or method 'isfinite' for input arguments of type 'sym'.

Error in ==> fzero at 317
    elseif ~isfinite(fx) || ~isreal(fx)

Error in ==> Untitled4 at 15
Ts= fzero(dL,5000);
 

et si j'utilise à la place solve :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
syms Ts lambda 
% Constantes
c = 3e8;
h = 6.6262e-34;
k = 1.38e-23;
% 
L = 1/(lambda^5*(exp(h*c/(lambda*k*Ts))-1));
%
dL =@(Ts)diff(L,'lambda');
% 
lambda = 600e-9;
%
Ts =  solve('dL=0');

je n'ai pas d'erreurs mais dans la grille des variables (workspace) il n'y a aucune valeur attribuée à Ts, juste <1x1 sym>.
Pourrais-tu m'expliquer?

[Invité]

pour la transformation en fonction symbolique:

Justement non PAS SYMBOLIQUE, tu peux t'aider si tu veux des variables symboliques pour avoir l'expression de la dérivée, mais la fonction passée à FSOLVE doit être numérique.

Code :

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

1
2
3
4
5
6
7
8
9
10
<s>syms Ts lambda </s>
% Constantes
c = 3e8;
h = 6.6262e-34;
k = 1.38e-23;
% 
 
dL = @(Ts,lambda) <expression de ta dérivée en fonction de c,h,k,Ts et lambda>;
 
Ts= fzero(@(Ts) dL(Ts, 600e-9),5000);

et si j'utilise à la place solve :

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
syms Ts lambda 
% Constantes
c = 3e8;
h = 6.6262e-34;
k = 1.38e-23;
% 
L = 1/(lambda^5*(exp(h*c/(lambda*k*Ts))-1));
%
dL = @(Ts)diff(L,'lambda');
% 
lambda = 600e-9;
%
Ts =  solve(subs(dL)); 
% subs: pour substituer la valeur de lambda
% Sans parenthèses: sinon tu résous seulement x = 0

[Esar001]

Ok, j'ai compris, merci beaucoup et bonne soirée