Телескопы покупают здесь


A A A A Автор Тема: пересчет  (Прочитано 2017 раз)

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

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
пересчет
« : 06 Мар 2010 [15:56:15] »
не могу найти как пересчитать из геоцентрической инерциальной СК в обычную долготу и широту.Может кто знает где или как

Tau

  • Гость
Re: пересчет
« Ответ #1 : 06 Мар 2010 [16:15:32] »
Поверните инерциальный вектор (X,Y,Z) вокруг оси Z на угол, равный -S (звездное гринвическое время) и разложите его на углы.
Или:
Долгота = arctg(Y/X) - S
Широта = arcsin(Z/sqrt(X^2+Y^2+Z^2))  ' широта в гринвической и инерциальной СК одинакова.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #2 : 06 Мар 2010 [16:31:47] »
А можно еще вопросик,как S расчитать...и вобще можно ли все это делать вариант для КА

Tau

  • Гость
Re: пересчет
« Ответ #3 : 06 Мар 2010 [17:55:48] »
Для КА это называется долготой и  широтой подспутниковой точки. Геодезическая широта вычисляет по-другому. Смотрите прицепляемый файл.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #4 : 07 Мар 2010 [13:28:41] »
спасибо вам огромное....только чуть не понятно с последней частью Фк и Фк-1

Tau

  • Гость
Re: пересчет
« Ответ #5 : 07 Мар 2010 [13:43:48] »
Вычисляете Ф0 - это геоцентрическая широта (если бы  Земля была шаром, то эта широта была бы широтой той точки на поверхности Земли, в которой спутник находится в зените). Но так как Земля - не шар, а эллипсоид, то нужно вычислить не геоцентрическую широту, а геодезическую, поэтому Ф0 - это первой приближение. Вычислили, теперь вычисляете C  и Ф1 (используя Ф0 в качестве Ф(k-1)) . Сравниваете, насколько отличаются Ф0 и Ф1. Если отличаются сильно, то вычисляете С  и Ф2 (используя Ф1 в качестве Ф(k-1)). Сравниваете Ф1 и Ф2. И так до тех пор, пока Ф(k) и Ф(k-1) не совпадут с нужной точностью. Вообщем, Ф(k-1) - это предыдущее вычисленное Ф, а Ф(k) - новое вычисленное Ф.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #6 : 07 Мар 2010 [14:05:04] »
а если считать землю сферичной тогда так грубо можно?

Tau

  • Гость
Re: пересчет
« Ответ #7 : 07 Мар 2010 [14:10:24] »
Максимальное расхождение между геодезической и геоцентрической широтами = 0.18 градуса, если такая точность устраивает, то можно :).
Арктангенсы-то не забудьте. Причем при вычислении долготы арктангенс должен быть с двумя аргументами - числителем и знаменателем. Возьмете простой арктангенс - будет ерунда.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #8 : 07 Мар 2010 [14:24:13] »
или  возможно брать ту же формулу для долготы но без -S

Tau

  • Гость
Re: пересчет
« Ответ #9 : 07 Мар 2010 [14:27:11] »
Если у вас вектор спутника действительно инерциальный, то без S долготу вы не получите. Вы получите геоцентрическое прямое восхождение спутника.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #10 : 07 Мар 2010 [14:33:06] »
эта S все портит :) она у мменя получается где то 2369478.2 так и должно или я считать не умею?

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #11 : 07 Мар 2010 [14:39:42] »
а если я просто Y\X  а потом от результата возьму арктангенс?

Tau

  • Гость
Re: пересчет
« Ответ #12 : 07 Мар 2010 [15:03:32] »
Цитата
она у меня получается где то 2369478.2 так и должно
Нет, так быть не должно. Число получится большое, но не настолько. Чтобы вогнать его в диапазон от 0 до 360 градусов, отнимайте 360 радусов, пока число не станет меньше 360.
Цитата
а если я просто Y\X  а потом от результата возьму арктангенс?
Без звездного времени вы получите ошибку, равную звездному времени.

Тест.
X = -11863.80783
Y = -6250.023262
Z = 21723.26637
4 марта 2010 года, всемирное время UT 23 часа 15 минут
MD = 55259
T = 0.101697467488022
S = 4111.35260997684 = 151.352609976844 градусов
Долгота =  56.4282914538887 градусов
Геоцентрическая широта =  58.313626370675 градусов

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #13 : 07 Мар 2010 [15:22:44] »
а где если не секрет UT yнайти или расчитать

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #14 : 07 Мар 2010 [15:27:46] »
MD это юлианская дата?а то у меня дрогое значение ее

Tau

  • Гость
Re: пересчет
« Ответ #15 : 07 Мар 2010 [15:36:00] »
Цитата
MD это юлианская дата?а то у меня дрогое значение ее
Да, MD - это юлианская дата. Проверяйте, как у вас происходит операция деления нацело - INT.
Цитата
а где если не секрет UT yнайти или расчитать
У вас дано какое-то время Time в дополнение к координатам спутника. Если это время местное зимнее, то UT = Time - n - 1
n - номер часового пояса.
Если это время местное летнее, то UT = Time - n - 2.
Например, я нахожусь во втором часовом поясе, сегодня время зимнее, на моих часах 15 часов 35 минут 20 секунд, следовательно всемирное время 12 часов 35 минут 20 секунд. Отняла 3 часа.

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #16 : 16 Мар 2010 [07:55:40] »
По этим формулам получается если брать х положительное то график функции от х возрастает,а если отрицательное то резко падает вниз как этого избежать

Tau

  • Гость
Re: пересчет
« Ответ #17 : 16 Мар 2010 [20:35:08] »
График какой функции от Х?

Оффлайн SicАвтор темы

  • Новичок
  • *
  • Сообщений: 24
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от Sic
Re: пересчет
« Ответ #18 : 16 Мар 2010 [22:55:07] »
Arctan(y/x)-s вот с ней проблема...трасса ка бежит а потом резко по этой координате уходит в сторону...может это быть связанно что для восточной и западной  долготы разные формулы

Tau

  • Гость
Re: пересчет
« Ответ #19 : 16 Мар 2010 [23:29:25] »
Долгота вычисляется одинаково: как западная,так и восточная. Долгота - она вообще одна, западная или восточная - это просто способ выражения относительно гринвического меридиана. Долгота измеряется в диапазоне от 0 до 360 градусов. Если перевести долготу в диапазон от -180 до 180 градусов, то отрицательная долгота считается западной, положительная - восточной.
Вы рисуете долготу как функцию от Х? А Y неужели константа?
Если не разберетесь, то конкретный пример покажите - что там куда прыгает :)