Друк

[Гість]

Ласкаво просимо, Гість. Будь ласка, увійдіть або зареєструйтеся.
19 Квітня 2024, 19:21:24

1 година 1 день 1 тиждень 1 місяць Назавжди

Увійти

[Новини]

•    Новини
•      Астрономія
•      Астрозаходи
•    Життя клуба
•      Новини
•      Проекти
•      Заходи
•      Звіти про спостереження
•      Документи
•    Статті
•      Обладнання та телескопи
•      Спостереження
•      Астрофотографія
•      Астрософт
•      Астрономія в Україні
•      Астрономи шуткують
•      Інше
•     
•      Файловий архів

Наш канал

Наші сайти

Астро-сайти

Програми

Наш канал

Наші сайти

Астро-сайти

Програми

[uncle]

Кто знает подскажите формулу расчета положенния солнца.
Задача следующая: известно расчетное время, известен угол стояния солнца, известен азимут. требуется определить угол стояния и азимут через Х часов.

[Олег Малый]

Кто знает подскажите формулу расчета положенния солнца.
Задача следующая: известно расчетное время, известен угол стояния солнца, известен азимут. требуется определить угол стояния и азимут через Х часов.

Если мне объясните, что сие значит - угол стояния солнца, то я тогда смогу помочь с решением данном задачи.

А лучше сформулируйте задачу целиком.

[SP]

Если мне объясните, что сие значит - угол стояния солнца

По идее, это  угловая высота Солнца над горизонтом в определённый момент времени.
Но зачем морочиться над решением этой задачи, если есть программы-планетарии? Разве что из любви к процессу...

[Олег Малый]

Может быть 

К примеру,  20 июня 2008 г., т.е. в день летнего солнцестояния , в момент киевского времени 11h41m00s, для г. Запорожье, азимут Солнца равен  0h 00m 06s, а его топоцентрическая высота равна  +65° 36' 42".

Также смотрите темы:
1. Вычисление экваториальных координат Солнца (с невысокой точностью (порядка 0.01 градуса = 6’) http://astroclub.biz/forum/index.php?topic=213.0
2. Вычисление экваториальных  топоцентрических координат из геоцентрических http://astroclub.biz/forum/index.php?topic=418.0
3. Формулы перехода между системами координат http://astroclub.biz/forum/index.php?topic=216.0

[uncle]

2 SP  нет не из любви к процессу
2 Олег Малый положение солнца в конкретное время меня не интересует.

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

[KMM]

2 Олег Малый положение солнца в конкретное время меня не интересует.

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

Не находите противоречия самому себе?

[uncle]

Не находите противоречия самому себе?

не нахожу, у меня есть координаты положения солнца полько в определенный час, а я хочу одновременно иметь графическое построение лучей за весь световой день

[Олег Малый]

Не находите противоречия самому себе?

не нахожу, у меня есть координаты положения солнца полько в определенный час, а я хочу одновременно иметь графическое построение лучей за весь световой день

Т.е. вам нужны горизонтальные координаты (высота и азимут) Солнца для любого момента времени светового дня?

[uncle]

Т.е. вам нужны горизонтальные координаты (высота и азимут) Солнца для любого момента времени светового дня?

Совершенно верно.
Извиняюсь за корявость формулировки просто я не спец в данной области.

[Олег Малый]

В таком случае, вам необходимо ознакомится с содержанием вот этих страниц:

1. Вычисление экваториальных координат Солнца (с невысокой точностью (порядка 0.01 градуса = 6’) http://astroclub.biz/forum/index.php?topic=213.0
2. Вычисление экваториальных  топоцентрических координат из геоцентрических http://astroclub.biz/forum/index.php?topic=418.0
3. Формулы перехода между системами координат http://astroclub.biz/forum/index.php?topic=216.0

[KMM]

Вот и я ж об чём!

[Олег Малый]

В данном случае, п. 2 Вычисление экваториальных  топоцентрических координат из геоцентрических можно пропустить, по причине условно-низких требований к точности.

[uncle]

в конечном итоге я должен получить вот такой график но интерактивно подстраиваемый под любой город

[tlgleonid]

Вы, конечно, понимаете, что такая картинка будет справедлива 1-2 дня в году?
Какой тогда в ней практический смысл?

[uncle]

1

Вы, конечно, понимаете, что такая картинка будет справедлива 1-2 дня в году?

только 1 день

2

Какой тогда в ней практический смысл?

а это уже специфика моей работы

СанПиН 2.2.1/2.1.1.1076-01
2.4. Нормативная продолжительность инсоляции устанавливается на определенные календарные периоды с учетом географической широты местности:
северная зона (севернее 58° с. ш.) -  22 апреля
центральная зона (58° с. ш. - 48° с. ш.) -  22 марта
южная зона (южнее 48° с. ш.) - 22 февраля

[Олег Малый]

Чтобы вычислить горизонтальные координаты (высота и азимут) Солнца для любого момента времени светового дня, сперва необходимо вычислить экваториальные координаты (прямое восхождение RA и склонение DE) Солнца для любого момента времени светового дня.

Когда экваториальные координаты Солнца определены, то тогда их можно перевести горизонтальные координаты.

Итак, прямое восхождение RA Солнца вычисляем по следующей формуле:

tan(RA) = cos(eps)*tan(theta)  (4.1)

а склонение DE Солнца вычисляем по формуле:

sin(DE) = sin(eps)*sin(theta)  (4.2)

Во всех приведённых формулах, T =  (JD - 2451545.0)/36525  (4.3)

Примечание: о вычислении JD смотрите тему "Вычисление юлианской даты JD".

Вычисляем наклон плоскости земного экватора eps (в градусах!)  к плоскости эклиптики

eps = 23+26/60+21.448/3600 - 46.8150/3600*T - 0.00059/3600*T*T + 0.001813*T*T*T  (4.4)

Вычисляем истинную долготу Солнца theta (в градусах!)  по формуле

 theta = Lo + Co (4.5)

где Lo - геометрическая средняя долгота Солнца (в градусах!) , приведённая к среднему текущему равноденствию, вычисляется по формуле

Lo = 280.46645 + 36000.76983*T +0.0003032*T*T  (4.6)

где Co - положение центра диска Солнца (в градусах!) , вычисляется по формуле

Co = (1.914600 - 0.004817*T - 0.00014*T*T)* sin( Mo) + (0.019993 - 0.000101*T)* sin(2* Mo) + 0.000290* sin(3* Mo)  (4.7)

а Mo - средняя аномалия Солнца (в градусах!) , вычисляемая по формуле

 Mo = 357.52910 +35999.05030*T -0.0001559*T*T -0.00000048*T*T*T.

[Олег Малый]

Как вычислить JD?

JD нам необходимо вычислить для T =  (JD - 2451545.0)/36525  (4.3).

Для вычисление юлианской даты JD для любого момента текущего времени (и места наблюдения, а не только для Гринвича!) Hour:Minute:Second для календарной даты Date.Month.Year, предлагается следующая формула:

JD = Int(365.25*(Year + 4716)) + Int(30.6001*(Month + 1)) + Date + B - 1525.5 - zone/24 + (Hour  +  Minute/60 + Second/3600)/24  (2.1)

где B = 2 - A +Int(A/4) (2.2)
      A = Int(Year/100)   (2.3)
      zone - временная зона места наблюдения (с учётом летнего времени).

Если Month = 1 или 2 (январь или февраль), то для формулы (2.1) Year = Year - 1, а Month = Month + 12.

Функция Int(x) определяет целую часть числа x.

Например: Int(3.56847325) = 3.

[Олег Малый]

Как перейти от экваториальных координат к горизонтальным?

1. находим местное звёздное S время для момента Т местного поясного времени (см. http://astroclub.biz/forum/index.php?topic=208.0

2. вычисляем часовой угол t по формуле t=S-RA интересующего нас небесного светила.

3. вычисляем зенитное расстояние z небесного светила по формуле

             cos(z)=sin(:fi:)*sin(DE)+cos(:fi:)*cos(DE)*cos(t)


4.  азимут A получим из формулы sin(z)*sin(A)=cos(DE)*sin(t)

где RA- прямое восхождение небесного светила
а DE - его склонение.

:fi: - широта места наблюдения.

[Олег Малый]

Как вычислить местное звёздное S время для момента Т местного поясного времени?

. Среднее звёздное время и его вычисление.

В 1982 г. для вычисления среднего звёздного времени для Гринвича на момент времени 0h UT, Международным Астрономическим Союзом, была принята следующая формула:

So = 6h41m50s.54841 + 8 640 184s.812866*T + 0s.093 104*T*T - 0s.000 006 2*T*T*T (1.1)

где T = (JD - 2451 545.0)/36525 (1.2)
это время в юлианских столетиях (по 36525 эфемеридных суток), прошедших от начала эпохи 2000,0 (1 января, 12h00m UT),
а JD - юлианская дата.

Для вычисления среднего звёздного времени для Гринвича на любой момент времени t, прошедший с момента времени 0h UT, к So следует прибавить величину t*1.002 737 909 35.   

Формулы для вычисления среднего звёздного времени для места наблюдения на любой момент tm местного времени:

А) для зимнего времени по Киеву:

St = So + (tm - 2)*1.00273790935 + long/15 (1.4a)

Б) для летнего времени по Киеву:

St = So + (tm - 3)*1.00273790935 + long/15 (1.4б)

где long - долгота места наблюдения в градусной мере,
а So из формулы (1.1).   

[Олег Малый]

Данные алгортмы изложены на форуме "Альтаира" и находятся в темах по ссылке http://astroclub.biz/forum/index.php?board=22.0