Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Importer

Bienvenue sur Scribd !

  • 56 vues

    tp1 01

    Transféré par

    L'HungrY
    Ce document présente trois exercices sur la résolution numérique d'équations différentielles du premier ordre par des schémas explicites. Le premier exercice compare la solution exacte et approchée du modèle de Gompertz et d'une autre EDO. Le deuxième exercice étudie si une solution atteint un équilibre ou non. Le troisième exercice compare trois schémas explicites sur un exemple.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    tp1 01

    Transféré par

    L'HungrY
    56 vues2 pages
    Ce document présente trois exercices sur la résolution numérique d'équations différentielles du premier ordre par des schémas explicites. Le premier exercice compare la solution exacte et approchée du modèle de Gompertz et d'une autre EDO. Le deuxième exercice étudie si une solution atteint un équilibre ou non. Le troisième exercice compare trois schémas explicites sur un exemple.

    Titre original

    tp1-01

    Copyright

    © © All Rights Reserved

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce document présente trois exercices sur la résolution numérique d'équations différentielles du premier ordre par des schémas explicites. Le premier exercice compare la solution exacte et approchée du modèle de Gompertz et d'une autre EDO. Le deuxième exercice étudie si une solution atteint un équilibre ou non. Le troisième exercice compare trois schémas explicites sur un exemple.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    56 vues2 pages

    tp1 01

    Transféré par

    L'HungrY
    Ce document présente trois exercices sur la résolution numérique d'équations différentielles du premier ordre par des schémas explicites. Le premier exercice compare la solution exacte et approchée du modèle de Gompertz et d'une autre EDO. Le deuxième exercice étudie si une solution atteint un équilibre ou non. Le troisième exercice compare trois schémas explicites sur un exemple.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    sur 2

    U NIVERSITÉ D ’A IX -M ARSEILLE L3 de Mathématiques

    TP 1, équations différentielles du 1er ordre, schémas explicites

    Exercice 1 (Schéma d’Euler explicite).


    1. Modèle de Gompertz. On reprend ici le modèle de Gompertz (TD2, Exercice 1). Tracer sur un même graphe les
    solutions exactes et approchées par le schéma d’Euler explicite, avec différents pas de temps entre 0.1 et 1, de
    l’EDO  
    ′ K
    y (t) = y(t) ln
    y(t)
    pour K = 10, 0 ≤ t ≤ 10 et pour les conditions initiales suivantes
    (a) y(0) = 0.1,
    (b) y(0) = 2.
    On peut vérifier que les solutions sont globales et sont de la forme
     e−t
    y(0)
    y(t) = K .
    K

    Dans les deux cas demandés, que vaut limt→+∞ y(t) ?


    2. EDO étudiée en cours. Tracer sur un même graphe les solutions exactes et approchées par le schéma d’Euler
    explicite, avec différents pas de temps entre 0.1 et 0.001, de l’EDO y ′ (t) = (y(t))2 pour 0 ≤ t ≤ T et pour les
    conditions initiales suivantes
    (a) T = 2, y(0) = 0.2,
    (b) T = 0.49, y(0) = 2.
    Vérifier que les solutions sont, sur l’intervalle considéré de la forme
    1
    y(t) = 1 .
    y(0) −t

    Exercice 2 (Equilibre atteint ou non atteint).


    1. On étudie l’équation différentielle
    y ′ (t) = −y(t)2 , t ∈ IR + ,
    avec une condition initiale en t = 0 notée y0 . Dans la suite, on prendra y0 = 1 ou 2 et on s’intéressera à la
    solution pour t ∈ [0, T ], avec T = 1, 5 ou T = 10.
    (a) Vérifier (ou admettre) que les solutions sont bien de la forme :
    1
    y(t) = 1 .
    y0 +t

    (b) Représenter maintenant (pour différentes valeurs de y0 ) sur un même graphique la solution exacte et la
    solution approchée obtenue avec le schéma d’Euler explicite et un pas de temps de 0.01. Commenter.

    1
    2. On étudie maintenant l’équation p
    y ′ (t) = − |y(t)|, t ∈ IR + ,
    avec une condition initiale en t = 0 notée y0 . Dans la suite, on prendra aussi y0 = 1 ou 2 et on s’intéressera à la
    √
    solution pour t ∈ [0, T ], avec T = 2 y0 .
    √
    (a) Vérifier que les fonctions définies pour t ≤ 2 y0 par
    √ t 2
    y(t) = y0 −
    2
    sont bien des solutions.
    (b) Représenter maintenant (pour différentes valeurs de y0 ) sur un même graphique la solution excacte et
    la solution approchée obtenue avec le schéma d’Euler explicite et un pas de temps de 0.01. La solution
    approchée est elle strictement positive pour tout t<T ?

    Exercice 3 (Comparaison de trois schémas explicites).


    On s’intéresse ici au calcul pour t = 1 de la solution du problème

    y ′ (t) = y(t)(1 + e−y(t) ) + e2t ,


    y(0) = 0.

    1. Pour un pas de temps h = 1/n, n ≥ 10, programmer les trois méthodes suivantes :
    (a) Schéma d’Euler explicite,
    (b) Schéma de Heun,
    (c) Schéma de Runge-Kutta d’ordre 4 (noté RK4).
    2. En prenant comme valeur (presque) exacte la solution donnée par le schéma RK4 avec n = 1280, donner l’ordre
    de convergence de chacune des 3 méthodes considérées à la question 1 (en prenant n ≤ 640).
    3. Pour chacune des méthodes considérées à la question 1, en prenant des valeurs de n de la forme 2p 5, pour
    quelle valeur de n a-t-on la solution avec 2 décimales exactes ? pour quelle valeur de n a-t-on la solution avec 6
    décimales exactes ?

    Vous aimerez peut-être aussi