Обзор работ по определению скорости гравитации

[Petrovich_Tot]

Проверь модель вывода результатов и опубликуй здесь

Код: [Выделить]

Reap[NDSolve[{Table[eq1[[k]] == 0, {k, 1, 3}],
    Table[eq2[[k]] == 0, {k, 1, 3}], eq0}, {x1[t], y1[t], z1[t],
    x2[t], y2[t], z2[t], x1'[t], y1'[t], x2'[t], y2'[t]},
  {t, 0, 10*t0},
   Method -> {"EventLocator", "Event" -> y1[t], "Direction" -> 1,
    "EventAction" :> Sow[{t, x1[t], y1[t], x2[t], y2[t]}]}]]

[BJIaquMup]

Проверь модель вывода результатов и опубликуй здесь

Код: [Выделить]

Reap[NDSolve[{Table[eq1[[k]] == 0, {k, 1, 3}],
    Table[eq2[[k]] == 0, {k, 1, 3}], eq0}, {x1[t], y1[t], z1[t],
    x2[t], y2[t], z2[t], x1'[t], y1'[t], x2'[t], y2'[t]},
  {t, 0, 10*t0},
   Method -> {"EventLocator", "Event" -> y1[t], "Direction" -> 1,
    "EventAction" :> Sow[{t, x1[t], y1[t], x2[t], y2[t]}]}]]

Чё понаписал? Чё проверять. Ничего это не работает. Это откуда, куда и к чему?
Пиши нормальное описание и вывод результата.

[Petrovich_Tot]

Чё понаписал? Чё проверять. Ничего это не работает. Это откуда, куда и к чему?
Пиши нормальное описание и вывод результата.

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

Код: [Выделить]

c0 = 299792458.; G = 6.67384*10^(-11); GM =
 G*1.9891*10^(30); a = 1950100000; T1 = 7.7519387738647*3600; t0 =
 T1*Sqrt[GM/a^3]; e0 = 0.6171334; m2 = 1.387; m1 = 1.441; v0 =
 450000/Sqrt[GM/a]*2.003857; c =
 c0/Sqrt[GM/a]; x10 = (1 - e0)*m2/(m1 + m2); x20 = -(1 - e0)*
  m1/(m1 + m2); y10 = 0.; y20 = 0; v1y =
 v0*m2/(m1 + m2); v2y = -v0*m1/(m1 + m2); v1x = 0.; v2x = 0;
r1 = {x1[t], y1[t], z1[t]};
r2 = {x2[t], y2[t], z2[t]};
v1 = D[r1, t];
v2 = D[r2, t];
r = r1 - r2;
ra = Norm[r];
rn = r/ra;
eq1 = -D[v1, t] -
   m2*rn/ra^2 + (m2/(c*ra)^2)*(rn*((5*m1 + 4*m2)/ra - v1.v1 -
         2*v2.v2 + 4*v1.v2 + (3/2)*(v2.rn)^2) + (v1 -
         v2)*(rn.(4*v1 - 3*v2)));
eq2 = -D[v2, t] +
   m1*rn/ra^2 - (m1/(c*ra)^2)*(rn*((5*m2 + 4*m1)/ra - v2.v2 -
         2*v1.v1 + 4*v1.v2 + (3/2)*(v1.rn)^2) - (v1 -
         v2)*(rn.(4*v2 - 3*v1)));
eq0 = {x1[0] == x10, y1[0] == 0, z1[0] == 0, x1'[0] == 0,
   y1'[0] == v1y, z1'[0] == 0, x2[0] == x20, y2[0] == 0, z2[0] == 0,
   x2'[0] == 0, y2'[0] == v2y, z2'[0] == 0};
Reap[NDSolve[{Table[eq1[[k]] == 0, {k, 1, 3}],
    Table[eq2[[k]] == 0, {k, 1, 3}], eq0}, {x1[t], y1[t], z1[t],
    x2[t], y2[t], z2[t], x1'[t], y1'[t], x2'[t], y2'[t]},
  {t, 0, 10*t0},
   Method -> {"EventLocator", "Event" -> y1[t], "Direction" -> 1,
    "EventAction" :> Sow[{t, x1[t], y1[t], x2[t], y2[t]}]}]]


[BJIaquMup]

Перепиши от начала до конца, что должно быть.

[Petrovich_Tot]

что переписать?

[BJIaquMup]

Всю программу. От и до.

[Petrovich_Tot]

для чего? у меня все работает, код опубликовал.

[BJIaquMup]

для чего? у меня все работает, код опубликовал.

А у меня не работает. Потому что возможно у меня в математике №5 нет искуственново интеллекту. Или у меня нету интелекту.
Вообще-то желательно в тексте программы писать комментарии.
Вот это я не понял, чё к чему добавлять?

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

[Petrovich_Tot]

А у меня не работает. Потому что возможно у меня в математике №5 нет искуственново интеллекту. Или у меня нету интелекту.
Вообще-то желательно в тексте программы писать комментарии.
Вот это я не понял, чё к чему добавлять?

не надо добавлять, используй этот код

Код: [Выделить]

c0 = 299792458.; G = 6.67384*10^(-11); GM =
 G*1.9891*10^(30); a = 1950100000; T1 = 7.7519387738647*3600; t0 =
 T1*Sqrt[GM/a^3]; e0 = 0.6171334; m2 = 1.387; m1 = 1.441; v0 =
 450000/Sqrt[GM/a]*2.003857; c =
 c0/Sqrt[GM/a]; x10 = (1 - e0)*m2/(m1 + m2); x20 = -(1 - e0)*
  m1/(m1 + m2); y10 = 0.; y20 = 0; v1y =
 v0*m2/(m1 + m2); v2y = -v0*m1/(m1 + m2); v1x = 0.; v2x = 0;
r1 = {x1[t], y1[t], z1[t]};
r2 = {x2[t], y2[t], z2[t]};
v1 = D[r1, t];
v2 = D[r2, t];
r = r1 - r2;
ra = Norm[r];
rn = r/ra;
eq1 = -D[v1, t] -
   m2*rn/ra^2 + (m2/(c*ra)^2)*(rn*((5*m1 + 4*m2)/ra - v1.v1 -
         2*v2.v2 + 4*v1.v2 + (3/2)*(v2.rn)^2) + (v1 -
         v2)*(rn.(4*v1 - 3*v2)));
eq2 = -D[v2, t] +
   m1*rn/ra^2 - (m1/(c*ra)^2)*(rn*((5*m2 + 4*m1)/ra - v2.v2 -
         2*v1.v1 + 4*v1.v2 + (3/2)*(v1.rn)^2) - (v1 -
         v2)*(rn.(4*v2 - 3*v1)));
eq0 = {x1[0] == x10, y1[0] == 0, z1[0] == 0, x1'[0] == 0,
   y1'[0] == v1y, z1'[0] == 0, x2[0] == x20, y2[0] == 0, z2[0] == 0,
   x2'[0] == 0, y2'[0] == v2y, z2'[0] == 0};
Reap[NDSolve[{Table[eq1[[k]] == 0, {k, 1, 3}],
    Table[eq2[[k]] == 0, {k, 1, 3}], eq0}, {x1[t], y1[t], z1[t],
    x2[t], y2[t], z2[t], x1'[t], y1'[t], x2'[t], y2'[t]},
  {t, 0, 10*t0},
   Method -> {"EventLocator", "Event" -> y1[t], "Direction" -> 1,
    "EventAction" :> Sow[{t, x1[t], y1[t], x2[t], y2[t]}]}]]


[BJIaquMup]

NDSolve::bdmtd: The value of the option Method->\[LeftSkeleton]1\
\[RightSkeleton] is not a known built-in method, a symbol which could be a \
user-defined method, or a list with a name followed by method options

Это результат исполнения.

[Petrovich_Tot]

NDSolve::bdmtd: The value of the option Method->\[LeftSkeleton]1\
\[RightSkeleton] is not a known built-in method, a symbol which could be a \
user-defined method, or a list with a name followed by method options

Это результат исполнения.

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

[BJIaquMup]

Да ладно, пиши как есть, я на слово поверю. В принципе, у меня есть чем заняться. 
Просто, хотел вникнуть.

[Petrovich_Tot]

так вникай - моделируем смещение орбиты двойного пульсара с учетом релятивистских поправок. В ОТО задача в общем виде решена в 1938 (Robertson), ответ

Код: [Выделить]

teor = 6*Pi*GM*(m1 + m2)/(a*(1 - e0^2)*c0^2)
0.0000652111
Можешь проверить, все параметры заданы в программе.
Считаем орбиты с опцией локатор событий

Код: [Выделить]

Reap[NDSolve[{Table[eq1[[k]] == 0, {k, 1, 3}],
   Table[eq2[[k]] == 0, {k, 1, 3}], eq0}, {x1[t], y1[t], z1[t], x2[t],
    y2[t], z2[t], x1'[t], y1'[t], x2'[t], y2'[t]}, {t, 0, 10*t0},
  Method -> {"EventLocator",
    "Event" -> (x1'[t] - x2'[t])*(x1[t] - x2[t]) + (y1'[t] -
         y2'[t])*(y1[t] - y2[t]), "Direction" -> 1,
    "EventAction" :>
     Sow[{t/t0, ArcTan[ x1[t] - x2[t], y1[t] - y2[t]]/teor}]}]]Получаем ответ в формате {время/t0, угол/teor}

Код: [Выделить]

{{{1.00091, 1.0109}, {2.00182, 2.01899}, {3.00273, 3.02878}, {4.00363,
     4.03746}, {5.00454, 5.04621}, {6.00544, 6.05499}, {7.00634,
    7.06502}, {8.00724, 8.07369}, {9.00814, 9.08244}}}}Т.е. период чуть больше измеренного в эксперименте и угол смещения периастра чуть больше теоретического. 

[BJIaquMup]

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

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

Давайте, я отвечу в вашей теме, ибо, по теме.

Ну, если вы заметили, что данные (полученные, стесно, не вами) приводят к каким-то забавным результатам, то и давите на эти точки. Только опишите проблему попонятнее и поинтереснее для простого обывателя (типа меня).

Вот я ещё раз попытался перечитать ваши первые сообщения в теме -- "ниасилил". "Многа букафф".
Я почему и настаиваю вот на этом самом, почему от вас всех и добиваюсь, чтобы вы в одном кратком сообщении чётко и ярко рассказали всю суть вашего открытия.

[Petrovich_Tot]

При этом параметры системы были следующие - масса пульсара 1,42*Мсолнца, масса компаньона 1,41*Мсолнца, большая полуось 971,437 тыс. км, эксцентриситет 0,617, долгота периастра 178,9 градуса, угол наклона орбиты к расчетной плоскости 42,9 градуса.


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


Варианты расчета___________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

Сергей, нет ли ошибки в данных dAlfaP? 1522189,5 - Это за 100 или за 1000 лет?

[Ser100]

Ну, если вы заметили, что данные (полученные, стесно, не вами) приводят к каким-то забавным результатам, то и давите на эти точки. Только опишите проблему попонятнее и поинтереснее для простого обывателя (типа меня).

Повторяю еще раз для самых понятливых. Все данные, которые приведены в моих статьях, получены именно мною, хотя для сравнения я привожу и данные других авторов. Вот, например, данные по вековым смещениям параметров орбит четырех внутренних планет, т.е. от Меркурия (1) до Марса (4), полученные мною на программах Solsys6 и Solsys7 с использованием данных полученных мною по имитаторам (астрономическим теориям) Солнечной системы New0, APC2 и JPL2, а также взятые непосредственно из моей астрономической теории Ser0, которая получена мною при обработке десятков тысяч данных оптических наблюдений за планетами и Солнцем различных обсерваторий мира в интервале примерно с 1800 по 1980 год, где учтены так же данные полученные мною по таблицам Птолемея, АльБаттани, АльХорезми и других древних астрономов. А так же я привожу и данные полученные мною на математических моделях Солнечной системы с использованием тех же моих программ Solsys6 и Solsys7.

Здесь JPL2 (эфемериды НАСА DE405), New0 (теория Ньюкома без периодических возмущений по данным книги Справочное руководство по небесной механике и астродинамике, под редакцией Г.Н. Дубошина, издание второе, дополненное и переработанное - М.: Наука, 1976. – 864 с. //см. стр. 487-494//), АРС2 (теория аналогичная теории Ньюкома с периодическими возмущениями из книги Астрономия с персональным компьютером, под редакцией О.Монтенбрук, Т.Пфлегер, пер. с нем. М.: Мир, 1993 – 279 с.), Ser0 (моя теория описанная в статье Кинематическая теория планет http://modsys.narod.ru/Stat/Stat_Est/Kinematik/Kinematik1.html ), а dAlfaP - смещение перигелия, dAlfaU - смещение узла восхождения, dBetta - смещение угла наклона, dEks - смещение (изменение) эксцентриситета орбит. При этом, все смещения получены мною в фиксированной эклиптике эпохи J2000.
   
__________Сглаженные теориями данные наблюдений________Модели Ньютона
Параметр __New0______APC2______Ser0______JPL2________/R^2_______/R^n*
dAlfaP1___+570,72____+574,02____+578,00____+572,20_____+529,64_____+573,38
dAlfaU1___-452,17_____-450,34____-433,20____-449,95______-450,37______-450,31
dBetta1____-21,43______-21,44_____-19,84_____-21,44_______-21,45______-21,45
dEks1_____+20,55_____+20,52_____+20,10_____+20,51______+20,43______+20,43

dAlfaP2____+54,46_____+54,74_____+30,30_____+40,54______+41,63______+57,32
dAlfaU2____-999,88____-1011,04____-996,50____-998,37______-998,97_____-998,92
dBetta2_____-2,50______-2,49_______-2,83______-2,50________-2,53_______-2,53
dEks2______-47,98_____-48,09_____-45,33______-48,37_______-49,46______-49,47

dAlfaP3___+1160,54___+1185,68___+1141,60___+1156,74_____+1140,88___+1151,45
dAlfaU3____-869,86____-872,65_____-868,00____-858,83______-660,48_____-660,26
dBetta3_____-47,18_____-47,14______-47,45_____-47,19_______-47,25______-47,25
dEks3______-41,42_____-41,64______-41,57_____-41,81_______-43,38______-43,36

dAlfaP4___+1601,07___+1599,64___+1582,70___+1599,96_____+1596,82___+1602,43
dAlfaU4___-1054,62___-1058,58____-1028,00____-1050,28_____-1052,33___-1052,34
dBetta4____-29,11_____-29,16______-29,39______-28,93_______-29,00_____-29,00
dEks4_____+92,10_____+92,07_____+95,53______+91,42______+96,21_____+96,17

* - показатель степени n во второй модели в формуле Ньютона F=G*m*M/R^n брался таким, каким он был у Ньюкома, т.е. n=2,0000001612 (у Холла было n=2,0000001574).

Вот я ещё раз попытался перечитать ваши первые сообщения в теме -- "ниасилил". "Многа букафф".
Я почему и настаиваю вот на этом самом, почему от вас всех и добиваюсь, чтобы вы в одном кратком сообщении чётко и ярко рассказали всю суть вашего открытия.

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

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

[Ser100]

Сергей, нет ли ошибки в данных dAlfaP? 1522189,5 - Это за 100 или за 1000 лет?

Это за 100 лет, но в угловых секундах, а в градусах будет 422.

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

[BJIaquMup]

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

Меня интересует, почему именно 100. А не 10 и не 1000.
Разумеется, ничего не понял.

Слушайте, по тутошней ссылке на "кинематик1" html, у вас проделана чудовищная работа, огромная просто. Почему вы так далеко запрятали ссылку на эту страницу с главного своего сайта? Я вам разрешаю открыть здесь столько тем, сколько вы захотите. Думается, надо открыть специальную тему, где бы было чётко и ясно всё разрисовано, вот эта таблица, её надо сделать красиво. Я вам помогу в этом, пожалуй.. Чуть позже. Сейчас надо самому разобраться кое с чем.  У меня тут тоже небольшой затык - небольшая правка основной формулы.  Не так-то просто всё.
Короче, вот эту всю таблицу надо представить более красиво. Почему у вас там ВАШИ данные где-то посредине? По-моему, данные Ser0 должны быть выделены в отдельную колонку. А нето всё как каша перемешано. Я не вижу, чем ваши данные отличаются. Вот так, на первый взгляд. Ну должны же вы как-то объяснить эту вашу скорость 100. Вы ведь этого хотите, я правильно понял? Чтобы из ваших рассчётов, из таблицы, было именно это видно?

[Petrovich_Tot]

Это за 100 лет, но в угловых секундах, а в градусах будет 422.

т.е. в угловых секундах за секунду - 0.000482344
период обращения чему равен и как определялись массы?

[Ser100]

Меня интересует, почему именно 100. А не 10 и не 1000.
Разумеется, ничего не понял.

За 100 лет потому что у астрономов принято давать смещения параметров орбит планет за век и называются эти смещения вековые. А, если ничего не поняли, то начинайте внимательно читать все мои статьи по этому вопросу начиная со статьи О равноденствиях Гиппарха и Птолемея http://modsys.narod.ru/Stat/Stat_Est/Hipp_Ptolo/Hipp_Ptolo2.html .

Слушайте, по тутошней ссылке на "кинематик1" html, у вас проделана чудовищная работа, огромная просто. Почему вы так далеко запрятали ссылку на эту страницу с главного своего сайта? Я вам разрешаю открыть здесь столько тем, сколько вы захотите. Думается, надо открыть специальную тему, где бы было чётко и ясно всё разрисовано, вот эта таблица, её надо сделать красиво. Я вам помогу в этом, пожалуй..

Ссылка на эту статью находится там же, где и ссылки на все остальные статьи http://modsys.narod.ru/Stat/statii.html , т.е. никуда она не запрятана. Вот только обсуждать свою методику обработки данных оптических наблюдений или результатов моделирования Солнечной системы (или системы двойного пульсара, о чем спрашивает Петрович) я сейчас не собираюсь, т.к. для себя я этот вопрос решил давным-давно и все изложено в этой статье, а обсуждаю на форумах я только или какие то спорные для меня моменты или недавно законченные работы (особенно, если я их собираюсь продолжать, т.к. во время обсуждения всплывают неизвестные мне эксперименты).

Короче, вот эту всю таблицу надо представить более красиво. Почему у вас там ВАШИ данные где-то посредине? По-моему, данные Ser0 должны быть выделены в отдельную колонку. А нето всё как каша перемешано. Я не вижу, чем ваши данные отличаются. Вот так, на первый взгляд. Ну должны же вы как-то объяснить эту вашу скорость 100. Вы ведь этого хотите, я правильно понял? Чтобы из ваших рассчётов, из таблицы, было именно это видно?

Отличие видно из других таблиц. Так в статье Аномальные смещения орбит планет http://modsys.narod.ru/Stat/Stat_Est/Anomal/Anomal1.html я привожу таблицы, где сравниваю аномальные смещения полученные Ньюкомом по его данным и аномальные смещения полученные мною. А затем уже в статье Влияние скорости гравитации на изменения параметров орбит планет http://modsys.narod.ru/Stat/Stat_Est/Vlijanie3.html я и пытался найти эту скорость гравитации такой, чтобы с применением моей теории потенциалов запаздывающих по координатам не было этих аномальных смещений. И в этой статье тоже куча интересных таблиц (одну из них сейчас изучает Петрович), но мне так и не удалось добиться, чтобы не было аномальных смещений сразу по всем параметрам и для всех планет. Так что тут тоже особо пока обсуждать нечего. Вот когда мне удастся провести наблюдения за двойным пульсаром В1913+16, вот тогда возможно и появится повод для нового обсуждения, а старое обсуждение уже было.

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