Главная > Лекции по небесной механике (К. ЗИГЕЛЬ, Ю. МОЗЕР)
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
След.
Макеты страниц

Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше

Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике

ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO

Пусть система Гамильтона w˙=JHw удовлетворяет тем же самым условиям, как и в предыдущем параграфе. Следовательно, H= =12wSw+ будет степенным рядом с действительными коэффициентами, начинающимся с квадратичных членов, причем этот ряд сходится в некоторой окрестности точки w=0. Пусть также JS имеет 2n различных собственных значений λ1,,λn,λ1,,λn. Докажем следующую теорему существования:

Пусть λ1 чисто мнимое и пусть ни одно из n1 отношений (λ2/λ1),,(λn/λ1) не является целым числом. Тогда существует семейство действительных периодических решений системы Гамильтона, зависящее аналитически от действительного параметра ρ и обращающееся при ρ=0 в равновесное решение. Период τ=τ(ρ) будет также аналитической функцией ρ и τ(0)=2π|λ1|.

Для доказательства этой теоремы будем искать решения в форме степенных рядов с неизвестными коэффициентами. Сначала будем строить формальные ряды, как уже делалось в §4, а доказательство сходимости этих рядов проведем в следующем параграфе. Переменными будут теперь неизвестные z1,,zm, не имеющие определенных численных значений, в то время как коэффициенты будут какими-то определенными комплексными числами. Если определить равенство, сумму и произведение по правилам, которые имеют место в случае сходимости, то степенные ряды образуют кольцо R(z) без делителей нуля. Если ввести новые переменные ζ1,,ζq посредством подстановки степенных рядов относительно z1,,zm, которые не содержат постоянных членов, то каждый степенной ряд относительно новой переменной ζ посредством перестановки членов перейдет в степенной ряд относительно старой переменной z, и поэтому кольцо R(ζ) изоморфно отобразится на часть кольца R(z). Можно определить частные производные по z1,,zm как результат почленного проведения этой операции, причем дифференцирование многочлена сводится к чисто алгебраическим операциям. Тогда в кольце R(z) получатся обычные правила для производной суммы и произведения, остается в силе и цепное правило.

Сначала рассмотрим опять вместо системы Гамильтона общую нелинейную систему (13;1) при том ограничении, что матрица A имеет два противоположных по знаку собственных значения λ1,λ2=λ1. Попытаемся найти частное решение, в котором x1,,xm разлагаются в степенные ряды по двум неизвестным функциям ξ=ξ(t),η=η(t). Система (13;1) переходит при этом в
xξξ˙+xηη˙=f(x)

Предполагая теперь, что функции ξ,η удовлетворяют двум дифференциальным уравнениям
ξ˙=αξ,η˙=βη,

причем α,β суть степенные ряды по ξ,η, мы нахождение каждого частного решения разделим на два этапа. Прежде всего рассмотрим линейное дифференциальное уравнение в частных производных
xξξα+xηηβ=f(x)

для x(ξ,η) при соответствующем выборе α(ξ,η) и β(ξ,η), а затем уже проинтегрируем систему (1). Преимущество этого способа заключается в том, что при рассмотрении уравнения (2) можно оставаться в кольце формальных степенных рядов и не исследовать вопрос об их сходимости. Уравнение (2) при линейной подстановке x=Cy переходит в уравнение
yξξα+yηηβLy=g(y),

где степенной ряд
g(y)=C1f(Cy)Ly

начинается с членов второго порядка и опять L=[λ1,λ2,,λm]. Чтобы из уравнения (3) путем сравнения коэффициентов можно было однозначно определить m+2 степенных рядов yk(ξ,η)(k=1,,m), α(ξ,η),β(ξ,η), введем еще следующие три условия: все ряды y1ξ, y2η и yk(k=3,,m) должны начинаться с членов второго порядка; y1ξ не должно иметь членов вида ξ(ξη)l и y2=η — членов вида η(ξη)l;α и β должны быть рядами только по степеням одного произведения ξη=ω.m собственных значений λ1,λ2=λ1,λ3,,λm, как и раньше, предполагаются различными; потребуем еще, чтобы никакое из m2 отношений λ3λ1,,λmλ1 не было целым числом.

Подставим m+2 степенные ряда yk,α,β с неопределенными коэффициентами в уравнение (3) и приравняем прежде всего линейные члены. При этом получится, что α,β имеют постоянные члены λ1,λ1. Пусть теперь s — натуральное число, и предположим, что приравниванием членов порядков s однозначно найдено y до s-го порядка и α,β до порядка s1. Приравняем теперь в уравнении (3) коэффициенты членов вида ξpηq(p+q=s+1). Так как g(y) начинается с членов второго порядка, то соответствующий коэффициент в g(y) есть многочлен относительно уже известных коэффициентов в y1,,ym. Пусть γ искомый коэффициент при ξpηq в yk. Этот член дает в соответствующем коэффициенте левой части уравнения (3) величину {λ1(pq)λk}γ. Если теперь не будут одновременно иметь место равенства k=1 и p=q+1 или k=2 и q=p+1, то слева прибавляется еще один многочлен относительно уже известных коэффициентов разложений функций y,α,β. Так как тогда множитель перед γ отличен от нуля, то γ определяется однозначно. Если все-таки или k=1,p=q+1, или k=2,q=p+1, то в этом случае, согласно второму упомянутому выше требованию, γ=0; но, с другой стороны, тогда слева прибавляется для k=1 неизвестный коэффициент при ωp в α и для k=2 неизвестный коэффициент при ωq в β, в то время как остальные члены уже известны. Поэтому предположенное для s доказано также для s+1, и так как оно выполняется для s=1, то отсюда следует по индукции разрешимость уравнения (3) при поставленных условиях.

Теперь допустим, что ряд f(x) имеет действительные коэффициенты. Пусть собственное значение λ1 является чисто мнимым, тогда λ2=λ¯1. Посмотрим, какие следствия выводятся отсюда для рядов y(ξ,η),α(ξ,η),β(ξ,η). Если положить C1C=T и T1y=y, то последняя подстановка в силу уравнений (13;6),(13;7) имеет простейший вид
yl=ρlyk(l=lk;k=1,,m),

и, в частности, y1=ρ¯1y2,y2=ρ¯2y1,ρ¯1ρ2=1. В тождестве (3) перейдем теперь к комплексно сопряженным коэффициентам, причем неопределенные ξ,η могут оставаться постоянными. Тогда, если положить y¯=y¯(ξ,η), то
C1f¯(Cy)=T1C1f(CTy¯),

и, в силу (4), уравнение (3) будет удовлетворяться, если в нем y,α, β заменить на Ty¯,α¯,β¯. Обе первые составляющие вектора Ty¯ суть ρ¯1y¯2=ρ¯1η+,ρ¯2y¯1=ρ¯2ξ+, а остальные начинаются опять с членов второго порядка. Если ξ,η заменить теперь на ρ1η,ρ2ξ, то, очевидно, получится, что Ty¯(ρ1η,ρ2ξ),β¯(ρ1η,ρ2ξ),α¯(ρ1η,ρ2ξ) также составят решение уравнения (3), для которого выполнены все три требования. Вследствие доказанной единственности, отсюда получается, что
Cy(ξ,η)=Cy¯(ρ1η,ρ2ξ),α(ξ,η)=β¯(ρ1η,ρ2ξ).

Таким образом, при ξ¯=ρ1η в случае сходимости значения функций α(ξ,η),β(ξ,η) будут комплексно сопряженными, а значение функции Cy(ξ,η)=x(ξ,η) будет действительным.

Аналогично можно рассмотреть тот случай, когда ряд f(x) имеет действительные коэффициенты и λ1 является действительной величиной, следовательно, λ¯1=λ1,λ2=λ1=λ¯2. Теперь l1=1,l2=2 и можно нормировать, принимая ρ1=ρ2=1, после этого получается y1=y1,y2=y2. Так как тогда обе первые компоненты Ty¯ имеют уже форму y¯1=ξ+,y¯2=η+, то в этом случае
Cy(ξ,η)=Cy¯(ξ,η),α(ξ,η)=α¯(ξ,η),β(ξ,η)=β¯(ξ,η),

и потому в случае сходимости для действительных ξ,η значения функций x(ξ,η),α(ξ,η),β(ξ,η) также будут действительными.

Не будем опять требовать вещественности функции f(x) и рассмотрим частный случай системы Гамильтона. Нужно показать, что тогда α(ξ,η)=β(ξ,η). Использовавшиеся до сих пор обозначения целесообразно изменить следующим образом. Вместо m переменных y1,,ym возьмем 2n переменных zk=xk,zk+n=yk(k=1,,n), причем, в частности, переменные x1,y1 будут играть роль переменных y1,y2 и, соответственно, λ2=λ1 заменим на λn+1=λ1. Предположим, что в функции Гамильтона H уже произведено линейное каноническое преобразование переменных, которое члены второго порядка приводит к нормальной форме, заданной равенством (13;21). Тогда вместо уравнений (3), (4) имеем
zξξα+zηηβ=JHz,

и потому для формально образованного дифференциала dH=Hξdξ+ +Hηdη получается выражение
dH=Hzdz=(αξzξ+βηzη)J(zξdξ+zηdη)=(αξdηβηdξ)zξJzη.

Если положить для сокращения еще zξJzη=Δ, то Hξ=βηΔ, Hη=αξΔ.
αξHξ+βηHη=0,

причем H теперь является степенным рядом по ξ,η.
С помощью уравнения (6) покажем, что H есть ряд только по ω=ξη. Пусть уже доказано, что члены в H порядка s образуют многочлен относительно ω. Для s=2 это справедливо, так как
H=k=1nλkxkyk+=λ1ω+

Пусть теперь γξpηq(p+q=s+1) — член порядка s+1. Так как α=λ1+,β=λ1+, то в тождестве
λ1(ηHηξHξ)=(αλ1)ξHξ+(β+λ1)ηHη

множители αλ1,β+λ1 будут степенными рядами по ω, не содержащими постоянных членов. Применяя здесь метод полной индукции, допустим, что стоящие справа члены порядков s+2 образуют многочлен относительно ω. Коэффициент при ξpηq в левой части равенства будет равен λ1(pq)γ, следовательно, γ=0 при peqq, и этим сформулированное выше утверждение доказывается для s+1.
Так как H зависит только от ω, то
Hξ=ηHω,Hη=ξHω,

и уравнение (6) переходит в
(α+β)ωHω=0.

Так как ωHω=λ1ω+ не будет степенным рядом, тождественно равным нулю, то получаем условие
α+β=0.

Возвратимся еще раз к общему случаю уравнения (2) и допустим, что для найденного решения выполняется также условие (7). В следующем параграфе будет показано, что тогда в случае сходимости функции f(x) ряды y(ξ,η),α,β будут также сходящимися для комплексных ξ,η, модули которых достаточно малы. Поэтому дифференциальные уравнения (1) получают некоторый смысл. Согласно условию (7), будем иметь теперь
ω˙=ξ˙η+ξη˙=(α+β)ξη=0,

следовательно, ω,α,β не зависят от t. Отсюда следует далее
ξ=ξ0eαt,η=η0eβt.

Пусть ряд f(x) имеет опять действительные коэффициенты и λ1 чисто мнимое; выберем начальные значения ξ0,η0 для переменных ξ, η при t=0 в соответствии с условием ξ¯0=ρ1η0 и потребуем, чтобы |ξ0| было достаточно мало. Тогда, согласно ранее полученному результату, числа α=α(ξ,η)=α(ξ0,η0) и β=β(ξ,η)=β(ξ0,η0) будут комплексно сопряженными, следовательно, согласно условию (7), они будут сопряженными и чисто мнимыми. В соответствии с (8), тогда также ξ¯=ρ1η для всех действительных t и |ξ|=|ξ0|; следовательно, согласно упомянутому уже результату, x(ξ,η) также будет действительным. Поэтому (8) представляет семейство действительных периодических решений системы дифференциальных уравнений (13;1), которое содержит комплексный параметр ξ0. Так как правые части дифференциальных уравнений не зависят явно от t, то каждая кривая, изображающая решение, переходит сама в себя, если t заменить на t+c при произвольном c. Поэтому достаточно выбрать ξ0=ρ0. Период имеет величину τ(ρ)=2π/|α|, где α=λ1+; следовательно, τ(0)=2π/|λ1|. Так как в соответствии с нашей заменой y1=ξ+,y2=η+, то получим также, что y1,y2 при достаточно малом ρ>0 действительно зависят от t, и тогда τ(ρ) есть примитивный период. В силу (8) найденные степенные ряды могут быть записаны в следующем параграфе как ряды Фурье по кратным |α|t. Поэтому, после того как будет проведено доказательство сходимости, огсюда будет следовать утверждение теоремы.

Пусть теперь функция f(x) имеет действительные коэффициенты и λ1 действительно, а следовательно, и λ2=λ1 действительно; выберем тогда начальные значения ξ0,η0 действительными и |ξ0|,|η0|
достаточно малыми. В этом случае числа α=α(ξ,η)=α(ξ0,η0) и β=α будут действительными, следовательно, в силу равенств (8) и ξ, η действительны для всех действительных t. Тогда x(ξ,η) сходится, во всяком случае, для достаточно малых |t| и тоже является действительным. Семейство решений (8), зависящее в плоскости ( ξ,η ) от действительных параметров ξ0,η0, есть семейство равносторонних гипербол или, точнее, если различать разные знаки ξ0,η0, то четыре семейства гипербол. Таким образом, мы получили четыре однопараметрических семейства действительных решений системы (13;1). Если теперь α>0 и β=α<0, то eαt (соответственно eβt ) стремится к при t (соответственно при t ). Следовательно, при ξ0η0eq0 в соответствии с равенствами (8) точка ξ,η остается в заданной ограниченной окрестности начал координат только в течение ограниченного времени, в то время как при ξ0=0,η0eq0 и t она, двигаясь по оси η, возвращается в начало координат; аналогичная картина имеет место по оси ξ при ξ0eq0,η0=0,t. Соответственно двум возможным знакам η0 и ξ0 получаются четыре рєшения x(t), которые при t (соответственно при t ) асимптотически приближаются к решению, соответствующему положению равновесия. Напротив, для случая ξ0η0eq0 найденное решение x(t) обладает тем свойством, что оно остается в достаточно малой окрестности начала координат x=0 только в течение ограниченного времени; значит, оно сначала будет некоторое время находиться в этой окрестности, а потом опять уйдет в бесконечность. Впрочем, нельзя заключить, что между этими решениями будут периодические с действительным периодом, так как невозможно установить периодичность, используя только локальное исследование. Если положить eαt=q, то xk(t)(k=1,,m) будут рядами Лорана относительно переменной q и будут иметь по t чисто мнимый период 2πi/α. Чтобы лучше разъяснить этот результат, рассмотрим для сравнения систему
x˙=x+x(xy)g,y˙=y+y(xy)g,

которая образована по аналогии с системой (13;24); здесь g — опять заданное натуральное число. Тогда
(xy)=2(xy)g+1,yx˙xy˙=2xy.

Рассмотрим сначала случай xyeq0, тогда
xy=(a2gt)1/g,x=bye2t

с двумя постоянными интегрирования aeq0,beq0; отсюда
x=b1/2et(a2gt)1/2g,y=b1/2et(a2gt)1/2g,

в то время как при xy=0 получается частное решение x=0,y=cet и x=cet,y=0 с постоянной c. В этом случае λ1=1,λ2=1, но все же нет решения из однопараметрического семейства (9) с комплексным периодом. Отсюда вытекает аналитичность решений системы Гамильтона и для действительного λ1, которая, вообе говоря, не имеет места для системы (13;1).

Для последовательного нахождения коэффициентов разложения в ряды функций y,α,β путем сравнения коэффициентов существенно, чтобы все m2 отношения λ3/λ1,,λm/λ1 не равнялись целым числам. В случае системы Гамильтона это соответствует тому, чтобы ни одно из n1 отношений λ2/λ1,,λn/λ1 не было целым. На примере покажем, что это предположение не является лишним для справедливости теоремы существования. Возьмем функцию Гамильтона в виде многочлена третьей степени
H=12(x12+y12)x22y22+x1y1x2+12(x12y12)y2.

Тогда соответствующая система
{x˙1=y1+x1x2y1y2,y˙1=x1y1x2x1y2,x˙2=2y2+12(x12y12),y˙2=2x2x1y1,

очевидно, имеет равновесное решение x1=x2=y1=y2=0; собственными значениями будут λ1=λ3=i,λ2=λ4=2i. Беря в качестве собственного значения, рассматриваемого в теореме существования, чисто мнимое значениє λ2, удовлетворим предположению, что λ1/λ2=1/2 не будет целым числом, откуда следует существование однопараметрического семейства периодических решений с периодом, приблизительно равным 2πi/λ2=π. Это решение легко получить непосредственно. Действительно, при начальных значениях x1=y1=0 получим в силу теоремы единственности решений дифференциальных уравнений общее решение
x1=0,y1=0,x2=αcos2tβsin2t,y2=αsin2t+βcos2t

с постоянными α,β, являющееся окружностью в плоскости (x2,y2), которая пробегается один раз за время π. При этом радиус остается произвольным. Это как раз то решение, которое получается из теоремы существования, и потому период здесь равен точно π и не только в первом приближении. Но если взять в качестве собственного значения, рассматриваемого в теореме существования, величину λ1, то можно показать, что фактически не существует периодических решений с периодом, близким к 2πi/λ1=2π, если не рассматривать тривиальное решение. Так как в этом случае нарушено предположение теоремы существования о том, что λ2/λ1 не должно равняться целому числу (здесь оно равно двум), то отсюда видно, что указанное предположение не является несущественным. Мы видели уже, что все решения с начальными значениями x1=y1=0 имеют период π. Покажем теперь, что все другие действительные решения не будут периодическими. По теореме о единственности решений дифференциальных уравнений для таких решений, во всяком случае, будет всегда справедливо неравенство p=x12+y12>0. Положим q=x22+y22, тогда из системы (10), если соответствующим образом сгруппировать члены, весьма просто получается дифференциальное уравнение
p¨=4pq+p2.

Вследствие условий p2>0,4pq0 отсюда получается, что p есть выпуклая функция t в строгом смысле, которая, следовательно, не является периодической. Итак, в рассмотренном случае фактически не существует других периодических решений, кроме кругового.

1
Оглавление
email@scask.ru