4.7 Fonction delta de Dirac (impulsion)

Dans les sections précédentes, nous avons exploré des problèmes de valeur initiale pour des équations différentielles du second ordre à coefficients constants, en nous concentrant sur les cas où la fonction de forçage, f(t), est soit continue, soit continue par morceaux sur l’intervalle [0, oo).

 ay''+by'+cy=f(t)\ ; y(0)=k_0, \ \ y'(0)=k_1

Intéressons-nous maintenant à un autre type de fonction de forçage : celle qui représente une force d’impulsion. Les forces d’impulsion se caractérisent par des amplitudes très importantes sur des intervalles de temps extrêmement courts, à la façon d’une brusque « secousse » ou d’un « pic » soudain dans le système. De telles impulsions se produisent dans divers contextes, notamment dans les circuits électriques lors de l’activation, dans les systèmes mécaniques lors d’une collision ou dans tout scénario où une force soudaine et importante est appliquée pendant une brève période.

A. Fonction delta de Dirac

Pour modéliser mathématiquement ces forces d’impulsion, on utilise la fonction delta de Dirac, à savoir delta(t). La fonction delta de Dirac n’est pas une fonction au sens traditionnel du terme, mais plutôt une fonction généralisée ou une distribution dotée des propriétés suivantes.

1. Nulle partout sauf en 0 :

  delta(t)={(0 \ si\ t!=0),(oo \ si \ t=0):}

2. Intégrale égale à 1 :

int_-oo^oo delta(t) = 1

3. Propriété de criblage :

int_-oo^oo f(t) delta(t) = f(0) pour toute fonction f(t) qui est continue sur l’intervalle contenant t=0

En déplaçant l’argument t dans delta(t), il est possible de modéliser les impulsions qui se produisent à des instants autres que t=0. La fonction delta de Dirac déplacée, delta(t-a), marque un pic à t=a et est définie comme suit :

 delta(t-a)={(0 \si \ t!=a),(oo \si \ t=a):}

La propriété de criblage s’étend donc à

4. Criblage en t=a :

 int_-oo^oo f(t) delta(t-a) = f(a) pour toute fonction f(t) qui est continue sur l’intervalle contenant t=a

B. Transformée de Laplace de la fonction delta de Dirac

La transformée de Laplace offre une moyen pratique traiter la fonction Delta de Dirac dans le cadre de la résolution d’équations différentielles. La transformée d’une fonction delta de Dirac déplacée est donnée par

 \mathcal{L}{delta(t-a)}=e^(-as)(4.7.1)

Il importe de comprendre la fonction Delta de Dirac et ses propriétés pour modéliser et analyser des systèmes soumis à des forces d’impulsion.

 

Exemple 4.7.1 : Résoudre le PVI avec une fonction de forçage impulsive

Trouver la solution du problème de valeur initiale

 y''+16y=4delta(t-pi),   y(0)=1,   y'(0)=0

Afficher/Masquer la solution

 

En prenant la transformée de Laplace de l’équation et en appliquant l’équation de la section 4.7.1 avec a=pi à la fonction delta, on obtient

 s^2 Y(s)-s+16Y(s)=4e^(-pis)

La résolution pour Y(s) donne

 Y(s)=(4e^(-pis) +s)/(s^2+16)

 =e^(-pis) 4/(s^2+16)+s/(s^2+16)

En calculant la transformée inverse de Laplace, on obtient

 y(t)=u(t-pi) sin(4(t-pi)) +cos(4t)

Ce qui est équivalent à

 image

   image

La figure ci-dessous montre y(t). La force d’impulsion est appliquée et ajoute de l’élan au système en t=pi. Pour comparaison, la ligne en pointillés représente le système non perturbé.

 

Prenons un exemple

 

Section 4.7 Exercices

  1. Résous le problème de valeur initiale

     y'' + 25 y = delta(t-4), \quad \ y(0) = y'(0) = 0.

    Afficher/Masquer la réponse

     image

  2. Résous le problème de valeur initiale

     y'' + 4 y = 80 e^(4 t) + delta(t-7), \quad \ y(0) = 11,\ y'(0) = 32

    Afficher/Masquer la réponse

     image

Licence

Équations différentielles© par AMIR TAVANGAR. Tous droits réservés.

Partagez ce livre