Discussion :

[romane333]

Bonjour bonjour,

Tout d'abord merci beaucoup à tous ceux qui vont répondre

J'ai ici une équation à résoudre pour sigma, que je n'arrive plus à simplifier, où les x_i sont des données d'un vecteur x (taille 100, donc j'ai pas vraiment envie de réécrire la formule à la main), et n est donné :



J'ai tenté d'utiliser solve:

Code :

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

solve('b+sum(x.*exp(-x/b))/sum(exp(-x/b))=2', 'b')

--> not a valid expression.
malheureusement, il me semble qu'écrire dans autre chose qu'une expression la division par un nombre inconnu est impossible...

Je suppose donc que solve n'est pas habilité à faire ce genre de calculs ? Y a-t-il une autre solution?

Merci! (encore )

(ma version matlab: 2011)

[procrastination]

Bonjour,

Code :

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

1
2
syms b
solve(b+(sum(x.*exp(-x/b))/sum(exp(-x/b)))==2, b)

Je te renvoie sur le liens du site mathworks spécifiant les paramètres de la fonction solve.
http://nl.mathworks.com/help/symbolic/solve.html


Mehdi

[romane333]

ça ne fonctionne pas non plus!

Je me suis mal exprimé, j'ai essayé tout, avec les syms ou avec des guillements autour de mes équations, dans tout les cas ça ne fonctionne pas. Matlab me renvoie "is not a valid expression or equation."


Est-ce que ce que vous m'avez envoyé fonctionne chez vous? c'est peut être un problème qui a été corrigé pour matlab 2014 mais ça ne fonctionne pas dans matlab 2011...

[Invité]

Bonsoir,

Déjà, vu que tu souhaites utiliser les valeurs de ton vecteur x, oublie les guillemets sinon le 'x' n'a rien à voir avec celui déjà défini.

Ceci fait, quelle erreur obtiens-tu sans avec ce que t'a fourni procrastination ?

Code :

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

1
2
syms b
solve(b+(sum(x.*exp(-x/b))/sum(exp(-x/b)))==2, b)

As-tu aussi pensé à la résoudre numériquement avec fzero ?

[romane333]

Bonsoir!

bon à savoir, pour les guillemets

Voici les erreurs avec la solution de procrastination, désolée de ne l'avoir pas bien écrit:

??? Error using ==> char
Conversion to char from logical is not possible.

Error in ==> solve>getEqns at 189
   vc = char(v);

Error in ==> solve at 67
[eqns,vars] = getEqns(varargin{:});

Sinon sur ton conseil j'ai utilisé fzero (pas bête!), mais ça me pose problème aussi:
- le fait de devoir mettre un guess c'est pas top, mais ça c'est un détail
- j'ai une fonction contenant un log, et la résolution de cette fonction avec fzero marche une fois sur cinq :/

Error using ==> fzero at 318
Function value at starting guess must be finite
and real.

voici la fonction à résoudre:

Code :

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

1
2
3
4
5
6
7
function y = muSol(mu, sample,n,beta)
    sum1=sum(exp(-sample/beta));
    y = -beta*log(sum1/n)-mu;
end
 
%appelée par:
 mu=fzero(@muSol, 10, [], sample,n,beta);

où sample est un vecteur, et n et beta des constantes.
Est-ce que tu saurais pourquoi ce problème?
(Je comprends si tu n'a pas envie de te pencher là-dessus, autant essayer d'utiliser plutot solve si c'est possible!)


Merci beaucoup de vos réponses!

[Invité]

Il me semble que sur les versions plus anciennes de MATLAB, il fallait rentrer une expression symbolique que l'on supposait nulle, donc ton cas on aurait donc :

Code :

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

solve(b+(sum(x.*exp(-x/b))/sum(exp(-x/b)))-2, b)

Pour ce qui est de fzero, sum1=sum(exp(-sample/beta)); doit être évalué à nul (ou tout du moins trop petit pour le type double). Ton log donne donc -Inf qui donne alors cette erreur.

[romane333]

Eh bien oui aha, c'était bien seulement ça le problème!
maintenant, ça fonctionne

et pour le -Inf, je vais faire avec..

Merci beaucoup!