Киевский клуб любителей астрономии "Астрополис"

astromagazin.net
* *
Добро пожаловать, Гость. Пожалуйста, войдите или зарегистрируйтесь.
11 Декабря 2017, 02:22:31


Автор Тема: подскажите формулу расчета  (Прочитано 6634 раз)

0 Пользователей и 1 Гость просматривают эту тему.

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
подскажите формулу расчета
« : 19 Января 2008, 14:23:53 »
. - .

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

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #1 : 19 Января 2008, 14:58:48 »
. - .

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

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

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

SP

  • Клуб Астрополис, Модератор
  • Оффлайн Оффлайн
  • Сообщений: 3362
  • Благодарностей: 383
Re: подскажите формулу расчета
« Ответ #2 : 19 Января 2008, 19:02:06 »
. - .

Цитировать
Если мне объясните, что сие значит - угол стояния солнца
По идее, это  угловая высота Солнца над горизонтом в определённый момент времени.
Но зачем морочиться над решением этой задачи, если есть программы-планетарии? Разве что из любви к процессу...
Записан
бино 20х80, МТО-1000, "Мицар", самодельные 150, 265, 500-мм добсоны, SW-120/600, SW-25012 EQ6, QHY-5,-6, Canon 350Da
300 комет, >1000 дипов

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #3 : 19 Января 2008, 19:46:51 »
. - .

Может быть  ;D

К примеру,  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
« Последнее редактирование: 19 Января 2008, 19:54:08 от Олег Малый »
Записан

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
Re: подскажите формулу расчета
« Ответ #4 : 19 Января 2008, 22:23:05 »
. - .

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

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

KMM

  • Оффлайн Оффлайн
  • Сообщений: 2192
  • Благодарностей: 28
Re: подскажите формулу расчета
« Ответ #5 : 19 Января 2008, 22:56:31 »
. - .

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

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

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
Re: подскажите формулу расчета
« Ответ #6 : 19 Января 2008, 23:06:35 »
. - .

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

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #7 : 19 Января 2008, 23:18:41 »
. - .

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

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

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
Re: подскажите формулу расчета
« Ответ #8 : 19 Января 2008, 23:28:43 »
. - .

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

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #9 : 19 Января 2008, 23:37:28 »
. - .

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

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

  • Оффлайн Оффлайн
  • Сообщений: 2192
  • Благодарностей: 28
Re: подскажите формулу расчета
« Ответ #10 : 19 Января 2008, 23:44:31 »
. - .

Вот и я ж об чём!
Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #11 : 19 Января 2008, 23:49:42 »
. - .

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

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
Re: подскажите формулу расчета
« Ответ #12 : 20 Января 2008, 20:49:48 »
. - .

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

tlgleonid

  • Moderator
  • Оффлайн Оффлайн
  • Сообщений: 6618
  • Благодарностей: 404
Re: подскажите формулу расчета
« Ответ #13 : 20 Января 2008, 21:43:37 »
. - .

Вы, конечно, понимаете, что такая картинка будет справедлива 1-2 дня в году?
Какой тогда в ней практический смысл?
Записан
Задокументировано наблюдение  >1500 Deepsky объектов.
ТАЛ75R, GSO 6" 1:4 +HEQ5Pro SynScan+Canon1000D+QHY6+фильтры R,V+Юпитер21М+QHY5, 265мм и 415 мм Добсоны HandMade+o3+hb+uhc+...
FAQ по любительской астрономии

uncle

  • Оффлайн Оффлайн
  • Сообщений: 7
  • Благодарностей: 0
Re: подскажите формулу расчета
« Ответ #14 : 20 Января 2008, 22:07:23 »
. - .

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

2
Какой тогда в ней практический смысл?
а это уже специфика моей работы

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

« Последнее редактирование: 20 Января 2008, 22:11:53 от uncle »
Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #15 : 20 Января 2008, 22:31:51 »
. - .

Чтобы вычислить горизонтальные координаты (высота и азимут) Солнца для любого момента времени светового дня, сперва необходимо вычислить экваториальные координаты (прямое восхождение 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.
Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #16 : 20 Января 2008, 22:38:30 »
. - .

Как вычислить 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.



Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #17 : 20 Января 2008, 22:40:23 »
. - .

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

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: - широта места наблюдения.
Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #18 : 20 Января 2008, 22:42:42 »
. - .

Как вычислить местное звёздное 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).   
Записан

Олег Малый

  • Гость
Re: подскажите формулу расчета
« Ответ #19 : 20 Января 2008, 22:44:44 »
. - .

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