3.4 Méthode des coefficients indéterminés

La méthode des coefficients indéterminés est une technique permettant de trouver des solutions particulières, y_p​à des équations différentielles linéaires non homogènes à coefficients constants

 ay'' +by' + cy = f(x)

Pour appliquer cette méthode, il faut d’abord identifier la forme de la fonction de forçage f(x), puis faire une supposition raisonnée de y_p​ avec des coefficients indéterminés. Cette supposition est ensuite remplacée dans l’équation pour trouver la valeur de ces coefficients. Cette méthode est utile lorsque la fonction de forçage, f(x), est une fonction relativement simple, comme un polynôme, une exponentielle, un sinus ou un cosinus, ou une combinaison de ces fonctions.

 

Exemple 3.4.1 : Forme de la supposition d’une solution particulière

Fonctions de forçage polynomiales : Pour y'' +y' - 3y = 9x^2 + 7x + 5, nous ne connaissons pas de solution particulière. Mais, en regardant f(x), on se demande quel type de fonction laisserait un polynôme. Supposons Y_p = AX^2+Bx+C et trouvons la valeur de A, B, C.

 

Fonctions de forçage exponentielles : Pour y'' - 3y' + 2y = 5e^(4x), nous supposons Y_p = Ae^(4x). Si f(x) était 5e^(2x), il faudrait multiplier notre supposition par x: Y_p = Axe^(2x).

 

Ajuster la supposition en fonction des solutions d’équations complémentaires : Si l’équation complémentaire a une solution correspondant partiellement à f(x), il faut ajuster la supposition en conséquence. Par exemple, si f(x)= (3x-2)e^(4x), il faut commencer par Y_p = (Ax+B)e^(4x). Si e^(4x) est une solution à l’équation non homogène, il faut utiliser Y_p=(Ax+B)xe^(4x). Pour une racine répétée, utiliser Y_p = (Ax+B)x^2e^(4x).

 

Remarque : nous transcrivons Y_p avec un « Y » majuscule pour indiquer qu’il s’agit de notre supposition initiale pour la solution particulière. En revanche, y_p avec un « y » minuscule indique la solution particulière après détermination des coefficients.

 

Exemple 3.4.2 : Résoudre une équation avec une fonction de forçage exponentielle

Trouver la solution générale de l’équation suivante.

 (d^2y)/(dx^2) + 10(dy)/(dx) +25y = 3e^(-2 x)
Afficher/Masquer la solution

 

Trouver la solution complémentaire :

Bien qu’il ne soit pas nécessaire de connaître la solution complémentaire pour trouver la solution particulière, il est utile de la connaître. Comprendre la solution complémentaire permet de faire de meilleures suppositions initiales pour la solution particulière et de les ajuster en conséquence avant de procéder à l’algèbre nécessaire pour déterminer les coefficients indéterminés.

L’équation auxiliaire associée à l’équation complémentaire est r^2+10r+25=0, qui a une racine répétée r=-5. Ainsi, {e^(-5x),xe^(-5x)} est un ensemble fondamental de solutions de l’équation complémentaire.

Supposer la forme de la solution particulière :

Étant donné que f(x) est une fonction exponentielle et que les fonctions exponentielles ne changent jamais d’exposant ou ne disparaissent jamais par différenciation, nous supposons que la solution particulière aura une forme similaire à la composante exponentielle dans f(x). De même, l’exposant dans f(x) diffère de l’exposant dans la solution complémentaire, de sorte qu’aucun ajustement n’est nécessaire.

 Y_p=Ae^(-2x)

Introduire la supposition dans l’équation pour trouver A :

Il faut ensuite introduire la supposition et ses dérivées dans l’équation différentielle afin de déterminer le coefficient indéterminé A.

 Y_p=Ae^(-2x) ,  Y'_p=-2Ae^(-2x) ,  Y''_p=4Ae^(-2x)

 4Ae^(-2x) - 20Ae^(-2x) +25Ae^(-2x) = 3e^(-2 x)

(4A-20A+25A)e^(-2x)=3e^(-2x)

9Ae^(-2x)=3e^(-2x)

9A=3

A=1/3

La solution particulière de l’équation différentielle est donc

 y_p=1/3e^(-2x)

Trouver la solution générale :

La solution générale d’une équation non homogène est

 y = y_p + c_1 y_1 + c_2 y_2.

y_1 et y_2 sont les solutions de l’équation complémentaire et y_p est la solution particulière de l’équation non homogène.

 y=1/3e^(-2x)+c_1e^(-5x)+c_2xe^(-5x)

 

Prenons un exemple

 

 

Exemple 3.4.3 : Fonction de forçage similaire à la solution complémentaire avec racine répétée

Trouver la solution générale de l’équation suivante.

 (d^2y)/(dx^2) + 10(dy)/(dx) +25y = e^(-5 x)
Afficher/Masquer la solution

 

Trouver la solution complémentaire :

L’équation complémentaire est similaire à celle de l’exemple 3.4.2. Ainsi, {e^(-5x),xe^(-5x)} est un ensemble fondamental de solutions de l’équation complémentaire et la solution complémentaire est y_c=c_1e^(-5x)+c_2xe^(-5x).

Supposer la forme de la solution particulière :

Notre supposition initiale est Y_p=Ae^(-5x). Cependant, comme e^(-5x) est aussi la solution complémentaire, il faut ajuster notre supposition. Étant donné que r=-5 est une racine répétée, nous multiplions notre suppositions originale par x^2.

 Y_p=Ax^2e^(-5x)

Introduire la supposition dans l’équation pour trouver A :

Il faut ensuite introduire la supposition et ses dérivées dans l’équation différentielle afin de déterminer le coefficient indéterminé A.

 Y_p=Ax^2e^(-5x),

 Y'_p=Ae^(-5x) (2x-5x^2),

 Y''_p=-5Ae^(-5x)(2x-5x^2)+Ae^(-5x)(2-10x)

 Y''_p=Ae^(-5x)(25x^2-20x+2)

 (d^2y)/(dx^2) + 10(dy)/(dx) +25y = e^(-5 x)

 Ae^(-5x)(25x^2-20x+2)+ 10Ae^(-5x) (2x-5x^2)+ 25Ax^2e^(-5x)=e^(-5x)

En factorisant le terme exponentiel et en rassemblant les termes semblables, on obtient

 Ae^(-5x)(25x^2-20x+2+20x-50x^2+25x^2)=e^(-5x)

Ae^(-5x)(2)=e^(-5x)

2A=1

 A=1/2

La solution particulière de l’équation différentielle est donc

 y_p=1/2x^2e^(-5x)

Trouver la solution générale :

La solution générale est

y=y_p+y_c

 y=1/2x^2e^(-5x)+c_1e^(-5x)+c_2xe^(-5x)

 

Prenons un exemple

 

 

Exemple 3.4.4 : Résoudre le PVI avec une équation non homogène

Résoudre le problème de valeur initiale suivant.

 (d^2y)/(dx^2) + 10(dy)/(dx) +25y = e^(-5 x),       y(0)=2,\ y'(0)=-3
Afficher/Masquer la solution

 

Trouver la solution générale :

L’équation est similaire à celle de l’exemple 3.4.3. La solution générale est donc

 y=1/2x^2e^(-5x)+c_1e^(-5x)+c_2xe^(-5x)

Appliquer les conditions initiales :

Appliquer la condition initiale à y :

 y(0)=2

 c_1e^(0)=2

 c_1=2

Appliquer la condition initiale à y' :

 y'=xe^(-5x)-5/2x^2e^(-5x)-5c_1e^(-5x)+c_2(e^(-5x)-5xe^(-5x))

y'(0)=-3

-5c_1e^0+c_2(e^0)=-3

-5c_1+c_2=-3

En introduisant cela dans c_1=2, on obtient c_2=7.

La solution du problème de valeur initiale est donc

 y=1/2x^2e^(-5x)+2e^(-5x)+7xe^(-5x)

Remarque : les conditions initiales doivent satisfaire la solution tout entière de l’équation non homogène, et non pas seulement la partie complémentaire. Il faut donc appliquer les conditions initiales directement à la solution générale de l’équation non homogène donnée pour déterminer les constantes.

La section suivante synthétise les formes appropriées de supposition pour différents types de fonctions de forçage et explique comment modifier ces suppositions si une partie quelconque de la fonction de forçage f(x) correspond à des solutions de l’équation complémentaire.

 

Méthode des coefficients indéterminés (supposition de Y_p)

Trouver une solution particulière à l’équation différentielle

 ay''+by'+cy=f(x)

 f(x)   Y_p Supposition
 n^(e) degré polynomial  A_nx^n+A_(n-1)x^(n-1)+...+A_1x+A_0
 ae^(rx)  Ae^(rx)
 acos(betax)  Acos(betax)+Bsin(betax)
 bsin(betax)  Acos(betax)+Bsin(betax)
 acos(betax)+bsin(betax)  Acos(betax)+Bsin(betax)

Remarques

1. Produits exponentiels et polynomiaux : Sif(x) ne contient que des fonctions exponentielles ou des produits de fonction exponentielle et des polynômes et si e^(rx) est aussi la solution de l’équation complémentaire associée, il faut alors multiplier la partie exponentielle de Y_p par x pour une racine simple x^2 pour une racine répétée.

2. Racines complexes : Si f(x) a trait à la racine complexe d’une équation complémentaire, c’est-à-dire si  alpha+beta est une racine complexe de l’équation auxiliaire associée, il faut alors multiplier la supposition Y_p par x.

3. Produits exponentiels et trigonométriques/polynomiaux : Si f(x) comprend des produits d’une fonction exponentielle et une fonction polynomiale ou trigonométrique, il ne faut retenir que la partie polynomiale ou trigonométrique de la supposition initiale, puis multiplier par la partie exponentielle de f(x).

4. Produits polynomiaux et trigonométriques : Si f(x) comprend des produits de fonctions polynomiales et trigonométriques, il faut d’abord noter la supposition pour la fonction polynomiale uniquement, puis multiplier par le cosinus approprié, puis ajouter une autre fonction polynomiale supposée avec différents coefficients et multiplier par le sinus approprié.

 

Exemple 3.4.5 : Trouver la forme de la solution particulière

Trouver la forme d’une solution particulière à

 y''+y'-6y=f(x)

où f(x) est

a) 5cos(4x)      b) 3x^2sin(pix)       c) 7xe^(2x) cos(8x)       d) 2e^(-3x)       e) (9x^2+3)e^(2x)

Afficher/Masquer la solution

 

L’équation auxiliaire associée à l’équation est r^2-r-6=0, qui a pour racines r_1=-3 et r_2=2.

a)        Y_p=Acos(4x)+Bsin(4x)

 

b) Cette fonction comprend le produit de fonctions polynomiales (second degré) et trigonométriques. Suivant la remarque 4, il faut d’abord supposer le polynôme et le multiplier par le bon cosinus, puis l’ajouter au produit d’un autre polynôme supposé avec différents coefficients et un sinus.

 Y_p=(A_2x^2+A_1x+A_0)cos(pix)+(B_2x^2+B_1x+B_0)sin(pix)

 

b) Cette fonction comprend le produit de fonctions exponentielles, polynomiales (premier degré) et trigonométriques. Suivant les remarques 3 et 4, il faut d’abord supposer le polynôme et le multiplier par le bon cosinus, puis l’ajouter au produit d’un autre polynôme supposé avec différents coefficients et un sinus. Enfin, il faut multiplier la partie exponentielle

 Y_p=e^(2x)((A_1x+A_0)cos(8x)+(B_1x+B_0)sin(8x))

 

d) Étant donné que r=-3 est la racine de l’équation auxiliaire et donc que e^(-3x) est une solution de l’ensemble fondamental, Ae^(-3x) n’est pas une supposition correcte. Suivant la remarque 1, il faut le multiplier par x. Par conséquent,

 Y_p=Axe^(-3x)

 

b) Cette fonction comprend le produit de fonctions exponentielles et polynomiales (second degré). Suivant la remarque 3, il faut d’abord supposer le polynôme et multiplier la partie exponentielle. La supposition polynomiale sera A_2x^2+A_1x+A_0. La partie exponentielle e^(2x) doit être multipliée par x car e^(2x) est dans l’ensemble fondamental de solutions (remarque 1). Par conséquent,

Y_p=xe^(2x)(A_2x^2+A_1x+A_0)

 

Prenons un exemple

 

Section 3.4 Exercices

  1. Trouve la solution particulière de l’EDO
     y'' + 8y'+ 16y = -12e^(-4 x)
    Afficher/Masquer la réponse

     y_p=-6x^2 e^(-4x)

  2. Trouve la solution générale de l’EDO

     (d^2y)/(dt^2) + 2(dy)/(dt) + 17y = 5e^(-4 t)

    Afficher/Masquer la réponse

    y(t)=c_1e^(- t)cos(4t)+c_2 e^(- t)sin(4t)+ 1/5e^(-4t)

  3. Trouve la solution particulière de l’EDO
     y'' + 12y'+ 40y = 4cos(10 x)
    Afficher/Masquer la réponse

    y_p=2/75sin(10x)-1/75cos(10x)

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

     y'' -3y' -10y = 7e^(-2x),\ \ y(0) = 1,\ y'(0) = -17

    Afficher/Masquer la réponse

    y(x)=3 e^(-2x)-2 e^(5x)-xe^(-2x)

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

     y''+y'-6y = 12t-62,   \ \ y(0) = 4,\ y'(0) = 1

    Afficher/Masquer la réponse

     y(t)=-3e^(-3t)-3e^(2t)-2t+10

Licence

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

Partagez ce livre