Université Pierre et Marie Curie Travaux dirigés

Licence de Mathématiques LM83 Année universitaire 2005-2006

Equations diérentielles. Méthodes de résolution numérique.

TD no 3 : Intégration numérique

Exercice 1. (extrait de l'examen partiel, mars 2000)

Soit f ∈ C 1 ([0, 1]. On considère la formule de quadrature élémentaire :

Z 1
f (x) dx ∼ w0 f (0) + w1 f 0 (0) + w2 f 0 (ξ)
0

où ξ ∈]0, 1] et w0 , w1 , w2 sont des réels. On pose

Z 1
E(f ) = f (x) dx − [w0 f (0) + w1 f 0 (0) + w2 f 0 (ξ)]
0

a) Déterminer les paramètres ξ, w0 , w1 , w2 pour que la formule de quadrature soit exacte si f

est un polynôme de degré inférieur ou égal à 3.

b) Les paramètres ξ, w0 , w1 , w2 étant ainsi xés, calculer E(x → x4 ) et en déduire l'ordre de la

méthode. Déterminer le noyau de Peano K associé à la méthode.

c) En déduire une expression de l'erreur E(f ) lorsque f ∈ C 4 ([0, 1]).

d) A l'aide d'un changement de variable, construire une méthode de quadrature élémentaire

sur un intervalle [a, b] et donner la valeur de l'erreur.

Exercice 2. Interpolation/intégration numérique utilisant les polynômes de Hermite

a) Soit αi , i = 1 . . . 4 quatre constantes réelles. Montrer qu'il existe un seul polynôme p ∈ P3

vériant les égalités :

p(0) = α1 , p0 (0) = α2 , p(1) = α3 , p0 (1) = α4 . (1)

b) Calculer les quatre polynômes pi ∈ P3 , i = 1 . . . 4 vériant les relations (1), avec respective-
ment :
(α1 , α2 , α3 , α4 )i = (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1). (2)
Montrer que le polynôme de la question a) s'écrit comme une combinaison linéaire de pi :
4
X
p(x) = αi pi (x). (3)
i=1

1
c) Soit f une fonction de classe C 4 ([0, 1]) et soit pf le polynôme de la question a) avec :
α1 = f (0), α2 = f 0 (0), α3 = f (1), α4 = f 0 (1). (4)

Montrer que pour chaque point x ∈]0, 1[, il existe un point ξx ∈]0, 1[ où on a l'égalité :

π(x) (4)
f (x) − pf (x) = f (ξx ), avec π(x) = x2 (x − 1)2 . (5)
4!

d) On approxime l'intégrale d'une fonction f à l'aide de la formule de quadrature

Z 1 Z 1
f (x)dx ∼ pf (x)dx. (6)
0 0

Montrer
Z 1 que cette formule est exacte pour les polynômes de P3 . Calculer les poids wi =
pi (x)dx, i = 1, . . . , 4 et en déduire une forme explicite de la formule de quadrature.
0

e) Montrer que la fonction x 7→ f (4) (ξx ) est prolongeable par Rcontinuité sur
R 1[0, 1]. En déduire
1
qu'il existe un point η ∈]0, 1[ où l'erreur d'intégration E(f ) = 0 f (x)dx − 0 pf (x)dx s'écrit :

1 (4)
E(f ) = f (η). (7)
720

f) Calculer le noyau de Peano de la méthode élémentaire. Retrouver le résultat démontré au

point e).

g) Soit maintenant [a, b] un intervalle fermé, borné de R et f une fonction de classe C 4 ([a, b]).
En utilisant les résultats précédents, construire une formule de quadrature de la forme :

J(f ) = γ1 f (a) + γ2 f 0 (a) + γ3 f (b) + γ4 f 0 (b), (8)

qui soit exacte pour les polynômes de P3 . Calculer l'erreur

Z b
E(f ) = f (x) dx − J(f ). (9)
a

h) Soit une discrétisation de l'intervalle [a, b] avec des points équidistants

xi = a + ih, h = (b − a)/n, i = 0, 1, . . . , n.

Construire une formule de quadrature composée qui utilise les valeurs de f (xi ) et f 0 (xi ) pour
i = 0, 1, . . . , n. Estimer l'erreur de cette formule de quadrature.

2
Exercice 3. Calcul des coecients de la formule de Newton-Cotes

Soit (NC` ) la formule de quadrature élémentaire de Newton-Cotes de rang ` (avec ` entier ≥ 1) :

Z 1 X̀ (`)
f (x) dx ∼ 2 ωj f (xj ),
−1 j=0

j
où les points xj sont équidistants, d'abscisses xj = −1 + 2 , j = 0, 1, . . . `. On note E` (f )
R1 P `
(`)
l'erreur d'approximation −1 f (x) dx − 2 `j=0 ωj f (xj ).

a) Montrer que :
Z
(`) (−1)(`−j) C`j `
π` (s) Ỳ
ωj = ds, j = 0, 1, . . . , ` avec π` (s) = (s − j).
` · `! 0 s−j
j=0

b) Si le polynôme π` (s) est écrit sous la forme :

`+1
X (`+1) i
π` (s) = Si s,
i=0

(`+1)
montrer que les coecients Si (appelés les nombres de Stirling de 1ère espèce) vérient :
 (`+1)
( 
(1)  S0 =0
S0 = 0 (`+1) (`) (`)
(1) et Si = Si−1 − ` · Si , i = 0, 1, . . . `
S1 = 1 
 (`+1)
S`+1 = 1

c) Si
π` (s) X̀ (`) i
= aj,i s ,
s−j
i=0
(`+1) (`)
en déduire un algorithme de calcul des coecients Sj et aj,i . Trouver les expressions des
(`)
poids ωj par cet algorithme.

d) Construire les formules de Newton-Cotes pour ` = 1, 2, 3.

d) Calculer le noyau de Peano associé aux formules de Newton-Cotes pour ` = 1, 2. Evaluer les
erreurs d'interpolation.

3
Exercice 4. Intégration numérique aux points de Tchebichev.

Nous nous proposons d'établir quelques résultats sur la formule approchée :

Z 1 Z 1
f (x) dx ∼ Pn (x) dx,
−1 −1

où Pn (x) est le polynôme d'interpolation de degré n de f aux points de Tchebichev

µ ¶
2i + 1
xi = cos θi = cos π , i = 0, 1, . . . n.
2n + 2

(a) Pour tout x ∈ [−1, 1] et n ∈ N, on note

Z 1
cos(n arccos x) − cos(n arccos y)
an (x) = dy.
−1 x−y

Montrer que :
Z 1 n
X (−1)i sin θi
Pn (x) dx = ωi f (xi ), avec ωi = an+1 (xi ).
−1 n+1
i=0

(b) Calculer an+1 (x) + an−1 (x) ; en déduire la valeur de

an+1 (x) − 2xan (x) + an−1 (x).

(c) En distinguant deux cas suivant la parité de n, montrer l'égalité :

X 1
sin θ an (cos θ) = 2 sin(nθ) − 4 sin(n − 2q)θ.
4q 2 − 1
1≤ q ≤n/2

(d) En déduire l'expression de ωi :

 
1  X 1
ωi = 2−4 cos(2qθi )
n+1 4q 2−1
1≤ q ≤(n+1)/2

Montrer que ωi > 0 pour i = 0, 1, . . . , n.