Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

  1. #1
    Nouveau membre du Club
    Résolution algébrique d'un système d'équations
    Bonjour à tous!

    Voila j'ai les deux éqautions suivantes:

    A=f+d
    B=f-3d


    avec 2 inconnues f et d.

    Je voudrais que Matlab puisse trouver les inconnues en fonction de A et B de telle sorte que j'obtienne la solution suivante:

    d=(A-B)/4
    f=A-(A-B)/4


    Dois-je passer par le calcul matriciel?
    J'avais pour idée d'utiliser la fonction "evala" mais je n'ai pas trouver les arguments pour obtenir la solution!

    Merci de votre aide!!!!

    JB

  2. #2
    Rédacteur/Modérateur

    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  3. #3
    Nouveau membre du Club
    J'ai regarder le site. Cela marche pour l'exemple que j'avais donner précedement. Mais pas pour celui que je veux résoudre:

    K = solve('m=g/v-h*x/v-y*i/w','n=-x*g/v+h/v-y*i/w','0=-z*g/v-z*h/v+i/w','r=(1+x)*l/v','0=y-z*w/v');
    K.v
    K.w
    K.x
    K.y
    K.z

    j'obtiens l'erreur suivante:

    Warning: Explicit solution could not be found.
    > In solve at 140
    In Stiff at 11
    ??? Access to an object's fields is only permitted within its methods.

    Error in ==> Stiff at 12
    K.v
    Pourtant j'ai 5 équations:

    -----------------------------
    m=g/v-h*x/v-y*i/w

    n=-x*g/v+h/v-y*i/w

    0=-z*g/v-z*h/v+i/w

    r=(1+x)*l/v

    0=y-z*w/v
    ----------------------------

    et je veux connaitre les 5 inconnues v, w, x, y, z en fonction de m, n, r, g, h, i, et l.


    Est-il possoble de le résoudre ce problème,
    Merci d'avance!

    JB

  4. #4
    Rédacteur/Modérateur

    Que vaut K à la sortie de SOLVE ?
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  5. #5
    Membre régulier
    Tu veux absolument résoudre ton problème en symbolique ?
    Sinon tu peux le faire en numérique, comme ca tu auras pour une solution (fk,dk) pour chaque couple (Ak,Bk).

  6. #6
    Nouveau membre du Club
    Citation Envoyé par Dut Voir le message
    Que vaut K à la sortie de SOLVE ?
    K est vide !

    Citation Envoyé par nant44 Voir le message
    Tu veux absolument résoudre ton problème en symbolique ?
    Sinon tu peux le faire en numérique, comme ca tu auras pour une solution (fk,dk) pour chaque couple (Ak,Bk).
    Je voudrais bien le faire en symbolique: cela me simplifierai très largement la tache pour les futur calculs!

    -------------------------

    Je viens aussi de récrire le systeme et il y avait une équation qui était redondante. Le nouveau systeme est:

    K = solve('m=g/v-h*x/v-y*i/w','n=-x*g/v+h/v-y*i/w','0=-y*g/w-y*h/w+i/w','r=(1+x)*l/v');
    K.v
    K.w
    K.x
    K.y

    Donc 4 équations:
    -------------------------
    m=g/v-h*x/v-y*i/w

    n=-x*g/v+h/v-y*i/w

    0=-y*g/w-y*h/w+i/w

    r=(1+x)*l/v

    -------------------------

    et 4 inconnues : v, w, x et y (à exprimer en fonction de g, h, i, l, m, n et r)

    J'obtiens toujours malheureusement la meme erreur que précedement...

    Savez-vous d'ou peut provenir cette erreur?

    Merci d'avance!!!!

    JB

  7. #7
    Rédacteur/Modérateur

    Je ne possède pas cette Toolbox mais as-tu essayé d'utiliser SOLVE comme ceci ?

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    [v,w,x,y] = solve('m=g/v-h*x/v-y*i/w','n=-x*g/v+h/v-y*i/w','0=-y*g/w-y*h/w+i/w','r=(1+x)*l/v');
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  8. #8
    Nouveau membre du Club
    J'obtiens l'erreur suivante:

    Warning: Explicit solution could not be found.
    > In solve at 140
    In Stiff at 22


    A noter:
    si au lieu de la premiere equation (m=g/v-h*x/v-y*i/w) j'ai:
    m=-g/v-h*x/v-y*i/w

    (- devant le premier terme)

    Matlab arrive a trouver une solution!!!

    Cela veut-il dire que pour mon systeme aucune solution existe????

    Merci d'avance!

    JB

  9. #9
    Expert confirmé
    en effet cela veut dire que matlab et maple ne trouvent aucune solution analytique comme on peut le vérifier en faisant également appel directement au noyau maple:
    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    maple('solve({-m+g/v-h*x/v-y*i/w,-n-x*g/v+h/v-y*i/w,-y*g/w-y*h/w+i/w,r=(1+x)*l/v},{v,w,x,y})')

    cela retourne une chaine vide.
    MATLAB 8.2 (R2013b)

    le blog des MATLAB geeks

  10. #10
    Nouveau membre du Club
    Merci pour la réponse meme si ca ne m'arange pas beaucoup!
    Il est quand meme sans doute possible de résoudre numériquement cette equation, non?

  11. #11
    Rédacteur/Modérateur

    Citation Envoyé par jbollagnier Voir le message
    Il est quand meme sans doute possible de résoudre numériquement cette equation, non?
    Essaie peut être de poser ta question sur le forum Algorithmes > Math en mentionnant cette discussion.

    Certain "spécialistes" sauront peut être mieux t'orienter vers des méthodes de résolutions...
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  12. #12
    Nouveau membre du Club
    Merci pour tout!

  13. #13
    Nouveau membre du Club
    Bonjour,

    Me revoila!

    J'ai un nouveau systeme:

    [v,w,x,y] = solve('m=g*v-h*y*w-x*z*v','n=-y*g*w+h*w-y*z*w','0=-x*g*v-y*h*w+z*v','r=(w+w*w*y/v)*l')

    avec toujours 4 équations

    ---------------------------
    m=g*v-h*y*w-x*z*v'

    n=-y*g*w+h*w-y*z*w

    0=-x*g*v-y*h*w+z*v

    r=(w+w*w*y/v)*l

    ------------------------------

    et 4 inconnues v, w, x et y.


    Si j'essaie de le résoudre avec matlab j'obtiens une nouvelle erreur:

    ??? Error using ==> solve
    Unable to find closed form solution.

    Error in ==> Stiff at 3
    [v,w,x,y] = solve('m=g*v-h*y*w-x*z*v','n=-y*g*w+h*w-y*z*w','0=-x*g*v-y*h*w+z*v','r=(w+w*w*y/v)*l')
    J'aimerais cependant trouver une solution numérique (g, h, i, l, m, n et r peuvent etre remplacés numériquement).

    Comment puis-je faire ?

    Merci!

    JB

  14. #14
    Expert confirmé
    Un autre utilisateur a eu ce problème sur le ML Central, jette un oeil, ça pourrait peut etre t'aider.
    MATLAB 8.2 (R2013b)

    le blog des MATLAB geeks

  15. #15
    Nouveau membre du Club
    Merci j'ai regardé et j'ai vu aussi un conseil dans l'aide: fsolve (optimization Toolbox) (http://www.mathworks.com/access/help...uche&meta=)

    J'ai essaye d'écrire un programme semblable à l'aide:


    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
     
    function F = fun(x) 
    F = [-2.8*10^-3-7.21*10^6*x(1)-1.19*10^9*x(4)*x(2)-x(3)*2.36*10^8*x(1);
        -1.13*10^-3-x(4)*-7.21*10^6*x(2)+1.19*10^9*x(2)-x(4)*2.36*10^8*x(2);
        -x(3)*-7.21*10^6*x(1)-x(4)*1.19*10^9*x(2)+2.36*10^8*x(1);
        4.38*10^-6+(x(2)+x(2)*x(2)*x(4)/x(1))*1.89*10^6];
     
     
    %Next, call an optimization routine. 
     
    x0 = [-5; -5];           % Make a starting guess at the solution
    options=optimset('Display','iter');   % Option to display output
    [x,fval] = fsolve(@fun,x0,options)  % Call optimizer


    [QUOTE]??? Input argument "x" is undefined.

    Error in ==> Stiff at 2
    F = [-2.8*10^-3-7.21*10^6*x(1)-1.19*10^9*x(4)*x(2)-x(3)*2.36*10^8*x(1);[/QUOTE]



    D'ou provient cette erreur?

    J'aimerais aussi pouvoir donner des valeurs de départ différentes à x(1), x(2), x(3) et x(4)!

    Quelqu'un a-t-il des idées?

    Merci!

    JB

  16. #16
    Expert confirmé
    Pour que ton code marche il faut que x0 doit ait 4 valeurs vu que dans fun tu as x(1),x(2),x(3) et x(4).
    En prenant par exemple
    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    x0=[-5;1;2;3]

    ton code fonctionne sans problème
    MATLAB 8.2 (R2013b)

    le blog des MATLAB geeks

  17. #17
    Nouveau membre du Club
    Oups !
    Je m'en suis apercu mais j'obtiens toujours la meme erreur !

    Code :Sélectionner tout -Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
     
    function F = fun(x) 
    F = [-2.8*10^-3-7.21*10^6*x(1)-1.19*10^9*x(4)*x(2)-x(3)*2.36*10^8*x(1);
        -1.13*10^-3-x(4)*-7.21*10^6*x(2)+1.19*10^9*x(2)-x(4)*2.36*10^8*x(2);
        -x(3)*-7.21*10^6*x(1)-x(4)*1.19*10^9*x(2)+2.36*10^8*x(1);
        4.38*10^-6+(x(2)+x(2)*x(2)*x(4)/x(1))*1.89*10^6];
     
     
    %Next, call an optimization routine. 
     
    x0 = [0.3;0.3;100e9;100e9];           % Make a starting guess at the solution
    options=optimset('Display','iter');   % Option to display output
    [x,fval] = fsolve(@fun,x0,options)  % Call optimizer



    Erreur:
    [QUOTE]??? Input argument "x" is undefined.

    Error in ==> Stiff at 2
    F = [-2.8*10^-3-7.21*10^6*x(1)-1.19*10^9*x(4)*x(2)-x(3)*2.36*10^8*x(1);[/QUOTE]

    Ais-je fais une erreur dans la syntaxe?

    Merci!

    JB

  18. #18
    Rédacteur/Modérateur

    Citation Envoyé par jbollagnier Voir le message
    Ais-je fais une erreur dans la syntaxe?
    Comment fais-tu appel à la fonction fun ?
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  19. #19
    Nouveau membre du Club
    J'ai ecris le programme tel quel.
    Sans une ligne de plus ou de moins !

    JB

  20. #20
    Rédacteur/Modérateur

    Et où est-ce que tu définis le paramètre d'entrée x de la fonction fun ?
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et Arduino)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)