الرئيسية

مسابقة دكتوراه 2021Université Ferhat Abbas - Sétif 1

مسابقة تخصص · EDP · المدة: 3سا

إضافة يدوية — Université Ferhat Abbas - Sétif 1 2021

التمرين 1

تمرين 1

Soit l'équation différentielle ordinaire

(E):y2y=2x25.(E) : y' - 2y = 2x^2 - 5.
  1. Montrer que l'équation (E)(E) admet une unique solution polynômiale.
  2. En déduire l'ensemble des solutions de (E)(E) dans R\mathbb{R}.
  3. Déterminer la solution φ\varphi de (E)(E) qui vérifie la condition initiale φ(0)=1\varphi(0) = 1.
الحل

Soit l'équation différentielle : (E):y2y=2x25(E) : y' - 2y = 2x^2 - 5

1) Unique solution polynômiale.

Une fonction φ0\varphi_0 est solution de (E)(E) si :

xR:φ0(x)2φ0(x)=2x25(0,25p)\forall x \in \mathbb{R} : \varphi_0'(x) - 2\varphi_0(x) = 2x^2 - 5 \quad \text{(0,25p)}

Si φ0\varphi_0 est polynômiale, alors :

dφ0<dφ0(0,25p)d^\circ \varphi_0' < d^\circ \varphi_0 \quad \text{(0,25p)}

Donc

d(φ02φ0)=dφ0(0,25p)d^\circ(\varphi_0' - 2\varphi_0) = d^\circ \varphi_0 \quad \text{(0,25p)}

Donc il faut que dφ0=2d^\circ \varphi_0 = 2 (0,25p), et par conséquent φ0\varphi_0 est de la forme :

φ0(x)=ax2+bx+c(0,25p)\varphi_0(x) = ax^2 + bx + c \quad \text{(0,25p)}

Donc

xR:(2ax+b)2(ax2+bx+c)=2x252ax2+(2a2b)x+b2c=2x25(0,25p)\forall x \in \mathbb{R} : (2ax + b) - 2(ax^2 + bx + c) = 2x^2 - 5 \Leftrightarrow -2ax^2 + (2a - 2b)x + b - 2c = 2x^2 - 5 \quad \text{(0,25p)}

Ce qui implique que

{2a=22a2b=0b2c=5{a=1b=1c=2\begin{cases} -2a = 2 \\ 2a - 2b = 0 \\ b - 2c = -5 \end{cases} \Rightarrow \begin{cases} a = -1 \\ b = -1 \\ c = 2 \end{cases}

Donc, la fonction xx2x+2x \mapsto -x^2 - x + 2 est l'unique solution polynômiale de l'équation. (0,25p)

2) Ensemble des solutions de (E)(E) dans R\mathbb{R}.

On résout l'équation homogène associée (EH):y2y=0(E_H) : y' - 2y = 0.

dydx=2ydyy=2dxdyy=2dxln(y)=2x+Cy=e2xeCy=Ke2x, K=eCR(0,75p)\frac{dy}{dx} = 2y \Rightarrow \frac{dy}{y} = 2dx \Rightarrow \int \frac{dy}{y} = \int 2dx \Rightarrow \ln(|y|) = 2x + C \Rightarrow |y| = e^{2x}e^{C} \Rightarrow y = Ke^{2x},\ K = e^C \in \mathbb{R} \quad \text{(0,75p)}

Les solutions générales de (EH)(E_H) sont les fonctions xKe2xx \mapsto Ke^{2x}KRK \in \mathbb{R}.

Donc, la fonction φ\varphi cherchée est de la forme :

φ(x)=x2x+2+Ke2x(0,5p)\varphi(x) = -x^2 - x + 2 + Ke^{2x} \quad \text{(0,5p)}

3) Solution vérifiant φ(0)=1\varphi(0) = 1.

La solution générale de (E)(E) est φ(x)=x2x+2+Ke2x\varphi(x) = -x^2 - x + 2 + Ke^{2x}, donc

φ(0)=2+K=1K=1(0,25p)\varphi(0) = 2 + K = 1 \Leftrightarrow K = -1 \quad \text{(0,25p)}

Par suite

xR:φ(x)=x2x+2e2x(1p)\forall x \in \mathbb{R} : \varphi(x) = -x^2 - x + 2 - e^{2x} \quad \text{(1p)}

التمرين 2

تمرين 2

(Oscillations libres d'une corde vibrante fixée aux extrémités avec 00 déplacement)

En utilisant la méthode de séparation des variables, trouver la solution formelle du problème

a22ux22ut2=0(0<x<, t>0), aRa^2 \frac{\partial^2 u}{\partial x^2} - \frac{\partial^2 u}{\partial t^2} = 0 \quad (0 < x < \ell,\ t > 0),\ a \in \mathbb{R}

avec les conditions aux limites

u(0,t)=0,u(,t)=0,t>0u(0,t) = 0,\quad u(\ell,t) = 0, \quad t > 0

et les conditions initiales

u(x,0)=f(x),u(x,0)t=g(x)(0<x<).u(x,0) = f(x),\quad \frac{\partial u(x,0)}{\partial t} = g(x) \quad (0 < x < \ell).
الحل

On utilise la méthode de séparation des variables : soit u(x,t)=X(x)T(t)u(x,t) = X(x)\,T(t), on obtient

X(x)X(x)=λetT(t)a2T(t)=λ(0,5p)\frac{X''(x)}{X(x)} = -\lambda \quad \text{et} \quad \frac{T''(t)}{a^2 T(t)} = -\lambda \quad \text{(0,5p)}

or

X(x)+λX(x)=0(3)X''(x) + \lambda X(x) = 0 \qquad (3) T(t)+λa2T(t)=0(4)(0,5p)T''(t) + \lambda a^2 T(t) = 0 \qquad (4) \quad \text{(0,5p)}

D'après les conditions aux limites nous avons :

X(0)T(t)=0,X()T(t)=0 i.e. X(0)=0, X()=0(5)(0,5p)X(0)T(t) = 0, \quad X(\ell)T(t) = 0 \ \text{i.e.}\ X(0) = 0,\ X(\ell) = 0 \qquad (5) \quad \text{(0,5p)}

Résolvons (3) :

(i) Si λ=k2<0\lambda = -k^2 < 0, alors X(x)=C1ekx+C2ekxX(x) = C_1 e^{kx} + C_2 e^{-kx} (0,25p)

Utilisons les conditions aux limites :

X(0)=C1+C2,C1ek+C2ek=0 ce qui donne C1=C2=0.(0,25p)X(0) = C_1 + C_2, \quad C_1 e^{k\ell} + C_2 e^{-k\ell} = 0 \ \text{ce qui donne}\ C_1 = C_2 = 0. \quad \text{(0,25p)}

Alors X(x)=0 xX(x) = 0\ \forall x, donc u(x,t)=0u(x,t) = 0. Mais alors le problème donné n'a pas de solution si au moins une des fonctions ff et gg est non nulle. (0,5p)

(ii) Si λ=0\lambda = 0 alors X(x)=C1+C2xX(x) = C_1 + C_2 x (0,25p) — d'après (5) nous obtenons la même conclusion. (0,25p)

(iii) Si λ=k2\lambda = k^2 alors

X(x)=C1coskx+C2sinkx(0,5p)X(x) = C_1 \cos kx + C_2 \sin kx \quad \text{(0,5p)}

(5) donne

X(0)=C1=0,X()=C2sink=0(0,5p)X(0) = C_1 = 0, \quad X(\ell) = C_2 \sin k\ell = 0 \quad \text{(0,5p)}

Si C2=0C_2 = 0 on obtient la solution triviale. Supposons C20C_2 \neq 0 :

sink=0 est vraie pour k=nπ, n=1,2,(0,5p)\sin k\ell = 0 \ \text{est vraie pour}\ k = \frac{n\pi}{\ell},\ n = 1, 2, \dots \quad \text{(0,5p)}

Les valeurs propres du problème considéré sont :

λ=λn=n2π22(nN)(6)(0,25p)\lambda = \lambda_n = \frac{n^2\pi^2}{\ell^2} \quad (n \in \mathbb{N}) \qquad (6) \quad \text{(0,25p)}

Les fonctions propres associées sont

Xn(x)=sinnπx(nN, C2=1)(7)(0,25p)X_n(x) = \sin\frac{n\pi x}{\ell} \quad (n \in \mathbb{N},\ C_2 = 1) \qquad (7) \quad \text{(0,25p)}

Pour λ\lambda donné par (6), la solution de l'équation différentielle ordinaire (4)

T(t)+λa2T(t)=0(0,5p)T''(t) + \lambda a^2 T(t) = 0 \quad \text{(0,5p)}

est sous la forme

T(t)=Ancosnaπt+Bnsinnaπt(nN)(8)(0,5p)T(t) = A_n \cos\frac{na\pi t}{\ell} + B_n \sin\frac{na\pi t}{\ell} \quad (n \in \mathbb{N}) \qquad (8) \quad \text{(0,5p)}

An,BnA_n, B_n sont des constantes arbitraires.

Multiplions (7) et (8), on obtient la solution du problème :

un(x,t)=X(x)T(t)=(Ancosnaπt+Bnsinnaπt)sinnπx,nN(0,25p)u_n(x,t) = X(x)\,T(t) = \left(A_n \cos\frac{na\pi t}{\ell} + B_n \sin\frac{na\pi t}{\ell}\right)\sin\frac{n\pi x}{\ell}, \quad n \in \mathbb{N} \quad \text{(0,25p)}

Utilisons la superposition sur l'ensemble des solutions, on obtient :

u(x,t)=n=1(Ancosnaπt+Bnsinnaπt)sinnπx(0,25p)u(x,t) = \sum_{n=1}^{\infty}\left(A_n \cos\frac{na\pi t}{\ell} + B_n \sin\frac{na\pi t}{\ell}\right)\sin\frac{n\pi x}{\ell} \quad \text{(0,25p)}

Cette série converge uniformément si n=1(An+Bn)\sum_{n=1}^{\infty}(|A_n| + |B_n|), où An,BnA_n, B_n sont les coefficients de Fourier, converge. (0,5p)

التمرين 3

تمرين 3

Considérons le système différentiel non linéaire

{x=2yy=2x3x2\begin{cases} x' = 2y \\ y' = 2x - 3x^2 \end{cases}
  1. Déterminer les points d'équilibre.
  2. Linéariser ce système au voisinage de chaque point d'équilibre.
  3. Déterminer la nature de chaque point d'équilibre.
  4. Considérons la fonction H(x,y)=y2+x3x2H(x,y) = y^2 + x^3 - x^2 : a) Calculer dHdt\dfrac{dH}{dt}. b) Que peut-on déduire de la relation entre HH et le système ?
الحل

Considérons le système différentiel non linéaire

{x=2yy=2x3x2\begin{cases} x' = 2y \\ y' = 2x - 3x^2 \end{cases}

1) Les points d'équilibre.

{2y=02x3x2=0sontM=(x=0, y=0),N=(x=23, y=0)(1p)\begin{cases} 2y = 0 \\ 2x - 3x^2 = 0 \end{cases} \quad \text{sont} \quad M = (x = 0,\ y = 0), \quad N = \left(x = \frac{2}{3},\ y = 0\right) \quad \text{(1p)}

2) Linéarisation.

La matrice jacobienne du système est la suivante :

J(x,y)=(x(2y)y(2y)x(2x3x2)y(2x3x2))=(0226x0)(1p)J(x,y) = \begin{pmatrix} \frac{\partial}{\partial x}(2y) & \frac{\partial}{\partial y}(2y) \\ \frac{\partial}{\partial x}(2x - 3x^2) & \frac{\partial}{\partial y}(2x - 3x^2) \end{pmatrix} = \begin{pmatrix} 0 & 2 \\ 2 - 6x & 0 \end{pmatrix} \quad \text{(1p)}

Donc le linéarisé du système non linéaire au voisinage d'un point d'équilibre (a,b)(a,b) est le système linéaire suivant :

(xy)=(0226a0)(xy){x=2yy=(26a)x(2p)\begin{pmatrix} x' \\ y' \end{pmatrix} = \begin{pmatrix} 0 & 2 \\ 2 - 6a & 0 \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} \quad \Rightarrow \quad \begin{cases} x' = 2y \\ y' = (2 - 6a)x \end{cases} \quad \text{(2p)}

3) Nature de chaque point d'équilibre.

  • Le point d'équilibre (0,0)(0,0) :
J(0,0)=(0220),valeurs propres : 2, 2J(0,0) = \begin{pmatrix} 0 & 2 \\ 2 & 0 \end{pmatrix}, \quad \text{valeurs propres : } -2,\ 2

donc (0,0)(0,0) est un point selle. (1p)

  • Le point d'équilibre (23,0)\left(\frac{2}{3}, 0\right) :
J(23,0)=(0220),valeurs propres : 2i, 2iJ\left(\tfrac{2}{3}, 0\right) = \begin{pmatrix} 0 & 2 \\ -2 & 0 \end{pmatrix}, \quad \text{valeurs propres : } -2i,\ 2i

donc (23,0)\left(\frac{2}{3}, 0\right) est un centre ou foyer faible. (1p)

4) Considérons la fonction H(x,y)=y2+x3x2H(x,y) = y^2 + x^3 - x^2 :

a) On a

dHdt=Hxx+Hyy=(2y)x(y2+x3x2)+(2x3x2)y(y2+x3x2)=0(1p)\frac{dH}{dt} = H_x x' + H_y y' = (2y)\frac{\partial}{\partial x}(y^2 + x^3 - x^2) + (2x - 3x^2)\frac{\partial}{\partial y}(y^2 + x^3 - x^2) = 0 \quad \text{(1p)}

b) On en déduit que la fonction H(x,y)=y2+x3x2=kH(x,y) = y^2 + x^3 - x^2 = k, kRk \in \mathbb{R}, est une intégrale première du système. (1p)