Équation différentielle linéaire/Équations différentielles linéaires d'ordre deux

Une page de Wikiversité.

**Équations différentielles linéaires d'ordre deux**
Leçon : Équation différentielle linéaire

Chapitre 5
Chap. préc. :	Changements de variable

En raison de limitations techniques, la typographie souhaitable du titre, « Équation différentielle linéaire : Équations différentielles linéaires d'ordre deux
Équation différentielle linéaire/Équations différentielles linéaires d'ordre deux », n'a pu être restituée correctement ci-dessus.

[modifier] Introduction

Nous dédions un chapitre à l'étude spécifique des équations d'ordre deux, qui apparaissent souvent en physique et jouissent ainsi d'un formalisme très général. En particulier, ses solutions existent toujours dans le cas d'équations linéaires à coefficients constants.

Pour tracer le lien avec l'étude du traitement du signal, nous présentons également une étude fréquentielle de ces équations, qui amène aux phénomènes de filtres et de résonance. Si des notions seront introduites, ce chapitre n'abordera pas la théorie générale des filtres linéaires.

[modifier] Formalisme général

Soit une équation différentielle linéaire d'ordre deux à coefficient constants. On qualifie chacun de ces cas de « régime » :

$af'' + bf' +cf = d\,$

On suppose bien sûr que a n'est pas nul. En ramenant le coefficient devant la dérivée de plus grand ordre à l'unité :

$f'' + \frac{b}{a}f' + \frac{c}{a}f = \frac{d}{a}$

On introduit la notation suivante :

Notations et définitions

Une équation différentielle linéaire d'ordre deux à coefficients constants peut toujours être écrite de la manière suivante :

$f'' + \frac{\omega_0}{Q}f' + \omega_0^2 f = D$

La quantité $ω 0$ est homogène à une pulsation, elle est appelée pulsation propre de f ;
La quantité $Q$ est adimensionnée, elle est appelée facteur de qualité.

Comme c'est le cas dans de nombreux problèmes physiques, on supposera ici que les équations ne font intervenir que des quantités réelles.

[modifier] Résolution de l'équation

Comme démontré dans les chapitres précédents, la résolution de l'équation peut se faire en étudiant le polynôme caractéristique associé à l'équation homogène :

$P \left(X \right) = X^2 + \frac{\omega_0}{Q}X + \omega_0^2$

Le déterminant de ce polynôme est :

$\Delta_P = \frac{\omega_0^2}{Q^2} - 4 \omega_0^2 = 4 \omega_0^2 \left[\frac{1}{4Q^2} - 1 \right]$

Se présentent alors trois cas :

Soit $Δ P > 0$ , auquel cas P admet deux racines réelles ;
Soit $Δ P = 0$ , alors P admet une racine double ;
Soit $Δ P < 0$ , auquel cas P admet deux racines complexes conjuguées.

Dans ce dernier cas, on peut réécrire le discriminant :

$\Delta_P =\left( 2 i \omega_0 \right)^2 \cdot \left[1 - \frac{1}{4Q^2} \right]$ .

[modifier] Cas où Δ_P > 0

Dans le cas où P possède deux racines réelles, :

$\begin{align} r_1 & = - \frac{\omega_0}{2Q} - \frac{\sqrt{\Delta}}{2} \\ \ & = - \frac{\omega_0}{2Q} - \omega_0 \cdot \sqrt{\frac{1}{4Q^2}-1} \\ \ & = - \frac{\omega_0}{2Q} \left[ 1 + 2Q \sqrt{\frac{1}{4Q^2}-1} \right] \\ \ & = - \frac{\omega_0}{2Q} \left[ 1 + \sqrt{1-4Q^2} \right] \end{align}$

De même, l'autre racine est :

$r_2 = - \frac{\omega_0}{2Q} \left[ 1 - \sqrt{1-4Q^2} \right]$

La solution homogène est ainsi de la forme :

$f(t) = A e^{r_1t} + B e^{r_2t}$

La solution générale est somme de la solution homogène ci-dessus et d'une constante (solution particulière).

[modifier] Cas où Δ_P = 0

Dans ce cas, $Δ$ s'annule, la racine est racine double :

$r = - \frac{\omega_0}{2Q}$

La solution homogène est donc de la forme :

$f(t) = \left( At + B \right) e^{rt}$

La solution générale est somme de la solution homogène ci-dessus et d'une constante (solution particulière).

[modifier] Cas où Δ_P < 0

On utilise ici la seconde forme du discriminant. Elle permet en effet d'écrire plus facilement les racines de P :

$r_{\pm} = -\frac{\omega_0}{2Q} \pm i \omega_0 \sqrt{1 - \frac{1}{4Q^2}}$

que l'on réécrit :

$r_{\pm} = -\frac{\omega_0}{2Q} \pm i \Omega$ avec

Ω

pseudo-pulsation du système.

La solution (réelle) homogène est donc de la forme :

$\left[ A \cos \left( \Omega t \right) + B \sin \left( \Omega t \right) \right] e^{-\frac{\omega_0 t}{2Q}}$

ou encore :

$C e^{-\frac{\omega_0 t}{2Q}} \cos \left( \Omega t + \phi \right)$ (en posant

A = C cosφ

et

B = - C sinφ

)

La solution générale est somme de la solution homogène ci-dessus et d'une constante (solution particulière).

[modifier] Lien avec le facteur de qualité

Régimes et facteur de qualité

Les trois cas étudiés ci-dessus sont liés à la valeur du facteur de qualité :

Si Q = 1/2, le discriminant s'annule : c'est le régime critique ;
Si Q > 1/2, le discriminant est négatif : c'est le régime pseudo-périodique ;
Si Q < 1/2, le discriminant est positif : c'est le régime amorti.

[modifier] Principe de superposition

Ce principe est un théorème important des équations linéaires. Nous le rappelons ici :

Principe de superposition

Soit l'équation différentielle :

f'' + b f' + c f = d 1 + d 2

Alors une solution est donnée par :

f = f 1 + f 2

où ces deux fonctions sont respectivement solutions de :

f'' + b f' + c f = d 1

f'' + b f' + c f = d 2

Ce principe se généralise à des sommes infinies.

[modifier] Analyse fréquentielle

[modifier] Mise en forme

On suppose que le système répond à une excitation sinusoïdale. Physiquement, une telle étude se justifie par le théorème de Fourier, d'après lequel tout signal (suffisamment régulier) peut être décomposé en somme (possiblement infinie) de sinusoïdes. En vertu du principe de superposition, on peut étudier indépendamment chaque sinusoïde.

Pour plus de simplicité, on utilise la notation complexe. L'équation différentielle :

$f'' + \frac{\omega_0}{Q}f' + \omega_0^2f = D$

S'écrit alors :

$-\omega^2f + i \omega \frac{\omega_0}{Q} f + \omega_0^2f = D$

$\left[\omega_0^2-\omega^2 + i \frac{\omega \omega_0}{Q} \right] f = D$

Ce qui donne :

$f \left( i \omega \right) = \frac{D}{\left( \omega_0^2 - \omega^2 \right) + i \frac{\omega \omega_0}{Q}}$

Le module de cette fonction est :

$F \left( i \omega \right) = \frac{D}{ \sqrt{\left(\omega_0^2 - \omega^2 \right)^2 + \frac{\omega^2 \omega_0^2}{Q^2}}}$

[modifier] Résonance

Tout d'abord, on voit que le dénominateur de F peut atteindre un minimum. cela signifie seulement qu'il y a une forte augmentation de f, car des effets non-linéaires apparaissent qui viennent toujours atténuer les phénomènes physiques, et qui ne sont pas pris en compte ici.

Résonance

On appelle résonance un maximum du module de la fonction de transfert d'un système.

Cela se produit lorsque la pulsation égale la pulsation propre. Alors :

$F \left( i \omega_0 \right) = \frac{DQ}{\omega_0^2} \in \mathbb R$

[modifier] Comportement asymptotique

Lorsque $\omega \gg \omega_0$ , on a :

$F \left( i \omega \right) \approx \frac{D}{\omega^2}$

En passant au logarithme, le gain en décibels est :

$G_{\rm dB} \left( \omega \right) \approx 20 \log D - 40 \log \omega$

Lorsque $\omega \to 0$ , on a :

$F \left( i \omega \right) \approx \frac{D}{\omega_0^2}$

$G_{\rm dB} \left( \omega \right) \approx 20 \log D - 40 \log \omega_0$

Dans ce cas, le système décrit par l'équation différentielle est un passe-bande : certaines fréquences l'excitent — en dehors de celles-ci, les excitation sont atténuées. Notons que le cas général dépend de D : ici, c'est une constante.

[modifier] Différents filtres

Les différents filtres d'ordre deux

Les différents comportements possibles sont :

Passe-bande (le cas étudié ici) ;
Coupe-bande, ou réjecteur (le cas opposé) ;
Passe-bas (les seules basses fréquences passent) ;
Passe-haut (les seules hautes fréquences passent) ;

[modifier] Équations à coefficients non-constants

[modifier] Cas homogène

Définition

Une relation de la forme :

$a(t)f''(t) + b(t)f'(t) + c(t)f(t) = 0\,$

où a, b et c sont des fonctions au moins continues, a ne s'annulant pas, est appelée équation différentielle linéaire homogène d'ordre deux à coefficients non-constants.

Disons-le d'entrée de jeu : il n'existe pas de solution générale à un tel problème. Cela a été évoqué dans le chapitre 4 sur la théorie générale. Certaines fonctions (comme les fonctions de Bessel) qui ne peuvent pas être explicitées, sont justement définies comme solutions d'une équation différentielle d'ordre deux.

Néanmoins, restent valables les théorèmes généraux :

l'ensemble des solutions forme un espace vectoriel de dimension deux ;
une condition initiale étant donnée, il existe une unique solution maximale sur un ouvert de $\scriptstyle \mathbb R$ .

Il n'est donc pas toujours possible d'exhiber une solution complète à un tel problème. Parfois, un changement de variable, comme évoqué dans un chapitre précédent, est utile et permet de se ramener à un cas plus simple — sinon... que peut-on faire ?

Le théorème du wronksien, décrit dans la section suivante, est un outil précieux : connaissant une solution à cette équation différentielle, il permet d'en trouver une seconde, linéairement indépendante. On obtient ainsi une base de l'ensemble des solutions, ce qui revient à résoudre complètement le problème. Rappelons toutefois qu'il n'existe pas de méthode pour trouver la première solution.

[modifier] Wronskien

Matrice wronskienne et wronskien

Soient u et v deux fonctions dérivables. La matrice wronskienne associée au couple (u, v) est la matrice :

$\mathcal W_{u,v} = \begin{pmatrix} u & v \\ u' & v' \end{pmatrix}$

On appelle wronskien du couple (u, v) le déterminant de $\scriptstyle \mathcal W_{u,v}$ . Il est noté W_u,v. Autrement dit :

$W_{u,v} = uv' - u'v\,$

Supposons que u et v sont deux solutions d'une équation telle que celle définie à la section précédente. Elles vérifient ainsi :

$au'' + bu' + cu = 0\,$

$av'' + bv' + cv = 0\,$

Dérivons leur wronskien :

$\begin{align} W_{u,v}' = \left( uv' - u'v \right)' & = u'v' + uv'' -u''v - u'v' \\ \ & = uv'' - u''v \\ \ & = u \left( - \frac{b}{a}v' - \frac{c}{a}v \right) - v \left( - \frac{b}{a}u' - \frac{c}{a}u \right) \\ \ & = -uv' \frac{b}{a} - uv \frac{c}{a} +u'v \frac{b}{a} + uv\frac{c}{a} \\ \ & = -\frac{b}{a} \left( uv' - v'u\right)\end{align}$

Soit en conclusion :

$W_{u,v}' + \frac{b}{a} W_{u,v} = 0$

Il s'agit d'une équation différentielle linéaire homogène à coefficients non-constants d'ordre un, que l'on sait complètement résoudre : en posant $Ω$ une primitive de la fonction b/a, la solution s'écrit :

$W_{u,v} = W_0 \exp \left( - \Omega t \right)$

On peut alors utiliser le wronskien dans deux buts :

trouver une seconde solution linéairement indépendante connaissant une première solution de l'équation ;
vérifier que deux solutions de l'équation sont linéairement indépendantes et forment une base de l'espace des solutions.

Rappelons le théorème utile :

Théorème du wronskien

Le wronskien de deux fonctions est soit :

toujours nul, auquel cas elles sont liées ;
toujours non-nul, auquel cas elles sont linéairement indépendantes.

Tout cela mérite bien un exemple simple :

Utilisation du wronskien

On cherche à résoudre l'équation suivante :

x 2 y'' - x y' + y = 0

On travaille sur l'ouvert $]0; + \infty[$ , de sorte que $x 2$ ne s'annule pas.

Une solution immédiate est la fonction $y : x\mapsto x$ . Nous allons nous en servir pour déterminer une autre solution avec le wronskien.

On note z la seconde solution de l'équation, alors le wronskien s'écrit : xz' - z.

Le wronskien vérifie l'équation : $W' - \frac{1}{x}W = 0$ .

On sait résoudre cela avec les outils du premier chapitre, et on obtient finalement que toutes les solutions de l'équation sur cet intervalle sont de la forme :

$y : x \mapsto Ax + W_0 x \ln x$

[modifier] Cas général

La résolution du cas général se résume, en fait, au gros problème de trouver une solution particulière. En effet, connaissant l'ensemble des solutions homogènes et une seule solution particulière, on trouve l'ensemble des solutions générales.

On peut tenter plusieurs approches pour trouver une telle solution :

Utiliser le principe de superposition : on trouve plusieurs solutions non-homogènes éventuellement plus simples, et on les somme ;
Rechercher des solutions sous une forme particulière : polynômes, puissances, exponentielles, fonctions trigonométriques, séries entières... cela est parfois efficace ;
Utiliser la méthode de « variation de la constante » évoquée dans le premier chapitre.

Équation différentielle linéaire

Changements de variable

Équation différentielle linéaire/Équations différentielles linéaires d'ordre deux

Une page de Wikiversité.

Sommaire

[modifier] Introduction

[modifier] Formalisme général

[modifier] Résolution de l'équation

[modifier] Cas où Δ_P > 0

[modifier] Cas où Δ_P = 0

[modifier] Cas où Δ_P < 0

[modifier] Lien avec le facteur de qualité

[modifier] Principe de superposition

[modifier] Analyse fréquentielle

[modifier] Mise en forme

[modifier] Résonance

[modifier] Comportement asymptotique

[modifier] Différents filtres

[modifier] Équations à coefficients non-constants

[modifier] Cas homogène

[modifier] Wronskien

[modifier] Cas général

Affichages

Outils personnels

Rechercher

Navigation

Contribuer

Imprimer / exporter

Boîte à outils

Équation différentielle linéaire/Équations différentielles linéaires d'ordre deux

Une page de Wikiversité.

Sommaire

[modifier] Introduction

[modifier] Formalisme général

[modifier] Résolution de l'équation

[modifier] Cas où ΔP > 0

[modifier] Cas où ΔP = 0

[modifier] Cas où ΔP < 0

[modifier] Lien avec le facteur de qualité

[modifier] Principe de superposition

[modifier] Analyse fréquentielle

[modifier] Mise en forme

[modifier] Résonance

[modifier] Comportement asymptotique

[modifier] Différents filtres

[modifier] Équations à coefficients non-constants

[modifier] Cas homogène

[modifier] Wronskien

[modifier] Cas général

Affichages

Outils personnels

Rechercher

Navigation

Contribuer

Imprimer / exporter

Boîte à outils

[modifier] Cas où Δ_P > 0

[modifier] Cas où Δ_P = 0

[modifier] Cas où Δ_P < 0