![](/user_photo/_userpic.png)
книги из ГПНТБ / Митропольский Ю.А. Интегральные многообразия в нелинейной механике
.pdf20 |
ГЛ. I. |
В В Е Д Е Н И Е |
|
|
можно ввести нормальные координаты гь г%.......ZN , для ко |
||||
торых |
|
|
|
|
|
Z É=1 |
1/ - 4 - |
Е юМ |
(1.18) |
|
Z |
A=1 |
|
|
(фf \ |
®k — соответственно нормальные функции и собствен |
ные частоты, соответствующие квадратичным формам (1.16)), и уравнения Лагранжа для невозмущенного движе
ния примут следующий |
вид: |
|
|
= |
0 |
(ft = 1, 2, . . . , N). |
(1.19) |
Допустим теперь, что на рассматриваемую колебатель |
|||
ную систему действуют малые возмущающие силы |
|
||
zQk = е {Qi0) (q{, qc) + 2 |
[Qw? (<7,-, ?;) cos QJ -f |
|
|
a |
|
|
|
+ Q*2)(</i, ^)sinQ efl} |
( t = l , 2 ..........N), |
(1.20) |
где Qa — частоты возмущающих сил, e — малый положи тельный параметр. Тогда, переходя в выражениях (1.20) также к нормальным координатам, получаем следующую систему нелинейных дифференциальных уравнений:
-~pL' + |
- ?zdU zh *i) (ft = 1, 2.......... |
N), (1.21) |
где Zk определяются из условия эквивалентности работ согласно формулам
N |
|
|
( 1.22) |
2 * = = S Q /ФІЛ |
(ft = |
1,2 |
|
/=і |
|
|
|
Уравнения (1.21) путем замены переменных |
|||
+ |
x-bfi-1"**, |
J |
|
Д = |
— mkx_ke~1^ |
(1.23) |
|
j |
|||
легко приводятся к уравнениям в стандартной форме |
|||
-~Ж~ — EX k(t, Xi) ( f t = ± l , |
± 2 , . . . , |
± N), (1.24) |
|
в которых |
—шJ |
|
|
|
|
|
|
|
|
|
(1.25) |
и, кроме того, —со-А = со*, |
= Zk. |
§ 1. И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я |
21 |
В дальнейшем для упрощения систему дифференциаль ных уравнений (1.24) целесообразно записывать в вектор
ной форме |
|
- § - = е * (/,*), |
(1.26) |
где X , X — n-векторы (в нашем случае п — 2N). Правые части системы (1.26) являются периодическими функциями или почти периодическими в зависимости от того, как зави сели от t функции eQk.
К уравнениям типа (1.26) могут быть приведены также задачи о воздействии на колебательные системы внешних сил высокой частоты. Так, например, пусть движение та кой системы характеризуется уравнениями вида
Ц |
Чь 4 г ) |
(s = |
1, 2, |
... . N), (1.27) |
где Fs (wf, qlt |
— периодические |
(или |
почти периоди |
ческие) функции t и о — «большой» параметр.
Тогда, вводя новые неизвестные согласно формулам
* = |
qs = xs, |
-jf- = y s (s= |
1, 2, |
. . . , N), (1.28) |
вместо системы (1.27) получаем систему уравнений |
||||
" $ Г |
= еу*’ |
= s f ^ x’ Хі’ Уі> |
(s = |
Ь 2...........N)’ |
|
1 |
|
|
(1.29) |
где e = |
параметр. |
|
|
|
—---- малый |
|
|
К уравнениям вида (1.26) может быть приведен и более общий случай, чем указанный выше, когда состояние дина мической системы характеризуется угловой переменной а и п переменными xlf х2.......хп и описывается следующей
системой |
уравнений: |
|
|
|
|
|
dx |
xs(об, х^у |
. • • t хп) |
(s |
1> 2, |
. . . j |
|
dt |
||||||
da |
|
|
|
|
(1.30) |
|
= Ы (xlf . . . |
, хп) + |
А (а, |
х1г |
, хп), |
||
dt |
||||||
где X — большой параметр, X s (а, |
хіу |
хп), А (а, хи ... |
...,хп) — периодические (или почти-периодические) функции
22 |
ГЛ. I. В В Е Д Е Н И Е |
угловой |
переменной а с периодом 2л. (Системы типа (1.30) |
встречаются при изучении гироскопических систем, в тео рии ускорителей и т. д.)
Вводя новую независимую переменную т = Kt, можем
представить систему (1.30) в виде |
|
|
|
|
|
|
|
|
|||||||
-fö- = |
eXs{a, |
xlt |
, хп) |
(s = |
1, 2......... п), |
|
(1.31) |
||||||||
~ |
= |
со(х1, |
... , хп) + еЛ(а, |
|
. . . |
, |
хп), |
|
|
||||||
Ху, |
|
|
|
||||||||||||
где е = |
1---- малый |
параметр, |
или, |
исключая |
т, |
в |
виде |
||||||||
dx s |
= |
eX*s (а, |
Ху, |
. . . , хп, |
е) |
(s — |
1, |
2, |
... |
, п), |
|||||
da. |
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.32) |
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
X s (a, |
Ху, |
. . . |
, |
х п) |
|
|
|
|||
К (а, Ху, |
. . . , хп, е) |
|
|
|
|
||||||||||
Ö) (хх, |
. . . . |
х п) + 8Д (а, |
Ху, |
. . . , |
х п) ' |
||||||||||
|
|
|
|
|
К уравнениям вида (1.26) могут быть сведены также уравнения более общего вида, чем (1.21).
Как известно, нестационарные процессы в нелинейных колебательных системах во многих случаях приводятся к рассмотрению систем дифференциальных уравнений типа
-4 " { 2 ац 00 <7;} + |
2 |
Ьц (т) <7 ( = |
|
|
= eQj (т, Ѳ, qy, ... |
, qN, qt, ... , qj4, |
e) |
( / = 1 , 2 , . . . , M), |
|
|
|
|
|
(1.33) |
где qt — обобщенные |
координаты, |
eQj |
— внешние возму |
|
щающие силы, т — et — медленное |
время (медленное по |
сравнению с естественной единицей времени — порядка пе риода собственных колебаний), е — малый положительный параметр, а;/ (т) = щі (т), Ьц (т) = Ьц (т) — инерционный и квазиупругий коэффициенты, медленно изменяющиеся со
временем, |
dQ |
— |
, |
, |
|
|
|
ѵ (т). |
|
|
|||
Система типа (1.33) с помощью замены |
переменных со |
|||||
гласно формулам |
|
|
|
|
||
|
|
N |
|
|
|
|
|
<7і = |
2 |
Ф |
(t' = 1, 2, ... |
, N), |
(1.34) |
k—\
§ ! И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я |
23 |
где tp£fe>(т) — нормальные функции, зависящие от парамет ра т и удовлетворяющие условию ортогональности, может быть приведена к следующему виду:
-ZjjT + |
Zk(т’ ѳ- 2i> • ■■» г‘Ѵ, Zj......... zN, e) |
|||
|
( * = 1 , 2 ......... N), |
(1.35) |
||
после чего, произведя замену согласно формулам |
|
|||
г* = Хке1Ы |
^ ‘ + |
х . ке-^ |
(1.36) |
|
г, = «о*(т) Ѵ Ь |
(Т,Л - |
|||
|
||||
приходим к системе уравнений вида |
|
|||
|
|
= еХ (т, Ѳ, X, е) |
(1.37) |
|- ^ - = ѵ(т), т = eij , в которой правые части зависят от
медленного времени х и являются периодическими, почтипериодическими или квазипериодическими по Ѳ в зависи мости от того, как зависят от Ѳфункции Q,-.
Заметим, что к системе уравнений вида (1.37) может быть приведена и более общая система дифференциальных урав нений, описывающая колебания системы со многими степе нями свободы, характеризующиеся в невозмущенном со стоянии (т. е. при е = 0, х — const) лагранжевой функ
цией
j Г N N N 1
£ = -Ö- | 2 ■аіі(х)йіЯі + 2 |
g iiW q â — 2 |
bij(x)qiqi\, |
||
1 |
U ',/= i |
> ./= 1 |
i , j = i |
) |
где gij |
(г) = |
gji (T). |
|
(1.38) |
|
|
Кроме рассмотренных классов дифференциальных урав нений, содержащих малый параметр (в основном приводя щихся к уравнениям в стандартной форме), с помощью метода интегральных многообразий можно эффективно исследовать несколько более общие дифференциальные си стемы, описывающие динамические процессы в сложных колебательных системах.
В связи с этим мы будем рассматривать также системы, близкие к автономным, вида
~ = Х(х) + |
X, е), |
(1 .39) |
24 |
ГЛ. 1. В В Е Д Е Н И Е |
|
и системы, близкие к неавтономным: |
|
|
|
= X (/, х ) + еХ* (yt, X , е). |
(1.40) |
К системам типа (1.39) и (1.40) при ряде предположений приводятся многочисленные задачи теории колебаний (од ночастотные колебания систем со многими степенями сво боды, не близкие к гармоническим, колебания релаксацион ного типа и т. д.).
Будем также рассматривать и более общие системы диф
ференциальных уравнений вида |
|
|
|
= |
X (т, х) -f еХ* (т, 0, X, |
е), |
(1.41) |
= |
X (т, Ѳ, х ) + еХ* (т, 0, |
X, г ). |
(1.42) |
Многие задачи приводят к рассмотрению колебательных систем, описываемых дифференциальными уравнениями, в коіорых одни переменные изменяются быстро, а другие медленно. К уравнениям такого типа относятся, например, указанные выше уравнения (1.30).
В общем виде такие уравнения можно записать следую щим образом:
~ ~ = X (х , у) + |
(t, X, у), |
(1.43)
=(*’ х’ У)-
где X , у — соответственно п- и m-векторы, е — малый по ложительный параметр.
Проблемы, приводящие к рассмотрению уравнений типа (1.43), встречаются в динамике спутников и искусственных небесных тел. Частным случаем системы (1.43), когда пе ременная X отсутствует, являются уравнения в стандартной форме. При Y (t, X , у) £= 1 приходим к рассмотрению урав нений типа (1.41).
Как известно, ряд актуальных задач радиотехники, ав томатического регулирования, химической кинетики и др. приводит к исследованию систем дифференциальных урав нений с малым параметром при старших производных. В связи с этим нами будут рассмотрены системы нелиней-
§ 1. И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я |
25 |
ных нерегулярно-возмущенных дифференциальных уравне ний вида
= fix, г, i), е |
= F (х, z, t), |
(1.44) |
где е — малый параметр, х, / — «-векторы, г, F — т-век- торы, а также системы линейных нерегулярно-возмущенных дифференциальных уравнений
■^r = A(t)x + B(t)y + h(t),
(1.45)
е= С (t) у + eF (t) X + еН (t),
где Ху h — m-векторы, у, Н — п-векторы, А, В, С, F — со ответственно (m X т)-, (т X п)-, (п X «)-, (п X т)-мат- рицы.
В связи с большим интересом, возникшим в последнее время к теории дифференциальных уравнений с запазды вающим аргументом, будет уделено внимание также рассмот рению интегральных многообразий для уравнений с запаз дыванием вида
|
j £ - = ALt)x(t) + A1( f ) x ( t - A ) |
+ |
|
||||
|
+ X(t,g(t), |
x(t), |
X (t — A), у (t), |
y(t — Д), e), |
|||
|
(ef) {/(() + |
(eO y(t — b) + |
(1.46) |
||||
|
+ eY (t, g (t), |
X (t), |
X (t — A), у (t), у (t A), e), |
||||
|
= © (t) + |
G (t, |
g (t), |
X (t), |
у (t), e), |
|
|
где X , |
у, g — соответственно n - , m - , |
/г-векторы, А |
(/), A t (/), |
||||
В (et), |
Bi (et) — ограниченные |
квадратные |
матрицы, |
© (t) — fe-вектор-функция, е — малый положительный па раметр.
Будут рассмотрены также нерегулярно-возмущенные уравнения с запаздыванием
~= f(t, х,у, е),
8 ЧГ = А V' х>Ха) у + В У>х>Ха) У* + § ((>х’ А'Д' У' Уь>е)’
(1.47)
26 |
; |
|
ГЛ. I. В В Е Д Е Н И Е |
где |
х — т-вектор, |
у — п-вектор, хд — х (t — еД), //д = |
|
= |
у |
(t — еД), е > |
0 — малый параметр. |
Как уже указывалось, рассмотрение исходной системы уравнений на многообразии сводится к рассмотрению урав нений, число которых равно размерности многообразия, в результате чего исследование решений исходной системы уравнений значительно упрощается.
В связи с этим большое значение приобретают развитые методы для исследования нелинейных дифференциальных уравнений в бесконечномерных пространствах. Использо вание спектральной теории линейных операторов позволи ло развить теорию интегральных многообразий для различ ных классов нелинейных дифференциальных уравнений в бесконечномерном банаховом пространстве.
§2. Вспомогательные сведения из линейной алгебры
ианализа
Приведем некоторые вспомогательные сведения из линейной алгеб ры и анализа, которые играют важную роль в теории дифференциаль ных уравнений.
1. Матричные обозначения. Рассмотрим квадратную матрицу А порядка п с элементами аи- (і, / = 1, .... п), ко торые, вообще говоря, являются комплексными. Матрица (кі — А), где к — независимое переменное, / — единичная матрица порядка п, называется характеристической мат рицей для А. Ее определитель
Ф (Я) ~ det (kl — А) |
(2.1) |
является многочленом относительно к и называется харак теристическим многочленом матрицы А. Корни характери стического многочлена матрицы называются ее характери стическими числами или собственными значениями. Обоз
начив эти корни kt |
(і = 1, ..., л), |
запишем: det (kl — А) |
= |
П |
|
|
|
= п(к —к). |
|
|
|
1 |
собственных |
значений) матрицы |
А |
Спектр (спектр |
есть множество всех ее собственных значений. Он состоит
из тех значений Я,, для которых матрица (kl — А) |
необрати |
ма. Те значения к, для которых матрица (kJ — А) |
обратима, |
называются регулярными, а матрица (kJ — Л)-1 называется
§ 2. С В Е Д Е Н И Я ИЗ Л И Н Е Й Н О Й А Л Г Е Б Р Ы И А Н А Л И З А |
27 |
резольвентой матрицы Л. Спектр п собственных значений матрицы А порядка п совпадает с множеством корней ал гебраического уравнения
det (Я/ — А) =-- 0 |
(2.2) |
(характеристического уравнения матрицы А), или в развер
нутом виде: |
— а, |
|
|
'X — ап |
Мп |
|
|
|
М2 |
|
|
det |
|
— Я, |
(2.3) |
|
= 0. |
||
С1цЛ |
Сіп |
|
|
Кратность каждого корня Я/ этого уравнения равна его
алгебраической кратности т,- как собственного значения,
так что т\ + т-2 + ... = п.
Для любой квадратной матрицы А с собственными зна чениями Я/ матрица аА имеет собственные значения аЯ;;
матрица А р — собственные значения Я/ (р — 0, 1, 2,...); каждый многочлен или аналитическая функция / (Л) име ет собственные значения / (Я/).
Если две квадратные |
матрицы Л |
и ß имеют |
собствен |
||
ные значения Я/ и |
то |
спектром |
прямого произведения |
||
Л 0 |
В является |
множество всевозможных |
произведе |
||
ний |
Я/рй. |
|
|
|
|
Имеет место следующий фундаментальный результат о |
|||||
канонической форме матрицы. |
|
|
|||
Каждая квадратная матрица Л порядка п с помощью |
|||||
преобразования подобия |
|
|
|
||
|
|
J = 5Л5~' |
|
|
(det S ф 0) может быть приведена к нормальной жордановой форме
У = diag [Ух (ЯА), . . . , Jrn(Ят)] (т < п), |
(2.4) |
где |
|
К |
1 |
0 |
•• • |
0 |
|
о |
яр 1 |
.. • |
0 |
|
|
Jр(Яр) |
|
|
|
(р) |
(2.5) |
|
|
|
|
||
0 |
0 |
0 |
• |
1 |
|
0 |
0 |
0 |
.. • |
я„ |
|
28 |
ГЛ. |
I. В В Е Д Е Н И Е |
(р ■= 1, |
m; li > 1, |
lm > 1), N\p) — нильпотентная |
матрица, причем каждому характеристическому корню Хр
кратности а р соответствует |
одна или |
несколько клеток |
||
Жордана размерами 1{р\ |
.... Ір], |
таких, |
что |
|
4 *4" |
• • • |
+ |
Ір*— |
|
Каждой клетке Жордана |
J р (Хр) порядка /р с точностью |
до ненулевого скалярного множителя отвечает один и толь ко один собственный вектор матрицы Л, причем различным кьеткам Жордана соответствуют линейно-независимые соб ственные векторы.
Степень вырождения г собственного значения Хр пред ставляет собой максимальное число линейно-независимых собственных векторов матрицы А , соответствующихХр. В об щем случае г С ар. Если степень вырождения г собствен ного значения равна его кратности, т. е. г = ар, то, оче
видно, /р*= ... = t'p = 1. В этом случае все соответствую щие клетки Жордана будут содержать по одному элементу (простые клетки).
Пусть Хі,..., Хт— характеристические числа матрицы Л, соответствующие различным клеткам Жордана (не обяза тельно различные между собой). Тогда, так как XI —
= S~lXS, А — S~lJS, |
то |
характеристический полином |
||
Д (X) матрицы Л можем записать в виде |
|
|||
Д (X) = det (XI — J) = (X — XJ1* . . . (X — X jm |
|
|||
(h + |
’ ' • |
+ |
Ал “ n)‘ |
(2.6) |
Множители (X — XP)1P (p = |
1, .... m) называются |
эле |
ментарными делителями матрицы Л, а натуральные числа Ір, т. е. размеры (порядки) клеток Жордана,— показателя ми элементарных делителей, соответствующих числу Хр.
Если все характеристические числа Хр имеют простые |
|
элементарные делители (Ір — 1), то матрица Жордана J бу |
|
дет чисто диагонального вида |
|
Хг |
О |
J = |
К |
= diag (Xv . . . , Ха), (2.7) |
ОХп
причем числа Хр (р = 1.......п) не обязательно различны.
§ 2. С В Е Д Е Н И Я ИЗ Л И Н Е Й Н О Й А Л Г Е Б Р Ы И А Н А Л И З А |
29 |
2. Линейные преобразования. Каждой квадратной мат рице А порядка я соответствует линейное преобразование А
векторного координатного пространства Сп размерности п,
т. е. вектору х = (хи ..., хп) пространства Сп ставится в соответствие вектор
А х = у = ( У ! .......... |
У п ) , |
(2.8) |
определяемый соотношением
п
У і = 2 a ikx k.
k=i
Нулевой матрице О соответствует нулевое преобразова ние О, переводящее каждый вектор в нуль; единичной матрице соответствует тождественное преобразование. Пре образование является примером линейного оператора, дей
ствующего в Сп.
Говорят, что подпространство Сп линейного простран ства Сп инвариантно относительно преобразования А , если каждый вектор из Сп преобразование А снова переводит в некоторый вектор из Сп , т. е. А С п сг Сп .
Если пространство Сп распадается в прямую сумму под пространств, инвариантных относительно линейного пре образования А , то в надлежащей координатной системе матрица преобразования А принимает клеточно-диагональ ный вид и диагональные клетки являются матрицами пре образований, индуцированных преобразованием А в ин вариантных подпространствах.
Характеристическим многочленом преобразования А на зывается характеристический многочлен <р (Я) = det (XI —
— А) матрицы А.
Степень характеристического многочлена преобразова ния Л равна порядку матрицы А, а порядок матрицы равен
размерности пространства Сп. Поэтому степень характерис тического многочлена преобразования А равна размернос ти пространства, в котором действует это преобразование.
Так же, как и квадратная матрица, каждое линейное преобразование является корнем своего характеристиче ского многочлена, т. е. q> (А ) = 0.
Число % называется собственным значением линейного преобразования А , если в рассматриваемом пространстве