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

Bienvenue sur Scribd !

  • 5 vues

    SMA-S5 Projet3

    Transféré par

    groupestheorie
    0pmp0

    Droits d'auteur :

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

    SMA-S5 Projet3

    Transféré par

    groupestheorie
    5 vues2 pages
    0pmp0

    Titre original

    SMA-S5_Projet3

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

    0pmp0

    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
    5 vues2 pages

    SMA-S5 Projet3

    Transféré par

    groupestheorie
    0pmp0

    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

    UNIVERSITE CADI AYYAD Année universitaire 2018-2019

    Faculté des Sciences Semlalia Filère SMA Semestre S5


    Département de Mathématiques Automne 2018
    Analyse Numérique II
    Projet de Travaux Pratiques n◦ 3

    Écrire les scripts des fonctions dans des chiers : · · · .sci, et les tests dans des chiers : · · · .sce
    Mettre votre Nom en commentaire à l'en-tête de chaque chier. Respecter la syntaxe du nom
    des fonctions et variables demandées et commenter vos résultats.
    Travail à faire sur le compte rendu :
    On considère l'équation diérentielle suivante :
    −u”(x) = 1, pour x ∈]0, 1[
    (
    (1)
    u(0) = u(1) = 0

    1. Montrer que (1) admet une solution unique qui peut s'écrire sous la forme intégral suivante :
    Z 1
    u(x) = G(x, t)dt
    0
    où G est une fonction à décrire.
    2. On eectue une subdivision régulière de l'intervalle [0, 1] en N intervalles [xi , xi+1 ] où xi = Ni .
    Montrer que pour une valeur approchée de la solution u au point xi notée ui , on a :
    −ui−1 + 2ui − ui+1 = N −2 pour 1 ≤ i < N
    (

    avec u0 = uN = 0

    3. Écrire la matrice A de ce système linéaire et son second membre B . Quelle est sa dimension?
    Est-elle symétrique dénie positive? justier.
    4. Réécrire et répondre aux questions précédentes dans le cas général:
    −u”(x) = f (x), pour x ∈]a, b[
    (
    (2)
    u(a) = u(b) = c
    où a, b et c sont dans IR, avec b > a.
    5. Calculer la solution exacte u(x) pour les cas suivants:
    (a) f (x) = 1., [a, b] = [0, 1] et c = 0.
    (b) f (x) = sin(π x−a
    b−a
    ), [a, b] = [0, 2] et c = 0.

    (c) f (x) = cos(4π x−a


    b−a
    ), [a, b] = [−4, 4] et c = 1.

    (d) f (x) = 1 + cos(4π x−a


    b−a
    ), [a, b] = [0, 4] et c = 2.

    1
    Scripts des fonctions : chier : Fonctions.sci

    Dans un chier "Fonctions.sci", écrire les fonctions :


    • [y]=Fun(x,a,b) qui correspond à la fonction f (x) du problème (2),
    • [z]=uExact(x,a,b,c) qui correspond à la solution exacte du problème (2)
    calculée dans la question 5.
    où x, y et z sont des vecteurs colonnes.
    Script des diérences nies : chier: DifFinies.sci

    Dans un chier "DifFinies.sci", écrire une fonction [A, B]=DifFinies(a, b, c, N, Fun)


    programmant la discrétisation par diérence nie de (2) avec :
    1. Les arguments d'appel sont :
    • [a, b] : l'intervalle [a, b],
    • c : la condition aux bords c,
    • N : le nombre de subdivisions N ,
    • Fun : le nom de la fonction f (x) de (2).

    2. Les arguments de sortie sont :


    • A : la matrice carrée de discrétisation A,
    • B : le vecteur colonne du second membre B .

    Test n◦ 1 : chier : test1.sce

    Pour des input : a, b, c, N .


    • Écrire le script "test1.sce" pour tracer le graphique de uExa(.) solution exacte et la solution
    approchée uApp(.) obtenue par diérences nies aux points x = [x1 = a; x1 ; ...; xN +1 = b]
    pour les diérentes valeurs de N = 10 et N = 50.
    • Ajouter au graphique un titre et une légende et exporter le graphique vers un chier image.
    • Acher la valeur de ErDif f = kuExa − uAppk∞ pour chaque valeur de N donnée par
    l'utilisateur.
    • Noter les valeurs d'entrées et de sorties dans un tableau sur le compte rendu et commenter
    vos résultats.
    D'autres précisions peuvent être ajoutées au moment du TP.

    Test n◦ 2 : chier : test2.sce

    D'autres testes pour d'autres fonctions f peuvent être ajoutées au moment du TP.

    ◦
    Projet de Travaux Pratiques n 3 SMA-S5 2018

    Vous aimerez peut-être aussi