Скорость гравитации

………………….                      

                                                    

 

                                                            СКОРОСТЬ ГРАВИТАЦИИ

 

                                                                             Часть 4

ВЛИЯНИЕ СКОРОСТИ ГРАВИТАЦИИ НА СМЕЩЕНИЯ ПАРАМЕТРОВ ОРБИТ ПЛАНЕТ

 

…................................... третья редакция (переработанная и дополненная) 08.09.2014      

…................................... вторая редакция (переработанная и дополненная) 01.12.2013      

                                                  …            .первая редакция 07.07.2013                                                                                 

 

С. Ю. Юдин          http://modsys.narod.ru  .                                                                     . modsys@yandex.ru

                               https://googledrive.com/host/0BwnV2Ac6zvhMalpOVGktQ1Jic1U

 

 

 

 

                                                                 навигация по статье

стр.2 – исторический обзор  

стр.4 –  эффект запаздывания потенциалов по координатам

стр.11 – эффект динамического давления поля

стр.32 – учет эффекта динамического давления поля в формуле учитывающей запаздывание потенциалов по координатам 

стр.36 – практическая проверка различных формул, учитывающих конечность скорости гравитации, на примере Солнечной системы.

стр.46 – выводы

стр.48 – список литературы

 

стр.50 – приложение 1 (экспериментальные данные при учете эффекта запаздывания по координатам)

стр.63 – приложение 2 (модификация метода Рунге-Кутта для расчета ускорений по ускорениям, которая необходима при использовании формул Вебера, Гербера и Лиенара-Вихерта)

стр.72 – приложение 3 (экспериментальные данные при учете и эффекта запаздывания по координатам и эффекта динамического давления)

стр.76 – приложение 4 (расчеты скорости гравитации Лапласа и Ван Фландерна, а также данные по двойным пульсарам и методика их получения на примере пульсара PSR 1913+16)

 

 

 

В статье рассматриваются различные варианты усовершенствования закона тяготения Ньютона, где авторы пытаются учесть конечность скорости распространения гравитации. Это и эффект сопротивления гравитационной жидкости у Лапласа и потенциалы Вебера, Гербера и Лиенара-Вихерта, а также геометрическое решение этой задачи в ОТО. При этом я предлагаю и свой вариант решения этой задачи, а именно использование потенциалов запаздывающих по координатам, где предлагаю и свое решение учета динамического давления гравитационного поля. Только не надо путать мой вариант запаздывания потенциалов по координатам с тем, что дается в учебниках и который использовал Ван Фландерн в своих расчетах, т.к. там дается смесь моего вариант запаздывания потенциалов по координатам и варианта Лапласа с эффектом сопротивления гравитационной жидкости, что является элементарной ошибкой и подробно об этой ошибке я пишу в приложении 4. А в приложении 2 я привожу свою модификацию метода Рунге-Кутта для расчета ускорений по ускорениям, которая необходима при использовании формул Вебера, Гербера и Лиенара-Вихерта.

 

 

Эта третья редакция статьи отличается от второй тем, что я теперь, найдя первоисточники, не называю ошибочно потенциалы Лиенара-Вихерта потенциалами запаздывающими по координатам, а рассматриваю их как совершенно другие потенциалы. А также, кроме потенциалов Лиенара-Вихерта я теперь рассматриваю и решение с использованием ОТО. При этом из статьи убран материал по различным подходам для описания явлений Природы и по эффекту Доплера, которые рассмотрены мною в отдельных статьях (см. ссылки в списке литературы), а также я убрал материал, где я рассматриваю нарушение законов сохранения энергии и момента количества движения при использовании потенциалов запаздывающих по координатам, т.к. и по этому вопросу планирую написать отдельную статью, где планирую провести более точные вычислительные эксперименты с использованием третьей версии программы Solsys7. Таким образом, в статье остался только материал непосредственно относящийся к влиянию скорости гравитации при различных модернизациях закона тяготения Ньютона на вековые смещения параметров орбит планет, что позволило значительно сократить ее объем и не распылять внимание читателей на другие вопросы. 

 

 

 

 

 

                                                                        Исторический обзор

 

 О скорости распространения гравитации ученые стали задумываться буквально сразу после открытия Гуком и Ньютоном закона тяготения, но и сейчас сказать, что-то определенное о величине этой скорости никто не может. Лаплас, исходя из данных наблюдений за Луной, сделал прикидочный расчет и получил, что для этого скорость гравитации должна быть значительно больше скорости света. Конкретно он писал [41] "Мгновенно ли передается притяжение от одного тела к другому? Продолжительность его передачи, если бы она была для нас ощутима, обнаружилось бы главным образом в вековом ускорении движения Луны. Я предлагал таким способом объяснить наблюдаемое ускорение этого движения и нашел, что удовлетворить наблюдениям можно, лишь приписав силам притяжения скорость, в 7 000 000 раз большую скорости светового луча. Так как причина векового уравнения Луны в настоящее время хорошо известна, мы можем утверждать, что тяготение передается, по крайней мере, в 50 000 000 раз быстрее света. Поэтому, не боясь внести заметную ошибку, можно считать его распространение мгновенным".

 

А уже современные исследования Ван Фландерна, который определял скорость гравитации по периоду обращения двойного пульсара, дали еще более значительный результат. Конкретно он писал [42] "Using the same equation with binary pulsar PSR1534+12 and the parameters in Table I, we can place the most stringent limit yet from the observed uncertainty in dP/dt: Vgr => 2*10^10 c". Отсюда можно сделать пессимистический вывод о том, что определить скорость гравитации мы никогда не сможем, т.к. приведенные Лапласом и Ван Фландерном значения скорости гравитации для наших приборов означают практически бесконечность, т.е. никакими нашими приборами зафиксировать такую скорость не возможно. И получается, что наука пришла к тому с чего она и начиналась после открытия закона всемирного тяготения (1), а именно к тому, что Ньютон не стал измышлять гипотез и принял в своей теории тяготения скорость гравитации равной бесконечности. К счастью, это не совсем так и результаты Лапласа и Ван Фландерна можно смело проигнорировать, т.к. они не имеют никакого отношения к действительности.

 

F= G*M*m/R^2                 (1)

 

Так Лаплас пытался объяснять с помощью конечности скорости гравитации увеличение орбитальной угловой скорости Луны, но в математическом (эфемеридном) времени, в котором Лаплас делал расчеты, эта скорость не увеличивается, а уменьшается. При этом, когда Лаплас отказался от своей ошибочной оценки скорости гравитации, которая давала превышение скорости света в 7 000 000 раз, и выдал (неизвестно из каких туманных соображений) оценку в 50 000 000 раз, то он даже это сделал ошибочно, т.к. он ошибочно определил и то, что ускорение орбитальной угловой скорости Луны объясняется не конечностью скорости гравитации, а изменением эксцентриситета орбиты Земли. А Ван Фландерн неизвестно зачем искал с учетом скорости гравитации теоретическое уменьшение угловой скорости у пульсаров, хотя по приведенным им якобы экспериментальным данным эта скорость у них увеличивается. К тому же, как я выяснил, надежность этих якобы наблюдаемых данных равна нулю, а используемое им наблюдаемое значение увеличения радиуса орбиты Земли является очень не надежным, о чем я писал в [45]. А, чтобы выяснить достоверность якобы наблюдаемого уменьшения периода обращения двойных пульсаров я даже провел целое исследование в этой статье (приложение 4), которое показало, что это значение вообще не является наблюдаемым.

 

 

При этом, я рассмотрел на примере самого известного двойного пульсара PSR 1913+16 методику определения всех параметров систем двойных пульсаров, где обнаружил массу различных ошибок или, что еще хуже, подтасовок. Например, интересующее нас якобы наблюдаемое уменьшение периода обращения пульсара по орбите из тех же данных наблюдений может получиться и увеличением его периода обращения. А в таком случае уравнение ОТО, которое приводит к уменьшению периода обращения из-за излучения гравитационных волн, использовать при определении "наблюдаемых" параметров орбит двойных пульсаров нельзя. Причем, его нельзя использовать даже в том случае, если выяснится, что период обращения двойных пульсаров действительно уменьшается, т.к. это может происходить и по другим причинам. В общем, использованные Ван Фландерном якобы наблюдаемые параметры орбит двойных пульсаров таковыми никак не являются.

 

 

А если к этому добавить еще и то, что в расчетах Лапласа и Ван Фландерна имеются и ошибки, то говорить о том, что они определяли скорость гравитации, можно только в историческом аспекте. Если коротко, то Лаплас учел только динамическое давление гравитационной жидкости при движении планет по круговым орбитам, которое дает тангенциальную составляющую силы взаимодействия между планетой и центральным телом, а при движение по эллипсам появится еще и дополнительная радиальная составляющая от давления гравитационной жидкости при удалении планеты от центрального тела и при приближении. А Ван Фландерн запутался в звездной аберрации и в планетной (он ее называет транзитной задержкой) и в результате саму величину тангенциальной силы определил как Лаплас по аналогии со звездной аберрацией, а знак взял тот, который получается от эффекта планетной аберрации, т.е. от запаздывания потенциалов по координатам, и в результате получил ту же по величине тангенциальную силу, что и Лаплас, но с другим знаком и поэтому у него большая полуось орбит не уменьшается, как у Лапласа, а увеличивается.

 

 

Кого интересуют подробности расчетов Лапласа и Ван Фландерна с разбором их ошибок, а также подробности по определению параметров орбит двойных пульсаров я отсылаю к Приложению 4, а здесь перейду к изложению основного содержания статьи. Так вот, я считаю, что определить экспериментально скорость гравитации все же можно и такой эксперимент уже поставила сама Природа, а нам надо только правильно обработать экспериментальные данные, которые мы зафиксировали в этом эксперименте. А конкретно правильно определить по данным оптических наблюдений вековые смещения параметров орбит планет и подобрать скорость гравитации и скорость Солнечной системы так, чтобы расчетные значения смещений, полученные на корректной математической модели Солнечной системы, совпали с наблюдаемыми значениями.

 

 

Как я писал в предыдущей статье этого цикла статей [1], в середине 19-го века Леверье, при создание своей теории планет, обнаружил, что наблюдаемые значения положений планет отличаются от рассчитанных им по теории Ньютона (аналогичные отклонения обнаружил в 1895 году, при создании своей теории планет, и Ньюком). И после открытия Леверье все ученые кинулись создавать свои физические теории, объясняющие аномальные вековые смещения перигелиев планет, и совершенно забыли, что, согласно Ньюкому, и по другим параметрам орбит планет имеются аномальные отклонения. Некоторые из этих теорий я уже приводил в предыдущей статье (см. таблицу 1). Но все они (кроме астрономической теории Зеелингера), как я уже писал ранее, никак не объясняют аномальные смещения линии узлов, углов наклона орбит и их эксцентриситетов, которые, по приводившимся мною данным Ньюкома, тоже имеются. Это относится и к ОТО Эйнштейна, которая, не смотря на это, почему то считается верной теорией. Более того, для обоснования ее справедливости, как экспериментальное доказательство, приводят именно данные по вековым смещениям перигелиев планет, но умалчивают об имеющихся аномальных вековых смещениях других параметров орбит планет, которые она никак не объясняет.

 

 

Таблица 1. Аномальное смещение перигелия внутренних планет (в угл.сек за 100 лет).

 

 

Меркурий

Венера

Земля

Марс

Наблюдения  (Ньюком [6])

575,06

42,52

1162,92

1602,69

теория Ньютона (Ньюком [6])

533,82

49,85

1156,95

1594,65

Аномальный остаток (Ньюком [6])

41,24

-7,33

5,97

8,04

теория Холла (Ньюком [7])

43,37

16,98

10,45

5,55

теория Эйнштейна (Субботин [7])

43,03

8,62

3,83

1,35

теория Гербера (Хайдаров [3])

43,03

8,62

3,83

1,35

теория Ритца k=6,4 (Роузвер [9])

41,0

8,0

3,4

-

теория Маха (Зайцев [10])

43,0

23,0

17,0

11,0

теория Зеелингера (Роузвер [9])

41,3

7,3

4,2

6,3

 

 

 

Но, как будет показано мною ниже, для объяснения аномальных смещений  перигелиев планет можно было бы и не создавать сложные физические теории, а объяснить эти смещения с использованием феноменологической формулы Ньютона (1), но при этом грамотно учесть эффекты, вытекающие из конечности скорости распространения гравитации в какой то среде, что при разных скоростях всей Солнечной системы позволило бы получить самые разнообразные значения этих смещений. При этом от конечности скорости распространения гравитации у нас может возникнуть два эффекта. Первый это запаздывание потенциалов по координатам, а второй это что-то похожее на запаздывание потенциалов у Вебера. И хотя названия у этих эффектов одинаковые, но они совершенно разные. Первый эффект вызван действительно именно запаздыванием потенциалов, а во втором эффекте никакого запаздывания как такового нет, а учитывается динамическое взаимодействие движущегося потенциала поля с движущимся пробным зарядом (массой), т.е. как бы учитывается примитивное динамическое давление Лапласа, но возникающее не при корпускулярной природе гравитации, а уже на более высоком теоретическом уровне, т.е. с привлечением потенциалов поля. При этом первый эффект у нас наблюдается, когда поле  создается движущимся зарядом (массой), а пробный заряд (масса), на которое действует это поле, может быть и покоящимся, а второй эффект у нас проявляется при движении любого из зарядов (масс). Для простоты я эти эффекты называю запаздыванием по координатам и динамическим давлением. А, судя по тому, что термин «запаздывание по координатам», который я часто употреблял в дискуссиях на форумах, уже прижился и используется в общении на форумах, надеюсь, что приживется и термин «динамическое давление».

 

 

                                              Эффект запаздывания потенциалов по координатам

 

 

Основная идея в запаздывание потенциалов по координатам  состоит в том, что из-за конечности скорости гравитации одно тело воздействует на другое не из той точки, где оно находится в данный момент времени, а из той точки, где оно находилось несколько минут (секунд) назад, которые необходимы, чтобы гравитационное поле из той точки достигло приемника этого поля. Я эти точки, где источник поля находился несколько минут тому назад, называю следами планет. А на приведенном рисунке 1, когда у нас движутся и источник поля и приемник, например, Меркурий в конкретный момент времени t3 находиться в точке 1-3, а Венера в точке 2-3. Если пренебречь (для простоты изложения) разностью расстояний между точками 1-3 и 2-1 а также 2-3 и 1-1 и принять их равными L, то можно считать, что в данный момент времени t3 на Меркурий со стороны Венеры будет действовать сила притяжения, рассчитанная по формуле (1), из точки 2-1, т.е. из точки, где Венера находилась некоторое время тому назад dt= t3-t1= L/Vgr, где Vgr - это скорость распространения гравитации. Аналогично и на Венеру со стороны Меркурия будет действовать сила притяжения из точки 1-1, т.е. из точки, где находился Меркурий dt секунд тому назад.

 

 

 

Рис.1 Взаимодействие Меркурия и Венеры при конечной скорости гравитации демонстрирующее запаздывание потенциалов по координатам.

 

Ведь, если принять, что гравитация распространяется с конечной скоростью, то ясно, что, когда Меркурий будет находиться в точке 1-3, гравитационное поле, создаваемое Венерой, достигнет его только из точки 2-1, где Венера находилась dt секунд тому назад, за которые это поле преодолело расстояние L2-1 равное L и достигло Меркурия. А, когда Венера находилась в точке 2-2, ее гравитационное поле преодолело только расстояние L2-2 равное L/2 и еще не достигло Меркурия. Естественно, из точки 2-3 ее поле даже еще не успело вылететь, т.к. Венера только что прибыла в эту точку. И, следовательно, со стороны Венеры на Меркурий будет в данный момент времени действовать сила F2-1 из точки 2-1, а на Венеру со стороны Меркурия будет действовать сила F1-2 из точки 1-1. Таким образом, точки 2-1 и 1-1 у нас будут следами Венеры и Меркурия для данного момента времени при расчете сил притяжения между Венерой и Меркурием. То же самое мы с Вами наблюдали и рассматривая в предыдущей части этого цикла статей планетную аберрацию, т.е. аберрацию вызванную конечностью скорости распространения света и движением источника, т.е. можно сказать, что она тоже обусловлена запаздыванием по координатам от движения источника. А, проводя аналогию с эффектами аберрации, можно вполне предположить, что при взаимодействии двух тел у нас будет присутствовать и звездная аберрация, т.е. вызванная движением приемника поля, которая тоже будет немного изменять силу притяжения между двумя массами, т.к. планеты являются не точками, а имеют определенную протяженность, но рассматривать и этот эффект мы не будем. Только немного уточню, что данный эффект от звездной аберрации не имеет ничего общего с эффектами Лапласа и Ван Фландерна, хотя они для расчета тангенциальной силы и использовали аналогию со звездной аберрацией.

 

 

При этом, учитывая эффект запаздывания по координатам, надо четко понимать, что гравитационное поле, создаваемое какой то массой, будет присутствовать во всех точках пространства всегда, но вот величина напряженности этого поля при движении этой массы будет во всех точках пространства со временем изменяться. Поэтому для расчета силы притяжения пробной массы при взаимодействии ее с этим полем надо всегда определять место, из которого это гравитационное поле достигло пробной массы в данный момент времени. При этом силы F2-1 и F1-2 не только будут разными по модулю, но и направлены будут не по одной прямой. Да, при этом будет нарушаться 3-ий закон Ньютона, но он справедлив только при непосредственном контакте двух тел, а здесь у нас тела взаимодействуют не непосредственно, а через создаваемые ими поля, поэтому и не требуется соблюдение этого закона. Но сейчас в науке распространено и иное мнение, согласно которому при равномерном движении заряда (массы) его поле, распространяющееся на всю вселенную, является стационарным и являясь какой-то сплошной средой движется вместе с телом как приклеенное к нему и только при ускорении заряда (массы) происходящие при этом возмущения начинают распространяться от заряда (массы) в этой сплошной среде до самых окраин вселенной с какой то скоростью. Но, т.к. при этом на самом краю вселенной напряженность поля изменяется синхронно с движением массы, то все же получается, что здесь скорость распространения гравитации принимается равной бесконечности, т.е. ни о какой скорости здесь вообще не может быть и речи.

 

 

Я не знаю кто первым попытался применить потенциалы запаздывающие по координатам на практике, хотя, многие авторы пишут, что Лаплас, но это не правильно, т.к. он рассмотрел совсем другой эффект. Возможно, это был Гаусс, но наверняка я в этом не уверен, поэтому в дальнейшем я различные варианты описания взаимодействия между планетами буду именовать фамилиями авторов этих расчетов, например, расчет по Ньютону или по Веберу, а вот вариант расчета с эффектом запаздывания потенциалов по координатам я так и буду именовать или буду писать, что это потенциалы Юдина, когда буду рассматривать их комбинацию с динамическим давлением. Давайте теперь посмотрим, как мы можем реализовать этот эффект на практике, т.е. какие мы будем использовать формулы и для этого воспользуемся рисунком 2. Здесь у нас в момент времени  t в точке Р расположено пробное тело массой m1 (или заряд q1) и нам надо вычислить силу с которой это тело будет притягиваться массой m2 (или зарядом q2), которая находится сейчас в точке 2. При этом тело m2 движется со скоростью V вдоль оси X и текущее расстояние между двумя массами равно R. А, т.к. мы принимаем гипотезу о том, что любое взаимодействие между телами происходит с конечной скоростью, то, например, гравитационный потенциал массы m2 в момент времени t достигнет точки Р только из точки 2', где эта масса находилась в момент времени t', который определится как (2), где Vgr скорость распространения гравитации. И таким образом и потенциал в точке Р в момент времени t создаваемый массой m2 определится как (3), а сила притяжения массы m1 к следу тела m2 (точка 2'), т.е. F2'1 определится по формуле (4), где E2' напряженность поля создаваемая в точке Р телом m2 при его нахождении в точке 2'.

 

 

Рис. 2. Схема к расчету потенциалов запаздывающих по координатам.

 

t'= t - ∆t = t - R' / Vgr                                                              (2)

 

φ = G*m2 / R’                                                                         (3)

 

F2'1= m1*E2'= m1*m2 / R'^2                                                (4)

 

 

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

 

R’ ≈ R + V*R / VgrR + Vr * ∆t                            (5)

 

то становится ясно, что время изменения радиуса ∆t=R/Vgr дает приблизительный результат, который отличается от точного ∆t=R'/Vgr и радиальная скорость изменения радиуса Vr=V*cos(a) тоже будет справедлива только вблизи точки 2, а вблизи точки 2' у нас будет Vr=V*cos(a') и то с учетом того, что у нас будет равномерное движение, т.е. будет V'=V. Можно, конечно, получить и точное решение, если масса 2 движется прямолинейно и равномерно. Для этого надо в уравнение (2) подставить R', который определится по теореме Пифагора из условия, что координаты точки Р заданы, а абсцисса точки 2' будет V*t'. Получаем квадратное уравнение и находим его два корня для ∆t, а потом оставляем одно решение. Да, в своей программе Potencial2 я так и делаю, рассматривая прямолинейное и равномерное движение, но для решаемой нами задачи (движение планет Солнечной системы) возможен только вариант расчета, который напоминает расчет по формуле (5), но при этом расчеты ведутся за несколько итераций.

 

 

Конкретно в своей программе Solsys7 я делал две итераций для расчета координат следов планет, что при существующих скоростях планет в Солнечной системе дает почти точное решение. При первой итерации я (согласно рис.2, где у нас масса 2 движется равномерно вдоль оси X) по известному времени ∆t1=R/Vgr находил новую (запаздывающую) координату массы m2 как X2(1)= V*t - V*∆t1, а затем по новым координатам вычислял первое приближение расстояния R1 и вычислял новое время ∆t2=R1/Vgr для второй итерации, при которой опять находил запаздывающую координату X2(2)= V*t - V*∆t2, по которой находил новое расстояние R2 и уже по нему потом определял силу притяжения между массой 1 и следом массы 2. А конкретно в программе Solsys7, запаздывающие координаты планет, т.е. координаты из которых они через созданное ими в прошлом времени гравитационное поле будут воздействовать на планету, ускорения которой мы определяем в дифференциальных уравнениях по силам, действующим на планеты, я рассчитывал так (смотрите нижеприведенный код - это один из вариантов расчета - простейший, который мною назван DX1). Сначала я определял время T(i, j) необходимое гравитации, чтобы долететь от планеты i до планеты j, на которую оно будет воздействовать, по текущим координатам планет, а потом откатывал планету i немного назад пропорционально времени запаздывания и ее скоростям по осям координат (за две итерации). И уже из этой точки (следа планеты) вычислял силу притяжения планетой i планеты j. Причем здесь дан алгоритм расчета запаздывающих координат планеты воздействующей своим полем на исследуемую планету не абсолютных, а сразу относительных. При этом первое значение времени запаздывания T(i, j) определялось в программе до этого расчета по текущим координатам планет.

 

 

For k = 1 To 2 '                                      количество итераций равно 2

For i = 0 To Ne '  NE=10

For j = 0 To Ne

If i = j Then GoTo 300 '                        расчет координат следов планет

DX(i, j) = X(j) - X(i) + VX(i) * T(i, j) 'расстояние между следом i-го тела воздействующим на j-е тело по оси X

DY(i, j) = Y(j) - Y(i) + VY(i) * T(i, j)

DZ(i, j) = Z(j) - Z(i) + VZ(i) * T(i, j)

Rxy(i, j) = Sqr(DX(i, j) * DX(i, j) + DY(i, j) * DY(i, j)) ' расстояния в плоскости XY между следом i-го тела воздействующим на j-е тело

R(i, j) = Sqr(Rxy(i, j) * Rxy(i, j) + DZ(i, j) * DZ(i, j)) '    общее расстояния между следом i-го тела воздействующим на j-е тело

T(i, j) = R(i, j) / Vgr '                                                          время необходимое для распространения гравитации

300: Next j

Next i

Next k

 

 

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

 

 

А, чтобы визуально увидеть погрешность от не точного определения запаздывающих координат по приведенному выше коду программы нам придется задать очень маленькую скорость гравитации, т.к. при скорости гравитации равной или больше скорости света эта погрешность на рисунке будет не видна. Поэтому на рис.3 я привожу расчет запаздывающих координат для Меркурия и Венеры, который можно выполнить в программе Solsys7 на форме для тестового расчета влияния скорости гравитации на силы притяжения между планетами, при скорости гравитации равной одной тысячной скорости света. На рисунке до первого приближения силы притяжения одинаковые F0, а после первого и второго приближений  (F21>F12) и кружками показаны положения планет в текущий момент времени и их следы при первом и втором приближении. Кому-то может показаться, что при таком расчете запаздывающих координат, когда учитывается только прямолинейное и равномерное движение (в программе это вариант расчета DX1), будет большая погрешность в их определении, но это не так. При скорости гравитации равной или больше скорости света запаздывание будет такими маленькими, что результат будет практически таким же, как и в варианте расчета, когда я учитываю еще и ускорения планет (в программе это вариант расчета DX2).

 

 

Рис.3 Расчет запаздывающих потенциалов, т.е. сил притяжения планет следами других планет после двух приближений расчета времени запаздывания потенциала. Скриншот программы Solsys7.

 

 

 

А сейчас я покажу на экспериментальных данных, полученных при проведении вычислительных экспериментов на математической модели Солнечной системы, что в равномерно и прямолинейно движущейся ИСО, при учете запаздывания по координатам, получаются совсем другие смещения параметров орбит планет, чем в покоящейся системе. При этом в модели скорость распространения гравитации будет задаваться как скорость распространения в какой то среде, которую можно назвать хоть эфиром, как это было в классике, хоть физическим вакуумом, как это сейчас в ОТО, т.к. для нас будет важно только то, что эта среда есть и находится в абсолютном покое, а уже в ней распространяется гравитация и движутся планеты и Солнце. А то, что любое взаимодействие между телами распространяется с конечной скоростью (в теории относительности даже есть ограничение на максимальное значение этой скорости), сейчас не оспаривает никто и, следовательно, эффект запаздывания потенциалов по координатам обязательно должен учитываться в модели Солнечной системы. Да, без компьютеров решить эту задачу ученые, которые работали над этой проблемой 100 лет назад, не могли, но ведь компьютеры уже пол века как в употреблении, но людей так оболванили СТО и ОТО, что никому и в голову не приходило проверить практически - а действительно ли в инерциальных системах все процессы протекают одинаково. А теперь получается, что мне, полученные на моей модели Солнечной системы смещения параметров орбит планет в движущейся Солнечной системе с учетом запаздывания потенциалов по координатам, и сравнить то не с чем. Нет таких данных. Поэтому, приведу только свои данные по смещениям параметров орбит планет, получающиеся при учете запаздывания по координатам потенциалов планет, которые у меня получились при использовании программы Solsys7m.

 

 

Вычислительные эксперименты проводились мною при задании начальных параметров орбит планет и их долготы в плоскости фиксированной эклиптики, т.е. в эпохе J2000, на 12 часов 1 января 1801 года, по моей теории планет Ser0. Затем по стандартным формулам определялись начальные координаты планет и их скорости при условии, что эти параметры орбит получены при движении планет относительно неподвижного Солнца, находящегося в начале абсолютно неподвижной системы координат. А потом начальные координаты и скорости Солнца с учетом воздействия планет не уточнялись, хотя вообще-то желательно бы было это сделать (программа позволяет). После этого я задавал произвольную скорость всей Солнечной системе (от -300 до +300 км/с) по одной из осей координат (по остальным осям скорость принималась равной нулю), т.е. увеличивал на эту величину начальные скорости всех планет и Солнца и начинал вычислительный эксперимент, т.е. моделировал движение объектов Солнечной системы. Во время эксперимента я определял до 31 декабря 2000 года, при каждом обороте планеты вокруг Солнца, значения получающихся значений долготы перигелия и восходящего узла, угла наклона орбиты и эксцентриситета и записывал их в файл. Затем, полученный массив данных, я подвергал статистической обработке и находил вековые смещения параметров орбит планет. Подробности о примененной мною методике статистической обработке этих данных смотрите в предыдущей статье [1].

 

 

 

При этом решение системы дифференциальных уравнений выполнялось с основным шагом 3600 секунд, а вблизи перигелия или узла восхождения с уменьшенным шагом решения. Для Меркурия этот уменьшенный шаг составлял 3,6 сек, для Венеры 7,2 сек, для Земли 18 сек и для Марса 36 сек. А область, где уменьшенный шаг решения применялся, определялась в интервале +/-1 градус от полученного значения перигелия или узла восхождения при предыдущем обороте планеты. При этом т.к. у Венеры наблюдается очень большой разброс в мгновенных значениях долготы перигелия, вызванный очень маленьким эксцентриситетом орбиты, когда воздействия от других планет могут значительно сместить его мгновенное значение, я для нее определял перигелий и узел восхождения в увеличенном интервале, т.е. +/-2 градуса от предыдущих значений.

 

 

А в связи с тем, что при проведении поисковых экспериментов по выявлению влияния скорости гравитации на смещения параметров орбит планет, я даже ориентировочно не знал какова должна быть эта скорость (данные Лапласа и Фландерна использовать нельзя), то я решил начать с одной скорости света (10^0) и ограничиться значением 1000 скоростей света (10^3). Полученные при этом мною значения по вековым смещениям параметров орбит четырех планет (от Меркурия до Марса, т.к. для других планет вычислять их не имеет пока смысла) я привожу в таблицах 1-1-X…4-6-Z. При этом в таблицах в первой колонке указаны показатели степени у 10 для задания скорости гравитации (10^0…10^3), а звездочками отмечены вековые смещения, которые рассчитаны только приблизительно из-за нелинейностей графиков смещений. А при очень больших нелинейностях, вызванных малой скоростью гравитации, например, когда 100 лет смещение уменьшается, а потом 100 лет увеличивается или изменяется очень не устойчиво, я вообще не привожу полученные значения. Данные по влиянию скорости гравитации и скорости Солнечной системы на все параметры орбит для четырех внутренних планет я привожу в приложение 1. Здесь же я, как пример, могу привести данные по смещению перигелия Меркурия при изменении скорости Солнечной системы по оси Y и смещения перигелия Марса при изменении скорости Солнечной системы по оси Z.

 

 

Таблица 1-1-Y. AlfaP1-VYsys  - вековые смещения перигелия Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-15612,46*

-9466,43

-4237,82

529,23

4990,88

9117,31

13154,15

1

-859,33

-392,79

+69,79

529,63

986,07

1438,82

1966,32

2

+392,20

+438,03

483,90

529,68

575,67

621,36

667,33

3

+515,89

+520,55

525,11

529,69

534,26

538,83

543,48

 

 

Таблица 4-1-Z. AlfaP4-VZsys  - вековые смещения перигелия Марса (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

1694,43

1653,10

1617,09

1587,93

1565,70

1551,53

1546,45

1

1595,85

1593,11

1590,58

1587,98

1585,45

1582,94

1580,57

2

1588,85

1588,59

1588,17

1587,95

1587,71

1587,45

1587,15

3

1588,03

1587,95

1587,94

1587,99

1587,95

1587,99

1587,85

 

 

Как видим, в первом случае (Меркурий) смещение перигелия очень сильно зависит от изменения, как скорости гравитации, так и скорости Солнечной системы, а во втором случае (Марс) положение перигелия остается почти неизменным. Таким образом, влияние, как скорости гравитации, так и скорости Солнечной системы, на различные параметры орбит различных планет существенно различно, но однозначно это влияние на параметры орбиты Меркурия является самым заметным. И при скорости гравитации равной скорости света у нас даже при небольшой скорости Солнечной системы получается смещение перигелия Меркурия в тысячи угловых секунд, а не какие-то жалкие 43 угловые секунды аномального смещения сверх классических 529, которые якобы объясняет ОТО и где скорость гравитации принята равной скорости света. Таким образом, мы ясно видим, что скорость гравитации при распространении ее в какой то среде должна быть гораздо больше скорости света, чтобы не было таких больших смещений. При этом, если смещение таких параметров орбит планет как перигелия, узла восхождения и угла наклона непосредственно не наблюдаются, а являются вычисляемыми при обработке данных наблюдений параметрами и тут сторонники ОТО всегда могут заявить, что мы не правильно обработали данные, то вот такие параметры как период обращения, эксцентриситет и большая полуось можно получить и из непосредственных наблюдений и сравнить их с данными древних астрономов, которые в таких простейших вычислениях уж точно не допустили ошибки. И ниже я привожу таблицы изменения большой полуоси Меркурия и его угловой скорости обращения вокруг Солнца при различной скорости Солнечной системы по оси Y.

 

 

Таблица 1-5-Y. Rsr1-VYsys  - вековые смещения большой полуоси Меркурия (тыс. км)

 

n

-300

-200

-100

0

100

200

300

0

-9287,97*

-5579,17*

-2545,99

1,16

2159,49

3995,73

5520,32*

1

-720,01

-476,02

-235,99

0,12

232,38

460,89

685,30

2

-70,42

-46,90

-23,42

0,01

23,40

46,78

70,10

3

-7,02

-4,68

-2,34

0,00

2,34

4,70

7,02

 

Таблица 1-6-Y. Wsr1-VYsys  - вековые смещения орбитальной угловой скорости Меркурия (рад/век)

 

n

-300

-200

-100

0

100

200

300

0

-

459,0944*

187,3206

-0,0781

136,3162

-238,5178*

-314,5845*

1

49,8034

32,6638

16,0672

-0,0077

-15,5838

-30,6835

-45,2992

2

4,7691

3,1739

1,5840

-0,0007

-1,5799

-3,1557

-4,7260

3

0,4741

0,3166

0,1579

-0.0000

-0,1581

-0,3174

-0,4740

 

 

По крайней мере, за последние 2000 лет у нас имеются более-менее точные значения этих параметров и они говорят о том, что эти параметры изменились очень незначительно. А по данным вычислительных экспериментов получается, что большая полуось, например, при скорости Солнечной системы 200 км/с увеличивается на 4 млн. км за сто лет и за 20-ть веков она должна была увеличиться на 80 млн. км. Но большая полуось Меркурия и сейчас составляет только 58 млн. км и, следовательно, скорость гравитации никак не может быть равна скорости света. А, если учесть то, что уже во времена Гиппарха и Птолемея период обращения Земли вокруг Солнца был известен очень точно и споры шли об одной сотой дня, то и данные по угловой скорости обращения Земли вокруг Солнца тоже говорят о том, что скорость гравитации должна быть, как минимум, на несколько порядков больше скорости света. Ведь, согласно нижеприведенной таблице, угловая скорость обращения Земли вокруг Солнца при скорости Солнечной системы по оси X, т.е. VXsys=200 км/с, только за 100 лет должна при этом увеличиться на 6,4 рад/век, а за 20-ть веков это будет 20-ть оборотов, т.е. одна пятая от ста оборотов, т.е. за век, но по данным наблюдений период обращения Земли за это время остался неизменным (в пределах точности измерений).

 

 

Таблица 3-6-X. Wsr3-VXsys  - вековые смещения орбитальной угловой скорости Земли (рад/век)

 

n

-300

-200

-100

0

100

200

300

0

-9,7175

-6,4737

-3,2339

-0,0016

3,2231

6,4381

9,6422

1

-0,9688

-0,6456

-0,3229

-0,0001

0,3227

0,6454

0,9680

2

-0,0968

-0,0646

-0,0323

+0,0000

0,0323

0,0646

0,0969

3

-0,0097

-0,0064

-0,0032

+0,0000

0,0033

0,0065

0,0097

 

 

И по полученным данным можно сделать уже сейчас вполне определенные выводы и один из них о том, что инерциальных систем отсчета, т.е. систем движущихся прямолинейно и равномерно, где все явления описываются точно так же, как и в покоящейся системе, не существует. Ведь как бы ни была велика скорость гравитации, но все равно у нас получатся данные по смещениям параметров орбит планет хоть на немного, но отличные от данных полученных в абсолютно неподвижной системе отсчета. Таким образом, вся СТО основанная на замене одной ИСО другой ИСО является ошибочной теорией. Точно также, и ОТО является ошибочной теорией, т.к. в ней принята скорость гравитации равная скорости света, а это явно не так. К тому же ОТО, хоть и декларирует конечность скорости гравитации, но совершенно не учитывает запаздывание по координатам, т.е. фактически использует бесконечную скорость гравитации. Более того, как показано в моей работе [1] и как я покажу ниже, ОТО позволяет учесть только аномальные смещения перигелиев планет в покоящейся Солнечной системе, а на остальные параметры орбит планет она совершенно не влияет, но по приведенным выше данным мы видим, что, если наша Солнечная система движется, то уж точно надо объяснять и отклонения других параметров орбит от чисто Ньютоновского расчета.

 

 

Причем, приведенные мною тут данные о различном поведении систем при разной поступательной скорости движения этих систем не являются единственными. Например, в натурных экспериментах, проведенных с вращающимися маятниками [39], у автора тоже получается, что в движущихся поступательно и равномерно системах вращающиеся маятники ведут себя по-разному. Но, даже при наличии этих экспериментальных данных найти скорость гравитации будет не просто, т.к. нам не известна даже абсолютная скорость Солнечной системы, хотя ориентировочные данные есть VXsys= -358 км/с,  VYsys= 57 км/с, VZsys= -73 км/с, но особо доверять им не стоит. Итого мы имеем задачу с четырьмя неизвестными, если не считать того, что я не уверен и какова должна быть механическая динамическая математическая модель Солнечной системы, хотя в том, что она должна учитывать запаздывание по координатам, я уверен.

 

 

                                             Эффект динамического давления поля

 

 

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

 

 

Сейчас считается, что Пауль Гербер [2] объяснил аномальное вековое смещение перигелия Меркурия именно запаздыванием потенциалов, вызванным конечностью скорости гравитации, и здесь Гербер попытался учесть не только запаздывание потенциалов по координатам, рассмотренное нами ранее, но и   динамическое давление. А сама идея запаздывания потенциала при движении заряженных частиц, была высказана Гауссом незадолго до смерти в его письме Веберу в 1835 году, которое его друг Вебер скрыл от научной общественности, а потом изложил эту идею как свою, но в искаженном виде, хотя позже работа Гаусса все же была опубликована. А сейчас аналогичную теорию запаздывания потенциалов в виде давления эфира использует Карим Хайдаров [3], который тоже объясняет с ее помощью аномальные отклонения перигелиев планет. Но началось все официально с Вебера, который предложил исходя из закона Ампера свою формулу (6-1), если ее записать не для зарядов, как она была первоначально опубликована, а для масс, а Гербер в 1898 году получил уравнение (6-2), которое подробно расписал в 1902 году.

 

Fw = - (G*M*m/R^2) * (1 – 0,5*Vr^2/c^2 + R*Ar/c^2)     (6-1)

 

Fg = - (G*M*m/R^2) * (1 – 3*Vr^2/c^2 + 6*R*Ar /c^2)     (6-2)

 

В принципе ничего нового в этой формуле Гербера нет, т.к. аналогичные формулы предлагали и Тессеран (1872) и Ньюком (1882), и отличается она у Гербера только коэффициентами при втором и третьем членах для выражения потенциала у Вебера. Но вот идеи, заложенные Гербером, при выводе его формулы, принципиально отличаются от идей Вебера, т.к. у Вебера в формуле величина c - это всего-навсего коэффициент перехода от электростатического взаимодействия к электродинамическому, т.е. ни о каком запаздывании потенциалов он не говорил и скорость взаимодействия он принимал равной бесконечности. Это уже потом, когда ученые ознакомились с работой Гаусса и нашли, что коэффициент перехода Вебера равен корню квадратному из скорости света, то его формула тоже стала восприниматься как учет в ней запаздывания потенциала, поэтому формулу (6-1) я записал так, как это будет при коэффициенте с равном скорости света.

 

 

В этих формулах, т.к. они достаточно большие, а нам, при рассмотрении вывода Гербера, придется часто их записывать, я ввел свои обозначения, которые буду использовать и при рассмотрении далее других вопросов. Скорость изменения расстояния R между двумя массами будет Vr = dR/dt, а ускорение, с которым это расстояние будет изменяться, будет Ar = dVr/dt = d2R/dt2. Таким образом, как видите, идеи СТО заложены в этих формулах задолго до создания СТО, т.к. и скорости и ускорения здесь используются относительные, т.е. одной массы относительно другой. В этих формулах (6-1) и (6-2), как сейчас объясняют, второй член, зависящий от скорости сближения или удаления двух тел Vr, отвечает за динамическое давление гравитационного поля, а третий член, зависящий от ускорения с которым изменяется расстояние между двумя гравитирующими массами Ar, отвечает за образование при этом гравитационных волн, т.е. за излучение. Вот только с такими относительными величинами получаются абсурдные результаты. Если обе скорости будут равны, то не будет запаздывания по координатам, а, если оба ускорения равны, то не будет никакого излучения. А в том случае, если основной заряд, создающий поле, движется без ускорения, а с ускорением движется пробный заряд, то в этом случае должен будет излучать основной заряд.

 

 

А говорю я, рассматривая динамическое давление, чаще всего именно о теории Гербера по тому, что он и весь вывод строит из идеи запаздывания потенциала и ему удалось получить конкретное решение задачи, а именно формулу для расчета аномального смещения перигелия Меркурия, которая давала необходимые 43 угл.сек. При этом формула Вебера, как пишет [9], дает только 7,2 угл.сек, хотя, по выполненным ранее расчетам Тиссерана, считалось, что она дает 14 угл.сек. Правда Гербер получил формулу для определения скорости гравитации по заданному смещению перигелия Меркурия и эта скорость у него при этом получилась равной скорости света, но из этой формулы легко определить и смещение перигелия при заданной скорости гравитации (7), где с- скорость гравитации, равная скорости света, а- большая полуось эллипса, е- эксцентриситет, а Т- период обращения планеты вокруг Солнца. Сейчас считается, что эта формула (7) получена Эйнштейном из его ОТО, но как видим она была получена до Эйнштейна и без всякой ОТО.

 

dAlfaP= 24*pi^3*a^2 / [c^2*(1-e^2)*T^2]                            (7)

 

 

Считается так же, что эффект смещения перигелия в теории Гербера получается из-за того, что планета при движении вокруг Солнца по эллипсу то удаляется от Солнца (двигаясь к афелию), то приближается (двигаясь к перигелию) и вследствие этого на планету половину пути действует сила меньшая, чем положено по формуле Ньютона, а половину пути большая, что вызвано или уменьшением гравитационного давления на планету или увеличением. И во всех работах пишут, что поворот перигелия объясняется воздействием, вызванным именно вторым членом в формулах (6-1) и (6-2), но это не так и смещение перигелия обеспечивает третий член в этих формулах, который отвечает у авторов всех этих работ за излучение. И это легко проверить с помощью моей программы Solsys7m. Выбираем на форме 23 расчет динамического давление по формуле Гербера, задаем значение k1=3, k2=0 и рассчитываем параметры орбиты Меркурия за период с 1801 по 2000 год. Обрабатываем данные и находим, что смещение перигелия получилось равное нулю. Теперь задаем k1=0, k2=6 и получаем 43 угл.сек. Естественно, тот же эффект будет и по формуле Вебера.

 

 

 

Рис.4. Скриншот программы Solsys7mmm с вариантами расчета динамического давления гравитации и запаздывания по координатам.

 

 

А вот как так получилось, что у Вебера, Гербера и других авторов при учете вроде бы одного и того же эффекта получились разные формулы – это интересно. Но вопрос этот не такой простой, т.е. одной математикой мы здесь не отделаемся, и надо будет рассматривать представления ученых того времени о распространении взаимодействия в пространстве. Более того, в этих формулах получается, что сила зависит от силы или иными словами ускорение зависит от ускорения, т.к. сила и ускорение это суть одно и тоже и, здесь тоже надо будет рассматривать методологические подходы к описанию процессов протекающих в Природе. Вот, например, Герц [19] создал всю свою механику без единого упоминания о силе, т.е. используя в своих уравнениях только ускорения, которые ему с успехом заменяют силы.

 

 

Интересно в этой связи заметить, что и в формуле Гаусса ускорения не было и зависимость ускорения (силы) от ускорения (силы) это уже изобретение именно Вебера. К сожалению, работа Гаусса, после того, как она все же была опубликована, осталась как-то в стороне от бурных событий того времени и о ней быстро забыли. А Максвелл в своем трактате [27] прямо заявил, что формула Вебера правильная, а Гаусса не правильная. Но и сейчас некоторые ученые возвращаются к этому вопросу о связи силы и ускорения. И, например, автор работы [11] прямо заявляет «Итак, сила воздействия на тело есть не что иное, как его ускорение, но выраженное в других единицах. Поэтому, естественно, сила воздействия на тело не может зависеть ни от ускорения тела, ни от производных ускорения по времени». А, если предположить, что в формуле Вебера ускорение действительно отвечает за излучение и все с этим соглашались, то мне не понятно, например, почему же при рассмотрении простейших колебаний груза на пружинке никто не учитывал излучения груза, т.е. дополнительной силы, вызванной его ускорением, и экспериментальные данные при этом совпадали с расчетными.

 

 

Так откуда же все-таки взялась эта глупость с ускорениями в уравнениях Вебера и Гербера и почему ее не было у Гаусса. Чтобы выяснить, где конкретно у Вебера и Гербера появились ускорения, надо посмотреть их вывод этих уравнений, но мне это сделать практически не возможно, т.к. я нигде не могу найти работ Гаусса и Вебера по этому вопросу. Все авторы, которые касаются рассмотрения этого вопроса, ссылаются на книгу Роузвера [9] и трактат Максвелла [27], поэтому и я буду ссылаться только на тот материал, который есть у них, но про вывод уравнений у них, к сожалению, практически ничего нет. А вот вывод Гербера я рассмотрю подробно, т.к. его работа имеется у меня в распоряжении, да и Роузвер частично рассматривает его вывод. Но, прежде, чем мы приступим к рассмотрению этого вопроса будет очень полезно ознакомиться с механикой Герца, хотя бы в общих чертах, чтобы понимать как мыслили ученые того времени и какие у них были проблемы. Обосновывая свой подход к рассмотрению вопросов механики, Герц дает и обзор двух существовавших до него подходов. Ниже я перечисляю эти подходы с основной их характеристикой и с указанием их основных признаков, а также сразу привожу и предлагаемый мною 4-ый подход и добавлю нулевой подход Ньютона, который Герц почему то пропустил, т.к. отождествил его с силовым подходом Эйлера.

 

Нулевой – импульсный (пространство – время – масса – импульс)      - Ньютон

Первый – силовой (пространство – время – масса – сила)                      - Эйлер

Второй – энергетический (пространство – время – масса – энергия)     - Лагранж

Третий – кинематический (пространство – время – масса)                      - Герц

Четвертый – мощностной (пространство – время – масса – мощность) - Юдин

 

 

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

 

 

А Лагранж добавил в своем подходе (в конце 18 века, начале 19-го) к подходу Эйлера принцип виртуальных скоростей, но применил его не правильно, т.к. использовал при этом принципы вариационного исчисления, и поэтому получился принцип виртуальных перемещений, который и дал энергетический подход, хотя Лагранж этого и не желал и даже не знал не только о законе сохранения энергии, но даже слова энергия тогда не было. Ну, а потом, когда он связался с принципом наименьшего действия, то его энергетический подход превратился в энергетическо-геометрический и сейчас используется только в таком виде. Герц же в своем подходе объединил движение тел по инерции и дифференциальный принцип наименьшего принуждения Гаусса, а я в своем подходе к движению тел по инерции добавил принцип Даламбера в своем мощностном прочтении.

 

 

Вот так, вроде бы из одних и тех же теоретических предпосылок получилось пять различных подходов для описания явлений Природы. Но четкое и однозначное понимание всех этих подходов мы сможем получить только после того, когда мы дадим точные и однозначные формулировки пространству, времени, массе, силе и энергии, а также сможем перевести словесное описание этих подходов на современный язык математики, т.е. язык формул или, иначе говоря, создадим «аналитическую механику». Например,  Ньютон, по сути, только собрал в одном месте, т.е. в своих Началах [30] словесные формулировки уже известных законов (кроме третьего) и попытался их однозначно идентифицировать, но опять-таки словесно и у него это не очень удачно получилось. Первым, кто смог это сделать был Эйлер, который не только создал силовой подход, но и дал четкие определения всем используемым им величинам, а также изложил свой подход на языке математики. А до Эйлера все задачи решались геометрическими методами с пространными словесными доказательствами, как, например, рассмотренный мною в работе [28] вывод Бернулли для брахистохроны или решения задач Ньютоном в его Началах.

 

 

Здесь необходимо заметить, что официальная наука настойчиво вдалбливает в головы ученым, что «аналитическая механика» это именно геометрическая механика, т.е. аналитическая механика с энергетическо-геометрическим подходом Лагранжа, который получился у него при использовании вариационного исчисления. При этом, все остальные «механики» (Ньютона, Эйлера, Герца) это тоже «механики», но второго сорта, хотя у Герца тоже используется вариационный принцип, но не интегральный, как у Лагранжа, а дифференциальный. Но, как я показал в своей статье [8], именно механика Лагранжа является механикой 2-го сорта, и ее использовать можно только для решения учебных задач. Я бы даже назвал ее механикой 3-го сорта, т.к. даже аналитическая механика Герца тянет только на 2-ой сорт, а механика Ньютона вообще не аналитическая. И, если мы хотим говорить именно об аналитической механике, то всегда надо указывать о какой именно, например, об «аналитической силовой механике», т.е. о «механике Эйлера» и т.д. А такая низкая оценка «аналитической энергетическо-геометрической механике», т.е. «механике Лагранжа» связана с тем, что вариационное исчисление создавалось Эйлером и Лагранжем для решения именно геометрических задач на максимум и минимум, хотя уже в древней Греции такие простейшие задачи решали и без вариационного исчисления, а применено оно было Лагранжем для задач динамики, чего делать нельзя (подробнее об этом читайте в статье [8]). К тому же, ни одной сложной, т.е. реальной, задачи с использованием  вариационного исчисления решить нельзя и сейчас, например, в задачах управления, хотя для этих целей существует уже и динамическое программирование и принцип максимума Понтрягина.

 

 

Таким образом, различных подходов для описания явлений Природы существует много, но, в любом случае мой подход и, как его частный случай, подход Эйлера являются все же самыми продуктивными, понятными и корректными. Но при этом я пока не буду утверждать, что мой подход поможет разрешить все проблемы, например, существующие в электродинамике или поможет нам учесть динамическое давление гравитационного поля. А вот в связи с наличием разных подходов у меня возникает вопрос. Почему официальная наука нам преподносит законы Ньютона в силовой интерпретации Эйлера, а сама пропагандирует энергетическо-геометрический подход рассмотрения законов Природы Лагранжа-Гамильтона и при этом ссылается на авторитет Ньютона, хотя мы ясно видим, что у Ньютона просматривается все-таки подход, который ближе к мощностному. При этом  фраза «согласно законам Ньютона» стала просто штампом, т.к. при этом подразумеваются законы Ньютона в интерпретации силового подхода Эйлера, т.е. уже совсем другие законы. Но многие люди, считающие себя серьезными учеными, даже не догадываются об этом, считая, что энергетическо-геометрический подход Лагранжа-Гамильтона является просто логическим развитием подхода Ньютона, а других подходов для описания явлений Природы вообще не существует. Кого интересуют  подробности различных подходов описания явлений Природы я отсылаю к своей статье [67], а мы продолжим рассмотрение вопроса запаздывания потенциалов Вебера и Гербера.

 

 

Но, прежде чем приступить к рассмотрению этого вопроса, необходимо сделать еще и такое уточнение. Сейчас у нас потенциалом называется удельная потенциальная энергия, т.е. отношение потенциальной энергии пробного тела, помещенного в наше поле, к массе или заряду этого пробного тела. Но, ученые 19-го и начала 20-го века под потенциалами понимали ту же самую потенциальную энергию, но когда говорили о ней при передаче на расстояние и за конечное время. Вот, например, открываем трактат Максвелла [27] и читаем на стр. 373 «Это возражение не применимо по отношению к формуле Вебера, ибо им было показано, что если принять в качестве потенциальной энергии системы, состоящей из двух электрических частиц, величину φ= e*e’/r *[1 – Vr^2/2/c^2] …..». А уже на стр. 379 мы читаем «С другой стороны, электрическая частица посылает потенциал, величина которого  e*e’/r  зависит не только от заряда e излучающей частицы, но также от заряда e’ принимающей частицы и от расстояния r между частицами в момент испускания». Т.е. ни о какой удельной потенциальной энергии тут нет и речи, поэтому, говоря о потенциалах Вебера или Гербера, мы должны ясно понимать, что речь там идет именно о потенциальной энергии, а не о современных понятиях потенциалов. Не понятно только почему у Вебера, как мы это видим у Максвелла, потенциальная энергия получается положительной при отталкивании частиц, а, как мы увидим дальше, в потенциале самого Вебера она, почему-то, получается положительной и при рассмотрении взаимодействия двух притягивающихся масс.

 

 

Конечно же, было бы интересно рассмотреть подход Вебера к определению им запаздывающих потенциалов, но, как я писал, в моем распоряжении нет его двух работ по этому вопросу, а к тому же свою формулу (8) Вебер получил уже после того, как его раскритиковали за его формулу (6-1) потому, что он не понял того, что своровал у Гаусса. Таким образом, велика вероятность, что он, как и все математико-физики, просто подогнал решение под готовый ответ, т.е. под свою формулу (6-1) и, быстрее всего, ничего достойного вниманию там нет. Поэтому, давайте рассмотрим вывод формулы для запаздывающих потенциалов (9) и для силы взаимодействия между массами (6-2) при конечной скорости гравитации у Гербера. Я считаю, что этот разбор будет очень полезен в плане рассмотрения взглядов ученых того времени на природу вещей и раскрывает их методы познания Природы. К тому же, очень многие ученые писали, что формула Гербера (6-2) не правильная, а формула Вебера (6-1) правильная, но за все сто лет так никто и не удосужился найти в выводе Гербера хоть одну ошибку. Вот и Роузвер [9] пишет только, что «Вывод его формулы довольно запутанный, но приведенные ниже выкладки дают требуемый результат». Вот давайте и рассмотрим его выкладки на предмет наличия в них ошибок.

 

Pw= G*M*m * [1 - Vr^2/(2*c^2)] / R            (8)

 

Pg = G*M*m / [R*(1 – Vr/c)^2]                     (9)

 

 

Сразу надо заметить, что общий подход Гербера был верным и, по задумке, формула должна была учитывать не только динамическое давление, но и запаздывание потенциала по координатам. И, если бы Гербер более четко изложил свой вывод (а может быть, и сам понимал, что он делает) и не сделал в нем логические, методологические и математические ошибки, то цены бы его работе не было. Итак, свой вывод Гербер начинает с учета запаздывания потенциала по координатам, но при этом ни на кого не ссылается, т.е. такая точка зрения тогда считалась общепринятой. Но вывод у него очень, очень, очень запутанный, поэтому я дам свою интерпретацию его вывода и после каждого этапа буду сравнивать, что у него должно было получиться и что получилось. Итак, Гербер рассматривает две массы – притягивающую M и притягиваемую m, которые движутся с какими то скоростями так, что расстояние между ними R постоянно изменяется. А, т.к. потенциалу, чтобы передаться от массы M к притягиваемой массе m нужно какое то время ∆t, то, следовательно, он дойдет от массы M до массы m не из того положения, где сейчас находиться масса M (на рис.5 положение 2), а из того, где она была ∆t минут тому назад (положение 2’). И тогда, согласно Герберу, надо в расчетах учитывать потенциал, который определится по формуле

 

Pg1 = G*M*m / (R - R)                                   (10)

 

 

 

 

Рис.5 Пояснение к расчету запаздывания по координатам у Гербера.

 

 

Изменение расстояния ∆R у Гербера положительно при увеличении расстояния между массами. Следовательно, если скорости этих масс у нас будут направлены так, как показано на рис. 5 то это расстояние ∆R за время ∆t определиться как ∆R= ∆R2= Vr2*∆t и нас совершенно не интересует где в прошлом была масса m, т.е. нам не надо учитывать ∆R1= Vr1*∆t, но Гербер здесь допускает ошибку и использует значение ∆R= ∆R2 - ∆R1. А в таком случае при равных скоростях масс ∆R= 0 и у него не будет никакого запаздывания потенциала по координатам. Но, Гербер этого не замечает, а вот проблему определения ∆t он не мог не заметить, поэтому он очень долго и путано заметает следы пока не обнародует, что все таки ∆t=R/c, где c – это у него скорость распространения потенциала. Но, как мы выяснили при рассмотрении потенциалов запаздывающих по координатам эту величину ∆t= (R-∆R2)/c определить точно нельзя и можно найти только приблизительное ее значение по текущим координатам, то и Гербер вынужден написать, что ∆t=R/c, но он не пишет, что это приблизительное значение. Причем, здесь он еще и противоречит самому себе, т.к. у него скорость распространения потенциала равна непосредственно этой скорости c плюс скорость источника, т.е. будет c-Vr2, а он использует только c, но зато, наконец-то, находит расстояние R - ∆R в нужном ему виде

 

R - ∆R = R*[1 – (∆R/∆t) / c]                             (11)

 

И, хотя у него здесь неизвестно не только ∆R, но и ∆t, и если подставить в формулу найденное им значение ∆t=R/c, то мы опять придем к тому, что R - ∆R= R - ∆R, но тут получается, что при малых ∆R можно записать, что ∆R/∆t это будет dR/dt = Vr, а последнюю легко найти, зная Vr1 и Vr2. Вот так Гербер получил ошибочную формулу для учета запаздывания потенциала по координатам (12), но, если бы он учитывал только запаздывание по координатам потенциала источника (масса М) и не заигрывал бы с суммарной скоростью распространения потенциала, то он бы мог получить более правильную формулу (12’).

 

Pg1  = G*M*m / [R*(1 – Vr/c)]         (12)

 

Pg1’ = G*M*m / [R*(1 – Vr2/c)]       (12’)

 

А вот на втором этапе, т.е. с динамическим давлением потенциала, Гербер расправился двумя строчками, написав, что «Так как далее скорость, с которой движущиеся потенциалы минуют друг друга, имеет значение  c - ∆R/∆t потенциал из-за расхода времени на свою передачу m, становится пропорциональным c / (c - ∆R/∆t) таким образом, находим» и дальше записывает окончательную формулу, которая после замены ∆R/∆t = dR/dt = Vr будет выглядеть как (9), а логически должна выглядеть как (9’)

 

Pg  = G*M*m / [R*(1 – Vr/c)^2]                          (9)

 

Pg’ = G*M*m*(1 – Vr/c) / [R * (1 – Vr2/c)]       (9’)

 

Ведь, когда он пишет, что «потенциал из-за расхода времени на свою передачу», то здесь, как я понимаю, он в скрытом виде говорит именно о динамическом давлении потенциала, т.е. чем больше времени потенциал будет передаваться массе m, тем меньшее значение его передастся в единицу времени, т.е. тем меньше будет мощность поля и, следовательно, сила притяжения. По крайней мере, Роузвер понял эти его две строчки также как и я, т.к. он пишет

«Получение потенциала составляет существенную часть теории Гербера, поскольку его формула (6.24) выведена при рассмотрении влияния, которое конечное время распространения оказывает, во-первых, на расстояние и, во-вторых, на количество полученного потенциала. Во втором случае возможна интерпретация в рамках представлений Ритца о гипотетических частицах: чем быстрее поглощающее тело движется по направлению к излучающему, тем больше частиц «собирает» оно, и, аналогично, чем быстрее движется излучающее тело, тем быстрее движутся излучаемые им частицы и тем больше поглощается их в единицу времени. Однако детали этого процесса представить трудно».

 

 

Да, детали представить трудно, но логику нарушать нельзя в любом случае. А у Гербера в формуле (9) получается все наоборот. Ведь у него Vr положительна при удалении тел друг от друга и получается, что потенциал будет увеличиваться в этом случае. А логически надо было эту составляющую записать в числителе, чтобы при удалении тел передаваемый потенциал уменьшался, т.е. уменьшалось динамическое давление поля. И тогда формула (9’) будет соответствовать логике рассуждений Гербера для двух эффектов от запаздывания потенциала. Хотя, может быть, мы с Роузвером неправильно поняли Гербера и он рассуждал как раз наоборот, что чем дольше будет передаваться потенциал, тем больше его передастся. Ведь он писал, что передаваемый потенциал будет "пропорциональным c / (c - ∆R/∆t)", а тогда, действительно, логически правильно будет записать формулу (9). И в этом Гербера поддерживает  автор русского перевода его работы вот здесь [76].

 

 

Но получение так называемых запаздывающих потенциалов Вебера или Гербера это только пол дела, т.к. для описания поведения систем нам нужны ускорения масс системы, а их мы определяем через силы, действующие на тело с массой m. И тут мы видим во всей красе нелепости, вытекающие из формального энергетическо-геометрического подхода Лагранжа в механике. Ведь Гербер, не моргнув глазом, имея ввиду свое выражение (9), где у него в тексте Pg обозначено как V, пишет, что «Здесь в этом выражении для V содержатся не только R  но и производная от R по времени. Поэтому благодаря общим уравнениям движения по Лагранжу для ускорения массы m получаем» и применяет колдовские манипуляции Лагранжа для своего потенциала точно так же, как и для кинетической энергии (14). Правда перед этим Гербер переписывает свою формулу (9) заявляя «откуда с помощью биноминального разложения с сохранением членов до 2-го порядка малости следует» (13) и потом уже с этой формулой проделывает манипуляции (14).

 

Pg = G*M*m*(1 + 2*Vr/c + 3*Vr^2/c^2)  / R                   (13)

 

(dT/dR)/m – d(dT/dVr)/dt /m = dP/dR – d(dP/dVr)/dt        (14)

 

 

Здесь сразу бросается в глаза нелепость того, что Гербер ускорение в левой части приравнял силе в правой части, но тут он, можно сказать, не виноват, т.к. это у Лагранжа в Аналитической механике  получается путаница с массами, когда он выводит свою основную формулу динамики (подробности по этой ошибке Лагранжа смотрите в моей статье [67]). А сейчас давайте закроем на это глаза и выполним манипуляции предложенные Гербером. Что касается правой части, то здесь действительно получается выражение для силы (6-2), которое и получил Гербер. А вот, проделав указанные манипуляции с кинетической энергией T= m*Vr^2/2 в левой части, мы получим минус d2R/dt2, т.е. и ускорение в левой части и сила притяжения в правой части у нас получились со знаком минус или можно сказать, что обе величины будут с плюсом. Таким образом, уравнение, полученное с использованием выражения (14) будет отражать у Гербера не закон тяготения, а закон отталкивания. Странно только, что даже эту элементарную ошибку никто у Гербера не заметил, т.е. мы видим, что оболванивание ученых уравнениями Лагранжа уже в то время было на очень высоком уровне. К тому же я ума не приложу какое отношение к кинетической энергии тела m в формуле (14) имеет скорость тела M относительно него, т.к. Vr это у него относительная скорость двух тел.

 

 

Но и это еще не все. Выше я писал, что Гербер делал приближенное аналитическое решение, т.к. дифференцировал не свое точное выражение для потенциала (9), а приблизительное, т.е. после его разложение в ряд Тейлора с сохранением трех членов (13). А, если проделать его манипуляции (14) с выражением (9), то мы получим силу (6-2’). Кого интересуют подробности всех этих математических преобразований могут посмотреть их в файле Lagrange_Gerber.mws для математического пакета Maple, который я выложу вместе со статьей.

 

Fg=  - (G*M*m / R^2) * (1 – 3*Vr^2 / c^2 + 6*R*Ar /c^2)                                        (6-2)

 

Fg2= - G*M*m*c^2 * (c^2 – 4*c*Vr + 3*Vr^2 + 6*R*Ar) / [R^2 * (c -Vr)^4]         (6-2’)

 

 

Причем, здесь надо отметить, что у Гербера c=Vgr+Vi, т.е. скорость распространения гравитации складывается из непосредственно скорости гравитации и скорости источника. Но аналитически решить такую задачу Гербер бы не смог, поэтому в решении он принимает c=Vgr, что, впрочем, практически никак не может сказаться на величине смещения перигелия Меркурия, если рассматривать как Гербер его взаимодействие только с Солнцем, т.к. скорость последнего будет ничтожно мала по сравнению со скоростью Меркурия. Впрочем, и сама скорость Меркурия в его формулах (6-2) и (6-2’) практически никак не влияет на величину этого смещения и результат будет в обоих случаях 43 угл.сек (результаты я приведу позже) не смотря на то, что в одной формуле 3*Vr^2 с плюсом, а в другой с минусом, т.к. в этих формулах смещение обеспечивает ускорение, а не скорость, а оно в обеих формулах одинаковое. Хотя, конечно же, сама сила взаимодействия, рассчитанная по приближенной формуле (6-2) и по точной (6-2’) будет разной. И при решении других задач с этими формулами мы получим существенно другие результаты, но останавливаться на этом сейчас я не буду, т.к. в работе [67] я специально рассмотрел вопрос о том, что аналитическое решение реальных задач всегда является очень приблизительным. 

 

 

А сейчас мы видим, что уравнение для смещения перигелия Меркурия Гербер получил, проделав массу математических фокусов, которые никто, почему-то, за сто лет так и не разоблачил. Но, нас должны интересовать не математические фокусы, а принципы механики, которые позволяют получить уравнения описывающие поведение системы и подробно ознакомиться с этим Вы сможете в работе [67]. А тут мне стало интересно – а где это Гербер нашел у Лагранжа эти колдовские манипуляции с потенциальной энергией (15), которые позволили ему получить из потенциала (9) силу притяжения (6-2), и специально просмотрел еще раз Аналитическую механику Лагранжа [21], но ничего подобного не нашел. Да я даже не нашел широко известных манипуляций с выражением L=T-U (16), которое называется лагранжианом и где у нас T это кинетическая энергия, U потенциальная, R это обобщенная координата, а Vr это первая производная от этой координаты.

 

F= dP/dR – d(dP/dVr)/dt                 (15)

 

d(dL/dVr)/dt – dL/dR =0                 (16)        

 

d(dT/dVr)/dt – dT/dR + dV/dR =0   (17)

 

 

Тогда я стал смотреть более внимательно, чтобы найти хоть что-то похожее на манипуляции с лагранжианом, но нашел только его манипуляции с функциями T и V (17), где у него T это половина живой силы, а V это силовая функция, которая в принципе может быть и не потенциальной энергией, но и подробно об этом тоже смотрите в [67]. Тогда я вспомнил, что это у Гамильтона в его принципе наименьшего действия (ПНД) функционал равен T - U, а у Лагранжа в его ПНД функционал равен T и т.к. все эти манипуляции с функционалом вытекают из уравнения Эйлера для поиска экстремалей в вариационном исчислении, то по идее формулу (16) должен был придумать Гамильтон. Стал смотреть его работы, что у меня были, но тоже никакого упоминания о лагранжиане не нашел. Кругом он пишет только о своей главной функции, которая равна интегралу от T + U. Но вот в статье Полака [29] я нашел такую фразу «где L – функция Лагранжа (в частном случае разность кинетической и потенциальной энергии). Она введена для случая, когда H = H(t), Гамильтоном и обобщена в форме (33) Остроградским». Таким образом, мы хотя бы установили, что не только выражения (15) Гербер у Лагранжа найти не мог, но и все ссылки кого бы то ни было на Лагранжа, когда используется лагранжиан, неправомерны.

 

 

Здесь приходиться вспомнить критику Герцем энергетического подхода, где потенциальная энергия не только не локализована, но даже не говориться от чего она может зависеть. И в свое время многие ученые были категорически против того, что потенциальная энергия у Вебера стала зависеть от скорости. Например, это Гельмгольц и Зеелингер, а последний даже задал головоломку математико-физикам заявив, что если к потенциалам Вебера применить манипуляции (15), которые позволили найти силу взаимодействия Герберу, то у Вебера получится сила взаимодействия (6-1), где у второго и третьего членов будут другие знаки. А в случае с Меркурием это означает, что у него перигелий должен смещаться в другую сторону. Но тут математко-физики, как всегда придумали очередной хитрый трюк и Роузвер пишет, что к потенциалам Вебера эти манипуляции (15) применять нельзя. Их можно применить к лагранжиану Вебера, который будет выглядеть так

 

Lw= - G*M*m * [1 + Vr^2 / (2*c^2) / R                (18)

 

 

Выражение получилось странное, а откуда оно взялось это вообще тайна, но давайте попробуем применить манипуляции, которые рекомендует Роузвер. В результате мы действительно получим, что знаки у 2-го и 3-го членов станут такими как надо, но сила то получится уже со знаком плюс, т.е. будет не притяжение, а отталкивание. Тогда давайте попробуем применить к этому выражению манипуляции как к лагранжиану, т.е. для получения дифференциального уравнения. Действительно получаем выражение, совпадающее с силой Вебера, но мы должны получить дифференциальное уравнение, а не одну только силу. Но для этого в лагранжиане должна быть и кинетическая энергия, а ее в уравнении (18) нет и поэтому не понятно, что это выражение вообще из себя представляет. Более того, манипуляции (14) Гербера, которые он взял у Лагранжа, соответствуют манипуляциям Гамильтона (20), если все знаки заменить на противоположные, а значит и у ускорения тоже получится отрицательный знак, т.е. мы опять таки получим закон отталкивания, а не притяжения.

 

(dT/dR)/md(dT/dVr)/dt /m = dP/dRd(dP/dVr)/dt        (14)

 

d(dT/dVr)/dt - d(dU/dVr)/dt – dT/dR + dU/dR =0     (20)

 

 

В общем, мы видим полное не понимание математико-физиками механики и все их способности в механике заключаются только в манипуляциях Лагранжа-Гамильтона с произвольными формулами для энергий, что напоминает элементарное передергивание карт. А кто в этом виноват. Конечно же, Лагранж со своей Аналитической механикой, которая, по моему мнению, нанесла вред науке даже больше, чем потом принцип наименьшего действия Гамильтона и курсы физики Ландау и Фейнмана. Хотя, в оправдание Фейнмана могу сказать, что хоть и редко, но он задумывался о том, что пишет. А основное уравнение динамики Лагранжа все забыли очень быстро, т.к. там и забывать то было нечего, и поэтому, говоря о его механике, говорят только об особом математическом приеме, позволяющем простым способом сократить число переменных, т.е. говорят о его формуле (17) вытекающей из его ПНД, но, почему-то, применяют формулу (16) вытекающую из ПНД Гамильтона. О различиях между этими двумя принципами и о других вариационных принципах я подробно написал в своей статье [8] и останавливаться на этом не буду. Вообще то, число переменных можно сократить и другими способами, но Лагранж пишет, что так проще. А в настоящее время этот способ, предложенный Лагранжем, не имеет никакого практического значения, т.к. при наличии компьютеров все реальные задачи механики решаются численными методами и здесь нам практически безразлично, сколько у нас переменных.

 

 

Я уже писал выше, что, например, я решаю сейчас таким методом 33 дифференциальных уравнения второго порядка, описывающих поведение планет в Солнечной системе, а в то время, чтобы хоть как-то решить эту задачу, т.к. задача трех тел уже не решается аналитически, вынуждены были изыскивать различные приемы хотя бы для приблизительного решения этой задачи. То есть речь здесь совершенно не идет о механике и это чисто математическая задача. И надо сказать, Лагранж, как математик, с ней справился, введя свою возмущающую функцию, а Гамильтон использовал для этого свою характеристическую функцию. Хотя, конечно же, без теории Луны Эйлера (1753 год) они бы ничего не сделали. Мало того, что теория Луны на порядок сложнее теорий планет, поэтому Лагранж и Гамильтон, как математики, справились только с задачей по теории планет, так вдобавок Эйлер и показал им как надо решать эту задачу, а ведь он писал эту работу уже будучи слепым на оба глаза. При этом замечу, что теорий планет потом было много, а следующая теория Луны Хилла появилась только через 104 года после теории Эйлера. Много для решения задач в астрономии сделали и Пуассон и Лаплас, и именно поэтому астрономы смогли найти аномальные смещения параметров орбит, которые я рассматривал в предыдущей статье, но я их ищу уже другими методами.

 

 

Вообще-то, тут надо сказать и то, что весь 2-ой том Аналитической механики Лагранжа, посвященный математическому решению конкретных задач, является очень полезным, а вот от его 1-го тома науке только вред. Да и собственно механики в его Аналитической механике ничуть не больше чем в квантовой механике или релятивистской механике, а ту, что была, он добил, связавшись с формализмом ПНД. А вот Эйлер, как только узнал, что на другой стороне шара ПНД не работает, так сразу о нем забыл, хотя именно Эйлер и дал математическую формулировку этому принципу, т.е., по сути, он его и создал. Да и сам Лагранж ведь писал о ПНД, что «Таков тот принцип, которому, хотя и не вполне точно, я даю здесь название принципа наименьшего действия и на который я смотрю не как на метафизический принцип, а как на простой и общий вывод из законов механики». Но нет, наверное, это черт его дернул изложить более простой способ уменьшения переменных и теперь его механику и всю современную физику связывают только с ПНД, а он считал главным в ней свою основную формулу динамики.

 

 

Вот только я не понимаю (вернее догадываюсь) почему представители официальной науки так расхваливают его Аналитическую механику, т.к. там вообще нет механики, а его основная формула, которая, вроде бы должна быть механикой, является просто недоразумением. Тем более, что эти глупости Лагранж писал в 1788 году, когда был знаком уже не только с двумя томами Аналитической механики Эйлера (1736 год), которая называлась «Механика, или наука о движении, в аналитическом изложении» и с третьим томом 1765 года, который назывался «Теория движения твердых тел», а также с теорией Луны 1772 года. Более того, в моем распоряжении находится 2-ое издание Аналитической механики Лагранжа 1813 года, т.е. года, когда Лагранж умер и, следовательно, он мог до переиздания внести поправки в свой первый том, но он сделал только добавления к 1-у изданию, которые ничего не изменили. А, чтобы Вам было более понятно, о чем я говорил, критикуя Лагранжа, я отправляю Вас к своей статье [67], где я на трех страничках изложил всю его Аналитическую механику, а затем указал на его элементарные ошибки, которые говорят о его полном не понимании механики.

 

 

Ну, а теперь давайте вернемся непосредственно к различным теориям, учитывающим конечность скорости гравитации и поговорим о потенциалах Лиенара-Вихерта. Да, они так же, как и потенциалы Вебера, создавались для взаимодействия зарядов, но их точно так же можно применить и к массам, если не учитывать эффектов, возникающих от магнитных полей. Впервые выражения для этих потенциалов получили Лиенар в 1898 году [70] и Вихерт в 1901 году [71], т.е. получены они были до формальной даты создания СТО, т.е. до 1905 года, но сейчас считается, что они полностью соответствуют СТО, хотя это не так. Да, влияние идей СТО просматривается в них, т.к. работы Лоренца в это время активно обсуждались и, конечно же, и Лиенар и Вихерт не хотели быть в стороне от мэйн стрима того времени, поэтому и отразили эти идеи в этих потенциалах. Но при этом не менее важным для их создателей было и стремление получить выражения для потенциалов, которые зависят от скорости, т.к. уже стало ясно, что потенциалы Вебера не отвечают модным идеям СТО. При этом потенциалы Лиенара-Вихерта часто называют еще запаздывающими потенциалами, но как такового запаздывания, т.е. запаздывания потенциалов по координатам в них нет, а есть только учет запаздывающих скоростей и ускорений. Таким образом, получается, что сам потенциал определяется по законам Ньютона, т.е. при скорости гравитации равной бесконечности, а динамическое давление создаваемое этим потенциалом определяется исходя из конечности скорости гравитации.

 

 

Мне это сочетание кажется странным, но не менее странным и мутным является и сам способ получения формул для этих потенциалов. Сейчас наиболее известны два вывода - Ландау [4] и Фейнмана [5]. При этом Ландау дает вывод для точечных зарядов, а Фейнман строит весь свой вывод именно на том, что рассматривает движущийся заряд размазанным по объему, но в конце вывода пишет, что так как в полученной им формуле отсутствует размер заряда, то эти потенциалы справедливы и для точечного заряда. При этом Фейнман вычисляет странный интеграл, который у него пролетает над зарядом и в результате он суммирует потенциалы от разных частей движущегося объемного заряда так как будто они долетают до пробного заряда в один и тот же момент времени, что является явным нонсенсом. А Ландау тот вообще не очень то вникает в тонкости вывода, хотя начинает все вроде бы нормально. Сначала он записывает уравнение (63,1*), которое соответствует моей формуле (2), а потом уравнение (63,2*), которое соответствует моей формуле (3), и которые вытекают из элементарного анализа запаздывания потенциалов согласно рис.2.

 

 

t’ + R(t’) / c = t               (63,1*)

 

φ = e / R(t’)                     (63,2*)                                                     

 

t'= t - ∆t = t - R' / Vgr                                                              (2)

 

φ = G*m2 / R’                                                                         (3)

 

 

Рис.2 Запаздывание потенциала в точке Р от движущегося заряда 2.

 

 

 

 

Но вот тут Ландау делает хитрое заявление о том, что оказывается в его формуле (63,2*) приведено выражение для потенциала заряда в момент времени t при условии, что он покоился в момент времени t' и, чтобы получить выражение потенциала для движущегося заряда, ему надо срочно уйти в 4-х векторы, где он не особенно это расписывая делает некоторые преобразования, но от туда он возвращаемся с уже готовыми формулами для запаздывающих скалярного φ и векторного A потенциалов (63,5*), где все величины теперь должны быть взяты в запаздывающий момент времени. Выполним эту замену и получим формулы (63,5')

 

φ = e / (RV*R / c

                       

A  = e * V / c  / (R – V*R / c)           (63,5*)

               

φ = e / Rf

 

A = φ * V' / c                                     (63,5')                     

 

где я сделал замену Rf= R' – V' * R' / c . Здесь мы видим, что не только скалярный потенциал, т.е. по сути потенциальная энергия, стала зависеть от скорости, что и требовалось от этих потенциалов, но еще и появился векторный потенциал A, т.е. по сути теперь потенциальная энергия еще и перестала быть скаляром (тут надо отметить и то, что у Фейнмана векторный потенциал определяется как φ * V' / c^2).  При этом у нас получились выражения, где нет никакого запаздывания самих потенциалов, т.к. RfR а это возможно только при скорости электрического взаимодействия равной бесконечности. Ну, а в связи с тем, что у нас теперь появился еще и векторный потенциал, которого не было ни у Вебера, ни у Гербера, то, чтобы найти теперь напряженность поля, создаваемую этими потенциалами, уже недостаточно или просто взять производную от потенциала по Rf или проделать манипуляции Лагранжа с потенциальной энергией, чтобы найти силу взаимодействия. Поэтому для этого придумали другие хитрые манипуляции для потенциальной энергии не только зависящей от скорости, но и являющейся вектором, чтобы рассчитать напряженность поля.

 

E = - (dA/dt) / c - grad φ

 

 

Ну, а результатом таких манипуляций является формула (63,8*) у Ландау, но я ее запишу в том виде, как она дана у [72] - это его формула (5.17*), которая преобразуется потом в формулу (5.20*) идентичную формуле (63,8*) у Ландау и, где a' это вектор ускорения в запаздывающий момент времени.

 

E = (e / Rf^3 / c^2) * [(R' - V' * R' / c) * (R' * a' + c^2 - V'^2) - R' * Rf * a']                 (5.17')

 

 

Как нетрудно заметить, выражение R' - V' * R' / c это тот же самый фиктивный радиус RfR и поэтому при малых скоростях и ускорениях мы получим просто напряженность поля покоящегося заряда, но при больших скоростях и ускорениях у нас уже получится существенно другой результат. Вот только здесь получается одна неувязочка с требованиями СТО. Ландау утверждает, что эта формула получена для произвольной системы отсчета, но у нас не может быть ИСО связанных с обоими зарядами, т.к. они могут двигаться и с ускорением, а в этом случае это будут уже не ИСО. К тому же, скорость пробного заряда может быть абсолютно любой и это никак не повлияет на результат, т.к. эта скорость не входит в формулу (5.17'), поэтому утверждение Ландау об относительных скоростях в этой формуле видится надуманным. А вот то, что эта формула не учитывает именно запаздывание потенциалов по координатам говорит уже сама СТО, т.к. напряженность поля движущегося заряда, полученная с использованием преобразований Лоренца, дает тот же эллипсоид Хэвисайда, что и формула (5.17') и Ландау об этом сам пишет в конце 63 параграфа. А этот эллипсоид явно не учитывает запаздывание потенциалов по координатам.

 

 

 

И на то, что преобразования Лоренца не учитывают запаздывание потенциалов, т.е. скорость гравитации в них равна бесконечности, неоднократно указывали, например, авторы работы [13]. Да и в работе [14] это тоже подтверждается, когда речь идет об эллипсоиде Хевисайда, который, как утверждают сторонники СТО как раз и демонстрирует запаздывание потенциалов. Здесь необходимо заметить, что этот эллипсоид (см. рис 6.b) является просто сплюснытым по оси X и растянутым по оси Y полем статического заряда, т.е. никакое запаздывание потенциалов он не отражает. Но в этом случае скорость движения источника V создающего потенциал должна быть равна нулю (см. рис 6.а) и, следовательно, вообще никаких преобразований СТО не надо делать. А, если принять, что на рис. 4.a мы имеем движущийся заряд, то такая картина линий равной напряженности поля может наблюдаться только при условии, что это статическое поле является упругой субстанцией связанной с зарядом и простирающейся до бесконечности или для динамического поля при скорости его распространения  равной бесконечности, и при бесконечно большой частоте кадров, но изобразить это на рисунке в процессе движения заряда не возможно, поэтому сторонники СТО просто рисуют картинку статического поля.

 

Рис. 6 Поверхности равных напряженностей поля (эквипотенциали) при разных текущих расстояниях от движущегося заряда (или диаграммы напряженностей поля при одинаковых текущих расстояниях от движущегося заряда) при скорости распространения потенциала: a – бесконечной, c – конечной и искажение этих линий напряженностей от преобразований Лоренца, соответственно, b и d.

 

 

Здесь надо также заметить, что Хэвисайд получил свой эллипсоид в 1882 году, т.е. не только до СТО и формул Лиенара и Вихерта, но даже до преобразований Лоренца, хотя размеры для эквипотенциалей у него изменяются согласно этим преобразованиям точно так же, т.е. по оси X сжимаются в sqrt(1 – V^2 / c^2) раз, а по оси Y во столько же раз увеличиваются, хотя я считаю, что из логики СТО должно следовать увеличение напряженности поля именно вдоль движения заряда, т.к. линейные размеры при этом должны уменьшаться. Интересно так же заметить, что и выражение для силы Лоренца, Хэвисайд получил в 1889 году, т.е. за три года до Лоренца. А вот с учетом именно запаздывания потенциала по координатам, т.е. при учете конечности скорости распространения взаимодействия, мы должны получить для движущегося в неподвижной системе отсчета заряда картинку 6.с. И, если теперь, вследствие движения заряда, сжать и растянуть эту картинку согласно требованиям СТО для движущейся системы, то мы должны получить картинку 6.d, которая никак не эквивалентна картинке 6.b. Но сторонники СТО упорно доказывают, что эллипсоид Хэвисайда как раз и отражает запаздывание потенциалов, но не уточняют, что он отражает запаздывание потенциалов не по координатам, а по скоростям и ускорениям, т.е. учитывает только динамическое воздействие поля и то в их симметричной интерпретации.

 

 

И еще. По поводу эллипсов. Я тут написал программу Potencial2 и построил на ней диаграммы напряженности поля, вернее сил притяжения пробного заряда к движущемуся заряду, с использованием формулы (5.17') и у меня получаются эллипсы при равномерном движении заряда только при его относительно небольших скоростях. Например, на рис.7 слева даны силы притяжения пробных зарядов к движущемуся, который находиться в центре, а справа эти силы спроецированы в центр этого движущегося заряда, т.е. получился так называемый ёжик напряженности при скорости заряда 0,5*c, а вот на рисунке 8 дано тоже самое, но при скорости заряда 0,9*c. Как видим, вместо эллипса мы получили восьмерку. На всякий случай даю и кусок кода программы, где приведен аналогичный расчет и для потенциалов Вебера, Гербера и моих, т.е. запаздывающих по координатам, а также дан расчет преобразований Лоренца для поля движущегося заряда.

     

 

 

Рис. 7 Диаграмма напряженностей поля (сил) при одинаковых текущих расстояниях от заряда движущегося со скоростью 0,5*с.

 

 

Рис. 8 Диаграмма напряженностей поля (сил) при одинаковых текущих расстояниях от заряда движущегося со скоростью 0,9*с.

 

 

For i = 0 To 360 / N

An(1) = i * N  '                                       угол луча зрения в текущий момент от оси Х

X(1) = Xn2(1) + R0(1) * Cos(An(1) / kUgol) 'находим текущую координату Х(1)

Y(1) = Yn2(1) + R0(1) * Sin(An(1) / kUgol) 'находим текущую координату y(1)

Picture1.Circle (X(1) / ML, Y(1) / ML), 0.2, RGB(0, 0, 250) 'рисуем координаты приемника

 

If Option1.Value = True Or Option2.Value = True Then '         если Лиенар-Вихерт или Юдин

dX = X(1) - Xn2(1): dY = Y(1) - Yn2(1)

a = VX(2) ^ 2 + VY(2) ^ 2 - Vs ^ 2

b = -2 * dX * VX(2) - 2 * dY * VY(2)

c = dX ^ 2 + dY ^ 2

T1 = -(-b + Sqr(b ^ 2 - 4 * a * c)) / 2 / a '   первый корень

T2 = -(-b - Sqr(b ^ 2 - 4 * a * c)) / 2 / a '   второй корень

T00 = T1: If T1 < 0 Then T00 = T2

X(2) = Xn2(1) - VX(2) * T00 '                     находим запаздывающую координату Х(2)

Y(2) = Yn2(1) - VY(2) * T00 '                     находим запаздывающую координату y(2)

dX = X(1) - X(2): dY = Y(1) - Y(2): R = Sqr(dX ^ 2 + dY ^ 2):

An(1) = ArcTangens(dX, dY) '                      находим новый угол An(1)

End If

 

Fp(1) = Kq * Q2 * Q1 / R ^ 2 '                                                       Кулон, Юдин-окончательно

Qn1(1) = An(1) - QV1: Qn2(1) = An(1) - QV2 '                                         угол между радиус-вектором и скоростью

dR = V(1) * Cos(Qn1(1) / kUgol) - V(2) * Cos(Qn2(1) / kUgol) '                     это для Вебер, Гербер

 

If Option2.Value = True And Check2.Value = 1 Then '                  Юдин если с давлением

Fp(1) = Fp(1) * (1 - V(1) * Cos(Qn1(1) / kUgol) / Vs) / (1 - V(2) * Cos(Qn2(1) / kUgol) / Vs)

End If

 

If Option12.Value = True Then '                                                      Вебер

Fp(1) = Fp(1) * (1 - 0.5 * dR ^ 2 / Vs ^ 2 + 1 * R * A2 * Cos(Qn2(1) / kUgol) / Vs ^ 2)

End If

 

If Option13.Value = True Then '                                                      Гербер

Fp(1) = Fp(1) * (1 - 3 * dR ^ 2 / Vs ^ 2 + 6 * R * A2 * Cos(Qn2(1) / kUgol) / Vs ^ 2)

End If

 

If Option14.Value = True Then '                       преобразования Лоренца (аналог Лиенара-Вихерта)

Fp(1) = Fp(1) * (1 - V(2) ^ 2 / Vs ^ 2) / (1 - V(2) ^ 2 * (Sin(Qn2(1) / kUgol)) ^ 2 / Vs ^ 2) ^ 1.5

End If

 

Fx = Fp(1) * Cos(An(1) / kUgol): Fy = Fp(1) * Sin(An(1) / kUgol)

 

If Option1.Value = True Then '                                     Лиенар-Вихерт

Rfik = R - R * V(2) * Cos(Qn2(1) / kUgol) / Vs ' фиктивный радиус

If Check3.Value = 0 Then '                                    формула Болюсина 5.20 (до упрощения)

Fx = (Kq * Q1 * Q2 / Rfik ^ 3 / Vs ^ 2) * ((Vs ^ 2 - V(2) ^ 2) * (dX - VX(2) * R / Vs) + _

((dX - VX(2) * R / Vs) * R * A2 * Cos(Qn2(1) / kUgol) - AX2n * R * (R - R * V(2) * Cos(Qn2(1) / kUgol) / Vs)))

Fy = (Kq * Q1 * Q2 / Rfik ^ 3 / Vs ^ 2) * ((Vs ^ 2 - V(2) ^ 2) * (dY - VY(2) * R / Vs) + _

((dY - VY(2) * R / Vs) * R * A2 * Cos(Qn2(1) / kUgol) - AY2n * R * (R - R * V(2) * Cos(Qn2(1) / kUgol) / Vs)))

End If

If Check3.Value = 1 Then '                                     формула Болюсина 5.17

Fx = (Kq * Q1 * Q2 / Rfik ^ 3 / Vs ^ 2) * ((dX - R * VX(2) / Vs) * (R * A2 * Cos(Qn2(1) / kUgol) + Vs ^ 2 - V(2) ^ 2) - R * Rfik * AX2n)

Fy = (Kq * Q1 * Q2 / Rfik ^ 3 / Vs ^ 2) * ((dY - R * VY(2) / Vs) * (R * A2 * Cos(Qn2(1) / kUgol) + Vs ^ 2 - V(2) ^ 2) - R * Rfik * AY2n)

End If

Fp(1) = Sqr(Fx ^ 2 + Fy ^ 2)

End If

 

 

Таким образом, ни потенциалы Вебера, ни потенциалы Лиенара-Вихерта не учитывают запаздывание потенциалов по координатам, а вот, если бы создатели законов электродинамики и теории поля при их создании учитывали это, то электродинамика могла бы получиться совсем другой и в ней не было бы столько различных парадоксов. Например, современная электродинамика совершенно не может объяснить работу EH антенн [15], но эти антенны не только работают, но и хорошо работают в воде и земле. А вот, если в законах электродинамики учесть запаздывание по координатам, то возможно, что нам не только не потребуется сила Лоренца, но и вообще магнитное поле. Рассмотрим рис.9, где у нас имеются два движущихся заряда (1 заряжен отрицательно и 2 заряжен положительно), которые движутся со скоростями V1 и V2. Давайте найдем силу, с которой заряд 2 будет воздействовать на заряд 1. Сейчас эту силу вычисляют так. Находят силу электростатического притяжения Fq, которая действует без всякого запаздывания из точки 2, где в данный момент находится заряд 2. Затем находят силу Лоренца Fm, которая в нашем случае будет направлена вверх, т.к. магнитное поле, создаваемое зарядом 2, по правилу правого винта будет направлено к нам, и по правилу левой руки для положительного заряда сила Fm будет направлена вниз, а для отрицательного в другую сторону. Теперь находим суммарную силу Fsum1 и видим, что она очень похожа на силу электростатического притяжения с учетом запаздывания по координатам, т.е. когда заряд 1 будет притягиваться следом заряда 2, т.е. из точки 2’.    

Рис. 9 Схема к расчету силы притяжения движущимся зарядом 2 неподвижного заряда 1 с использованием силы Лоренца и запаздывания потенциала по координатам.

 

А, если мы сначала рассчитаем силу электростатического притяжения с учетом запаздывания, т.е. из точки 2’, а потом просуммируем ее с силой Лоренца, то получим суммарную силу Fsum2 и не понятно какая же из этих сил Fsum1 или Fsum2 является именно суммарной силой. Да, при разных направлениях скоростей и разных знаках зарядов получаются самые разнообразные комбинации суммирования сил и я  не утверждаю, что уже отменил силу Лоренца. Но ведь эта сила и другие законы электродинамики были получены при обработке экспериментальных данных без учета запаздывания по координатам, поэтому надо заново проанализировать экспериментальные данные с учетом эффекта запаздывания потенциала по координатам. Возможно, некоторые законы изменятся или появятся новые, а в результате и теория поля будет отличаться от теории поля Максвелла, у которого скорость распространения электромагнитного взаимодействия была равна бесконечности [69]. Вообще то, сегодняшнее изложение теории поля Максвелла имеет мало общего с там, что писал сам Максвелл, но именно запаздывание потенциалов по координатам и там не учтено. И при этом возможно, что в учебниках перестанут писать, что сила Лоренца вычисляется по известной формуле только при скоростях намного меньше скорости света, но не пишут почему. Я, конечно, могу предположить, что частично это связано с синхротронным излучением, но хотелось бы иметь четкие объяснения.

 

 

 

А то, что вопрос переработки электродинамики для математико-физиков действительно сложный, доказывает то, что пока даже те, кто старается устранить ошибки СТО и современной электродинамики, т.е. альтернативные ученые, путаются в вопросе о запаздывании потенциала по координатам. Так автор [14] правильно указывает на многие ошибки при выводе потенциалов Лиенара-Вихерта у Фейнмана, которые могут привести к формуле без запаздывания потенциалов, но на стр.51 у него формула (4.4*) записана без штрихов, а первоначально на стр.21 была записана со штрихами. Но, не смотря на это, он, учитывая запаздывание потенциалов,  все же приходит к выводу о том, что специфического магнитного поля не существует, а есть лишь специфическое проявление обычных радиальных электрических полей. И даже приводит цитату Ампера, который сказал «магнетизма нет, есть одно электричество». Но вот с формулой для запаздывания потенциалов он все же путается. Очень интересен и подход автора работы [34], который, применив преобразования Лоренца в своей трактовке, тоже приходит к выводу о том, что существующая сейчас формула для расчета силы Лоренца справедлива только при скорости электромагнитных взаимодействий равной бесконечности.

 

 

Но, не так просто учесть даже запаздывание скоростей и ускорений в аналитических расчетах и автор [40], ссылаясь на Парселла, говорит, что при взаимодействие двух зарядов движущихся не параллельно (или двух токов), надо бы учитывать потенциалы Лиенара-Вихерта. Однако, опасаясь больших сложностей такого решения, рассматривает этот вопрос взаимодействия двух токов с учетом продольных сил Николаева (только я думаю, что это еще одни силы Лоренца, но вид сбоку). Но я не хочу углубляться во все эти математические тонкости по которым желающие могут поспорить с математико-физиками Ландау и Фейнманом или авторами работ [13, 14]. Тем более, что со мной спорить по этим математическим тонкостям бесполезно, т.к. я в своих исследованиях придерживаюсь исключительно инженерной точки зрения решения задач, т.е. решаю их практически. Поэтому, у меня потенциалы  рассчитываются только в абсолютной системе отсчета, т.е. никаких преобразований Лоренца и прочих заумных действий у меня нет. При этом запаздывание потенциалов я тупо вычисляю методом приближений за несколько итераций, а решение системы дифференциальных уравнений, описывающих движение планет Солнечной системы, я осуществляю самым надежным методом численного решения, т.е. методом Рунге-Кутта по четырем коэффициентам. А кого интересуют вопросы точности решения дифференциальных уравнений, описывающих поведение реальных систем, численными методами, которые дают почти точное решение, и аналитическими, которые дают приблизительное или ошибочное решение, я отсылаю к своей статье [67].

 

 

А сейчас давайте рассмотрим один из современных альтернативных подходов учета динамического давления поля. Конкретно рассмотрим подход, который применяет очень известный и уважаемый альтернативщик Хайдаров, который пропагандирует теорию эфира. И, исходя из нее, он предложил идею уменьшения или увеличения динамического давления гравитационного поля на движущуюся в нем массу [3]. При этом он получил формулу (22) для давления эфира на массу m, где у него с- скорость света, а Vr =dR/dt – скорость относительного движения тел (положительная при удалении тел и отрицательная при сближении).

 

Hooke= m*c^2* (1 – 2*Vr/c +Vr^2/c^2)                (22)   

 

 

Рис. 10. Относительное движение внутренней планеты (Pinn) и внешней планеты (Pout) относительно средней планеты P, для которой ведется расчет силы тяготения, действующей на нее со стороны внутренней и внешней планет (воспроизведено из работы [3]).

 

 

А далее он, используя формулу Гербера, рассчитал аномальные смещения перигелиев планет, вызванные притяжением Солнца, а потом, используя формулу (22), посчитал смещения перигелиев от воздействия на исследуемые планеты других планет и потом уже их суммарные смещения (см. таблицу 2). При этом у него получается, что согласно формуле (22) на рис. 10 в нижних квадратах планета Pout удаляется от исследуемой планеты P и ее воздействие на нее меньше, чем дает формула Ньютона, а в верхних квадратах она приближается и ее воздействие больше, чем по формуле Ньютона. И наоборот планета Pinn в нижних квадратах приближается к планете P, а в верхних удаляется. В связи с этим, как пишет Хайдаров, тангенциальная составляющая круговых интегралов гравитационного воздействия будет не нулевой и воздействие от внешней планеты приведет к положительному движению перигелия планеты P, а воздействие внутренней планеты к отрицательному. И именно такое воздействие на исследуемую планету мы и видим в его таблице 2.

 

 

К сожалению, я не могу проверить данные, полученные Хайдаровым, т.к. он отказался сообщить мне формулу для расчета силы притяжения между двумя массами вытекающую из его эфирного подхода, т.е. давления эфира по формуле (22). А пытаться получить ее из его формулы давления эфира и с использованием тех расплывчатых рекомендаций, что он мне дал, не имеет смысла, т.к. он всегда сможет заявить, что я что-то не понял и поэтому получил неправильную формулу. А, как следует из его ответа на мою просьбу, правильную формулу он мне не сообщит, т.к. это коммерческая тайна. И, чтобы у читателей не было никаких вопросов, я свое последнее письмо ему и его ответ привожу полностью (посмотреть можно здесь http://bourabai.kz/guest/index.php – сообщение № 267).

 

 

 

Сергей Юдин (04.04.2013 23:08:41)

Добрый день Карим Аменович. Я извиняюсь за настойчивость, но Вы так и не ответили на мой вопрос о том какое выражение для силы притяжения между планетами Вы использовали при расчете смещений перигелиев планет, данные по которым Вы привели в таблице 1 своей статьи «Эфир великий часовщик». Вопрос был мною задан почти два месяца назад, но ответа я так и не получил. Время пока немного терпит, но процесс написания моей статьи, где я рассматриваю теории различных авторов, объясняющих смещения перигелиев планет, близится к концу, а Ваши данные я так и не проверил, хотя уже разместил и рисунок и таблицу из Вашей статьи в своей статье.
А проблема в том, что данные всех авторов я не только привожу в своей статье, но и провожу свой расчет, полученных ими значений по приведенным ими формулам. Вот только Ваши данные я проверить не могу, т.к. в своей статье Вы пишите об уменьшающейся или увеличивающейся силе притяжения между планетами, но приводите формулу только давления эфира, которое без Ваших пояснений ни о чем не говорит и трансформировать его в силу притяжения я не могу. А Вы в Ваших ответах пишите «Дорогой сергей, это зависит от конкретного случая, который Вас интересует, так как в отличие от потенциала, сила "привязана" к конкретной механической конфигурации». Но ведь данные, приведенные Вами в таблице 1, получены Вами по какой то конкретной формуле для конкретной ситуации. Вот я и писал Вам «Карим Аменович, меня интересует именно тот случай, где сила привязана к механической конфигурации, которую Вы использовали при расчете смещений перигелиев планет, данные по которым приведены в Вашей статье в таблице 1». И сейчас у меня та же убедительная просьба сообщить мне эту секретную формулу, чтобы я мог проверить сам Ваши расчеты.

В общем, то и без Вашей теории существует множество других теорий, объясняющих аномальные смещения перигелиев планет, но Ваша теория вроде бы самая передовая, т.к. используется давление эфира, и в тоже время Ваша теория очень разрекламирована в Интернете. По этому, я думаю, что мои читатели меня не поймут, если я ничего не скажу о Вашей теории. Но, пока я сам не проверю Ваши данные, я ничего определенного о Вашей теории сказать не могу, а, следовательно, я вынужден буду объяснить моим читателям почему я не могу этого сделать, т.е. я вынужден буду написать, что Вы отказались от проверки Ваших данных и поэтому они, как минимум, не могут считаться достоверными, а все Ваши исследования являются антинаучными, т.к. не поддаются проверке. Вообще то, я и так вижу, что Ваши данные ошибочные, но лучше бы было, чтобы Вы указали Вашу ошибочную формулу и я после ее проверки просто написал, что, к сожалению, и Ваша теория не позволяет адекватно описать процесс смещения перигелиев планет. Но, когда Вы отказываетесь от такой проверки, то это просто смертельно для Вашей репутации. Надеюсь, что Ваша репутация то для Вас не пустой звук, т.к. в одном из сообщений Вы написали «Если серьезно, наш форум, на мой взгляд, посвящен главной проблеме науки: проблеме совести».

С наилучшими пожеланиями Сергей Юдин.


 

Дорогой Сергей,
из Вашего длинного послания я вижу, что Вы совершенно отчетливо понимаете, о чем идет речь, но пытаетесь педалировать на совесть.
Это педалирование действует только на тех, кто не умеет мыслить.
Меня это не проймет, - как это видно из того, что я для читателей сохранил, а не удалил этот Ваш пост, где Вы уже "поливаете" меня, - лишь только бы достичь своей цели - получить бесплатно технологию моих вычислений.
Стоимость этой новой технологии существенно превышает стоимость существующих технологий расчета полетов, которые, думаю, известно или по крайней мере понятно читателям, - стоят очень и очень дорого и по существу во многом определяют первенство той или иной страны в космосе...
Не зря пословица гласит:
"Что бы ни говорил человек - он говорит о себе"
У меня нет ложной гордыни: что там будет судить о мне г-н Юдин или еще кто-то. Можете и дальше продолжать в том же духе, но этот Ваш пост дезавуирует Вас самого.
Судить будет История, ибо, как известно с античных времен, - Истина - дочь времени, а сейчас судит здравый смысл. У кого его нет, тот проигрывает.
Здравый смысл советует мне придерживаться всем известного "принципа нераспространения" технологий, а того объема, в котором мой эфирный подход уже изложен публично, достаточно для понимания его чисто научных принципов.
Искренне, КХ

 

 

Ну, что же попробуем хотя бы проанализировать данные, приведенные самим Хайдаровым, если мы их не можем проверить. Во-первых, сразу видим методологическую ошибку, т.к. он для одного и того же природного явления – гравитационного взаимодействия между телами использует одновременно две разные формулы. Основное смещение от воздействия планет на исследуемую планету он вычисляет по формуле, которая вытекает из его формулы (22), а от воздействия Солнца на эту планету появляется добавка, которую он получает по формуле Гербера. При этом, что интересно, суммарные значения смещений у меня в таблице 3, взятой из работы [1], и у Хайдарова в таблице 2, как ни странно, получились очень похожими, но вот составляющие от воздействия других планет совершенно другие. Не менее интересно и то, что составляющие перигелия Меркурия от воздействия других планет у Хайдарова практически совпали с полученными мною значениями, которые в свою очередь у меня практически совпадают с данными по всем планетам полученными другими астрономами. И, если отталкиваться от данных по Меркурию, то надо сделать вывод о том, что при расчете воздействия на Меркурий других планет у Хайдарова, с использованием его формулы, полностью отсутствует динамическое давление, т.к. и я и другие астрономы рассчитывали это воздействие с использованием только законов Ньютона, т.е. без учета динамического давления.

 

 

Таблица 2. Рассчитанные по теории Гербера-Хайдарова со скоростью гравитации равной скорости света значения смещений перигелиев орбит планет (в угл.сек за 100 лет) от действия на них отдельных планет по Хайдарову и Солнца по Герберу (воспроизведено из работы [3]).

 

 

 

 

Таблица 3. Рассчитанные численными методами за период времени с 1801 по 2000 год по теории Ньютона, т.е. с бесконечной скоростью гравитации, значения смещений перигелиев орбит планет (в угл.сек за 100 лет) от действия на них отдельных планет с системой масс JPL.

 

 

Меркурий

Венера

Земля

Марс

Меркурий

+0,006+/-0,003

-147,945+/-0,204

-13,585+/-0,099

+0,763+/-0,028

Венера

+275,950+/-0,091

+0,002+/-0,003

+341,837+/-5,898

+49,881+/-0,942

Земля

+90,196+/-0,023

-558,944+/-11,995

+0,003+/-0,006

+228,374+/-1,460

Марс

+2,467+/-0,003

+74,696+/-6,729

+99,476+/-3,495

+0,002+/-0,007

Юпитер

+153,545+/-0,174

+674,666+/-11,611

+705,877+/-15,863

+1247,544+/-4,50

Сатурн

+7,322+/-0,058

+8,184+/-0,936

+19,361+/-1,590

+66,301+/-0,863

Уран

+0,152+/-0,009

+0,236+/-0,061

+0,562+/-0,019

+1,082+/-0,081

Нептун

+0,036+/-0,004

+0,047+/-0,175

+0,078+/-0,062

+0,493+/-0,030

Плутон

+0,006+/-0,002

+0,002+/-0,002

+0,003+/-0,006

+0,002+/-0,007

итого

+529,680+/-0,376

+50,944+/-31,716

+1153,612+/-27,04

+1594,442+/-7,92

 

 

Таким образом, логично сделать вывод, что у него и при расчете смещения перигелиев других планет от воздействия на них других планет по его формуле, которая даст результат как по формуле Ньютона, будут те же данные, что и у меня и только будет добавка от давления поля создаваемого Солнцем по Герберу.  Но, нет. Здесь, почему-то, данные по воздействию планет на другую исследуемую планету получились у Хайдарова очень отличающимися от моих данных. Вот, например, Земля при расчете ее воздействия на Венеру по закону тяготения Ньютона смещает перигелий Венеры на – 559 угл.сек, а у Хайдарова она смещает перигелий Венеры на + 892 угл.сек. Или, например, для Земли смещение перигелия от Венеры по законам Ньютона будет +342 угл.сек, а у Хайдарова получается -299 угл.сек. Таким образом, мы видим, что динамическое давление у Хайдарова просто огромное, т.к. оно даже изменяет знак смещений в рассмотренных примерах.

 

 

Но почему же оно у него тогда никак не проявляется при воздействии планет на Меркурий. Ведь это воздействие по законам Ньютона вполне приличное и Венера с Землей смещают его перигелий, соответственно, на 276 и 90 угл.сек. К тому же скорость у Меркурия самая большая из всех планет и он наиболее чувствителен к эффекту динамического давления. Ответа на эти вопросы нет и, следовательно, все данные, приведенные Хайдаровым, мягко говоря, сомнительные, хотя может быть в его теории эфира и есть что то здравое, а вот в расчете динамического давления гравитации мы здравого смысла не видим. Поэтому, данные, приведенные Хайдаровым, нет даже смысла проверять и, следовательно, отсутствие у меня его формулы для расчета силы притяжения между планетами никак не может повлиять на мои исследования по поиску скорости гравитации с использованием данных по смещениям параметров орбит планет. Следует также заметить, что я не знаю откуда Хайдаров взял такие данные НАСА, т.к. их данные, выложенные вот здесь  http://ssd.jpl.nasa.gov/txt/p_elem_t1.txt немного отличаются от того, что привел он.

 

 

 

Учет эффекта динамического давления поля в формуле учитывающей запаздывание потенциалов по координатам 

 

А сейчас давайте подведем предварительные итоги по сделанному обзору различных потенциалов и наметим пути дальнейшего исследования. Вообще то, я идею о запаздывании потенциала понимаю именно как его запаздывание по координатам, а не как возникающее при этом динамическое давление, но не отрицаю и возможности существования такого давления. И хотя этот эффект и будет второстепенным, я все же предусматриваю в своей математической модели Солнечной системы использование и этого эффекта. Но вот только как его учесть не ясно, т.к. предварительный теоретический анализ показал, что формулы Вебера, Гербера, Лиенара-Вихерта и Хайдарова для этого явно не подходят, но я все же продолжу исследование этих потенциалов уже практически. Поэтому, давайте попробуем сами получить хоть какую-то более-менее правдоподобную формулу для учета этого эффекта, чтобы посмотреть, как он себя проявит в изменении параметров орбит планет. А здесь мы с Вами в самом начале сталкиваемся с большой проблемой, т.к. мы даже не знаем какова природа гравитации – полевая, волновая или корпускулярная. Ведь если мы имеем дело с потоковыми явлениями (полевая и корпускулярная природа гравитации), то можно применять аналогии, например, с давлением ветра в паруса, а если волновая, то надо искать аналогии в радиационном давлении света и звука на препятствия. При этом, если радиационное давление света на препятствия ничтожно мало по сравнению с мощностью светового луча, то в звуковых волнах оно может быть уже заметной величиной.

 

 

Например, в [20] приводятся такие данные «в поле стоячей волны при частоте f=20 кГц и интенсивности звука 10^-4 Вт/м^2 в условиях земной гравитации может левитировать сфера с радиусом 0,4 см и массой =2 г». Да, здесь речь идет об акустических левитаторах, где используются стоячие волны, которые заставляют частицы занимать строго определенные положения или в узлах скорости или в пучностях и такой эффект, больше подходит для аналогии размещения атомов в веществе в строго заданных положениях, но и при бегущей волне тоже будет довольно заметная сила, чтобы ее легко зафиксировать. Опять таки, в работе [20] приводятся такие данные «для интенсивности звука в воде 10^-3 Вт/м^2 амплитуда звукового давления p=5*10^3 Па, тогда как радиационное давление на полностью поглощающую стенку составляет всего  p=10^2 Па». Таким образом, аналогии с давлением ветра в паруса и с давлением электромагнитных и звуковых волн говорят о том, что динамическая составляющая гравитационного взаимодействия вполне возможна и вопрос только в том, как ее рассчитать.

 

 

При этом ясно только то, что скорость распространения гравитационного взаимодействия должна быть конечной величиной. Так же я предполагаю, что должен быть и какой-то носитель этого взаимодействия между телами, т.е. здесь я четко ориентируюсь на наличие эфира или чего-то подобного, что должно переносить гравитационное взаимодействие, но меня пока совершенно не интересует каковы свойства этого эфира. Для меня важно, что есть какая-то среда, в которой распространяется гравитация. А для обоснования своей формулы для этого давления я сейчас попробую применить свой мощностной подход для описания явлений Природы [67]. Вот только, применение моего мощностного подхода обязательное при наличие проскальзывания при взаимодействие двух тел не вызывает затруднений, когда у нас имеется и непосредственный контакт взаимодействующих тел, а при гравитационном взаимодействии все гораздо сложнее, т.к. здесь у нас и гравитационное поле будет со скоростью Vgr проскальзывать относительно неподвижного тела k и это тело, при наличии скорости Vk, будет с этой скоростью дополнительно проскальзывать относительно движущегося поля.

 

 

Рассмотрение этого вопроса значительно упростилось бы, если бы мы рассмотрели взаимодействие гравитации с каким-то неподвижным объектом, т.е. с эфиром, а потом уже рассмотрели взаимодействие пробного тела движущегося в этом эфире, свойства которого искажены прошедшими через него гравитационными полями или волнами. Но такой подход требует создания еще и теории эфира, а это в мои планы, при нахождении сейчас скорости гравитации, не входит. Поэтому, никаких теорий я сейчас создавать не буду, а только на качественном уровне, используя методологию мощностного подхода, попробую создать какую ни будь правдоподобную формулу для учета эффекта динамического давления гравитации и при этом я буду стараться делать так, чтобы логика формулы не противоречила ни потоковому типу взаимодействия, ни волновому. Но, даже такую более-менее правдоподобную формулу для учета эффекта динамического давления гравитации, будет создать не просто, т.к. мой мощностной подход требует обязательного движения тела, даже если внешнее воздействие проскальзывает относительно него. Так при определении силы Fik = Nik / Vk, действующей от тела i на тело k и где Nik это мощность переданная от тела i телу k, скорость Vk стоит в знаменателе и не может быть равна нулю, а при гравитационном взаимодействие у нас должна быть сила, действующая на тело k, даже при скорости Vk равной нулю. Правда, здесь можно применить принцип виртуальных скоростей, т.к. тело покоится и значит это будет уже задача статики, а не динамики, но это еще больше запутает решение нашей задачи при тех наших познаниях о Природе, что у нас есть.

 

 

Поэтому, выводя формулу для учета динамического давления гравитации, я буду только на качественном уровне применять свой мощностной подход и при этом аналогию я буду проводить не с притяжением тел друг к другу, а с отталкиванием. В таком случае, если предположить, что скорости взаимодействующих тел влияют на силу их притяжения, то при этом должна изменяться и мощность воздействия одного тела на другое. При абсолютно неподвижных телах мощность воздействия будет такова, что мы получим силу воздействия Fn, которая определиться по формуле закона тяготения Ньютона (1) и единичная сила воздействия Fv, т.е. сила при единичной скорости взаимодействия определится по формуле (23). Тогда, если пробное тело покоится, то скорость прохождения поля сквозь пробное тело будет равна скорости гравитации Vgr и, действующая на него сила будет равна Fv*Vgr = Fn, а если пробное тело движется со скоростью Vk, которая положительна при увеличении расстояния между телами, то потребляемая им мощность будет больше при движение навстречу телу i и меньше при удалении от тела i и поэтому силу притяжения можно выразить формулой (24).

 

FvFn / Vgr           (23)

 

F = Fv*(VgrVk)  =  Fn*(VgrVk) /Vgr = Fn*(1 – Vk/Vgr)       (24)

 

F(i,k) = Fn(i,k)*(1 – Vk/Vgr) / (1 – Vi/Vgr)      (25)

 

Здесь получается, что чем чаще, например, гравитационные волны или гравитоны будут воздействовать на тело, тем будет больше сила притяжения (отталкивания), а при удалении k-го тела со скоростью равной скорости гравитации сила притяжения будет равна нулю и получается вроде бы все логично при потоковом походе, т.е. если рассматривать этот процесс, как давление ветра в паруса. Но, если природа гравитации не потоковая, а волновая, то перед Vk/Vgr должен стоять маленький коэффициент, который позволяет изменяться силе притяжения в статике (1) только на небольшую величину от этого значения. И потом, по формуле (24) получается, что если Vk=Vgr, то силы притяжения не будет, а, если Vk > Vgr, то получится уже сила отталкивания, но по идее такого не должно быть, если тело будет взаимодействовать с эфиром, который все равно уже был до этого искажен гравитационным полем источника в тот момент, как пробное тело попало в этот участок пространства. Да, даже, если рассматривать взаимодействие пробного тела с неискаженным эфиром, а непосредственно с полем, то и в этом случае оно будет присутствовать в этой части пространства. Только добиралось оно в эту точку пространства гораздо дольше, чем пробное тело. Так что пока у меня самого много вопросов к моей формуле (24), но я считаю, что на качественном уровне с ней можно работать.

 

 

А вот вопрос как скорость тела создающего это поле, повлияет на мощность этого поля, распространяемую в разные стороны, для меня остается пока совсем не ясным. Вроде бы, логично предположить, что, если гравитация объясняется действием гравитонов, как у Ритца, то, чем больше скорость, с которой тело i приближается к покоящемуся телу k, тем больше будет скорость гравитонов и тем больше будет мощность поля в этом направлении. Аналогично, если гравитация носит волновой характер, то тогда в направление к покоящемуся телу k волны будут распространяться с большей частотой и логично предположить, что и мощность в этом направлении будет больше. Тогда с учетом этой поправки формула (24) запишется как (25). Здесь у меня Fn(i,k) это сила притяжения i-ой массой k-ой массы, полученная по формуле Ньютона с учетом запаздывания потенциала по координатам согласно рис.1 и рис.3, а Vi и Vk это проекции скоростей i-ой и k-ой масс на вектор скорости гравитации Vgr при ее распространении от i-ой массы к k-ой массе. При этом, если проекции скоростей i-ой и k-ой масс совпадают с направлением вектора Vgr, то они будут положительные, а если направлены в другую сторону, то отрицательные. Необходимо также отметить, что скорость Vk будет соответствовать текущему моменту времени, а скорость Vi будет соответствовать моменту времени, когда гравитация начала распространяться от i-ой массы к k-ой массе. При этом в формулах Вебера и Гербера dR/dt = VkVi это тоже будет векторная разность проекций скоростей приемника и источника на вектор распространения Vgr.

 

 

Сомневаясь в своей формуле (24) я тем более не собираюсь доказывать, что и полученная мною формула (25) объективно отражает гравитационное взаимодействие между телами и у меня остаются большие сомнения по поводу ее справедливости и особенно в наличии в ней знаменателя. Ведь, например, скорость света в абсолютной системе отсчета, быстрее всего, остается неизменной при любой скорости источника света. Да и скорость звука в воздухе остается постоянной вне зависимости от того в какую сторону движется источник звука и тут аналогия с потоковым подходом для знаменателя формулы пропадает. Поэтому, я сделаю расчеты и по формуле (25) и по формуле (24), т.е. по формуле (25), когда в знаменателе там будет просто единица, т.е. не буду учитывать скорость источника. Более того, я в своей программе Solsys7m буду использовать для расчета давления гравитации формулу (25), где у меня перед Vk/Vgr и Vi/Vgr будут стоять коэффициенты k1 и k2, как это видно на рис. 4. При этом, когда k1=1 и k2=1, то это будет явно потоковый подход, а если мы уменьшим эти коэффициенты на несколько порядков, то это будет соответствовать волновому подходу, где мощность радиационного давления значительно меньше мощности самого излучения, а при k2=0 у нас формула (25) переходит в (24) и при k1=0 и k2=0 у нас вообще не будет динамического давления, а будет учитываться только запаздывание потенциалов по координатам.

 

 

А, если мы внимательно присмотримся к полученной нами формуле (25), то мы легко угадаем в ней формулу для продольного эффекта Доплера. А это уже интересно и, наверное, стоит более подробно ознакомиться с эффектом Доплера в общем виде на предмет обнаружения каких-то новых эффектов, которые мы могли пропустить при выводе формулы (25). Сейчас для эффекта Доплера в общем виде, а не для частного случая продольного эффекта, который рассмотрел сам Доплер, известно очень много формул и самая известная из них это формула Лоренца (26), где v0 - частота передатчика на движущемся источнике сигнала, а v - частота сигнала, принимаемая на движущемся приемнике. При этом b1=V1/c и b2=V2/c, V1 это скорость приемника сигнала и V2 это скорость источника сигнала, а Q1 и Q2 это углы между этими скоростями и радиус-вектором, соединяющим источник и приемник. Эту формулу частенько критикуют за то, что она дает неверные результаты, но если в этой формуле учесть запаздывание по координатам, то ошибки не будет. Таким образом, надо использовать углы Q1 и Q2 между векторами скорости приемника и источника и радиус-вектором распространения сигнала из той точки, где был источник, когда сигнал его покинул, в точку, где его принял приемник.

 

v =v0*(1 – b1*cos(Q1)) / (1 – b2*cos(Q2))                  (26)

 

v =v0*(1 – b*cos(Q)) / sqrt(1 – b^2)                              (27)

 

 

А вот согласно теории относительности не может быть абсолютной скорости и есть только относительные скорости, поэтому формула (26) упрощается в формулу (27) и в ней присутствует только коэффициент b, который равен скорости приемника относительно источника деленной на скорость света. Если использовать только числитель этой формулы, то получающийся результат будет незначительно отличаться от того, что дает формула Лоренца с моей поправкой, но здесь появился дополнительный релятивистский множитель, который дает, так называемый, поперечный эффект Доплера, которого нет в классической формуле (26). Но, альтернативные ученые пытаются доказать, что такой эффект существует и в их классических формулах для эффекта Доплера. Я же в их формулах кроме ошибок никакого поперечного эффекта Доплера не нашел (см. мою статью [68]). А откуда же он тогда взялся в  формуле (27). Оказывается, эта формула учитывает, согласно СТО, замедление времени при увеличении относительной скорости источника сигнала и в результате уменьшается частота передатчика, расположенного на источнике. Причем, она уменьшается и при движении источника вдоль линии соединяющей источник и приемник и в поперечном направлении, т.е. этот эффект никак не является именно поперечным и здесь надо говорить о другом названии для этого эффекта. И по аналогии с гравитационным красным смещением, которое дает ОТО, его можно назвать скоростным красным смещением СТО, т.к. и там и там этот эффект уменьшения частоты передатчика расположенного на источнике объясняется замедлением времени и эти эффекты не зависят от направления движения источника.

 

 

Таким образом, релятивистская формула для эффекта Доплера (27) с методологической точки зрения явно ошибочна, т.к. должна отражать оба релятивистских эффекта замедления времени. Более того, в учебниках по ТО дается и другое объяснение гравитационного красного смещения, которое объясняет сдвиг спектра в красную сторону, т.е. уменьшение частоты излучения, исходя из корпускулярной природы света. Причем эти два объяснения считаются равнозначными (наверное, потому, что результат получается один и тот же), но ведь это два совершенно разных физических эффекта (если их можно назвать физическими) и если оба эти объяснения верны, то учитывать надо оба этих эффекта, т.е. результат от гравитационного красного смещения согласно ТО должен быть в два раза больше. Но, даже, если сторонники ТО договорятся между собою какова должна быть формула для эффекта Доплера с учетом всех релятивистских эффектов, формула (27) в наше случае для учета динамического давления гравитации уже сейчас означает дополнительное увеличение силы притяжения между массами как при их удалении друг от друга, так и при их приближении, т.к. относительная скорость b в знаменателе будет в квадрате. А это для нашего случая является явным нонсенсом и поэтому эта формула нам явно не подходит. Ну и, естественно, эта формула нам не подходит потому, что здесь используется только  угол Q отражающий текущие положения приемника и источника, т.е. формула не отражает запаздывание по координатам. К тому же, у меня возникло много вопросов к этой формуле и чисто по эффекту Доплера, которые я изложил в своей статье [68], поэтому я от нее и отказался.

 

 

Но вот, как я недавно узнал, идея с эффектом Доплера, т.е. с формулой аналогичной моей формуле (25), но якобы учитывающая и поперечный эффект Доплера,  использовалась авторами работы [12] для объяснения аномального смещения перигелия Меркурия. Формула у них получилась очень странная и очень длинная, чтобы ее воспроизвести, это формула (30*) в текущих координатах, а с учетом запаздывания координат источника формула (31*), которую они и рекомендуют для учета давления гравитационного поля как формулу (66*). А, т.к. моя формула аналогична формуле для продольного эффекта Доплера, то необходимо разобраться более подробно почему у Сухоруковых получилась формула (66*) и возможно, что поперечный эффект Доплера все же существует, а я не учел этот эффект при выводе своей формулы (25). Главной формулой у Сухоруковых, от которой начинаются все преобразования, является формула для расчета скорости сближения волн с приемником. Эта формула в их работе находиться под их формулой (29*). К сожалению, эта формула не только не пронумерована, но авторы и совершенно не объясняют из какого рукава они ее достали. Но в своей статье по эффекту Доплера я все же выяснил откуда у них взялась эта формула и в результате показал, что она ошибочна. Таким образом, получается, что я, быстрее всего, рассуждал правильно, используя свой мощностной подход, для вычисления динамического давления гравитационного поля, когда выводил формулу (25), т.к. формула получилась аналогичной формуле Лоренца для эффекта Доплера (с моими поправками), а это хоть и не доказывает ее правильность, но позволяет надеяться на это.

 

 

 

Практическая проверка различных формул, учитывающих конечность скорости гравитации, на примере Солнечной системы.

 

А теперь давайте посмотрим, что получается по моей формуле (25), а заодно и проверим формулы Вебера, Гербера и Лиенара с Вихертом. При этом, хотя формула Сухоруковых [12] явно ошибочна, давайте проверим и, что получается у них с использованием их формулы (31*), которая аналогична их формуле (66*), и где вместо частоты у них используются силы. А для сравнения я приведу и данные, которые получаются с использованием теории Ньютона и ОТО. Ну, а начнем мы с практического исследования формул Вебера, Гербера, моей и Сухоруковых, т.е. с тех формул, которые воспроизводила программа Solsys7mm и которые я рассматривал в предыдущих редакциях статьи. А потом я в этой третьей редакции статьи приведу и данные, которые я получил на программе Solsys7mmm, где я добавил рассмотрение формул для потенциалов Лиенара-Вихерта и ОТО, а также повторю данные по своей формуле.

 

Хотя формально у Вебера и Гербера учитывается скорость распространения гравитации, но получается, что только для учета динамического давления, и при этом координаты самих планет у них используются только текущие и расстояния между планетами и скорости, используемые в расчетах, вычисляются относительные, т.е. одной планеты относительно другой. И при этом у них силы, действующие между телами, направлены по одной прямой соединяющей координаты тел на текущий момент времени, т.е. здесь принимается, что скорость гравитации равна бесконечности и получается, что нет запаздывания по координатам, но для Меркурия это не критично. Ведь запаздывания по координатам от силы притяжения Солнцем Меркурия практически не будет, т.к. Солнце все равно будет стоять практически на месте, а запаздывание от силы притяжения Меркурием Солнца можно не принимать в расчет, т.к. все равно он Солнце значительно не сдвинет с места и, следовательно, расчеты, выполненные как с запаздыванием по координатам, так и без запаздывания практически совпадут. То же самое относится и к формуле Сухоруковых, но я выполню расчет по формулам Вебера, Гербера и Сухоруковых и с учетом запаздывания по координатам, используя при этом расчет самого запаздывания по своим формулам.

 

 

Ну и заодно давайте и протестируем программу Solsys7m на правильность учета сил давления гравитационного поля при его скорости равной скорости света, т.к. все приведенные мною авторы использовали именно эту скорость, и сравним результаты, полученные мною при численном решении дифференциальных уравнений, с аналитическими решениями, полученными различными авторами. При этом, то, как мне пришлось модернизировать метод Рунге-Кутта, чтобы с его применением решить дифференциальные уравнения, где у нас используются силы Вебера и Гербера, которые у них зависят и от ускорений, т.к. в стандартном методе Рунге-Кутта не предусмотрен расчет ускорений по ускорениям, вы можете посмотреть в приложении 2. Результаты тестирования я привожу в таблице 4, где данные получены при шаге численного решения уравнений равном 3600 сек, с уменьшением его в 1000 раз вблизи положений перигелия (афелия) и восходящего (нисходящего) узла орбиты, чтобы точнее зафиксировать их оскулирующие положения.  В таблице смещения перигелия даны в угловых секундах за 100 лет и приняты следующие обозначения - 1pl- расчет ведется только для взаимодействия Меркурия с Солнцем, 9pl – расчет ведется при взаимодействии Меркурия и с Солнцем и с 8-ю другими планетами, +Pgr – давление гравитационного поля учитывается, а -Pgr – не учитывается, +DX - учитывается запаздывание по координатам, по моим формулам, которые я приводил выше (вариант расчета DX1), а -DX  - не учитывается запаздывание по координатам. При этом в колонках Гербер1 и Гербер2 будут приведены расчеты по формулам (6-2) и (6-2’), а Юдин1 и Юдин2 данные при расчете по формуле (25), когда, соответственно, k1=1, k2=1, т.е. используются и числитель и знаменатель формулы, и k1=1, k2=0, т.е. используется только числитель формулы. А в двух последних строках здесь отражен только эффект от динамического давления при наличии или отсутстствии запаздывания по координатам в системах с 9-ю планетами, т.е. приведена разница между получающимся значением с заданными условиями и значением при таких же условиях, но без учета динамического давления, т.е. при -Pgr.

 

 

Таблица 4. Смещение перигелия Меркурия в покоящейся Солнечной системе и при скорости гравитации равной скорости света.

 

Варианты расчета

Вебер

Гербер1

Гербер2

Юдин1

Юдин2

Сухоруков

1pl   +Pgr  -DX

7,1684

42,9926

42,9937

1,7742

1,7777

8,7294

1pl   +Pgr +DX

7,1685

42,9779

42,9781

1,7295

1,7711

8,7346

9pl +Pgr  -DX

536,8530

572,6796

572,6795

618,0663

654,7382

625,0600

9pl +Pgr +DX

536,3834

572,1960

572,1961

614,2655

650,8803

621,2551

9pl +Pgr  -DX      -529,7076 (-Pgr)

7,1454

42,9720

42,9719

88,3587

125,0306

95,3524

9pl +Pgr +DX      -529,2340 (-Pgr)

7,1497

42,9620

42,9621

85,0315

121,6463

92,0211

 

 

Как видим по данным Вебера и Гербера программа Solsys7m правильно определила смещения перигелия Меркурия, когда рассматривается его взаимодействие только с Солнцем, т.к. по данным аналитического решения, которые приводил Роузвер, у Вебера должно быть 7,2 угловых секунд, а сам Гербер приводит значение 41 угл.сек (для полученной им скорости света 305,5 тыс.км/с). В общем, можно считать, что наше численное решение почти точно совпадает с аналитическим. Здесь, правда, еще не ясно, что с чем совпадает, т.к. аналитические решения могут быть в разы менее точные, чем численные, да и у Гербера явно приблизительное решение. Этот вопрос я подробно рассмотрел в статье [67] и останавливаться на нем не буду, а сейчас просто будем считать, что приведенные данные аналитического решения довольно таки точные и поэтому моя программа правильно учитывает эффект динамического давления при гравитационном взаимодействии, а из приведенных в таблице данных можно сделать и некоторые выводы.

 

 

Первое это то, что на результат по формулам Вебера и Гербера в покоящейся Солнечной системе совершенно не оказывает влияния то, сколько у нас планет в системе и имеется ли при этом запаздывание потенциала по координатам или нет, т.к. в двух последних строчках таблицы я привожу смещение перигелия, которое получается при вычете из общего эффекта (от статического взаимодействия планет между собой и динамического давления) этого же эффекта, но без учета динамического давления, как при отсутствие запаздывания по координатам, так и при наличии. А в результате получаем те же значения смещений, как если бы у нас в системе были только Солнце и Меркурий. Данный результат является явно противоестественным, что позволяет сделать вывод о неадекватности формул Вебера и Гербера. А второе, это то, что моя формула и формула Сухоруковых оказываются чувствительны, как к количеству планет в системе, так и к наличию или отсутствию запаздывания по координатам, хотя у самих Сухоруковых этого запаздывания нет, что позволяет надеяться на то, что они адекватно отражают реальность, хотя наблюдательным данным пока больше соответствуют данные Гербера. А теперь давайте посмотрим, как учет по этим формулам динамического давления гравитации, отразится на вековых смещениях других параметров орбиты Меркурия. Ведь не должно же быть такого, чтобы этот эффект влиял на смещения перигелия dAlfaP1 и не отражался на смещениях узла восхождения dAlfaU1, угла наклона плоскости орбиты к эклиптике dBetta1, эксцентриситета dEks1, большой полуоси эллипса dR1 и угловой скорости dW1. Здесь и далее  dAlfaP, dAlfaU и dBetta в угл.сек за век, dEks в абсолютных единицах за век увеличенное в 1 000 000 раз, dRsr в тыс.км. за век и dW в рад/век за век. И в каждой третьей строке здесь отражен только эффект от динамического давления при заданных условиях.

 

 

Таблица 5. Смещение всех параметров орбиты Меркурия в покоящейся Солнечной системе и при скорости гравитации равной скорости света.

 

Варианты расчета

Вебер

Гербер1

Гербер2

Юдин1

Юдин2

Сухоруков

1pl +Pgr  +DX                              dAlfaP1

7,1685

42,9779

42,9781

1,7295

1,7711

8,7346

9pl +Pgr +DX

536,3834

572,1960

572,1961

614,2655

650,8803

621,2551

9pl +Pgr +DX      -529,2340 (-Pgr)

7,1497

42,9620

42,9621

85,0315

121,6463

92,0211

1pl +Pgr  +DX                             dAlfaU1

+0,0044

-0,0101

-0,0092

+0,0204

-0,0251

+0,020

9pl +Pgr +DX

-451,4799

-451,4746

-451,4728

-506,5038

-506,8767

-506,5227

9pl +Pgr +DX      +451,4677 (-Pgr)

-0,0122

-0,0069

-0,0051

-55,0361

-55,4090

-55,0550

1pl +Pgr  +DX                              dBetta1

+0,0005

+0,0021

+0,0021

-0,0002

-0,0002

-0,0002

9pl +Pgr +DX

-21,4717

-21,4704

-21,4704

-29,9754

-29,9116

-29,9761

9pl +Pgr +DX      +21,4720 (-Pgr)

+0,0003

+0,0016

+0,0016

-8,5034

-8,4396

-8,5041

1pl +Pgr  +DX                                 dEks1

+0,0091

-0,0213

-0,0216

49902,1632

49902,1370

49902,1649

9pl +Pgr +DX

20,9037

20,8731

20,8731

49815,0833

49862,0270

49815,0858

9pl +Pgr +DX      -20,9098 (-Pgr)

-0,0061

-0,0367

-0,0367

49794,1735

49841,1172

49794,1760

1pl +Pgr  +DX                                   dR1

0,0164

0,0156

0,0156

1510,7438

1510,7430

1510,7439

9pl +Pgr +DX

1,1555

1,1549

1,1549

1513,5620

1515,4341

1513,5617

9pl +Pgr +DX      -1,1557 (-Pgr)

-0,0002

-0,0008

-0,0008

1512,4063

1514,2784

1512,4060

1pl +Pgr  +DX                                   dW1

-0,0011

-0,0011

-0,0011

-97,5800

-97,5799

-97,5800

9pl +Pgr +DX

-0,0781

-0,0780

-0,0780

-97,7554

-97,8708

-97,7554

9pl +Pgr +DX      +0,0781 (-Pgr)

0,0000

+0,0001

+0,0001

-97,6773

-97,7927

-96,6773

 

 

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

 

 

Но давайте все же попробуем «растормошить» формулы Вебера и Гербера, чтобы они хоть как-то проявили себя в изменение параметров орбиты Меркурия. Для этого давайте проведем вычислительные эксперименты не в покоящейся Солнечной системе, а в движущейся со скоростью 100 км/с по оси Y. При этом в таблице 6-1 не будем учитывать запаздывание по координатам, а в таблице 6-2 будем. Здесь, так же, как и в таблице 5, в первой строчке приводится результат, получающийся при учете динамического давления, а во второй разность при вычете результата, получающегося при тех же условиях, но без учета динамического давления, т.е. в каждой второй строке отражен только эффект от динамического давления. А вот данные Юдин2 я здесь не привожу, т.к. вследствие очень маленькой скорости гравитации (равной скорости света) получаются очень не стабильные результаты, по которым нельзя определить смещения параметров орбит.

 

 

Таблица 6-1. Смещение всех параметров орбиты Меркурия в Солнечной системе движущейся со скоростью VYsys=100 км/с при отсутствии запаздывания по координатам и при скорости гравитации равной скорости света.

 

Варианты расчета

Вебер

Гербер1

Гербер2

Юдин1

Юдин2

Сухоруков

9pl +Pgr -DX                              dAlfaP1

536,94

572,75

572,78

618,30

-

632,73

9pl +Pgr -DX      -529,79 (-Pgr)

7,15

42,96

42,99

88,51

-

102,94

9pl +Pgr -DX                             dAlfaU1

-451,37

-451,41

-451,42

-506,26

-

-506,75

9pl +Pgr -DX      +451,39 (-Pgr)

+0,02

-0,02

-0,03

-54,87

-

-55,36

9pl +Pgr -DX                              dBetta1

-21,45

-21,45

-21,45

-30,03

-

-29,94

9pl +Pgr -DX      +21,45 (-Pgr)

0,00

0,00

0,00

-8,59

-

-8,49

9pl +Pgr -DX                                 dEks1

20,46

20,43

20,43

49823,98

-

49895,41

9pl +Pgr -DX      -20,46 (-Pgr)

0,00

-0,03

-0,03

49803,52

-

49874,95

9pl +Pgr -DX                                   dR1

0,00

0,00

0,00

1512,52

-

1515,32

9pl +Pgr -DX      -0,00 (-Pgr)

0,00

0,00

0,00

1512,52

-

1515,32

9pl +Pgr -DX                                   dW1

0,00

0,00

0,00

-97,69

-

-97,86

9pl +Pgr -DX      -0,00 (-Pgr)

0,00

0,00

0,00

-97,69

-

-97,86

 

 

Таблица 6-2. Смещение всех параметров орбиты Меркурия в Солнечной системе движущейся со скоростью VYsys=100 км/с при наличии запаздывания по координатам и при скорости гравитации равной скорости света.

 

Варианты расчета

Вебер

Гербер1

Гербер2

Юдин1

Юдин2

Сухоруков

9pl +Pgr +DX                              dAlfaP1

4997,54

5030,97

5031,29

5598,88

-

5615,42

9pl +Pgr +DX      -4990,88 (-Pgr)

6,66

40,09

40,41

608,00

-

624,54

9pl +Pgr +DX                             dAlfaU1

-6886,86

-6888,26

-6889,68

-7914,31

-

-7917,05

9pl +Pgr +DX      +6886,56 (-Pgr)

-0,30

-1,70

-3,12

-1027,75

-

-1030,49

9pl +Pgr +DX                              dBetta1

-1278,79

-1278,72

-1278,40

-1458,33

-

-1458,70

9pl +Pgr +DX      +1278,81 (-Pgr)

+0,02

+0,09

+0,41

-179,72

-

-179,89

9pl +Pgr +DX                                 dEks1

928,46

928,03

927,87

47911,41

-

47979,92

9pl +Pgr +DX      -928,55 (-Pgr)

-0,09

-0,52

-0,68

46982,86

-

47051,37

9pl +Pgr +DX                                   dR1

2159,24

2157,96

2157,25

3947,09

-

3949,52

9pl +Pgr +DX      -2159,49 (-Pgr)

-0,25

-1,53

-2,24

1787,60

-

1790,03

9pl +Pgr +DX                                   dW1

-136,30

-136,22

-136,15

-237,56

-

-237,69

9pl +Pgr +DX      +136,32 (-Pgr)

+0,02

+0,10

+0,17

-101,24

-

-101,37

 

 

Как видим, и это воздействие на Солнечную систему по формулам Вебера и Гербера никак не отражается, т.е. эти формулы являются «мертвыми» формулами никак не связанными с реальностью. А вот моя формула и формула Сухорукова чувствительны к любым изменениям условий функционирования системы и в этом смысле они являются «живыми» формулами, которые отражают эффект динамического давления гравитации. Правда, здесь у нас, даже без учета запаздывания по координатам (см. табл. 6-1), получаются просто сумасшедшие значения смещений эксцентриситета, большой полуоси и угловой скорости обращения вокруг Солнца, чего в действительности точно не наблюдается, но судить об их ошибочности или вообще об отсутствии эффекта динамического давления пока рано. 

 

 

Например, в таблицах 1-11-X … 1-61-Z, которые я даю в приложении 3, у меня приведены данные по смещению параметров орбиты Меркурия, которые получаются при учете запаздывания по координатам и динамического давления по моей формуле (25) с k1=1 и k2=1 при различных скоростях Солнечной системы и гравитации. Две из этих таблиц приведены ниже и из них следует, что уже при скорости гравитации равной 1000 скоростей света (n=3) у нас смещения перигелия при разных скоростях Солнечной системы отличаются не значительно, а смещения большой полуоси, хоть и продолжают различаться, но значения получаются уже не такими большими, как были при скорости гравитации равной скорости света. А в таком случае, если скорость гравитации окажется значительно больше 1000 скоростей света, то и смещения эксцентриситета, большой полуоси и угловой скорости будут вполне приемлемыми для длительного существования Солнечной системы. Но, если у нас скорость гравитации получиться больше 1000 скоростей света, то тогда мы, возможно, не сможем объяснить смещение параметров орбит планет эффектом запаздывания по координатам, т.к. из данных, приведенных в приложении 1, видно, что в этом случае этот эффект становиться еле заметным. Поэтому, если эффект динамического давления и имеет место быть, то быстрее всего это давление носит не потоковый характер по аналогии с давлением ветра в паруса, а волновой, т.е. имеется только радиационное давление, т.е. в формуле (25) значения k1 и k2 должны быть много меньше единицы.

 

 

Таблица 1-11-Y. AlfaP1-VYsys  - вековые смещения перигелия Меркурия (угл. cекунды)

 

n

-300

-200

-100

0

100

200

300

0

-

-11548,11

-5003,69

614,27

5598,88

10245,42

14617,83

1

-878,09

-404,21

+65,89

531,81

994,96

1472,98

1855,50

2

+392,06

+438,18

483,97

529,98

576,00

621,67

667,55

3

515,97

520,51

525,08

529,74

534,31

539,09

543,41

 

 

Таблица 1-51-X. Rsr1-VXsys  - вековые смещения большой полуоси Меркурия (тыс. км)

 

n

-300

-200

-100

0

100

200

300

0

32080,03

22778,72

12687,02

+1513,56*

-12608,47*

-

-

1

3261,48

2222,77

1173,38

+112,61

-960,03

-2045,56

-3145,01

2

325,95

221,05

116,06

+10,94

-94,29

-199,63

-305,06

3

32,59

22,10

11,59

+1,09

-9,41

-19,91

-30,41

 

 

 

А теперь давайте рассмотрим как повлияют на вековые смещения параметров орбит планет и потенциалы Лиенара-Вихерта, а также рассмотрим и принципиально другое описание гравитационного взаимодействия, которое отражено в ОТО. При этом я приведу и данные, которые получаются по моей формуле (25), но они будут немного отличаться от приведенных выше в лучшую сторону при движение всей солнечной Системы. Дело в том, что в программах Solsys7m и Solsys7mm у меня решение выполнялось в абсолютной системе отсчета и при наличии скорости всей Солнечной системы очень быстро увеличивались значения координат планет и Солнца в направлении этой скорости. В результате увеличивалось число цифр до запятой и уменьшалось после запятой, т.к. в сумме восьми битное представление данных, которое используется в программе, обеспечивает в сумме только 15 цифр. Но сейчас я в программе Solsys7mmm, кроме того, что добавил туда расчеты с использованием потенциалов Лиенара-Вихерта и по уравнениям ОТО, изменил действие чекбоксов <Xsol=0> и <Xzem=0> на 2-ой форме программы, где выполняется моделирование, а также добавил и чекбокс <Xцм=0>. Если раньше при отмеченном одном из чекбоксов все координаты планет оставались неизменными, а только на рисунок выводились координаты планет относительно Солнца <Xsol=0>,  относительно Земли <Xzem=0> или относительно центра масс системы <Xцм=0>, что позволяло видеть относительные траектории планет при заданной скорости всей системы, и расчет всех параметров орбит велся или относительно Солнца или относительно Земли или относительно центра масс системы, то сейчас после каждого шага решения уже сами координаты планет пересчитываются в относительные, а координаты Солнца, Земли или центра масс системы принимаются равными нулю. Это позволяет значительно повысить точность вычислений при задании всей системе поступательной скорости и получать более стабильные данные по вековым смещениям параметров орбит планет.

 

 

Кроме этого я в программе Solsys7mmm добавил и еще один вариант расчета ускорений, которые используются в формулах Вебера, Гербера и Лиенара-Вихерта, т.к. стандартный метод Рунге-Кутта не рассчитан на такую неожиданность. А во всем остальном программа Solsys7mmm полностью соответствует программе Solsys7mm. При этом, как я писал выше, силы притяжения между объектами с использованием потенциалов Лиенара-Вихерта я вычислял с использованием формулы (5.17'), но предусмотрел два варианта расчета. В первом варианте расчет идет по полной формуле (5.17'), а во втором по укороченной, которая получается при взятии производной только от скалярного потенциала по Rf, т.е. получаем формулу Ньютона, но не с реальным радиусом, а с фиктивным, который вычисляется из запаздывающего радиуса Rf= R' – V' * R' / Vgr. И ниже приведен код программы, где отражены эти два варианта расчета. Если на форме 23 (см. рис. 4) отмечен вариант расчета по упрощенной формуле, то будет kodLV0 = 1 и расчет сил будет идти по последним трем уравнениям, а, если этот вариант расчета не отмечен, то расчет будет идти по полной формуле. При этом так же, как и для других формул, на форме 23 можно задать значения kPgr1 и kPgr2 в окошках k1 и k2, чтобы учесть или не учитывать воздействие ускорения.

 

 

Rfik(i, j) = R(i, j) - R(i, j) * dVi(i, j) / Vgr ' фиктивный радиус

F(i, j) = -kodF(i, j) * m(i) * m(j) * gamma / Rfik(i, j) ^ 3 / Vgr ^ 2 'находим условную силу притяжения между следом i-го тела воздействующим на j тело

 

If kodLV0 = 0 Then 'полная формула Лиенара-Вихерта

FX(i, j) = F(i, j) * ((DX(i, j) - R(i, j) * VXold(i, j) / Vgr) * (R(i, j) * kPgr1 * Ai(i, j) + Vgr ^ 2 - Vold(i, j) ^ 2) - R(i, j) * Rfik(i, j) * kPgr2 * AXzap(i, j))

FY(i, j) = F(i, j) * ((DY(i, j) - R(i, j) * VYold(i, j) / Vgr) * (R(i, j) * kPgr1 * Ai(i, j) + Vgr ^ 2 - Vold(i, j) ^ 2) - R(i, j) * Rfik(i, j) * kPgr2 * AYzap(i, j))

FZ(i, j) = F(i, j) * ((DZ(i, j) - R(i, j) * VZold(i, j) / Vgr) * (R(i, j) * kPgr1 * Ai(i, j) + Vgr ^ 2 - Vold(i, j) ^ 2) - R(i, j) * Rfik(i, j) * kPgr2 * AZzap(i, j))

End If

 

If kodLV0 = 1 Then 'упрощенная формула Лиенара-Вихерта

FX(i, j) = F(i, j) * Vgr ^ 2 * (DX(i, j) - R(i, j) * VXold(i, j) / Vgr)

FY(i, j) = F(i, j) * Vgr ^ 2 * (DY(i, j) - R(i, j) * VYold(i, j) / Vgr)

FZ(i, j) = F(i, j) * Vgr ^ 2 * (DZ(i, j) - R(i, j) * VZold(i, j) / Vgr)

End If

 

 

А формулы для расчета по ОТО я использовал те, которые получаются у Ландау из лагранжиана (106,17) [4] и в этом мне помог Дробышев на форуме sciteclibrary (http://www.sciteclibrary.ru/cgi-bin/yabb2/YaBB.pl?num=1400127682/700  ), где я обсуждал потенциалы Лиенара-Вихерта. Он сначала получил из этого лагранжиана дифференциальные уравнения в векторном виде, а потом расписал их мне в проекциях по трем осям координат. Кого интересуют получившиеся у него формулы я отсылаю на этот форум, а здесь приведу уже сам код программы для расчета по этим формулам.

 

 

 

Рис. 11. Скриншот сообщения Дробышева с форума SciTecLibrary.

 

 

k1 = kodF(i, j) * m(i) * gamma / R(i, j) ^ 3

k2 = gamma * (5 * m(j) + 4 * m(i)) / R(i, j) - VX(j) ^ 2 - VY(j) ^ 2 - VZ(j) ^ 2 - 2 * VX(i) ^ 2 - 2 * VY(i) ^ 2 - 2 * VZ(i) ^ 2 + 4 * VX(j) * VX(i) + 4 * VY(j) * VY(i) + 4 * VZ(j) * VZ(i) + 3 * (VX(i) * DX(i, j) + VY(i) * DY(i, j) + VZ(i) * DZ(i, j)) ^ 2 / 2 / R(i, j) ^ 2

k3 = (4 * VX(j) - 3 * VX(i)) * DX(i, j) + (4 * VY(j) - 3 * VY(i)) * DY(i, j) + (4 * VZ(j) - 3 * VZ(i)) * DZ(i, j)

 

FX(i, j) = -k1 * DX(i, j) + k1 / Vgr ^ 2 * (DX(i, j) * k2 + (VX(j) - VX(i)) * k3)

FY(i, j) = -k1 * DY(i, j) + k1 / Vgr ^ 2 * (DY(i, j) * k2 + (VY(j) - VY(i)) * k3)

FZ(i, j) = -k1 * DZ(i, j) + k1 / Vgr ^ 2 * (DZ(i, j) * k2 + (VZ(j) - VZ(i)) * k3)

 

 

И ниже в таблицах 7-1 и 7-2 я привожу данные, которые у меня получаются с использованием потенциалов Лиенара-Вихерта (Л-В) как по полной формуле, так и по упрощенной и по формуле ОТО. Причем данные расчета /ОТО Ландау/ это те, которые получились у меня по вышеприведенным формулам, а /ОТО JPL2/ это данные, которые получились в лаборатории реактивного движения (JPL) на их более сложной модели ОТО и которые я получил статистически обработав результаты их расчетов, которые они записали в эфемеридах DE 405. Также я повторяю и данные, которые получаются у меня при учете эффекта запаздывания потенциалов по координатам при расчете по формуле (25), как с учетом динамического давления (k1=1 k2=1), так и без учета (k1=0 k2=0), т.е. с учетом только запаздывания по координатам. А кроме этого я привожу и данные, которые получаются по теории Ньютона, а также для сравнения привожу здесь и наблюдательные данные по вековым смещениям параметров орбит, которые получились у Ньюкома и у меня при обработке данных оптических наблюдений за планетами и Солнцем. При этом, т.к. расчетные данные по большой полуоси и угловой скорости взаимосвязаны, я привожу только данные по большой полуоси, а наблюдаемые данные по этим параметрам отсутствуют, т.к. мы с Ньюкомом не смогли обнаружить их изменение. В таблице 7-1 приведены расчетные данные в покоящейся Солнечной системе, а в таблице 7-2 в движущейся со скоростью VY=100 км/с. При этом все расчеты выполнены, когда в системе присутствуют все 9 планет и при скорости гравитации равной скорости света, а данные даны при 95% доверительной вероятности, когда доверительный интервал составляет +/- 2 среднеквадратичных отклонения.

 

 

 

 

 

Таблица 7-1. Смещение параметров орбиты Меркурия наблюдаемые и расчетные в покоящейся Солнечной системе при скорости гравитации равной скорости света.

 

Варианты расчета

dAlfaP1

dAlfaU1

dBetta1

dEks1

dR1

Наблюдения Ньюком

570,73+/-3,90

-452,18+/-7,38

-21,43+/-1,60

20,55+/-4,84

-

Наблюдения Юдин

578,04+/-7,25

-433,15+/-8,21

-19,84+/-0,51

20,10+/-3,98

-

Ньютон

529,79+/-0,36

-451,40+/-0,02

-21,45+/-0,01

20,49+/-0,19

0

ОТО Ландау

572,76+/-0,32

-451,42+/-0,02

-21,45+/-0,01

20,49+/-0,19

-0,001+/-0,001

ОТО JPL2

572,20+/-0,07

-449,95+/-0,42

-21,44+/-0,00

20,50+/-0,06

0

Л-В полный

529,75+/-0,35

-451,40+/-0,02

-21,45+/-0,01

20,49+/-0,19

-0,001+/-0,001

Л-В скалярный

529,79+/-0,36

-451,41+/-0,02

-21,46+/-0,01

20,49+/-0,19

-0,001+/-0,001

Юдин k1=0 k2=0

529,30+/-0,30

-451,47+/-0,01

-21,47+/-0,01

20,91+/-0,18

1,155+/-0,004

Юдин k1=1 k2=1

614,50+/-0,05

-506,57+/-11,08

-29,99+/-2,08

49819,75+/-

1513,88+/-12,7

 

 

Таблица 7-2. Смещение параметров орбиты Меркурия наблюдаемые и расчетные в Солнечной системе движущейся со скоростью VY=100 км/с и при скорости гравитации равной скорости света.

 

Варианты расчета

dAlfaP1

dAlfaU1

dBetta1

dEks1

dR1

Наблюдения Ньюком

570,73+/-3,90

-452,18+/-7,38

-21,43+/-1,60

20,55+/-4,84

-

Наблюдения Юдин

578,04+/-7,25

-433,15+/-8,21

-19,84+/-0,51

20,10+/-3,98

-

Ньютон

529,79+/-0,36

-451,40+/-0,02

-21,45+/-0,01

20,49+/-0,19

0

ОТО Ландау

572,74+/-0,33

-451,42+/-0,02

-21,45+/-0,01

20,49+/-0,19

-0,001+/-0,001

ОТО JPL2

-

-

-

-

-

Л-В полный

529,80+/-0,34

-451,37+/-0,04

-21,45+/-0,01

20,48+/-0,19

-0,001+/-0,001

Л-В скалярный

529,81+/-0,32

-451,39+/-0,03

-21,45+/-0,01

20,48+/-0,19

-0,001+/-0,001

Юдин k1=0 k2=0

4991,2+/-57,6

-6887,9+/-113,1

-1278,5+/-24,9

928,4+/-13,8

2158,8+/-55,4

Юдин k1=1 k2=1

5600,6+/-138,3

-7919,9+/-413,7

-1458,8+/-20,6

47922,1+/-792

3948,9+/-133,5

 

 

Как видим, здесь ОТО, как и потенциалы Гербера, тоже дает нужное значение смещения перигелия Меркурия, но по всем остальным параметрам она дает те же результаты, что и теория Ньютона и никак не отражает то, что изменились условия функционирования системы, т.е. то, что система движется со скоростью 100 км/с. А вот потенциалы Лиенара-Вихерта, как по полной формуле, так и по сокращенной,  по всем параметрам дают тот же самый результат, что и расчеты по Ньютону, т.е. при тех скоростях и ускорениях, что у нас имеются у планет, потенциалы Лиенара-Вихерта полностью идентичны потенциалам Ньютона, где скорость гравитации равна бесконечности и нет никакого запаздывания потенциалов. Поэтому, по крайней мере для моего исследования, где я пытаюсь найти по наблюдаемым данным смещений параметров орбит планет скорость гравитации и абсолютную скорость Солнечной системы и потенциалы Лиенара-Вихерта и теория ОТО так же, как и потенциалы Вебера и Гербера, точно не подходят. А, что касается моих потенциалов запаздывающих по координатам и учитывающих динамическое давление, то можно констатировать, что они реагируют как на изменение скорости Солнечной системы так и на изменение скорости гравитации и есть надежда, что они адекватно отражают это воздействие.

 

 

Но, возможно, что при скоростях близких к скорости света и при больших ускорениях и потенциалы Лиенара-Вихерта будут давать отличные от теории Ньютона результаты, поэтому давайте проверим, что они дадут для двойных пульсаров, а заодно посмотрим и то, что дадут потенциалы Гербера и мои, а так же геометрическая теория Эйнштейна, т.е. ОТО. Для этого на форме 25 загрузим из файла параметры двойного пульсара PSR 1913+16 и рассчитаем необходимые нам для начала моделирования на форме 2 начальные данные, т.е. координаты и скорости пульсара и его компаньона. А, т.к. алгоритм поиска перигелиев и афелиев орбит планет на форме 2 удовлетворительно работает только при небольших углах наклона орбит, что мы и имеем у планет Солнечной системы, а у пульсара угол наклона составляет 42,9 градуса, то изменим угол наклона на форме 25 перед расчетом начальных данных на 2,9 градуса. При этом и положение периастра изменим со 178,9 градуса на 185 градусов, чтобы во время работы программы апоастр не проходил через 0 градусов. Вообще то в программе предусмотрено безболезненное прохождение 0 градусов, т.е. когда угол с 360 изменяется до 0, но на всякий случай я решил сместить периастр. Здесь можно было бы вообще не задавать и угол наклона орбиты и положение восходящего узла, т.к. в задаче 2-х тел они все равно не должны изменяться, т.е. можно было бы рассмотреть движение звезд только в расчетной плоскости, но для большей наглядности я решил рассмотреть движение в трех плоскостях.

 

 

Теперь после того, как мы задали параметры орбиты пульсара и начальные данные, переходим к форме 1 программы Solsys7mmm. Сейчас параметры орбиты пульсара и начальные данные находятся в памяти компьютера как текущие значения объектов 0 и 1. Это у нас на форме 1 Солнце и Меркурий. Но нам надо их задать именно как начальные параметры и данные. Если мы нажмем на переключатели <параметры орбит объектов> и <начальные данные>, то мы в соответствующих окошках у Солнца и Меркурия увидим параметры орбиты пульсара или начальные данные. Теперь выбираем переключатель <параметры или данные с экрана> и, нажимая два раза на кнопку <Загрузить параметры на>, загружаем сначала параметры орбиты, а потом начальные данные. Теперь убираем галочки у всех объектов, кроме Солнца и Меркурия, т.е. оставляем в системе только два этих объекта и, нажимая на кнопку <К моделированию>, переходим на форму 2, где мы и будем выполнять вычислительные эксперименты на математической модели системы состоящей из двух звезд, где у нас в модели будут использованы различные законы описывающие гравитационное взаимодействие между объектами.

 

 

При этом определять параметры орбиты пульсара можно не только относительно его компаньона (надо отметить на форме 2 чекбокс <Xsol=0>), как мы это делали для планет, но и относительно центра масс системы (надо отметить чекбокс <Xцм=0>). Все параметры орбиты получатся в обоих случаях те же самые, кроме размера большой полуоси, которая в первом случае получится в два раза больше, чем было задано на форме 25, что, естественно, и должно быть. Потом данные, записанные в различные файлы, подвергнем как обычно статистической обработке на форме 6 и получившиеся при этом значения смещений параметров орбиты пульсара я привожу в таблице 7-3. А подробно о задании начальных данных для двойного пульсара я пишу потому, что эта процедура является не совсем типичной для программы Solsys7mmm поэтому из описания программы может быть не совсем понятно как это сделать для тех, кто решит повторить мои вычислительные эксперименты. Да, при этом надо отметить, что время проведения вычислительных экспериментов здесь у меня не регламентировалось и составляло примерно один или два года, поэтому при проведении повторных вычислительных экспериментов, получающиеся данные, при другой продолжительности вычислительных экспериментов могут незначительно отличаться от приведенных мною.

 

 

Таблица 7-3. Смещения параметров орбиты двойного пульсара PSR 1913+16 рассчитанные при скорости гравитации равной скорости света (Vgr=Vsv) для потенциалов Гербера, Лиенара-Вихерта (Л-В) и для ОТО и при Vgr=100*Vsv для потенциалов Юдина. Все эксперименты выполнены в покоящейся системе двойного пульсара (кроме нескольких при скорости системы VY=100 км/с) и при расчете запаздывающих координат для потенциалов Лиенара-Вихерта и Юдина по методике DX2, кроме нескольких, где было DX1, а также с основным шагом решения P0=3,6 с, который уменьшался в 100 раз вблизи периастра, узла восхождения и при переходе 0 градусов долготы. При этом расчет предварительных ускорений для потенциалов Гербера и Лиенара-Вихерта велся по ускорениям ускорений кроме двух вариантов A=f(V), где расчет велся по скоростям.

 

Варианты расчета

dAlfaP

dAlfaU

dBetta

dEks

dR

ОТО Ландау 

1522189,5+/-16,1

4,6+/-2,0

-0,01+/-0,00

0,0+/-0,0

0,0+/-0,0

ОТО Ландау      VY=100

1522190,5+/-19,7

3,2+/-2,1

-0,00+/-0,00

-0,1+/-0,1

0,0+/-0,0

Гербер1

1522169,2+/-19,2

2,0+/-6,7

-0,00+/-0,00

403,3+/-0,1

1,6+/-0,0

Гербер1              A=f(V)

1522284,1+/-27,4

7,1+/-5,1

-0,00+/-0,00

-2749,6+/-0,2

-10,7+/-0,0

Гербер1              VY=100

1522172,6+/-22,3

2,6+/-3,6

-0,00+/-0,00

402,5+/-0,1

1,6+/-0,0

Гербер2

1522177,8+/-29,4

-11,4+/-7,2

-0,00+/-0,00

403,4+/-0,1

1,6+/-0,0

Л-В полный            DX1

-308784,1+/-127,6

19,6+/-25,3

-0,04+/-0,05

352285,8+/-193,1

1372,9+/-0,5

Л-В полный 

-182326,8+/-79,6

4,7+/-30,4

-0,09+/-0,08

349474,3+/-164,8

1356,7+/-0,6

Л-В полный       A=f(V)

-182471,1+/-87,9

-5,9+/-13,6

-0,02+/-0,04

350485,8+/-243,9

1373,9+/-0,8

Л-В полный       VY=100

-180993,2+/-68,8

8,5+/-26,1

-0,17+/-0,10

293913,9+/-105,6

1138,3+/-0,5

Л-В скалярный 

64134,5+/-60,5

-15,4+/-17,7

0,04+/-0,07

1472,9+/-21,9

7,6+/-0,1

Юдин k1=0 k2=0     DX1

-24,4+/-11,5

2,8+/-7,9

0,04+/-0,01

1771690,8+/-35623,6

31688,4+/-277,6

Юдин k1=0 k2=0

0,4+/-16,9

-15,9+/-14,8

0,04+/-0,01

1787409,1+/-33697,7

31811,7+/-261,0

Юдин k1=0 k2=0  VY=100

21,2+/-6,9

-77,0+/-7,4

-32,3+/-0,8

1691640,5+/-45463,3

31064,9+/-363,6

Юдин k1=1 k2=1

-21,7+/-4,5

-5,2+/-8,3

0,03+/-0,00

4747236,3+/-188749,9

48419,1+/-132,9

 

 

Как видим и здесь ОТО и потенциалы Гербера дают, как и для Меркурия, по смещению периастра результат совпадающий с тем, что получается по формуле (7), т.е. 422,8 градуса за век и это значение совпадает с якобы наблюдаемым значением 423,0 (1522800,0 угл.сек за век). Кстати, это смещение можно рассчитать и по другой формуле (7'), которую получает Ландау [4] в задаче 3 к параграфу 106. При этом мы так же видим, что мои потенциалы, как и для Меркурия, дают большие смещения эксцентриситета и большой полуоси. А изменения положения восходящего узла dAlfaU и угла наклона dBetta во всех вариантах расчета получаются практически в пределах статистической ошибки проведения вычислительных экспериментов (теоретически они все должны быть равны нулю). Но при этом у потенциалов Гербера появились некоторые изменения в смещениях эксцентриситета и большой полуоси, а потенциалы Лиенара-Вихерта дают принципиально другие результаты и здесь я пока не могу объяснить однозначно чем вызваны все эти отличия.

 

dAlfaP= 24*pi^3*a^2 / [c^2*(1-e^2)*T^2]                            (7)

 

dAlfaP= 6*pi*G*(m1+m2) / [c^2*a*(1-e^2)]                         (7')

 

 

Что касается потенциалов Лиенара-Вихерта, то здесь вроде бы все понятно и смещение периастра при методике расчета запаздывания по координатам DX1, естественно, должно отличаться от значения полученного по более точной методике DX2. И, естественно, результаты, получающиеся по упрощенной формуле Лиенара-Вихерта, когда не учитываются ускорения, тоже должны давать результаты отличные от расчета по полной формуле. А вот с потенциалами Гербера мне не понятно почему при расчете предварительных ускорений по скоростям A=f(V) получаются другие значения изменения эксцентриситета и большой полуоси, т.к. смещение периастра, которое, как мы выяснили выше, и вызвано именно ускорением, получается такое же, как и при расчете предварительных ускорений по ускорениям. Да и у потенциалов Лиенара-Вихерта при расчете предварительных ускорений по обоим этим вариантам результаты получаются примерно одинаковые, т.е. оба этих варианта расчета ускорений дают примерно одинаковые результаты. Прояснить этот вопрос могло бы аналитическое решение задачи двух тел с использованием потенциалов Гербера дающее выражения для смещения эксцентриситета и большой полуоси, но, к сожалению, в работе [73] не приводится такого решения для потенциалов Гербера, поэтому придется подождать, когда это сделает кто-то из математиков, заинтересовавшийся этим вопросом.  

 

 

Интересно так же было бы узнать, как официальная наука может трактовать тот факт, что она признает и ОТО и потенциалы Лиенара-Вихерта, но они дают принципиально разные результаты по смещениям параметров орбиты двойного пульсара. Здесь получается примерно такой же парадокс, как и с красным смещением, который я подробно рассмотрел в статье [68] и где получается, что то же самое значение смещения получается и, если мы рассматриваем гравитационное замедление времени, и, если мы рассматриваем покраснение фотона, т.е. потерю им энергии при вылете из гравитационного поля. Но это два совершенно разных физических эффекта и, если они оба имеют место быть, то получается красное смещение в два раза больше, чем мы наблюдаем. Поэтому официальная наука не нашла ничего лучше, чем заявить, что это два разных объяснения одного и того же физического эффекта.

 

 

А в нашем случае наоборот и ОТО и потенциалы Лиенара-Вихерта (при применении их для расчета напряженности гравитационного поля) дают разные результаты, но обе эти теории признаются правильными. Мне могут возразить, что потенциалы Лиенара-Вихерта не применимы к гравитационным полям, но тут возникает вопрос - а почему. Ведь формулы закона тяготения Ньютона для масс и закона Кулона для зарядов полностью идентичны, т.е. наблюдаем полную их эквивалентность. Да и против применения потенциалов Вебера для взаимодействия масс, хотя они были получены для зарядов, у ученых никогда не было никаких возражений. Так почему же тогда нельзя применить и формулу потенциалов Лиенара-Вихерта (без члена магнитной напряженности поля) для взаимодействия масс. А, если их можно применить, то почему они дают результат, который отличается от того, что дает ОТО.

 

 

При этом было бы очень заманчиво использовать наблюдаемые данные по смещениям параметров орбиты двойного пульсара для оптимизации скорости гравитации и скорости всей системы пульсара при использовании моделей с различными формулами описания взаимодействия, т.к. у пульсара на небольших интервалах времени мы можем фиксировать очень большие смещения параметров его орбиты. Но, к сожалению, как это следует из моего исследования (см. приложение 4), так называемые, наблюдаемые значения изменения параметров орбиты двойного пульсара PSR 1913+16, которые получились у Тейлора и Вайсберга, являются очень ненадежными, если не сказать, что просто сфальсифицированными, поэтому я не могу их использовать в своем исследовании по определению скорости гравитации. А, что касается надежности данных по Солнечной системе, то тут, к сожалению, как я выяснил в [1], бесспорно надежными данными можно считать только значение по аномальному смещению перигелия Меркурия, т.е. получается не густо.

 

 

А в заключение могу констатировать, что да, мои потенциалы запаздывающие по координатам и учитывающие динамическое давление при разных скоростях Солнечной системы и разных значениях скорости гравитации дают самые разнообразные значения вековых смещений параметров орбит планет, т.е. есть вероятность найти такие скорости, когда наблюдаемые смещения совпадут с расчетными, но вероятность эта из-за не очень надежных данных по аномальным смещениям всех параметров орбит планет очень мала. Однако, я все же попробую это сделать. При этом у нас для одних вековых смещений оптимальными получатся одни скорости, а для других другие и наша задача будет состоять в том, чтобы найти оптимальные скорости, которые дадут нужные значения вековых смещений для всех параметров и для всех планет. Естественно, делать мы это будем не наугад, а с использованием планов многофакторного планирования, а т.к. при этом у нас должен быть единый критерий оптимизации, то будем использовать комплексный критерий оптимизации (см. формулу (4-2), которую я рассматривал в предыдущей статье [1]), а, что у нас при этом получится, об этом Вы узнаете уже из следующей статьи.

 

 

 

 

Выводы

 

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

 

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

 

3. Формулы Вебера и Гербера для учета запаздывания потенциалов являются не только ошибочными, т.к. их запаздывание потенциалов совершенно не влияет на смещения параметров орбит планет (кроме перигелия), но и являются принципиально антинаучными, т.к. согласно этим формулам ускорение тела (действующая на него сила) зависит от его ускорения (силы), чего не должно быть в принципе. То же самое относится и к формуле Лиенара-Вихерта в которой тоже ускорения зависят от ускорений, а получающееся запаздывание потенциалов не влияет на смещения всех параметров орбит планет, хотя для двойных пульсаров эти потенциалы и дают заметное смещение его параметров орбиты.

 

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

 

5. Наблюдаемое значение уменьшения периода обращения двойного пульсара PSR 1913+16, которое приводят Тейлор и Вайсберг является не столько наблюдаемым значением сколько расчетным и  полученным при этом с использованием уравнений ОТО. Таким образом, и оценка скорости гравитации Ван Фландерном по изменению этого параметра является некорректной, а, если учесть, что в расчетах у него имеется и явная ошибка, то следует заключить, что на сегодняшний день не существует объективных оценок скорости гравитации.

 

 

 

P.S. Я, конечно же, много критиковал в этой статье математиков за то, что они из-за отсутствия математических способностей переквалифицировались в математико-физиков и может быть, где-то я был и не прав, но и Вы поймите меня. Я вот при написании предыдущей статьи «Опять о принципе наименьшего действия» столкнулся с проблемой, которую математики не могут решить уже лет 300, т.е. разрешить спор двух создателей вариационного исчисления Эйлера и Лагранжа о том варьируется там время или нет. И тут получается как в анекдоте. Я их спрашиваю - кто прав, а они отвечают и Эйлер прав и Лагранж прав. Я говорю, что такого не может быть, а они отвечают, что и Вы правы. Вот и все на что способны современные математики, которые в основной массе являются математико-физиками (не путать с физико-математиками, коих считанные единицы). Я, конечно, имею, как механик, свое мнение по этому вопросу и считаю, что время вообще не может присутствовать в вариационном исчислении, но хотелось бы увидеть чисто математические выкладки по этому вопросу.

 

 

А сейчас вот после написания этой статьи занялся непосредственно поиском скорости гравитации с использованием планов многофакторного планирования и опять у меня возникают сложности из-за математиков, т.к. все эти планы показывают отличные результаты, когда критерием оптимизации является сам отклик системы. А мне надо оптимизировать параметры по отклонению (причем по абсолютной величине) отклика системы от оптимального значения. Я уже писал об этой проблеме в статье «Аномальные смещения параметров орбит планет», где на конкретных тестах показал, что по этому критерию оптимизации все планы (ортогональные, рототабельные и т.д.) показывают только удовлетворительные результаты. Но математики упорно лезут в физику, где они сто лет не нужны, а вот заняться своими математическими проблемами никак не хотят. А я был бы им признателен, если бы они разработали планы многофакторного планирования, где можно производить оптимизацию параметров систем по критерию отклонения отклика системы от оптимального значения (как по абсолютной величине, так и этого отклонения в квадрате). И, как бы не была печальна ситуация в современной науке, я все же надеюсь, что, после этого моего крика вопиющего в пустыне, хоть один математик займется реальной математикой.

 

 

 

                                                                 Список литературы

 

1.- С.Ю.Юдин //Аномальные смещения параметров орбит планет//  2012, - 32 с http://modsys.narod.ru/Stat/Stat_Est/Anomal/Anomal1.html

2.- П. Гербер //Пространственное и временное распространение гравитации// http://bourabai.kz/articles/gerber/gerber-rus.htm

3. – Карим Хайдаров //Эфир великий часовщик// http://www.bourabai.kz/watchmaker.htm

4.- Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: Учеб. пособие в 10-ти т.Т.2. Теория поля. - 7-е изд. М.: Физматлит, 1988. - 512 с.

5.- Р. Фейнман, Р. Лейтон, М. Сэндс Фейнмановские лекции по физике. Том. 6. Электродинамика - 339 с.

6. – S. Newcomb //Four inner planets and the fundamental constans of astronomy// - Washington, Government printing office, 1895, 202 p.

7. – М.Ф.Субботин //Введение в теоретическую астрономию// - М.: Наука, 1968, - 800 с.

8.- С.Ю.Юдин //Опять о принципе наименьшего действия//  2010, - 43 с http://modsys.narod.ru/Stat/Stat_Est/Princip2/princip21.html

9. – Н.Т.Роузвер //Перигелий Меркурия. От Леверье до Эйнштейна//: пер. с англ. – М.: Мир, 1985 - 246 с.

10. – О.В.Зайцев //Принцип Маха и орбитальная прецессия планет //    http://www.physical-congress.spb.ru/russian/pmoop/pmoop.asp#_Toc479550256

11. И.И.Смульский /Теория взаимодействия/ - Новосибирск: Из-во Новосиб. Ун-та, НИЦ ОИГГМ СО РАН, 1999, 294 с. http://www.ikz.ru/~smulski/TVfulA5_2.pdf

12. Сухоруков Г.И., Сухоруков В.И., Сухоруков Э.Г, Сухоруков Р.Г. //Реальный физический мир без парадоксов// Братск. БрГТУ, 2001 – 229 с.    http://suxorucov.narod.ru/Kniga.rar

13. – М.Корнева, В.Кулигин, Г.Кулигина //Вы очень жаждете иметь новый Чернобыль?// http://www.sciteclibrary.ru/rus/catalog/pages/9448.html

14. – В.В.Ерохин //Основы конструктивной электродинамики. Часть 1. Магнитное поле в релятивистском приближении// - Торез: , 2002. – 92 с. http://ru.convdocs.org/download/docs-107777/107777.doc

15. – В.И.Коробейников //Мифы и реальность EH антенн// 2005   http://www.qrz.ru/articles/detail.phtml?id=282

16. – А.Г. Замятин //Принцип близкодействия// Свердловск, 1988, 153 с. http://yadi.sk/d/mpuZzuId1xkL5

17. – Н.В.Купряев //Классический эффект Доплера// http://www.sciteclibrary.ru/rus/catalog/pages/8803.html

18. – О.Е.Акимов //О формуле, описывающей классический эффект Доплера// http://sceptic-ratio.narod.ru/fi.htm#Doppler

19. - Г. Герц //Принципы механики, изложенные в новой связи// М.: Академия наук СССР, 1959, 386 с.

20. – В.А. Красильников, В.В. Крылов. Введение в физическую акустику. М.: Наука, 1984. 403 с.

21. - Лагранж //Аналитическая механика// Т.1 М-Л.: ОНТИ-НКТП, 1938, 348 с.  

Лагранж //Аналитическая механика// Т.2 М-Л.: ОНТИ-НКТП, 1950, 440 с.  

22. – А.Т. Григорьян, Л.С. Полак //Основные идеи механики Генриха Герца// статья размещена в книге Г.Герц  //Принципы механики, изложенные в новой связи// (см п.19).

23. - Л.Эйлер //Основы динамики точки// М-Л.: ОНТИ-НКТП, 1938, 500 с.   http://narod.ru/disk/3486349000.68bf093210ba8df860be54c7ad004003/Euler_1938_4.pdf.html

24. - Брайс Ф. Девитт. Динамическая теория групп и полей. Москва, Наука, 1987 г. 

 http://bookfi.org/dl/453421/8fcfca

25. - Юдин С.Ю. //Две меры механической формы движения материи//. 2005, 16 с.      http://modsys.narod.ru/Stat/Stat_Est/Udar/udar.html

26. - Юдин С.Ю. //Методы оценки эффективности применения МТА и оптимизации их параметров с использованием математических моделей/ диссертация на соискание ученой степени кандидата технических наук, Волгоград.: 1991, 270 с.

27. - Дж.К. Максвелл. Трактат об электричестве и магнетизме. В двух томах. Т.II. М.: Наука, 1989, 436 с.

28.- Юдин С.Ю. //О принципах кратчайшего времени и наименьшего действия// 2005, 20 с.  http://modsys.narod.ru/Stat/Stat_Est/Princip/princip.html

29. - Полак Л.С. Вариационные принципы механики // Вариационные принципы механики// Под ред. Л.С. Полака, М.: Физматгиз, 1959. с. 780-879. http://ftp.kinetics.nsc.ru/chichinin/pmlic.htm#P

30. – Ньютон Исаак. Математические начала натуральной философии. - М.: Наука, 1989, 688 с.

31. - Юдин С.Ю. //О равноденствиях Гиппарха и Птолемея// 2010, 55 с.  http://modsys.narod.ru/Stat/Stat_Est/Hipp_Ptolo/Hipp_Ptolo2.html

32. Юдин С.Ю. //О формуле Планка и кванте действия// 2005, 26 с. http://modsys.narod.ru/Stat/Stat_Est/Plank/plank.html 

33. Юдин С.Ю. //Выбор языка программирования для научных работников// 2007, 24 с. http://modsys.narod.ru/Stat/Stat_Prog/Vibor/Vibor2.html 

34. В.С.Вепринцев //Об алгоритмах преобразований Галилея и Лоренца и динамике// http://vvs-ya.narod.ru/about12.htm

35. Л.Ландау, Л.Пятигорский, Механика. – М-Л.: Техтеорлит, 1940. - 200 с.

36. – А.Эйнштейн //Объяснение движения перигелия Меркурия в общей теории относительности// Собрание научных трудов. Том 1. М.: Наука, 1965, 701 с.

37.- E. Myles Standish and James G. Williams. Orbital Ephemerides of the Sun, Moon, and Planets http://iau-comm4.jpl.nasa.gov/XSChap8.pdf 

38. - Е.Майлс Стэндиш, Джеймс Г.Вильямс. Орбитальные эфемериды Солнца, Луны и планет  http://vadimchazov.narod.ru/text_htm/xsru00.htm  )

39.- Шмидт С.Н. Способ определения скорости «корабля Галилея» и ускорения «лифта Эйнштейна»

http://www.dot2008-shmidt.narod.ru/Exhibits/Exhibit9.htm

40.- Томилин А.К. //Основы обобщенной электродинамики// 2009, 126 с.  http://vev50.narod.ru/Tomilin.html

41.- П.С. Лаплас Изложение системы мира Л.: Наука, 1982, 376 с.

42.- T. Van Flandern The Speed of Gravity What the Experiments Say http://www.metaresearch.org/cosmology/speed_of_gravity.asp

43.- А. Лайтман, В. Пресс, Р. Прайс, С. Тюкольски Сборник задач по теории относительности и гравитации. Пер. с англ. А.П.Бондарева и Ю.А.Данилова. М.: Мир, 1979, 536 с.

44.- Горелик И. Устойчивость солнечной системы http://darkenergy.narod.ru/solarru.html

45.- Юдин С.Ю. Кинематическая теория планет.  2012, 72 с http://modsys.narod.ru/Stat/Stat_Est/Kinematik/Kinematik1.zip

46.- Дж.Х. Тэйлор (мл.) Двойные пульсары и релятивистская гравитация (Нобелевская лекция. Стокгольм, 8 декабря 1993 г.) Успехи физических наук, Июль 1994 г. Том 164, №7

47.- J.H. Taylor Binary Pulsars and Relativistic Gravity, Nobel Lecture, December 8, 1993. http://www.nobelprize.org/nobel_prizes/physics/laureates/1993/taylor-lecture.pdf

48.- Липунов В.М. В мире двойных звезд - М.: Наука. Гл. ред. физ.-мат. лит., 1986 - 208 с

49.- Засов А.В., Постнов К.А. Общая астрофизика. - Фрязино, 2006. - 496 с.  http://alexandr4784.narod.ru/zasow.htm

50.- R.A. Hulse The Discovery of the Binary Pulsar, Nobel Lecture, December 8, 1993.     http://www.nobelprize.org/nobel_prizes/physics/laureates/1993/hulse-lecture.pdf

51.- Родин А.Е. Прецизионная астрометрия пульсаров в присутствии низкочастотных шумов. Диссертация на соискание ученой степени кандидата физико-математических наук. М. 2000 г.

52.- Weisberg, J.M., Taylor, J.H. and L.A. Fowler Gravitational Waves from an Orbiting Pulsar, Scientific American 1981 http://www.nature.com/scientificamerican/journal/v245/n4/pdf/scientificamerican1081-74.pdf

53.- R.A. Hulse and J.H. Taylor Discovery of a Pulsar in Binary System The Astrophysical Journal, 195. 15 January 1975. p. 51-53.

54.- Бисноватый-Коган Г.С. Двойные и подкрученные радиопульсары: через 30 лет после наблюдательного открытия УФН 176 59–75 c. (2006)  http://ufn.ru/ufn06/ufn06_1/Russian/r061c.pdf

55.- I. H. Stairs, S. E. Thorsett, J. H. Taylor and A. Wolszczan 

Studies of the Relativistic Binary Pulsar PSR B1534+12. I. Timing Analysis The Astrophysical Journal 581, 10 December 2002, p. 501-508.

56.- I.H. Stairs Testing General Relativity with Pulsar Timing   http://www.livingreviews.org/lrr-2003-5

57.- M. Maggiore Gravitation Waves: Volume 1: Theory and Experiments, Chapter 6   http://www.oxfordscholarship.com/view/10.1093/acprof:oso/9780198570745.001.0001/acprof-9780198570745   

58.- Д. тер Хаар Пульсары. Успехи физических наук, Июль 1976 г. Том 119, вып. 3.

59.- J.H. Taylor, J.M. Weisberg A New Test of General Relativity: Gravitational Radiation and the Binary Pulsar PSR 1913+16, The Astrophysical Journal, 253. 15 February 1982, p. 908-920.

60.- J.M. Weisberg, J.H. Taylor The Binary Pulsar B1913+16, meaning in China of August 2002 on Radio Pulsars Conference.

61.- Бутяев Е.А. Загадки природы Гл. 4 О противоречивости различных интерпретаций красного смещения http://butjaev.narod.ru/BIB/gl4.pdf

62.- T. Damour and J.H. Taylor On the Orbital Period Change of the Binary Pulsar PSR 1913+16, The Astrophysical Journal, 366. 10 January 1991, p. 501-511.

63.- S.A. Balbus and K. Brecher Tidal Friction in the Binary Pulsar System PSR 1913+16, The Astrophysical Journal, 203. 1 January 1976, p. 202-205.

64.- M. Bails and other Transformation of a Star into a Planet in a Millisecond Pulsar Binary  http://arxiv.org/PS_cache/arxiv/pdf/1108/1108.5201v1.pdf

65.- Open Exoplanet Catalogue  http://www.openexoplanetcatalogue.com/system.html?id=PSR%20J1719-1438%20b              

66.- В. Воробьева Алмазная планета у нейтронной звезды  http://allplanets.ru/novosti_2011.htm

67.- Юдин С.Ю. // Математическое описание явлений Природы// 2014, 70 с.     http://modsys.narod.ru/Stat/Stat_Est/Matopisanie1.html

68.- Юдин С.Ю. //Эффект Доплера// 2014, 20 с.  http://modsys.narod.ru/Stat/Stat_Est/Dopler2.html

69.- Салль С.А. Истоки и заблуждения релятивизма http://www.spbs.rusphysics.ru/files/Istoki.pdf

70.- A. Lienard Produit par une electrique concentree en un point et anemee d'un mouveement quelconque. L'Eclairage Electrique. Tome XVI. 1898.

71.- E. Wiechert Elektrodinamishe Elementargesetze. 1901.

72.- С.В. Беллюстин Классическая электронная теория М.: Высшая школа, 1971, 352 с.

73.- А.Ф. Богородский Всемирное тяготение. Киев: Наукова думка, 1971, 351 с.

74.- Laplace, P., Mechanique Celeste. English translation. Volume IV. Boston, 1839, 1018 p.

75.- А.Ф. Богородский Уравнения поля Эйнштейна и их применение в астрономии. Киев: Киевский университет, 1962, 196 с.

76.- В. Орлов. Теория гравитации Пауля Гербера https://sites.google.com/site/testsofphysicaltheories/russian/gerber  

 

 

 

P.S. Предыдущую 2-ю версию этой статьи можно скачать отсюда http://modsys.narod.ru/Stat/Stat_Est/Vlijanie/Vlijanie2doc.zip или отсюда https://googledrive.com/host/0BwnV2Ac6zvhMalpOVGktQ1Jic1U/Stat/Stat_Est/Vlijanie/Vlijanie2doc.zip .

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

                                                                                                                                                 ПРИЛОЖЕНИЕ 1

 

 

Влияние скорости гравитации и скорости Солнечной системы на вековые смещения параметров орбит планет при учете запаздывания потенциалов по координатам при разных значениях скорости гравитации и скорости Солнечной системы.

 

 

 

 

 

 

Вычислительные эксперименты проводились на математической модели Солнечной системы, где учитывалось запаздывания по координатам. В начале эксперимента задавались начальные параметры орбит планет и их долгота в плоскости фиксированной эклиптики, т.е. в эпохе J2000, на 12 часов 1 января 1801 года, по моей теории планет Ser0. Затем по стандартным формулам определялись координаты планет и их скорости при условии, что эти параметры орбит получены при движении планет относительно неподвижного Солнца, находящегося в начале абсолютно неподвижной системы координат. А потом начальные координаты и скорости Солнца с учетом воздействия планет не уточнялись, хотя вообще-то желательно бы было это сделать (программа позволяет). После этого я задавал произвольную скорость всей Солнечной системе (от -300 до +300 км/с) по одной из осей координат (по остальным осям скорость принималась равной нулю), т.е. увеличивал на эту величину начальные скорости всех планет и Солнца и начинал вычислительный эксперимент, т.е. моделировал движение объектов Солнечной системы.

 

 

Во время эксперимента я определял до 31 декабря 2000 года, при каждом обороте планеты вокруг Солнца, значения получающихся значений долготы перигелия AlfaP и восходящего узла AlfaU, угла наклона орбиты Betta, эксцентриситета Eks, большой полуоси Rsr, угловой скорости Wsr и записывал их в файл. Затем, полученный массив данных, я подвергал статистической обработке и находил вековые смещения параметров орбит планет. Подробности о примененной мною методике статистической обработке этих данных смотрите в предыдущей статье [1].

 

 

При этом решение системы дифференциальных уравнений выполнялось с основным шагом 3600 секунд, а вблизи перигелия или узла восхождения с уменьшенным шагом решения. Для Меркурия этот уменьшенный шаг составлял 3,6 сек, для Венеры 7,2 сек, для Земли 18 сек и для Марса 36 сек. А область, где уменьшенный шаг решения применялся, определялась в интервале +/-1 градус от полученного значения перигелия или узла восхождения при предыдущем обороте планеты. При этом, т.к. у Венеры наблюдается очень большой разброс в мгновенных значениях долготы перигелия, вызванный очень маленьким эксцентриситетом орбиты, когда воздействия от других планет его могут значительно сместить, я для нее определял перигелий и узел восхождения в увеличенном интервале, т.е. +/-2 градуса от предыдущих значений.

 

 

Скорость гравитации в экспериментах изменялась от одной скорости света (10^0) до 1000 скоростей света (10^3), а скорость Солнечной системы по осям координат от -300 до +300 км/с. При этом в таблицах в первой колонке указаны показатели степени у 10 для задания скорости гравитации (10^0…10^3), а звездочками отмечены вековые смещения, которые рассчитаны только приблизительно из-за нелинейностей графиков смещений. А при очень больших нелинейностях, вызванных малой скоростью гравитации и уменьшением точности расчета одной из координат при ее большом увеличении от постоянной скорости всей системы, например, когда 100 лет смещение уменьшается, а потом 100 лет увеличивается или изменяется очень не устойчиво, я вообще не привожу полученные значения.

 

 

 

 

 

 

Таблица 1-1-X. AlfaP1-VXsys  - вековые смещения перигелия Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

529,23

1787,01

-

-

1

307,47

369,66

410,41

529,63

624,74

721,49

832,64

2

502,16

511,37

520,65

529,68

539,10

548,90

557,46

3

526,91

528,10

528,82

529,69

530,58

531,49

532,40

 

Таблица 1-1-Y. AlfaP1-VYsys  - вековые смещения перигелия Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-15612,46*

-9466,43

-4237,82

529,23

4990,88

9117,31

13154,15

1

-859,33

-392,79

+69,79

529,63

986,07

1438,82

1966,32

2

+392,20

+438,03

483,90

529,68

575,67

621,36

667,33

3

+515,89

+520,55

525,11

529,69

534,26

538,83

543,48

 

Таблица 1-1-Z. AlfaP1-VZsys  - вековые смещения перигелия Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-1706,05*

-1033,93

-296,99

529,23

1435,84

2235,04*

3089,96*

1

+275,92

+360,20

+444,56

529,63

615,33

701,34

786,61

2

504,12

512,59

521,16

529,68

538,32

547,04

555,68

3

527,10

527,93

528,80

529,69

530,54

531,43

532,29

 

Таблица 1-2-X. AlfaU1-VXsys  - вековые смещения узла восхождения Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-16640,93*

-11741,57

-6560,55

-451,47

8293,04

-

-

1

-2464,39

-1811,98

-1152,11

-451,39

+262,54

+1002,29

+1768,87

2

-660,88

-591,69

-521,57

-451,40

-381,05

-310,46

-239,70

3

-472,50

-465,48

-458,42

-451,40

-444,32

-437,33

-430,28

 

Таблица 1-2-Y. AlfaU1-VYsys  - вековые смещения узла восхождения Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

14947,73

10519,26

5352,15

-451,47

-6886,54

-13952,81

-21722,00

1

+1355,45

+759,52

+157,22

-451,39

-1066,22

-1687,08

-2315,12

2

-268,13

-329,17

-390,16

-451,40

-512,65

-573,93

-635,45

3

-433,04

-439,14

-445,25

-451,40

-457,52

-463,62

-469,74

 

Таблица 1-2-Z. AlfaU1-VZsys  - вековые смещения узла восхождения Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

-451,47

-

-

-

1

-31721,79*

-18967,90*

-8732,18

-451,39

6304,91

11867,56

16477,37*

2

-2757,47

-1972,93

-1204,59

-451,40

+286,60

+1010,06

+1719,17

3

-675,95

-601,20

-526,17

-451,40

-376,96

-302,65

-228,60

Таблица 1-3-X. Betta1-VXsys  - вековые смещения угла наклона Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-2898,74*

-2078,27*

-1254,58

-21,47

2183,42*

-

-

1

-450,35

-313,94

-171,04

-21,46

+136,51

303,73

481,17

2

-67,24

-52,04

-36,79

-21,45

-6,05

+9,44

+25,03

3

-26,07

-24,53

-23,00

-21,45

-19,91

-18,38

-16,84

 

Таблица 1-3-Y. Betta1-VYsys  - вековые смещения угла наклона Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

5542,41*

3238,99*

1452,55

-21,47

-1278,81

-2368,88

-3321,74*

1

394,83

253,82

+115,07

-21,46

-155,87

-287,80

-419,94

2

+12,27

+5,67

-7,89

-21,45

-34,99

-48,51

-61,99

3

-17,39

-18,74

-20,10

-21,45

-22,81

-24,17

-25,52

 

Таблица 1-3-Z. Betta1-VZsys  - вековые смещения угла наклона Меркурия (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

-21,47

17155,05

33360,16

47529,64*

1

-4730,89

-3225,31

-1657,26

-21,46

1654,66

3357,96

5079,10

2

-517,17

-352,36

-187,12

-21,45

+144,58

311,02

477,82

3

-71,19

-54,63

-38,04

-21,45

-4,87

11,71

28,31

 

Таблица 1-4-X. Eks1-VXsys  - вековые смещения эксцентриситета Меркурия (увеличено в 1000000 раз)

 

n

-300

-200

-100

0

100

200

300

0

10174,14*

7302,69*

3556,98*

20,91

-5677,12

-

-

1

1320,84

928,96

388,11

20,53

-445,67

-924,33

-1430,77

2

186,87

111,61

66,13

20,49

-25,28

-71,23

-117,29

3

34,15

29,69

25,09

20,49

+15,92

+11,31

+6,71

 

Таблица 1-4-Y. Eks1-VYsys  - вековые смещения эксцентриситета Меркурия (увеличено в 1000000 раз)

 

n

-300

-200

-100

0

100

200

300

0

-4470,99*

-2503,59*

-1117,04

20,91

928,55

1673,37*

-

1

-295,24

-187,60

-82,48

20,53

121,37

220,64

334,29

2

-10,03

+0,04

+10,35

20,49

30,72

40,86

50,94

3

+17,46

18,42

19,48

20,49

21,53

22,54

23,46

 

Таблица 1-4-Z. Eks1-VZsys  - вековые смещения эксцентриситета Меркурия (увеличено в 1000000 раз)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

20,91

874,35*

2086,37*

-

1

-107,35*

-71,15

-28,19

20,53

75,89

137,22

204,91*

2

+5,08

+10,16

+15,31

20,49

25,80

31,18

36,58

3

18,93

19,45

19,97

20,49

21,01

21,15

22,04

Таблица 1-5-X. Rsr1-VXsys  - вековые смещения большой полуоси Меркурия (тыс. км)

 

n

-300

-200

-100

0

100

200

300

0

28359,33

19419,66

10057,49

1,16

-11204,99

-25301,75*

-38035,28*

1

3104,44

2079,32

1045,98

0,12

-1055,43

-2122,55

-3201,85

2

314,53

209,80

104,96

0,01

-105,05

-210,21

-315,49

3

31,50

21,00

10,50

0,00

-10,50

-21,00

-31,49

 

Таблица 1-5-Y. Rsr1-VYsys  - вековые смещения большой полуоси Меркурия (тыс. км)

 

n

-300

-200

-100

0

100

200

300

0

-9287,97*

-5579,17*

-2545,99

1,16

2159,49

3995,73

5520,32*

1

-720,01

-476,02

-235,99

0,12

232,38

460,89

685,30

2

-70,42

-46,90

-23,42

0,01

23,40

46,78

70,10

3

-7,02

-4,68

-2,34

0,00

2,34

4,70

7,02

 

Таблица 1-5-Z. Rsr1-VZsys  - вековые смещения большой полуоси Меркурия (тыс. км)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

1,16

1853,23*

4869,78*

8839,65*

1

-280,14

-201,74

-108,27

0,12

123,24

260,89

412,91

2

-34,06

-22,84

-11,49

0,01

11,67

23,46

35,42

3

-3,46

-2,31

-1,16

0,00

1,16

2,32

3,48

 

Таблица 1-6-X. Wsr1-VXsys  - вековые смещения орбитальной угловой скорости Меркурия (рад/век)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-515,7168*

-0,0781

1221,9867*

-

-

1

-190,6637

-131,6263

-68,3515

-0,0077

73,8607

154,1369

241,8354

2

-21,0359

-14,0791

-7,0673

-0,0007

7,1224

14,3009

21,5379

3

-2,1259

-1,4176

-0,7094

-0.0000

0,7098

1,4197

2,1304

 

Таблица 1-6-Y. Wsr1-VYsys  - вековые смещения орбитальной угловой скорости Меркурия (рад/век)

 

n

-300

-200

-100

0

100

200

300

0

-

459,0944*

187,3206

-0,0781

-136,3162

-238,5178*

-314,5845*

1

49,8034

32,6638

16,0672

-0,0077

-15,5838

-30,6835

-45,2992

2

4,7691

3,1739

1,5840

-0,0007

-1,5799

-3,1557

-4,7260

3

0,4741

0,3166

0,1579

-0.0000

-0,1581

-0,3174

-0,4740

 

Таблица 1-6-Z. Wsr1-VZsys  - вековые смещения орбитальной угловой скорости Меркурия (рад/век)

 

n

-300

-200

-100

0

100

200

300

0

-

-

-

-0,0781

-118,6362

-289,4285*

-501,5755*

1

19,1215

13,7283

7,3426

-0,0077

-8,2932

-17,4785

-27,5325

2

2,3039

1,5451

0,7769

-0,0007

-0,7877

-1,5838

-2,3895

3

0,2341

0,1563

0,0782

-0.0000

-0,0783

-0,1566

-0,2351

Таблица 2-1-X. AlfaP2-VXsys  - вековые смещения перигелия Венеры (угл. секунды)

 

n

-300

-200

-100

0

100

200

300

0

-

978,46*

226,64

58,64

72,00

-237,01

-338,66

1

67,58

97,64

-52,20

57,75

38,13

54,06

23,02

2

60,36

58,00

58,23