ch11 12 rr uklady

background image

UkÃlady r´

owna´

n r´

o˙zniczkowych liniowych pierwszego rze

,

du

Poprawiony 26 maja, 22:07

Pokazywali´smy jak mo˙zna rozwia

,

za´c r´ownanie x

00

+ x = 2 cos t + 12t sin t uzmien-

niaja

,

c obie sta le, kt´ore pojawiaja

,

sie

,

w rozwia

,

zaniu og´olnym r´ownania jednorodnego

x

00

+ x = 0 . W trakcie rozwia

,

zywanie ta

,

metoda

,

troche

,

niespodziewanie pojawi lo sie

,

r´ownanie c

0

1

(t) cos t + c

0

2

(t) sin t = 0 , po prostu zosta lo dopisane. Mo˙zna jednak wyja´sni´c

ska

,

d ono sie

,

bierze. Trzeba jednak w tym celu rozwa˙zy´c nieco og´olniejszy problem.

Mamy rozwia

,

za´c r´ownanie x

00

+x = 2 cos t+12t sin t . Wprowad´zmy oznaczenie x

0

= y .

Wtedy y

0

= x

00

= −x + 2 cos t + 12t sin t . Mo˙zemy wie

,

c zamiast r´ownania r´o˙zniczkowego

drugiego rze

,

du rozwia

,

zywa´c uk lad dw´och r´owna´

n pierwszego rze

,

du:



x

0

= y ,

y

0

= −x + 2 cos t + 12t sin t .

Niech A =



0

1

1 0



, ~g(t) =

0

2 cos t+12t sin t



oraz ~x =

x
y



. A jest wie

,

c macierza

,

, a ~x

niewiadoma

,

funkcja

,

zmiennej t , kt´orej warto´sciami sa

,

wektory, r´ownie˙z ~g jest funkcja

,

o

warto´sciach wektorowych. Uk lad r´owna´

n mo˙zemy zapisa´c tak:

~x

0

(t) = A~x(t) + ~g(t) .

(ur)

To bardzo przypomina r´ownanie liniowe

x

0

(t) = λx(t) + b(t) .

Rozwia

,

zaniem r´ownania liniowego jednorodnego x

0

= λx jest funkcja e

λt

· c , gdzie c ∈ R

lub c ∈ C , je´sli interesuja

,

nas rozwia

,

zania zespolone. Wyka˙zemy, ˙ze teraz jest tak samo,

no mo˙ze nieomal tak samo. W przypadku uk ladu funkcja powinna mie´c warto´sci wekto-

rowe. Chodzi te˙z o to ska

,

d wzia

,

´c liczbe

,

λ . To bardzo proste: nale˙zy znale´z´c warto´sci

w lasne i odpowiadaja

,

ce im wektory w lasne. Niech A~v = λ~v i ~x(t) = ~ve

λt

. Wtedy

~x

0

(t) = λ~ve

λt

= A~ve

λt

, a to oznacza, ˙ze funkcja ~x jest rozwia

,

zaniem jednorodnego

uk ladu r´owna´

n

~x

0

(t) = A~x(t) .

(jur)

Je´sli mamy troche

,

szcze

,

´scia i uda lo nam sie

,

znale´z´c wektory ~v

1

, ~v

2

, . . . , ~v

k

odpowia-

daja

,

ce warto´sciom w lasnym λ

1

, λ

2

, . . . , λ

k

, przy czym wektory ~v

1

, ~v

2

, . . . , ~v

k

sa

,

liniowo niezale˙zne,* to rozwia

,

zaniem og´olnym uk ladu jednorodnego jest funkcja

*

Oznacza to, ˙ze je´sli ~

v

1

c

1

+~

v

2

c

2

+···+~

v

k

c

k

=~

0 , to c

1

=c

2

=...=c

k

=0 ; je´sli k=1 , to ~

v

1

6=~

0 ; je´sli

k=2 , to wektory ~

v

1

,~

v

2

nie le˙za, na jednej prostej; je´sli k=3 , to wektory ~v

1

,~

v

2

,~

v

3

nie le˙za, na jednej

p laszczy´

znie.

1

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

~x(t) = ~v

1

e

λ

1

t

· c

1

+ ~v

2

e

λ

2

t

· c

2

+ · · · + ~v

k

e

λ

1

t

· c

k

.

Niech X(t) be

,

dzie macierza

,

, kt´orej kolumnami sa

,

kolejno funkcje o warto´sciach wektoro-

wych: ~v

1

e

λ

1

t

, ~v

2

e

λ

2

t

, . . . , ~v

k

e

λ

k

t

. Wtedy zachodzi r´owno´s´c ~x(t) = X(t) · ~c , kolejnymi

wsp´o lrze

,

dnymi wektora c sa

,

liczby c

1

, c

2

. . . , c

k

. Zauwa˙zmy przy okazji, ˙ze zachodzi te˙z

r´owno´s´c

X

0

(t) = AX(t) .

Wynika ona sta

,

d, ˙ze analogiczna r´owno´s´c jest spe lniona dla ka˙zdej kolumny macierzy X

i ze sposobu mno˙zenia macierzy. Macierz X zwana jest rozwia

,

zaniem fundamentalnym

uk ladu jednorodnego, o ile jej kolumny sa

,

wektorami liniowo niezale˙znymi, a to o niej

za lo˙zyli´smy.

Rozwia

,

zanie r´ownania liniowego niejednorodnego mogli´smy znale´z´c w postaci e

λt

·c(t) .

Rozwia

,

zania uk ladu r´owna´

n niejednorodnych mo˙zemy poszuka´c w analogiczny spos´ob, tj.

w postaci ~x(t) = X(t) · ~c(t) . Mamy wtedy

~x

0

(t) = X

0

(t)~c(t) + X(t)~c

0

(t) = AX(t)~c(t) + X(t)~c

0

(t) = A~x(t) + X(t)~c

0

(t) .

Wystarczy wie

,

c, by spe lniona by la r´owno´s´c

Ax(t) + ~g(t) = x

0

(t) = A~x(t) + X(t)~c

0

(t) ,

czyli ~g(t) = X(t)~c

0

(t) , by funkcja ~

X by la rozwia

,

zaniem uk ladu niejednorodnego. Ot´o˙z

mo˙zna udowodni´c, ˙ze det X(t)



6= 0 , a z tego wynika, ˙ze istnieje

X(t)



1

. Mo˙zemy

wie

,

c napisa´c ~c

0

(t) = X(t)



1

~g(t) . Sprowadzili´smy wie

,

c rozwia

,

zanie uk ladu r´owna´

n do

sca lkowania pewnej funkcji, kt´orej warto´sciami sa

,

wektory. Prze´sledzimy te

,

procedure

,

na

przyk ladzie, od kt´orego rozpocze

,

li´smy omawianie uk lad´ow r´owna´

n r´o˙zniczkowych.

Mamy wie

,

c k = 2 , A =



0

1

1 0



, ~g(t) =



0

2 cos t + 12t sin t



.

Uk lad jednorodny wygla

,

da wie

,

c tak:



x

0

1

(t)

x

0

2

(t)



=



0

1

1 0



·



x

1

(t)

x

2

(t)



.

R´ownanie charakterystyczne macierzy A wygla

,

da tak 0 =

−λ

1

1 −λ

= λ

2

+ 1 , zatem

jej warto´sciami w lasnymi sa

,

liczby ±i . Znajdziemy wektor w lasny

v

1

v

2



odpowiadaja

,

cy

warto´sci w lasnej i . Musi by´c spe lniony uk lad r´owna´

n (algebraicznych):



0 · v

1

+ 1 · v

2

= iv

1

,

1 · v

1

+ 0 · v

2

= iv

2

.

Oba r´ownania sa

,

r´ownowa˙zne, zatem musi by´c spe lniona r´owno´s´c v

2

= iv

1

. Niech v

1

= 1 .

Widzimy, ˙ze

1

i



jest wektorem w lasnym odpowiadaja

,

cym warto´sci w lasnej i . Poniewa˙z

macierz A jest rzeczywista, wie

,

c warto´sci w lasnej −i odpowiada wektor w lasny

1

−i



.

2

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

Wobec tego funkcje

1

i



e

it

oraz

1

−i



e

−it

sa

,

rozwia

,

zaniami wektorowymi jednorodnego

uk ladu r´owna´

n (r´o˙zniczkowych). Rozwia

,

zanie og´olne wygla

,

da wie

,

c tak:

1

i



e

it

c

1

+

1

−i



e

−it

c

2

=



e

it

e

−it

ie

it

−ie

−it



c

1

c

2



.

Mo˙zemy wie

,

c przyja

,

´c X(t) =



e

it

e

−it

ie

it

−ie

−it



. Szukamy rozwia

,

zania uk ladu niejednorod-

nego w postaci X(t)·~c(t) =



e

it

e

−it

ie

it

−ie

−it



·



c

1

(t)

c

2

(t)



. Podstawiaja

,

c do wyj´sciowego uk ladu

i redukuja

,

c wszystko, co tylko mo˙zna, otrzymujemy:



e

it

e

−it

ie

it

−ie

−it



·



c

0

1

(t)

c

0

2

(t)



=



0

2 cos t + 12t sin t



,

co mo˙zna zapisa´c te˙z nie u˙zywaja

,

c macierzy



c

0

1

(t)e

it

+ c

0

2

(t)e

−it

= 0,

c

0

1

(t)ie

it

− c

0

2

(t)ie

−it

= 2 cos t + 12t sin t.

Otrzymali´smy wie

,

c uk lad r´owna´

n z niewiadomymi funkcjami c

0

1

, c

0

2

taki sam, jak w me-

todzie uzmienniania sta lych stosowanej do r´ownania r´o˙zniczkowego drugiego rze

,

du. Tym

razem nie trzeba by lo wprowadza´c sztucznie pierwszego r´ownania, pojawi lo sie

,

samo. Teraz

nale˙zy rozwia

,

za´c otrzymany uk lad r´owna´

n. Mno˙za

,

c pierwsze przez e

it

, drugie przez −ie

it

i dodaja

,

c r´ownania stronami otrzymujemy

2c

0

1

(t)e

2it

= −i(2 cos t + 12t sin t)e

it

.

Sta

,

d wynika, ˙ze

c

0

1

(t) = −i(cos t + 6t sin t)e

−it

= −i(cos t + 6t sin t)(cos t − i sin t) =

= sin t cos t − 6t sin

2

t − i(cos

2

t + 6t sin t cos t) .

Z pierwszego r´ownania wynika, ˙ze

c

0

2

(t) = c

0

1

(t)e

2it

= i(cos t + 6t sin t)e

it

= i(cos t + 6t sin t)(cos t + i sin t) =

= sin t cos t − 6t sin

2

t + i(cos

2

t + 6t sin t cos t) .

Pozosta lo sca lkowa´c otrzymane pochodne poszukiwanych funkcji:

R

c

0

1

(t)dt =

R

sin t cos t − 6t sin

2

t − i(cos

2

t + 6t sin t cos t)



dt =

= 2 cos

2

t + 3 sin t cos t −

3
2

t

2

− i 2 sin t cos t + 2t − 3t cos

2

t



+ d

1

,

gdzie d

1

jest pewna

,

liczba

,

, by´c mo˙ze nierzeczywista

,

. Analogicznie

R

c

0

2

(t)dt =

R

sin t cos t − 6t sin

2

t + i(cos

2

t + 6t sin t cos t)



dt =

= 2 cos

2

t + 3 sin t cos t −

3
2

t

2

+ i 2 sin t cos t + 2t − 3t cos

2

t



+ d

2

,

gdzie d

2

jest pewna

,

liczba

,

by´c mo˙ze nierzeczywista

,

. Wobec tego rozwia

,

zania uk ladu nie-

jednorodnego wygla

,

daja

,

w tym przypadku tak:

x

1

(t) = 2 cos

2

t + 3t sin t cos t −

3
2

t

2

− i(2 sin t cos t + 2t − 3t cos

2

t) + d

1



e

it

+

+ 2 cos

2

t + 3t sin t cos t −

3
2

t

2

+ i(2 sin t cos t + 2t − 3t cos

2

t) + d

2



e

−it

=

3

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

= 2 cos

2

t + 3t sin t cos t −

3
2

t

2

− i(2 sin t cos t + 2t − 3t cos

2

t) + d

1



cos t + i sin t



+

+ 2 cos

2

t + 3t sin t cos t −

3
2

t

2

+ i(2 sin t cos t + 2t − 3t cos

2

t) + d

2



cos t − i sin t



=

= 4 cos t − 3t

2

cos t + 2t sin t + (d

1

+ d

2

) cos t + i(d

1

− d

2

) sin t .

Analogicznie

x

2

(t) = i 2 cos

2

t + 3t sin t cos t −

3
2

t

2

− i(2 sin t cos t + 2t − 3t cos

2

t) + d

1



e

it

− i 2 cos

2

t + 3t sin t cos t −

3
2

t

2

+ i(2 sin t cos t + 2t − 3t cos

2

t) + d

2



e

−it

=

= i 2 cos

2

t + 3t sin t cos t −

3
2

t

2

− i(2 sin t cos t + 2t − 3t cos

2

t) + d

1



cos t + i sin t



− i 2 cos

2

t + 3t sin t cos t −

3
2

t

2

+ i(2 sin t cos t + 2t − 3t cos

2

t) + d

2



cos t − i sin t



=

= 3t

2

sin t − 2t cos t − (d

1

+ d

2

) sin t + i(d

1

− d

2

) cos t .

Przyjmimy d

1

= α + βi , α, β ∈ R . Je´sli chcemy, by funkcja x

1

przyjmowa la jedynie

rzeczywiste warto´sci, nale˙zy za lo˙zy´c, ˙ze d

1

+ d

2

jest liczba

,

rzeczywista

,

, a d

1

− d

2

— liczba

,

czysto urojona

,

. Oznacza to, ˙ze Imd

2

= Imd

1

= −β i Red

2

= Red

1

= α . Rozwia

,

zania

przyjmuja

,

posta´c

x

1

(t) = 4 cos t − 3t

2

cos t + 2t sin t + 2α cos t − 2β sin t ,

x

2

(t) = 3t

2

sin t − 2t cos t − 2α sin t − 2β cos t .

Oczywi´scie ˙zadnych ogranicze´

n na liczby rzeczywiste α, β nie ma.

Doda´c trzeba, ˙ze mogli´smy rozwia

,

zanie og´olne uk ladu jednorodnego zapisa´c w inny

spos´ob. Np. zamiast rozwia

,

za´

n wektorowych

1

i



e

it

i

1

−i



e

−it

mogli´smy rozwa˙za´c funkcje

Re



1

i



e

it



=

1
2



1

i



e

it

+

1

−i



e

−it

] =

cos t

sin t



oraz

Im



1

i



e

it



=

1

2i



1

i



e

it

1

−i



e

−it

] =

sin t

cos t



.

Prowadzi do wzoru

X(t) =



cos t

sin t

sin t cos t



.

Zache

,

cam student´ow do przeprowadzenia oblicze´

n startuja

,

c z z w la´snie wskazanego roz-

wia

,

zania fundamentalnego. Obliczenia wcale nie be

,

da

,

trudniejsze ni˙z pokazane przeze

mnie.

Zajmiemy sie

,

teraz dwoma uk ladami trzech r´owna´

n liniowych z trzema niewiadomymi

funkcjami.

Przyk lad 15.1

Niech A =

5 2

4

2

0

2

2

1 1

. Rozwa˙zymy uk lad r´owna´n r´o˙zniczkowych

jednorodnych (zapisany w postaci wektorowej)

x

0

1

(t)

x

0

2

(t)

x

0

3

(t)

 = ~x

0

(t) = A~x(t) =

5 2

4

2

0

2

2

1 1

x

1

(t)

x

2

(t)

x

3

(t)

 =

5x

1

(t) 2x

2

(t) + 4x

3

(t)

2x

1

(t) + 2x

3

(t)

2x

1

(t) + x

2

(t) − x

3

(t)

.

Zaczniemy od znalezienia warto´sci w lasnych, potem zajmiemy sie

,

wektorami w lasnymi tej

4

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

macierzy. R´ownanie charakterystyczne wygla

,

da tak

0 =

5 − λ

2

4

2

0 − λ

2

2

1

1 − λ

= −λ

3

+ 4λ

2

5λ + 2 = (1 − λ)

2

(2 − λ) .

Wobec tego warto´sciami w lasnymi sa

,

liczby 1, 1, 2 .* Je´sli wektor ~v o wsp´o lrze

,

dnych

v

1

, v

2

, v

3

odpowiada warto´sci w lasnej 2 , to musza

,

by´c spe lnione r´owno´sci

( 5v

1

2v

2

+ 4v

3

= 2v

1

,

2v

1

+ 0v

2

+ 2v

3

= 2v

2

,

2v

1

+ 1v

2

1v

3

= 2v

3

,

czyli

( 3v

1

2v

2

+ 4v

3

= 0,

2v

1

2v

2

+ 2v

3

= 0,

2v

1

+ v

2

3v

3

= 0.

Odejmuja

,

c drugie r´ownanie od pierwszego otrzymujemy v

1

+2v

3

= 0 , a zatem v

1

= 2v

3

.

Po podstawieniu tego wyniku do dowolnego z trzech r´owna´

n otrzymujemy v

2

= −v

3

.

Przyjmuja

,

c np. v

3

= 1 (mo˙zemy te˙z v

3

= 1683 ·

1410

1525 ) stwierdzamy, ˙ze wektor

~v

1

=



2

1

1



odpowiada warto´sci w lasnej 2 . Teraz znajdziemy wektory odpowiadaja

,

ce

warto´sci w lasnej 1 . Je´sli v

1

, v

2

, v

3

sa

,

wsp´o lrze

,

dnymi takiego wektora, to

( 5v

1

2v

2

+ 4v

3

= v

1

,

2v

1

+ 0v

2

+ 2v

3

= v

2

,

2v

1

+ 1v

2

1v

3

= v

3

,

czyli

( 4v

1

2v

2

+ 4v

3

= 0,

2v

1

− v

2

+ 2v

3

= 0,

2v

1

+ v

2

2v

3

= 0.

Otrzymali´smy trzy r´ownowa˙zne r´ownania. Oznacza to, ˙ze mo˙zemy wybra´c dowolnie np.

warto´sci v

1

i v

3

i wyznaczy´c v

2

z wzoru v

2

= 2v

1

+ 2v

3

. Wybieraja

,

c np. v

1

= 1 . v

3

= 0

otrzymujemy v

2

= 2 . Podobnie z r´owno´sci v

1

= 0 i v

3

= 1 wynika r´owno´s´c v

2

= 2 .

Znale´zli´smy wie

,

c dwa wektory w lasne ~v

2

=



1

2
0



i ~v

3

=



0

2
1



, kt´ore nie le˙za

,

na jednej

prostej, bowiem z r´owno´sci ~0 = c

2

~v

2

+ c

3

~v

3

=



c

2

2c

2

+2c

3

c

3



wynikaja

,

oczywi´scie r´owno´sci

c

2

= 0 = c

3

. Mo˙zna bez trudu stwierdzi´c, ˙ze je´sli ~0 = c

1

~v

1

+ c

2

~v

2

+ c

3

~v

3

=



2c

1

+c

2

c

1

+2c

2

+2c

3

−c

1

+c

3



,

to c

1

= c

2

= c

3

= 0 . Wynika sta

,

d latwo, ˙ze dla ka˙zdego wektora ~v R

3

istnieje i to

dok ladnie jedna taka tr´ojka liczb c

1

, c

2

, c

3

, ˙ze ~v = c

1

~v

1

+ c

2

~v

2

+ c

3

~v

3

. Mo˙zemy wie

,

c

stwierdzi´c, ˙ze funkcja o warto´sciach wektorowych

~x(t) = c

1

~v

1

e

2t

+ c

2

~v

2

e

t

+ c

3

~v

3

e

t

=



2c

1

e

2t

+c

2

e

t

c

1

e2t+2c

2

e

t

+2c

3

e

t

−c

1

e

2t

+c

3

e

t



jest rozwia

,

zaniem og´olnym jednorodnego uk ladu r´owna´

n, tzn. ka˙zde rozwia

,

zanie mo˙zemy

zapisa´c w tej postaci dobieraja

,

c sta le c

1

, c

2

, c

3

w odpowiedni spos´ob.

Przyk lad 15.2

Niech A =

5 2 1

3 0

1

3 2

1

. Rozwa˙zymy wie

,

c uk lad r´owna´

n:

*

Liczba 1 jest powt´

orzona, bo jest pierwiastkiem podw´

ojnym r´

ownania charakterystycznego, czyli jest

podw´

ojna, warto´scia, w lasna,.

5

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

x

0

1

(t)

x

0

2

(t)

x

0

3

(t)

 = ~x

0

(t) = A~x(t) =

5

2 1

3

0

1

3

2

1

x

1

(t)

x

2

(t)

x

3

(t)

 =

5x

1

(t) + 2x

2

(t) − x

3

(t)

3x

1

(t) + x

3

(t)

3x

1

(t) + 2x

2

(t) + x

3

(t)

.

Zaczniemy, jak poprzednio, od znalezienia warto´sci w lasnych macierzy A . Szukamy takich

liczb λ , dla kt´orych spe lniona jest r´owno´s´c

0 =

5 − λ

2

1

3

0 − λ

1

3

2

1 − λ

= −λ

3

+ 6λ

2

12λ + 8 = (2 − λ)

3

.

Tym razem macierz ma jedna

,

, za to trzykrotna

,

, warto´s´c w lasna

,

. Je´sli ~v =



v

1

v

2

v

3



jest

wektorem w lasnym macierzy A , to musza

,

by´c spe lnione r´owno´sci:

( 5v

1

+ 2v

2

− v

3

= 2v

1

,

3v

1

+ 0v

2

+ v

3

= 2v

2

,

3v

1

+ 2v

2

+ v

3

= 2v

3

,

czyli

( 3v

1

+ 2v

2

− v

3

= 0,

3v

1

2v

2

+ v

3

= 0,

3v

1

+ 2v

2

− v

3

= 0.

Oznacza to, ˙ze ~v jest wektorem w lasnym macierzy A wtedy i tylko wtedy, gdy v

3

=

3v

1

+ 2v

2

. Niech ~v

1

=



1

0
3



, ~v

2

=



0

1
2



. Jest jasne, ˙ze je´sli ~v =



v

1

v

2

v

3



i v

3

= 3v

1

+ 2v

2

,

to ~v = v

1

· v

1

+ v

2

· v

2

. Wynika sta

,

d, ˙ze u˙zywaja

,

c rozwia

,

za´

n postaci ~ve

2t

nie znajdziemy

nigdy rozwia

,

zania ~x , dla kt´orego spe lniona by laby r´owno´s´c ~x(0) =



0

0
1



. Wobec tego

funkcja c

1

~v

1

e

2t

+ c

2

~v

2

e

2t

nie jest rozwia

,

zaniem og´olnym uk ladu.

W tej sytuacji nale˙zy znale´z´c tzw. uog´olniony wektor w lasny, tj. wektor ~

w , dla kt´orego

istnieje taki wektor w lasny ~v , ˙ze spe lniona jest r´owno´s´c

A~

w = λ~

w + ~v

bo

====

λ=2

2~

w + ~v .

Wtedy, jak to sprawdzimy za chwile

,

, funkcja (~

w + t~v)e

λt

bo

====

λ=2

(~

w + t~v)e

2t

. W ten spos´ob

znajdziemy brakuja

,

ce rozwia

,

zania.

Najpierw wyka˙ze

,

, ˙ze znalezienie wektora ~

w rozwia

,

zuje

,

problem:

(~

w + t~v)e

λt



0

= ~ve

λt

+ λ(~

w + t~v)e

λt

= (λ~

w + ~v)e

λt

+ tλ~ve

λt

= (A~

w)e

λt

+ t(A~v)e

λt

=

= A (~

w + t~v)e

λt



.

Teraz znajdziemy wektor ~

w . R´ownanie A~

w = λ~

w + ~v mo˙zna przepisa´c w postaci

(A − λI)~

w =

3

2 1

3 2

1

3

2 1

~w = ~v . Spr´obujmy przyja

,

´c ~

w =



0

0
1



. Wtedy

~v =

3

2 1

3 2

1

3

2 1

0
0
1

 =

1

1

1

~v jest wektorem w lasnym, bowiem

Przypominamy, ˙ze I to macierz jednostkowa, ma jedynki na g l´

ownej przeka,tnej, poza nia, — zera.

6

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

3

2 1

3 2

1

3

2 1

1

1

1

 =

0
0
0

.

Jest jasne , ˙ze je´sli ~

w =



w

1

w

2

w

3



, to

3

2 1

3 2

1

3

2 1

~w = (3w

1

2w

2

+ w

3

)

1

1

1

,

a wie

,

c zasta

,

pienie wektora ~

w innym mo˙ze jedynie w nieistotny spos´ob zmieni´c wektor

~v — zasta

,

pimy go r´ownoleg lym do niego, pomimo tego, ˙ze mamy do dyspozycji ca la

,

p laszczyzne

,

!

Mo˙zemy napisa´c w ko´

ncu

~x(t) = c

1

~v

1

e

2t

+ c

2

~v

2

e

2t

+ c

3

(~

w + t~v)e

2t

=

(c

1

− c

3

t)e

2t

(c

2

+ c

3

t)e

2t

(3c

1

+ 2c

2

+ c

3

− c

3

t)e

2t

.

Oczywi´scie dla ka˙zdego wektora ~u =



u

1

u

2

u

3



mo˙zna znale´z´c takie sta le c

1

, c

2

, c

3

, ˙ze ~u =

c

1

~v

1

+ c

2

~v

2

+ c

3

~

w = ~x(0) ,* a to oznacza, ˙ze funkcja ~x jest rozwia

,

zaniem og´olnym

uk ladu.

Obejrzymy jeszcze jeden uk lad jednorodny.

Przyk lad 15.3

Spr´obujemy rozwia

,

za´c uk lad r´owna´

n

x

0

1

(t) = 2x

1

(t) + x

2

(t),

x

0

2

(t) = 2x

2

(t) + x

3

(t),

x

0

3

(t) = 2x

3

(t).

Zapiszemy go w postaci wektorowej. Niech A =

2 1 0
0 2 1
0 0 2

. Bez trudu sprawdzamy,

˙ze liczba 2 jest potr´ojna

,

warto´scia

,

w lasna

,

macierzy A oraz ˙ze wektorami w lasnymi jej

odpowiadaja

,

cymi sa

,

wektory postaci ~v =



0

0

v

3



= v

3



0

0
1



. Zbi´or z nich z lo˙zony jest

wie

,

c jednowymiarowa

,

przestrzenia

,

liniowa

,

. Niech ~v

3

=



0

0
1



. Niech ~v

2

=



0

1
0



. Mamy

(A − 2I)~v

2

= ~v

3

. Mamy wie

,

c rozwia

,

zanie postaci c

2

(~v

2

+ t~v

3

)e

2t

+ c

3

~v

3

e

2t

. Nie jest

to niestety rozwia

,

zanie og´olne, bo niezale˙znie od tego, co uczynimy ze sta lymi c

2

, c

3

,

to podana funkcja w punkcie 0 nie przyjmie warto´sci



1

0
0



. Zauwa˙zmy jednak, ˙ze je´sli

~v

1

=



1

0
0



, to (A − 2I)~v

1

= ~v

2

, czyli A~v

1

= 2~v

1

+ ~v

2

. Niech

*

Wystarczy przyja,´c c

1

=u

1

, c

2

=u

2

i c

3

=u

3

2u

2

3u

1

.

7

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

~x

1

(t) = ~v

1

+ t~v

2

+

1
2

t

2

~v

3



e

2t

.

Mamy teraz

~x

0

1

(t) = ~v

2

+ t~v

3



e

2t

+ 2 ~v

1

+ t~v

2

+

1
2

t

2

~v

3



e

2t

=

= (2~v

1

+ ~v

2

)e

2t

+ t(2~v

2

+ ~v

3

)e

2t

+ t

2

~v

3

e

2t

= A~v

1

e

2t

+ tA~v

2

e

2t

+

1
2

A~v

3

e

2t

=

= A ~v

1

+ t~v

2

+

1
2

t

2

~v

3



e

2t

= A~x

1

(t) .

Wskazali´smy wie

,

c naste

,

pne rozwia

,

zanie. Niech

~x(t) = c

1

~v

1

+ t~v

2

+

1
2

t

2

~v

3



e

2t

+ c

2

~v

2

+ t~v

3



e

2t

+ c

3

~v

3

e

2t

=

c

1

e

2t

(c

1

t + c

2

)e

2t

(

1
2

c

1

t

2

+ c

2

t + c

3

)e

2t

.

Funkcja ta jest rozwia

,

zaniem og´olnym uk ladu, bo manipuluja

,

c wsp´o lczynnikami c

1

, c

2

, c

3

mo˙zna uzyska´c dowolna

,

warto´s´c ~x(0) .

Pojawi l sie

,

tu nowy spos´ob znajdowania. rozwia

,

zania og´olnego. Nie be

,

dziemy dowodzi´c

og´olnego twierdzenia, kt´ore kryje sie

,

za przedstawianymi rozumowaniami, bo nie mamy na

to czasu.

Definicja 15.1 (serii)

Cia

,

g wektor´ow ~v

1

, ~v

2

, . . . , ~v

m

nazywamy seria

,

Jordana d lugo´sci m odpowiadaja

,

ca

,

warto´sci w lasnej λ macierzy A wtedy i tylko wtedy, gdy

A~v

1

= λ~v

1

+ ~v

2

, A~v

2

= λ~v

2

+ ~v

3

, . . . , A~v

m−1

= λ~v

m−1

+ ~v

m

i A~v

m

= λ~v

m

.

W ostatnim przyk ladzie wysta

,

pi la jedna seria d lugo´sci 3 ( ~v

1

, ~v

2

, ~v

3

) i dzie

,

ki niej

uda lo sie

,

napisa´c rozwia

,

zanie uk ladu w postaci og´olnej. W poprzednim przyk ladzie wy-

ste

,

powa la seria d lugo´sci 2 . Wektor w lasny to seria d lugo´sci 1 .

Twierdzenie Jordana, kt´orego nie be

,

dziemy dowodzi´c, m´owi, ˙ze je´sli mamy

n –krotna

,

warto´s´c w lasna

,

, to mo˙zemy znale´z´c serie Jordana je odpowiadaja

,

ce w taki

spos´ob, by ich suma d lugo´sci by la r´owna krotno´sci warto´sci w lasnej n i by wszystkie

wektory wyste

,

puja

,

ce w tych seriach by ly liniowo niezale˙zne. Nie m´owimy te˙z jak mo˙zna

to zrobi´c. W przyk ladach, kt´ore rozpatrujemy, be

,

dzie to proste.

Maja

,

c dana

,

serie

,

d lugo´sci m mo˙zemy napisa´c odpowiadaja

,

ce jej rozwia

,

zania ~v

1

+

t~v

2

+

1
2

t

2

~v

3

+ · · · +

1

(m−1)!

t

m−1



e

λt

, ~v

2

+ t~v

3

+

1
2

t

2

~v

4

+ · · · +

1

(m−2)!

t

m−2



e

λt

,

. . . , ~v

m−1

+ t~v

m



e

λt

, ~v

m

e

λt

.

W rozwia

,

zaniu og´olnym pojawi sie

,

sk ladnik

c

1

~v

1

+ t~v

2

+

1
2

t

2

~v

3

+ · · · +

1

(m−1)!

t

m−1



e

λt

+ c

2

~v

2

+ t~v

3

+

1
2

t

2

~v

4

+ · · · +

1

(m−1)!

t

m−1



e

λt

+

+ · · · + c

m−1

(~v

m−1

+ t~v

m

)e

λt

+ c

m

~v

m

e

λt

.

Wynika z tego, ˙ze wsp´o lrze

,

dne rozwia

,

zania og´olnego sa

,

quasiwielomianami, kt´orych

stopie´

n jest o co najmniej jeden mniejszy od d lugo´sci ka˙zdej serii zwia

,

zanej z dana

,

war-

to´scia

,

w lasna

,

, kt´ora oczywi´scie pojawia sie

,

w wyk ladniku.

8

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

Mo˙zna wie

,

c nie szuka´c serii, lecz za lo˙zy´c ˙ze wsp´o lrze

,

dne sa

,

quasiwielomianami, kt´orych

stopie´

n jest o co najmniej jeden mniejszy od krotno´sci warto´sci w lasnej (to na og´o l za du˙zo)

podstawi´c taka

,

og´olna

,

funkcje

,

wektorowa

,

w miejsce poszukiwanej do uk ladu i znale´z´c

warunki na wsp´o lczynniki; oczywi´scie na wste

,

pie nie zak ladamy ˙zadnych relacji mie

,

dzy

wsp´o lczynnikami wyste

,

puja

,

cym na r´o˙znych wsp´o lrze

,

dnych.

Zadania

15. 01 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= x − 3y,

y

0

= 3x + y.

15. 02 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= x + y,

y

0

= 3y − 2x.

15. 03 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

+ x + 5y = 0,

y

0

− x − y = 0.

15. 04 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= 2y − 3x = 0,

y

0

= y − 2x.

15. 05 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= 4x − y,

y

0

= 5x + 2y.

15. 06 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= 2x + y,

y

0

= 4y − x.

15. 07 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

= 3x − y,

y

0

= 4x − y.

15. 08 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:



x

0

5x − 3y = 0,

y

0

+ 3x + y = 0.

15. 09 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= 2x + y + z

y

0

= 2x − z,

z

0

= 2x + y + 2z.

15. 10 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

(

x

0

= x − y − z,

y

0

= x + y,

z

0

= 3x + z.

15. 11 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= 4x − y − z,

y

0

= x + 2y − z,

z

0

= 2y + 3z − x.

15. 12 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= y − 2z − x,

y

0

= 4x + y,

z

0

= 2x + y − z.

15. 13 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= 2x − y − z,

y

0

= 2x − y − 2z,

z

0

= 2z − x + y.

15. 14 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= 2x − y + z,

y

0

= x + 2y − z,

z

0

= x − 2y + 2z.

9

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

15. 15 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= 2x + 2z − y,

y

0

= x + 2z,

z

0

= y − 2x − z.

15. 16 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

x

0

= x − y + z,

y

0

= x + y − z,

z

0

= 2z − y.

15. 17 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

(

x

0

= 2x + y,

y

0

= 2y + 4z,

z

0

= x − z.

15. 18 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

(

x

0

= 4x − y,

y

0

= 3x + y − z,

z

0

= x + z.

15. 19 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

w

0

= 7w − 2x + 5y − 10z,

x

0

= 5w − 2x − 5z,

y

0

= 2w + 3y − 4z,

z

0

= 5w − x + 5y − 8z.

15. 20 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n:

w

0

= 7w − 4x + 5y − 10z,

x

0

= 3w − 2x − 3z,

y

0

= 2w + 3y − 4z,

z

0

= 5w − 2x + 5y − 8z.

15. 21 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= y + tg

2

t − 1,

y

0

= −x + tg t ;



x

0

= 2y − x,

y

0

= 4y − 3x +

e

3t

1+e

2t

.

15. 22 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= 4x − 2y +

2

e

t

1

,

y

0

= 6x + 3y −

3

e

t

1

;



x

0

= x − y +

1

cos t

,

y

0

= 2x − y.

15. 23 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= 3x − 2y,

y

0

= 2x − y + 15e

t

t.

15. 24 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= y + tg

2

t − 1,

y

0

= −x + tg t.

15. 25 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= 2y − x,

y

0

= 4y − 3x +

e

3t

1+e

2t

.

15. 26 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= 4x − 2y +

2

e

t

1

,

y

0

= 6x + 3y −

3

e

t

1

.

15. 27 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= x − y +

1

cos t

,

y

0

= 2x − y.

10

background image

Uk lady r´owna´

n r´o˙zniczkowych liniowych pierwszego rze

,

du

15. 28 Rozwia

,

za´c uk lad r´owna´

n (mo˙zna ewentualnie uzmiennia´c sta le):



x

0

= 3x − 2y,

y

0

= 2x − y + 15e

t

t.

15. 29 Niech M =

0 1 1
1 0 1
0 1

0

, v =

2

1

1

.

Znale´z´c iloczyn M · v .

Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n x

0

(t) = M · x(t) .

Znale´z´c rozwia

,

zanie uk ladu r´owna´

n x

0

(t) = M · x(t) spe lniaja

,

ce warunek

x(0) =

2

1

1

.

Znale´z´c rozwia

,

zanie uk ladu r´owna´

n x

0

(t) = M · x(t) spe lniaja

,

ce warunek

x(0) =

1
1
2

.

15. 30 Znale´z´c rozwia

,

zanie og´olne uk ladu r´owna´

n r´o˙zniczkowych



x

0

(t) = −x(t) + y(t),

y

0

(t) = −x(t) 3y(t).

Znale´z´c rozwia

,

zanie szczeg´olne spe lniaja

,

ce warunek pocza

,

tkowy x(0) = 0 , y(0) = 7 .

15. 31 Rozwia

,

za´c uk lad r´owna´

n



x

0

(t) = 3x(t) − y(t) + e

4t

· sin t,

y

0

(t) = x(t) + 5y(t) − e

4t

· sin t.

(4)

Znale´z´c rozwia

,

zanie uk ladu (4) spe lniaja

,

ce warunki x(0) = 2 i y(0) = 2 .

Znale´z´c rozwia

,

zanie uk ladu (4) spe lniaja

,

ce warunki x(0) = 0 = y(0) .

15. 32 Znale´z´c warto´sci i wektory w lasne macierzy



1 2
2

5



.

Rozwia

,

za´c uk lad r´owna´

n



x

0

(t) = x(t) 2y(t) + 3e

3t

· cos(3t),

y

0

(t) = 2x(t) + 5y(t) 3e

3t

· cos(3t).

15. 33 Rozwia

,

za´c uk lad r´owna´

n



x

0

(t) = 6x(t) 3y(t),

y

0

(t) = 8x(t) 5y(t).

15. 34 Znale´z´c rozwia

,

zanie uk ladu r´owna´

n



x

0

(t) = 13x(t) + 25y(t),

y

0

(t) = 9x(t) + 17y(t),

kt´ore spe lnia warunek x(0) = 2 , y(0) = 1 .

11


Wyszukiwarka

Podobne podstrony:

więcej podobnych podstron