În analiza matematică , metoda variației constantelor sau metoda Lagrange este o procedură generală care permite determinarea integralei generale a unei ecuații diferențiale liniare de orice ordin și oricare ar fi funcția continuă {\ displaystyle f (t)} care constituie termenul cunoscut. Această metodă este aplicabilă acolo unde este posibil să se determine n soluții independente ale ecuației omogene asociate și a primitivelor funcțiilor adecvate care oferă soluția unui sistem.
Metoda este ilustrată aici inițial pentru ecuații de primul și al doilea ordin, și apoi generalizată la ecuații de ordin n arbitrar. Variabila de care depinde funcția necunoscută {\ displaystyle y} se numește în toate exemplele {\ displaystyle t} .
Ecuații de ordinul întâi
O ecuație diferențială de ordinul întâi în formă normală arată astfel:
- {\ displaystyle y '(t) + a (t) y (t) = f (t)}
Metoda de variație a constantelor constă în căutarea soluțiilor de tipul:
- {\ displaystyle \ displaystyle {\ tilde {y}} (t) = c (t) e ^ {- A (t)}}
obținut pornind de la o soluție a ecuației omogene asociate:
- {\ displaystyle y '(t) + a (t) y (t) = 0}
de tipul:
- {\ displaystyle y (t) = ce ^ {- A (t)}}
unde este {\ displaystyle A (t)} este o primitivă a {\ displaystyle a (t)} Și {\ displaystyle c} este o constantă arbitrară. Motivul pentru care metoda este numită aceasta se datorează constantei {\ displaystyle c} se transformă în funcție {\ displaystyle c (t)} a fi determinat.
Metoda constă în esență în înlocuirea {\ displaystyle {\ tilde {y}}} în ecuația diferențială originală. Pentru a efectua substituția, este necesar să se calculeze mai întâi prima derivată , folosind regula lui Leibniz :
- {\ displaystyle {\ tilde {y}} '(t) = c' (t) y (t) + c (t) y '(t)}
Înlocuind ceea ce tocmai am obținut în ecuația de pornire, obținem:
- {\ displaystyle c '(t) y (t) + c (t) y' (t) + a (t) {\ tilde {y}} = f (t)}
prin urmare, înlocuind:
- {\ displaystyle c '(t) e ^ {- A (t)} - c (t) a (t) e ^ {- A (t)} + c (t) a (t) e ^ {- A (t)} = f (t)}
Simplificând, obținem:
- {\ displaystyle c '(t) e ^ {- A (t)} = f (t)}
Prin izolarea a ceea ce ne interesează și integrarea ambilor membri:
- {\ displaystyle c (t) = \ int f (t) e ^ {A (t)} dt}
din care integralul general al ecuației complete este:
- {\ displaystyle {\ tilde {y}} = e ^ {- A (t)} \ int f (t) e ^ {A (t)} dt}
În acest moment, singura dificultate este calcularea unei integrale care poate să nu fie imediată, sau chiar să nu se rezolve cu metode analitice.
Ecuații de ordinul doi
O ecuație diferențială de ordinul doi arată astfel:
- {\ displaystyle y '' (t) + a (t) y '(t) + b (t) y (t) = f (t)}
Metoda de variație a constantelor în acest caz constă în căutarea soluțiilor de tipul:
- {\ displaystyle \ displaystyle {\ tilde {y}} = c_ {1} (t) y_ {1} (t) + c_ {2} (t) y_ {2} (t)}
construit pornind de la două soluții {\ displaystyle y_ {1} (t)} Și {\ displaystyle y_ {2} (t)} a ecuației omogene asociate:
- {\ displaystyle y '' (t) + a (t) y '(t) + b (t) y (t) = 0}
Deoarece adesea ecuația omogenă asociată este mai ușor de rezolvat, această metodă este utilă în multe cazuri concrete.
Metoda constă în esență în înlocuirea {\ displaystyle {\ tilde {y}}} în ecuația diferențială originală. Pentru a efectua substituția, este necesar să se calculeze mai întâi prima derivată , folosind regula lui Leibniz :
- {\ displaystyle {\ tilde {y}} '(t) = c_ {1}' (t) y_ {1} (t) + c_ {2} '(t) y_ {2} (t) + c_ {1 } (t) y_ {1} '(t) + c_ {2} (t) y_ {2}' (t)}
Pentru a simplifica calculele, se impune următoarea condiție:
- {\ displaystyle c_ {1} '(t) y_ {1} (t) + c_ {2}' (t) y_ {2} (t) = 0}
Acest lucru are ca rezultat:
- {\ displaystyle {\ tilde {y}} '= c_ {1} (t) y_ {1}' (t) + c_ {2} (t) y_ {2} '(t)}
si in consecinta:
- {\ displaystyle {\ tilde {y}} '' = c_ {1} '(t) y_ {1}' (t) + c_ {2} '(t) y_ {2}' (t) + c_ {1 } (t) y_ {1} '' (t) + c_ {2} (t) y_ {2} '' (t)}
Înlocuind ceea ce tocmai am obținut în ecuația de pornire, obținem:
- {\ displaystyle (c_ {1} '(t) y_ {1}' (t) + c_ {2} '(t) y_ {2}' (t) + c_ {1} (t) y_ {1} ' '(t) + c_ {2} (t) y_ {2}' '(t)) + a (c_ {1} (t) y_ {1}' (t) + c_ {2} (t) y_ { 2} '(t)) + b (c_ {1} (t) y_ {1} (t) + c_ {2} (t) y_ {2} (t)) = f (t)}
prin urmare:
- {\ displaystyle c_ {1} (t) (y_ {1} '' (t) + ay_ {1} '(t) + by_ {1} (t)) + c_ {2} (t) (y_ {2 } '' (t) + ay_ {2} '(t) + by_ {2} (t)) + (c_ {1}' (t) y_ {1} '(t) + c_ {2}' (t ) y_ {2} '(t)) = f (t)}
Primele două completări sunt identice nule, deoarece {\ displaystyle y_ {1} (t)} Și {\ displaystyle y_ {2} (t)} sunt soluții ale ecuației omogene, deci totul se reduce la:
- {\ displaystyle c_ {1} '(t) y_ {1}' (t) + c_ {2} '(t) y_ {2}' (t) = f (t)}
Toate acestea duc la studiul sistemului liniar al două ecuații în necunoscute {\ displaystyle c_ {1} '} Și {\ displaystyle c_ {2} '} :
- {\ displaystyle \ left \ {{\ begin {matrix} c_ {1} '(t) y_ {1} (t) + c_ {2}' (t) y_ {2} (t) = 0 \\ c_ { 1} '(t) y_ {1}' (t) + c_ {2} '(t) y_ {2}' (t) = f (t) \ end {matrix}} \ right.}
Determinantul matricei :
- {\ displaystyle \ left ({\ begin {matrix} y_ {1} (t) & y_ {2} (t) \\ y_ {1} '(t) & y_ {2}' (t) \ end {matrix }} \ dreapta)}
este Wronskianul din {\ displaystyle y_ {1} (t)} Și {\ displaystyle y_ {2} (t)} : acest lucru este nul dacă și numai dacă cele două soluții sunt dependente. Rezultă că în acest caz nu este niciodată nul, iar sistemul are întotdeauna o soluție, dată de:
- {\ displaystyle c_ {1} '(t) = {\ frac {-y_ {2} (t) f (t)} {y_ {2}' (t) y_ {1} (t) -y_ {1} '(t) y_ {2} (t)}} \ qquad \ qquad c_ {2}' (t) = {\ frac {y_ {1} (t) f (t)} {y_ {2} '(t ) y_ {1} (t) -y_ {1} '(t) y_ {2} (t)}}}
Integrarea {\ displaystyle c_ {1} '(t)} Și {\ displaystyle c_ {2} '(t)} este posibil să se obțină fie o soluție specială a ecuației de pornire (integrarea definitivă), fie integrala generală a ecuației de pornire (integrarea la nesfârșit).
Ecuații de ordine n
În cazul ecuațiilor de ordine n :
- {\ displaystyle y ^ {(n)} (t) + a_ {n-1} (t) y ^ {(n-1)} (t) + \ cdots + a_ {0} (t) y (t) = f (t)}
sunt luate în considerare n soluțiile independente ale ecuației omogene și se caută o soluție specială a ecuației sub forma:
- {\ displaystyle {\ tilde {y}} = c_ {1} (t) y_ {1} (t) + c_ {2} (t) y_ {2} (t) + \ cdots + c_ {n} (t ) y_ {n} (t)}
Sistemul liniar este apoi rezolvat în n necunoscute {\ displaystyle c_ {i} '(t)} :
- {\ displaystyle {\ begin {cases} c_ {1} '(t) y_ {1} (t) + c_ {2}' (t) y_ {2} (t) + \ cdots + c_ {n} '( t) y_ {n} (t) = 0 \\ c_ {1} '(t) y_ {1}' (t) + c_ {2} '(t) y_ {2}' (t) + \ cdots + c_ {n} '(t) y_ {n}' (t) = 0 \\\ cdots \ cdots \ cdots \ cdots \\ c_ {1} '(t) y_ {1} ^ {(n-2)} (t) + c_ {2} '(t) y_ {2} ^ {(n-2)} (t) + \ cdots + c_ {n}' (t) y_ {n} ^ {(n-2) } (t) = 0 \\ c_ {1} '(t) y_ {1} ^ {(n-1)} (t) + c_ {2}' (t) y_ {2} ^ {(n-1 )} (t) + \ cdots + c_ {n} '(t) y_ {n} ^ {(n-1)} (t) = f (t) \ end {cases}}}
Determinantul acestui sistem se numește determinant Wronskian și, ca mai sus, se poate arăta că este întotdeauna diferit de zero pornind de la independența soluțiilor ecuației omogene. Funcțiile necunoscute sunt determinate prin integrarea celor n termeni de soluție ai sistemului de mai sus, pentru a obține integralul general al ecuației.
În mod specific, având în vedere o ecuație liniară neomogenă obișnuită:
- {\ displaystyle y ^ {(n)} (t) + \ sum _ {i = 0} ^ {n-1} a_ {i} (t) y ^ {(i)} (t) = b (t) }
este {\ displaystyle y_ {1} (t), \ ldots, y_ {n} (t)} un sistem fundamental de soluții ale ecuației omogene corespunzătoare:
- {\ displaystyle y ^ {(n)} (t) + \ sum _ {i = 0} ^ {n-1} a_ {i} (t) y ^ {(i)} (t) = 0}
Apoi, o soluție specială a ecuației neomogene este dată de:
- {\ displaystyle y_ {p} (t) = \ sum _ {i = 1} ^ {n} c_ {i} (t) y_ {i} (t)}
unde este {\ displaystyle c_ {i} (t)} sunt funcții diferențiate care se presupune că îndeplinesc condițiile:
- {\ displaystyle \ sum _ {i = 1} ^ {n} c_ {i} '(t) y_ {i} ^ {(j)} (t) = 0 \ qquad j = 0, \ ldots, n-2 }
Luând în considerare soluția particulară a ecuației omogene, diferențierea în mod repetat și utilizarea condițiilor anterioare:
- {\ displaystyle y_ {p} ^ {(j)} (t) = \ sum _ {i = 1} ^ {n} c_ {i} (t) y_ {i} ^ {(j)} (t) \ qquad j = 0, \ ldots, n-1}
Cu o diferențiere finală avem:
- {\ displaystyle y_ {p} ^ {(n)} (t) = \ sum _ {i = 1} ^ {n} c_ {i} '(t) y_ {i} ^ {(n-1)} ( t) + \ sum _ {i = 1} ^ {n} c_ {i} (t) y_ {i} ^ {(n)} (t)}
Înlocuind deci soluția particulară din ecuația de pornire și aplicând ultimele două relații obținem:
- {\ displaystyle \ sum _ {i = 1} ^ {n} c_ {i} '(t) y_ {i} ^ {(n-1)} (t) = b (t)}
Această ecuație și cea anterioară sunt sisteme liniare care pot fi rezolvate cu regula lui Cramer :
- {\ displaystyle c_ {i} '(t) = {\ frac {W_ {i} (t)} {W (t)}} \ qquad i = 1, \ ldots, n}
unde este {\ displaystyle W (t)} este Wronskianul sistemului fundamental de soluții și {\ displaystyle W_ {i} (t)} este Wronskianul sistemului fundamental cu coloana a I-a înlocuită cu {\ displaystyle (0,0, \ ldots, b (t))} .
Soluția particulară a ecuației neomogene poate fi scrisă ca:
- {\ displaystyle y_ {p} (t) = \ sum _ {i = 1} ^ {n} y_ {i} (t) \, \ int {\ frac {W_ {i} (t)} {W (t )}} dt}
Bibliografie
- (EN) Earl A. Coddington și Norman Levinson, Theory of Ordinary Differential Equations , New York, McGraw-Hill, 1955.
- ( EN ) WE Boyce și RC DiPrima, ecuații diferențiale elementare și probleme de valoare limită Ediția a 8-a , Wiley Interscience, 1965. , paginile 186-192, 237-241
- ( EN ) Gerald Teschl , Ecuații diferențiale ordinare și sisteme dinamice , Providence , Societatea Americană de Matematică .
Elemente conexe
linkuri externe