СИДЕРИЧЕСКИЕ ВАРИАЦИИ ПАРАМЕТРОВ ОРБИТ GPS: ЕЩЁ ОДНО ПОДТВЕРЖДЕНИЕ НОВОЙ КОНЦЕПЦИИ ТЯГОТЕНИЯ

 

А.А.Гришаев,  независимый исследователь

 

 

Введение.

Анализ траекторий искусственных спутников Земли (ИСЗ) может служить мощным средством верификации тех или иных моделей тяготения. Так, параметры орбит спутников GPS испытывают ярко выраженные вариации с периодом в сидерический месяц. Эти вариации не находят объяснения в рамках традиционной концепции всемирного тяготения – тогда как они прямо следуют из наших представлений о тяготении.

Как мы излагали ранее [1], движение космических аппаратов в окололунном пространстве убедительно свидетельствует о малости области действия лунного тяготения. Граница этой области отстоит примерно на 10000 км от поверхности Луны – и, таким образом, лунное тяготение далеко не достаёт до Земли. Поэтому Земля, вопреки распространённому заблуждению, не обращается, в противофазе с Луной, около их общего центра масс. Вместо этого имеет место модуляция, с периодом в синодический месяц, скорости движения земной частотной воронки [2], т.е. области земного тяготения, по орбите вокруг Солнца. В результате, в системе Земля-Луна, Земля совершает одномерные колебания около условного центра, вокруг которого Луна выписывает двумерную траекторию. Такой подход даёт немедленное, чисто геометрическое, объяснение одного из главных неравенств в движении Луны – т.н. вариации [3]. Кроме того, факт одномерности колебаний Земли в кинематике пары Земля-Луна [4] обнаруживается, например, при радиолокации планет, при приёме импульсов пульсаров, при слежении за межпланетными космическими аппаратами. Мы полагаем, что назначение модуляции орбитальной скорости у земной частотной воронки – быть синхронизатором обращения Луны вокруг Земли [3], т.е. быть параметрическим задатчиком периода этого обращения.

Действительно: земная частотная воронка, на наш взгляд, задаёт околоземный «инерциальный фон» для механических и электромагнитных процессов [5]. Ускорение этого «инерциального фона» земной частотной воронки – вперёд-назад вдоль местного участка околосолнечной орбиты – даёт соответствующие поправки в локальные векторы силы тяжести. Именно эти поправки, как мы постарались показать, являются истинной причиной лунных приливов в океанах [4]. Что касается свободного движения тел в области земного тяготения, то дополнительное ускорение «инерциального пространства», в котором они движутся, должно порождать «инерциальный снос», возмущающий их движение. Это должно быть справедливо не только в случае с Луной [3], но и в случае с ИСЗ. В данной статье мы дадим теоретические предсказания эволюций больших полуосей орбит GPS из-за колебаний, вперёд-назад, земной частотной воронки при её орбитальном движении вокруг Солнца, а также покажем, как эти предсказания согласуются с экспериментальными данными, полученными при слежении за спутниками с помощью системы NORAD.

 

Теоретические предсказания.

Первая оценка массы Луны, на основе закона всемирного тяготения, была получена через периодическую компоненту видимой долготы Солнца, с периодом в синодический месяц (29.53 сут), что интерпретировали как колебания гелиоцентрической долготы Земли – т.н. лунное неравенство [6]. Заметим: лунное неравенство как раз свидетельствует о том, что Земля, при своём орбитальном движении, совершает колебания вперёд-назад, но отнюдь не вправо-влево – которые тоже имели бы место, если бы Земля обращалась около общего с Луной центра масс. Из амплитуды лунного неравенства, которая составляет, в линейном исчислении, около 4670 км, и периода в синодический месяц, для амплитуды ускорения при колебаниях земной частотной воронки вперёд-назад мы получаем величину 2.83×10-5 м/с2.

Мгновенному вектору этого ускорения противоположен мгновенный вектор ускорения «инерциального сноса» малых тел, свободно движущихся в пределах земной частотной воронки. Движение спутников GPS возмущается этим ускорением «инерциального сноса», которое на одной половине синодического месяца направлено по вектору орбитальной скорости Земли, а на другой половине – против этого вектора. Результирующие эволюции параметров орбит GPS должны иметь периодические компоненты – которые, на первый взгляд, должны иметь период в синодический месяц. Однако, при движении пары Земля-Луна вокруг Солнца, орбиты GPS почти не изменяют своих ориентаций относительно неподвижных звёзд – поэтому, как и в случае с Луной, геоцентрическая долгота которой повторяется с периодом в сидерический месяц (27.32 сут), период вариаций орбит GPS также должен быть равен сидерическому месяцу.

Для оценки амплитуд этих вариаций можно приближённо считать, что, на протяжении сидерического месяца, не изменяется ориентация прямой, вдоль которой действует ускорение «инерциального сноса», относительно розетки орбит GPS. Рассмотрим случай, когда плоскость орбиты параллельна линии действия «инерциального сноса». Мы имеем задачу об эволюциях почти круговой 12-часовой орбиты в условиях постоянно действующего малого возмущающего ускорения, вектор которого медленно изменяется с периодом в сидерический месяц, действуя вдоль зафиксированной в пространстве линии.

Такая задача, насколько нам известно, не имеет аналитического решения. Поэтому мы решали её с помощью численного интегрирования орбиты, причём, использовали предложенный нами в 2001 году упрощённый метод [7], работающий для почти круговых орбит – в котором учитывается не только центростремительное ускорение к «силовому центру», но и центробежное ускорение, а также, в рассматриваемом случае, ещё и возмущающее ускорение. При этом задача решается для плоской орбиты в полярных координатах, и текущими параметрами являются: модуль радиус-вектора R, угол q поворота радиус-вектора (отсчитываемый от точки перигея), модуль вектора скорости V и угол b между вектором скорости и перпендикуляром к радиус-вектору. При шаге интегрирования T, алгоритм имеет вид:

;                                                                      (1)

;                                             (2)

;                                                                      (3)

,                                         (4)

где Vt=Vcosb и Vr=Vsinb - тангенциальная и радиальная компоненты скорости, at=a×cos(q+j0) и ar=a×sin(q+j0) – аналогичные компоненты возмущающего ускорения, причём a=a0sin(2pt/23.32 сут), a0=2.83×10-5 м/с2, t – текущее время, j0 – стартовый угол между вектором возмущающего ускорения и линией «перигей-апогей», GM – произведение гравитационной постоянной на массу Земли. Численное интегрирование по алгоритму (1)-(4) даёт зависимость длины радиус-вектора от времени – которая позволяет судить об эволюции большой полуоси орбиты. Для выполнения этих расчётов и визуализации эволюции большой полуоси мы написали программу GPS-num-integrator [8]; конечный результат при шаге интегрирования T=2 с для случая j0=1.57 рад и эксцентриситета 0.005 приведён на Рис.1. Расчёты проводились для различных значений эксцентриситетов, в диапазоне от 0.002 до 0.012, характерном для орбит GPS (см. ниже ссылку на данные NORAD). При этом стартовой являлась точка перигея с идеально соответствующим ей

 

           

            Рис.1. Сидерическая вариация большой полуоси орбиты GPS, рассчитанная

по программе GPS-num-integrator. Параметры см. в тексте.

 

вектором перигейной скорости – и по модулю, и по направлению (b=0). Результирующие амплитуды сидерических вариаций большой полуоси, как функции угла j0 при различных значениях эксцентриситета, приведены на Рис.2. Можно видеть, что сидерическая вариация, практически, линейно растёт при увеличении эксцентриситета, и что при заданном эксцентриситете она максимальна, когда линия действия «инерциального сноса» перпендикулярна линии «перигей-апогей».

           

           

            Рис.2. Расчётные амплитуды сидерических вариаций больших полуосей орбит GPS.

 

Теперь вспомним, что расчёты выполнялись для упрощённого частного случая, когда вектор «инерциального сноса» на протяжении всего сидерического месяца находился в плоскости орбиты. Но космический сегмент GPS имеет шесть орбитальных плоскостей. Все орбиты GPS имеют одинаковые номинальные наклонения, около 55о, номинальные же долготы восходящих узлов для орбит из орбитальных плоскостей A,B,C,D,E,F составляют, соответственно, 30о,90о,150о,210о,270о,330о. Причём, такая геометрия орбит GPS имеет место в геоцентрической системе координат – у которой ось OZ, проходящая через северный полюс Земли, наклонена на 23о.5 к плоскости эклиптики, в которой всегда лежит вектор «инерциального сноса». В общем случае, этот вектор будет иметь две компоненты: параллельную плоскости орбиты и ортогональную к ней. Но, как известно, малая тяга, ортогональная к плоскости орбиты, не приводит к эволюции её большой полуоси [9]. Поэтому для нашей задачи представляет интерес лишь величина проекции вектора «инерциального сноса» на ту или иную орбитальную плоскость. После стандартного преобразования декартовых координат при переходе от системы, у которой оси OX и OY лежат в плоскости эклиптики, причём ось OX проходит через точку весеннего равноденствия, к геоцентрической системе, получаемой поворотом вокруг оси OX эклиптической системы на 23о.5, можно для каждой орбитальной плоскости получить коэффициент, умножение на который «инерциального сноса» даёт его проекцию на эту орбитальную плоскость. Временные зависимости этих коэффициентов для каждой из шести орбитальных плоскостей GPS показаны на Рис.3 на годичном интервале – от одного весеннего равноденствия до следующего.

 

Рис.3. Коэффициенты проекций «инерциального сноса» на орбитальные плоскости

A,B,C,D,E,F – на протяжении одного года.

 

Сопоставление диаграмм на Рис.2 и Рис.3 позволяет сформулировать предсказания нашего подхода: для характерных параметров орбит GPS можно ожидать амплитуду сидерической вариации большой полуоси примерно от 20 метров и менее – в зависимости от эксцентриситета орбиты, её принадлежности к той или иной из шести орбитальных плоскостей, а также фазы обращения спутника. Посмотрим, как эти предсказания согласуются с экспериментальными реалиями.

 

Экспериментальные реалии.

Параметры орбитального движения спутников GPS, по результатам работы NORAD, мы брали из базы данных официального портала [10]. В выдаваемом там наборе параметров орбитального движения – для каждого момента из дискретного ряда отсчётов – отсутствует величина большой полуоси, но она однозначно рассчитывается на основе периода обращения (фактически – на основе количества витков за сутки). На построенных таким образом диаграммах, отображающих зависимости больших полуосей от времени, месячная волна легко заметна невооружённым глазом, но точно определить период и амплитуду этой волны – затруднительно. Поэтому мы написали программу NORAD-Spectr [11], которая производит спектральный анализ вариаций большой полуоси орбиты ИСЗ, получаемых на основе данных портала [10]. При этом в программе NORAD-Spectr применён метод спектрального анализа, аналога которого нам не удалось отыскать.

Здесь камнем преткновения для известных методов спектрального анализа является то, что обрабатываемый дискретный ряд данных имеет, помимо зашумлённой периодической компоненты, ещё и низкочастотный тренд. В этом случае плохо работает традиционная методика умножения исследуемого ряда F(t) на гармоническую функцию sin(wt+x) и, для исследуемого временного интервала t2-t1, нахождения определённого интеграла

,                                                                     (5)

который, при подходящем значении x, максимален для спектральной компоненты с частотой w. Проблема в том, что низкочастотный тренд даёт паразитные вклады в амплитуды спектральных компонент, вычисляемых через интеграл (5) – отчего в искомый спектр добавляются искажения, которые проблематично учесть. Впрочем, известны методы устранения этого влияния тренда. Так, в методе спектрального анализа, основанном на использовании т.н. вариации Аллана [12], низкочастотный тренд игнорируется благодаря тому, что анализируется не ряд разностей между значениями функции F(t) и её средним на интервале t2-t1, а ряд разностей между последовательными значениями функции F(t). Однако, имеется серьезное ограничение на применимость такого метода: он работает лишь для рядов данных, в которых последовательные моменты отсчётов значений F(t) разделены одинаковыми дискретами – т.е. отсчёты делаются раз в секунду, раз в час, раз в сутки и т.п. Это же ограничение на применимость имеет и такой известный метод спектрального анализа, как быстрое преобразование Фурье. Поэтому названные методы не подходят для обработки рядов данных NORAD – где моменты отсчётов параметров орбиты разделены неодинаковыми промежутками. Хуже того, в рядах данных NORAD встречаются скачки параметров орбиты – например, из-за импульсных манёвров космического аппарата. Насколько нам известно, традиционные методы спектрального анализа не способны корректно обрабатывать подобные скачки: приходится искусственно, с некоторой степенью произвола, делать «сшивку» ряда данных. Метод же спектрального анализа, применённый в программе NORAD-Spectr, позволяет обрабатывать дискретные ряды данных со стохастическими промежутками между моментами последовательных отсчётов – причём в этом методе изящно игнорируются как низкочастотный тренд, так и скачки в потоке данных.

Идея метода основана на том, что анализируется не исходная функция F(t), а её производная. Ясно, что каждая гармоническая компонента, присутствующая в исходной функции, будет присутствовать и в её производной – правда, с соответственно изменёнными амплитудой и фазой; но эти изменения легко учесть. Преимущества же анализа не исходной функции, а её производной, заключаются в следующем. Во-первых, низкочастотный тренд при взятии производной превращается в постоянную составляющую, которую легко удалить – благодаря чему можно устранить паразитные вклады или, по крайней мере, минимизировать их. Во-вторых, скачок в потоке данных даёт единичный аномальный выброс значения производной, который легко отбросить – без ущерба для корректности последующего анализа, ведь «сшивка» здесь не требуется.

Эти принципы и реализованы в программе NORAD-Spectr, которая производит следующие операции.

1. Визуализация исходного ряда значений большой полуоси, полученных пересчётом из значений количества витков за сутки.

2. Взятие производной от исходного ряда и её визуализация.

3. Отбрасывание аномальных выбросов значений производной, лежащих за пределами утроенного средне-квадратического отклонения. Эту процедуру можно выполнять неоднократно – до тех пор, пока производная не приобретёт вид без аномальных выбросов.

4. Сдвиг графика производной по оси ординат таким образом, чтобы среднее значение производной стало равно нулю.

5. Спектральный анализ результирующей производной. Для каждой спектральной компоненты в полосе анализа, находится набор интегралов (5) при различных фазах x, и выбирается максимальное значение интеграла, дающее максимальную амплитуду спектральной компоненты.

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

На Рис.4-9 показаны некоторые результаты – исходные временные зависимости большой полуоси орбиты и их спектры в окрестностях компоненты с периодом в сидерический месяц. Указаны орбитальная плоскость, номер спутника по каталогу NORAD, эксцентриситет и годичный интервал исходного массива данных. Как можно видеть, с полученными спектрами весьма неплохо согласуются наши теоретические предсказания, изложенные выше.

 

 

 

Рис.4, Рис.5. Большая полуось, по данным NORAD, и спектр её вариаций.

Плоскость B, № 27663, эксц. 0.005-0.006, 21 марта 2009 – 21 марта 2010.

 

 

 

 

Рис.6, Рис.7. Большая полуось, по данным NORAD, и спектр её вариаций.

Плоскость E, № 28129, эксц. 0.005, 21 марта 2009 – 21 марта 2010.

 

 

 

 

Рис.8, Рис.9. Большая полуось, по данным NORAD, и спектр её вариаций.

Плоскость A, № 32711, эксц. 0.001-0.002, 21 марта 2008 – 21 марта 2009.

 

 

Небольшое обсуждение.

Налицо факт переменных составляющих в вариациях больших полуосей орбит GPS – с периодом в сидерический месяц и амплитудами порядка или менее 20 метров. Этот факт мы считаем прямым следствием того, что спутники GPS испытывают «инерциальный снос» из-за модуляции орбитальной скорости земной частотной воронки с периодом в синодический месяц, что эквивалентно колебаниям земной частотной воронки вперёд-назад с амплитудой около 4670 км, наложенным на её орбитальное движение вокруг Солнца.

В рамках же традиционных представлений о тяготении, периодические вариации больших полуосей орбит GPS, о которых идёт речь, не находят разумного объяснения. Действительно, поскольку период этих вариаций равен сидерическому месяцу, то их единственной причиной, в рамках концепции всемирного тяготения, могли бы быть возмущения движения спутников GPS из-за тяготения Луны. Но, даже если тяготение Луны доставало бы до Земли, то оно сообщало бы, в первом приближении, одинаковые ускорения и Земле, и спутникам GPS. Возмущения в движение спутников были бы обусловлены дифференциальной разностью ускорений, сообщаемых Луной Земле и спутникам. Величина этой разности была бы максимальна, когда спутник находился бы либо со стороны Луны от Земли, либо с противоположной стороны – но даже тогда эта величина была бы на порядок меньше, чем значение 2.83×10-5 м/с2, на основе которого получены вышеприведённые результаты. Хуже того, возмущающее ускорение – при движении спутника как по ближнему к Луне сегменту орбиты, так и по дальнему от неё сегменту орбиты – было бы направлено наружу от орбиты. При этом приращения большой полуоси, набираемые на этих сегментах, имели бы противоположные знаки [9], и результирующее приращение большой полуоси на одном витке было бы, практически, нулевое. Поэтому сидерическим вариациям большой полуоси с амплитудой ~20 м было бы неоткуда взяться. Вот почему специалисты помалкивают об этих вариациях.

Как можно видеть, особенности движения спутников GPS в очередной раз указывают на ограниченность применимости концепции т.н. всемирного тяготения и, в очередной раз, подтверждают не только работоспособность нашей модели тяготения [13], но и её предсказательную силу.

 

 

Ссылки.

 

1.                         А.А.Гришаев. Граница области тяготения Луны: анализ полётов в окололунном пространстве. – Доступна на данном сайте.

2.                         А.А.Гришаев. Взаимное тяготение звёзд и планет обусловлено… алгоритмически? – Доступна на данном сайте.

3.                         А.А.Гришаев. Синхронизатор орбитального движения Луны. – Доступна на данном сайте.

4.                         А.А.Гришаев. Свидетельства об одномерности колебаний Земли в кинематике пары Земля-Луна. – Доступна на данном сайте.

5.                         А.А.Гришаев. Иерархия частотных склонов в роли «светоносного эфира». – Доступна на данном сайте.

6.                         М.Ф.Субботин. Введение в теоретическую астрономию. «Наука», М., 1968.

7.                         А.А.Гришаев. Частотно-градиентная природа центробежных сил. – Доступна на данном сайте.

8.                         Программа GPS-num-integrator и описание к ней доступны на данном сайте:  http://newfiz.narod.ru/soft/soft.htm

9.                         К.Б.Алексеев, Г.Г.Бебенин, В.А.Ярошевский. Маневрирование космических аппаратов. «Машиностроение», М., 1970.

10.                      http://www.space-track.org

11.                      Программа NORAD-Spectr, описание к ней и демонстрационные файлы доступны на данном сайте: http://newfiz.narod.ru/soft/soft.htm

12.                      D.W.Allan. Time and Frequency (Time-Domain) Characterization, Estimation and Prediction of Precision Clocks and Oscillators. IEEE Trans. on Ultrasonics, Ferroelectrics and Freq. Control. Vol. UFFC-34, N6, November 1987, p.647.

13.                      А.А.Гришаев. Книга «Этот «цифровой» физический мир», 2010. – Доступна на данном сайте.

 

Источник:  http://newfiz.info

Поступило на сайт: 15 декабря 2010.