On débute par la résolution d'une équation différentielle simple. Le grand avantage de cette méthode est qu'elle englobe, dans sa procédure, les conditions initiales. C'est particulièrement appréciables dans les problèmes d'évolutions issus de des différents champs de la Physique.
Question 1
Soit f une fonction univariée et causale qui doit satisfaire aux deux conditions initiales suivantes : f(t=0)=1 et f′(t=0)=dtdf)t=0=0. En outre, cette fonction numérique f vérifie l'équation différentielle ordinaire (EDO) suivante : dt2d2f(t)+4f(t)=cos(3t).
Déterminer la fonction numérique f sui satisfait aux conditions énoncées ci-avant.
Correction
Comme on recherche une fonction causale, introduisons l'échelon unité dans l'équation différentielle proposée : dt2d2f(t)+4f(t)=cos(3t)U(t) Prenons maintenant la transformée de Laplace de chaque membre de cette équation. On a alors : TLp[dt2d2f(t)+4f(t)]=TLp[cos(3t)U(t)] Par linéarité de la transformée de Laplace, on obtient alors : TLp[dt2d2f(t)]+4TLp[f(t)]=TLp[cos(3t)U(t)] On va adopter la notation usuelle TLp[f(t)]=F(p), avec p∈C. D'après les propriétés de la transformée de Laplace, on a alors : (p2F(p)−pf(0)−f′(0))+4F(p)=p2+9p En faisant usage des conditions initiales, on a alors : (p2F(p)−p−0)+4F(p)=p2+9p Soit : p2F(p)−p+4F(p)=p2+9p En isolant F(p), on obtient : (p2+4)F(p)=p+p2+9p Ainsi : F(p)=p2+4p+(p2+9)(p2+4)p En décomposant en élément simples, on peut écrire que : (p2+9)(p2+4)p=p2+9Ap+B+p2+4Cp+D Dans cette relation, les quatre quantités A, B, C et D sont des nombres réels. Si p=0 alors on constate que : (02+9)(02+4)0=02+9A0+B+02+4C0+D⟺0=9B+4D⟺4B=−9D⟺D=−94B Si maintenant on pose p=1 alors on constate que : (12+9)(12+4)1=12+9A1+B+12+4C1−94B⟺501=10A+B+5C−94B⟺101=2A+B+C−94B Soit : 1=5A+5B+10C−940B⟺1=5A+10C+945B−940B⟺1=5A+95B+10C Soit encore : 1−5A−10C=95B⟺A+2C−51=−91B⟺4A+8C−54=−94B De plus, ceci nous permet d'écrire que : B=−9A−18C+59 Ainsi, on a : (p2+9)(p2+4)p=p2+9Ap−9A−18C+59+p2+4Cp+4A+8C−54 En réduisant au même dénominateur, on trouve que : (p2+9)(p2+4)p=(p2+9)(p2+4)(p2+4)(Ap−9A−18C+59)+(p2+9)(Cp+4A+8C−54) D'où, en développant puis en regroupant après, on obtient : (p2+9)(p2+4)0p3+0p2+1p+0=(p2+9)(p2+4)(A+C)p3+(1−5A−10C)p2+(4A+9C)p+0 Mais, on sait que 1−5A−10C=95B, d'où : (p2+9)(p2+4)0p3+0p2+1p+0=(p2+9)(p2+4)(A+C)p3+95Bp2+(4A+9C)p+0 Par identification, on trouve que : ⎩⎨⎧A+C95B4A+9C===001 Donc C=−A et B=0. De fait D=0 également et on a donc 4A−9A=1. Ainsi −5A=1 et on en déduit immédiatement que : ⎩⎨⎧ABCD====−510510 Ainsi : (p2+9)(p2+4)p=−51p2+91+51p2+41 Ce qui implique que : F(p)=p2+4p−51p2+91+51p2+41⟺F(p)=56p2+4p−51p2+91 Nous allons écrire ceci comme : F(p)=56p2+22p−51p2+321 D'où : TLp[f(t)]=56TLp[cos(2t)U(t)]−51TLp[cos(3t)U(t)] Soit : TLp[f(t)]=TLp[56cos(2t)U(t)]+TLp[−51cos(3t)U(t)] Soit encore : TLp[f(t)]=TLp[56cos(2t)U(t)−51cos(3t)U(t)] En factorisant : TLp[f(t)]=TLp[51(6cos(2t)−cos(3t))U(t)]µ Finalement : f(t)=51(6cos(2t)−cos(3t))U(t) On vérifie graphiquement que les deux conditions initiales sont bien simultanément satisfaites :
Signaler une erreur
Aide-nous à améliorer nos contenus en signalant les erreurs ou problèmes que tu penses avoir trouvés.
Connecte-toi ou crée un compte pour signaler une erreur.