📚 الرئيسية

مسابقة دكتوراه 2025Source inconnue — الموضوع 02

مسابقة تخصص · Analyse Numérique & Optimisation

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\displaystyle\int_{-1}^{1} g(t)\, dt sous la forme :

Q(g)=αg(1)+βg(0)+αg(1).Q(g) = \alpha g(-1) + \beta g(0) + \alpha g(1).

  1. Calculer α\alpha et β\beta pour que QQ soit exacte sur R2[X]\mathbb{R}_2[X] (polynômes de degré 2\leq 2).
  2. Montrer que QQ est en fait exacte sur R3[X]\mathbb{R}_3[X], mais pas sur R4[X]\mathbb{R}_4[X].
  3. Par changement de variable affine, en déduire la formule de quadrature sur un intervalle quelconque [xi,xi+1][x_i, x_{i+1}] de milieu mim_i et de longueur h=xi+1xih = x_{i+1}-x_i.
  4. Écrire la formule composite FF obtenue en subdivisant [a,b][a,b] en NN sous-intervalles de même longueur h=(ba)/Nh = (b-a)/N, et montrer qu'elle est stable (coefficients positifs).
  5. En utilisant un développement de Taylor à l'ordre 33 au point mim_i, montrer que l'erreur locale vérifie RiM4h52880|R_i| \leq \dfrac{M_4 h^5}{2880}M4=supg(4)M_4 = \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 α,β\alpha, \beta

On impose l'exactitude sur 1,t,t21, t, t^2 :

  • Pour g=1g=1 : 111dt=2=α+β+α=2α+β\displaystyle\int_{-1}^1 1\,dt = 2 = \alpha + \beta + \alpha = 2\alpha+\beta.
  • Pour g=tg=t : 11tdt=0=α+0+α\displaystyle\int_{-1}^1 t\,dt = 0 = -\alpha + 0 + \alpha, automatiquement vérifiée (symétrie).
  • Pour g=t2g=t^2 : 11t2dt=23=α+0+α=2α    α=13\displaystyle\int_{-1}^1 t^2\,dt = \frac{2}{3} = \alpha + 0 + \alpha = 2\alpha \implies \alpha = \frac13.

De 2α+β=22\alpha+\beta = 2 : β=223=43\beta = 2 - \frac23 = \frac43.

α=13,β=43.\boxed{\alpha = \frac{1}{3}, \qquad \beta = \frac{4}{3}.}

On reconnaît la formule de Simpson : Q(g)=13(g(1)+4g(0)+g(1))Q(g) = \dfrac{1}{3}\big(g(-1) + 4g(0) + g(1)\big).

2. Exactitude sur R3[X]\mathbb{R}_3[X] mais pas R4[X]\mathbb{R}_4[X]

Par linéarité, il suffit de tester g=t3g = t^3 et g=t4g=t^4 (les degrés 2\le 2 sont déjà acquis).

Pour g=t3g=t^3 : 11t3dt=0\displaystyle\int_{-1}^1 t^3\,dt = 0 (fonction impaire), et Q(t3)=13(1)+430+13(1)=0Q(t^3) = \frac13(-1) + \frac43\cdot0 + \frac13(1) = 0. Égalité : QQ est exacte sur les impairs, donc sur R3[X]\mathbb{R}_3[X].

Pour g=t4g=t^4 : 11t4dt=25\displaystyle\int_{-1}^1 t^4\,dt = \frac{2}{5}, et Q(t4)=13(1)+43(0)+13(1)=23Q(t^4) = \frac13(1) + \frac43(0) + \frac13(1) = \frac23.

2325    Q n’est pas exacte sur R4[X].\frac23 \neq \frac25 \implies \boxed{Q \text{ n'est pas exacte sur } \mathbb{R}_4[X].}

Donc QQ est exacte exactement jusqu'au degré 33.

3. Formule sur [xi,xi+1][x_i, x_{i+1}]

Soit φ(t)=mi+h2t\varphi(t) = m_i + \dfrac{h}{2}t le changement de variable affine envoyant [1,1][-1,1] sur [xi,xi+1][x_i,x_{i+1}] (dx=h2dtdx = \frac{h}{2}dt). Alors :

xixi+1g(x)dx=h211g(φ(t))dth213(g(xi)+4g(mi)+g(xi+1)),\int_{x_i}^{x_{i+1}} g(x)\,dx = \frac{h}{2}\int_{-1}^1 g(\varphi(t))\,dt \approx \frac{h}{2}\cdot\frac13\Big(g(x_i) + 4g(m_i) + g(x_{i+1})\Big),

xixi+1g(x)dxh6[g(xi)+4g(mi)+g(xi+1)].\boxed{\int_{x_i}^{x_{i+1}} g(x)\,dx \approx \frac{h}{6}\Big[g(x_i) + 4g(m_i) + g(x_{i+1})\Big].}

4. Formule composite et stabilité

En sommant sur les NN sous-intervalles [xi,xi+1][x_i,x_{i+1}], i=0,,N1i=0,\dots,N-1, avec xi=a+ihx_i = a+ih :

F(g)=h6i=0N1[g(xi)+4g(mi)+g(xi+1)]=h6[g(a)+g(b)+2i=1N1g(xi)+4i=0N1g(mi)].F(g) = \frac{h}{6}\sum_{i=0}^{N-1}\Big[g(x_i) + 4g(m_i) + g(x_{i+1})\Big] = \frac{h}{6}\left[g(a) + g(b) + 2\sum_{i=1}^{N-1}g(x_i) + 4\sum_{i=0}^{N-1}g(m_i)\right].

Tous les coefficients (1,2,41, 2, 4 multipliés par h/6>0h/6 > 0) sont positifs. Une formule de quadrature à coefficients positifs est stable : pour une perturbation gg~ε\|g-\tilde g\|_\infty \leq \varepsilon,

F(g)F(g~)εcoefficients=ε(ba),|F(g)-F(\tilde g)| \leq \varepsilon \sum |\text{coefficients}| = \varepsilon(b-a),

car la somme des coefficients vaut (ba)(b-a) (exactitude sur les constantes). D'où :

F est stable.\boxed{F \text{ est stable.}}

5. Estimation de l'erreur

Sur [xi,xi+1][x_i,x_{i+1}], en effectuant un DL de Taylor à l'ordre 33 de gg au point mim_i et en intégrant (calcul standard de l'erreur de Simpson), on obtient l'erreur locale :

Ri=xixi+1gh6[g(xi)+4g(mi)+g(xi+1)]=h52880g(4)(ξi),ξi[xi,xi+1],R_i = \int_{x_i}^{x_{i+1}} g - \frac{h}{6}[g(x_i)+4g(m_i)+g(x_{i+1})] = -\frac{h^5}{2880}g^{(4)}(\xi_i), \quad \xi_i \in [x_i,x_{i+1}],

d'où :

RiM4h52880,M4=sup[a,b]g(4).\boxed{|R_i| \leq \frac{M_4\, h^5}{2880}, \qquad M_4 = \sup_{[a,b]}|g^{(4)}|.}

Erreur globale : en sommant sur les N=bahN = \dfrac{b-a}{h} sous-intervalles :

E=i=0N1RiNM4h52880=bahM4h52880=(ba)M4h42880.|E| = \left|\sum_{i=0}^{N-1} R_i\right| \leq N\cdot\frac{M_4h^5}{2880} = \frac{b-a}{h}\cdot\frac{M_4h^5}{2880} = \boxed{\frac{(b-a)M_4 h^4}{2880}.}

La méthode composite est donc d'ordre 4 en hh.