[fsolve] résolution d'une équation polynomiale

Bonjour,

J'ai un problème itératif dont je dois résoudre une équation non linéaire, et qui dépend des paramètres et des variables de mon problème.

j'ai définie dans un fichier .m (version R2016a de MATLAB) la fonction:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 for i=1:n+1 if i==1 v(i)=0; r(i)=0; theta(i)=e/(l-a); Ro(i)=1/(l-b).*theta(i); w(i)=(((A+C)*w_r)/(A*Ra+C.*(Rs_0-Ro(i)))).*Ro(i); Nc(i)=Cste_3*m.*Ro(i)*w(i).^2; d(i)=2*pi*(Rs_0+r(i)); else v(i)=2*pi*kv/H.*Nc(i-1).*d(i-1); % r(i) solution de pi/3*(l-a)*(3*Rr_0*r(i)+r^2(i))=v(i) r(i)=fsolve(pi/3*(l-a)*(3*Rr_0*r(i)+r^2(i))-v(i)=0,0); theta(i)=(e+r(i))/(l-a); Ro(i)=1/(l-b).*theta(i); w(i)=(((A+C)*w_r)/(A*Ra+C.*(Rs-Ro(i)))).*Ro(i); Nc(i)=Cste_3*m.*Ro(i)*w(i).^2; d(i)=2*pi*(Rs_0+r(i)); end end
mais ce programme ne marche pas (un problème avec la fonction fsolve), en effet je sais pas comment je résous une telle équation et intégrer la solution dans le programme.
merci d'avance.

résolution d'une équation polynomiale

Bonjour,

plutôt :

Code:

r(i)=fsolve(pi/3*(l-a)*(3*Rr_0*r(i-1)+r(i-1)^2)-v(i),0);

Citation:

Envoyé par phryte

Bonjour,

plutôt :

Code:

r(i)=fsolve(pi/3*(l-a)*(3*Rr_0*r(i-1)+r(i-1)^2)-v(i),0);

Bonjour,

j'ai essayé ce code, et il m'a affiché cet erreur pour la fonction fsolve:
Code:

1 2 3 4 5 Error using lsqfcnchk (line 108) FUN must be a function, a valid string expression, or an inline function object. Error in fsolve (line 210) funfcn = lsqfcnchk(FUN,'fsolve',length(varargin),funValCheck,gradflag);
je pense qu'il y a un problème avec la fonction fsolve. que pensez vous?

salut

j'avoue que je comprend pas trop ce que tu veux faire, cad de quelle fonction (fonction de quelle(s) variable(x)) tu cherches une solution
phryte t'as surtout indiqué une erreur de syntaxe évidente dans l'écriture r^2(i)) qui ne veut rien dire

sinon pour utiliser fsolve correctement il faut écrire quelque chose comme :

Code:

Xsolution = fsolve(@(x) f(x, parametres(i)), x0);

je te laisse écrire correctement la fonction, soit directement dans le solve, soit dans un fichier séparé
(je te conseille la deuxième option car elle est "moins piégeuse" et que tu as pas l'aire de maîtriser le sujet)

Fabien

1 pièce(s) jointe(s)

Citation:

Envoyé par le fab

pour utiliser fsolve correctement il faut écrire quelque chose comme :

Code:

Xsolution = fsolve(@(x) f(x, parametres(i)), x0);

Fabien

Salut,

effectivement je maitrise pas MATLAB et le problème que j'ai, est combiné.
Voici une image descriptive du problème itératif:
Pièce jointe 319754

pour ce la, j'ai le code MATLAB suivant, sauf que il marche pas, j'ai encore un problème avec la fonction fsolve, j'ai essayé de faire la première solution que vous m'avez proposer:
Code:

1 2 f=@(x)(pi/3*(l-a)*(3*Rr_0*x+x^2)-v(i)); r(i) = fsolve(@(x) f(x,(i-1)), 0);
que pensez vous? comment puisse je résoudre ce problème?

plutôt quelque chose comme :
Code:

1 2 f=@(x,i)(pi/3*(l-a)*(3*Rr_0*x+x^2)-v(i)); r(i) = fsolve(@(x) f(x,(i-1)), 0);

Je vous remercie pour votre aide.