IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
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

MATLAB Discussion :

Problème avec fsolve


Sujet :

MATLAB

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Nouveau candidat au Club
    Profil pro
    Inscrit en
    Décembre 2007
    Messages
    1
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Décembre 2007
    Messages : 1
    Par défaut Problème avec fsolve
    Bonjour à tous


    Navré de vous importuner avec mon petit problème, j'espère que pour vous ça sera plus lumineux que pour moi

    Je dois résoudre un système 5x5 (5 eq & 5 inconnus) faiblement non-linéaire (2 seulement équations le sont) ; le solveur d'Excel est hélas inefficace car les résultats ne sont pas dans le bon intervalle de recherche. Je suis donc passé sous Matlab où j'espère m'en sortir avec la fonction Fsolve

    Pour ceux qui sont curieux le problème à résoudre provient des équations de la combustion de l'ethane suivant un volume d'air variant : outre la bonne vieille réaction de combustion il y a en plus deux réactions de dissociation des molécules (ce sont ces deux équations qui sont non-linéaires) ; l'objectif est de déterminer les compositions molaires ce qui permettra de remonter à la température

    Les équations sont les suivantes :



    "e" représente l'excès (ou le défaut) de dioxygène ; ici on prendra e = -0.5 ; toutes les solutions doivent être positives (ou nulles)

    Voici pour l'instant mon programme

    %Resolution des equations de combustion

    % x1 : x CO2
    % x2 : y CO
    % x3 : w H2O
    % x4 : v H2
    % x5 : q O2
    % i : i N2

    %Donnees
    P= 1;
    Kp4= 0.000194017;
    Kp8= 5.398915192;

    e= -0.5

    i= 18.8*(1+e) % mole de diazote

    F= inline('[x(1)+x(2)-2 ; x(3)+x(4)-3 ; 2*x(1)+x(2)+x(3)+2*x(5)-10*(1+e) ; x(2).^2*x(5)*P./(x(1).^2*(x(1)+x(2)+x(3)+x(4)+x(5)+i)-Kp4; x(2)*x(3)./x(1)-Kp8]','x');

    [x, fval] = fsolve(F, [0;0], optimset('fsolve'))
    Evidemment ça ne marche pas, et je n'ai pas trouvé de solutions grace à l'aide (sinon je ne viendrais pas vous ennuyer avec ça :pfff: )

    Si l'un ou plusieurs d'entre vous arrive à faire fonctionner tout ça merci beaucoup de me montrer ; évidemment si il y en a qui savent utiliser autre chose que "Fsolve" je prends aussi ^^

    Merci d'avance !!

  2. #2
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 317
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 317
    Par défaut
    J'essaierai comme ceci :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    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
    26
    27
    28
    29
    30
    31
    32
    33
    function combustion
     
    %Resolution des equations de combustion
     
    % x1 : x CO2
    % x2 : y CO
    % x3 : w H2O
    % x4 : v H2
    % x5 : q O2
    % i : i N2
     
    % F= inline('[x(1)+x(2)-2 ; x(3)+x(4)-3 ; 2*x(1)+x(2)+x(3)+2*x(5)-10*(1+e) ; x(2).^2*x(5)*P./(x(1).^2*(x(1)+x(2)+x(3)+x(4)+x(5)+i)-Kp4; x(2)*x(3)./x(1)-Kp8]','x');
     
    x0=[0;0;0;0;0];
     
    options=optimset('Display','iter'); 
    [x, fval] = fsolve(@myfun,x0,options)
     
    function F=myfun(x)
    %Donnees
    P= 1;
    Kp4= 0.000194017;
    Kp8= 5.398915192;
     
    e= -0.5;
     
    i= 18.8*(1+e); % mole de diazote
     
    F=[x(1)+x(2)-2 ; 
        x(3)+x(4)-3 ;
        2*x(1)+x(2)+x(3)+2*x(5)-10*(1+e) ; 
        x(2).^2*x(5)*P./(x(1).^2*(x(1)+x(2)+x(3)+x(4)+x(5)+i))-Kp4 ; 
        x(2)*x(3)./x(1)-Kp8 ];
    En changeant bien entendu les valeurs de départ de x0

Discussions similaires

  1. problème avec fsolve
    Par lilettre dans le forum MATLAB
    Réponses: 3
    Dernier message: 22/09/2012, 18h26
  2. Problème avec la fonction fsolve
    Par VegaLyra dans le forum Scilab
    Réponses: 2
    Dernier message: 01/03/2011, 13h59
  3. Problème avec le type 'Corba::Any_out'
    Par Steven dans le forum CORBA
    Réponses: 2
    Dernier message: 14/07/2002, 18h48
  4. Problème avec la mémoire virtuelle
    Par Anonymous dans le forum CORBA
    Réponses: 13
    Dernier message: 16/04/2002, 16h10

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo