Электронный журнал «Труды МАИ». Выпуск № 45
www. ша1. гц/ваепсеАгМу/
УДК: 521; 521.1; 521.1:523.03; 521.2; 521.3; 531.3
Исследование резонансных колебаний в фотогравитационной
задаче Ситникова
В. О . Калас Аннотация
Рассматривается задача о движении пассивно гравитирующей точки в поле притяжения двух одинаковых массивных тел, излучающих световую энергию. Считается, что масса гравитирующей точки много меньше массы притягивающих тел, поэтому ее влиянием на движение основных тел можно пренебречь. Траектории движения основных тел относительно общего центра масс - кеплеровские эллипсы.
В окрестности устойчивого тривиального равновесия исследуются одномерные колебания точки вдоль оси Z, перпендикулярной плоскости орбиты главных тел и проходящей через центр масс системы. Получены уравнения колебаний с точностью до членов четвертого порядка малости по отклонениям. Получены условия существования параметрического резонанса 1:2. Методом усреднения исследованы резонансные колебания точки как при строгом резонансе, так и в окрестности резонанса: выведены усредненные уравнения, показано, что они допускают первый интеграл, построен фазовый портрет колебаний в окрестности резонанса (и при строгом резонансе). Результаты исследований дублируются в разных системах координат.
Примером описываемой задачи является система из двух одинаковых по массе и излучению звезд, между которыми сосредоточены облака из пылевых частиц, подверженных влиянию как сил светового давления (парусный эффект), так и сил гравитации.
Ключевые слова:фотогравитационная задача ситникова; равновесие; параметрический резонанс; окрестность резонанса; фазовый портрет; бифуркационное значение параметра.
1. Общая постановка задачи и основные определения
Рассмотрим задачу о движении пассивно гравитирующей точки С, находящейся в поле притяжения двух массивных тел одинаковой массы щ = т2, излучающих световую энергию. В этом случае точка С испытывает, помимо гравитационного притяжения, световое давление со стороны тел т1, т2. Предполагается, что точкаС массы т движется вдоль оси Ог, проходящей через центр масс притягивающих тел, перпендикулярно плоскости их движения. Одномерное движение вдоль оси Ог возможно в силу симметрии задачи (см. рис. 1). Считается, что масса гравитирующей точки много меньше массы притягивающих тел, т.е. т«щ, поэтому ее влиянием на движение основных тел можно пренебречь (ограниченная постановка задачи). Траектории движения основных тел относительно общего центра масс - кеплеровские эллипсы эксцентриситета е .
В классическом случае, когда световое давление равно нулю, задачу о движении пассивно гравитирующей точки Свдоль оси Ъ в поле притяжения двух массивных тел одинаковой массы щ = щ принято называть задачей Ситникова.
Известно, что задача Ситникова связана с проблемой классификации финальных движений в задаче трех тел. Полную классификацию различных типов финальных движений дал Шази [1], постулируя наличие осциллирующих решений, которым отвечают неограниченные колебания координаты г, при условии, что 2 не стремится к бесконечности со временем. В 1954 году А.Н. Колмогоров предложил изучить частный случай задачи трех тел на предмет исследования топологии некоторых подмножеств фазового пространства, порождающих разные типы финальных движений. К.А. Ситников [2] доказал для этого случая существование осциллирующих движений, для которых координата 2 испытывает
гги
Рис. 1.1. Система трех тел
бесконечное число выбросов на произвольно большие расстояния, однако всегда возвращается в начало координат. Позднее, В.М. Алексеев [3] исследовал хаотические движения в задаче Ситникова методами символьной динамики и показал, что при определенных условиях в этой задаче реализуются все возможные комбинации финальных движений по Шази, доказал также существование осциллирующих и гиперболических решений, меняющих свой финальный тип за счет явления «полного захвата». MozeгJ. познакомил западных исследователей с задачей Ситникова, читая лекции по небесной механике в Принстонском университете в начале 70-ых и публикуя результаты своих исследований этой задачи в виде монографии [4]. Впоследствии появилось множество работ, посвященных задаче Ситникова, в основном зарубежных авторов. Например, в работах [5-9] исследовались регулярные (периодические) орбиты, работы [10-13] посвящены вопросам хаотической динамики.
Если световое давление отлично от нуля, задачу о движении точкиС будем называть фотогравитационной задачей Ситникова.
Получим уравнения движения точки массы т вдоль оси Ог. Пусть ^, ^ - силы,
действующие со стороны масс т, т соответственно. Эти силы будут иметь следующий вид:
дтщ дтт2
1 п2 ' Г2 / „2 , К\ К2
где = Я2 = V г2 + г2 . Здесь ^ - коэффициент редукции массы частицы. Он меньше
единицы, так как обусловлен воздействием светового потока на частицуС, а это воздействие направлено в сторону, прямо противоположную силе притяжения (см. [14]).
Тогда уравнение движения точкиС будет таким та = ^ + Г2
Проектируем это уравнение на ось Ог и, учитывая, что
л I I 1 1 1 2•> 22 + г2
имеем уравнение второго порядка
г = --
V 3/2 •
('2 + г' )3
Здесь л = /(т + т) , 2г - расстояние между притягивающими телами, зависящее от эксцентрической аномалии Е по формуле
2r = a(1 — e cos E) (1.1)
Параметр a - большая полуось орбиты. Аномалия E зависит от времени t в силу уравнения Кеплера
E — e sin E = n(t — т),
1/2 —3/2
где n = /л a - среднее движение притягивающих тел, т - один из кеплеровских элементов орбиты (1.1), представляющий собой время прохождения перицентра.
Будем считать, без ограничения общности, что f = 1, щ + щ = 1, a = 1. Этого можно добиться с помощью выбора единиц измерения времени, расстояния и массы притягивающих тел. В этом случае / = 1, n = 1, поэтому период обращения основных тел будет равен 2ж . Уравнения движения примут следующий вид:
q •z z-i z+--ш = 0
( z 2 + r 2(t) ) (1.2)
Очевидно, что z = 0 - частное решение этого уравнения.
Ниже приводится известная в небесной механике формула решения E(t) уравнения Кеплера в виде ряда по e, сходящегося при e < 0.662...: » 1 dn—1 (sinn g)
E(e,g) = д + У---—г—^e", g = (t — т) (1.3)
V W Ь t n! dgn—1 Ь '
2. Нахождение параметрического резонанса для системы первого приближения
Разложим второе слагаемое уравнения (1.2) в ряд Тейлора относительно z
z + ±z — ^z3 + O( z5) = 0 r3 2r
С точностью до членов четвертого порядка относительно z имеем
q г 3q -з
- 1 —-z--
r3 2r
z +X z — z3 = 0 (2.1)
Используя уравнение (1.1), разложим сомножители -1,-1 в ряд относительно
3 ' 5
r r
эксцентриситета e
1 =1 (8 + (24 cos E) e + ( 48 cos2 E) e2 + (80 cos3 E) e3 +...) r5 a
— = — (32 + (160 cos E) e + (480 cos2 E) e2 +(1120 cos3 E) e3 +...)
Используя разложение в ряд (1.3), получим следующие выражения для 1 1
3 ' 5
r r
1 = \ (8 + 24e cos g + e2 (12 + 36 cos 2g) + e3 (27 cos g + 53 cos 3g)) (2.2)
\ =1 (32 + 160e cos g + e2 (160 + 320 cos 2g) + e3 (540 cos g + 580 cos 3g)) (2.3)
Учитывая, что параметр a=1 и предполагая, не ограничивая общности, что г = 0, запишем уравнение (2.1) в следующем виде
if + q -a2 z = q • g(t) z + q • h(t )z3 (2.4)
где
a = ^(8 + 12e2),
g (t) = - (24e cos t + e2 (36 cos 2t) + e3 (27 cos t + 53 cos 3t)),
h(t) = (48 + 240e cos t + e2 (240 + 480 cos 2t) + e3 (810 cos t + 870 cos 3t))
Оставляя члены до второго порядка относительно z и относительно e, из (2.4) получим следующее уравнение
z + 8q(1 + 3e cos t) z = 0 (2.5)
Положим w2 = 8q . Тогда уравнение (2.5) примет вид z + w2(1 + 3e cos t) z = 0 (2.6)
Рассмотрим задачу о нахождении резонанса в уравнении (2.6). Здесь w - собственная частота колебаний, а параметрическая частота колебаний равна единице (сомножитель перед переменной t). С помощью стандартной замены переменных z = х cos y, z = -xw sin y
это уравнение сводится к одночастотной неавтономной системе. Вводим новую быструю переменную y2 = t, исключая время t из правой части. Тогда уравнения колебаний предстанут в виде двухчастотной автономной системы х = 3ewx sin y cos y cos y
y = w + 3ew cos2 y cos y (2.7)
& = 1
Вычислим временное среднее от правой части уравнения для медленной переменной x вдоль быстрых решений y = wt + y0, y = t невозмущенной задачи
X(x, y10) = lim — j*3wx sin(wt + y10)cos(wt + y10)cos t dt
TT 0
Проводя интегрирование, получим
X (x, У10) = Jim
3wx Г 4w cos 2 y10 + (1 - 2w)cos ( 2(wT + y10) + T
t 4t
- lim
t 4t
, (2w + 1)(2w -1)
4 ' (2 8) 3wx ( (1 + 2w) cos (2( wT + y10) - T ) ^ '
(2w + 1)(2w -1)
Учитывая положительность w, можно сделать вывод, что при w = 1 присутствует
неопределенность типа 0/0, которая ведет к разрыву среднего в окрестности резонанса. Этот резонанс представляет собой главный параметрический резонанс. Других резонансов, проявляющих себя в первом приближении метода усреднения, нет. Найдем параметр q, соответствующий данному резонансу.
2 1
W 1
q = — = — . 8 32
Итак, мы исследовали уравнение на наличие резонансов в первом приближении и нашли коэффициент редукции массы частицы q, при котором присутствует резонанс.
3. Нахождение параметрического резонанса для системы в нелинейном
приближении
Рассмотрим уравнение (2.4): z + q • со2z = q • g(t)z + q • h(t)z3.
Предположим, что величина z является малой величиной порядка s. Делая замену переменных z = sy, z = sy,
гдеу~1, после чего сокращая левую и правую часть уравнения на s, получим y + qC y = qg (t) y + qh(t )s2 y3
y + qC y = q(-24e cos t - 36e2 cos 2t) y + q(48 + 240e cos t )s2 y3
Положим s2 = e. С точностью до членов второго порядка малости по ебудем иметь: y + 8qy = e ((-24q cos t)y + 48qy3) .
Пусть w2 = 8q . Тогда уравнение будет иметь следующий вид:
y + w2 y = e ( (-24q cos t) y + 48qy3) (3.1)
Исследуем параметрические колебания этого уравнения. Преобразуем (3.1), используя замену переменных Ван-дер-Поля: y = x cos wt - x sin wt
. 2 (3.2)
y = XjW Sin wt - x2w cos wt
В новых переменных x1, x2 уравнение (3.1) примет вид стандартной по Боголюбову системы
e sin wt,
x =-"
- ((-24q cos t)(x cos wt - x2 sin wt) + 48q(x cos wt - x2 sin wt)3)
w (3.3)
e cos wt 3
x2 =--((-24q cos t)(x cos wt - x2 sin wt) + 48q(x cos wt - x2 sin wt) )
w v '
Эту систему можно рассматривать как систему с двумя быстрыми переменными t и
wt. Таким образом, уравнение (3.1) эквивалентно следующей системе уравнений
x =--y1 ((-24q cos y2)(x cos y - x2 sin y ) + 48q(x cos y - x2 sin y )3)
w v '
e cos y í \
x =--— ((-24q cos y2)(x cos y - x2 sin y ) + 48q(x cos y - x2 sin y )3) (3 .4)
yi = w y2 = 7
при условии, что начальные значения быстрых переменных y1, y2 равны нулю y1 (0) = 0, y2(0) = 0 . Здесь 7 = 1. Вектор частот будет таким: w = (w,7) = (w,1). Система (3.4)
содержит две медленные и две быстрые переменные.
Вычисляя временное среднее от правой части уравнений (3.3), получим следующую систему:
18qv2(4w2 -1) 2 2ч V = e * 2'-rr(v2 + v22)
w(2w - 1)(2w +1) 18qy (4w2 -1) w(2w - 1)(2w +1)
(3.5)
V = -e '8^ -1> (v2 + v2)
Отсюда следует, что при w = 0, ± 1 числители и знаменатели правых частей системы
(3.5) обращаются в ноль, хотя прохождение параметра через значения ± 1 не ведет к разрыву временного среднего, поскольку 4w2 -1 = (2w -l)(2w + 1) . Тем не менее, частоты
<
w = ± — являются резонансными, так как вид усредненных уравнений резко меняется, когда
н принимает указанные значения. Отбрасывая нулевые и отрицательные значения частоты н, получим резонансные значения частоты собственных колебаний: 1
w = —. 2
Отсюда д = -1. Заметим, что коэффициент ц получился таким же, как и в случае линейного приближения.
4. Исследование системы в нелинейном приближении при строгом резонансе и в
окрестности резонанса
Исследуем колебания системы при резонансе w «1. В этом случае разность частот
А = w -1 у есть величина порядка е, если 8 «1: А = ед.
Вводим новую медленную переменную 6 = у -1 у . Выражая отсюда у = 6 +1 у, получим новый вид уравнений (3.4):
х, =-е- 1
| 6 +1 у I х
(-24д с08 у2 ) Х1 С0Б [6 + ^ у2 У Х2 ^ |6 + ^ у2
2
- е
1
1
1 48дI х1с08[6+-У2у Х2 [6+-у2
2 + А
С0§|б+ 2 У2 - | | 1 I I 1 ^
х^ = -е-Ч--(-24дсобу ) х СОБI 6 +—у I-х2 Б1ПI 6 + — у2
1 2
+ А
V
-е
С0Б [6 + 1 У2 ^ А
+ А
48д х соб |б +1 у I - х бш 1б +1 у2
6 = ед у2 = 1
(4.1)
3
<
3
1
Система уравнений (4.1) замкнута относительно переменных х, X,0, У2 и не имеет резонансов, поскольку вектор частот и = (у) = (1) не является резонансным. Она содержит три медленных и одну быструю переменные. Это значит, что мы можем усреднить ее как в нерезонансном случае с обычной гарантией точности приближения.
Усредняем уравнения (4.1) по быстрой переменной у2 на периоде 2ж :
= \2ед ¡зи2 + 3 _„ + 2м, б1П0СОБ0 + 2щ(СОБ0)2)
1 1 + 2Д1 12221 ^ ;
12ед /о • з ? \
щ =-(2м, СОБ 0 _ 2щ Б1п 0 СОБ 0 _ м, _ Зщ _ Зим,, )
2 1 + 2Д1 1 2 1 1 1 '
(4.2)
где 0 = еЛ. Таким образом, колебания в окрестности резонанса и«1 описываются
уравнениями (4.2), содержащими медленно меняющийся со временем параметр 0 .
Исследуем поведение решений системы уравнений (4.2), когда частота собственных
колебаний и = 1 строго. Тогда величина 8 равна нулю и 0 = еЛ = 0. Усредненные
уравнения (4.2) примут более простой вид: щ = 12ед (Зм2 щ + Зи\ + и )
(4.3)
М2 = _12ед(Змщ\ + ЗиЗ _щ)
Для простоты введем обозначение ед = £ . Тогда система (4.3) примет вид: щ = 12^ (Зм2 щ + ЗиЗ + щ)
/ 2 3 \ (44)
| щ = _12£(Змщ2 + Зм1 _ щ)
Легко видеть, что уравнения (4.4) можно записать в виде уравнений Гамильтона
дИ . дИ
д = е — , р = _£— (4.5)
др дд
если переменную и1 считать обобщенной координатой д, а и2 - обобщенным импульсом р.
Заметим, что в уравнениях Гамильтона переменная д выступает в качестве обобщенной координаты и не имеет ничего общего с коэффициентом редукции массы частицы.
Функция Гамильтона Н(д, р)будет такой: И = 9(р2 + д2)2 + 6(р2 _д2)
Учитывая, что функция Гамильтона не зависит явно от времени I, приходим к выводу, что уравнения (4.5) допускают первый интеграл
<
H (q, p) = const
В исходных переменных этот интеграл имеет форму
3(u2 + u2)2 + 2(u2 — u2) = C
(4.6)
Для получения параметрического представления фазовых кривых, описываемых интегралом (4.6), положим щ = ±>/у, и2 = .
Подставим эти выражения в уравнение (4.6). Получим квадратное уравнение относительно V.
2 ) 2 2 1
2 w--V + w + — w — С = О
3) 3 3
Решение уравнения будет таким:
w —1 |± — 12w + 3C
На рисунке 4.1 представлен фазовый портрет резонансных колебаний, рассчитанный на основе параметрического представления кривых.
Рис.4.1. Фазовый портрет резонансных колебаний Сепаратриса колебаний выделена жирным черным цветом, она проходит через неустойчивое положение равновесия щ = и2 = 0, отвечающее константе С=0. В окрестности
■4
стационарных значений ^ =±/Щ = 0, соответствующих константе С = -1,
изображающая точка описывает замкнутые кривые. Амплитуда колебаний растет с приближением константыС к нулю, поэтому замкнутые кривые заполняют всю область, ограниченную сепаратрисой. Дальнейшее возрастание константыС, когда постоянная
v
V1,2 = —
интегрирования принимает положительные значения, ведет к появлению замкнутых периодических решений, охватывающих начало координат, сепаратрису и движения вблизи
равновесия их = и2 = 0 •
Заметим, что положение равновесия их = и2 = 0 является устойчивым при s = 0.
Однако при сколь угодно малом s, отличным от нуля, это равновесие неустойчиво, что можно легко проверить, используя теорему об устойчивости по первому приближению. Это значит, что в системе наблюдается явление параметрического резонанса, то есть нарастание амплитуды колебаний при малых возмущениях в начальных условиях, выводящих систему из равновесия. Амплитуда таких колебаний - величина порядка единицы, независимо от величины малого параметра s .
Теперь проведем исследования в окрестности резонанса w = 1. Уравнения (4.2)
допускают неавтономный первый интеграл. Но пользоваться им неудобно, поэтому
исследование окрестности резонанса w = 1 будем проводить в переменных амплитуда (а ) -
фаза (Р). В этом случае фазовые переменные уравнения (3.1) имеют следующее
представление:
y = acosР, y = -awsinР
Для простоты в уравнении (3.1) положим s = eq . Тогда оно примет вид:
y + w2y = s((-24cost)y + 48y3) (4.7)
Уравнение движения в новых переменных будет таким: s
(4.8)
á =--((-24 cos t )á cos P + 48a3 cos3 p) sin p
wv '
P = w - — ((-24 cos t) cos P + 48á2 cos3 p) cos p w
Вводя новую быструю переменную y1 = t, получаем многочастотную систему á = - — ((-24 cos y )á cos P + 48á3 cos3 p) sin P
P = w -—((-24 cos y) cos p + 48á2 cos3 p) cos p (4.9)
— w
y =1
<
Для преобразования этой системы в окрестности резонанса и = 1 введем медленную фазу щ = Р_ 1 у . После подстановки быстрой переменной /3 = щ +1 у в систему (4.9)
получим
(
а =--
w
1
1
(-24cosy)acosу+2УJ + 48а cos у+2У
f f
sin \V + - У1
у/ = A-— (-24cosy)cos+1 y j + 48а2 cos3 +1 y j cos+1 y
(4.10)
У = 1
После последующего усреднения по времени вдоль быстрых движений у = I. получим усредненные уравнения:
v = 6—v sin2v2 w
— 9
v2 =A + 6—(cos2v2 -3v ) w
(411)
Здесь v,, v, имеют смысл переменных а, у соответственно, A = w -1 - резонансная
расстройка, порядок малости которой равен s : A = sS.
Перейдем к полярным координатам u = v cos v2, v = v sin v2
Усредненные уравнения (4.11) примут вид:
. 6s II = — v w
3 (u2 + v2) +1 A v ' 6s
6s
V =--u
w
3 (u2 + v2)-1 A v ' 6s
(4.12)
Несложно видеть, что эти уравнения допускают автономный первый интеграл
u
1 Aw j 2 (л Aw j 3 / 2 2\2 ^ 1 +- I-v 1--I--(u + v I = C ,
V 6s J V 6s J 2V ' 1
(4.13)
который гораздо проще неавтономного интеграла усредненной задачи, представленного в
8и
переменных щ, м2. Учитывая, что А = 8 и полагая а = , представим интеграл (4.13) в
следующем виде:
s
<
<
3 (и2 + V2) + 2 [у2 — и2)-а (и2 + V2 ) = С
2
(4.14)
Варьируя значения параметра а, можно проследить за качественным изменением фазового портрета в зависимости от выбора данного параметра. Ниже представлены фазовые портреты для различных областей параметра а. 1) а <-2.
В этом случае коэффициенты (2 — а) > 0, (2 + а) < 0 при членах V2, и2 соответственно. На рис. 4.2 представлен фазовый портрет для этого случая, а = —4 .
Как видно из графика, здесь есть одна устойчивая особая точка в начале координат. 2) —2 < а < 2.
В этом случае (2 — а) > 0, (2 + а) > 0. Ниже представлен фазовый портрет при а = 1. Сепаратриса проходит через неустойчивое положение равновесия и = V = 0, отвечающее константе С = 0. Также есть две устойчивые особые точки.
Рис.4.2. Фазовый портрет при а < —2
Рис.4.3. Фазовый портрет при —2 < а < 2
3) а > 2.
Коэффициенты (2 — а) < 0, (2 + а) > 0.
Рис.4.4. Фазовый портрет при а > 2 Рисунок 4.4 соответствует значению а = 10. В этом случае имеем три устойчивых положения равновесия и два неустойчивых, через которые проходит сепаратриса.
Итак, мы имеем три качественно различных фазовых портрета, соответствующих интервалам а <—2, —2< а < 2 и а > 2. Значения параметра а = ±2 являются
бифуркационными, при которых наблюдается изменение топологического типа фазового портрета колебаний.
5. Выводы
Основным результатом работы является описание параметрических резонансных колебаний в нелинейном приближении. Получены условия параметрического резонанса 1:2 для системы первого приближения и для системы в нелинейном приближении. Методом усреднения исследованы резонансные колебания точки как при строгом резонансе, так и в окрестности резонанса: выведены усредненные уравнения, показано, что они допускают первый интеграл, построен фазовый портрет колебаний. Найдено бифуркационное значение параметра a (a = 2), при котором наблюдается изменение топологического типа фазового портрета колебаний. Результаты исследований дублируются в разных системах координат.
Примером описываемой задачи является система из двух одинаковых по массе и излучению звезд, между которыми сосредоточены облака из пылевых частиц, подверженных влиянию как сил светового давления (парусный эффект), так и сил гравитации.
Библиографический список
1. J. Chazy. Surl'allurefinaldumouvementdansleproblemedestrioscorpsquant le temps croitindéfiniment. Annales de l'Ecole Norm. Sup., 3 ser., 39,1922, P. 29-130
2. К.А. Ситников Существование осциллирующих движений в задаче трех тел. ДАН, 1960, Т. 133, № 2, С. 303-306
3. В.М. Алексеев Квазислучайные динамические системыI, II, III. Матем. сб.,1968, 76:1, с. 72-134; 1968, 77:4, с. 545-601; 1969, 78:1, с. 3-50
4. J. Moser, Stable and Random Motions in Dynamical Systems, Princeton University Press, Princeton, NJ, 1973
5. E. Belbruno, J. Llibre, and M. Olle, On the families of periodic orbits which bifurcate from the circular Sitnikov motions, Celestial Mech. and Dynam. Astronom., 1994, 60 , p. 99-129.
6. M. Corbera and J. Llibre, Periodic orbits of the Sitnikov problem via a Poincare map, Celestial Mech. Dynam. Astronom., 2000, 77 , р. 273-303
7. L. Jiménez-Lara, A. Escalona-Buendía, Symmetries and bifurcations in the Sitnikov problem, Celestial Mech. Dynam. Astronom., 2001, 79 97-117
8. J. Llibre, R. Ortega, On the families of periodic orbits of the Sitnikov problem, SIAM J. Applied Dynamical Systems, 2008, 7, p. 561-576
9. J Hagel A new analytic approach to the Sitnikovproblem, Celes. Mech., 1992,53, p. 267-292.
10. T. Kovacsand B. Erdi The structure of the extended phase space of the Sitnikov problem Astron. Nachr., 2007, AN 328, No. 8, р. 801 - 804
11. Jie Liu, Yi-Sui Sun On the Sitnikov problem, Celes. Mech., 1990,49, р. 285-302.
12. M. A. Jalali S. H. PourtakdoustRegular and Chaotic Solutions of the Sitnikov Problem near the 3/2 CommensurabilityCelestial Mechanics and Dynamical Astronomy, 1997, 68, 151-162
13. Clark Robinson Uniform subharmonic orbits for Sitnikov problemDiscrete and continuous dynamical systems, Series S, 2008 Vol. 1, No 4, p. 647-652
14. Радзиевский В.В. Астрономический журнал, Т. 27, 1950, стр. 250
15. Д.Р. Меркин, Введение в теорию устойчивости движения (М., «Наука», 1976), с. 319
16. В.Н. ТхайПММ, 70, № 5, 813 (2006)
17. А.М. Ляпунов Общая задача об устойчивости движения (М.-Л., ГИТТЛ, 1950), с. 471
Сведения об авторах
Калас Вячеслав Олегович, студент Московского авиационного института (государственного
технического университета);
тел.: 7-906-100-20-30; e-mail: sss@mail.ru