Лекции по дифференциальным уравнениям (Абрамов А.А
.).pdf10. Матричные формулы решений обыкновенной системы линейных дифференциальных уравнений с п
Теорема: Пусть 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: Теорема доказана.