Concours de Doctorat 2025 — Épreuve commune, Analyse Numérique (15/02/2025)
التمرين 1
Formule de quadrature de type Simpson sur $[-1,1]$ et formule composite
#analyse numérique#intégration numérique#formule de Simpson#ordre de convergence#erreur de quadrature
On cherche une formule de quadrature approchant ∫−11g(t)dt sous la forme :
Q(g)=αg(−1)+βg(0)+αg(1).
Calculer α et β pour que Q soit exacte sur R2[X] (polynômes de degré ≤2).
Montrer que Q est en fait exacte sur R3[X], mais pas sur R4[X].
Par changement de variable affine, en déduire la formule de quadrature sur un intervalle quelconque [xi,xi+1] de milieu mi et de longueur h=xi+1−xi.
Écrire la formule composite F obtenue en subdivisant [a,b] en N sous-intervalles de même longueur h=(b−a)/N, et montrer qu'elle est stable (coefficients positifs).
En utilisant un développement de Taylor à l'ordre 3 au point mi, montrer que l'erreur locale vérifie ∣Ri∣≤2880M4h5 où M4=sup∣g(4)∣, puis en déduire une majoration de l'erreur globale.
Remarque : L'énoncé original de cette épreuve (première page, contexte général) n'est pas présent sur le scan disponible (page 2 uniquement) ; l'exercice a été reconstitué fidèlement à partir des questions visibles sur la formule de quadrature de Simpson et son erreur composite.
◀الحل
1. Détermination de α,β
On impose l'exactitude sur 1,t,t2 :
Pour g=1 : ∫−111dt=2=α+β+α=2α+β.
Pour g=t : ∫−11tdt=0=−α+0+α, automatiquement vérifiée (symétrie).
Pour g=t2 : ∫−11t2dt=32=α+0+α=2α⟹α=31.
De 2α+β=2 : β=2−32=34.
α=31,β=34.
On reconnaît la formule de Simpson : Q(g)=31(g(−1)+4g(0)+g(1)).
2. Exactitude sur R3[X] mais pas R4[X]
Par linéarité, il suffit de tester g=t3 et g=t4 (les degrés ≤2 sont déjà acquis).
Pour g=t3 :∫−11t3dt=0 (fonction impaire), et Q(t3)=31(−1)+34⋅0+31(1)=0. Égalité : Q est exacte sur les impairs, donc sur R3[X].
Pour g=t4 :∫−11t4dt=52, et Q(t4)=31(1)+34(0)+31(1)=32.
32=52⟹Q n’est pas exacte sur R4[X].
Donc Q est exacte exactement jusqu'au degré 3.
3. Formule sur [xi,xi+1]
Soit φ(t)=mi+2ht le changement de variable affine envoyant [−1,1] sur [xi,xi+1] (dx=2hdt). Alors :
Tous les coefficients (1,2,4 multipliés par h/6>0) sont positifs. Une formule de quadrature à coefficients positifs est stable : pour une perturbation ∥g−g~∥∞≤ε,
∣F(g)−F(g~)∣≤ε∑∣coefficients∣=ε(b−a),
car la somme des coefficients vaut (b−a) (exactitude sur les constantes). D'où :
F est stable.
5. Estimation de l'erreur
Sur [xi,xi+1], en effectuant un DL de Taylor à l'ordre 3 de g au point mi et en intégrant (calcul standard de l'erreur de Simpson), on obtient l'erreur locale :