Equation non linéaire

Bonjour je ne sais pas si je trouverai quelqu'un sur ce site qui puisse m'aider car je souhaite coder sous MATLAB...
Le titre est plutôt clair je souhaite résoudre un systéme de 4 équations non linéaires. J'ai un bon début de code mais j'ai un message d'erreur que je ne comprends pas du tout...

Voici le code :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

function F=myfun(x) V1=400000 %V2=x(2) %P1=x(1) P2=100000000 %Q1=x(4) Q2=100000000 R=1.575 X=22.275 Z=sqrt(R*R+X*X) teta1=0 %teta2=x(3) d=atan(R/X) dprime=((pi/2)-d) Y=0.00014361   F=[(x(1))-(V1*x(2)/Z)*(sin(teta1-(x(3))-dprime))+(V1*V1/Z)*(sin(dprime)); P2-(V1*x(2)/Z)*(sin(x(3)-teta1-dprime))+(x(2)*x(2)/Z)*(sin(dprime)); x(4)-(V1*x(2)/Z)*(cos(teta1-x(3)-dprime))+(V1*V1)*((cos(dprime)/Z)-Y); Q2-(V1*x(2)/Z)*(cos(x(3)-teta1-dprime))+(x(2)*x(2))*((cos(dprime)/Z)-Y)]     x0=[-5;-5;-5;-5] options=optimset('Display','iter'); [x,fval]=fsolve(@myfun,x0,options,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y)

---

Et voici le message d'erreur:

??? Input argument "x" is undefined. 

Error in ==> myfun at 23 
F=[(x(1))-(V1*x(2)/Z)*(sin(teta1-(x(3))-dprime))+(V1*V1/Z)*(sin(dprime));

HELP !!

Merci d'avance.

Bonjour,

Il faut faire attention que l'appel à la fonction fsolve et la définition de ta fonction myfun doivent se faire dans deux fonctions distinctes.
De la façon dont tu effectues l'appel:

Code:

---

[x,fval]=fsolve(@myfun,x0,options,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y)

---

Tu dois définir ta fonction et tes différentes variables comme ceci:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

function resolution_equation V1=400000 %V2=x(2) %P1=x(1) P2=100000000 %Q1=x(4) Q2=100000000 R=1.575 X=22.275 Z=sqrt(R*R+X*X) teta1=0 %teta2=x(3) d=atan(R/X) dprime=((pi/2)-d) Y=0.00014361   x0=[-5;-5;-5;-5]; options=optimset('Display','iter'); [x,fval]=fsolve(@myfun,x0,options,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y) % ou [x,fval]=fsolve(@(x) myfun(x,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y),x0,options)

---

Puis

Code:

---

1 2 3 4 5 6

function F=myfun(x,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y)   F=[(x(1))-(V1*x(2)/Z)*(sin(teta1-(x(3))-dprime))+(V1*V1/Z)*(sin(dprime)); P2-(V1*x(2)/Z)*(sin(x(3)-teta1-dprime))+(x(2)*x(2)/Z)*(sin(dprime)); x(4)-(V1*x(2)/Z)*(cos(teta1-x(3)-dprime))+(V1*V1)*((cos(dprime)/Z)-Y); Q2-(V1*x(2)/Z)*(cos(x(3)-teta1-dprime))+(x(2)*x(2))*((cos(dprime)/Z)-Y)];

---

PS: Pense à mettre des ';' à la fin de tes lignes, ce qui empêche l'affichage de la valeur dans la ligne de commande.

Merci ça marche,

Puis-je poser une question de plus?

Peux-tu m'expliquer ce que l'ont doit renseigner dans les valeur de x0.
Je me rends compte que c'est super important pour la résolution et je n'ai pas compris entièrement son utilité.

Merci !!

La valeur x0 est le point de départ de la résolution, c'est à partir de celle-ci que l'algorithme commence sa convergence vers une solution. Et selon sa valeur, il se peut qu'il converge vers une "mauvaise" solution (Exemple: minimum local et non global), ou n'en trouve pas du tout.
Donc, si tu as une idée de la valeur à trouver, il est bon de fixer une valeur x0 approchante.

Juste pour être certains d'avoir bien compris.
Le premier terme du tableau de x0 je présume qu'il correspond à x(1) et ainsi de suite ?

Puis-je prendre des valeurs différentes pour chacune de mes inconnues parce que je n'arrive pas à faire en sorte que mon vecteur f(x) soit proche de zéro avec un vecteur de gradient stable.

Encore merci,

Je dis oui à tes deux propositions :) x0 correspond en tous points à x, et tu peux très bien lui mettre n'importe quelles valeurs, aussi différentes quelles soient.

Ok,
Mais j'ai l'impression de pas faire converger la résolution...

Je dois mettre dans les valeurs de x(1) les valeurs les plus proches du résultats comment fait on ?

Comme je te l'ai dit, si tu as une idée des valeurs finales rentre des valeurs proches. Mais si tu n'en as aucune, il n'y a pas vraiment de valeurs miracles :aie:, essaie avec plusieurs... à tatons :P

Merci beaucoup pour ton aide.

:ccool:

Ce forum est hyper réactif !! :mrgreen: