Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Все в одном

.pdf
Скачиваний:
39
Добавлен:
13.02.2015
Размер:
23.06 Mб
Скачать

(t)-решение этой З.К. Покажем, что можно подобрать c1..n в (6): (6) задает функцию, которая будет решением З.К. По теореме ! это будет (t). Подставим (6) в нач. условия З.К.:

c1 c2

... cn

0

 

 

 

 

c2 2

...

cn n

'0

 

c1 1

.Эта система лин-х неоднор-х ур-ий для определения

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

c n 11

c

2

n 12

... c

n 1n (n 1)

0

1

 

 

 

 

n

 

 

c1..n. Если она однозначно разрешима, то из нее однозначным обоазом наход-ся набор

 

1

1 ...

1

 

 

 

 

 

1

2 ...

n

 

( i j )

i j

 

c1..n

= /определитель Вандермонда/ =

0

. =>

 

...

... ... ...

 

1 i j n

 

 

 

 

 

 

 

 

 

 

 

n 11

n 12

... n 1n

 

 

 

 

все нах-ся однозначно.

 

 

 

 

 

Вещественное решение уравнений с вещественными коэффициентами.

Пусть L(p):=a0pn+a1pn-1+...+an, аi , a00 L(p)z=0. (1).

Лемма: Пусть λ1..nÎC, λi≠ λj,i≠j => функции zi (t) e it ,i=1..n, образуют ЛН систему на R.

Теорема: все корни характеристического полинома L() простые и среди них ровно p вещественных корней 1..р и q пар комплексно-сопряженных корней 1..q, 1..q. p+2q=n. s= s+i s, s 0. Все корни различные для того, чтобы

p

q

z(t) ak e k t (bs e st d s e

 

st ) (5), где a,b,d – произвольные константы, была вещ-ой

 

k 1

s 1

<=> ak ak .bs ds

1) Пусть z(t) вещ-е реш-е, т.е. z(t)

_ p

свойству em em . Тогда z(t) z(t)

k 1

 

 

 

p

 

q

 

z(t) 0 . Из (5): z(t) ak e k t (

 

 

 

 

 

 

 

bs e

 

st ds e st )

по

 

 

 

k 1

 

s 1

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ak ak )e k t ((d s

 

 

 

 

st

 

 

 

s )e st ) (*).

 

 

bs )e

 

(bs

d

По

s 1

 

 

 

 

 

 

ak

 

 

 

k

 

 

 

 

 

 

 

a

.

 

 

 

 

 

 

лемме e k t , e st , e st - ЛНЗ на , т.е. все коэф-ты в (*) =0. Итак, z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bs

ds bs

ds

ak

ak

, то очевидно z(t) z(t) .

 

 

 

 

 

 

 

 

 

2) Если

 

 

 

 

 

 

 

 

 

 

 

 

ds

 

 

 

 

 

 

bs

 

 

 

 

 

 

 

 

 

 

 

 

Теорема: все корни характеристического полинома L() простые и среди них ровно

p вещественных корней 1..р и q пар комплексно-сопряженных корней 1..q,

 

1..q.

 

 

 

 

 

 

 

 

p

q

 

 

 

 

 

 

 

 

 

 

p+2q=n. s= s+i s, s 0. z(t) Ak e k t [Bs cos( s t) Ds sin( s t)] (7) задаёт

 

 

 

 

 

 

k 1

s 1

 

 

 

 

 

 

 

 

 

 

вещественное решение уравнения (1) при Ak,Bs,Ds . Всякое вещественное решение

(1) записывается ! образом в виде (7).

Пример 1.

x a 2 x 0, a 0 . 2-a2=0. =a. x(t)=C1eat+C2e-at, Ci или x(t)=Aeat+Be-at, A,B .

2

Пример 2. Ур-ие гармонич.осциллятора

x a 2 x 0, a 0 . 2+a2=0. =ia. x(t)=C1eiat+C2e-iat, Ci - обобщенное комплексное

_

решение или x(t)=Ceiat+ C e-iat, C - обобщенное вещественное решение в

комплексной форме или x(t)=Acos(at)+Bsin(at), A,B - обобщенное вещественное решение в вещественной форме.

Пример 3. Осциллятор с трением

 

 

2

x

 

 

 

 

 

 

2

2

x

2hx

 

0 , h>0, >0-частота колебаний, h-трение. Уравнение +2h+ =0 имеет

 

 

 

 

 

 

 

 

 

 

 

 

 

2 корня

 

 

h

h2 2 .

 

 

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

Если трение сильное и h> то , 1,2 отрицательны, различны

 

 

x(t) A e 1t A e 2t 0 , A . x(t) e 1t ( A

A e( 2 1 )t

) .

 

 

 

 

1

 

2

t

i

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если трение слабое и h< то =-h iq где q 2

 

h2 >0. Тогда x(t) Ce 1t

C

e 2t

 

 

 

R

t

 

x(t) e ht (Ceiqt

C

e iqt ) e ht R cos(qt ) ,где C

ei . x(t) 0

(маятник соверш-т

2

 

 

 

 

 

затух.колеб.).

 

 

 

 

 

 

 

 

 

3

~ ~ ~
z (t) z1 (t) z2 (t)

3. Линейное неоднородное уравнение n-го порядка с постоянными коэффициентами.

L(p)z=f(t) (1) и f(t): -> . Ему отвечает однородное L(p)z=0 (2).

Теорема:

Общее решение (1) имеет вид: z(t) z

~

~

 

0 (t) z (t) где z (t) - некоторое частное решение (1)

а z0(t) – общее решение (2).

 

 

 

Если ak общее решение (2) можно найти. Нужно найти

~

z (t) .

Свойство частного решения (1).

Если f(t)=f1(t)+ f2(t), ~1,2 (t) – частные решения уравнений L(p)z=f1(t) и L(p)z=f2(t)

z

- частное решение (1).

Лемма:(формула смещения) L(p)( e tz(t))=e tL(p+ )z(t), где ÎC, z(t)-некоторая ф-ия от t.

k

Опр: f(t) – квазимногочлен если f (t) pm (t)e mt , где m , mт, m≠n, pm

m 1

многочлен с коэффициентами.

Свойства (квазимногочлена).

илиΠ (квазимногочленi) = квазимногочлен. Const*(квазимногочлен) = квазимногочлен. (квазимногочлен)` = квазимногочлен.

(квазимногочлен) = квазимногочлен. L(p) (квазимногочлен) = квазимногочлен.

Теорема:

Частные решения уравнения L(p)=f(t)e t (5), где f(t)многоч-н degf(t)=r, ÎC,

представимы в виде ~ k t (6), где deg g(t) = deg f(t) = r. k=0,if L( )0, else, k -

z t g(t)e

кратность корня .

Подставим выр-е для реш-я (6) в ур-е (5): L(p)(tkg(t)e t)=f(t)e t. Восп-ся формулой смещ-я: e t L(p+ )(tkg(t))=f(t)e t.

L(p+ )={разложим в ряд Тейлора в окр-ти тчк

 

 

 

 

 

 

 

 

 

 

 

(n)

( )

 

 

 

 

 

 

 

t=0}= L( )

L ( )

p

 

L ( )

p2

...

L

pn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

2!

 

 

 

n!

 

 

 

 

 

 

 

Если корень хар-го многоч-на кратности k, то L( )=(-)kM( ), M( )≠0.

L’( )=k(-)k-1M( )+(-)kM’( ) => L’( )=0. =>

L’( )=…=L(k-1)( )=0.

Пусть f(t)=b0tr+b1tr-1+…+br-1t+br,

bi-известны;

 

 

 

 

g(t)=c0tr+c1tr-1+…+cr-1t+cr.

 

 

 

 

 

 

 

 

 

 

 

 

L ( )

 

L ( )

 

2

 

 

L(n) ( )

 

n

r

 

r-1

r

 

r-1

 

L( )

 

p

 

 

 

p

 

...

 

 

p

( c0t

+c1t

+…+cr-1t+cr)= b0t

+b1t

 

+…+br-1t+br

1!

2!

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В левой и правой частях стоят многочлены.Поэтому они равны, если равны коэф-ты при одинаковых степенях.

1

(i 1)t

t

r

 

 

 

L(k ) ( )

pk c0

(k r)(k r 1)...(r 1) b0

c0

k!b

 

 

 

 

 

 

 

 

0

 

 

 

 

 

k!

L(k ) ( )(k r)(k r 1)...(r 1)

 

 

 

 

 

 

 

 

 

 

 

t r 1

 

L(k 1) ( )

pk 1c0

(k r)(k r 1)...r

L(k ) ( )

pk c1 (k r 1)...r b1

 

c1 ...

 

(k 1)!

k!

 

 

 

 

 

 

 

 

 

 

 

 

И т.д. …

Т.о. мы приходим к системе для опред-я ci с нижнетреугольной матр-ей:

 

 

 

 

A11c0 b0

 

 

 

 

A11

0

...

0

 

 

 

 

 

 

A21c0 A22c1

b1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

A21

A22

...

0

 

Т.о.все коэф-ты однозначно определ-ся.

 

 

 

 

.......... .......... ..........

 

A

 

 

 

 

 

 

 

 

 

 

 

...

...

...

...

 

 

A

c

 

A

c ... A c

 

br

 

 

Ar1

Ar 2

...

Arr

 

 

 

r1

 

0

r 2

1

rr

r

 

 

 

 

 

Метод комплексных амплитуд.

Он примен-ся для поиска частн. реш-й ур-я L(p)z=F(t),L – полином с коэффициентами.

Теорема:

 

 

 

n

n-1

 

~

~

~

L(p)z=F(t), L(p):=a0p +a1p

+...+an, аi , a0 0 имеет частое решение

z (t) x (t) iy(t) ,

где

~

~

вещ-ные функции

~

 

~

 

x (t) ,

y (t)

x (t) - частное решение L(p)x=Re(F(t)),

y (t) -

 

частное решение L(p)y=Im(F(t)).

L(

~

F (t) L( p)(z (t))

~

p)x Re(F(t)),L(

L(

~

p) y

~ ~ p)(x (t) iy(t))

Im(F(t)) .

в _ силу _ лин ти

оп ра _ L( p)

~ ~

L( p)x (t) iL( p) y(t) Re(F(t)) i Im(F (t))

Пусть L(p)y=f(t), где L(p) оп-р с вещ.коэф-ми, f(t)-зад-ая вещ.ф-ция.

Рассм-м ур-е L(p)z=F(t) (1), где F(t): f(t)=Re(F(t)) or Im(F(t)). Находим частное решение

(1)

~

~

~

~

z (t) . Тогда реш-е исход-го вещ.ур-я

y(t) Re(z (t))or Im(z (t)).

Пример. x x et sin(t) .

Так как sin(t)=Im(eit) etsin(t)=Im(e(1+i)t) решение x(t)=x0(t)+xч(t). x0(t) – решение однородного x x 0 , 2-=0 1=0, 2=1 x0(t)=C1+C2et. Найдем частное решение:

z z e . =1+i – не корень 2-=0, f(t)=1. zч=Ae(1+i)t. A(1+i)2e(1+i)t- A(1+i)e(1+i)t=e(1+i)t (сокращаем) A(1+i)2- A(1+i)=1, A(1+2i-1)- A(1+i)=1, Ai-A=1.

A

1

 

 

i 1

zч i 1 e(1 i)t .

 

 

 

 

i 1

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

xч=Im(zч)=

et

Im(1 i)eit

 

 

et

Im(1 i)(cos t i sin(t))

 

1

t

(cos(t) sin(t)) . Ответ:

 

 

 

 

 

 

 

 

e

 

2

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) C

C

et

1

et (cos(t) sin(t)) .

 

 

 

 

 

 

 

 

 

1

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4. Нормальные линейные системы ДУ с переменными коэффициентами...

...имеют вид

a (t)

11

A(t) ...

an1 (t)

 

 

 

 

 

x1

 

 

 

 

F1 (t)

n

 

 

,

 

 

x

A(t)x

F (t) (1)

в , где

x

...

F (t)

... .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

Fn (t)

... a1n (t)

...

... . aij(t), Fj(t) – известные функции, определены и непрерывны на

...

 

ann (t)

t (r1,r2) (r1,r2 могут быть бесконечными), x – неизвестен. Система неоднородная если

 

 

t . Координатная форма

xi

n

(t)xi

fi (t) , i=1..n, (2).

 

 

 

 

F (t) 0

aij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

Свойства однородной системы

 

 

(3)

 

 

 

 

 

 

 

 

x

A(t)x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t) - некоторое решение (3), определенное на t (r1,r2). Если такая точка

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

t0 (r1,r2), что (t0 )

(t) 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

От противного: условия выполнены, но (t) 0 . Но оно – удовл-т З.К.:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x A(t)x

имеет решение x(t)=0, а в силу ед-ти з.к. др. реш-й нет.

 

 

 

 

 

 

 

 

x(t0 )

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

1..k (t) - набор из k решений (3) при Ci (t)

Ck k (t) - решение (3).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t) - некоторое решение

▲ (=>из св-в лин-ти системы). По предположению теоремы i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)=> i (t) A(t) i . Умножим обе части на ci и просуммируем:

 

 

 

 

 

k

 

 

k

 

k

 

 

 

 

- решение.

 

 

 

 

(t) Ci i

Ci A(t) i A(t) Ci i A (t)

(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

i 1

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1..k (t) - набор из k решений (3), определенный на (r1,r2) и такая точка t0 (r1,r2),

 

 

(t0 ) - линейно зависимый набор

 

 

 

 

 

 

 

 

 

что 1..k

1..k

(t)

линейно зависимый набор,

 

 

определенный на (r1,r2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ненулевой набор С1..k такой, что Ci i (t0 ) 0 ,

причем C 2 i 0 . Т.к. 1..k (t) -реш-

 

 

 

 

 

 

 

 

 

i 1

 

 

 

i 1

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

k

 

 

 

0 t (r1,r2),

 

 

я (3), то (t) Ci i (t) тоже реш-е.

Но (t0 )

Ci i (t0 )

0 (t)

 

 

 

 

 

 

i 1

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

что и означает ЛЗ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Опр: Фундаментальной системой решений x

A(t)x (3) называется набор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{ 1 (t),...,

n (t)} решений системы (3) ЛНЗ на (r1,r2). Число решений в наборе =

 

 

размерность пространства.

 

 

 

 

 

 

 

 

 

 

 

 

Теорема: (об общем решении).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t)} - ФСР (3)

общее решение (3)

 

 

 

 

n

 

 

 

 

 

 

 

 

 

{ 1

(t),..., n

записывается в виде x(t) Ci i

(t) (4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) фиксированного набора {C

}n

x(t) – решение системы.

 

 

 

 

 

 

 

 

 

 

i

 

i 1

 

 

 

 

 

 

 

 

 

 

2) всякое решение (3) представимо в (4) ! образом

1

1){ 1 (t),..., n (t)} - реш-я (3) => (по 2му св-ву) x(t) -тоже реш-е (3).

2)По То и!ЗК.

 

 

 

 

 

x(t) - некоторое решение (3) на (r1,r2), удовлетв-ее нач.усл-ю

 

 

 

n

 

 

 

 

Любой в-р можно разл-ть по выбранному базису:

x0

(t) Ci ai .

i 1

x0

x(t

0 ) при t0(r1,r2).

 

 

 

(*)

 

 

 

 

 

 

 

 

 

 

 

n

 

 

Рассм-м (4). В кач-ве const возьмем {Ci }in 1

 

 

 

из (*).

(t) Ci i

(t) -реш-е (3).

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

Рассмотрим это реш-е в тчк t0: (t0 ) Ci i (t0 ) .

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

Рассмотрим набор векторов

 

 

 

 

n

 

 

{ 1

(t0 ),..., n (t0 )} - ЛНЗ в . Если в кач-ве этого набора

 

 

 

 

 

 

 

 

 

 

 

 

взять {a1

,..., an } , то получим, что реш-е (t) и x(t)

удовл-т одним и тем же нач.усл-ям

 

 

 

 

 

 

 

 

 

 

 

 

x

A(t)x

 

 

 

 

x(t0 ) x0

(t0 ) . (т.е. З.К.

 

 

)

 

 

 

 

 

 

 

 

 

 

x(t0 ) x0

 

 

 

 

В силу ! З.К. они совп-т.

 

 

 

 

 

 

Теорема: (Критерий фундаментальности набора решений)

 

 

 

 

 

 

 

 

 

 

 

 

 

Набор решений (3) { 1 (t),..., n (t)} является ФСР вронскиан

 

 

 

 

 

 

11 (t)

...

n1 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отличен от нуля хотя бы в 1 точке t0 (r1,r2).

W (t) W ( 1

(t),..., n (t))

 

...

 

...

...

 

 

 

 

 

1n (t)

...

nn (t)

 

 

 

 

 

(=> тогда можно показать, что W≠0 t (r1,r2))

 

 

 

 

 

 

(t

 

 

 

 

(t

)}n

 

-ЛН => (по 3му св-ву реш-ий

W ( (t

),...,

n

)) 0 => система вект-в {

 

 

 

1 0

 

 

0

 

 

 

i

0

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

лин.систем) { i (t)}in 1 t (r1,r2)

 

 

 

 

 

 

 

 

 

 

 

 

(t)}n

 

 

 

 

 

 

)}n

-ЛН

 

(t)) 0

t0 (r1,r2).

{

i

 

- ФСР t (r1,r2). {

(t

W(t0) W (

(t),...,

 

i 1

 

 

 

i

0

i 1

 

 

 

 

1

n

 

 

Формула Лиувилля.

Лемма (о дифференцировании определителя)

uij – известные непрерывно-дифференцируемые функции на (r1,r2) и

u(t)

u11 (t)

... u1n (t)

ф-я u(t) диф-ма: u(t) u1 (t) ... un (t) , где ui(t) – определитель,

...

... ...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

un1 (t) ... unn (t)

 

 

 

 

 

 

 

 

 

 

 

полученный из u(t) дифференцированием i-й строки.

 

 

 

 

 

u(t) sgn u1 (1)

u1 (1) (π-перестановка, sgn π-четность перестановки).

 

 

Sn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

нечетная;

 

 

 

 

 

Sn – группа всех перестановок. sgn

четная.

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

 

 

 

 

 

 

 

 

 

 

 

 

 

u1 (t) ... un (t)

 

 

 

 

 

u(t) sgn u1 (1) u1 (1) ... sgn u1 (1)

u1 (1)

 

 

 

 

Sn

 

 

Sn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

ФСР (1) W ( (t),...,

 

 

 

sp A( )d

 

Теорема: набор функций {

i

(t)}n

-

n

(t)) W (t

0

)et0

(2)

 

 

 

 

 

i 1

 

 

1

 

 

 

(spA - след м-цы А).

2

 

 

 

 

k

(t) ( 1k

(t),..., nk

(t))

11 ...

... ...

Поэтому W (t) i1 ...

... ...

n1 ...

 

 

 

n

и k

A(t) k

записывается в координатах как ik

ais (t) sk (3) .

 

 

 

 

 

 

 

 

s 1

 

 

 

 

1n

 

 

 

 

...

 

по лемме W (t) W1 (t) ... Wn (t) , где

 

in

 

 

 

 

 

 

 

...

 

 

 

 

nn

 

 

 

 

 

11

...

1n

 

 

11

 

...

 

1n

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... ... ...

 

n

 

...

 

...

n

 

...

 

 

...

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wi (t)

 

...

 

 

 

a

is

(t)

s1

...

 

a

is

(t)

sn

ais (t)

i1

i1

in

 

 

 

 

 

 

 

 

 

 

s 1

 

 

 

 

 

s 1

 

 

 

 

s 1

 

 

... ... ...

 

 

 

...

 

...

 

 

...

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

...

nn

 

 

n1

 

...

 

nn

 

 

n1

 

 

 

 

 

 

 

 

...

1n

 

 

 

 

...

...

i s,det 0

...

in

 

aiiW (t) , поэтому

... ...

 

 

...

nn

 

 

n

W (t) W1 (t) ... Wn (t) aiiW (t) sp A(t)W (t) - уравнение с разделяющимися

i 1

t

 

sp A( )d

, где С - , в частности и С=W(t0).

переменными. W (t) Cet0

3

1