Линия "Терминатора"

[Tau]

от его положения, тоже много зависит,

Что вы подразумеваете под положением Солнца?

Ааа.. да.. я как-то послушала местный эфир радиолюбителей. Ничего так мужики чатились на радиоволнах   

[RN3QGA]

Закат , рассвет, в зените
Да увлекательная весчь радио

[RN3QGA]

Уважаемая Tau, а не подскажите как для луны расчет терминатора сделать. Буду весьма признателен.
На УКВ многие работают через луну, думаю будет не лишним добавить хотя бы её положение (но лучше все)  в свою прогу.
С уважением RN3QGA.

[Tau]

Угловой радиус терминатора Луны 89 градусов, а не 90, как у большого круга. Проигнорируем этот факт и будем считать терминатор Луны большим кругом.
При вычислении терминатора Солнца проводились следующие вычисления:
1) Вычисление юлианской даты
2) Вычисление вектора положения Солнца в эклиптической системе координат
3) Вычисление вектора положения Солнца в экваториальной системе координат
4) Вычисление звездного гринвического времени
5) Поворот вектора Солнца в гринвическую систему координат
6) Разложение вектора Солнца на углы – полюс терминатора
7) Вычисление географических координат точек, лежащих на терминаторе.

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

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

T = (MD+UT/24 - 51544.5) / 36525
 q = pi / 180 
' эклиптическая долгота
 L = 218.32 + 481267.883 * T +
+6.29 * Sin(q * (134.9 + 477198.85 * T)) -
-1.27 * Sin(q * (259.2 - 413335.38 * T)) +
+ 0.66 * Sin(q * (235.7 + 890534.23 * T)) +
+0.21 * Sin(q * (269.9 + 954397.7 * T)) -
-0.19 * Sin(q * (357.5 + 35999.05 * T)) - 
-0.11 * Sin(q * (186.6 + 966404.05 * T))
' Эклиптическая широта
 F = 5.13 * Sin(q * (93.3 + 483202.03 * T)) +
+0.28 * Sin(q * (228.2 + 960400.87 * T)) -
-0.28 * Sin(q * (318.3 + 6003.18 * T)) - 
-0.17 * Sin(q * (217.6 - 407332.2 * T))
' Параллакс
 p = 0.9508 + 0.0518 * Cos(q * (134.9 + 477198.85 * T)) +
+0.0095 * Cos(q * (259.2 - 413335.38 * T)) +
+0.0078 * Cos(q * (235.7 + 890534.23 * T)) + 
+0.0028 * Cos(q * (269.9 + 954397.7 * T))
' Расстояние от Земли до Луны в радиусах Земли
 r = 1 / Sin(p * q)
 ER = 6378.136
 X = r * ER * Cos(q*L) * Cos(q*F)
 Y = r * ER * Sin(q*L) * Cos(q*F)
 Z = r * ER * Sin(q*F)
Это не самое точное вычисление, оно дает ошибку ~ 10 угловых минут, зато оно самое короткое.
А дальше пункты 3-7: вычисление экваториального вектора Луны, звездное время (если вы его уже вычислили для солнца, то второй раз на данный момент времени вычислять не нужно) , поворот вектора Луны в гринвическую систему координат, разложение на углы, вычисление географических координат терминатора Луны.
Тест:
11.02.2010 UT=5.033
X =  144273.846684896
Y = -378584.242770559
Z = -140.334392420613
Дальнейшие числа не привожу. Если у вас для Солнца правильно вычисляется, то и для Луны вычислится правильно. Картинка ниже - результат вышеприведенного теста.

Что вам еще необходимо? Восход,заход, кульминация Солнца и Луны для конкретного географического пункта?

[RN3QGA]

  Ну просто хочется крикнуь, - А давайте всеееее.  но думаю пока и это хватит. Понадобиться я обязательно спрошу.
Огромное Вам спасибо уважаемая Tau.

[kress]

RN3QGA, поздравляю
Правда, я так и не поняла, что человеку дает лицезрение терминатора? Информативность-то какая?

Дело не в информативности, а в.....не знаю. Tau, если Вас не затруднит дайте , пожалуйста, формулу нахождения азимута точки весеннего равнодействия.

[Tau]

kress, азимут точки с координатами RA=0, DEC=0 (точки весеннего равноденствия):
Az = 180/pi * arctg(tgS/sinF)
S - местное звездное время, S = S0 + L
S0 - гринвическое звездное время на нужный момент
L - географическая долгота места наблюдения
F - широта места наблюдения
На экваторе (при F=0) формула работать не будет.

[Tau]

Но если я вас неправильно поняла и вы хотите нарисовать точку весеннего равноденствия на карте, то рисуйте точку с географическими координатами: долгота = -S, широта = 0.
S - гринвическое звездное время

[kress]

Но если я вас неправильно поняла и вы хотите нарисовать точку весеннего равноденствия на карте, то рисуйте точку с географическими координатами: долгота = -S, широта = 0.
S - гринвическое звездное время

Большое спасибо.

[kress]

kress, азимут точки с координатами RA=0, DEC=0 (точки весеннего равноденствия):
Az = 180/pi * arctg(tgS/sinF)
S - местное звездное время, S = S0 + L
S0 - гринвическое звездное время на нужный момент
L - географическая долгота места наблюдения
F - широта места наблюдения
На экваторе (при F=0) формула работать не будет.

Tau, я рассчитал азимут по Вашей формуле, но полученное значение отличается от истинного(в StarCalcke) на 180 градусов.

[Tau]

Действительно, в половине случаев ерунда какая-то.
Так лучше:
Az = 180/pi * atn2(sinS,cosS*sinF)
atn2 - арктангенс с двумя аргументами

[kress]

Действительно, в половине случаев ерунда какая-то.
Так лучше:
Az = 180/pi * atn2(sinS,cosS*sinF)
atn2 - арктангенс с двумя аргументами

Хм, у меня результат аналогичный, что и по старой формуле и надо еще следить за знаком.

[Tau]

17 февраля 2010, Долгота = 30 градусов, Широта = 60 градусов

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

UT       Азимут
 0ч   176.378971250693
 1ч   193.680886596857
 2ч  210.411171232754
 3ч   226.10115615483
 4ч   240.672008915473
 5ч   254.344674002029
 6ч   267.496192748073
 7ч   280.562016197894
 8ч   293.985054426742
 9ч   308.173161460715
 10ч   323.422658434349
 11ч   339.787143933509
 12ч  356.947551901653
 13ч   14.2416680986013
 14ч   30.9432018306744
 15ч   46.5957397712078
 16ч   61.1324513111718
 17ч   74.7815519144397
 18ч   87.923218431595
 19ч   100.993789652243
 20ч   114.435777884126
 21ч   128.654882840984
 22ч   143.941818273455
 23ч   160.339665852684
 24ч   178.654083244087
Что не так?

[RN3QGA]

Доброго времени .
Уважаемая Tau,я рано радовался , начал двигать по месяцам и обнаружил, что сам терминатор ОК , а вот зарисовочка все время сверху карты, нашел вроде ошибку, что нужно менять, но теперь загвоздка другая, никак не могу разобраться с вычислением U (Ваш док по терминатору раздел Зарисовка ночной области),  определить либо градус либо знак cosU, что то получается все время 1, вроде не округляю, все в радианах, переменные в Double, а вот если просто U вычисляю без Cos,  то и знак появляется, но почему - то тогда, в одном месяце может быть и "-"  и "+";
Может наставите на путь истинный.
С уважением RN3QGA.

[Tau]

RN3QGA, когда вы зарисовываете ночную область, то для каждой точки карты (поверхности Земли) определяется угол между вектором этой точки и вектором солнца. Если угол больше 90 градусов (=косинус угла отрицательный), то данная точка поверхности Земли не освещена солнцем. Не поняла, в чем у вас трудности. Примерчик какой-нибудь конкретный приведите.

[RN3QGA]

Вот кусочек кода вычисляющий cosU:
   
    ао3=cos(A-Ao);
    uu3:=sin(F)*sin(Fo);
    uu2:=cos(F)*cos(Fo)*ao3;
    uu1:=uu3+uu2;
    Uu:=cos(uu1);


    Form_Terminator.Label1.Caption:=FloatToStr(Uu);

Uu все время равно 1.

[Tau]

И получился у вас такой ужас:
Uu = cos(sinF*sinFo+cosF*cosFo*cos(A-Ao))

У вас uu1 это и есть cosU ! Не надо вычислять косинус косинуса. Смотрите на знак uu1.
Если вы хотите вычислить угол U, то U = 180/PI*arccos(uu1)

Непонятно, как вы закрасили ночную область в ответе№ 57, если пользовались таким кодом 

[RN3QGA]

  Да его тогда еще и не было этого кода красилось ошибочно для одного параметра переменной.
А когда стал проверять - "Опа  яма ....".
В Вашем описании так:   cosU=sinF* ......  ну я из U cos  и вычислял 
но вот мне непонятно, в одном месяце могут быть разные знаки"-" и "+" -- это правильно???, или чего то опять на так?
С уважением RN3QGA.

[Tau]

в одном месяце могут быть разные знаки"-" и "+" -- это правильно???, или чего то опять на так?

Разные знаки для большей части точек земной поверхности должны быть не то что в разные месяцы, а в разное время одного дня. Только в областях полярного дня и ночи знак будет постоянен в течение нескольких месяцев, пока полярный день/ночь не закончится. На большей части Земли Солнце восходит и заходит каждый день, пока длится день в конкретной точке Земли - знак будет положительный, когда наступит ночь - знак станет отрицательным.

[RN3QGA]

Вот посмотрите пожалста по месяцам, будет ясно, что у меня не получается.