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

Лекции по дифференциальным уравнениям (Абрамов А.А

.).pdf
Скачиваний:
142
Добавлен:
12.03.2016
Размер:
3.56 Mб
Скачать

10. Матричные формулы решений обыкновенной системы линейных дифференциальных уравнений с п

Теорема: Пусть B произвольная комплексная n n-матрица. Тогда ряд E + B + B2 + : : :

1! 2!

для каждого матричного элемента сходится и притом абсолютно.

Доказательство:

B =

 

 

 

 

 

 

 

, B2

(2)

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:b11: : : ::::::: :b:1:n:

= b: 11: : : : :: :: :: : :b:1:n:

brs

j

6

 

 

 

 

 

 

b

n1

: : :

b

nn

 

(2)

 

 

 

(2)

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bn1

: : : bnn

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

brs

= brqbqs

 

 

brs

 

6 n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

qP

 

 

 

 

 

j

(2)

j

 

(3)

 

n

 

(2)

 

(3)

 

(m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценим матричные элементы степеней:

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

brs =

=1 brq bqs; jbrs j 6 n2 3, jbrs j 6 nm 1 m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qP

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rs + 1!

 

+ 2! + : : :; rs = (0; r = s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

brs

 

 

brs

 

 

 

1; r = s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

n 2

 

n2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Этот ряд мажорируется рядом: 1 +

 

+

 

 

 

+

 

 

 

 

 

+ : : : сходится по признаку Д'Аламбера.

 

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

3!

 

 

 

 

 

 

 

 

 

Определение: Пусть B произвольная комплексная n n матрица.

 

 

 

 

 

 

 

 

 

 

 

eB

def

 

 

 

B

 

B2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= E +

 

 

 

 

+

 

 

 

+ : : :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

 

 

 

 

Свойства:

1. BC = CB ) eB eC = eB+C

Доказательство: Ряд для eB è eC формально перемножить, собрать члены при одинаковых степенях. Так можно делать, потому что ряды сходятся абсолютно.

2.(eB) 1 = e B

Доказательство:

B ( B) = B B ) eB e B = eB B = e0 = E

Функция etA определяется следующим образом: etA = E +

tA

tA2

 

+

 

 

+ : : :.

1!

2!

 

 

 

1)et1A et2A = e(t1+t2)A, t1A è t2A перестановочны;

2)detA = AetA = etAA, степенной ряд почленно дифференцируемый. dt

2.Теорема: Рассмотрим уравнение (1). Пусть 6 t0 6 . Тогда функция

~y(t) = ~y0e(t t0)A + tZ0

t

 

e(t )Af~( )d

(3)

является решением (1), удовлетворяющим

~y(t0) = ~y0:

(4)

32

ГЛАВА 1. МЕТОДЫ РЕШЕНИЯ НЕКОТОРЫХ УРАВНЕНИЙ.

V

Каждое решение задачи (1) (4) совпадает на сво¼м отрезке определения с функцией, задаваемой формулой (3); решение задаваемое формулой (3) определено на вс¼м отрезке [ ; ].

Äëÿ (2):

z(t) = ~z0eA(t t0); ãäå ~z0 = ~z(t0)

Доказательство: y(t) = etA~u(t) ~u(t) = e tA~y(t)

 

n

d(etA~u)

= AetA~u + f~o

 

n

etA

 

d~u

= AetA~u + f~o

 

 

 

 

)

d

~u + etA

 

 

 

)

 

dt

dt

dt

 

 

 

 

 

d~u

 

 

 

 

 

d~u

 

 

 

d~u

 

)

nAetA~u + etA

 

= AetA~u + f~o

)

 

netA

 

 

= f~o

)

n

 

= e tAf~o

dt

 

dt

dt

 

 

~u(t) = tZ0

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e Af( )d + C;~

C~ произвольная постоянная

 

Отсюда ~y(t) = etAC~ + tZ0

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

etAe Af~( )d , òàê êàê tA ( A) = ( A) tA.

 

 

 

~y(t0) = ~y0 = et0AC;~ C~ = ~y0e t0A: Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~y(t) = ~y0e(t t0)A + tZ0

e(t )Af~( )d

 

 

 

 

3. Алгебраический способ вычисления etA.

Вспомогательные утверждения.

1.A = CBC 1 ) etA = CetBc 1

Доказательство:

 

etA = E + t

CBC 1

+ t2

CBC 1CBC 1

+ : : :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

B2

2!

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

= C(E + t

 

+ t2

 

 

+ : : :)C 1

= CetBC 1

 

 

 

2!

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. A = CLC 1; L = diag( 1; : : : ; n) Тогда etL = diag(e 1t; : : : ; e nt)

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

2 L2

 

 

 

 

 

 

 

 

 

tL

 

 

 

 

L

 

 

 

 

 

 

 

 

 

e

 

= E + t

 

 

+ t

 

 

+ : : : ;

 

ãäå

 

 

 

 

1!

2!

 

 

 

 

 

 

L2 = diag( 12; : : : ; n2 ) è ò.ä.

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

12

 

 

 

 

 

 

 

 

 

 

n

 

n2

 

 

etL = diag 1 + t

 

+ t2

 

+ : : : ; : : : ; 1 + t

 

 

+ t2

 

+ : : : = diag e 1t; : : : ; e nt

1!

2!

1!

2!

10. Матричные формулы решений обыкновенной системы линейных дифференциальных уравнений с п

4. Жорданова форма матрицы

Определение: Упорядоченный набор векторов ~g1; : : : ;~gk называется цепочкой, соответству- ющей собственному значению линейного преобразования , если:

8 A~g2

= ~g2 + g16

~

, где векторы ~g

; : : : ;~g

присоединенные к вектору ~g .

>

A~g1

= ~g1; ~g1 = 0

2

k

1

: : : : : : : : : : : : : : : : : :

>

 

 

 

 

 

 

<

>

>

: A~gk = ~gk + ~gk 1

Теорема: (Жордана) Пусть А линейное преобразование в комплексном линейном пространстве. Тогда существует базис этого пространства, составленный из цепочек ( Жорданов базис).

A =

 

J1 ...

0

, где каждая клетка Jr соответствует отдельным цепочкам.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

JN

 

 

 

 

 

 

r 1 0

Jr

def

 

... 1

 

Жорданова клетка размера kr kr.

=

 

 

 

 

 

 

 

 

0

 

 

r

 

 

 

 

 

 

 

J1

 

 

 

0

 

 

 

 

 

 

J1

 

 

0

 

 

e 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

etA =

 

tJ

..

 

 

, òàê êàê An =

 

 

 

..

 

 

 

 

 

 

=

 

 

n

..

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

.

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

e

tJN

 

 

 

 

 

 

 

0

 

 

 

 

JN

 

 

 

 

0

 

n

 

 

 

 

 

 

 

 

 

J1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

JN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 J12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2!

 

 

 

 

 

 

 

 

 

 

 

E

 

 

0

 

 

t 1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

etA =

 

 

..

.

 

 

+

 

 

...

0

 

+

 

 

 

 

 

 

..

.

 

0

 

 

+ : : :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

N

 

 

 

 

 

0

 

 

 

 

 

2 JN

 

 

 

 

 

 

 

 

 

E

 

 

 

 

t

1!

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Íàéä¼ì etJ : J = E + Q; ãäå Q =

 

 

...

 

1

; t E и tQ перестановочны, поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

et( E+Q) = et E etQ = et EetQ = et etQ:

При каждом умножении матрицы Q на себя, ¾диагональ¿, содержащая единицы, смещается на

единицу в сторону правого левого угла матрицы. Тогда Qk 1 =

0

: : : 0

0

, Q

нулевая

 

 

 

 

 

 

 

 

 

0

: : : 0

1

 

k

 

 

 

 

 

 

 

 

: : : : : : : : : : : :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

: : : 0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

матрица, также Qs = 0 ïðè s

 

k. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

etQ = E + t

Q

 

+ t2

Q2

+ : : : + tk 1

Qk 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

2!

(k 1)!

 

 

 

 

34

ГЛАВА 1. МЕТОДЫ РЕШЕНИЯ НЕКОТОРЫХ УРАВНЕНИЙ.

конечная сумма и

 

 

 

 

(ktk 1)!

 

 

 

1 t=1! : : :

 

 

 

 

 

1

 

 

 

etA = et

..

.

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

t=1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11. Элементы операционного исчисления

1. Определение преобразования Лапласа.

Определение: Комплексная функция f(t) называется оригиналом, если:

1.f определена при 0 6 t 6 1;

2.f кусочно-непрерывна (может иметь разрывы только первого рода)

3.

f

 

функция конечного порядка роста, то есть

9

вещественные числа

; M

:

 

 

Me

t

, äëÿ âñåõ 0 6 t 6

1

 

 

 

jf(t)j

 

 

 

 

 

Теорема: Пусть f(t) оригинал, jf(t)j 6 Me t, p комплексное число, такое что Re p > .

1

 

 

 

Тогда Z

f(t)eptdt существует.

0

 

 

 

Доказательство:

 

 

 

 

jf(t)e ptj = jf(t)jje ptj 6 Me te tRe p = Me (Re p )t;

 

 

1

 

Re p > 0 )

Z0

f(t)e (Re p )tdt сходится, а следовательно по признаку сравнения су-

 

 

1

ществует и интеграл Z0

f(t)eptdt.

Определение: Пусть f(t) оригинал, jf(t)j Me t. Функция f~(p), определ¼нная для комплексных p при Re p > 0 формулой

1

Z

f~(p) = f(t)e ptdt

0

называется изображением по Лапласу оригинала f(t). Определение: Переход от f(t) к f~(p) преобразование Лапласа.

11. Элементы операционного исчисления

35

Далее Re p достаточно велико.

f(p) соответствующее изображение

 

f(t) L! f~(p) ,

 

 

f(t) оригинал

 

2. Свойства преобразования Лапласа

Теорема: Пусть f(t) ! f~(p), комплексное число. Тогда

L

f(t) ! f~(p)

L

Теорема: Пусть f1(t) ! f~1(p), f2(t) ! f~2(p). Тогда

L L

f1(t) + f2(t) ! f~1(p) + f~2(p)

L

Доказательство: Пусть выполнены условия

jf1(t)j 6 M1e 1t; jf2(t)j 6 M2e 2t:

Тогда по неравенству треугольника

jf1(t) + f2(t)j 6 Me t; = max( 1; 2); M = M1 + M2

Теорема: Пусть f(t) оригинал,

t

Z

F (t) = f( ) d :

0

Тогда F (t) оригинал.

Доказательство: Пусть jf(t)j 6 Me t. Бер¼м > 0. Тогда

tt

jF (t)j = Z

f( ) d 6 Z

Me = (e t 1) 6

e t:

 

 

 

M

M

 

 

 

 

 

00

Теорема: Пусть f(t) ! f~(p). Пусть также f(t) непрерывно дифференцируемая,

L

f0(t) = g(t); g(t) ! g~(p):

L

Тогда

g~(p) = pf~(p) f(0):

Доказательство:

 

 

 

+1

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

Z

 

 

 

 

 

 

g~(p) =

 

e ptg(t) dt =

e ptf0(t) dt =

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

+1

 

+

 

 

e pt

f(

t)

0+1

Z

pe ptf(t) dt = f(0) + p Z0

 

1 e ptf(t) dt

 

 

 

 

 

 

 

 

=0

 

 

0

 

 

 

 

|

 

{z

 

 

 

 

 

 

 

 

}

 

 

 

 

 

 

Теорема: Пусть f(t) ! f~(p), и выполнены следующие условия:

L

36

ГЛАВА 1. МЕТОДЫ РЕШЕНИЯ НЕКОТОРЫХ УРАВНЕНИЙ.

1.f(t) m раз непрерывно дифференцируемая, f(m)(t) = g(t).

2.g(t) ! g~(p)

L

Тогда

g~(p) = pmf~(0) pm 1f(0) pm 2f0(0) : : : f(m 1)(0)

Доказательство:

Z t

f(m 1)(t) = f(m)( ) d + f(m 1)(0) оригинал

0

f(m 2); : : : ; f0(t) оригиналы.

Утверждение доказывается по индукции.

Теорема: Пусть f(t) ! f~(p), и выполнены следующие условия:

L

1.f(t) непрерывно дифференцируемая.

2.jf(t)j 6 Me t.

3.> .

Тогда при t > 0 имеет место

 

 

T

 

 

 

 

( ) = 2 T !+1 Z

( +i )t

~

 

 

1

 

 

 

f t

 

lim e

 

 

f( + i ) d

(3)

T

При t = 0 правая часто да¼т не f(0), а 12f(0).

Теорема: Пусть f1(t) ! f~1(p), f2(t) ! f~2(p). Пусть кроме того, f~1(p) = f~2(p).

L L

Тогда f1(t) = f2(t)

Доказательство:

 

 

 

 

 

 

 

 

F1(t) = Z0t f1( ) d ; F2(t) = Z0t f2( ) d

 

~

 

 

 

 

 

~

 

 

F1(t) ! F1(p); F2

(t) ! F2(p)

 

 

L

 

 

 

 

L

 

 

 

f1(t) = F10(t); f~1(p) = pF1(p) f1(0) = pF1(p)

~

f~1(p)

 

f~2(p)

~

 

)

f2(p) = pF2(p); F2(p) =

 

 

=

 

 

= F2

(p)

p

 

p

F1(t) = F2(t); f1(t) = f2(t)

3.Пусть выполнены следующие условия:

11. Элементы операционного исчисления

37

1.'(t) = tke t, комплексное, 0 6 k, k 2 Z.

2.> Re

Найд¼м преобразование Лапласа функции '.

'(t)e t ! 0 при t ! +1. Кроме того, j'(t)e tj 6 M, j'(t)j 6 Me t. Рассматривая p: Re p > , Re ( p + ), вычисляем

Z +1

e pttke t dt:

0

 

 

def

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k = Z0

 

tke zt dt; Re z > 0

 

 

 

 

 

 

 

 

 

 

 

+1

 

 

1

 

 

 

 

 

 

 

 

 

 

0(z) = Z0

e zt

=

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

6 k :

 

+1

 

zt

+1

 

 

+1 k 1e

zt

 

k

 

 

 

 

 

 

 

 

k(z) =

 

tke zt dt = tk

e

0

 

 

 

 

 

kt

 

 

dt =

 

k 1

(z);

0

z

0

 

 

z

 

 

z

 

Z

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

tke t !

 

 

k!

где 0 6 k целое, комплексное:

 

 

 

 

 

 

;

 

 

(p

 

)k+1

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

Теорема: Пусть f(t) ! f~(p), f(t) квазимногочлен.

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

q(p)

 

Тогда это равносильно тому, что f~(p) правильная рациональная дробь, то есть f~(p) =

,

r(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где q; r многочлены, deg q < deg r.

 

 

 

 

 

 

 

 

Доказательство: Очевидно, f~(p) = Ps

 

rs

 

 

(p s)ks+1

 

 

 

 

Рассмотрим функцию

 

 

 

 

 

 

 

 

ks s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

def

 

 

 

 

 

 

 

 

h(t) = Xs

st e

 

 

 

Тогда h(t) ! f~(p)

 

 

 

 

 

ks!

 

 

)

f(t) = h(t)

 

 

 

 

 

 

 

 

L

8

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

y(n) + a y(n 1) + : : : + any = f(t)

 

 

 

 

 

 

0 6 t <1+1; ai комплексные

(5)

 

< y(0) = C0; y0(0) = C1; : : : ; y(n 1)(0) = Cn 1

 

 

 

предположение:

f(t)

оригинал.

 

 

 

Дополнительное

:

 

 

 

 

 

 

 

 

 

 

 

Тогда y(t); y0(t); : : : ; y(n)(t) оригиналы.

 

 

 

 

 

 

Ниже предлагается метод решения такого уравнения:

 

 

 

 

 

 

y(t) ! y~(p); y0(t) ! py~(p) y(0)

 

 

 

 

 

 

 

L

 

 

 

L

 

 

|{z}

=C0

38

ГЛАВА 1. МЕТОДЫ РЕШЕНИЯ НЕКОТОРЫХ УРАВНЕНИЙ.

y00(t) ! p2y~(p) C0p C1; : : :

L

8 f~(p) =

<

(5) , (50):

:

L(p) = pn + a1pn 1 + : : : + an;

y~(p)pn

C

pn 1

 

C

pn 2

 

: : :

 

C

+

+ :1: :

n 1

0

 

)

n

21

0

 

n 1

 

 

~(

 

p

 

 

 

 

 

 

+a p

 

y p

 

 

C

: : : +

 

 

L(p)~y(p) = F (p)

 

 

 

 

 

 

 

 

(500)

F (p) = f~(p) + C0pn 1 + (C0 + a1C0)pn 2 + : : :

Èç (500): y~(p) = FL((pp)). Находя обратное преобразование Лапласа функции решение.

Рассмотрим теперь такую систему:

8

 

d~y

 

= A~y + f~(t)

 

 

 

 

 

 

 

dt

 

a11 : : : a1n

 

 

y1

 

>

 

 

 

 

f1

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

>

~

 

.

; A =

 

: : : : : : : : : : : : :

 

; ~y = .

 

> f =

 

 

>

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

an1 : : : ann

 

 

<

 

 

 

fn

 

yn

 

 

 

 

 

 

~

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

~y(0) = C

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

Дополнительное предположение: âñå fr(t) оригиналы.

fr(t) ! f~r(p)

L

y~(p), получаем общее

(6)

Тогда все yr(t) оригиналы (см. 9)

 

.

= f~(p);

.

= ~y~(p);

 

 

f~1(p)

 

 

 

y~1(p)

 

 

f~n(p)

 

~

y~n(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Èç (6):

 

 

 

 

 

 

 

 

yr0 (t) ! py~r(p) Cr

 

 

 

L

 

 

 

~

~

~

~

(6) ,

~

py~(p) C = Ay~(p) + f(p);

 

 

 

~

 

~y~(p) = (pE A) 1 f~(p) + C~

Такой при¼м решения задач и называется операционным методом.

5. Заключительные замечания:

Обратить внимание на ограничения. Множества рассматриваются на (0; +1). Если не все fr оригиналы, то метод неприменим.

39

Глава 2. Элементы вариационного исчисления

В этой главе все величины вещественные.

Фиксируется класс функций K, рассматриваются функции F (y), где y 2 K, F (y) функционал.

Рассматривается такая задача: найти y : F min. Это и составляет предмет вариационного исчисления.

1. Простейшая задача вариационного исчисления

1.

Фиксируем отрезок t1 6 t 6 t2, y(t),

F (y) =

t2

f

t; y(t); y0(t) dt;

Z

def

 

 

 

 

t1

f задана при t1 6 t 6 t2, y, y0 любые.

Предположение: f непрерывна вместе со своими частными производными до второго порядка включительно.

Определение: Заданы числа y0; y1. Функция y(t) называется допустимой, если она удовлетворяет двум требованиям:

y0(t) непрерывна на отрезке [t1; t2]

y(t1) = y1; y(t2) = y2

Для того, чтобы подставить функцию в интеграл, достаточно первого условия, но мы сужаем класс функций, на которых рассматривается данный интеграл.

Определение: Допустимая функция y(t) дает минимум F , если существует положительное число ", такое, что для любой допустимой функции y~(t), удовлетворяющей условию

maxjy(t) y~(t)j + maxjy0(t) y~0(t)j < "

t t

имеет место F (y) 6 F (~y)

(1) Найти допустимую функцию, дающую минимум F . (Простейшая задача вариационного исчисления)

2. Возьмем y(t) любая непрерывно дифференцируемая функция, на отрезке [t1; t2], не обязательно подчиненная граничным условиям. Составим семейство функций

y(t) + (t);

40

 

 

ГЛАВА 2. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

ãäå (t) íåïð. äèôô. íà [t1; t2];

2 R:

 

!( )

def

 

 

 

= F (y(t) + (t)) = !(0) + !0(0) + o( )

Легко сосчитать, что !(0) = F (y),

 

(t) + 0(t))dt3

!0(0) =

2d Z

f(t; y(t) + (t); y0

 

 

 

 

 

t2

 

 

4

 

d

 

50

 

 

 

t1

 

t2

Z

=(fy + fy0 0)dt (fy = @y@ f(t; y; y0)) è ò.ä.

t1

Определения и обозначения:

(t) = y вариация y(t). (Аналог дифференциала независимого переменного)

!0(0) = F первая вариация F . (Линейная часть приращения функционала, аналог первого дифференциала функции)

t2

Z

F (y) = [fy y + fy0( y)0] dt

t1

Так вычисляется первая вариация функционала по первому дифференциалу независимого переменного (исходной функции).

Это играет такую же роль, как при изучении экстремума в математическом анализе. Вариацией называется малое изменение , малое число.

Определение: Вариация y называется допустимой, если

y t1 = y t2 = 0:

Имеет место: если y(t) допустимая функция, то для того, чтобы вариация y(t) была допустимая, необходимо и достаточно, чтобы y(t) + y(t) была допустимая.

Теорема: Пусть y(t) дает минимум F . Тогда F = 0 для любой допустимой вариации y.

Доказательство: Рассмотрим функции y~(t) = y(t) + (t), где (t1) = (t2) = 0. Тогда для фиксированной функции имеем:

maxt

y~(t) y(t)

+ maxt

y~0

(t) y0(t)

=

maxt

j (t)j + maxt

j 0(t)j

 

 

 

 

 

 

 

 

 

 

Поэтому !( ) достигает минимума при = 0. (Для достаточно малых получаем малое отклонение функции одна от другой.)

Поэтому !0(0) = 0, !0(0) = F (y) = 0: Теорема доказана.