O ecuație diferențială liniară de ordinul doi este un tip special de ecuație diferențială liniară .
Definiție
O ecuație diferențială liniară obișnuită de ordinul doi are forma: [1]
- {\ displaystyle y '' + a (x) \ cdot y '+ b (x) \ cdot y = f (x)}
unde este {\ displaystyle a (x)} Și {\ displaystyle b (x)} sunt funcții continue într-un interval real .
Pentru a o rezolva, luați în considerare ecuația diferențială omogenă asociată:
- {\ displaystyle y '' + a (x) \ cdot y '+ b (x) \ cdot y = 0}
care are ca o soluție banală {\ displaystyle y = 0} . Pentru a obține soluții non-banale, dacă {\ displaystyle y_ {1}} Și {\ displaystyle y_ {2}} sunt două soluții liniar independente ale acestei ecuații, de asemenea:
- {\ displaystyle y = c_ {1} y_ {1} + c_ {2} y_ {2}}
este soluția pentru fiecare valoare a constantelor {\ displaystyle c_ {1}} Și {\ displaystyle c_ {2}} . Mai exact, toate soluțiile ecuației omogene asociate sunt de această formă. Deoarece diferența oricăror două soluții ale ecuației neomogene trebuie să fie soluția ecuației omogene, pentru a găsi soluția generală a ecuației neomogene este suficient să găsim o soluție specială {\ displaystyle u (x)} și adăugați soluția generică a ecuației omogene asociate:
- {\ displaystyle y = u (x) + c_ {1} y_ {1} + c_ {2} y_ {2}}
În loc să indicați familia parametrică a tuturor soluțiilor ecuației neomogene, este posibil să vi se solicite să rezolvați ecuația cu valorile inițiale atribuite. Problema Cauchy astfel delimitată are forma:
- {\ displaystyle {\ begin {cases} y (x_ {0}) = y_ {0} \\ y '(x_ {0}) = y_ {1} \ end {cases}}}
iar aceste două condiții servesc la determinarea valorilor constantelor arbitrare asociate cu soluția anterioară pentru ecuația neomogenă, pentru a avea o soluție specială care testează problema la valorile inițiale.
Ecuații cu coeficienți constanți
Ecuația omogenă asociată
Ecuația omogenă asociată are forma:
- {\ displaystyle y '' + a \ cdot y '+ b \ cdot y = 0}
unde este {\ displaystyle a} Și {\ displaystyle b} sunt coeficienți constanți dați. Rezoluția sa constă în căutarea unei soluții precum:
- {\ displaystyle y = e ^ {\ lambda x}}
Înlocuind această expresie în ecuația omogenă anterioară, derivând și evidențiind{\ displaystyle e ^ {\ lambda x}} :
- {\ displaystyle e ^ {\ lambda x} \ left (\ lambda ^ {2} + a \ lambda + b \ right) = 0}
Deoarece exponențialul nu dispare niciodată, această ecuație dispare dacă și numai dacă:
- {\ displaystyle \ lambda ^ {2} + a \ lambda + b = 0}
Dacă rădăcinile sale sunt reale și distincte, asta este {\ displaystyle \ lambda _ {1} \ neq \ lambda _ {2} \ in \ mathbb {R}} , atunci soluția este așa:
- {\ displaystyle y = c_ {1} e ^ {\ lambda _ {1} x} + c_ {2} e ^ {\ lambda _ {2} x}}
dacă sunt reale și coincidente, adică {\ displaystyle \ lambda _ {1} = \ lambda _ {2}} , atunci soluția este așa:
- {\ displaystyle y = \ left (c_ {1} + c_ {2} \ cdot x \ right) și ^ {\ lambda _ {1} x}}
în timp ce dacă sunt complexe și conjugate, adică {\ displaystyle \ lambda _ {1,2} = \ alpha \ pm i \ beta \ in \ mathbb {C}} , atunci partea reală și imaginară pot fi considerate separat:
- {\ displaystyle y = e ^ {\ alpha x} \ left (c_ {1} \ cos \ beta x + c_ {2} \ sin \ beta x \ right)}
Ecuația completă
Ecuația completă are forma:
- {\ displaystyle y '' + a \ cdot y '+ b \ cdot y = f (x)}
Pentru a determina soluțiile, este suficient să adăugați o soluție specială a ecuației neomogene la soluția generică a ecuației omogene asociate. O astfel de soluție particulară poate fi găsită cu metoda variației constante sau luând în considerare unele cazuri particulare:
- De sine {\ displaystyle f (x) = P (x)} , unde este {\ displaystyle P (x)} este un polinom de grad {\ displaystyle m} , se caută o soluție specială de acest tip{\ displaystyle u (x) = P_ {1} (x)} , unde este {\ displaystyle P_ {1} (x)} este un polinom formal de același grad {\ displaystyle m} . Dacă totuși {\ displaystyle {\ lambda} = 0} este o soluție (simplă) a ecuației omogene asociate, atunci trebuie să căutăm o soluție de acest tip {\ displaystyle u (x) = x \ cdot P_ {1} (x)} .
- Considera {\ displaystyle f (x) = P (x) \ cdot e ^ {\ alpha x}} , unde este {\ displaystyle P (x)} este un polinom de grad {\ displaystyle m} . De sine {\ displaystyle \ alpha} nu este o rădăcină a ecuației {\ displaystyle \ lambda ^ {2} + a \ lambda + b = 0} apoi căutăm o soluție specială, cum ar fi:
- {\ displaystyle u (x) = P_ {1} (x) \ cdot e ^ {\ alpha x}}
- unde este {\ displaystyle P_ {1}} este un polinom formal de același grad {\ displaystyle m} . Dacă în schimb {\ displaystyle \ alpha} coincide cu una dintre rădăcini, apoi se caută o soluție specială, cum ar fi:
- {\ displaystyle u (x) = x ^ {r} \ cdot P_ {1} (x) \ cdot e ^ {\ alpha x}}
- unde este {\ displaystyle r} este multiplicitatea rădăcinii {\ displaystyle \ alpha} .
- Consideră-se pe ei înșiși {\ displaystyle f (x) = A \ cos (\ beta x)} , sau {\ displaystyle f (x) = A \ sin (\ beta x)} sau din nou {\ displaystyle f (x) = A \ cos (\ beta x) + B \ sin (\ beta x)} , unde este {\ displaystyle A, B} li se dau constante. În acest caz, dacă {\ displaystyle i \ beta} nu este o rădăcină a ecuației {\ displaystyle \ lambda ^ {2} + a \ lambda + b = 0} apoi căutăm o soluție specială, cum ar fi:
- {\ displaystyle u (x) = C \ cos (\ beta x) + D \ sin (\ beta x)}
- unde este {\ displaystyle C} Și {\ displaystyle D} sunt constante de determinat. Dacă nu, căutăm o soluție precum:
- {\ displaystyle u (x) = x \ cdot \ left (C \ cos (\ beta x) + D \ sin (\ beta x) \ right)}
- Consideră-se pe ei înșiși {\ displaystyle f (x) = P (x) \ cdot \ cos (\ beta x)} , sau {\ displaystyle f (x) = Q (x) \ cdot \ sin (\ beta x)} sau din nou {\ displaystyle f (x) = P (x) \ cdot \ cos (\ beta x) + Q (x) \ cdot \ sin (\ beta x)} , unde este {\ displaystyle P (x), Q (x)} sunt polinoame. În acest caz, dacă {\ displaystyle i \ beta} nu este o rădăcină a ecuației {\ displaystyle \ lambda ^ {2} + a \ lambda + b = 0} apoi căutăm o soluție specială, cum ar fi:
- {\ displaystyle u (x) = P_ {1} (x) \ cdot \ cos (\ beta x) + Q_ {1} (x) \ cdot \ sin (\ beta x)}
- unde este {\ displaystyle P_ {1} (x)} Și {\ displaystyle Q_ {1} (x)} sunt polinoame de grad respectiv egale cu cele ale lui {\ displaystyle P (x)} Și {\ displaystyle Q (x)} . Dacă nu ( {\ displaystyle i \ beta} este rădăcina multiplicității {\ displaystyle r} ecuației caracteristice), căutăm o soluție de tipul:
- {\ displaystyle u (x) = x ^ {r} \ cdot \ left (P_ {1} (x) \ cdot \ cos (\ beta x) + Q_ {1} (x) \ cdot \ sin (\ beta x ) \ dreapta)}
- De sine {\ displaystyle f (x) = f_ {1} (x) + f_ {2} (x) + ... + f_ {m} (x)} , pentru liniaritatea ecuației poate fi rezolvată separat:
- {\ displaystyle y '' + ay '+ by = f_ {i} (x) \ qquad i = 1, ..., m}
- și apoi adăugați {\ displaystyle u_ {i} (x)} soluții:
- {\ displaystyle u = u_ {1} + u_ {2} + \ cdots + u_ {m}}
Ecuații cu coeficienți variabili
Ecuația omogenă asociată are forma:
- {\ displaystyle y '' + a (x) \ cdot y '+ b (x) \ cdot y = 0}
unde este {\ displaystyle a (x)} Și {\ displaystyle b (x)} sunt funcții continue într-un interval al axei reale. Rezoluția sa constă în căutarea unei soluții de formă:
- {\ displaystyle y = c_ {1} y_ {1} + c_ {2} y_ {2}}
Având în vedere ecuația completă:
- {\ displaystyle y '' + a (x) \ cdot y '+ b (x) \ cdot y = f (x)}
în acest caz se poate folosi metoda variației constante . Căutăm o soluție de același tip cu cea a omogenului, considerând constantele ca funcții:
- {\ displaystyle u (x) = c_ {1} (x) y_ {1} + c_ {2} (x) y_ {2}}
unde este {\ displaystyle y_ {1}} Și {\ displaystyle y_ {2}} sunt două soluții independente ale ecuației omogene asociate (două soluții sunt independente una de cealaltă dacă raportul lor nu este constant). De cand {\ displaystyle y_ {1}} Și {\ displaystyle y_ {2}} sunt cunoscute și funcții {\ displaystyle c_ {1}} Și {\ displaystyle c_ {2}} sunt necunoscute, acestea din urmă trebuie determinate în așa fel încât {\ displaystyle u (x)} satisfaceți ecuația completă. Mai mult, deoarece există două funcții care urmează să fie determinate, o a doua condiție poate fi impusă {\ displaystyle c_ {1}} Și {\ displaystyle c_ {2}} după plac. Alege:
- {\ displaystyle c_ {1} ^ {'} (x) y_ {1} + c_ {2} ^ {'} (x) y_ {2} = 0}
astfel încât derivarea {\ displaystyle u (x)} de două ori și folosind această relație avem:
- {\ displaystyle {\ begin {cases} u '(x) = c_ {1} (x) y_ {1} ^ {'} + c_ {2} (x) y_ {2} ^ {'} \\ u' '(x) = c_ {1} (x) y_ {1} ^ {' '} + c_ {2} (x) y_ {2} ^ {' '} + c_ {1} ^ {'} (x) y_ {1} ^ {'} + c_ {2} ^ {'} (x) y_ {2} ^ {'} \ end {cases}}}
Înlocuind în ecuația completă obținem:
- {\ displaystyle c_ {1} (x) \ cdot \ left [y_ {1} ^ {''} + a (x) y_ {1} ^ {'} + b (x) y_ {1} \ right] + c_ {2} (x) \ left [y_ {2} ^ {''} + a (x) y_ {2} ^ {'} + b (x) y_ {2} \ right] + c_ {1} ^ {'} (x) y_ {1} ^ {'} + c_ {2} ^ {'} (x) y_ {2} ^ {'} = f (x)}
Este un sistem de necunoscute {\ displaystyle c_ {1} ^ {'} (x)} Și {\ displaystyle c_ {2} ^ {'} (x)} :
- {\ displaystyle {\ begin {cases} c_ {1} ^ {'} (x) y_ {1} + c_ {2} ^ {'} (x) y_ {2} = 0 \\ c_ {1} ^ { '} (x) y_ {1} ^ {'} + c_ {2} ^ {'} (x) y_ {2} ^ {'} = f (x) \ end {cases}}}
Odată obținut {\ displaystyle c_ {1} ^ {'} (x)} Și {\ displaystyle c_ {2} ^ {'} (x)} (se arată că acest lucru este întotdeauna fezabil având în vedere independența soluțiilor {\ displaystyle y_ {1}} Și {\ displaystyle y_ {2}} ), sunt obținute {\ displaystyle c_ {1} (x)} Și {\ displaystyle c_ {2} (x)} . În cele din urmă, soluția este:
- {\ displaystyle u (x) = c_ {1} (x) y_ {1} + c_ {2} (x) y_ {2}}
iar cea completă este:
- {\ displaystyle Y (x) = y + u (x)}
Notă
Bibliografie
- Arfken, G. „O a doua soluție”. §8.6 în Metode matematice pentru fizicieni , ed. A III-a. Orlando, FL: Academic Press, pp. 467-480, 1985.
- Boyce, WE și DiPrima, RC Ecuații diferențiale elementare și probleme cu valoarea limită , ediția a IV-a. New York: Wiley, 1986.
- Morse, PM și Feshbach, H. Metode de fizică teoretică, Partea I. New York: McGraw-Hill, pp. 667–674, 1953.
Elemente conexe