Главная > Лекции по небесной механике (К. ЗИГЕЛЬ, Ю. МОЗЕР)
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
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

В первом параграфе этой главы рассматривался метод определения периодических решений системы Гамильтона с помощью степенных рядов. В двух предыдущих параграфах были получены этим методом периодические решения плоской задачи трех тел, которые имеют значение в теории движения Луны. В этом параграфе будет рассмотрен третий метод определения периодических решений системы дифференциальных уравнений. Бо́льшая часть изложенных ниже результатов имеет место не только при регулярности, но и при более слабых предположениях; все же ради простоты предположение о регулярности будет в дальнейшем сохранено.
Рассмотрим систему дифференциальных уравнений
x˙k=fk(x,α)(k=1,,m),

которая зависит от одного параметра α. Пусть правые части fk при
|xlξl|<r(l=1,,m),αG

являются регулярными функциями m+1 комплексных переменных xl и α, причем G есть область комплексной α-плоскости. Далее, пусть при условиях (2) выполняются неравенства
|fk(x,α)|M(k=1,,m).

Прежде чем применять метод малого параметра, исследуем зависимость решений системы (1) от параметра α и начальных значений ξ1,,ξm.

Пусть ξ1,,ξm — какие-нибудь комплексные величины, удовлетворяющие условиям
|ξlξl|<r2(l=1,,m).

Тогда при
|xlξl|<r2,αG
fk(x,α) будут регулярными функциями переменных xl и α, и в этой области справедлива оценка (3). Согласно теореме существования Коши (§4), система (1) имеет единственное решение x(t,ξ,α), для которого x(0,ξ,α)=ξ и xk(t,ξ,α)(k=1,,m) суть регулярные аналитические функции комплексной переменной t в круге
|t|<r2(m+1)M=ρ.

Последнее справедливо для любого значения ξ из области (4) и для любого α из G. Покажем теперь, что xk(t,ξ,α) в области
|t|<ρ,|ξlξl|<r2(l=1,,m),αG

будут регулярными функциями всех m+2 независимых комплексных переменных t,ξl,α. Это следует из доказательства теоремы Коши, данного в §4. Коэффициенты αkn разложения в ряд xk(t,ξ,α) по степеням t при сравнении коэффициентов оказываются многочленами относительно коэффициентов разложения Тейлора функции fl(x,α) ( l= =1,,m) по степеням x1ξ1,,xmξm; эти последние коэффициенты по формуле Тейлора будут аналитическими функциями от ξ1,,ξm,α в области |ξhξh|<r2(h=1,,m),αG. Так как, с другой стороны, для разложений функций xk(t,ξ,α) по степеням t можно указать мажорирующие функции, коэффициенты которых зависят только от M и r, то эти ряды сходятся равномерно по ξh и α в каждом открытом отрезке |t|<ρ. Следовательно, xk(t,ξ,α) по известной теореме Вейерштрасса будут регулярными по всем m+2 переменным в области (5).

Если значения xk(t,ξ,α) при каком-либо t из интервала 0<t<ρ, например при t=ρ/2, считать опять начальными значениями, то решение можно продолжить аналитически и за точку t=ρ. Пусть решение x(t,ξ,α) при закрепленных ξ=ξ,α=α продолжено, как функция t, на весь интервал 0t<t1. Если тогда кривая x(t,ξ,α) при 0tt1 вся лежит в области регулярности f1(x,α),,fm(x,α), то, как это следует из теоремы о покрытии, из последовательного применения вышеупомянутых операций следует, что в достаточно малой окрестности U точки ξ=ξ,α=α решение x(t,ξ,α) может быть продолжено на интервал 0tt1 и там будет оставаться регулярной функцией всех переменных t,ξ,α. Но чем больше будет выбрано t1, тем меньше будет вообще окрестность U, в которой регулярность сохранится при аналитическом продолжении. Нужно заметить, что это рассуждение можно провести и для таких дифференциальных уравнений
x˙k=fk(x,t,α)(k=1,,m),

в правые части которых входит явно независимая переменная t. Именно, если ввести новую неизвестную x0 и заменить систему (6) системой
x˙0=1,x˙k=fk(x,x0,α)(k=1,,m),

то правые части этих m+1 дифференциальных уравнений не содержат более переменной t.

Так как xk(t,ξ,α) будет при 0tt1 регулярной функцией ξl в окрестности ξ=ξ, то там, в частности, существуют частные производные xkξl(t,ξ,α). Ввиду того что дальнейшее рассмотрение не зависит от α, ибо будет предполагаться, что α имеет постоянное значение, мы не будем писать α в качестве аргумента функций. Если известно теперь решение x(t,ξ) для системы фиксированных начальных значений ξl(l=1,,m), то частные производные xkξl=xkξl(t,ξ) определяются следующим образом из так называемых уравнений в вариациях. Так как ξl,t можно рассматривать по отношению к xk(t,ξ) как независимые переменные, то из дифференциальных уравнений (1) дифференцированием по ξl получаем
x˙kξl=r=1mfkxrxrξl,fkxr=fkxr[x(t,ξ)](k=1,,m;l=1,,m).

Если ввести матрицы m-го порядка X=(xkξl)F=(fkxl), то для X получается уравнение в вариациях
X˙=FX,

причем F известно. Так как x(0,ξ)=ξ, то при t=0X=E. Итак, матрицу X=X(t) можно получить интегрированием линейного дифференциального уравнения (7) при начальном условии X(0)=E. Это интегрирование можно провести последовательными приближениями с использованием интегрального уравнения
X=E+0tFXdt

для решения которого можно использовать матричную последовательность
X0=E,Xn=E+0tFXn1dt(n=1,2);

можно также использовать метод сравнения коэффициентов, как это было сделано при доказательстве теоремы существования Коши.
Для определителя Δ=|X| получается
Δ˙=k,l=1mx˙kξlXlk

причем Xkl есть алгебраическое дополнение элемента xkξl в матрице X, т. е. минор элемента xxξl, взятый со знаком (1)k+l. Если ввести еще матрицу Y=(Xkl), составленную из алгебраических дополнений элементов X, то равенство (8) можно записать в виде
Δ˙=σ(X˙Y),

где σ обозначает след матрицы X˙Y. С помощью уравнения (7) имеем
Δ˙=σ(FYY).

С другой стороны, XY=ΔE, откуда
Δ˙=Δσ(F)=Δσ,

где
σ=σ(F)=k=1mfkxk(x)

и x=x(t,ξ). Принимая во внимание, что начальное значение Δ(0) функции Δ=Δ(t,ξ)=Δ(t) вследствие X(0)=E равно единице, путем интегрирования уравнения (9) получим
lnΔ=0tσdt.

Можно считать, что система (1) есть система дифференциальных уравнений движения потока жидкости, причем xk(k=1,,m) суть координаты частиц жидкости. Так как правые части не содержат явно независимой переменной, то имеет место установившееся движение. При t=0 положение частиц жидкости определяется координатами ξk. По прошествии времени t частицы перейдут из ξ в x(t,ξ), чем устанавливается отображение ξ на x. Функциональная матрица этого отображения есть (xkξl)=X(t,ξ), и функциональный определитель равен Δ. Если Δ тождественно равно единице, получается отображение, сохраняющее объем, что соответствует несжимаемому потоку. В соответствии с уравнением (9) это означает, что
σ=k=1mfkxk=0

Для системы Гамильтона
x˙k=Eyk,y˙k=Exk(k=1,,n)

имеем
σ=k=1n[(Eyk)xk+(Exk)yk]=0,

так что в этом случае равенство (10) выполняется.

Метод малого параметра, предложенный Пуанкаре [1], возник из следующей задачи. Рассмотрим решение x(t,ξ,α) системы (1) опять в зависимости от ξ и α и допустим, что при α=α система имеет периодическое решение. Пусть этому решению соответствуют начальные значения ξ=ξ, тогда x=x(t,ξ,α). Предположим при этом, что речь идет не о равновесном решении. Пусть τ>0 будет периодом x(t,ξ,α) по t, причем это необязательно наименьший положительный период, и пусть вся кривая x(t,ξ,α) лежит при 0tτ в области регулярности функций f1,,fm по x и α. Тогда эти утверждения справедливы для всякого действительного t, так как
x(t+τ,ξ,α)=x(t,ξ,α).

По теореме о единственности решений дифференциальных уравнений равенство (11) справедливо при всех t, если оно верно хотя бы при одном значении, например при t=0. Поставим вопрос, имеет ли система (1) периодические решения для несколько измененных начальных значений ξ,α.

Будем искать сначала периодические решения с тем же самым периодом τ. Чтобы x(t,ξ,α) имело период τ, по теореме единственности необходимо и достаточно, чтобы x(τ,ξ,α)=x(0,ξ,α)=ξ. Если положить
φk(ξ,α)=xk(τ,ξ,α)ξk,

то, следовательно, нужно удовлетворить m аналитическим уравнениям
φk(ξ,α)=0(k=1,,m).

Это будет система неявных уравнений, которые удовлетворяются при ξ=ξ,α=α в силу периодичности исходного решения. Следовательно, если функциональный определитель |φkξl| порядка m отличен от нуля при ξ=ξ,α=α, то можно найти решение системы (13) вблизи α=α, и по известной теореме существования для неявных функций m разностей ξkξk получаются в виде степенных рядов по αα, причем эти ряды не содержат постоянных членов. Но этого здесь как раз и не может быть, так как именно определитель |φkξl| обязательно равен нулю. Эту трудность, однако, можно обойти небольшим видоизменением рассуждения. Исследуем, почему определитель должен обращаться в нуль. Пусть α=α. Если опять положить X=[xkξl(t,ξ)] и ввести матрицу C(t,ξ)=XE, то в силу равенства (12) имеем для функциональной матрицы
φkξl=C(τ,ξ)=C.

С другой стороны, если ξ будет какой-нибудь точкой на траектории x(t,ξ), то при соответствующем выборе t можно положить, что
ξ=x(t,ξ),
т. е. величина ξ будет функцией t. Так как правые части дифференциальных уравнений (1) явно не зависят от t, то
x(t+t,ξ)=x(t,ξ).

Дифференцируя равенство (16) по t, в соответствии с (1) и (15) получим уравнения
fk[x(t,ξ)]=l=1mxkξl(t,ξ)fl(ξ)(k=1,,m),

следовательно,
f(x)=Xf(ξ),f(x)f(ξ)=C(t,ξ)f(ξ),

где f(x) — вектор, имеющий составляющие f1(x),,fm(x), и где нужно положить x=x(t,ξ). Тогда, в частности, при ξ=ξ,t=τ будем иметь f(x)=f(ξ), и поэтому
Cf=0,f=f(ξ).

Так как рассматриваемое периодическое решение x(t,ξ) не является равновесным, то f(ξ) не будет нулевым вектором, следовательно, |C|=0. Поэтому основная причина обращения определителя C в нуль состоит в том, что при произвольном сдвиге начальных значений ξ на траектории x(t,ξ) получается опять периодическое решение, а именно, та же самая траектория, для которой только t увеличено на постоянную t. Этого можно избежать, если варьировать начальные значения ξ только в (m1)-мерной плоскости, которая не касается интегральной кривой в начальной точке ξ. Мы уже видели, что f\left(\xi^{*}\right)неестьнулевойвектор,поэтомуможновыбратьобозначениятак,чтобыпоследняясоставляющаяf_{m}\left(\xi^{*}\right)
eq 0.Тогдавсоответствииссистемой(1)x_{m}=\xi_{m}^{*}ибудеттакойплоскостью.Следовательно,можноположить\xi_{m}=\xi_{m}^{*}иварьироватьтолькоm-1начальныхзначений\xi_{1}, \ldots, \xi_{m-1}.Нотаккакдолжныбытьудовлетвореныmуравнений(13),мыбудемтеперьпредполагатьпериодискомогопериодическогорешения\tau$ также величиной переменной.
Если положить теперь
φk(τ,ξ,α)=xk(τ,ξ,α)ξk,

то m уравнений
φk(τ,ξ,α)=0(k=1,,m)

должны быть удовлетворены при дополнительном условии ξm=ξm. Они имеют очевидное решение τ=τ,ξ=ξ,α=α. Будем рассматривать ξ1,,ξm1 и τ в уравнениях (19) как неизвестные, а α как независимую переменную, и найдем соответствующую функциональную матрицу B порядка m, которая получается из C заменой столбца φξm, соответствующего ξm, столбцом φτ=φτ(τ,ξ,α). Теперь, в соответствии с системой (1) и соотношением (18),
φτ=x˙(τ,ξ,α)=f[x(τ,ξ,α),α],

следовательно, в точке τ=τ,ξ=ξ,α=α
φτ=f(ξ,α)=f

и
B=(φξ1φξm1f),

где первые m1 столбцов суть
φξk=φξk(τ,ξ,α)(k=1,,m1).

Если определитель |B| отличен от нуля, то систему уравнений (19) можно при ξm=ξm разрешить в окрестности α=α относительно τ, ξ1,,ξm1 и получить для разностей ττ,ξ1ξ1,ξm1ξm1 степенные ряды по αα, не содержащие постоянных членов. Следовательно, тогда для всех значений параметра α в достаточной близости от α можно определить такие начальные значения ξ1,,ξm1,ξm= =ξm и такой период τ, что соответствующее этим начальным значениям решение будет периодическим с периодом τ. Чтобы вычислить матрицу B, нужно в соответствии с равенством (14) проинтегрировать только систему линейных уравнений в вариациях (7), причем x=x(t,ξ,α) следует взять в качестве известного исходного периодического решения с начальными значениями ξ=ξ,α=α. Можно построить простые примеры, показывающие, что определитель |B| в отличие от определителя |C| не всегда равен нулю.

Пуанкаре распространил свой метод и на общий случай, когда правые части дифференциальных уравнений зависят явно от t; при этом правые части должны быть, однако, периодическими функциями t. Предполагается, что существует периодическое решение с тем же самым периодом; легко показать на примерах, что аналогично подсчитываемый определитель по крайней мере не всегда равен нулю. Последнее правдоподобно, так как для обоснования равенства нулю определителя существенную роль у нас играла стационарность потока. Мы не будем больше здесь и далее углубляться в важные и интересные вопросы, связанные с теорией дифференциальных уравнений с периодическими коэффициентами. Бо́льшую часть известных методов и результатов этой теории можно истолковать с помощью рассмотренных нами стационарных потоков; кроме того, при начальном рассмотрении не решенной еще задачи следует ограничиваться разбором простых нетривиальных случаев.

Мы покажем теперь, что можно найти периодическое решение для значения α, близкого к α, с тем же самым периодом τ=τ, как и в исходном решении, если известен не зависящий от t интеграл ψ(x,α), который при x=ξ,α=α не является стационарным. При этом предположим, что интеграл ψ(x,α) будет аналитическим относительно x, α в окрестности исходного периодического решения x(t,ξ,α) и значения α. Являясь интегралом системы (1), ψ=ψ(x,α) удовлетворяет уравнению в частных производных
k=1mψxkfk(x,α)=0

тождественно по x и α. Если ввести теперь вектор-строку ψx с составляющими ψxk(ξ,α), то, в частности,
ψxf=0.

Мы предположили, что ψ не является стационарным при x=ξ, α=α; поэтому ψx не является нулевым вектором. Так как, с другой стороны, fm(ξ,α)eq0, то не все величины ψxk(ξ,α) при k=1,,m1 равны нулю. Выберем теперь обозначения таким образом, чтобы ψxm1(ξ,α)eq0. Так как интеграл ψ(x,α) имеет постоянное значение на каждой траектории, то равенство
ψ[x(t,ξ,α),α]=ψ(ξ,α)

удовлетворяется тождественно по t,ξ,α. Отсюда, дифференцируя по ξl, получим
k=1mψxk(x,α)xkξl=ψxl(ξ,α)(l=1,,m)

на каждой траектории x=x(t,ξ,α). Если здесь положить t=τ,ξ= =ξ,α=α, т. е. если x=ξ, то в векторной форме получим
ψxX(τ,ξ)ψx=ψxC=0.

Из уравнений (14), (20), (21) и (23) прежде всего следует ψxB=0, так что в случае существования нестационарного интеграла определитель |B| равен нулю, и поэтому нельзя прямо применить метод, развитый выше. Условия того, что x(t,ξ,α) будет периодическим решением с периодом τ, даются m уравнениями (13). Положим опять ξm=ξm, и тогда нужно определить еще m1 неизвестных ξ1,,ξm1, удовлетворяющих этим m уравнениям. Разрешим прежде всего m1 уравнений
φk(ξ,α)=0(keqm1)

и выразим ξlξl(l=1,,m1) в виде степенных рядов по αα, не содержащих постоянных членов, причем предположим, что соответствующий функциональный определитель не равен нулю при ξ=ξ, α=α. Соответствующая функциональная матрица m1 порядка A получается из C вычеркиванием последнего столбца и предпоследней строки. Тогда при |A|eq0 в окрестности α=α уравнения (24) удовлетворяются. Остается показать, что вследствие существования интеграла ψ(x,α) выполняется также и последнее уравнение φm1(ξ,α)=0. Если образовать с найденными начальными значениями ξ1,,ξm1 и ξm=ξm решение x=x(t,ξ,α), то на этой траектории удовлетворяется уравнение (22). Положим здесь, в частности, t=τ и применим теорему о среднем значении из дифференциального исчисления к функции ψ(x,α) и переменной xm1. Тогда, принимая во внимание уравнения (24), получим
0=ψ[x(τ,ξ,α),α]ψ(ξ,α)=ψxm1(x~,α)φm1(ξ,α),

где x~k=ξk(keqm1), и x~m1 лежит между ξm1 и xm1(τ,ξ,α). Вследствие того что ψxm1(ξ,α)eq0, будем также иметь ψxm1(x~,α)eq0, если только α достаточно близко к α; отсюда и получится нужное нам уравнение φm1(ξ,α)=0. Этим самым доказано существование периодических решений с периодом τ в окрестности α в предположении, что |A|eq0.

Будем считать теперь τ параметром и положим α=α. Если m1 условий периодичности
φk(τ,ξ,α)=xk(τ,ξ,α)ξk=0(keqm1)

выполнены, то так же, как и выше, выполняется остающееся условие φm1(τ,ξ,α)=0. В качестве функционального определителя для τ=τ,ξ=ξ, очевидно, получим опять |A|. Следовательно, в предположении |A|eq0 существуют такжє периодические решения с фиксированным параметром α и с любым заданным периодом τ, достаточно близким к τ; начальные значения этого решения можно разложить по степеням ττ. Для некоторых исследований выгодно ввести вместо τ значение интеграла ψ(x,α)=γ на рассматриваемой замкнутой траектории как новую переменную. Пусть γ=γ для исходного решения x=x(t,ξ,α). Тогда к m1 уравнениям (25) прибавляется еще следующее:
ψ(ξ,α)γ=0.

Это будут m уравнений с m неизвестными ξ1,,ξm1,τ. Функциональная матрица этой системы при ξ=ξ,τ=τ получается из матрицы (m+1)-го порядка
D=Cfψx0

вычеркиванием m-го столбца и ( m1 )-й строки. Вследствие условий (17),(21) и (23) ( m1)-я строка в D не зависит от остальных; то же самое имеет место для m-го столбца. Следовательно, наш определитель обязательно отличен от нуля, если матрица D имеет ранг m. При этом предположении вышеприведенную систему уравнений можно разрешить в окрестности γ=γ разложением в ряды разностей ξkξk (k=1,,m1) и ττ по степеням γγ. Выполнение условия φm1(τ,ξ,α)=0 и периодичность обусловлены существованием интеграла.

Наконец, если α будет переменной, а γ=γ фиксировано, при тех же предположениях получаем существование соответствующих разложений по степеням αα. Тогда, следовательно, для каждого α вблизи α существует в окрестности исходного решения периодическое решение с одним и тем же значением постоянной интеграла γ=γ.

Для действительного определения тех разложений в степенные ряды, о которых шла речь, необходимо, разумеется, знать полное решение x(t,ξ) в окрестности ξ=ξ,t=τ; в то время как для того, чтобы проверить, отличен ли соответствующий функциональный определитель от нуля, требуется только интегрирование линейной системы (7) с использованием уже известного исходного периодического решения.

Если будет известно большее число интегралов, не зависящих от t, то метод можно соответствующим образом изменить, но не очень сильно.

Применим теперь метод малого параметра к ограниченной задаче трех тел. Пусть точки P1,P2,P3 имеют массы m1=μ,m2=1μ, m3=0, где 0<μ<1 и пусть P1,P2 вращаются с угловой скоростью, равной единице, вокруг общего центра инерции. Введем, как в §17, вращающуюся систему координат, относительно которой P1,P2 и P3 будут иметь координаты (1μ,0),(μ,0) и (x,y). Если положить x1=x,x2=y,x3=x˙ и x4=y˙, то система (17;3) даст следующие уравнения движения точки P3 :
x˙1=x3,x˙2=x4,x˙3=2x4+x1+Fx1,x˙4=2x3+x2+Fx2,

где
F=(1μ)[(x1+μ)2+x22]1/2+μ[(x1+μ1)2+x22]1/2.

Эта система имеет ту же форму, что и система (1), с параметром α=μ и с m=4. При μ=0 масса точки P1 обращается в нуль и формулы
x1=rc,x2=rs,x3=rωs,x4=rωc,c=cos(ωt),s=sin(ωt)

дают тогда при действительной постоянной ωeq0 периодическое решение, если r3(ω+1)2=1. Период этого решения τ=2π|ω|1. Это решение выберем за исходное, полагая α=0. При этом предположим, что req1; в противном случае точка P3 должна была бы пройти через место расположения (1,0) точки P1, а это невозможно, т. к. точка x1=1μ,x2=0 будет особой точкой системы (28) при μeq0, и эта точка стремится к P1 при μ0. Точно так же следует предположить, что ωeq2,1,0. Нам надо установить, существуют ли периодические решения системы (28) при достаточно малых положительных значениях μ.

Если через fk(x,μ)(k=1,,4) обозначить правые части системы (28) и через ξk обозначить начальные значения ξ1=r,ξ2=ξ3=0, ξ4=rω исходного решения при t=0, то f3(ξ,0)=rω2eq0. Следовательно, вместо fmeq0 будем иметь f3eq0. Для применения метода малого параметра необходимо найти решение уравнений в вариациях (7). В нашем случае они интегрируются в элементарных функциях, причем для интегрирования нужно сделать подстановку
y2k1=x2k1c+x2ks,y2k=x2k1s+x2kc(k=1,2).

Образуем теперь матрицу C=XE, и отсюда в соответствии с равенством (20) матрицу B, но при этом заменим на f третий столбец матрицы C вместо последнего столбца. Вычисление показывает, что |B|=0, поэтому первый метод применять нельзя. Причина этому существование для ограниченной задачи трех тел так называемого интеграла Якоби
ψ(x,μ)=12(x32+x42x12x22)F.

Далее, ψx4(ξ,0)=rωeq0, таким образом мы вместо ψxm1eq0 имеем ψx4eq0. Если образовать квадратную матрицу третьего порядка A вычеркиванием третьего столбца и четвертой строки матрицы C, то вычисление дает
|A|=24πsin2πω.

Чтобы этот определитель не был равен нулю, нужно потребовать, кроме ωeq2,1,0, также выполнения условия
ωeqg1(g=±1,±2,).

Для радиуса r исходного решения соответственно получаются значения (g1+1)2/3 с точкой накопления 1 , которая исключается из рассмотрения. Следовательно, при этих предположениях для достаточно малых положительных μ существуют периодические решения системы (28) с периодом τ=τ=2π|ω|1.

Далее, пусть μ=μ будет достаточно малым положительным числом, для которого имеется периодичєское решение системы (28) с периодом τ=τ=2π|ω|1 и пусть γ есть соответствующее значение постоянной интеграла Якоби (31). Покажем, используя метод Пуанкаpe, что для каждого достаточно близкого к γ значения γ будут существовать периодические решения с периодом, близким к τ. Для этого исследуем теперь ранг квадратной матрицы пятого порядка D, определенной равенством (27). Если вычеркнуть в D третий столбец и четвертую строку, то соответствующий минор при ξ=ξ,τ=τ,μ=0 имеет значение 4r2ω3sin2πω, которое не равно нулю, если выполняется условие (32). Если теперь при фиксированном ω выбрать положительное число μ достаточно малым, то ранг матрицы D при μ=μ на рассматриваемом периодическом решении равен 4 . Следовательно, для каждого такого значения μ существует семейство периодических решений, зависящее от γ; период τ этих решений может быть разложен в окрестности γ по степеням γγ, и при γ=γ он имеет исходное значение τ. Таким образом, исходя из μ=0 и кругового решения с периодом τ=2π|ω|1(ωeq0,2,g1), нам удается найти при малых положительных μ сначала периодические решения системы (28) с тем же самым периодом и затем после фиксирования μ удается найти семейство периодических решений, зависящих от параметра γ, причем период этих решений τ, вообще говоря, не равен τ.

Метод малого параметра дает периодические решения только для достаточно малой окрестности значений γ. Интересно было бы изучить поведение решений при аналитическом продолжении по γ. Рассмотрим только случай системы Гамильтона
x˙k=Eyk,y˙k=Exk(k=1,,n),

для которой m=2n; тогда ψ=E(x,y) будет интегралом этой системы. В частности, при n=2 и y1=x3x2,y2=x4+x1 получается система (28), для которой ψ определяется выражением (31). Пусть теперь G есть область действительного пространства (x,y), в которой функция Гамильтона E регулярна и не имеет стационарных точек. Будем исходить из периодического решения C, лежащего в G, которое к тому же не является равновесным. На этом решении интеграл ψ=E также не будет нигде стационарным. Пусть ранг соответствующей матрицы D, определенной равенством (27), есть m. Если обозначить через γ значение параметра E=γ для заданного решения, то метод малого параметра дает семейство периодических решений Cγ, зависящее от параметра γ, начальные значения которых x=ξ,y=η и период τ можно разложить по степеням γγ в окрестности γ. При этом Cγ будет исходной кривой семейства C, и для достаточно малых по абсолютной величине γγ кривая Cγ лежит целиком в G. Повторно применяя метод малого параметра, продолжим аналитически это решение вдоль действительной оси γ. При этом предположим, что все решения Cγ продолжены на весь интервал γγ<γ0, и все они лежат в G. Исследуем поведение Cγ при γγ0. Если для каждой замкнутой ограниченной подобласти H области G существует число ε>0, такое, что при условии γ0ε<γ<γ0 никакая Cγ уже не лежит целиком в H, то мы будем говорить, что Cγ покидает G при γγ0. Пусть рассматривается не этот случай. Тогда по теореме о накоплении можно найти такое H и такую последовательность γγ0, что все Cγ будут целиком лежать в H и соответствующие начальные значения ξ,η будут сходиться к точке ξ0,η0 области H. При этом возможен случай, когда для каждой такой последовательности период τ=τγ, соответствующий кривой Cγ, стремится к . Этот случай мы также исключаем из рассмотрения. Тогда можно найти такую подпоследовательность, для которой τγ стремится к конечному предельному значению τγ0. Величина τγ0 может быть нулем, так как в противном случае вследствие непрерывной зависимости решений от начальных данных точка ξ0,η0 соответствовала бы равновесному решению системы (33), в то время как было предположено, что из области G исключены стационарные точки интеграла E. Тогда вследствие тех же самых теорем непрерывности решения Cγ, принадлежащие названной последовательности, стремятся к решению Cγ0 с начальными значениями ξ0,η0 и периодом τγ0, и это решение, во всяком случае, лежит в H, а следовательно, и в G.

Если ранг матрицы D для решения Cγ0 равен опять m, то решение Cγ можно, очевидно, продолжить за γ0. Остается рассмотреть случай, когда ранг меньше m. Поэтому нужно исследовать при старых обозначениях m аналитических уравнений (25) и (26) вблизи γ=γ0,τ,ξk (keqm1), для которых функциональный определитель равен нулю тождественно относительно τ и ξk, но в то же самое время при γγ0 имеется однопараметрическое семейство действительных решений. Тогда путем использования леммы Вейерштрасса можно показать, что существует решение в виде рядов по степеням (γ0γ)1/p с действительными коэффициентами, где p есть выбранное подходящим образом наименьшее натуральное число. Следовательно, в этом случае имеется точка ветвления порядка p1, и решение можно продолжить и при γ=γ0. Если p нечетное, то для γ>γ0 получим опять действительные значения рядов для τ,ξk. Напротив, если p четное, то корень (γ0 γ)1/p имеет при γ<γ0 два различных действительных значения. Следовательно, если в последнем случае γ устремить к γ0 по другой действительной ветви (γ0γ)1/p, то получится второе семейство периодических решений, которое отлично от первоначального. Поэтому в последнем случае также можно построить аналитическое продолжение решений.

Аналогично можно рассмотреть случай, когда γ убывает от значения γ. С другой стороны, можно опять начать процесс продолжения при γ0 и получить в данном случае другие точки ветвления γ1,γ2, при действительном продолжении решения в интервале от γk1 до γk (k=1,2,). Такое продолжение будет возможным, если не встретится какой-нибудь исключенный из рассмотрения случай, т. е. либо если Cγ покинет область G, либо если τγ будет неограниченно возрастать.

Для ограниченной задачи трех тел в качестве G можно выбрать пространство всех действительных x1,x2,y1 и y2, из которого особые точки x1=μ и x2=0, а также пять стационарных точек функции E выброшены. Если траектория Cγ при γγ1 покидает область G, то это означает, что выброшенные точки являются точками накопления точек Cγ. Для стационарных точек функции E предельным переходом из Cγ получаем равновесные решения. Для особых точек известно соответствующее регуляризирующее преобразование, которое дает в пределе траектории столкновения и показывает, что и здесь можно построить аналитическое продолжение по γ. Процесс продолжения периодических решений ограниченной задачи трех тел Стремгреном и его сотрудниками был осуществлен численно. Встречающиеся при этом теоретические вопросы подробно разработаны Винтнером [2].

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