Вычисление элементов невозмущённой орбиты по наблюдениям спутника с пункта земной поверхности
Постановка задачи
Пусть с пункта земной поверхности, координаты которого известны, выполнены наблюдения искусственного спутника Земли и определены топоцентрические направления и расстояния до трёх его мгновенных положений. В результате вычислены геоёмного методов определения орбит. Наиболее широко используемым является классический метод Гаусса.
Исходные данные
Номер варианта Момент времени, с Координаты спутника в равноденственной системе координат
x, м y, м z, м
0 -9572849.511 9178349.804 20597123.981
3 120 -10016322.150 9082709.836 20428060.924
240 -10455884.323 8983523.894 20251022.560
Решение
1. Вычисляем геоцентрические прямые восхождения, склонения и модули векторов ИСЗ для трёх его положений
αS=arctgysxS
δS=arctgzsxS2+ys2
rS=xS2+ys2+zs2
α1=arctg -9572849.5119178349.804=136°12΄18,9˝
δ1=arctg20597123.981-9572849.5112+9178349.8042=57°13΄24,3˝
r1=-9572849.5112+9178349.8042+20597123.9812=24497409.438 м
α2=arctg -10016322.1509082709.836=137°47΄54,8˝
δ2=arctg20428060.924-10016322.1502+9082709.8362=56°29΄59,0˝
r2=-10016322.1502+9082709.8362+20428060.9242=24497510.088 м
α3=arctg -10455884.3238983523.894=139°19΄52,9˝
δ3=arctg20251022.560-10455884.3232+8983523.8942=55°45΄23,1˝
r3=-10455884.3232+8983523.8942+20251022.5602=24497614.848 м
2. Вычисляем угол наклона плоскости орбиты ИСЗ к плоскости экватора Земли по формуле
сos i=x1∙y3-x3∙y1z3∙y1-z1∙y32+x3∙z1-x1∙z32+x1∙y3-x3∙y12=AB
A=-9572849.511∙8983523.894--10455884.323∙9178349.804=
A=9969841510919.020
B=20251022.560∙9178349.804-20597123.981∙8983523.8942+
-10455884.323∙20597123.981--9572849.511∙20251022.5602+
(((-9572849.511)∙(8983523.894)-(-10455884.323)∙(9178349.8048))))^2
B=23714903118749.400
сos i=AB=9969841510919.02023714903118749.400=0.420404058
i=65°08΄23,6˝
3
. Вычисляем долготу восходящего узла Ω орбиты ИСЗ, используя формулы
Sin α1- Ω=ctg i·tgδ1=
= ctg65°.08΄23,6˝·tg57°.13΄24,3˝=0.719605902
α1-Ω=133°58΄41,0˝
Ω=136°12΄18,9˝-133°58΄41,0˝=2°13΄37,9˝
tgα1- Ω=tgδ1∙sinα3-α1tgδ3-tgδ1∙cosα3-α1=
=tg57°13΄24,3˝·sin139°19΄52,9˝-136°12΄18,9˝tg55°45΄23,1˝-tg57°13΄24,3˝·cos139°19΄52,9˝-136°12΄18,9˝
tgα1- Ω=-1.0363246
α1- Ω=133°58΄41,0˝
Ω=136°12΄18,9˝-133°58΄41,0˝=2°13΄37,9˝
4. Вычисляем значения аргументов широты для трёх положений ИСЗ
cosuS=cosaS-Ω∙cosδS
sinuS=sinδSsini∙
cosu1=cosa1-Ω∙cosδ1
cosu1=cos136°12΄18,9˝-2°13΄37,9˝·cos57°13΄24,3˝=-0.37591433
sinu1=sinδ1sini
sinu1=sin57°13΄24,3˝sin65°08΄23,6˝=0.926654421
u1=112°04΄51,0˝
cosu2=cosa2-Ω∙cosδ2
cosu2=cos137°47΄54,8˝-2°13΄37,9˝·cos56°29΄59,0˝=-0.394153595
sinu2=sinδ2sini
sinu2=sin56°29΄59,0˝sin65°08΄23,6˝=0.919044582
u2=113°12΄47,5˝
cosu3=cosa3-Ω∙cosδ3
cosu3=cos139°19΄52,9˝-2°13΄37,9˝·cos55°45΄23,1˝=-0.41223876
sinu3=sinδ3sini
sinu3=sin55°45΄23,1˝sin65°08΄23,6˝=0.91107585
u3=114°20΄44,0˝
5