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

Fortran Discussion :

Equation différentielle


Sujet :

Fortran

  1. #1
    Membre à l'essai
    Inscrit en
    Février 2009
    Messages
    5
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 5
    Par défaut Equation différentielle
    salut,
    j'ai un projet de resoudre ce systeme d' equation differentielle couplé
    V'= Y'/V - mV ; Y''=1/V - exp (-Y) . par la methode de kunta runge et difference fini
    et la programmer en fortran, donc j'ai trouvé des resultas pour tracé les courbes V=f(x),Y=f(x)..
    N.B la dirivé dans les equations c'est par rapport à x.
    mais j'ai arrivé pas à tracé des courbes en fonction de m.
    je souhaite de m'aider a trouver comment je devrai faire pour tracer les courbes V=f(m),Y=f(m)..

    svp reponds a mon message
    merci bien.

  2. #2
    Rédacteur

    Homme Profil pro
    Comme retraité, des masses
    Inscrit en
    Avril 2007
    Messages
    2 978
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 84
    Localisation : Suisse

    Informations professionnelles :
    Activité : Comme retraité, des masses
    Secteur : Industrie

    Informations forums :
    Inscription : Avril 2007
    Messages : 2 978
    Par défaut
    Salut!
    1. De deux choses l'une: ou bien tu as un problème avec conditions initiales et tu utilises la méthode de Runge-Kutta, ou bien un problème avec conditions aux limites et tu utilises la méthode des différences finies; mais ça ne peut pas être les deux à la fois.
    2. Dans un cas comme dans l'autre, tu dois écrire tes équations sous forme canonique, c'est-à-dire sans dérivées à droite du signe =.
    3. Ton problème est du troisième ordre. Tu dois donc avoir trois variables d'état; deux d'entre elles sont V et Y, mais où est la troisième. La variable indépendante est x, mais que représente m: une variable ou un paramètre?

    Jean-Marc Blanc

  3. #3
    Membre à l'essai
    Inscrit en
    Février 2009
    Messages
    5
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 5
    Par défaut equation differentielle sous fortra
    salut;
    tout d'abord merci bien pour votre reponse,
    par exemple selon l'algorithme d'Euler:

    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
       DO i=1, n
     
    !   do  j=1,10
     Y(i+1) = Y(i)+dY 
    ! alpha=0.1
    ! alpha=alpha*j                    
     
    C Le systeme differentiel descriptif du systeme
          f1 =(X2(i)/X1(i))-(alpha*(X1(i)**(1+gama)))
         f3= X2(i)
     
          f2 = -exp(-X3(i))+(v0/X1(i))
    C L'algorithme d'EULER
            X1(i+1) = X1(i)+f1*dY                  ! Equation (1) du systeme
            X2(i+1) = X2(i)+f2*dY              ! Equation (2) du systeme
     
            X3(i+1)=X3(i)+f3*dY              !Equation (3) du systeme
    end do
      ! write(*,*)  REAL(Y(i)),REAL(X1(i))   
     
      " alpha c'est une parametre variable".   
     
     "cest un systeme d'equation monodimensionell suivant y"
     
     
    C On sauvegarde les données pour tracer la courbe X1(i) = f(Y) mais on pourrait tracer d'autre courbes,en particulier X2(i)=f(Y),,.
    MOI je voudrai calculer et mémoriser
    des valeurs qui serviront à tracer les courbes de X1 et X2 en fonction de de differente valeur de alpha par exemple " X3(i)=f(alpha)" ????
    comment je peurai introduire des different valeur pour alpha dans une boucle + la boucle spacial de y"
    ou bien si je voudrai fixé des valeurs pour X3(i) pour afficher des resultas de Y=f(alpha).
    J'attend votre reponse car je ss totalment perdu de ce coté.
    Cordialement.
    merci bien à vous.

  4. #4
    Rédacteur

    Homme Profil pro
    Comme retraité, des masses
    Inscrit en
    Avril 2007
    Messages
    2 978
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 84
    Localisation : Suisse

    Informations professionnelles :
    Activité : Comme retraité, des masses
    Secteur : Industrie

    Informations forums :
    Inscription : Avril 2007
    Messages : 2 978
    Par défaut
    J'attend votre reponse car je ss totalment perdu de ce coté.
    Moi aussi à la suite de ton second message; après le premier, je croyais encore y voir clair: il n'y a pas de relation visible entre les deux messages. Par exemple, où est passée la variable m? Tu ne réponds pas à mes questions.
    Jean-Marc Blanc

  5. #5
    Membre à l'essai
    Inscrit en
    Février 2009
    Messages
    5
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 5
    Par défaut
    salut;
    a oui c'est vrai, c'est la question de notation seulement ben m c'est alpha dans second message c'est meme equation.
    V'= Y'/V - alpha*V ; Y''=1/V - exp (-Y)
    dans 2 message voila selon l'algoritme d'auler,
    X1 C'est V
    X3 c'est Y

  6. #6
    Membre à l'essai
    Inscrit en
    Février 2009
    Messages
    5
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 5
    Par défaut
    le probleme est resolu y a pas de problemes des conditions.., le code ça marche bien, j'ai obtenu des resulta des inconnu V et Y en fonction de x pour valeur fixé de alpha.
    mon diffeculté c'est que j'ai arrivé pas à faire des resultas lorsque je fais varié alpha "c'est une parametre variable" on le donne des valeurs".donc je voudrai calculé les inconus comme V et Y en fonction de alpha.
    cordialement

Discussions similaires

  1. [DSP] Equation différentielle
    Par electro23 dans le forum Autres architectures
    Réponses: 7
    Dernier message: 06/05/2009, 23h50
  2. Equation différentielle non lineaire
    Par polal2is dans le forum MATLAB
    Réponses: 7
    Dernier message: 04/12/2008, 17h41
  3. comment trouvé les racines d'une equation différentielle
    Par samplaid dans le forum Mathématiques
    Réponses: 4
    Dernier message: 05/11/2007, 16h53
  4. Equation différentielle du second ordre
    Par moustiqu3 dans le forum MATLAB
    Réponses: 1
    Dernier message: 21/05/2007, 09h38

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