Здравствуйте!Помогите с решением задачки. Дано: высота светила=h, азимут=A и экваториальные координаты светила склонение=dec, прямое восхождение=alp.Нужно найти мои географические координаты - долготу = lam, широту=phi.Широту, понятно, можно найти по полярной phi=d... Но мне нужен точный алгоритм для всех звезд. Добрался пока до такого решения:Из теоремы синусов для сферического треугольника:sint=cos(h)*sin(A)/cos(dec); t=asin(t); lam=alp+t-s; (поскольку, с другой стороны, t=s-alp+lam)Вроде бы программа правильно считает, но не для всех долгот, для долгот, меньших где-то 25 градусов начинаются глюки.. короче правильные ответы получаются где-то по половине окружности долгот (криво как-то сказанул).Теперь по широте. Из теоремы косинусов:sinh=sin(dec)*sin(phi)+cos(dec )*cos(phi)*cos(t); Все было бы хорошо, но, как видно, справа одновременно и синус и косинус широты, т.е. явно phi не выразить. Поэтому я попробовал решить это уравнение методом итерации, типа:G=0; do{pp=G; G=asin((sin(h) - cos(dec)*cos(pp)*cos(t)) / sin(dec)); }while(fabs(pp-G)>10e-6); Оно решается, но опять же не для все широт. Например, с одними звездами задача решается на любой широте, а с другими вообще не решается.Понятно, что тут проблемы с математикой. Типа неопределенности.Я пробовал вывести удобное выражение для получения аналитического выражения(а не численно итерацие), но ничего хорошего не вышло.В принципе, задача должна быть решаема, раз уж для каких-то звезд и для каких-то широт решения получаются.Искал в интернете, нигде алгоритма нет. Неужели все решают навигационную задачу, только по высоте звезды, во время ее кульминации?Спасибо.
 | 
Линейный вид
