Discussion :

[tpscience]

Bonjour à tous,

J'écris donc en ce moment un programme en langage matlab, et j'aurais besoin de résoudre une équation. Pour cela j'utilise actuellement la fonction "fzero", or un problème majeur inhérent à cette dernière est qu'elle nécessite de connaitre à l'avance autour de quelle valeur on cherche notre zéro, ce qui n'est pas mon cas ici !!
Je pensais donc passer à la fionction "solve", or avec celle-là, le problème vient du fait que j'utilise des paramètres dans mon équation, et les solutions qui apparaissent sont des solutions avec les paramètres, soit avec leur expression littérale et non numérique ! Comment pourrais-je garder mes paramètres tout en leur assignant une valeur pour les calculs. En gros, pour être clair rien ne vaut un petit exemple :

Si j'écris cela :

Code :

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

1
2
3
4
5
 
syms x ;
a=sym(2);
b=double(a);
solve('b*x + 1')

Matlab me renvoie :

ans =

-1/b

Or, vous l'aurez maintenant compris, j'attendais -1/2...!!!!!!

Si quelqu'un avait une solution pour un (voire les deux) de mes problèmes...!

Merci par avance.

[Jean Dumoncel]

Bonjour,

il faut utiliser la fonction subs, avec ton exemple :

Code :

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

solve(subs('b*x + 1','b',double(a)))

[tpscience]

Bonjour,

Avant tout, merci beaucoup pour la réponse, en effet cela fonctionne très bien !
En revanche, j'obtiens plusieurs solutions, or je ne voudrais garder que les solutions réelles et positives. Pour cela j'impose donc ma solution comme telle :

Code :

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

syms x real positive;

Or, le paramètre

Code :

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

real

semble bien fonctionner alors que le

Code :

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

positive

non. J'ai l'impression que cela est dû au fait que la fonction que je résous comprend des fonctions trigos, car quand je fais la même simulation avec de "simples" fonctions (style polynomiales) le paramètre

Code :

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

positive

est bien pris en compte...???!!

Merci par avance.

[Jean Dumoncel]

Nulle part dans la doc il est indiqué que l'on peut passer plusieurs paramètres (réel, positif,...) à syms.

De plus :

Code :

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

syms x  positive

est équivalent à :

Code :

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

sym('x', 'positive');

et d'après la doc :

sym('x', 'positive') creates the symbolic variable x and assumes that x is real and positive.

[tpscience]

En effet,

Code :

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

positive

implique bien réel et positif à la fois...!
En revanche, je ne m'explique toujours pas pourquoi ce n'est pas pris en compte avec des fonctions trigos...?

Merci encore.

[Jean Dumoncel]

Aurais-tu un exemple?

[tpscience]

Oui biensur, mais ce n'est rien d'exceptionnel. En gros, en imposant :

Code :

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

x = sym('x','positive')

et en prenant comme fonction :

Code :

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

cos(x) - sqrt(2)/2;

j'obtiens 2 solutions, positive et négative.
Alors qu'avec la fonction :

Code :

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

x^2 - 1

j'ai bien la solution négative filtrée...

[Jean Dumoncel]

Effectivement, et je ne vois pas comment faire à part tester le signe des solutions...

[tpscience]

Bonjour,

En effet, c'est ce que j'avais fais en attendant une meilleure solution...

Encore Merci