📚 الرئيسية

مسابقة دكتوراه 2017Université Amar Telidji - Laghouat — الموضوع 02

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

Faculté des Sciences, Département de Mathématiques, Concours d'accès au Doctorat 2017-2018, Épreuve : Analyse Mathématique, Date : 18/10/2017

التمرين 1

Exercice 1 (Laghouat 2017) — EDO $y''-(x^4+1)y=0$ et estimation

#équations différentielles#estimation énergétique#unicité

On considère l'équation différentielle y(x)(x4+1)y(x)=0y''(x)-(x^4+1)y(x)=0 pour x[0,+[x\in[0,+\infty[.

  1. Soit yC2([0,+[)y\in C^2([0,+\infty[) une solution telle que y(0)=y(0)=0y(0)=y'(0)=0. Montrer que y0y\equiv 0 sur [0,+[[0,+\infty[.

  2. En déduire l'unicité de la solution du problème de Cauchy associé à cette EDO.

Technique de l'énergie + Grönwall pour l'unicité. Cette méthode s'étend aux EDO linéaires à coefficients continus et aux EDP elliptiques/paraboliques.

الحل
  1. Posons E(x)=y(x)2+y(x)2E(x)=y(x)^2+y'(x)^2 (énergie). Alors E(x)=2yy+2yy=2y(y+y)=2y(y+(x4+1)y)=2yy(x4+2)E'(x)=2yy'+2y'y''=2y'(y+y'')=2y'(y+(x^4+1)y)=2y'y(x^4+2).

Par Cauchy-Schwarz : E(x)(x4+2)2yy(x4+2)(y2+(y)2)=(x4+2)E(x)|E'(x)|\le (x^4+2)\cdot 2|y y'|\le (x^4+2)(y^2+(y')^2)=(x^4+2)E(x).

Donc E(x)(x4+2)E(x)|E'(x)|\le (x^4+2)E(x). Par lemme de Grönwall (intégration) avec E(0)=0E(0)=0 : E(x)E(0)exp0x(t4+2)dt=0E(x)\le E(0)\exp\int_0^x(t^4+2)dt = 0. Donc E(x)=0E(x)=0 pour tout x0x\ge 0, ce qui donne y0y\equiv 0.

  1. Si y1,y2y_1,y_2 sont deux solutions du problème de Cauchy avec mêmes conditions initiales, alors y=y1y2y=y_1-y_2 est solution avec y(0)=y(0)=0y(0)=y'(0)=0. Par 1) y0y\equiv 0, donc y1=y2y_1=y_2.

التمرين 2

Exercice 2 (Laghouat 2017) — Polynôme d'interpolation par différences divisées de Newton

#interpolation#Newton#différences divisées

Soit la fonction f(x)=ln(x)f(x)=\ln(x) connue aux points x0=1x_0=1, x1=2x_1=2, x2=4x_2=4.

  1. Construire le polynôme d'interpolation de Newton P2(x)P_2(x) de ff passant par ces trois points.

  2. Calculer P2(3)P_2(3) et comparer avec f(3)=ln31,0986f(3)=\ln 3\approx 1{,}0986.

  3. Donner une majoration théorique de l'erreur f(3)P2(3)|f(3)-P_2(3)|.

Interpolation de Newton : facile à mettre à jour, pratique numériquement. Erreur 1/(n+1)!\propto 1/(n+1)! (Runge). Pour minimiser, choisir nœuds de Chebyshev sur [a,b][a,b] : xk=(a+b)/2+((ba)/2)cos((2k+1)π/(2(n+1)))x_k=(a+b)/2+((b-a)/2)\cos((2k+1)\pi/(2(n+1))).

الحل
  1. Table des différences divisées : f[x0]=ln1=0f[x_0]=\ln 1=0, f[x1]=ln20,6931f[x_1]=\ln 2\approx 0{,}6931, f[x2]=ln4=2ln21,3863f[x_2]=\ln 4=2\ln 2\approx 1{,}3863.

Ordre 1 : f[x0,x1]=(ln20)/(21)=ln2f[x_0,x_1]=(\ln 2-0)/(2-1)=\ln 2. f[x1,x2]=(ln4ln2)/(42)=ln2/2f[x_1,x_2]=(\ln 4-\ln 2)/(4-2)=\ln 2/2.

Ordre 2 : f[x0,x1,x2]=(ln2/2ln2)/(41)=(ln2/2)/3=ln2/6f[x_0,x_1,x_2]=(\ln 2/2 - \ln 2)/(4-1) = (-\ln 2/2)/3 = -\ln 2/6.

Donc P2(x)=0+ln2(x1)ln26(x1)(x2)P_2(x) = 0 + \ln 2\cdot(x-1) - \dfrac{\ln 2}{6}(x-1)(x-2).

  1. P2(3)=ln22ln2621=2ln2ln23=5ln231,1552P_2(3) = \ln 2\cdot 2 - \dfrac{\ln 2}{6}\cdot 2\cdot 1 = 2\ln 2 - \dfrac{\ln 2}{3} = \dfrac{5\ln 2}{3}\approx 1{,}1552.

Erreur : 1,15521,09860,057|1{,}1552-1{,}0986|\approx 0{,}057.

  1. Majoration : f(x)P2(x)M33!(xx0)(xx1)(xx2)|f(x)-P_2(x)|\le \dfrac{M_3}{3!}|(x-x_0)(x-x_1)(x-x_2)|M3=max[1,4]fM_3=\max_{[1,4]}|f'''|. Or f(x)=2/x3f'''(x)=2/x^3, M3=2M_3=2 (en x=1x=1). En x=3x=3 : (31)(32)(34)=2|(3-1)(3-2)(3-4)|=2. Majoration : f(3)P2(3)(2/6)2=2/30,667|f(3)-P_2(3)|\le (2/6)\cdot 2 = 2/3\approx 0{,}667. La majoration est nettement plus large que l'erreur réelle.