Scribd Logo
Importer

Bienvenue sur Scribd !

  • 354 vues

    Devoir Diff Finies

    Transféré par

    Lokmane xbox 360

    Droits d'auteur :

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

    Devoir Diff Finies

    Transféré par

    Lokmane xbox 360
    354 vues1 page

    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é ?

    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
    354 vues1 page

    Devoir Diff Finies

    Transféré par

    Lokmane xbox 360

    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
    Vous êtes sur la page 1sur 1

    Université de Versailles Saint-Quentin-en-Yvelines Année Universitaire 2013-2014

    Master 1 de Mathématiques EDP et Approximation Numérique

    Devoir N°1
    P. Gabriel

    TD2 : Différences finies

    Exercice 1 - Étude d’un schéma aux différences finies pour un problème 1D


    Soit f une fonction régulière. On souhaite approcher la solution du problème
    ( 0 (x)
    −u00 (x) + u1+x = f (x), 0 < x < 1,
    (P )
    u(0) = u(1) = 0.

    par la méthode des différences finies. Pour cela on considère une discrétisation (x0 , · · · , xN +1 ) de [0, 1]
    à pas h constant, ainsi que le schéma à trois points usuel pour la dérivée seconde :
    1
    (ui−1 − 2ui + ui+1 ) ' u00 (xi ), xi = ih.
    h2
    et le schéma centré suivant pour la dérivée première,
    1
    (ui+1 − ui−1 ) ' u0 (xi ).
    2h
    1. Ecrire le problème approché sous la forme d’un système linéaire Ah uh = bh en explicitant soi-
    gneusement la matrice Ah et le second membre bh . Quel est l’ordre de consistance du schéma ?
    2. Une matrice est dite "monotone" si elle est inversible et si sa matrice inverse est positive, au sens
    où tous ses coefficients sont positifs ou nuls. Montrer que la matrice Ah est monotone.
    3. On considère la fonction auxilliaire
    (1 + x)2 2
    θ : x ∈ [0, 1] 7→ − ln(1 + x) + (x2 + 2x) ln(2).
    2 3
    On note θi = θ(xi ), ∀i = 0, · · · , N + 1.
    (a) Montrer que θ est solution d’un problème aux limites de même type que (P ).
    (b) Montrer qu’il existe une constante C indépendante de h telle que

    θi+1 − 2θi + θi−1 θi+1 − θi−1
    − 1 6 Ch2 .

    max −

    2
    +
    16i6N h 2h(1 + ih)

    (c) On définit le vecteur Θ de composantes θi . Déduire de ce qui précède que

    (Ah Θ)i > 1 − Ch2 , ∀i = 1, · · · , N.

    4. (a) Montrer que pour la matrice Ah on a

    kA−1 −1
    h k∞ = kAh ek∞ ,

    où e est le vecteur dont toutes les composantes sont égales à 1.


    (b) Poser v = A−1
    h e et utiliser les résultats qui précèdent pour montrer qu’il existe une constante
    M indépendante de h telle que
    kA−1
    h k∞ 6 M.

    5. En précisant les hypothèses sur u solution de (P ), montrer que le schéma converge et préciser
    l’ordre de convergence obtenu.

    Vous aimerez peut-être aussi