Бортовые вычислительные комплексы навигации и самолетовождения. Шивринский В.Н. - 32 стр.

UptoLike

Составители: 

33
здесь ϕширота, λдолгота места ЛА; αпрямое восхождение, δ
склонение, t – часовой угол светила; S
гр
звёздное гринвичское время.
Для определения азимута светила по вычисленному значению A
1
, не-
обходимо воспользоваться таблицей 2.1.
Таблица 2.1
Определение квадранта азимута светила
Знаки
B
1
sin(t)
Азимут А
+ +
2⋅π – A
1
– +
π – A
1
– –
π – A
1
+ – A
1
Вычисление высоты и азимута светила производится по следующему
алгоритму:
а) после ввода исходных данных вычисляется часовой угол светила t по
формуле (2.7);
б) вычисляется значение B из уравнения (2.5)
B = sin(ϕ)·sin(δ) + cos(ϕ)·cos(δ)·cos(t);
из-за погрешности вычисления тригонометрических функций мо-
дуль B может быть чуть больше 1;
в) для исключения аварийного останова, предусмотреть округление
B
до 1, если его модуль больше 1;
г) произвести вычисления
h = arcsin(B); D = sin(t); C = sin(t)·cos(δ);
д) по формуле (2.8) вычислить B
1
;
е) прежде чем поделить C на B
1
, необходимо убедиться, что B
1
не равно
нулю; в противном случае B
1
нужно присвоить значение 10
-9
; эта не-
большая величина (меньше погрешности вычисления тригонометри-
ческих функций) практически не вносит погрешностей в результат,
однако не вызывает аварийного останова программы;
ж) заканчивается вычисление азимута определением квадранта согласно
таблице 2.1.