Równania różniczkowe, budownictwo, III semestr, Analiza matematyczna 3, Matematyka, Matma2odinnegozioma, Matma


1. RÓWNANIA RÓŻNICZKOWE RZĘDU PIERWSZEGO

1.1 POJĘCIA WSTĘPNE

Def. 1.1.1 (równanie różniczkowe zwyczajne pierwszego rzędu)

Równaniem różniczkowym zwyczajnym rzędu pierwszego nazywamy równanie postaci:

(R) 0x01 graphic
.

Uwaga. Będziemy się również posługiwali tzw. formą różniczkową równania różniczkowego, czyli równaniem postaci

0x01 graphic
.

Jednak najogólniejszą formą równania różniczkowego rzędu pierwszego jest równanie postaci

0x01 graphic
.

Inaczej mówiąc równanie różniczkowe rzędu pierwszego jest zależnością między funkcją niewiadomą, zmienną niezależną i pierwszą pochodną funkcji niewiadomej.

Def. 1.1.2 (rozwiązanie równania różniczkowego)

Funkcję y(t) nazywamy rozwiązaniem równania różniczkowego (R) na przedziale (a,b), jeżeli jest ona różniczkowalna na tym przedziale oraz zamienia to równanie w tożsamość

0x01 graphic
,

prawdziwą dla wszystkich t (a,b). Wykres (rys. 1) rozwiązania równania różniczkowego nazywamy jego krzywą całkową.

0x01 graphic

Rys. 1 Krzywa całkowa

Uwaga. Analogicznie określamy rozwiązania na przedziałach postaci: [a,b), (a,b], [a,b], (-∞,b), (-∞,b], [a,∞), (-∞,∞). Przy czym w przypadku, gdy rozwiązanie określone jest na przedziale domkniętym z jednego lub obu końców, przez jego pochodną na końcu przedziału rozumiemy odpowiednią pochodną jednostronną. Rozwiązanie równania różniczkowego zadane w postaci uwikłanej

0x01 graphic

nazywamy całką tego równania. Ponieważ każde rozwiązanie jest całką (niekoniecznie odwrotnie), więc często w odniesieniu do rozwiązań używa się także terminu całka. Stąd mówimy zwyczajowo scałkować równanie różniczkowe, zamiast rozwiązać równanie.

Def. 1.1.3 (zagadnienie początkowe)

Równanie różniczkowe (R) oraz warunek

(W) 0x01 graphic

nazywamy zagadnieniem początkowym lub zagadnieniem Cauchy'ego.

Uwaga. Zagadnienie początkowe będziemy zapisywali krótko

(RW) 0x01 graphic

Przy czym liczby t0 i y0 nazywamy wartościami początkowymi, a warunek (W) nazywamy warunkiem początkowym.

Def. 1.1.4 (rozwiązanie zagadnienia początkowego)

Funkcja y(t) jest rozwiązaniem zagadnienia początkowego (RW), jeżeli jest rozwiązaniem równania (R) na pewnym przedziale zawierającym t0 i spełnia warunek (W).

0x01 graphic

Rys. 2

Uwaga. W interpretacji geometrycznej, rozwiązanie zagadnienia początkowego, polega na znalezieniu wśród wszystkich krzywych całkowych równania (R) tej, która przechodzi przez punkt (t0,y0) (rys.2). Jednak zagadnienie to niekoniecznie musi mieć jednoznaczne rozwiązanie. Może istnieć więcej niż jedno rozwiązanie danego zagadnienia początkowego. Istnienie rozwiązań zagadnienia początkowego oraz ich jednoznaczność jest jednym z głównych problemów teorii równań różniczkowych zwyczajnych.

Tw. 1.1.5 (o istnieniu i jednoznaczności rozwiązań równania (R))

Niech funkcja f(t,y) oraz jej pochodna cząstkowa 0x01 graphic
będą określone i ciągłe na obszarze domkniętym D R2. Wtedy dla każdego punktu (t0,y0) D, zagadnienie początkowe (RW) ma dokładnie jedno rozwiązanie.

Uwaga. Inaczej mówiąc przy dowolnych wartościach początkowych wybranych z obszaru D istnieje zawsze rozwiązanie zagadnienia początkowego (RW). Co więcej, jeżeli dane są dwa rozwiązania o tych samych wartościach początkowych (W), przy czym każde z rozwiązań określone jest na pewnym przedziale zawierającym punkt t0, to rozwiązania te pokrywają się na wspólnej części rozważanych przedziałów.

Def. 1.1.6 (rozwiązanie ogólne i szczególne równania różniczkowego)

Rodzinę funkcji

0x01 graphic
,

zależnych od parametru rzeczywistego C, nazywamy rozwiązaniem ogólnym równania (R), jeżeli:

1. każda funkcja tej rodziny jest jego rozwiązaniem,

2. dla każdego warunku początkowego 0x01 graphic
, dla którego rozwiązanie istnieje i jest jednoznaczne można dobrać stałą C tak, aby 0x01 graphic
.

Każdą funkcję otrzymaną z rozwiązania ogólnego równania (R) przy ustalonej wartości parametru C nazywamy rozwiązaniem szczególnym tego równania.

Uwaga. Rozwiązanie zagadnienia początkowego, jeżeli istnieje i jest jednoznaczne, jest rozwiązaniem szczególnym. W praktyce znajomość rozwiązania ogólnego jest bardzo dogodna, gdyż przez odpowiedni dobór parametru C można otrzymać rozwiązanie zagadnienia początkowego. Rozwiązanie ogólne równania różniczkowego, określone w postaci uwikłanej

0x01 graphic
,

nazywamy całką ogólną tego równania.

1.2 RÓWNANIA RÓŻNICZKOWE O ZMIENNYCH ROZDZIELONYCH

Def. 1.2.1 (równanie różniczkowe o zmiennych rozdzielonych)

Równaniem różniczkowym o zmiennych rozdzielonych nazywamy równanie postaci

(S) 0x01 graphic
.

Uwaga. Zauważmy, że jeżeli h(y0) = 0 dla pewnego y0, to funkcja y(t) y0 jest jednym z rozwiązań powyższego równania. W formie różniczkowej o zmiennych rozdzielonych przyjmuje postać

0x01 graphic
.

Fakt 1.2.2 (całka ogólna równania (S))

Niech funkcje g(t) i h(y) będą ciągłe odpowiednio na przedziałach (a,b) i (c,d), przy czym h(y) 0 dla każdego y (c,d). Wtedy całka ogólna równania różniczkowego o zmiennych rozdzielonych (S) dana jest wzorem

0x01 graphic

Uwaga. Całki w powyższym wzorze rozumiane są jako dowolne, lecz ustalone funkcje pierwotne.

Tw. 1.2.3 (o istnieniu i jednoznaczności rozwiązań równania (S))

Niech funkcje g(t) i h(y) będą ciągłe odpowiednio na przedziałach (a,b) i (c,d), przy czym h(y) 0 dla każdego y (c,d). Wtedy dla każdego punktu (t0,y0), gdzie t0 (a,b), y0 (c,d), zagadnienie początkowe

0x01 graphic

ma dokładnie jedno rozwiązanie.

Uwaga. Inaczej mówiąc, przez każdy punkt (t0,y0) prostokąta (a,b)×(c,d) przechodzi dokładnie jedna krzywa całkowa (rys.3) równania y' = g(t)h(y).

0x01 graphic

Rys. 3

1.3 RÓWNANIA RÓŻNICZKOWE JEDNORODNE

Def. 1.3.1 (równanie różniczkowe jednorodne)

Równaniem różniczkowym jednorodnym nazywamy równanie postaci

(J) 0x01 graphic
.

Uwaga. Jeżeli f(u) u, to równanie jednorodne przyjmuje postać 0x01 graphic
i całkuje się przy pomocy metody rozdzielonych zmiennych. Jego rozwiązanie ogólne dane jest wtedy wzorem y(t) = Ct, gdzie C R. jeżeli f(u0) = u0 dla pewnego u0, to jedynym rozwiązaniem równania (J) jest funkcja y(t) = u0t.

Fakt 1.3.2 (o zamianie zmiennych w równaniu jednorodnym)

Równanie jednorodne (J) przez zamianę zmiennych

y=ut

sprowadza się do równania o zmiennych rozdzielonych postaci

tu' = f(u) - u.

Tw. 1.3.3 (o istnieniu i jednoznaczności rozwiązań równania (J))

Niech funkcja f(u) będzie ciągła na przedziale (a,b) i niech spełnia na nim warunek f(u) u. Wtedy dla każdego punktu (t0,y0) takiego, że 0x01 graphic
zagadnienie początkowe

0x01 graphic

ma dokładnie jedno rozwiązanie.

0x01 graphic

Rys. 4

Uwaga. Inaczej mówiąc, przez każdy punkt (t0,y0) obszaru

0x01 graphic
,

przechodzi dokładnie jedna krzywa całkowa (rys. 4) równania (J).

1.4 RÓWNANIA RÓŻNICZKOWE LINIOWE

Def. 1.4.1 (równanie różniczkowe liniowe)

Równaniem różniczkowym liniowym pierwszego rzędu nazywamy równanie postaci

(L) 0x01 graphic
.

Jeżeli q(t) 0, to równanie nazywamy niejednorodnym. W przypadku przeciwnym nazywamy je jednorodnym.

Uwaga. Równanie różniczkowe liniowe jednorodne jest szczególnym przypadkiem równania różniczkowego o zmiennych rozdzielonych y'=g(t)h(y), w którym przyjęto g(t) = - p(t), h(y) = y.

Tw. 1.4.2 (o istnieniu i jednoznaczności rozwiązań równania (L))

Niech funkcje p(t) i q(t) będą ciągłe na przedziale (a,b), gdzie - a < b . Wtedy dla każdego punktu (t0,y0), gdzie t0 (a,b) oraz y0 R, zagadnienie początkowe

0x01 graphic

ma dokładnie jedno rozwiązanie. Rozwiązanie to jest określone na przedziale (a,b).

0x01 graphic

Rys. 5

Uwaga. Inaczej mówiąc, przez każdy punkt pasa (a,b) × R przechodzi dokładnie jedna krzywa całkowa (rys. 5) równania różniczkowego liniowego.

1.5 RÓWNANIE RÓŻNICZKOWE BERNOULLIEGO

Def. 1.5.1 (równanie różniczkowe Bernoulliego)

Równaniem różniczkowym Bernoulliego nazywamy równanie postaci

0x01 graphic
,

gdzie rR-{0,1}.

Uwaga. Gdyby dopuścić r = 0, to równanie Bernoulliego byłoby równaniem różniczkowym liniowym niejednorodnym. Natomiast dla r = 1, równanie to byłoby równaniem różniczkowym liniowym jednorodnym. Zauważmy jeszcze, że dla r > 0 funkcja y(t) 0 jest zawsze jednym z rozwiązań równania Bernoilliego.

Fakt 1.5.2 (sprowadzanie równania Bernoulliego do równania liniowego)

Równanie Bernoulliego

0x01 graphic
,

gdzie r 0, 1, przez zamianę zmiennych 0x01 graphic
sprowadza się do równania liniowego niejednorodnego postaci

0x01 graphic
.

1. 6 KRZYWE ORTOGONALNE

Def. 1.6.1 (równanie rodziny krzywych)

Jeżeli równanie

0x01 graphic

dla każdej wartości parametru C z pewnego przedziału określa krzywą, to nazywamy je równaniem rodziny krzywych (rys. 6).

Def. 1.6.1 (rodzina krzywych ortogonalnych)

Mówimy, że rodziny krzywych (t,y,C) = 0, (t,y,C) = 0 są ortogonalne, jeżeli w każdym punkcie przecięcia krzywych z obu rodzin, krzywe te tworzą między sobą kąt prosty (rys. 7).

0x01 graphic

0x01 graphic

Rys. 6

Rys. 7

Fakt 1.6.2 (równanie różniczkowe rodziny krzywych ortogonalnych)

Jeżeli F(t,y,y') = 0 jest równaniem różniczkowym rodziny krzywych, to równanie różniczkowe rodziny krzywych ortogonal­nych ma postać

0x01 graphic
.

1.7 POJĘCIA WSTĘPNE DLA RÓWNAŃ RÓŻNICZKOWYCH WYŻSZYCH RZĘDÓW

Def. 1.7.1 (równanie różniczkowe zwyczajne n-tego rzędu)

Równaniem różniczkowym zwyczajnym rzędu n nazywamy równanie postaci

(R) 0x01 graphic
.

Uwaga. Najogólniejszą formą równania różniczkowego rzędu n jest wyrażenie postaci 0x01 graphic
.

Def. 1.7.2 (rozwiązanie równania różniczkowego n-tego rzędu)

Funkcję y(t), różniczkowalną n-krotnie na przedziale otwartym (a,b), nazywamy rozwiązaniem równania różniczkowego (R) na tym przedziale, jeżeli zamienia to równanie w tożsamość

0x01 graphic

prawdziwą dla wszystkich t należących do przedziału t (a,b). Wykres rozwiązania równania różniczkowego nazywamy jego krzywą całkową.

Def. 1.7.3 (zagadnienie początkowe)

Równanie różniczkowe (R) oraz warunki

(W) 0x01 graphic
,

nazywamy zagadnieniem początkowym lub zagadnieniem Cauchy'ego.

Uwaga. Zagadnienie początkowe będziemy zapisywali krótko

(RW) 0x01 graphic
,

przy czym liczby t0 i y0, y1, ..., yn-1 nazywamy wartościami początkowymi, a warunek (W) nazywamy warunkiem początkowym.

Def. 1.7.4 (rozwiązanie zagadnienia początkowego)

Funkcja y(t) jest rozwiązaniem zagadnienia początkowego (RW) jeżeli jest rozwiązaniem równania (R) na przedziale zawierającym punkt t0 i jeżeli spełnia warunki (W).

Def. 1.7.5 (rozwiązanie ogólne i szczególne równania różniczkowego)

Rodzinę funkcji

0x01 graphic

zależnych od n rzeczywistych parametrów C1, C2, ..., Cn nazywamy rozwiązaniem ogólnym równania 0x01 graphic
, jeżeli:

1. każda funkcja tej rodziny jest rozwiązaniem tego równania,

2. dla każdego układu warunków początkowych 0x01 graphic
, 0x01 graphic
, ..., 0x01 graphic
, dla którego rozwiązanie istnieje i jest jednoznaczne, można dobrać stałe C1, C2, ..., Cn tak, aby

0x01 graphic
.

Każdą funkcję otrzymaną z rozwiązania ogólnego równania (R) przy ustalonych wartościach parametrów C1, C2, ..., Cn nazywamy rozwiązaniem szczególnym tego równania.

Uwaga. Rozwiązanie zagadnienia początkowego, jeżeli istnieje i jest jednoznaczne, jest rozwiązaniem szczególnym. W praktyce znajomość rozwiązania ogólnego jest bardzo dogodna, gdyż przez odpowiedni dobór parametrów C1, C2, ..., Cn można otrzymać rozwiązanie zagadnienia początkowego.

1.8 RÓWNANIA RZĘDU DRUGIEGO SPROWADZALNE DO RÓWNAŃ RZĘDU PIERWSZEGO

Fakt 1.8.1 (równanie różniczkowe postaci y'' = f(t,y'))

Równanie różniczkowe drugiego rzędu postaci

0x01 graphic

przez podstawienie y' = u sprowadza się do równania różniczkowego rzędu pierwszego postaci

0x01 graphic
.

Fakt 1.8.2 (równanie różniczkowe postaci y'' = f(y,y'))

Równanie różniczkowe drugiego rzędu postaci

0x01 graphic

przez podstawienie y = q(y) sprowadza się do równania różniczkowego rzędu pierwszego postaci

0x01 graphic
.

2. RÓWNANIA RÓŻNICZKOWE LINIOWE WYŻSZYCH RZĘDÓW

2.1 POJĘCIA WSTĘPNE

Def. 2.1.1 (Równanie różniczkowe liniowe rzędu n)

Równaniem różniczkowym liniowym rzędu n nazywamy równanie postaci

(Ln) 0x01 graphic
.

Uwaga. Jeżeli n = 2, to będziemy pisali p(t) i q(t) zamiast odpowiednio p1(t) i p2(t). Równanie różniczkowe liniowe rzędu drugiego ma wtedy postać

(L2) 0x01 graphic
.

Tw. 2.1.2 (o istnieniu i jednoznaczności rozwiązań dla równania (Ln))

Niech funkcje p1(t), p2(t), ..., pn(t) i h(t) będą ciągłe na przedziale (a,b). Wtedy dla każdego punktu (t0, y0, y1, ..., yn-1) (a,b) × Rn zagadnienie początkowe

0x01 graphic
,

0x01 graphic

ma dokładnie jedno rozwiązanie. Rozwiązanie to jest określone na przedziale (a,b).

2.2 RÓWNANIA RÓŻNICZKOWE LINIOWE JEDNORODNE

Def. 2.2.1 (równanie różniczkowe liniowe jednorodne)

Równaniem różniczkowym jednorodnym rzędu n nazywamy równanie postaci

(LJn) 0x01 graphic
.

Uwaga. Dla równania drugiego rzędu w tym przypadku (LJ2) piszemy p(t) i q(t) zamiast odpowiednio p1(t) i p2(t), czyli

(LJ2) 0x01 graphic
.

Fakt 2.2.2 (o kombinacji liniowej rozwiązań równania jednorodnego)

Niech ϕ(t), (t) będą rozwiązaniami równania jednorodnego (LJn). Wtedy dla dowolnych stałych , funkcja 0x01 graphic
jest także rozwiązaniem tego równania.

Uwaga. Inaczej mówiąc, dowolna kombinacja liniowa rozwiązań równania różniczkowego liniowego jednorodnego jest również rozwiązaniem tego równania.

Def. 2.2.3 (układ fundamentalny równania (LJn))

Układ n rozwiązań (y1(t), y2(t), ..., yn(t)) równania jednorodnego (LJn) określonych na przedziale (a,b) nazywamy układem fundamentalnym tego równania na tym przedziale, jeżeli dla każdego t (a,b) spełniony jest warunek

0x01 graphic
.

Uwaga. Powyższy wyznacznik oznaczamy przez W(y1(t), y2(t), ..., yn(t)) i nazywamy wrońskianem układu funkcji (y1(t), y2(t), ..., yn(t)).

Fakt 2.2.4 (Wzór Liouville'a)

Niech (y1(t), y2(t), ..., yn(t)) będzie układem n rozwiązań równania jednorodnego (LJn) określonych na przedziale (a,b). Wtedy ich wrońskian W(t) = W(y1(t), y2(t), ..., yn(t)) spełnia warunek

0x01 graphic
,

gdzie t0 jest dowolnym punktem z przedziału (a,b).

Fakt 2.2.5 (o postaci rozwiązania ogólnego równania liniowego jednorodnego)

Niech (y1(t), y2(t), ..., yn(t)) będzie układem fundamentalnym równania jednorodnego (LJn). Wtedy rozwiązanie ogólne tego równania dane jest wzorem

0x01 graphic
,

gdzie C1, C2, ..., Cn są dowolnymi stałymi rzeczywistymi.

Uwaga. Przypominamy, że z rozwiązania ogólnego przez odpowiedni dobór stałych C1, C2, ..., Cn można otrzymać rozwiązanie każdego zagadnienia początkowego. Rozwiązanie ogólne równania liniowego jednorodnego będziemy oznaczali symbolem yRORLJ.

Fakt 2.2.6 (o obniżaniu rzędu równania liniowego jednorodnego)

Niech ϕ(t) będzie różnym od zera rozwiązaniem równania różniczkowego liniowego jednorodnego rzędu n. Wtedy przez podstawienie

0x01 graphic

równanie to sprowadza się do równania liniowego jednorodnego rzędu n - 1 (względem nowej zmiennej z).

Uwaga. Jeżeli znamy jedno rozwiązanie równania liniowego jednorodnego rzędu drugiego, to znalezienie rozwiązania ogólnego tego równania sprowadza się do rozwiązania, otrzymanego przez podstawienie, równania liniowego jednorodnego rzędu pierwszego.

2.3 RÓWNANIA RÓŻNICZKOWE LINIOWE O STAŁYCH WSPÓŁCZYNNIKACH

Def. 2.3.1 (równania różniczkowe liniowe o stałych współczynnikach)

Równaniem różniczkowym liniowym jednorodnym rzędu n o stałych współczynnikach nazywamy równanie postaci

(LSn) 0x01 graphic
,

gdzie p1, p2, ..., pn R.

Uwaga. Każde rozwiązanie równania różniczkowego liniowego o stałych współczynnikach (LSn) jest określone na R. Przypominamy, że zgodnie z poprzednią umową równanie różniczkowe jednorodne rzędu drugiego o stałych współczynnikach ma postać

(LS2) 0x01 graphic
.

Def. 2.3.2 (wielomian i równanie charakterystyczne)

Równanie postaci

0x01 graphic

nazywamy równaniem charakterystycznym równania różniczkowego liniowego o stałych współczynnikach (LSn). Natomiast wielomian

0x01 graphic

nazywamy wielomianem charakterystycznym tego równania.

Fakt 2.3.3 (o postaci układu fundamentalnego równania (LSn))

Niech 1, ..., s będą rzeczywistymi pierwiastkami o krotnościach odpowiednio k1, ..., ks i niech 0x01 graphic
, gdzie 0x01 graphic
będą zespolonymi pierwiastkami o krotnościach odpowiednio l1, ..., lm wielomianu charakterystycznego w() równania liniowego o stałych współczynnikach (LSn), przy czym k1+...+ks+2(l1+...+lm) = n. Wtedy układ fundamentalny tego równania tworzą funkcje:

0x01 graphic

Uwaga. W tym zestawieniu:

  1. pierwiastkowi rzeczywistemu jednokrotnemu odpowiada funkcja 0x01 graphic
    ;

  2. pierwiastkowi rzeczywistemu k-krotnemu odpowiada k funkcji 0x01 graphic
    ;

  3. parze jednokrotnych sprzężonych ze sobą pierwiastków zespolonych 0x01 graphic
    , 0x01 graphic
    odpowiada para funkcji

0x01 graphic
;

  1. parze l-krotnych sprzężonych ze sobą pierwiastków zespolonych 0x01 graphic
    , 0x01 graphic
    odpowiada l par funkcji

0x01 graphic
.

Fakt 2.3.4 (algorytm rozwiązywania równań o stałych współczynnikach)

Równanie jednorodne o stałych współczynnikach

(LSn) 0x01 graphic

0x01 graphic

Równanie charakterystyczne

0x01 graphic

0x01 graphic

Pierwiastki równania charakterystycznego

0x01 graphic

0x01 graphic

Układ fundamentalny równania (LSn)

0x01 graphic

0x01 graphic

Rozwiązanie ogólne równania (LSn)

0x01 graphic

2.4 RÓWNANIE EULERA

Def. 2.4.1 (równanie Eulera)

Równaniem różniczkowym Eulera rzędu n nazywamy równanie postaci

(En) 0x01 graphic
,

gdzie p1, p2, ..., pn R.

Uwaga. Do równania Eulera (En) sprowadza się równanie różniczkowe postaci

0x01 graphic
.

Aby tego dokonać należy podstawić 0x01 graphic
.

Fakt 2.4.2 (o sprowadzaniu równania Eulera do równania liniowego)

Podstawienie 0x01 graphic
przy t > 0 lub 0x01 graphic
przy t < 0 sprowadza równanie Eulera do równania różniczkowego liniowego jednorodnego o stałych współczynnikach.

2.5 RÓWNANIA RÓŻNICZKOWE LINIOWE NIEJEDNORODNE

Def. 2.5.1 (równanie różniczkowe liniowe niejednorodne)

Równaniem różniczkowym liniowym niejednorodnym rzędu n nazywamy równanie postaci

(LNn) 0x01 graphic

gdzie h(t)0.

Uwaga. Przypominamy, że zgodnie z poprzednią umową równanie różniczkowe liniowe niejednorodne rzędu drugiego ma postać

(LN2) 0x01 graphic
.

Fakt 2.5.2 (o postaci rozwiązania ogólnego równania niejednorodnego)

Niech (y1(t), y2(t), ..., yn(t)) będzie układem fundamentalnym równania jednorodnego (LJn) i niech ϕ(t) będzie dowolnym rozwiązaniem równania niejednorodnego (LNn). Wtedy rozwiązanie ogólne równania niejednorodnego dane jest wzorem

0x01 graphic
,

gdzie C1, C2, ..., Cn są dowolnymi stałymi rzeczywistymi.

Uwaga. Rozwiązanie ogólne równania niejednorodnego (LNn) będziemy oznaczali symbolerm yRORLN. Zauważmy, że rozwiązanie to jest sumą rozwiązania ogólnego yRORLJ równania jednorodnego i dowolnego rozwiązania ϕ równania niejednorodnego, co można zapisać w postaci

yRORLN = yRORLJ + ϕ.

Każde rozwiązanie równania liniowego jest rozwiązaniem szczególnym tego równania.

2.6 METODA UZMIENNIANIA STAŁYCH

Fakt 2.6.1 (metoda uzmienniania stałych)

Niech (y1(t), y2(t), ..., yn(t)) będzie układem fundamentalnym równania jednorodnego (LJn). Wtedy funkcja

0x01 graphic
,

gdzie c1(t), c2(t), ...,cn(t) są dowolnymi rozwiązaniami układu równań

0x01 graphic
,

jest rozwiązaniem szczególnym równania niejednorodnego (LNn).

Uwaga. Powyższy układ równań względem niewiadomych c1'(t), c2'(t), ...,cn'(t) ma jednoznaczne rozwiązanie, gdyż jego wyznacznik jest wrońskianem układu fundamentalnego (y1(t), y2(t), ..., yn(t)) równania jednorodnego (LJn), który jest różny od zera.

2.7 METODA WSPÓŁCZYNNIKÓW NIEOZNACZONYCH - METODA PRZEWIDYWANIA

Def. 2.7.1 (stała kontrolna)

Niech funkcja h(t) ma postać

0x01 graphic
.

Stałą kontrolną tej funkcji nazywamy liczbę 0x01 graphic
.

Fakt 2.7.2 (o rozwiązaniach szczególnych równania liniowego niejednorodnego)

Niech prawa strona równania liniowego niejednorodnego o stałych współczynnikach

0x01 graphic

ma postać

0x01 graphic

i niech w() będzie wielomianem charakterystycznym tego równania, a 0x01 graphic
stałą kontrolną funkcji h(t). Wtedy

1. jeżeli σ nie jest pierwiastkiem wielomianu w(), to funkcja postaci

0x01 graphic

2. jeżeli σ jest s-krotnym pierwiastkiem wielomianu w(), to funkcja postaci

0x01 graphic

gdzie m = max{k,l}, a A1, ..., Am, B1, ..., Bm są odpowiednio dobranymi współczynnikami rzeczywistymi, jest rozwiązaniem szczególnym tego równania.

Fakt 2.7.3 (o superpozycji rozwiązań)

Niech funkcje ψ(t) i ϕ(t) będą rozwiązaniami odpowiednio równań

0x01 graphic

0x01 graphic
.

Wtedy ich suma ψ(t) + ϕ(t) jest rozwiązaniem równania

0x01 graphic
.

3. UKŁADY RÓWNAŃ RÓŻNICZKOWYCH

3.1 POJĘCIA WSTĘPNE

Def. 3.1.1 (układ równań różniczkowych)

Układem równań różniczkowych rzędu pierwszego nazywamy układ równań postaci

(U) 0x01 graphic
.

Uwaga. Jeżeli n = 2, to będziemy pisali x, y zamiast y1, y2 oraz f, g w miejsce f1, f2. Podobnie, jeżeli n = 3, to będziemy pisali x, y, z zamiast y1, y2, y3 oraz f, g, h w miejsce f1, f2, f3.

Def. 3.1.2 (rozwiązanie układu równań)

Ciąg funkcji (y1(t), y2(t), ..., yn(t)) określonych i różniczkowalnych na przedziale (a,b) nazywamy rozwiązaniem układu równań (U) na tym przedziale, jeżeli zamienia on wszystkie równania tego układu w tożsamości

0x01 graphic

na przedziale (a,b).

Uwaga. W notacji wektorowej układ równań różniczkowych (U) można zapisać w postaci

0x01 graphic
,

gdzie

0x01 graphic
, 0x01 graphic
, 0x01 graphic
.

Wtedy rozwiązanie układu równań (U) jest funkcją wektorową

0x01 graphic
.

Def. 3.1.3 (zagadnienie początkowe dla układu równań)

Układ równań różniczkowych (U) oraz układ warunków

(W) 0x01 graphic
,

nazywamy zagadnieniem początkowym lub zagadnieniem Cauchy'ego.

Uwaga. Liczby 0x01 graphic
i 0x01 graphic
nazywamy wartościami początkowymi, a warunek (W) nazywamy warunkiem początkowym. Używając notacji wektorowej zagadnienie początkowe można zapisać w postaci

(UW) 0x01 graphic
.

Def. 3.1.4 (rozwiązanie zagadnienia początkowego)

Ciąg funkcji (y1(t), y2(t), ..., yn(t)) jest rozwiązaniem zagadnienia początkowego (UW), jeżeli jest rozwiązaniem układu równań (U) na pewnym przedziale zawierającym punkt t0 i spełnia warunki (W).

Tw. 3.1.5 (o istnieniu i jednoznaczności rozwiązań układu (U))

Niech funkcje fi(t,y1,y2,...,yn), gdzie 1 ≤ in, wraz ze swoimi pochodnymi cząstkowymi 0x01 graphic
, gdzie 1 ≤ i, jn, będą określone i ciągłe na obszarze D Rn+1. Wtedy dla dowolnego punktu 0x01 graphic
zagadnienie początkowe (UW) ma dokładnie jedno rozwiązanie. Rozwiązanie to jest określone na pewnym otoczeniu punktu t0.

Def. 3.1.6 (rozwiązania ogólne i szczególne układu równań)

Rodzinę funkcji wektorowych

0x01 graphic
,

zależnych od parametrów rzeczywistych C1, C2, ..., Cn nazywamy rozwiązaniem ogólnym układu równań (U) jeżeli:

1. każda funkcja wektorowa z tej rodziny jest rozwiązaniem układu,

2. dla każdego układu warunków początkowych 0x01 graphic
, dla którego rozwiązanie istnieje i jest jednoznaczne, można dobrać stałe C1, C2, ..., Cn tak, aby 0x01 graphic
.

Każdą funkcję wektorową otrzymaną z rozwiązania ogólnego układu (U) przy ustalonych wartościach parametrów C1, C2, ..., Cn nazywamy rozwiązaniem szczególnym tego układu.

Uwaga. Rozwiązanie zagadnienia początkowego, jeżeli istnieje i jest jednoznaczne, jest rozwiązaniem szczególnym. W praktyce znajomość rozwiązania ogólnego jest bardzo dogodna, gdyż przez odpowiedni dobór parametrów C1, C2, ..., Cn można otrzymać rozwiązanie zagadnienia początkowego.

3.2 UKŁADY RÓWNAŃ RÓŻNICZKOWYCH LINIOWYCH

Def. 3.2.1 (układ równań różniczkowych liniowych)

Układem równań różniczkowych liniowych rzędu pierwszego nazywamy układ równań postaci

(UL) 0x01 graphic
.

Uwaga. W notacji wektorowej układ równań (UL) przyjmuje postać

0x01 graphic

lub krótko

0x01 graphic
.

Tw. 3.2.2 (o istnieniu i jednoznaczności rozwiązań dla układu (UL))

Niech funkcje aij(t), gdzie i,j = 1, 2, ..., n, oraz hi(t), gdzie i = 1, 2, ..., n, będą ciągłe na przedziale (a,b). Wtedy dla każdego punktu 0x01 graphic
zagadnienie początkowe

0x01 graphic
, 0x01 graphic

ma dokładnie jedno rozwiązanie. Rozwiązanie to określone jest na przedziale (a,b).

3.3 UKŁADY JEDNORODNE RÓWNAŃ RÓŻNICZKOWYCH LINIOWYCH

Def. 3.3.1 (układ jednorodny równań różniczkowych jednorodnych)

Układem jednorodnym równań różniczkowych liniowych rzędu pierwszego nazywamy układ równań postaci

(UJ) 0x01 graphic

Uwaga. W notacji wektorowej układ jednorodny (UJ) przyjmuje postać

0x01 graphic
.

Fakt 3.3.2 (o kombinacji liniowej rozwiązań układu jednorodnego)

Niech 0x01 graphic
, 0x01 graphic
będą rozwiązaniami układu jednorodnego (UJ). Wtedy dla dowolnych stałych α,β R funkcja wektorowa

0x01 graphic

jest również rozwiązaniem tego układu.

Def. 3.3.3 (układ fundamentalny układu jednorodnego (UJ))

Układ rozwiązań 0x01 graphic
układu jednorodnego (UJ) określonych na przedziale (a,b) nazywamy układem fundamentalnym tego układu na tym przedziale, jeżeli dla każdego t (a,b) spełniony jest warunek

0x01 graphic

Uwaga. Powyższy wyznacznik oznaczamy przez 0x01 graphic
i nazywamy wrońskianem układu funkcji wektorowych 0x01 graphic
.

Fakt 3.3.4 (o postaci rozwiązania ogólnego układu jednorodnego)

Niech 0x01 graphic
będzie układem fundamentalnym układu jednorodnego równań różniczkowych (UJ). Wtedy rozwiązanie ogólne tego układu równań dane jest wzorem

0x01 graphic

gdzie C1, C2, ..., Cn są dowolnymi stałymi rzeczywistymi.

Uwaga. Przypomnijmy, że z rozwiązania ogólnego przez odpowiedni dobór stałych C1, C2, ..., Cn można otrzymać rozwiązanie każdego zagadnienia początkowego. Rozwiązanie ogólne układu jednorodnego (UJ) będziemy oznaczali symbolem 0x01 graphic
.

3.4 UKŁADY RÓWNAŃ RÓŻNICZKOWYCH LINIOWYCH O STAŁYCH WSPÓŁCZYNNIKACH

Def. 3.4.1 (układ jednorodny równań liniowych o stałych współczynnikach)

Układem jednorodnym równań różniczkowych liniowych o stałych współczynnikach nazywamy układ równań postaci

(US) 0x01 graphic
,

gdzie aij R dla 1 ≤ i, j n.

Uwaga. W notacji wektorowej układ jednorodny równań różniczkowych liniowych o stałych współczynnikach (US) przyjmuje postać

0x01 graphic
.

Każde rozwiązanie układu (US) określone jest na R. Układy (US) można rozwiązywać metodą eliminacji.

Def. 3.4.2 (wielomian i równanie charakterystyczne macierzy)

Niech 0x01 graphic
będzie macierzą kwadratową stopnia n. Wielomianem charakterystycznym macierzy A nazywamy wielomian

0x01 graphic
,

gdzie I oznacza macierz jednostkową stopnia n. Równaniem charakterystycznym macierzy A nazywamy równanie

0x01 graphic
.

Def. 3.4.3 (wartości własne i wektory własne macierzy)

Niech 0x01 graphic
będzie rzeczywistą macierzą kwadratową stopnia n. Wartością własną macierzy A nazywamy każdy (rzeczywisty lub zespolony) pierwiastek wielomianu charakterystycznego tej macierzy, tj. liczbę λ spełniającą równanie

0x01 graphic
.

Niezerowy wektor 0x01 graphic
(o rzeczywistych lub zespolonych współrzędnych) nazywamy wektorem własnym macierzy A odpowiadającym wartości własnej λ (rzeczywistej lub zespolonej) tej macierzy, jeżeli spełnia warunek

0x01 graphic
.

Uwaga. Ostatnią równość można zapisać w postaci układu równań

0x01 graphic
.

Fakt 3.4.4 (o postaci układu fundamentalnego układu (US) - metoda Eulera I)

Jeżeli macierz A układu jednorodnego równań o stałych współczynnikach (US) ma n różnych rzeczywistych wartości własnych λ1, λ2, ..., λn, to jego układ fundamentalny składa się z funkcji wektorowych:

0x01 graphic
, 0x01 graphic
, ..., 0x01 graphic
,

gdzie 0x01 graphic
są wektorami własnymi odpowiadającymi wartościom własnym λ1, λ2, ..., λn.

Fakt 3.4.5 (o postaci układu fundamentalnego układu (US) - metoda Eulera II)

Jeżeli macierz A układu jednorodnego równań o stałych współczynnikach (US) ma 2k=n parami różnych zespolonych wartości własnych 0x01 graphic
, 0x01 graphic
, 0x01 graphic
, 0x01 graphic
, …, 0x01 graphic
, 0x01 graphic
, gdzie βj0 dla 1 j k, to jego układ fundamentalny składa się z funkcji wektorowych:

0x01 graphic
, 0x01 graphic
, 0x01 graphic
,

0x01 graphic
, ..., 0x01 graphic
, 0x01 graphic
,

gdzie 0x01 graphic
oznaczają wektory własne odpowiadające wartościom własnym λ1, λ2, ..., λn.

Fakt 3.4.6 (o postaci układu fundamentalnego układu (US) - metoda Eulera III)

Jeżeli macierz A układu jednorodnego równań o stałych współczynnikach (US) ma s parami różnych zespolonych wartości własnych λ1, λ2, ..., λs oraz 2k parami różnych zespolonych wartości własnych 0x01 graphic
, 0x01 graphic
, 0x01 graphic
, 0x01 graphic
, …, 0x01 graphic
, 0x01 graphic
, gdzie s + 2k = n, βj ≠ 0 dla 1 ≤ j k, to jego układ fundamentalny składa się z funkcji wektorowych:

0x01 graphic
, 0x01 graphic
, …, 0x01 graphic
,

0x01 graphic
, 0x01 graphic
,

…, 0x01 graphic
, 0x01 graphic
,

gdzie 0x01 graphic
oznaczają wektory własne odpowiadające rzeczywistym wartościom własnym λ1, λ2, ..., λs, a 0x01 graphic
oznaczają wektory własne odpowiadające zespolonym wartościom własnym λs+1, λs+2, ..., λs+k.

3.5 UKŁADY NIEJEDNORODNE RÓWNAŃ RÓŻNICZKOWYCH LINIOWYCH

Def. 3.5.1 (układ niejednorodny równań różniczkowych liniowych)

Układem niejednorodnym równań różniczkowych liniowych pierwszego rzędu nazywamy układ równań postaci

(UN) 0x01 graphic
,

gdzie hi(t) ≠ 0 dla pewnego 1 ≤ i n.

Uwaga. W notacji wektorowej układ niejednorodny (UN) przyjmuje postać

0x01 graphic
.

Fakt 3.5.2 (o postaci rozwiązania ogólnego układu niejednorodnego)

Niech 0x01 graphic
będzie układem fundamentalnym układu jednorodnego (UJ) i niech 0x01 graphic
będzie dowolnym rozwiązaniem układu niejednorodnego (UN). Wtedy rozwiązanie ogólne układu niejednorodnego dane jest wzorem

0x01 graphic
,

gdzie C1, C2, ..., Cn są dowolnymi stałymi rzeczywistymi.

Uwaga. Przypomnijmy, że z rozwiązania ogólnego przez odpowiedni dobór stałych C1, C2, ..., Cn można otrzymać rozwiązanie każdego zagadnienia początkowego. Rozwiązanie ogólne układu niejednorodnego (UN) będziemy oznaczali symbolem 0x01 graphic
. Rozwiązanie to ma taką samą postać jak rozwiązanie ogólne równania liniowego niejednorodnego (LN)

0x01 graphic
.

3.6 METODA UZMIENNIANIA STAŁYCH

Def. 3.6.1 (metoda uzmienniania stałych)

Niech dany będzie układ fundamentalny

0x01 graphic
, 0x01 graphic
, ..., 0x01 graphic
,

układu jednorodnego (UJ). Wtedy funkcja wektorowa

0x01 graphic
,

gdzie c1(t), c2(t), ..., cn(t) są dowolnymi rozwiązaniami układu równań

0x01 graphic
,

jest rozwiązaniem szczególnym układu niejednorodnego (UN).

Uwaga. Powyższy układ równań względem niewiadomych 0x01 graphic
ma jednoznaczne rozwiązanie, gdyż jego wyznacznik jest wrońskianem układu fundamentalnego 0x01 graphic
układu (UJ), który jest różny od zera.

4. ELEMENTY RACHUNKU OPERATOROWEGO

4.1 PRZEKSZTAŁCENIE LAPLACE'A

Def. 4.1.1 (transformata Laplace'a)

Niech funkcja f będzie określona na przedziale [0,∞). Transformatę Laplace'a funkcji f oznaczamy symbolem F(s) lub L{f(t)} i definiujemy wzorem

0x01 graphic
L{f(t)}0x01 graphic
,

gdzie s jest zmienną rzeczywistą. Funkcję F(s) nazywamy także L-transformatą lub obrazem funkcji f(t).

Fakt 4.1.2 (transformaty ważniejszych funkcji)

Funkcja

Transformata

1

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Fakt 4.1.3 (warunki wystarczające istnienia transformaty Laplace'a)

Niech

1. funkcja 0x01 graphic
ma skończoną liczbę punktów nieciągłości pierwszego rodzaju na każdym przedziale [0,T], gdzie T > 0;

2. 0x01 graphic
.

Wtedy transformata Laplace'a L{f(t)} istnieje dla s > C.

Uwaga. Funkcję f spełniającą warunki 1. i 2. powyższego faktu będziemy nazywali oryginałem.

Fakt 4.1.4 (o liniowości przekształcenia Laplace'a)

Jeżeli istnieją transformaty Laplace'a funkcji f i g oraz c R, to

1. istnieje transformata Laplace'a funkcji f + g oraz

L{f(t)+g(t)} = L{f(t)} + L{g(t)};

2. istnieje transformata Laplace'a funkcji cf oraz

L{cf(t)} = cL{f(t)}.

Fakt 4.1.5 (o jednoznaczności transformaty Laplace'a)

Jeżeli funkcje ciągłe 0x01 graphic
mają takie same transformaty Laplace'a: F(s) = G(s), to są równe na przedziale [0,∞).

4.2 METODA OPERATOROWA ROZWIĄZYWANIA RÓWNAŃ RÓŻNICZKOWYCH

Fakt 4.2.1 (transformata n-tej pochodnej)

Jeżeli funkcja f(t) oraz jej pochodne f'(t), f''(t), ..., f(n-1)(t) są oryginałami, a ponadto funkcja ta ma na przedziale (0,) ciągłą n-tą pochodną, to istnieje transformata L{f(t)} oraz

L{f(n)(t)} =

= snL{f(t)} - sn-1f(0+) - sn-2f'(0+) + ... - sf(n-2)(0+) - f(n-1)(0+) =

= snF(s) - sn-1f(0+) - sn-2f'(0+) + ... - sf(n-2)(0+) - f(n-1)(0+),

gdzie F(s) = L{f(t)}, 0x01 graphic
, 0x01 graphic
, ..., 0x01 graphic
.

Uwaga. Jeżeli funkcje f(t), f'(t), ..., f(n-1)(t) są ciągłe prawostronnie w punkcie t0 = 0, to f(0+) = f(0), f'(0+) = f'(0), ..., f(n-1)(0+) = f(n-1)(0).

4.3 WŁASNOŚCI PRZEKSZTAŁCENIA LAPLACE'A

Fakt 4.3.1 (zmiana skali)

Jeżeli funkcja f(t) jest oryginałem, to dla dowolnej stałej α > 0

L{f(αt)} = 0x01 graphic
,

gdzie F(s) oznacza obraz funkcji f(t).

Fakt 4.3.2 (o różniczkowaniu obrazu)

Jeżeli funkcja f(t) jest oryginałem, to

L{tnf(t)} =(-1)nF(n)(s),

gdzie F(s) oznacza obraz funkcji f(t).

Fakt 4.3.3 (o przesunięciu argumentów obrazu)

Jeżeli funkcja f(t) jest oryginałem, to dla dowolnej stałej a R

L{eαtf(t)} = F(s - a),

gdzie F(s) oznacza obraz funkcji f(t).

Fakt 4.3.4 (o przesunięciu argumentów oryginału)

Jeżeli funkcja f(t) jest oryginałem, to dla dowolnej stałej τ > 0

L{1(t - τ)f(t - τ)} = e-sτF(s),

gdzie F(s) oznacza obraz funkcji f(t).

Fakt 4.3.4 (o całkowaniu oryginału)

Jeżeli funkcja f(t) jest oryginałem, to

L 0x01 graphic
,

gdzie F(s) oznacza obraz funkcji f(t).

Fakt 4.3.5 (transformaty ważniejszych funkcji c.d.)

Funkcja

Transformata

shαt

0x01 graphic

chαt

0x01 graphic

tneαt

0x01 graphic

eαtsinβt

0x01 graphic

eαtcosβt

0x01 graphic

4.4 SPLOT FUNKCJI

Def. 4.4.1 (splot funkcji)

Niech funkcje f(t) i g(t) będą całkowalne na każdym przedziale [0,T], gdzie T > 0. Splot funkcji f(t) i g(t) oznaczamy symbolem 0x01 graphic
i określamy wzorem

0x01 graphic
.

Fakt 4.4.2 (własności splotu funkcji)

Niech funkcje f(t), g(t), h(t) będą całkowalne na każdym przedziale [0,T], gdzie T > 0 i niech c R. Wtedy

1. 0x01 graphic
;

2. 0x01 graphic
;

3. 0x01 graphic
;

4. 0x01 graphic
.

Tw. 4.4.3 (wzór Borela)

Jeżeli funkcje f(t) i g(t) są oryginałami, to istniej transformata Laplace'a ich splotu oraz

L{0x01 graphic
} = L{f(t)} L{g(t)}.



Wyszukiwarka

Podobne podstrony:
Równania różniczkowe-ćwiczenia, budownictwo, III semestr, Analiza matematyczna 3, Matematyka, Matma2
egzamin analiza 2006, BUDOWNICTWO IL PW, SEMESTR I, Analiza Matematyczna I, Egzaminy
(3953) równania różniczkowe, Budownictwo 2, Budownictwo, Matematyka
Tabela6, Technologia INZ PWR, Semestr 2, Analiza Matematyczna 2.2, Tabele
am4 Szeregi liczbowe, Informatyka i Ekonometria SGGW, Semestr 1, Analiza Matematyczna, materialy od
POCHODNE I ICH ZASTOSOWANIA, ZiIP, Semestr I, Analiza matematyczna
ZALICZENIA I POPRAWY, ZiIP, Semestr I, Analiza matematyczna
Ekonomia matematyczna egz 30.01.2015, Ekonomia II stopień, UMK 2013-2015, III semestr, Ekonomia mate
zadania pochodne2 (dr R. Lizak), 2 Semestr, Analiza matematyczna i algebra liniowa, zad mat
27112009, Informatyka i Ekonometria SGGW, Semestr 1, Analiza Matematyczna, materialy od starszych ro
Egzamin ANA1 04092000, Informatyka i Ekonometria SGGW, Semestr 1, Analiza Matematyczna, materialy od
przykładowe pochodne, PWR, semestr I, analiza matematyczna
analiza (2), Informatyka i Ekonometria SGGW, Semestr 1, Analiza Matematyczna, materialy od starszych

więcej podobnych podstron