to całka szczególna równania niejednorodnego (7.1) jest postaci.
b e®*, gdy a nie jest pierwiastkiem równania charakterystycznego dla równania jednorodnego (7.2),
*1 'X ~ xkb e®*, gdy a jest pierwiastkiem o krotności k ^ 1 równania charakterystycznego dla równania jednorodnego (7.2),
gdzie współczynnik b należy wyznaczyć.
6° Jeżeli prawa strona f(x) równania (7.1) jest postaci:
gdzie: P„(x) i Q„(x) są wielomianami, z których jeden jest stopnia n, a drugi co najwyżej stopnia n, to całka szczególna yi(x) równania niejednorodnego (7.1) jest postaci:
yi (*) =
P„(x) cos + Sn (x) sin j3x, gdy liczba zespolona z = /?/ nie jest pierwiastkiem równania charakterystycznego dla równania (7.2), x? (Rn (x) cos Px + Sn(x) sin px), gdy liczba zespolona z = (U jest pierwiastkiem o krotności k ^ 1 równania charakterystycznego dla równania (7.2),
gdzie: R„(x) i Sn(x) są pewnymi wielomianami o własnościach analogicznych, jak Pn(x) i Q„(x), których współczynniki należy wyznaczyć.
" 'McTbilii druga (wariacji stałych dowolnych). Jeżeli prawa strona równania liniowego niejednorodnego (7.1) jest takiej postaci, że nie potrafimy odgadnąć postaci całki szczególnej równania niejednorodnego (7.1), stosujemy metodę ogólną wariacji stałych dowolnych. Metoda ta została podana w jednym z poprzednich punktów (por. 20.4) dla równania różniczkowego rzędu drugiego. Dla równania liniowego rc-tego (n ^ 3) rzędu (7.1) metoda ta przedstawia .się następująco:
1° Znajdujemy najpierw całkę ogólną równania jednorodnego (7.2) odpowiadającego równaniu (7.1) Przypuśćmy, że jest ona postaci:
(7-3) y = Ciyi(x) 4 C2y2(x) 4 ... 4 Cnyn(x).
2° Zakładamy następnie, że stałe Ck (k — 1,...,??) występujące po prawej stronie (7.3) są funkcjami zmiennej x, innymi słowy uzmienniamy stałe. Wzór (7.3) przyjmuje wówczas postać:
(7-4) y = C1 (x)y1 (x) 4 C2 (x) y2 (x) 4 ... 4 C„ (x)yn (x).
3° Nieznane funkcje Ck(x) (k = 1,...,»), a właściwie ich pochodne Ćk (x) (k = 1,..., n), wyznaczamy z następującego układu równań:
C[ (x) y i (x) 4 Q (x) y2 (x) 4 . • • 4 Q (*) yn (x) = 0 |
Ci W yi(*) + Q(*) yi(x) +... + c;(*) = o
(7.5) CJ (*) ^ (*) +C5(*) )4'(*)+ ... + «(*)< W = 0 }
gdzie: łk’(.x) oznacza wyznacznik charakterystyczny układu (7.5), zwany wyznacznikiem Wrońskiego (wrońskianem):
^(*) =
y i (*) |
y »(*) |
A W |
*(*) |
(*) | |
w |
(x) |
■ • • J^c*) |
4° Rozwiązując układ równań (7.5) względem C'k{x) (k = 1, . .n), otrzymujemy: (7.6)
Wk (x)
Ck (x) = T -rV , k =
łk(;t)
oraz gdzie (k=l,...,ri) oznaczają wyznaczniki, które otrzymujemy z wyznacznika W, zastępując w nim kolejno kolumnę pierwszą, drugą itd. kolumną wolnych wyrazów układu (7.5).
5° Całkujemy teraz równania (7.6). Mamy wtedy:
(7.7)
Wk(x)
W(x)
dx, k —
gdzie opuszczono stałe całkowania, gdyż chodzi o całkę szczególną.
6° Podstawiając (7.7) do (7.4), otrzymujemy szukaną całkę szczególną równania niejednorodnego (7.1).
Uwaga 1. Przy wyznaczaniu całki szczególnej równania liniowego niejednorodnego (7.1) bardzo praktyczne i pożyteczne okazuje się następujące twierdzenie.
Twierdzenie 2. Jeżeli ^(.y) oraz y2(*) są odpowiednio całkami szczególnymi równań:
y(n) + -f a2y(,,-2) + ... -f an-x y' + an y =fl (x),
oraz
y{n) + a1y(n-1) + a2yin~2) -f . .. -r an-x y' + an y — f2 (x), to funkcja }'i(x) + y2(A") jest całką szczególną równania:
y
,(n—i)
-f anmm 1 / + y = A (x) + f2 (x). dx
Uwaga 2. Jeżeli obok siły sprężystej — kx oraz siły — / • —j-, wynikającej z oporu
ośrodka, działa na punkt materialny m dodatkowa siła zewnętrzna /(?) zależna od czasu, to równanie różniczkowe ruchu przyjmuje postać:
d2x dx
m -kx -
czyli: (7.8)
dt2
d2x n dt2
dt
Równanie (7.8) obrazujące drgania mechaniczne wymuszone jest, jak widać, równaniem różniczkowym liniowym niejednorodnym rzędu drugiego o współczynnikach stałych.
265