Положение Солнца на небе

[xd]

Вы перепутали тангенс и арктангенс.
Возведение в квадрат с помощью Math.Pow - тоже то ещё варварство.
Вы сколько классов школы закончили?

[Sergey007]

давно, заканчивал и уже все забыл (

Возведение в квадрат с помощью Math.Pow - тоже то ещё варварство.

Как правильнее, в чем заключается варварство? ))
Да точно, но тоже результат неправильный.
Ra = Math.Atan(Y1 / X1);   результат 1,5406686819811 а должно быть 88.2738131051079
Dec = Math.Atan( Z1 / Math.Sqrt(Math.Pow(X1, 2) + Math.Pow(Y1, 2)) );  результат 0,408926998903626 а должно быть  23.4297911659417

[Sergey007]

Снова я все перепутал. ( результат был в радианнах его нужно перевести в градусы (мне нужен в градусах)

[Upsilon]

Молодец, сами догадались .
При вычислении RA нужно использовать арктангенс с двумя аргументами или анализировать знаки числителя и знаменателя. Арктангенс с одним аргументом даст результат в диапазоне от -90 до 90 градусов, а диапазон изменения прямого восхождения - 360 градусов.

[Sergey007]

Учимся потихоньку
Все сходится кроме азимута в чем то еще ошибся.

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

double Az = Math.Atan(Math.Sin(Th * Math.PI / 180) * Math.Cos(Dec * Math.PI / 180) * Math.Cos(Lat * Math.PI / 180) / (Math.Sin(H * Math.PI / 180) * Math.Sin(Lat * Math.PI / 180) - Math.Sin(Dec * Math.PI / 180))) * 180 / Math.PI;
Все переменные сходятся. результат -84.577900800592289 должно быть 95.4220991915847

[xd]

А если вспомнить школьную тригонометрию? Есть метод Math.Atan2

[Sergey007]

Спасибо считает.)

[Sergey007]

Точно начала отсчета углов берется по расчетам из этой темы, для азимута с юга 0 градусов до запада 90гр север 180 восток 270 и юг 360. Правильно Az будет менять от 0 до 360 ?
Th =   71.0833609149958
Az =  95.4220991915847
H =  29.3205618652633
На пальцах я стою смотрю по компасу на юг это 0 градусов если я повернусь на запад это буде 90 градусов дальше повернусь на север 180 на восток 270 обратно на юг 360 крутился я по часовой стрелке правильно мыслю ?

[Sergey007]

Набросал программку по расчетам Tau и проверочным данным (большое спасибо ). Часовой угол измеряется от 0 до 360 азимут от -180 до 180гр высота над горизонтом от -90 до +90 но только у меня больше 66 -66 не получилось получить какие я координаты задавал и время результат больше 66гр не получается.  Проверьте пожалуйста результат правильный получается ?

[xd]

Переместите точку наблюдения ближе к экватору - получите больше 66 градусов.

[Sergey007]

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

[xd]

Правдоподобно.

[Sergey007]

спасибо. сейчас достал компас посмотрел оказывается у меня азимут расчетный с магнитным не совпадает примерно на 20 градусов такая большая разница может быть?
Направил компас на север  и  нарисовал тень от солнца по расчетам угол должен быть 60 градусов а по компасу от севера тень падает на 80.

[xd]

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

[Sergey007]

Да вы правы все влияет на компас, ножи вилки линейки  .
А как узнать поправку магнитного полюса относительно географического ?

[xd]

Открыть формулы сферической тригонометрии и обсчитать треугольник "Северный географичесекий полюс" - "Северный магнитный полюс" - "Текущее положение". Угол, противолежащий стороне "Северный магнитный полюс" - "Северный географический полюс" и будет для данной местности искомой поправкой.

[Sergey007]

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

[xd]

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

[Sergey007]

Да совпадает действительно для моих координат магнитное склонение 18гр.
Калькулятор http://www.ngdc.noaa.gov/geomag-web/#declination
Как этот калькулятор вычисляет это склонение может кому встречались формулы?

[xd]

Открыть формулы сферической тригонометрии и обсчитать треугольник "Северный географичесекий полюс" - "Северный магнитный полюс" - "Текущее положение". Угол, противолежащий стороне "Северный магнитный полюс" - "Северный географический полюс" и будет для данной местности искомой поправкой.

Это Вы так игнорируете сказанное или издеваетесь?