Equations diferentielle du 2è degre

Je voudrais savoir comment prg -t-on sous matlab pour resoudre les equation differentielles avec seconde membre avec la methode des elements finis ! qelq'un peut m'aider ?
Configuration: Windows 2000
Firefox 1.0.7

7 réponses

  1. bonjour,

    ODE45 est une méthode pour équation différentielles de Matlab.
    Une résolution Runge-Kutta ou approchant si je me souviens bien.
    Il faut donner les différents coefficients et leux variations ensuite apres plusieurs itérations, tu peux ploter le résultats.

    J'espère t'avoir aidé un peu.
    Bonne chance ;)
    0
    1. merci d'avoir m'aider mais est ce que tu peut m'expliqué plus ou me donner un algo ? car je suis vraiment un débutant .J'ai juste commancé avec Matlab en L3 mais pendant les TDs le prof nous a laissé nous débrouiller seuls , elle dit si vous n'avez pas comment ça marche allez sur help ou si vous ne comprenez pas l'anglais allez consulter un dictionnaire , voilà quoi !!!!Donc du coup je nne pouvait pas faire grande chose !
      0
    2. Alors, je devrais avoir un bout de code qui traine quelque part.
      Je vais essayer de te trouver ca.

      Laisse-moi deux ou trois jours histoire de partir en quete du fichier perdu.
      Et ca risque d'être une grande expédition :p

      Je te tiens au courant
      0
    3. @khunlunhaha merci Khulun bcp , j'ai aussi fait des recherche sur net et j'aimerais savoir comment on peut resoudre l'eq (1+x)Y''+Y'=-x en methodes des elements finis ? mais pas en methde Rung-Kutta !?
      0
    4. @whitechips01Il n'y a pas de Y(x) dans ton équation.
      A moins que je dises une grosse grosse bétise, tu peux faire un changement de variable, Y' = K

      Ton équation te donne (1+x)K' + K = -x
      Et tu te retrouve avec une equation différentielle de premier ordre.
      Tu trouve la solution (facile ca ^^) et pis tu intègre la solution trouvée (car k=Y' et pas Y).

      Pour la résolution, une méthode d'Euler toute bete et hop :)
      Le secret de l'univers sera dans tes mains.
      0
    5. @khunlunj'ai l'ai vu le pb mais puisque je voudrais uiliser le met des elements finis pour le resoudre mais pas du methode d'Euler ou Runge-Kutta ! Est ce tu connais cette methode ou bien methode de Ritz ?
      0
  2. bonjour;
    j'ai un probleme pour resoudre sous matlab cette equation différentielle:

    d²(d)/dt²+((4.7+0.75/(d(x)/dt+d(y)/dt)^0.5)*d(d)/dt=Fd(t)

    Avec :
    d : c’est x ou y
    Fd(t) : c’est 1 ou bien -1
    On donc 2 équations différentielle, l'une en fonction de ' x 'et l'autre en fonction de ' y ' a résoudre ?
    0
    1. aller au site : www.les-mathmathique.net
      0
      1. slt est merci pour votre aide j'ai besoin les étapes pour la création du programme d'une equation du 2 éme degré.
        0
        1. mon amie moi aussi je cherche queleque cours sur la mèthode des èlèment fins parceque mon projet de fin d'etude baser sur cette mèthode si je trouve quelque document sur cette mèthode je envoyer a toi merci et bonne chance.
          0
          1. salut , bonjour a tous, je veux résoudre une équation différentielle de premier ordre en matlab par l'utilisation de la fonction "ode45" , et malheureusement j'ai trouvé des problèmes sur cette fonction ,j'ai essayé de respecter tous les conditions de fonction de "ode45" et j'ai pas trouvé la faute ou exacte , alors sa me fais plaisir vraiment quant vous allez m'aider SVP et merci en tout cas.
            voici le programme que je le construire:
            %%%% LA FONCTION%%%%%%%%%
            function dphi=rsprime(t,phi)
            Qx=0.5;
            phi=pi/12;
            gamma1=0.08;
            kphi=3;
            dphi=gamma1+kphi*Qx
            %%%%%%%%%%%%%%%
            clc;
            clear;
            Qx=0.5;
            phi0=pi/12;
            gamma1=0.08;
            t0=0;
            tf=10;
            kphi=3;
            [t,phi]=ode45(rsprime,t0,tf,phi0);
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            et l'erreur qui je le trouve c'est
            " Error using ==> exist"
            The first input to exist is a string.
            et "error in ==> funfun\private\odearguments at 79"
            if (exist(ode)==2) % M-file
            et"Error in ==> ode45 at 173
            [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

            Error in ==> exemple2ode45 at 9
            [t,phi]=ode45(rsprime,t0,tf,phi0);
            aide moi SVP
            0
            1. Organigramme de runge kutta d ordre 4
              Données Y0, n, f(y,t)
              Dt=(t2-t1)/2n
              i = 0
              t0 = 0
              ((test)) i inférieur ou égal à n
              si non : écrire y0,....., yn
              stop
              si oui : calcul de y chapeau (inconnue) i+1/2 =.....
              calcul de y double chapeau =....
              Dt = (t2-t1)/n
              i = i+1/2
              CALCUL DE y chapeau i+1 =.....
              Dt/6 = (t2-t1)/6
              calcul de y i+1+....
              0