Discussion :

[bilou_12]

Bonjour tout le monde,

Dans le code ci-dessous:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
        p = polyfit(abs, ord, degre_polyfit);
 
        syms 'x';
        %f(x) = p(1)*x^6 + p(2)*x^5 + p(3)*x^4 + p(4)*x^3 + p(5)*x^2 + p(6)*x^+ p(7); 
 
        dp=zeros(7);
        dp=polyder(p);
 
        eq = dp(1)*x^5 + dp(2)*x^4 + dp(3)*x^3 + dp(4)*x^2 + dp(5)*x + dp(6);
 
        solutions = solve(eq);
 
        if size(solutions,1)>0 && solutions(1,1)> 0

je cherche à résoudre l'équation

eq = 0

c'est à dire connaitre pour quels x l'eq s'annule.

Les dp() sont des "doubles". J'ai ajouté le mot clé syms pour créer une variable symbolique.

Pourtant, à la dernière ligne:

Code :

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

if size(solutions,1)>0 && solutions(1,1)> 0

, le code bugge.

Est-ce que quelqu'un pourrait m'expliquer comment régler ce problème svp?

Merci d'avance!

[Invité]

Bonsoir,

Il est toujours utile d'avoir le message d'erreur...
Je suppose toutefois qu'il s'agit de

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

solutions étant de type syms, tu ne peux pas effectuer une comparaison numérique dessus. Il te faut le repasser en type double :

Code :

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

subs(solutions(1,1))

[bilou_12]

Excuse-moi pour ne pas avoir écrit le message d'erreur.

Mais bien vu car c'est celui là!

Je vais tester ton bout de code! Merci beaucoup!!!