Scribd Logo
Importer

Bienvenue sur Scribd !

  • 129 vues

    Devoir 14 Analyse Numerique Corrige

    Transféré par

    idk04319
    Ce document présente des exercices sur le calcul numérique. Le premier exercice détaille le calcul de l'intégrale d'une fonction sur un intervalle en utilisant la méthode des trapèzes. Le deuxième exercice examine la convergence de la méthode de Newton-Raphson pour trouver les zéros d'une fonction.

    Droits d'auteur :

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

    Devoir 14 Analyse Numerique Corrige

    Transféré par

    idk04319
    129 vues2 pages
    Ce document présente des exercices sur le calcul numérique. Le premier exercice détaille le calcul de l'intégrale d'une fonction sur un intervalle en utilisant la méthode des trapèzes. Le deuxième exercice examine la convergence de la méthode de Newton-Raphson pour trouver les zéros d'une fonction.

    Description originale:

    Mo

    Titre original

    devoir-14-analyse-numerique-corrige

    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 des exercices sur le calcul numérique. Le premier exercice détaille le calcul de l'intégrale d'une fonction sur un intervalle en utilisant la méthode des trapèzes. Le deuxième exercice examine la convergence de la méthode de Newton-Raphson pour trouver les zéros d'une fonction.

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

    Devoir 14 Analyse Numerique Corrige

    Transféré par

    idk04319
    Ce document présente des exercices sur le calcul numérique. Le premier exercice détaille le calcul de l'intégrale d'une fonction sur un intervalle en utilisant la méthode des trapèzes. Le deuxième exercice examine la convergence de la méthode de Newton-Raphson pour trouver les zéros d'une fonction.

    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 2

    informatique commune

    Contrôle d’informatique

    Exercice 1
    a) Considérons tout d’abord la subdivision en 2n intervalles. Ces composantes x0 = a, x1 , . . . , x2n = b sont définies par :
    b−a
    xk = a + k et nous avons :
    2n
    2n−1
    b − a X f (xk ) + f (xk+1 )
    T2n (f ) = .
    2n 2
    k=0

    La subdivision du segment [a, b] en n intervalles correspond aux points x0 , x2 , x4 , . . . , x2n donc :


    n−1
    b − a X f (x2k ) + f (x2k+2 )
    Tn (f ) = .
    n 2
    k=0

    Notons enfin que pour tout k ∈ ~0, n − 1, x2k+1 est le milieu du segment [x2k , x2k+2 ] donc :
    n−1
    b−a X
    Mn (f ) = f (x2k+1 ).
    n
    k=0

    On en déduit :
     n−1 n n−1  n−1  f (b) − f (a) 
    b−a 1 X 1X b − a X
    X 
    Tn (f ) + Mn (f ) = f (x2k ) + f (x2k ) + f (x2k+1 ) = f (x2k ) + f (x2k+1 ) +
    n 2 2 n 2
    k=0 k=1 k=0 k=0
    2n−1
    b−a f (b) − f (a)
    X 
    = f (xk ) + = 2T2n (f ).
    n 2
    k=0

    T2p (f ) + M2p (f )
    De ceci il résulte que T2p+1 (f ) = .
    2
    b) On rédige une fonction qui calcule Mn (f ) :

    def milieu(f, a, b, n):


    s = 0
    h = (b − a) / n
    x = a + h / 2
    for k in range(n):
    s += f(x)
    x += h
    return h * s

    c) On calcule les termes de la suite (T2p ) à l’aide des formules :

    f (a) + f (b) T2k (f ) + M2k (f )


    T1 = (b − a) et T2k+1 = .
    2 2

    def trap_dicho(f, a, b, epsilon):


    n = 1
    t1 = (b − a) * (f(a) + f(b)) / 2
    t2 = (t1 + milieu(f, a, b, 1)) / 2
    while abs(t2 − t1) > epsilon:
    n *= 2
    t1, t2 = t2, (t2 + milieu(f, a, b, n)) / 2
    return t2

    page 1
    Exercice 2

    Question 1.
    (x − α)2 00 (x − α)2 00 f 00 (α)
    a) On a f (x) = f α) + (x − α)f 0 (α) + f (α) + o((x − α)2 ) = f (α) + o((x − α)2 ) donc h(x) = + o(1). On
    2 00
    2 2
    f (α)
    prolonge h par continuité en posant h(α) = .
    2
    f 0 (x) f (x)
    b) On calcule pour x , α : h0 (x) = −2 .
    (x − α)2 (x − α)3
    Sachant que f est de classe C 3 on peut écrire au voisinage de α :

    (x − α)2 00 (x − α)3 (3) (x − α)2 (3)


    f (x) = f (α) + f (α) + o((x − α)3 ) et f 0 (x) = (x − α)f 00 (α) + f (α) + o((x − α)2 )
    2 6 2

    f (3) (α)
    et en combinant ces deux développements on obtient h0 (x) = + o(1) : h0 possède une limite en α.
    6
    Ainsi, h est continue sur R, de classe C 1 sur R \ {α} et h0 possède une limite en α donc h est de classe C 1 sur R.

    Question 2.
    f (xn ) (xn − α)h(xn )
    a) La méthode de Newton-Raphson s’écrit : xn+1 = xn − 0
    = xn − donc :
    f (xn ) 2h(xn ) + (xn − α)h0 (xn )

    (xn − α)h(xn ) + (xn − α)2 h0 (xn ) en+1 h(xn ) + (xn − α)h0 (xn )
    xn+1 − α = et = .
    2h(xn ) + (xn − α)h0 (xn ) en 2h(xn ) + (xn − α)h0 (xn )

    f 00 (α) e 1
    Sachant que lim xn = α et que h est de classe C 1 , lim h(xn ) = h(α) = , 0 et lim h0 (xn ) = h0 (α) donc lim n+1 = .
    2 en 2
    Alors qu’en général la méthode de Newton-Raphson est d’ordre 2, cette méthode n’est que d’ordre 1 dans ce cas
    particulier (lorsque le zéro de f est d’ordre 2).
    f (x )
    b) Si l’on reprend ces calculs avec la nouvelle relation xn+1 = xn − p 0 n on obtient cette fois :
    f (xn )

    en+1 (2 − p)h(xn ) + (xn − α)h0 (xn )


    =
    en 2h(xn ) + (xn − α)h0 (xn )
    en+1 p
    ce qui montre que si p , 2 on a lim = 1 − , 0. La méthode est toujours d’ordre 1.
    en 2
    Il reste à examiner le cas où p = 2 :

    en+1 (xn − α)h0 (xn ) en+1 h0 (xn )


    = donc = .
    en 2h(xn ) + (xn − α)h0 (xn ) en2 2h(xn ) + (xn − α)h0 (xn )

    en+1 h0 (α) f (3) (α)


    Ainsi, lim = = et la méthode est donc au moins d’ordre 2 (et d’ordre 2 lorsque f (3) (α) , 0).
    en2 2h(α) 6f 00 (α)

    page 2

    Vous aimerez peut-être aussi