الرئيسية

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

مسابقة عامة · الرياضيات

Uploaded PDF 484646404_122188466804122683_3738905597730218888_n.pdf, page 1/7 — header not visible; metadata uncertain

التمرين 2

Erreur de quadrature par interpolation constante

#analyse numérique#intégration numérique#quadrature#estimation d'erreur

Exercice partiellement visible. Sur une subdivision x_i=a+ih, i=0,…,n, montrer notamment que ∫{x_i}^{x{i+1}} |x-x_i^| dx = ((x_i^-x_i)^2)/2 + ((h-(x_i^-x_i))^2)/2. En posant λ=x_i^-x_i∈[0,h] et g(λ)=λ²/2+(h-λ)²/2, étudier g sur [0,h] et en déduire ∫{x_i}^{x{i+1}} |x-x_i^|dx≤h²/2. En déduire une estimation |I(f)-R_n(f)|≤M(b-a)²/(2n). Si x_i^=(x_i+x_{i+1})/2, montrer que ∫{x_i}^{x{i+1}} |x-x_i^*|dx=h²/4 et en déduire |I(f)-R_n(f)|≤M(b-a)²/(4n).

الحل

1. Calcul de l'intégrale

On coupe l'intégrale au point xi[xi,xi+1]x_i^*\in[x_i,x_{i+1}] et on pose λ=xixi[0,h]\lambda=x_i^*-x_i\in[0,h] :

xixi+1xxidx=xixi(xix)dx+xixi+1(xxi)dx=λ22+(hλ)22.\int_{x_i}^{x_{i+1}}|x-x_i^*|\,dx =\int_{x_i}^{x_i^*}(x_i^*-x)\,dx+\int_{x_i^*}^{x_{i+1}}(x-x_i^*)\,dx =\frac{\lambda^2}{2}+\frac{(h-\lambda)^2}{2}.

D'où

xixi+1xxidx=(xixi)22+(h(xixi))22.\boxed{\int_{x_i}^{x_{i+1}}|x-x_i^*|\,dx=\frac{(x_i^*-x_i)^2}{2}+\frac{\big(h-(x_i^*-x_i)\big)^2}{2}.}

2. Étude de gg

g(λ)=λ22+(hλ)22,g(λ)=2λh,g(λ)=2>0.g(\lambda)=\frac{\lambda^2}{2}+\frac{(h-\lambda)^2}{2}, \qquad g'(\lambda)=2\lambda-h, \qquad g''(\lambda)=2>0.

gg décroît sur [0,h2][0,\tfrac{h}{2}] et croît sur [h2,h][\tfrac{h}{2},h] :

min[0,h]g=g ⁣(h2)=h24,max[0,h]g=g(0)=g(h)=h22.\min_{[0,h]}g=g\!\left(\frac{h}{2}\right)=\frac{h^2}{4}, \qquad \max_{[0,h]}g=g(0)=g(h)=\frac{h^2}{2}.

Donc

xixi+1xxidx  h22.\boxed{\int_{x_i}^{x_{i+1}}|x-x_i^*|\,dx\ \le\ \frac{h^2}{2}.}

3. Estimation générale

Si fM|f'|\le M, l'erreur locale de l'approximation par la constante f(xi)f(x_i^*) vérifie

xixi+1(f(x)f(xi))dxMxixi+1xxidxMh22.\left|\int_{x_i}^{x_{i+1}}\big(f(x)-f(x_i^*)\big)\,dx\right| \le M\int_{x_i}^{x_{i+1}}|x-x_i^*|\,dx\le\frac{Mh^2}{2}.

En sommant sur les nn sous-intervalles, avec h=banh=\dfrac{b-a}{n} :

I(f)Rn(f)  nMh22=M(ba)22n.\boxed{|I(f)-R_n(f)|\ \le\ n\cdot\frac{Mh^2}{2}=\frac{M(b-a)^2}{2n}.}

4. Cas du point milieu

Si xi=xi=xi+xi+12x_i^*=\overline{x_i}=\dfrac{x_i+x_{i+1}}{2}, alors λ=h2\lambda=\dfrac{h}{2} et

xixi+1xxidx=g ⁣(h2)=h24.\boxed{\int_{x_i}^{x_{i+1}}|x-x_i^*|\,dx=g\!\left(\frac{h}{2}\right)=\frac{h^2}{4}.}

D'où l'estimation améliorée

I(f)Rn(f)  M(ba)24n.\boxed{|I(f)-R_n(f)|\ \le\ \frac{M(b-a)^2}{4n}.}

التمرين 3

Série entière ∑ n z^{2n-1} et EDO associée

#analyse#séries entières#rayon de convergence#équations différentielles#récurrence
  1. Trouver le rayon de convergence R de la série entière ∑{n=1}^{∞} n z^{2n-1}. 2) Trouver la fonction somme f(x) pour x∈]-R,R[. 3) Démontrer que si y(x)=∑{n=0}^{∞} a_n x^n est une série entière solution de (x²-1)y''+6xy'+4y=0, alors a_{n+2}=((n+4)/(n+2))a_n, n∈N. 4) Trouver la solution y=φ(x) sur I telle que φ(0)=1 et φ'(0)=0.

Remarque : la récurrence a2k=k+1a_{2k}=k+1 reconstruit exactement la dérivée de la série géométrique en w=x2w=x^2 ; on peut vérifier que ϕ(x)=(1x2)2\phi(x)=(1-x^2)^{-2} satisfait bien l'équation (2)(2).

الحل

1. Rayon de convergence

On écrit

n=1+nz2n1=zn=1+n(z2)n1.\sum_{n=1}^{+\infty}n\,z^{2n-1}=z\sum_{n=1}^{+\infty}n\,(z^2)^{n-1}.

La série nwn1\sum n\,w^{n-1} a pour rayon 11 en ww, donc la condition est z2<1|z^2|<1 :

R=1.\boxed{R=1.}

2. Fonction somme

Pour w<1|w|<1 :

n=1+nwn1=1(1w)2.\sum_{n=1}^{+\infty}n\,w^{n-1}=\frac{1}{(1-w)^2}.

Avec w=x2w=x^2 :

f(x)=x(1x2)2,x]1,1[.\boxed{f(x)=\frac{x}{(1-x^2)^2},\qquad x\in\,]-1,1[.}

3. Relation de récurrence

Soit y(x)=n=0+anxny(x)=\sum_{n=0}^{+\infty}a_nx^n solution de

(x21)y+6xy+4y=0.(2)(x^2-1)\,y''+6x\,y'+4y=0. \tag{2}

On a

y=n=0+(n+2)(n+1)an+2xn,x2y=n=0+n(n1)anxn,6xy=n=0+6nanxn.y''=\sum_{n=0}^{+\infty}(n+2)(n+1)a_{n+2}x^n, \qquad x^2y''=\sum_{n=0}^{+\infty}n(n-1)a_nx^n, \qquad 6xy'=\sum_{n=0}^{+\infty}6n\,a_nx^n.

Le coefficient de xnx^n dans (2)(2) donne

(n+2)(n+1)an+2+[n(n1)+6n+4]an=0.-(n+2)(n+1)\,a_{n+2}+\big[n(n-1)+6n+4\big]a_n=0.

Or

n(n1)+6n+4=n2+5n+4=(n+1)(n+4),n(n-1)+6n+4=n^2+5n+4=(n+1)(n+4),

donc

(n+2)(n+1)an+2=(n+1)(n+4)an,(n+2)(n+1)\,a_{n+2}=(n+1)(n+4)\,a_n,

soit

an+2=n+4n+2an,nN.\boxed{a_{n+2}=\frac{n+4}{n+2}\,a_n,\qquad n\in\mathbb{N}.}

4. Solution avec ϕ(0)=1\phi(0)=1, ϕ(0)=0\phi'(0)=0

Les conditions donnent a0=1a_0=1, a1=0a_1=0 ; tous les coefficients impairs sont nuls. Pour n=2kn=2k :

a2k+2=2k+42k+2a2k=k+2k+1a2k  a2k=k+1.a_{2k+2}=\frac{2k+4}{2k+2}\,a_{2k}=\frac{k+2}{k+1}\,a_{2k} \ \Longrightarrow\ a_{2k}=k+1.

Donc

ϕ(x)=k=0+(k+1)x2k.\phi(x)=\sum_{k=0}^{+\infty}(k+1)\,x^{2k}.

Or, pour x<1|x|<1 :

k=0+(k+1)wk=1(1w)2,w=x2,\sum_{k=0}^{+\infty}(k+1)\,w^{k}=\frac{1}{(1-w)^2},\qquad w=x^2,

d'où

ϕ(x)=1(1x2)2,x]1,1[.\boxed{\phi(x)=\frac{1}{(1-x^2)^2},\qquad x\in\,]-1,1[.}