Scribd Logo
Importer

Bienvenue sur Scribd !

  • 13 vues

    Corrige Examen 2021

    Transféré par

    Ayoub Atifi
    1f ijmknv b b b b bj nn j. J nj. J. Jn jj njn jb bj bj. Nn bj. N. N.

    Droits d'auteur :

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

    Corrige Examen 2021

    Transféré par

    Ayoub Atifi
    13 vues4 pages
    1f ijmknv b b b b bj nn j. J nj. J. Jn jj njn jb bj bj. Nn bj. N. N.

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

    1f ijmknv b b b b bj nn j. J nj. J. Jn jj njn jb bj bj. Nn bj. N. N.

    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
    13 vues4 pages

    Corrige Examen 2021

    Transféré par

    Ayoub Atifi
    1f ijmknv b b b b bj nn j. J nj. J. Jn jj njn jb bj bj. Nn bj. N. N.

    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 4

    UNIVERSITE HASSAN II DE CASABLANCA ‫جامعة الحسن الثاني بالدار البيضاء‬

    --------- -------
    Faculté des Sciences Aïn Chock ‫كلية العلوم عين الشق‬

    Filière : Science des Matières Physique (SMP)


    Semestre 6 (S6)
    Corrigé Examen Projet Tutoré 1 : Analyse Numérique

    Exercice 1 (8 points)
    a ) Prenant les trois premiers points de collocation : (2 points)

    𝒕𝒊 V(𝒕𝒊 ) V[𝒕𝒊 , 𝒕𝒊+𝟏 ] V[𝒕𝒊 , 𝒕𝒊+𝟏 , 𝒕𝒊+𝟐 ]

    0 55 = 𝒂𝟎

    1 = 𝒂𝟏

    5 60 −𝟕⁄𝟓𝟎 = 𝒂𝟐

    −𝟐⁄𝟓

    10 58

    = −2 + 4𝑥
    𝑃2 (𝑡) = 𝑎0 + 𝑎1 (𝑡 − 𝑡0 ) + 𝑎2 (𝑡 − 𝑡0 )(𝑡 − 𝑡1 ) = 55 + 1(𝑡 − 0) + (−7⁄50)(𝑡 − 0)(𝑡 − 5)
    𝑃2 (𝑥) = 55 + 𝑡 + (−7⁄50)(𝑡 2 − 5𝑡) = 55 + (17⁄10)𝑡 + (−7⁄50)𝑡 2
    𝑃2 (2.5) = 58.375 𝑘𝑚⁄ℎ

    b ) Prenant les trois derniers points de collocation (2 points):

    𝒕𝒊 V(𝒕𝒊 ) V[𝒕𝒊 , 𝒕𝒊+𝟏 ] V[𝒕𝒊 , 𝒕𝒊+𝟏 , 𝒕𝒊+𝟐 ]

    35 57 = 𝒃𝟎

    -1 = 𝒃𝟏

    40 52 𝟏⁄𝟐𝟓 = 𝒃𝟐

    −𝟑⁄𝟓

    45 49

    = −2 + 4𝑥
    𝑄2 (𝑡) = 𝑏0 + 𝑏1 (𝑡 − 𝑡0 ) + 𝑏2 (𝑡 − 𝑡0 )(𝑡 − 𝑡1 ) = 57 + (−1)(𝑡 − 35) + (1⁄25)(𝑡 − 35)(𝑡 − 40)
    𝑄2 (𝑥) = 57 − 𝑡 + 35 + (1⁄25)(𝑡 2 − 75𝑡 + 1400) = 148 + −4𝑡 + (1⁄25)𝑡 2
    𝑄2 (42.5) = 50.25 𝑘𝑚⁄ℎ

    c) V = 69.248 Km/h pour t = 2.5 s et V = 27.02 Km/h pour t = 42.5 (1 point)


    Les valeurs accéptables sont celles obtenus aux questions a) et b), car :
    - Pour t = 2.5 les valeurs au-delà de 10 sont loin de 2.5 donc l’erreur d’interpolation augmente
    il en ets de même pour t = 42.5 les veleurs au-delà de 35 conduisent à une erreur
    d’interpolation plus grande.
    2 – Approximation d’ordre 8 de la distance 𝑑 (en métre ou en kilométre) parcourue pendant les
    40 premières secondes.(Algorithme de Romberg) (3 points):
    40
    𝑑 = ∫ 𝑉(𝑡)𝑑𝑡
    0
    Algorithme de Romberg :
    Prmière ligne méthode des trapèzes
    𝑏 𝑛−1
    ℎ (𝑏 − 𝑎)
    ∫ 𝑓(𝑥) ≊ [𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏)] , ℎ =
    𝑎 2 𝑛
    𝑖=1
    Deuxième ligne
    22 𝑇1,𝑖+1 − 𝑇1,𝑖
    𝑇2,𝑖 =
    22 − 1
    Troisième ligne
    24 𝑇2,𝑖+1 − 𝑇2,𝑖
    𝑇3,𝑖 =
    24 − 1
    Quatrième ligne
    26 𝑇3,𝑖+1 − 𝑇3,𝑖
    𝑇4,𝑖 =
    26 − 1
    Trapèze n = 1 Trapèze n = 2 Trapèze n = 4 Trapèze n =8 Ordre
    𝑇11 = 2140 𝑇12 = 2170 𝑇13 = 2205 𝑇14 = 2257.5 2
    = 594𝑚 = 602.7𝑚 = 612.5𝑚 = 627𝑚
    𝑇21 = 2180 𝑇22 = 2216.66 𝑇23 = 2275 4
    = 605.55 𝑚 = 615.74 = 631.94
    𝑇31 = 2219.104 𝑇32 = 2278.889 6
    = 616.417 𝑚 = 633.0248 𝑚
    𝑇41 = 2279.837 8
    = 633.288 𝑚
    Exercice 2 𝑓(𝑡, 𝑦) = 2𝑡 − 3𝑦 (7 points)
    1) h = 0.1 (2 points)

    n=1

    𝑦̂ = 𝑦0 + ℎ𝑓(𝑡0 , 𝑦0 ) = 1 + 0.1(2𝑥0 − 3𝑥1) = 0.7

    ℎ
    𝑦1 = 𝑦0 + [𝑓(𝑡0 , 𝑦0 ) + 𝑓(𝑡0 + ℎ, 𝑦̂)] = 1 + 0.05[(2𝑥0 − 3𝑥1 + 2𝑥0.1 − 3𝑥0.7)] = 0.755
    2

    𝑦(0.1) = −(2⁄9) + (2⁄3)0.1 + (11⁄9)𝑒−3𝑥0.1 = 0.749 888 936 4 ≅ 0.750

    𝑒𝑎 = 0.005

    n=2

    𝑦̂ = 𝑦1 + ℎ𝑓(𝑡1 , 𝑦1 ) = 0.755 + 0.1(2𝑥0.1 − 3𝑥0.755) = 0.5485

    ℎ
    𝑦2 = 𝑦1 + [𝑓(𝑡1 , 𝑦1 ) + 𝑓(𝑡1 + ℎ, 𝑦̂)] = 0.755 + 0.05[(2𝑥0.1 − 3𝑥0.755 + 2𝑥0.2 − 3𝑥0.5485)] = 0.589 475
    2

    𝑦(0.2) = −(2⁄9) + (2⁄3)0.2 + (11⁄9)𝑒−3𝑥0.2 = 0.581 880 888 6 ≅ 0.582

    𝑒𝑎 = 0.0076
    2) h = 0.05 (4 points)

    n=1

    𝑦̂ = 𝑦0 + ℎ𝑓(𝑡0 , 𝑦0 ) = 1 + 0.05(2𝑥0 − 3𝑥1) = 0.85

    ℎ
    𝑦1 = 𝑦0 + [𝑓(𝑡0 , 𝑦0 ) + 𝑓(𝑡0 + ℎ, 𝑦̂)] = 1 + 0.025[(2𝑥0 − 3𝑥1 + 2𝑥0.05 − 3𝑥0.85)] = 0.863 75
    2

    𝑦(0.05) = −(2⁄9) + (2⁄3)0.05 + (11⁄9)𝑒−3𝑥0.05 = 0.863 087 526 7 ≅ 0.8631

    𝑒𝑎 = 0.00066

    n=2

    𝑦̂ = 𝑦1 + ℎ𝑓(𝑡1 , 𝑦1 ) = 0.86375 + 0.05(2𝑥0.05 − 3𝑥0.86375) = 0.739 187 5

    ℎ
    𝑦2 = 𝑦1 + [𝑓(𝑡1 , 𝑦1 ) + 𝑓(𝑡1 + ℎ, 𝑦̂)] = 0.86375 + 0.025[(2𝑥0.05 − 3𝑥0.86375 + 2𝑥0.1 − 3𝑥0.7391875)]
    2
    = 0.751 029 687 5

    𝑦(0.1) = −(2⁄9) + (2⁄3)0.1 + (11⁄9)𝑒−3𝑥0.1 = 0.749 888 936 4 ≅ 0.750

    𝑒𝑎 = 0.001

    n=3

    𝑦̂ = 𝑦2 + ℎ𝑓(𝑡2 , 𝑦2 ) = 0.7510296875 + 0.05(2𝑥0.1 − 3𝑥0.7510296875) = 0.648 375 234 4

    ℎ
    𝑦3 = 𝑦2 + [𝑓(𝑡2 , 𝑦2 ) + 𝑓(𝑡2 + ℎ, 𝑦̂)]
    2
    = 0.751 029 687 5 + 0.025[(2𝑥0.1 − 3𝑥0.7510296875 + 2𝑥0.15 − 3𝑥0.6483752344]
    = 0.658 574 318 4

    𝑦(0.15) = −(2⁄9) + (2⁄3)0.15 + (11⁄9)𝑒−3𝑥0.15 = 0.657 101 074 2 ≅ 0.657

    𝑒𝑎 = 0.00157

    n =4

    𝑦̂ = 𝑦3 + ℎ𝑓(𝑡3 , 𝑦3 ) = 0.658 574 318 4 + 0.05(2𝑥0.15 − 3𝑥0.648 574318 4) = 0.566 288 170 6

    ℎ
    𝑦4 = 𝑦3 + [𝑓(𝑡3 , 𝑦3 ) + 𝑓(𝑡3 + ℎ, 𝑦̂)]
    2
    = 0.658 574 318 4 + 0.025[(2𝑥0.15 − 3𝑥0.658 574 318 4 + 2𝑥0.2 − 3𝑥0.566 288 170 6)]
    = 0.584 209 631 7

    𝑦(0.2) = −(2⁄9) + (2⁄3)0.2 + (11⁄9)𝑒−3𝑥0.2 = 0.581 880 888 6 ≅ 0.582

    𝑒𝑎 = 0.00233

    3) Extrapolation de RICHARDSON (1 point):

    22 𝑦4 − 𝑦2
    𝑦(0.2) = = 0.582 454 508 9 , 𝑒𝑎 = 0.00057
    22 − 1
    Exercice 3 (5 points)
    𝑉 ′ = 𝑔 − (𝑐 ⁄𝑚)𝑉(𝑡)
    {
    𝑉(0) = 0

    où la masse 𝑚 = 68 𝑘𝑔, le coefficient de résistance de l’air c = 14 kg/s et 𝑔 = 9.8 𝑚⁄𝑠 2 . On a résolu


    cette équation par la méthode de Range-Kutta d’ordre 4 et l’on a obtenu les résultats suivants :

    Solution
    𝑡 (𝑠) 𝑚
    𝑉 ( ⁄𝑠) 𝑡 (𝑠) 𝑉 (𝑚⁄𝑠)
    0 0.000 000
    1 8.855 846 6 33.760 030
    2 16.065 71 7 36.335 210
    3 21.933 240 8 38.431 230
    4 26.709 010 9 40.137 250
    5 30.596 160 10 ?
    a) Algorithme de Runge-Kutta ordre 4 (2 points)

    n = 10 pour calculer 𝑦10 , d’après le tableau on a h = 1

    𝑘1
    𝑘1 = ℎ𝑓(𝑡9 , 𝑦9 ) = 1𝑥(9.8 − (14⁄68)𝑥40.137250) = 1.536 448 529 4 , = 0.768 224 264 7
    2
    ℎ 𝑘1
    𝑘2 = ℎ𝑓 (𝑡9 + , 𝑦9 + ) = 1𝑥(9.8 − (14⁄68)𝑥40.905 474 426 7) = 1.378 284 676 9
    2 2
    𝑘2
    = 0.689 142 355 1
    2
    ℎ 𝑘2
    𝑘3 = ℎ𝑓 (𝑡9 + , 𝑦9 + ) = 1.394 566 279 8
    2 2
    𝑘4 = ℎ𝑓(𝑡9 + ℎ, 𝑦9 + 𝑘3 ) = 1.249 331 942 4

    1
    𝑦10 = 𝑦9 + (𝑘 + 2𝑘2 + 2𝑘3 + 𝑘4 ) = 41.525 830 397 5
    6 1
    b) Accélération du parachutiste à l’instant t = 4.5 s (1 point).

    𝒅𝑽 𝑽(𝒕+𝒉)−𝑽(𝒕−𝒉) 𝑽(𝟓)−𝑽(𝟒) 30.596 160−26.709 010


    𝑽′ = = = = = 𝟑. 𝟖𝟖𝟕 𝟏𝟓 𝒎⁄ 𝟐
    𝒅𝒕 𝟐𝒉 𝟐𝒙𝟎.𝟓 𝟐𝒉 𝒔
    c) Approximation d’ordre 4 (Simpson 1/3) de la distance verticale que le parachutiste parcourt
    pendant les 4 premières secondesAlgorithme de Simpson 1/3 (2 points) :

    𝑥𝑛 𝑛−1 𝑛−2
    ℎ
    ∫ 𝑓(𝑥)𝑑𝑥 = [𝑓(𝑥0 ) + 4 ∑ 𝑓(𝑥𝑖 ) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )]
    𝑥0 3
    𝑖=1,𝑖 𝐼𝑚𝑝𝑎𝑖𝑟𝑒 𝑖=2, 𝑖 𝑃𝑎𝑖𝑟𝑒

    4
    𝑑 = ∫0 𝑉(𝑡)𝑑𝑡 = 60.6656 𝑚.
    Approximation d’ordre 4 (Simpson 1/3) de la distance verticale que le parachutiste parcourt
    pendant les 8 premières soncondes.
    8
    𝑑 = ∫0 𝑉(𝑡)𝑑𝑡 = 194.127518 𝑚.

    Vous aimerez peut-être aussi