Теорема Ньютона о вращающихся орбитах - Newtons theorem of revolving orbits - Wikipedia

Рисунок 1: Сила притяжения F(р) заставляет голубую планету двигаться по голубому кругу. Зеленая планета движется в три раза быстрее и поэтому требует более сильного центростремительная сила, который обеспечивается добавлением притягивающей силы обратного куба. Красная планета неподвижна; сила F(р) уравновешивается отталкивающей силой обратного куба. А Гифка версия этой анимации найдена здесь.
Рисунок 2: Радиус р зеленых и синих планет одинаковы, но их угловая скорость отличается в разы k. Примеры таких орбит показаны на рисунках 1 и 3–5.

В классическая механика, Теорема Ньютона о вращающихся орбитах определяет тип центральная сила необходимо умножить угловая скорость частицы на фактор k не влияя на его радиальное движение (рисунки 1 и 2). Ньютон применил свою теорему к пониманию общего вращения орбит (апсидальная прецессия, Рис.3), что наблюдается для Луна и планеты. Термин «радиальное движение» означает движение к центру силы или от него, тогда как угловое движение перпендикулярно радиальному движению.

Исаак Ньютон вывел эту теорему в предложениях 43–45 книги I его Philosophiæ Naturalis Principia Mathematica, впервые опубликовано в 1687 году. В предложении 43 он показал, что добавленная сила должна быть центральной силой, величина которой зависит только от расстояния р между частицей и точкой, закрепленной в пространстве (центром). В предложении 44 он вывел формулу для силы, показав, что это сила обратного куба, изменяющаяся как обратный куб р. В предложении 45 Ньютон распространил свою теорему на произвольные центральные силы, предположив, что частица движется по почти круговой орбите.

Как отмечает астрофизик Субраманян Чандрасекар в его комментариях 1995 г. Principia, эта теорема оставалась в значительной степени неизвестной и неразработанной более трех столетий.[1] С 1997 г. теорему изучают Дональд Линден-Белл и сотрудники.[2][3] Его первое точное расширение произошло в 2000 году благодаря работе Магомеда и Вавды.[4]

Исторический контекст

Ретроградное движение Марс если смотреть с Земли.
Рисунок 3: Планеты, вращающиеся вокруг Солнца, движутся по эллиптическим (овальным) орбитам, которые постепенно со временем вращаются (апсидальная прецессия). В эксцентриситет этого эллипса увеличено для визуализации. Большинство орбит в Солнечной системе имеют гораздо меньший эксцентриситет, что делает их почти круглыми. А Гифка версия этой анимации найдена здесь.

Движение астрономических тел систематически изучается на протяжении тысяч лет. Было замечено, что звезды вращаются равномерно, всегда сохраняя одно и то же положение относительно друг друга. Однако было замечено, что другие тела бродить на фоне неподвижных звезд; большинство таких тел были названы планеты после греческого слова «πλανήτοι» (Planētoi) для "странников". Хотя обычно они движутся в одном и том же направлении по небу ( эклиптика ), отдельные планеты иногда на короткое время меняют свое направление, показывая ретроградное движение.[5]

Чтобы описать это движение вперед и назад, Аполлоний Пергский (c. 262 - с. 190 г. до н.э.) разработал концепцию семявыносящие протоки и эпициклы, согласно которому планеты перемещаются по вращающимся кругам, которые сами перемещаются по другим вращающимся кругам, и так далее. Любую орбиту можно описать с помощью достаточного количества разумно выбранных эпициклов, поскольку такой подход соответствует современному преобразование Фурье.[6] Примерно 350 лет спустя Клавдий Птолемей опубликовал свой Альмагест, в котором он разработал эту систему, чтобы соответствовать лучшим астрономическим наблюдениям его эпохи. Для объяснения эпициклов Птолемей принял геоцентрический космология Аристотель, согласно которому планеты были ограничены концентрическими вращающимися сферами. Эта модель вселенная был авторитетным почти 1500 лет.

Современное понимание движения планет возникло в результате совместных усилий астронома. Тихо Браге и физик Иоганн Кеплер в 16 веке. Тихо приписывают чрезвычайно точные измерения движения планет, из которых Кеплер смог получить свои законы движения планет.[7] По этим законам планеты движутся дальше эллипсы (нет эпициклы ) о солнце (не Земля). Второй и третий законы Кеплера делают конкретные количественные предсказания: планеты охватывают равные площади за равное время, и квадрат их орбитальные периоды равна фиксированной константе, умноженной на куб их большая полуось.[8] Последующие наблюдения планетных орбит показали, что длинная ось эллипса (так называемая линия апсид) постепенно вращается со временем; это вращение известно как апсидальная прецессия. В апсиды орбиты - это точки, в которых вращающееся тело наиболее близко или дальше всего от притягивающего центра; для планет, вращающихся вокруг Солнца, апсиды соответствуют перигелию (ближайшему) и афелию (самому дальнему).[9]

С публикацией его Principia примерно восемьдесят лет спустя (1687 г.), Исаак Ньютон предоставил физическую теорию, объясняющую все три закона Кеплера, теорию, основанную на Законы движения Ньютона и его закон всемирного тяготения. В частности, Ньютон предположил, что гравитационная сила между любыми двумя телами была центральная сила F(р), которые варьировались как обратный квадрат расстояния р между ними. Основываясь на своих законах движения, Ньютон показал, что орбита любой частицы, на которую действует одна такая сила, всегда коническая секция, особенно эллипс, если он не уходит в бесконечность. Однако этот вывод верен только при наличии двух тел ( проблема двух тел ); движение трех или более тел, действующих под действием их взаимного тяготения ( ппроблема тела ) оставался нерешенным в течение столетий после Ньютона,[10][11] хотя решения некоторых Особые случаи были обнаружены.[12] Ньютон предположил, что орбиты планет вокруг Солнца в основном эллиптические, потому что гравитация Солнца является доминирующей; к первое приближение, присутствие других планет можно не учитывать. По аналогии, эллиптическая орбита Луна около Земли преобладала сила тяжести Земли; в первом приближении гравитацией Солнца и других тел Солнечной системы можно пренебречь. Однако Ньютон заявил, что постепенная апсидальная прецессия планетных и лунных орбит была вызвана эффектами этих игнорируемых взаимодействий; в частности, он заявил, что прецессия орбиты Луны была вызвана возмущающими эффектами гравитационного взаимодействия с Солнцем.[13]

Теорема Ньютона о вращающихся орбитах была его первой попыткой количественно понять прецессию апсид. Согласно этой теореме, добавление особого типа центральной силы - силы обратного куба - может создать вращающуюся орбиту; угловая скорость умножается на коэффициент k, а радиальное движение не изменилось. Однако эта теорема ограничена определенным типом силы, которая может не иметь значения; несколько возмущающих взаимодействий обратных квадратов (например, на других планетах) кажутся маловероятными в сумме в точности с силой обратного куба. Чтобы применить свою теорему к другим типам сил, Ньютон нашел наилучшее приближение произвольной центральной силы. F(р) к потенциалу обратного куба в пределе почти круговых орбит, то есть эллиптических орбит с низким эксцентриситетом, что действительно верно для большинства орбит в Солнечной системе. Чтобы найти это приближение, Ньютон разработал бесконечный ряд, который можно рассматривать как предшественника Расширение Тейлора.[14] Это приближение позволило Ньютону оценить скорость прецессии для произвольных центральных сил. Ньютон применил это приближение к тестовым моделям силы, вызывающей апсидальную прецессию орбиты Луны. Однако проблема движения Луны пугающе сложна, и Ньютон никогда не публиковал точную гравитационную модель апсидальной прецессии Луны. После более точной модели Clairaut в 1747 г.,[15] аналитические модели движения Луны были разработаны в конце 19 века. холм,[16] Коричневый,[17] и Делоне.[18]

Однако теорема Ньютона является более общей, чем просто объяснение прецессии апсид. Он описывает эффекты добавления силы обратного куба к любой центральной силе. F(р), а не только к силам обратного квадрата, таким как Закон всемирного тяготения Ньютона и Закон Кулона. Теорема Ньютона упрощает орбитальные задачи в классическая механика исключив из рассмотрения силы обратного куба. Радиальные и угловые движения, р(т) и θ1(т), можно вычислить без силы обратного куба; впоследствии его эффект можно рассчитать, умножив угловую скорость частицы

Математическое утверждение

Рисунок 4: Все три планеты совершают одинаковое радиальное движение (голубой круг), но движутся с разными угловыми скоростями. Голубая планета чувствует только обратный квадрат сила и движется по эллипс (k = 1). Зеленая планета движется под углом в три раза быстрее, чем голубая планета (k = 3); он совершает три витка для каждой орбиты голубой планеты. Красная планета иллюстрирует чисто радиальное движение без углового движения (k = 0). Пути, по которым следуют зеленые и синие планеты, показаны на Рисунок 9. А Гифка версия этой анимации найдена здесь.
Рисунок 5: Зеленая планета движется под углом на одну треть быстрее, чем голубая планета (k = 1/3); он совершает один оборот на каждые три синих орбиты. Пути, по которым следуют зеленые и синие планеты, показаны на Рисунок 10. А Гифка версия этой анимации найдена здесь.

Рассмотрим частицу, движущуюся под произвольным центральная сила F1(р), величина которого зависит только от расстояния р между частицей и неподвижным центром. Поскольку движение частицы под действием центральной силы всегда лежит в плоскости, положение частицы можно описать следующим образом: полярные координаты (рθ1), радиус и угол частицы относительно центра силы (Рисунок 1). Обе эти координаты, р(т) и θ1(т), меняются со временем т как частица движется.

Представьте себе вторую частицу той же массы м и с таким же радиальным движением р(т), но с угловой скоростью k раз быстрее, чем у первой частицы. Другими словами, азимутальные углы двух частиц связаны уравнением θ2(т) = k θ1(т). Ньютон показал, что движение второй частицы может быть произведено добавлением центральной силы обратного куба к любой силе. F1(р) действует на первую частицу[19]

куда L1 величина первой частицы угловой момент, который является постоянная движения (сохранено) для центральных сил.

Если k2 больше единицы, F2 − F1 отрицательное число; таким образом, добавленная сила обратного куба равна привлекательный, как это наблюдается на зеленой планете на рисунках 1–4 и 9. Напротив, если k2 меньше единицы, F2F1 положительное число; добавленная сила обратного куба равна отталкивающий, как это наблюдается на зеленой планете на рисунках 5 и 10 и на красной планете на рисунках 4 и 5.

Изменение пути частицы

Добавление такой силы обратного куба также изменяет дорожка за которой следует частица. Траектория частицы игнорирует временные зависимости радиального и углового движений, такие как р(т) и θ1(т); скорее, он связывает переменные радиуса и угла друг с другом. Для этой цели угловая переменная не ограничена и может неограниченно увеличиваться, поскольку частица несколько раз обращается вокруг центральной точки. Например, если частица дважды оборачивается вокруг центральной точки и возвращается в исходное положение, ее конечный угол не совпадает с ее начальным углом; скорее, он увеличился на 2×360° = 720°. Формально угловая переменная определяется как интеграл от угловой скорости

Аналогичное определение справедливо для θ2, угол второй частицы.

Если путь первой частицы описан в виде р = грамм(θ1), путь второй частицы задается функцией р = грамм2/k), поскольку θ2 = k θ1. Например, пусть путь первой частицы будет эллипс

куда А и B константы; тогда путь второй частицы определяется выражением

Орбитальная прецессия

Если k близка, но не равна одной, вторая орбита похожа на первую, но постепенно вращается вокруг центра силы; это известно как орбитальная прецессия (Рисунок 3). Если k больше единицы, орбита прецессирует в том же направлении, что и орбита (рис. 3); если k меньше единицы, орбита прецессирует в обратном направлении.

Хотя может показаться, что орбита на рисунке 3 вращается равномерно, то есть с постоянной угловой скоростью, это верно только для круговых орбит.[2][3] Если орбита вращается с угловой скоростью Ω, угловая скорость второй частицы больше или меньше угловой скорости первой частицы на Ω; другими словами, угловые скорости удовлетворяли бы уравнению ω2 = ω1 + Ом. Однако теорема Ньютона о вращающихся орбитах утверждает, что угловые скорости связаны умножением: ω2 = 1, куда k является константой. Объединение этих двух уравнений показывает, что угловая скорость прецессии равна Ω = (k − 1)ω1. Следовательно, Ω постоянно, только если ω1 постоянно. Согласно сохранению углового момента, ω1 изменяется с радиусом р

куда м и L1 первые частицы масса и угловой момент соответственно, оба из которых постоянны. Следовательно, ω1 постоянно, только если радиус р постоянна, т.е. когда орбита представляет собой круг. Однако в этом случае орбита не изменяется при прецессии.

Наглядный пример: спирали Котеса

Рисунок 6: Для голубой частицы, движущейся по прямой линии, радиус р от данного центра изменяется с углом в соответствии с уравнением б = р cos (θ - θ0), куда б расстояние наибольшего сближения (прицельный параметр, показаны красным).

Простейшая иллюстрация теоремы Ньютона возникает, когда нет начальной силы, т.е. F1(р) = 0. В этом случае первая частица неподвижна или движется по прямой. Если он движется по прямой линии, которая не проходит через начало координат (желтая линия на рисунке 6), уравнение для такой линии может быть записано в полярных координатах (р, θ1) в качестве

куда θ0 - угол, при котором расстояние минимизируется (рисунок 6). Расстояние р начинается с бесконечности (когда θ1θ0 = −90°), и постепенно уменьшается до тех пор, пока θ1θ0 = 0°, когда расстояние достигает минимума, то постепенно снова увеличивается до бесконечности при θ1θ0 = 90°. Минимальное расстояние б это прицельный параметр, которая определяется как длина перпендикуляра от фиксированного центра до линии движения. Такое же радиальное движение возможно, когда добавляется центральная сила, обратная кубу.

Рисунок 7: Эписпирали, соответствующие k равно 2/3 (красный), 1,0 (черный), 1,5 (зеленый), 3,0 (голубой) и 6,0 (синий). Когда k меньше единицы, сила обратного куба является отталкивающей, тогда как когда k больше единицы, сила притяжения.

Центральная сила обратного куба F2(р) имеет вид

где числитель μ может быть положительным (отталкивающим) или отрицательным (притягивающим). Если ввести такую ​​силу обратного куба, теорема Ньютона говорит, что соответствующие решения имеют форму, называемую Спирали Котеса[требуется разъяснение ]. Это кривые, определяемые уравнением[20][21]

где постоянная k равно

Когда правая часть уравнения является положительной настоящий номер, решение соответствует эпспиральный.[22] Когда аргумент θ1θ0 равно ± 90 ° ×k, косинус стремится к нулю, а радиус стремится к бесконечности. Таким образом, когда k меньше единицы, диапазон допустимых углов становится небольшим, а сила отталкивающая (красная кривая справа на рисунке 7). С другой стороны, когда k больше единицы, диапазон допустимых углов увеличивается, что соответствует силе притяжения (зеленая, голубая и синяя кривые слева на рисунке 7); орбита частицы может даже несколько раз огибать центр. Возможные значения параметра k может находиться в диапазоне от нуля до бесконечности, что соответствует значениям μ в диапазоне от отрицательной бесконечности до положительного верхнего предела, L12/м. Таким образом, для всех сил притяжения обратного куба (отрицательные μ) существует соответствующая эпспиральная орбита, как и для некоторых отталкивающих (μ < L12/м), как показано на рисунке 7. Более сильные силы отталкивания соответствуют более быстрому линейному движению.

Фигура 8: Спирали Пуансо (шиш спирали), соответствующие значениям λ, равным 1,0 (зеленый), 3,0 (голубой) и 6,0 (синий).

Один из других типов решений представлен в терминах гиперболический косинус:

где постоянная λ удовлетворяет

Эта форма спиралей Котеса соответствует одной из двух Спирали Пуансо (Рисунок 8).[22] Возможные значения λ находятся в диапазоне от нуля до бесконечности, что соответствует значениям μ, превышающим положительное число. L12/м. Таким образом, спиральное движение Пуансо происходит только для отталкивающих центральных сил, обратных кубу, и применяется в том случае, если L не слишком велика для данного μ.

Принимая предел k или λ, стремящаяся к нулю, дает третью форму спирали Котеса, так называемую обратная спираль или же гиперболическая спираль, как решение[23]

куда А и ε - произвольные постоянные. Такие кривые возникают, когда сила отталкивания μ точно уравновешивает член угловой момент-масса

Замкнутые орбиты и центральные силы обратного куба

Рисунок 9: Гармонические орбиты с k = 1 (синий), 2 (пурпурный) и 3 (зеленый). Анимация синей и зеленой орбит показана на рисунке 4.

Два типа центральные силы - линейно увеличивающиеся с расстоянием, F = Cr, Такие как Закон Гука, и силы обратного квадрата, F = C/р2, Такие как Закон всемирного тяготения Ньютона и Закон Кулона - имеют очень необычное свойство. Частица, движущаяся под действием любого типа силы, всегда возвращается в исходное место со своей начальной скоростью, при условии, что ей не хватает энергии, чтобы уйти в бесконечность. Другими словами, путь связанной частицы всегда замкнут, и ее движение повторяется бесконечно, независимо от ее начального положения или скорости. Как показано Теорема Бертрана, это свойство неверно для других типов сил; в общем, частица не вернется в исходную точку с той же скоростью.

Однако теорема Ньютона показывает, что к частице, движущейся под действием линейной силы или силы обратного квадрата, может быть приложена обратная кубическая сила, так что ее орбита остается закрытым при условии, что k равно Рациональное число. (Число называется "рациональным", если его можно записать в виде дроби м/п, куда м и п являются целыми числами.) В таких случаях добавление обратно-кубической силы заставляет частицу завершить м вращение вокруг центра силы за то же время, что и исходная частица. п вращения. Этот метод создания замкнутых орбит не нарушает теорему Бертрана, поскольку добавленная обратно-кубическая сила зависит от начальной скорости частицы.

Рисунок 10: Субгармонические орбиты с k = 1 (синий), 1/2 (пурпурный) и 1/3 (зеленый). Анимация синей и зеленой орбит показана на рисунке 5.

Гармонические и субгармонические орбиты - особые типы таких замкнутых орбит. Замкнутая траектория называется гармоническая орбита если k является целым числом, т.е. если п = 1 в формуле k = м/п. Например, если k = 3 (зеленая планета на рисунках 1 и 4, зеленая орбита на рисунке 9), результирующая орбита является третьей гармоникой исходной орбиты. Наоборот, замкнутая траектория называется субгармоническая орбита если k это обратный целого числа, т.е. если м = 1 в формуле k = м/п. Например, если k = 1/3 (зеленая планета на рисунке 5, зеленая орбита на рисунке 10), получившаяся орбита называется третьей субгармоникой исходной орбиты. Хотя такие орбиты маловероятны в природе, они полезны для иллюстрации теоремы Ньютона.[2]

Предел почти круговых орбит

В предложении 45 его PrincipiaНьютон применяет свою теорему о вращающихся орбитах для разработки метода нахождения силовых законов, управляющих движением планет.[24] Иоганн Кеплер отметил, что орбиты большинства планет и Луны казались эллипсами, и длинная ось этих эллипсов может быть точно определена из астрономических измерений. Длинная ось определяется как линия, соединяющая позиции минимального и максимального расстояний до центральной точки, то есть линия, соединяющая две точки. апсиды. Для иллюстрации длинная ось планеты Меркурий определяется как линия, проходящая через ее последовательные положения перигелия и афелия. Со временем длинная ось большинства вращающихся тел постепенно поворачивается, обычно не более чем на несколько градусов за полный оборот, из-за гравитационных возмущений от других тел. сжатие в притягивающем теле, общие релятивистские эффекты, и другие эффекты. Метод Ньютона использует эту апсидальную прецессию в качестве чувствительного датчика типа силы, применяемой к планетам.[25]

Теорема Ньютона описывает только эффекты добавления центральной силы обратного куба. Однако Ньютон распространяет свою теорему на произвольную центральную силу. F(р), ограничивая свое внимание почти круглыми орбитами, такими как эллипсы с низкой орбитальный эксцентриситет (ε ≤ 0,1), что верно для семи из восьми планетных орбит в Солнечная система. Ньютон также применил свою теорему к планете Меркурий,[26] который имеет эксцентриситет ε примерно 0,21, и предположил, что это может относиться к Комета Галлея, чья орбита имеет эксцентриситет примерно 0,97.[25]

Качественное обоснование этой экстраполяции его метода было предложено Валлури, Уилсоном и Харпером.[25] Согласно их аргументам, Ньютон считал апсидальный угол прецессии α (угол между векторами последовательного минимального и максимального расстояния от центра) как гладкий, непрерывная функция эксцентриситета орбиты ε. Для силы, обратно пропорциональной квадрату, α равно 180 °; векторы к положениям минимального и максимального расстояний лежат на одной прямой. Если изначально α не равен 180 ° при малых ε (квазикруговые орбиты), то, как правило, α будет равняться 180 ° только для отдельных значений ε; случайно выбранное значение ε вряд ли даст α = 180 °. Таким образом, наблюдаемое медленное вращение апсид планетных орбит предполагает, что сила тяжести является законом обратных квадратов.

Количественная формула

Чтобы упростить уравнения, Ньютон пишет F(р) в терминах новой функции C(р)

куда р - средний радиус почти круговой орбиты. Ньютон расширяется C(р) в серии, теперь известной как Расширение Тейлораполномочия расстояния р, одно из первых появлений такого сериала.[27] Приравнивая результирующий член силы обратного куба к силе обратного куба для вращающихся орбит, Ньютон выводит эквивалентный угловой коэффициент масштабирования k для почти круговых орбит :[24]

Другими словами, приложение произвольной центральной силы F(р) на почти круговую эллиптическую орбиту может ускорить угловое движение в раз k без значительного влияния на радиальное движение. Если эллиптическая орбита неподвижна, частица вращается вокруг центра силы на 180 °, перемещаясь от одного конца длинной оси к другому (два апсиды ). Таким образом, соответствующий апсидальный угол α для общей центральной силы равна k× 180 °, используя общий закон θ2 = k θ1.

Примеры

Ньютон иллюстрирует свою формулу тремя примерами.В первых двух центральная сила - это сила закона, F(р) = рп−3, так C(р) пропорциональна рп. Приведенная выше формула показывает, что угловое движение умножается на коэффициент k = 1/п, так что апсидальный угол α равно 180 ° /п.

Этот угловой масштаб можно увидеть в апсидальной прецессии, то есть в постепенном вращении длинной оси эллипса (рис. 3). Как отмечалось выше, орбита в целом вращается со средней угловой скоростью Ω=(k−1)ω, куда ω равна средней угловой скорости частицы относительно стационарного эллипса. Если частица требует времени Т для перехода от одной апсиды к другой это означает, что при этом длинная ось будет вращаться на угол β = ΩТ = (k − 1)ωT = (k - 1) × 180 °. Для закон обратных квадратов Такие как Закон всемирного тяготения Ньютона, куда п равно 1, углового масштабирования нет (k = 1) апсидальный угол α равна 180 °, а эллиптическая орбита стационарна (Ω =β = 0).

В качестве последней иллюстрации Ньютон рассматривает сумму двух степенных законов

который умножает угловую скорость на коэффициент

Ньютон применяет обе эти формулы (степенной закон и сумму двух степенных законов) для изучения апсидальной прецессии орбиты Луны.

Прецессия орбиты Луны

Движение Луны более сложное, чем у планет, в основном из-за конкурирующего гравитационного притяжения Земли и Солнца.

Движение Луна можно точно измерить, и он заметно сложнее, чем у планет.[28] Древнегреческие астрономы, Гиппарх и Птолемей, отметил несколько периодических изменений орбиты Луны,[28] такие как небольшие колебания в его орбитальный эксцентриситет и склонность его орбиты к плоскости эклиптика. Эти колебания обычно происходят раз в месяц или два раза в месяц. Линия его апсиды прецессирует постепенно с периодом примерно 8,85 года, в то время как его линия узлов делает полный круг примерно вдвое больше, чем за 18,6 года.[29] Это объясняет примерно 18-летнюю периодичность затмения, так называемой Цикл Сароса. Однако обе линии испытывают небольшие колебания в своем движении, опять же в месячном масштабе.

В 1673 г. Джеремайя Хоррокс опубликовал достаточно точную модель движения Луны, в которой предполагалось, что Луна следует по прецессирующей эллиптической орбите.[30][31] Достаточно точный и простой метод предсказания движения Луны решил бы навигационную задачу определения движения корабля. долгота;[32] во времена Ньютона целью было предсказать положение Луны с точностью до 2 '(два угловые минуты ), что соответствовало бы ошибке земной долготы в 1 °.[33] Модель Хоррокса предсказывала положение Луны с ошибкой не более 10 угловых минут;[33] для сравнения, диаметр Луны составляет примерно 30 угловых минут.

Ньютон использовал свою теорему о вращающихся орбитах двумя способами для объяснения апсидальной прецессии Луны.[34] Во-первых, он показал, что наблюдаемую апсидальную прецессию Луны можно объяснить изменением силовой закон гравитации от закона обратных квадратов к сила закона в котором показатель был 2 + 4/243 (примерно 2,0165)[35]

В 1894 г. Асаф Холл принял этот подход, слегка изменив показатель степени в законе обратных квадратов, чтобы объяснить аномальный орбитальный прецессия из планета Меркурий,[36] который наблюдал в 1859 г. Урбен Леверье.[37] По иронии судьбы теория Холла была опровергнута тщательными астрономическими наблюдениями Луны.[38] В принятое в настоящее время объяснение поскольку эта прецессия включает в себя теорию общая теория относительности, который (к первое приближение ) добавляет силу обратной четвертичной степени, то есть силу, которая изменяется как обратная четвертая степень расстояния.[39]

В качестве второго подхода к объяснению прецессии Луны Ньютон предположил, что возмущающее влияние солнце на движение Луны может быть приблизительно эквивалентно дополнительной линейной силе

Первый член соответствует гравитационному притяжению между Луной и Землей, где р расстояние Луны от Земли. Второй член, как рассуждал Ньютон, может представлять собой среднюю возмущающую силу гравитации Солнца в системе Земля-Луна. Такой закон силы мог бы также возникнуть, если бы Земля была окружена сферическим пылевым облаком однородной плотности.[40] Используя формулу для k для почти круговых орбит и оценки А и B, Ньютон показал, что этот силовой закон не может объяснить прецессию Луны, поскольку предсказанный апсидный угол α составляла (≈ 180,76 °), а не наблюдаемая α (≈ 181,525 °). За каждый оборот длинная ось поворачивается на 1,5 °, примерно половину наблюдаемых 3,0 °.[34]

Обобщение

Исаак Ньютон впервые опубликовал свою теорему в 1687 году в качестве предложений 43–45 книги I его Philosophiæ Naturalis Principia Mathematica. Однако как астрофизик Субраманян Чандрасекар отметил в своем комментарии 1995 г. PrincipiaТеорема оставалась в значительной степени неизвестной и неразработанной более трех столетий.[1]

Первое обобщение теоремы Ньютона было открыто Магомедом и Вавда в 2000 году.[4] Как и Ньютон, они предположили, что угловое движение второй частицы k раз быстрее, чем у первой частицы, θ2 = k θ1. Однако, в отличие от Ньютона, Магомед и Вауда не требовали, чтобы радиальное движение двух частиц было одинаковым, р1 = р2. Скорее, они требовали, чтобы обратные радиусы были связаны линейным уравнением

Это преобразование переменных изменяет путь частицы. Если путь первой частицы записан р1 = грамм1), путь второй частицы можно записать как

Если движение первой частицы вызывается центральной силой F1(р), Магомед и Вауда показали, что движение второй частицы может быть вызвано следующей силой

Согласно этому уравнению вторая сила F2(р) получается путем масштабирования первой силы и изменения ее аргумента, а также путем добавления центральных сил обратного квадрата и обратного куба.

Для сравнения: теорема Ньютона о вращающихся орбитах соответствует случаю а = 1 и б = 0, так что р1 = р2. В этом случае исходная сила не масштабируется, и ее аргумент не изменяется; сила обратного куба добавляется, а член обратного квадрата - нет. Кроме того, путь второй частицы равен р2 = грамм2/k)в соответствии с приведенной выше формулой.

Производные

Вывод Ньютона

Вывод Ньютона находится в разделе IX его Principia, а именно предложения 43–45.[41] Его выводы этих предложений основаны в основном на геометрии.

Предложение 43; Проблема 30
Диаграмма, иллюстрирующая вывод Ньютона. Голубая планета движется по пунктирной эллиптической орбите, а зеленая планета движется по сплошной эллиптической орбите; два эллипса имеют общие фокус в момент C. Углы UCP и VCQ равны θ1, а черная дуга представляет угол UCQ, который равен θ2 = k θ1. Сплошной эллипс повернулся относительно пунктирного эллипса на угол UCV, равный (k−1) θ1. Все три планеты (красная, синяя и зеленая) находятся на одинаковом расстоянии р из центра силы C.
Требуется заставить тело двигаться по кривой, которая вращается вокруг центра силы так же, как другое тело по той же кривой в состоянии покоя.[42]

Вывод предложения 43 Ньютоном зависит от его предложения 2, полученного ранее в Principia.[43] Предложение 2 обеспечивает геометрическую проверку того, является ли результирующая сила, действующая на точечную массу (частицу) центральная сила. Ньютон показал, что сила является центральной тогда и только тогда, когда частица выметает равные площади за равное время, измеренное от центра.

Вывод Ньютона начинается с частицы, движущейся под произвольной центральной силой. F1(р); движение этой частицы под действием этой силы описывается ее радиусом р(т) от центра как функция времени, а также его угол θ1(т). В бесконечно малое время dt, частица выметает приближенный прямоугольный треугольник, площадь которого равна

Поскольку сила, действующая на частицу, считается центральной силой, частица выметает равные углы за равное время в соответствии с предложением Ньютона 2. Иначе говоря, ставка площади подметания постоянна

Эта постоянная площадная скорость можно рассчитать следующим образом. На апапсис и периапсис, позиции ближайшего и самого дальнего расстояния от центра притяжения, вектора скорости и радиус перпендикулярны; Следовательно угловой момент L1 на массу м частицы (записывается как час1) может быть связано со скоростью выметания участков

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

Поверхностная скорость второй частицы равна скорости первой частицы, умноженной на тот же коэффициент. k

С k - постоянная величина, вторая частица также заметает равные площади за равное время. Следовательно, согласно предложению 2, на вторую частицу также действует центральная сила F2(р). Это заключение предложения 43.

Предложение 44.
Разница сил, с помощью которых два тела могут двигаться одинаково, одно по фиксированной, а другое по одной и той же вращающейся орбите, изменяется обратно пропорционально кубу их общих высот.[44]

Чтобы узнать величину F2(р) от исходной центральной силы F1(р), Ньютон вычислил их разность F2(р) − F1(р) используя геометрию и определение центростремительное ускорение. В предложении 44 его Principia, он показал, что разница пропорциональна кубу, обратному радиусу, в частности, по приведенной выше формуле, которую Ньютон записывает через две постоянные поверхностные скорости: час1 и час2

Предложение 45; Проблема 31
Найти движение апсид на орбитах, очень близких к окружностям.[24]

В этом предложении Ньютон выводит следствия своей теоремы о вращающихся орбитах в пределе почти круговых орбит. Это приближение обычно справедливо для планетных орбит и орбиты Луны вокруг Земли. Это приближение также позволяет Ньютону рассматривать большое количество законов центральной силы, а не только законы обратных квадратов и обратных кубов.

Современное происхождение

Современные выводы теоремы Ньютона опубликованы Whittaker (1937)[45] и Чандрасекхар (1995).[42] По предположению вторая угловая скорость равна k раз быстрее, чем первый

Поскольку два радиуса ведут себя одинаково со временем, р(т) сохраняющиеся угловые моменты связаны тем же множителем k

Уравнение движения для радиуса р частицы массы м переезд в центральный потенциал V(р) дан кем-то Уравнения Лагранжа

Применение общей формулы к двум орбитам приводит к уравнению

который может быть преобразован в форму

Это уравнение, связывающее две радиальные силы, можно качественно понять следующим образом. Разница угловых скоростей (или, что то же самое, угловых моментов) вызывает разницу в центростремительная сила требование; чтобы компенсировать это, радиальная сила должна быть изменена на силу, обратную кубу.

Теорема Ньютона может быть выражена эквивалентным образом через потенциальная энергия, который определен для центральных сил

Уравнение радиальной силы можно записать в терминах двух потенциальная энергия

Интегрируя по расстоянию р, Теорема Ньютона утверждает, что k-кратное изменение угловой скорости является результатом прибавления потенциальной энергии в обратном квадрате к любой заданной потенциальной энергии V1(р)

Геометрическое доказательство Ньютона из принципов

Вращающиеся орбиты Principia 2

Упрощенное геометрическое доказательство предложения 44.

Хотя Ньютон утверждает, что проблема должна была быть решена с помощью предложения 6, он не использует его явно. В следующем упрощенном доказательстве предложение 6 используется, чтобы показать, как выводится результат.

Подробное доказательство Ньютона следует за этим, и, наконец, прилагается предложение 6, поскольку оно малоизвестно.

Предложение 44 использует предложение 6 для доказательства результата о вращающихся орбитах. В предложениях, следующих за предложением 6 раздела 2 Принципов, он применяет его к конкретным кривым, например к коническим сечениям. В случае предложения 44 оно применяется к любой орбите под действием произвольной силы, направленной к фиксированной точке, чтобы создать соответствующую вращающуюся орбиту.

На рис. 1 MN является частью этой орбиты. В точке P тело движется к Q под действием силы, направленной к S, как и раньше. Сила F (SP) определяется в каждой точке P кривой.

На рис. 2 соответствующая часть вращающейся орбиты - это mn с s в качестве центра силы. Предположим, что первоначально тело на статической орбите стартует под прямым углом к ​​радиусу со скоростью V. Тело на вращающейся орбите также должно начинаться под прямым углом и предполагать, что его скорость равна v. В случае, показанном на рис.1, и сила направлена ​​к S. Аргумент одинаково применим, если. Также сила может быть направлена ​​от центра.

Пусть SA будет начальным направлением статической орбиты, а sa - направлением вращающейся орбиты. Если через некоторое время тела на соответствующих орбитах окажутся в точках P и p, то отношения углов ; соотношения площадей; и радиусы, , .

Фигура pryx и дуга py на фиг. 2 - это фигура PRQT и дуга PQ на фиг. 1, растянутые линейно в горизонтальном направлении в соотношении , так что,, и . Прямые qt и QT действительно должны быть дугами окружности с центрами s и S и радиусами sq и SQ соответственно. В пределе их соотношение становится , будь то прямые или дуги.

Поскольку в пределе силы параллельны SP и sp, если бы на тело на фиг. 2 действовала та же сила, что и на фиг. 1, то тело прибыло бы в точку y, поскольку ry = RQ. Разница в горизонтальной скорости не влияет на вертикальные расстояния. Ньютон ссылается на следствие 2 законов движения, в котором движение тел разделяется на составляющую в радиальном направлении, на которую действует вся сила, и на другую составляющую, поперечную к ней, на которую не действует никакая сила.

Однако расстояние от y до центра, s, теперь больше, чем SQ, поэтому для перемещения тела к q требуется дополнительная сила, так что sq = SQ. Дополнительная сила представлена ​​yq, а f пропорциональна ry + yq, так же как F пропорционально RQ.

, .

Разница, , можно найти следующим образом:

,, так .

А в пределе, когда QT и qt стремятся к нулю, становится равным или 2SP так

.

Следовательно,.

Поскольку из предложения 6 (рис.1 и см. Ниже) сила равна.Поделить на , где k - постоянная величина, чтобы получить силы.

На рис. 3 в начальной точке A статической кривой проведите касательную AR, которая перпендикулярна SA, и окружность AQD, которая касается кривой в точке A. Пусть ρ - радиус этой окружности. Поскольку угол SAR является прямым углом, центр круга лежит на SA. Из свойства круга:, а в пределе, когда Q приближается к A, это становится.

Следовательно,.

А поскольку F (SA) задана, это определяет постоянную k. Однако Ньютон хочет, чтобы сила в точке A имела вид, где c - постоянная, так что,куда.

Выражение для f (sp) выше такое же, как у Ньютона в следствии 4 предложения 44, за исключением того, что он использует другие буквы. Он написал (где G и F не обязательно равны v и V соответственно) и использует букву «V» для константы, соответствующей «c», и букву «X» для функции F (sp).

Приведенное выше геометрическое доказательство очень ясно показывает, откуда возникает дополнительная сила, заставляющая орбиту вращаться по отношению к статической орбите.

Доказательство предложения 44 Ньютоном

Доказательство Ньютона сложно ввиду простоты приведенного выше доказательства. Например, его доказательство требует некоторой расшифровки, как показано в следующем предложении:

«И поэтому, если с центром C и любым радиусом CP или Cp описывается круговой сектор, равный общей площади VPC, которую тело P, вращающееся по неподвижной орбите, в любое время описывало радиусом, проведенным к центру, разница между силы, с помощью которых тело P на неподвижной орбите и тело p на подвижной орбите вращаются, будут равны центростремительной силе, с помощью которой какое-либо тело с радиусом, проведенным к центру, могло бы описать этот сектор однородно в том же самом время, когда область VPC описывалась как G2 - F2 к F2.”

Первоначально он рассматривает бесконечно малые при фиксированном значении площади SPQ и spq пропорциональны V и v соответственно; следовательно, и в каждой из точек P и p, и поэтому дополнительная сила изменяется обратно пропорционально кубу радиуса.

На рисунке 1 XQ представляет собой дугу окружности с центром S и радиусом SQ, пересекающую SP в точке X. Перпендикуляр XY пересекает RQ в точке Y и .

Позволять быть силой, необходимой для того, чтобы заставить тело двигаться по окружности радиуса SQ, если она имеет ту же скорость, что и поперечная скорость тела на статической орбите в Q.

в каждой точке P и, в частности, на верхней стороне A:

.

Но в точке A на рис. 3 отношение силы, заставляющей тело двигаться по статической кривой AE, к силе, необходимой для того, чтобы заставить его двигаться по окружности AB с радиусом SA, обратно пропорционально отношению их радиусов. кривизны, поскольку они оба движутся с одинаковой скоростью V перпендикулярно SA:.Из первой части доказательства .

Подставляя выражение Ньютона для F (SA), дает результат, полученный ранее.

Доказательство предложения 45 Ньютоном

«Чтобы найти движение апсидов по орбитам, приближающимся к кругам».

Предложение 44 было разработано специально для доказательства этого предложения.Ньютон хочет исследовать движение тела по почти круговой орбите, притягиваемого силой вида.

Он аппроксимирует статическую кривую эллипсом с силой, обратной квадрату, F (SP), направленной к одному из фокусов, заставленных вращаться добавлением силы, обратной кубу, согласно Предложению 44.

Для статического эллипса с силой, обратно пропорциональной квадрату SP,, поскольку c определено выше, так что .

Когда тело находится на статической орбите, начиная с верхней апсайды в точке A, оно достигнет нижней апсайды, точки, ближайшей к S, после перемещения на угол 180 градусов. Ньютону нужна соответствующая вращающаяся орбита, начинающаяся от апсайды a, вокруг точки s, с нижней апсайдерской стороной, смещенной на угол, α, где .

Начальная скорость V в точке A должна быть чуть меньше, чем требуется, чтобы тело двигалось по кругу. Тогда ρ можно принять равным SA или sa. Проблема состоит в том, чтобы определить v по значению n, чтобы можно было найти α или задать α, чтобы найти n.

Сдача ,

.

Затем «нашим методом сходящихся рядов»: плюс термины в X2 и выше, которые можно игнорировать, потому что орбита почти круглая, поэтому X мало по сравнению с sa.

Сравнивая 2 выражения для f (sp), следует, что.

Также, .

Отношение начальных сил при a определяется выражением.

Предложение 6 для доказательства предложения 44 выше.

На рис.1 тело движется по определенной кривой MN, на которую действует (центростремительная) сила, к фиксированной точке S. Сила зависит только от расстояния точки от S. Цель этого предложения состоит в том, чтобы определить как сила изменяется с радиусом, SP. Метод в равной степени применим к случаю, когда сила центробежная.

За короткое время , тело перемещается из P в ближайшую точку Q. Нарисуйте QR параллельно SP, пересекая касательную в точке R, и QT перпендикулярно SP, пересекаясь с ней в точке T.

Если бы не было силы, он двигался бы по касательной в точке P со скоростью, которую он имел в точке P, достигая точки R. Если бы сила, действующая на тело, движущееся из точки P в точку Q, была постоянной по величине и параллельной направлению в направлении SP дуга PQ будет параболической с PR в качестве касательной, а QR будет пропорционален этой постоянной силе и квадрату времени, .

И наоборот, если вместо того, чтобы достичь R, тело отклонилось на Q, то постоянная сила, параллельная SP, с величиной: заставил бы его достичь Q вместо R.

Однако, поскольку направление радиуса от S к точкам на дуге PQ, а также величина силы по направлению к S будут изменяться вдоль PQ, указанное выше соотношение не даст точной силы в P. Если Q достаточно близко к P, направление силы будет почти параллельно SP на всем протяжении PQ, и, если сила изменится незначительно, можно предположить, что PQ аппроксимируется параболической дугой с силой, указанной выше в терминах QR и .

Время, пропорциональна площади сектора SPQ. Это второй закон Кеплера. Доказательство представлено в предложении 1 книги 1 Принципов. Поскольку дугу PQ можно аппроксимировать прямой линией, площадь сектора SPQ и площадь треугольника SPQ можно принять равными, поэтому

, где k - постоянная величина.

Опять же, это не точно для конечных длин PQ. Закон силы получается, если предел приведенного выше выражения существует как функция SP, когда PQ приближается к нулю.

Фактически, со временем , тело без силы достигло бы точки W дальше от P, чем R. Однако в пределе QW становится параллельным SP. Точка W в доказательстве Ньютона игнорируется.

Кроме того, Ньютон описывает QR как проверенный синус дуги с P в центре и длиной вдвое QP. Хотя это не совсем то же самое, что QR, который он имеет на диаграмме (рис.1), в пределе они становятся равными.

Примечания:

Это предложение основано на анализе Галилеем тела, движущегося по параболической траектории под действием постоянного ускорения. В предложении 10 он описывает это как теорему Галилея и несколько раз упоминает Галилея в связи с ней в «Началах». Объединение его со вторым законом Кеплера дает простой и элегантный метод.

В исторически очень важном случае, когда MN на рис.1 был частью эллипса, а S - одним из его фокусов, Ньютон показал в предложении 11, что предел была постоянной в каждой точке кривой, так что сила, действующая на тело, направленная к фиксированной точке S, изменялась обратно пропорционально квадрату расстояния SP.

Помимо эллипса с центром в фокусе, Ньютон также применил предложение 6 к гиперболе (предложение 12), параболе (предложение 13), эллипсу с центром силы в центре эллипса (предложение 10), равноугольной спираль (Предложение 9), а также окружность с центром силы, не совпадающим с центром, и даже на окружности (Предложение 7).

Смотрите также

Рекомендации

  1. ^ а б Чандрасекхар, стр. 183.
  2. ^ а б c Линден-Белл, Д; Линден-Белл Р.М. (1997). «О формах вращающихся орбит Ньютона». Примечания и отчеты Лондонского королевского общества. 51 (2): 195–198. Дои:10.1098 / рснр.1997.0016. S2CID  73239002.
  3. ^ а б Линден-Белл Д., Джин С. (2008). «Аналитические центральные орбиты и их группа преобразований». Ежемесячные уведомления Королевского астрономического общества. 386 (1): 245–260. arXiv:0711.3491. Bibcode:2008МНРАС.386..245Л. Дои:10.1111 / j.1365-2966.2008.13018.x. S2CID  15451037.
  4. ^ а б Магомед FM, Vawda F (2000). «Применение симметрий к задачам центральной силы». Нелинейная динамика. 21 (4): 307–315. Дои:10.1023 / А: 1008317327402. S2CID  116319304.
  5. ^ Немирофф, Роберт (13 июня 2010 г.). «Ретроградный Марс». Астрономическая картина дня. НАСА. В архиве из оригинала 31 мая 2011 г.. Получено 31 октября 2016.
  6. ^ Sugon QM, Bragais S, McNamara DJ (2008) Эпициклы Коперника из закона гравитационной силы Ньютона через линейную теорию возмущений в геометрической алгебре В архиве 2016-10-29 на Wayback Machine.
  7. ^ Хейльброн 2005, стр.11
  8. ^ Фитцпартрик 2012, стр. 41–43
  9. ^ Ламбурн 2010, стр. 204–205
  10. ^ Уиттакер, стр. 339–385.
  11. ^ Сундман К.Ф. (1912). "Memoire sur le проблема де труа корпус". Acta Mathematica. 36 (1): 105–179. Дои:10.1007 / BF02422379.
  12. ^ Хильтебайтель А.М. (1911). «К проблеме двух неподвижных центров и некоторых ее обобщениях». Американский журнал математики. Издательство Университета Джона Хопкинса. 33 (1/4): 337–362. Дои:10.2307/2369997. JSTOR  2369997.
  13. ^ Хейльброн 2005, стр.139
  14. ^ Коэн, стр. 147.
  15. ^ Clairaut, AC (1745). "Du Système du Monde dans les Principes de la Gravitation Universelle". Histoire de l'Académie Royale des Sciences avec les mathématique et de Physique. 1749: 329–364. В архиве из оригинала 2011-06-07. Получено 2007-07-12.
  16. ^ Hill GW (1894 г.). «Буквальное выражение движения перигея Луны». Анна. Математика. 9 (1/6): 31–41. Дои:10.2307/1967502. JSTOR  1967502.
  17. ^ Браун EW (1891). «Неизвестный титул». Являюсь. J. Math. Издательство Университета Джона Хопкинса. 13 (2): 159–172. Дои:10.2307/2369812. JSTOR  2369812.
    Браун EW (1891). «Об определении определенного класса неравенств в движении Луны». Ежемесячные уведомления Королевского астрономического общества. 52 (2): 71. Bibcode:1891МНРАС..52 ... 71Б. Дои:10.1093 / mnras / 52.2.71.
  18. ^ Делоне С (1862 г.). «Неизвестный титул». Mémoires Acad. Imp. Sc.: 237.
    Делоне С (1867 г.). «Неизвестный титул». Mémoires Acad. Imp. Sc.: 451.
  19. ^ Ньютон, Principia, раздел IX книги I, предложения 43–45, стр. 135–147.
  20. ^ Гроссман 1996, стр. 33–34
  21. ^ Шикин 1995, стр. 139–140
  22. ^ а б Лоуренс 1972, стр. 192–194
  23. ^ Вайсштейн 2002, стр. 1427
  24. ^ а б c Чандрасекхар С 1995, стр. 192–194
  25. ^ а б c Валлури С. Р .; Wilson C .; Харпер В. (1997). "Апсидальная теорема Ньютона о прецессии и эксцентрические орбиты". Журнал истории астрономии. 28: 13–27. Bibcode:1997JHA .... 28 ... 13В. Дои:10.1177/002182869702800102. S2CID  117886193.
  26. ^ Ньютон, Principia, Книга III, Предложение 2, с. 406.
  27. ^ Коэн И.Б. (1990). "Два эссе Галлея о Ньютоне" Principia". В Нормане Троуэре (ред.). Стоя на плечах гигантов: более длинный взгляд на Ньютона и Галлея. Беркли, Калифорния: Калифорнийский университет Press. стр.91–108. ISBN  978-0-520-06589-5.
  28. ^ а б Повар А (2000). «Успех и неудача в лунной теории Ньютона». Астрономия и геофизика. 41 (6): 21–25. Bibcode:2000A и G .... 41f..21C. Дои:10.1046 / j.1468-4004.2000.41621.x.
  29. ^ Смит, стр. 252.
  30. ^ Хоррокс Дж. (1673). Иеремия Гороци, опера posthuma. Лондон: G Godbit для Дж. Мартина.
  31. ^ Уилсон C (1987). «О происхождении лунной теории Хоррока». Журнал истории астрономии. 18 (2): 77–94. Bibcode:1987JHA .... 18 ... 77 Вт. Дои:10.1177/002182868701800201. S2CID  115379870.
  32. ^ Коллерстрем Н (2000). Забытая лунная теория Ньютона: его вклад в поиски долготы. Green Lion Press. ISBN  978-1-888009-08-8.
  33. ^ а б Смит, стр. 254.
  34. ^ а б Ньютон, Principia, Книга I, раздел IX, предложение 45, стр. 141–147.
  35. ^ Чандрасекхар, стр. 198.
  36. ^ Зал А (1894 г.). «Предложение в теории Меркурия». Астрономический журнал. 14: 49–51. Bibcode:1894AJ ..... 14 ... 49H. Дои:10.1086/102055.
  37. ^ Le Verrier UJJ (1859 г.). "Теория движения Меркурия". Annales de l'Observatoire Impérial de Paris. 5: 1–196, особенно. 98–106. Bibcode:1859АнПар ... 5 .... 1л.
    Саймон Ньюкомб (1882 г.). «Обсуждение и результаты наблюдений за транзитами Меркурия с 1677 по 1881 год». Астрономические документы, подготовленные для использования в американских эфемеридах и морском альманахе. 1: 473. Bibcode:1882USNAO ... 1..363N.
  38. ^ Браун EW (1903). «О степени точности новой теории Луны». Ежемесячные уведомления Королевского астрономического общества. 64: 524–534. Bibcode:1904МНРАС..64..524.. Дои:10.1093 / минрас / 64.6.524.
  39. ^ Розевир Н. (1982). Перигелий Меркурия от Леверье до Эйнштейна. Оксфорд.
  40. ^ Симон К.Р. (1971). Механика (3-е изд.). Ридинг, Массачусетс: Аддисон – Уэсли. С. 267 (Глава 6, проблема 7). ISBN  0-201-07392-7.
  41. ^ Чандрасекар, стр. 183–192.
  42. ^ а б Чандрасекхар, стр. 184.
  43. ^ Чандрасекар, стр. 67–70.
  44. ^ Чандрасекхар, стр. 187.
  45. ^ Уиттакер, стр. 83.

Библиография

дальнейшее чтение

внешняя ссылка