Электронный журнал «Труды МАИ». Выпуск № 67
www.mai.ru/science/trudy/
УДК 629.78
Теория запаздывания сигналов применительно к ГЛОНАСС и
GPS
Вовасов В. Е.
Московский авиационный институт (национальный исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3, 125993, Россия
e-mail: vovasov@list.ru
Аннотация
Приведена теория запаздывания сигналов, базирующаяся на классическом представлении о пространстве и времени. Показано, что в отличие от специальной теории относительности в преобразованиях теории запаздывания сигналов нет зависимости пространства и времени от скорости движения тел, хотя полученные аналитические зависимости для неучтенных временных задержек и истинного значения радиальной дальности в момент предшествия достаточно близки соответствующему уходу шкал времени движущихся источников сигналов и преобразованиям координат в специальной теории относительности. Полученные в рамках теории запаздывания сигналов выражения позволяют объяснить многие эффекты, до настоящего времени не объясненные физикой.
Ключевые слова: теория запаздывания сигналов, специальная теория относительности, ГЛОНАСС и GPS
С развитием ГЛОНАСС и GPS в руках инженеров появился высокоточный измерительный инструмент, о котором ученые, жившие даже 50 лет тому назад, могли только мечтать [3,4]. Так эксперимент, описанный в [2], показал зависимость частоты Доплера, измеряемой навигационным приемником, от ускорения приемника. Такое поведение частоты Доплера невозможно описать ни специальной теорией относительности (СТО), ни общей теорией относительности (ОТО). В [2] дано объяснение этому эффекту с точки зрения запаздывания электромагнитного сигнала излучаемого источником относительно приемника. Так как это явление проявляет себя в различных областях, то появилась необходимость описать его отдельно под названием теория запаздывания сигналов (ТЗС).
Постулат Ньютона об абсолютном пространстве содержит идею об абсолютно неподвижной системе отсчета. Считалось, что среди множества движущихся относительно друг друга инерциальных систем, имеется одна, преимущественная, связанная с абсолютным пространством, которая действительно неподвижна. В отличие от постулата Ньютона основным постулатом Эйнштейна является принцип постоянства скорости света - c в вакууме во всех инерциальных системах отсчета по всем направлениям. Математический анализ явлений, происходящих в инерциальных системах, который Эйнштейн провел на основе постулатов, один из которых приведен выше, показал, что преобразования Галилея вида
х' = X-V • t (1)
несовместимы с этими постулатами и должны быть заменены другими , которые теперь называют преобразованиями Лоренца [1,5].
Здесь V -скорость источника сигнала, совпадающая с осью х , t - время. Следствием постулатов Эйнштейна является зависимость пространства, времени и массы от скорости движения тел и ряд других положений. Эти
X - V • t
(2)
X
выводы возможны лишь в том случае, если мировая среда - эфир не существует в природе.
ТЗС базируется на модифицированном постулате Ньютона, т.е. система отсчета должна быть неподвижна относительно эфира. Причем любая масса создает в непосредственной близости от себя неподвижный, относительно этой массы, эфир, так Земля вращается в неподвижном эфире. Этот эффект называется увлечение эфира Землей [6]. ТЗС, основываясь на классическом представлении о пространстве без постулатов Эйнштейна, с учетом запаздывания сигнала от источника до приемника, приводит к преобразованиям координат похожим на преобразования Лоренца. Покажем это.
Пусть электромагнитные колебания распространяются со скоростью света (с =299792458 м/с) относительно некоторой неподвижной среды, которую принято называть эфиром. Введем некоторую инерциальную систему координат, оси которой совпадают с геоцентрической экваториальной, а центр координат неподвижен относительно эфира. В геоцентрической экваториальной системе координат начало совпадает с центром Земли, опорная плоскость этой системы координат совпадает с плоскостью экватора. Первая ось - направлена из центра Земли в точку весеннего равноденствия или точку Овна у ( у -астрономический знак созвездия Овна) . Вторая ось проходит через северный полюс Земли. Третья ось дополняет систему до правой. Таким образом выбранная система практически неподвижна относительно звезд.
Введем в рассмотрение текущий момент времени ? и будем говорить, что в каждый текущий момент времени ? осуществляется наблюдение за сигналом в точке его приема.
Введем в рассмотрение момент времени предшествия
КР ()= I - г (пр(г)) (3)
текущему моменту наблюдения ? на время задержки г (г(г)) излучаемых колебаний в пространстве.
В связи с тем, что дальность между приемником и источником колебаний изменяется, будет изменяться задержка г (гпр (г)) и поэтому момент времени
г (г) является функцией текущего времени ? в точке наблюдения.
Из (3) получаем следующее выражение для задержки г (гпр (г))
г (гпр(г))= г - гпр (г) (4)
В принятой неподвижной системе координат координату источника обозначим как хй (г1д (г)), а координату приемника как Х1 (г), аналогично
скорость источника в этой системе координат Хй (г16 (г)), а скорость
приемника Х1 (г), скорость же волны относительно этой неподвижной среды можно обозначить следующим образом V (г1д (г)), где модуль этого вектора равен скорости света - с.
Тогда расстояние между приемником и источником сигнала равно
(5)
Радиальная скорость сближения источника сигнала и приемника будет равна
РЕ! (^д (г)) :
(г„ (г))- х, (г)
Хш (г„ (г))
йг
- (г)
Рш! (г,5 (г))
ет Ц
(г*д (г ))•
х^ш (г„ (г))
1 йЦ г1д (г)) йг
- (г)
(6)
Р
е (гд(г))-Р-1 (гд (г))-РЕ (г1д (г))
ёт( г1д (г))
йг
Здесь
Ре (¿¡б (г)) = ет (г1д (г))• (г1д (г)) радиальная скорость источника сигнала, рР (^ (г))= ет (г[д (г)) • • (г) радиальная скорость приемника сигнала,
т
ёт Ц
(^б (t)) =
(х))- X, (t)
так называемые направляющие косинусы.
Рй1 С*1))
Исходя из предположения о существовании эфира, будем считать, что сложение скоростей соответствует преобразованиям Галилея и тогда радиальная скорость сближения волны и приемника будет равна
Ра (0)=п+&(?й(t)) (7)
Для задержки т( (t)) справедливо следующее выражение
4б (t))=
Рш ((6 (х)) _ (Х))-X, (Х)]т [X, (t))-X, (t)
у* (t).
Ра, ЕДХ)}
П + Р, Чг6 (t
ш
(8)
т
Дифференцируя левую и правую части (8) получаем уравнение для
))
вычисления производной
ж
(t))
[X, (^ (t))-X, (t)]т
XX, (Хй (Х
ёт^б (t))'
ин-^')
ж
(t))
^6 (хРш (1д (х))-(хй (t^ (Хй (t))-(Хй (t,
ш
(tí6 (t ))Xi (х)
(п + Р, (х)))2
(п + Р, ^6 (/)))
(9)
Решая уравнение (9), получаем выражение для производной
йт{г1д {t) )
Ж
в виде
Рш (О)-р, {tгд (О)-
11д
{п + Р, к, {х)))
(10)
Пусть ускорение, как у источника сигнала, так и приемника отсутствует, а
для вычисления координат приемника сигнала используется алгоритм
получения координат источников сигналов на момент предшествия (3).
5
Используя разложение в ряд Тейлора, можно получить требуемое значение с помощью следующего выражения
*е (г,д (г)) = хш (г-т(г1д(г))) = х& (г)-т(г„ (г))• (г) (11)
Значение величины задержки т(г1д (г)) можно получить из (8) в виде
(г ^РТТПТЙ (12)
п + р1 (гй (г))
йт(г16 (г))
Значение величины —-—— с учетом (10) и принятого допущения об
йг
отсутствии ускорений будет равно йт(1,д (г)) ре {(¡3 (г))- р! (/,д (г))
йг п+р! (г,д (г))+ р(гй (г)) Из-за эффекта запаздывания сигнала, исходя из выражения (6), выражение (11) можно записать в виде
хш (г,д (г))-х1 (г ) = хш (г-т(г,д (г)))-х1 (г ) = хш (г)-х1 (г )-
т(ггд (г ))• ххш (г )•
' йт(г,д (г)) • 1--
йг
(14)
хЕ (г)-х (г)-т(г,д (г))• Йй (г) + т(г,д (г))• ^^• Йй (г)
Будем считать, что движение источника сигнала и приемника происходит на одной прямой, тогда ет (гй (г ))=ет (г).
Умножим каждый член выражения (14 ) на вектор направляющих косинусов ет (гй (г)) и запишем
РЕ1 (г,д(г)) = Рш (г)-т(г,д(г))• Рш (г) + ^(г,д(г))• ^^• РЕ (г) (15)
, ч (г))
Подставляя значение т( г1д (г)) и —-—— в (15) и производя перегруппировку
4 ' йг
членов, получим
(. /л\ Рш (г)- т(г ,д (г^ Ре (г) (л
(л ' Е \ ю \ // ' ! \ ю \ // ге V /
П + Р! (г,д (г))+ РЕ (г дд (г))' П + А" (г,д (г))
}
Если движение происходит по координате -к, то с учетом введенных ранее обозначений можно записать
х = Ре1 (1));
х = Ре1 (1); у=ре (1); и = Д (^ (1 )) = 0;
Для практических задач релятивистские эффекты рассматриваются за время равное задержке сигнала, т.е.
1 = (1)) • И тогда
х(1))-^ _ X-1 • V V - и V л у?2
{1 -----—} {1 -V—}
п + и + V п + и п + С • V
Если -□ 1 то
п
*'= Х-' •" 3 (18)
{1 - -г + }
п п
При этих условиях можно записать
X -1 • V X -1 • V X -1 • V
X' =
Л V2, V3, Г^ V2 V3
{1 + Л1 -^г--^ + ^ д1 --г
п2 п3 \1 п2 2 • п2 п V п2 Отличие полученного выражения от преобразования Лоренца заключается в
2 3
наличие слагаемого -—- + — в знаменателе. Очевидно, что при ~ □ 1 эти
2 • п п п
выражения практически идентичны.
Эйнштейн с помощью математического анализа явлений, происходящих в инерциальных системах, пришел к преобразованиям Лоренца, т.е. не используя физику явления. В случае ТЗС, явление объясняется с точки зрения физики. В нашем случае 1 не любая величина, а задержка сигнала из-за конечного времени его распространения. Величина X - г • V значение
координаты х в момент предшествия г (г) без учета запаздывания сигнала,
т.е. фиктивная величина, а х' - значение координаты х в момент предшествия г (г) с учетом запаздывания сигнала, т.е. реальная величина. Рассмотрим
также известное преобразование времени, которое трактуется в СТО как уход часов в одной системе координат относительно ухода в другой
V
г--2 •х
С = -гС=. (19)
{4
Перепишем это выражение в более удобном для дальнейшего сравнения виде
г - г 'Л - ^ = V • х. (20)
V п с
Из-за эффекта запаздывания сигнала, исходя из выражения (4), выражение (11) можно записать в виде
йг
V
йт{г1д (г))
хе (г,д (г )) = хе (г -т(г,д (г ))) = хе (г )-т(г,д (г ))• хе (г)
(21)
хе (г)-т(г,д(г))• хе (г) + т(г,д(г))• е (г) = хе (г)-[т(г,д(г))-д^]• хе (г)
Интервал времени
ДтТ.д (г ))• ^¡М. (22)
можно трактовать как не учитываемый интервал времени, приводящий к ошибке расчета координат источника сигнала. Этот интервал времени равен Дт = г - г' в принятых в СТО обозначениях. Подставляя в это выражение
/ / ч\ йт(^д (г)) vп
значения т(г[д (г)) и —-—и считая -□ 1, получим 4 ' йг п
! г_ Рй (г,д (г)) Рш Ц(г)) - Р1 Ц (г)) х • V ^ х • V ^
п + Р! (г,д (г)) п + р, (гй (г))+ Рш (г,д (г)) с2 + с • V ~ с2
Очевидно, что и выражения (20) и (23) при малых скоростях источника сигнала практически идентичны.
Заметим, что в отличие от СТО, в преобразованиях ТЗС нет зависимости пространства и времени от скорости движения тел, хотя аналитические зависимости достаточно близки. Таким образом, в ТЗС возвращаемся к тому, что ход времени никак не связан с системой отсчета и является абсолютным, т.е. к тому, как считал Ньютон.
Аналогичные рассуждения по запаздыванию сигнала могут быть произведены и для явления гравитации. В этом случае источником излучения является носитель передатчика, например спутник ГЛОНАСС, а приемником, Земля, как превалирующая масса для околоземного пространства. В этом случае выражения (3)-(10) справедливы, но скорость гравитационного взаимодействия будем считать неизвестной и равной п , так как в ОТО эту величину принимают равной скорости света без каких либо доказательств. Ввиду того, что приемником гравитационных взаимодействий является цент поля, в нашем случае центр Земли, то введем для этого случая новые обозначения.
Тогда расстояние между приемником и источником сигнала равно
р* ^ (1 )) = >/[^е ^ (1))! [^Е ^ {1))] (24)
Радиальная скорость сближения источника сигнала и приемника будет равна
Ра((1)):
((1))
Хе [11д (1))
1 (^ ('))
(1,6 (1 ))•
Хе (1,6 (1))
Ра (1,б (1)) 1 (^ (1))^
(25)
= Р а
Здесь
Ра (^ (1)) = еТ (^ (1 ))• (¿¡б (1)) радиальная скорость источника сигнала,
Хе (1,6 (*)Г
а \ г6
(1,6 (1 )) =
направляющие косинусы.
Ра(1,6 (1))
Поэтому для запаздывания гравитационной волны и его производной получим выражения
Т
,('„ (' )) = p%^) (26)
п~
dta (ts (t))_ P (ts (t))
(27)
dt ñ+ps(ti6 (0) Отсюда по аналогии с (22) получим
л„ = „(,„ (,))•£%«>=рМ.тАМС (28)
dt п п + ps(*ta (*))
Переходя к обозначениям, принятым в [7], запишем
Pa(ts (t)) = p = a • (1-e • cos E) (29)
P^s (t)) = pP (30)
v ' V a 1 - e •cos E Здесь a - большая полуось орбиты источника гравитационного сигнала, e - эксцентриситет орбиты источника гравитационного сигнала, л - геоцентрическая константа гравитационного поля Земли, E -эксцентрическая аномалия источника гравитационного сигнала. Окончательно по аналогии с (23), получаем
t - = • e • sin E (31)
Очевидно, что интервал времени Лг2 можно трактовать, как не учитываемый интервал времени, приводящий к ошибке расчета координат источника гравитационного сигнала по известным формулам [3].
d2 x dü _Kx dt2 dx m
^ + dü = S (32)
dt2 dy m d2z dü-K
dt2 dz m
2 , ,2 , _2
где p = ^J x2 + y2 + z
и = — - потенциал гравитационного поля Земли, Р
т - масса источника гравитационных волн,
т
Kx, Ky, Kz -вектор возмущающих сил, вызванных различными причинами
помимо гравитационного притяжения источника центральным полем Земли. Указанный интервал (31) зависит от U - потенциала гравитационного поля Земли следующим образом
t"-1 = Ar2=VU•4•e J1"e•cosE • sinE (33)
Физически это не уход часов, в движущейся системе координат вызванный потенциалом гравитационного поля, как утверждается в ОТО, а не учтенный интервал времени, из-за игнорирования ТЗС.
Некоторые задачи, такие как расчет координат навигационных спутников ГЛОНАСС и GPS, по принятой эфемеридной информации, требуют учета как запаздывания гравитационного взаимодействия, так и запаздывания электромагнитного сигнала излучаемого спутником.
Введем в рассмотрение момент времени предшествия гравитационного взаимодействия предшествию электромагнитного сигнала излучаемого спутником
w ()= *ПРа)-1 hp,(*)) (34)
Обозначение момента времени предшествия электромагнитного сигнала излучаемого спутником оставим прежним.
С учетом выражения (14) для явления гравитационного взаимодействия запишем
Xé (W (t)) =
(з ^^
Xé (tid (t ))-Та (tidn„p (t ))• Xé (t„ (t )) + Ta (tidn p (t ))• dTS (t'£ (t )) • Xé ^ (t ))
С учетом того же выражения (14) для явления запаздывания электромагнитного сигнала подставим значение XÉ (tid (t)) в выражение (35) и
получим
ХЕ (к ПР (1 )) =
Хе (1 )-г( 1,6 (1 ))• ХХе (1 ) + т{116 (1 ))• ХХе (1)
( (36)
Т (^ (1 ))• Хе (1,6 (1 Ж (1 ))• • ХХе (1,6 (1 )) =
Хе (1)-г(1,6 (1 ))• ХХе (1 )-Га(1,6^ (1 ))• ХГе (1,6 (1)) + [Аг2 • Хе (1,6 (1 )) + Аг1 • ХХЕ (1)]
Здесь
А,2 {l}}•-т%{г}}--¿кМ(1 (37)
Р -1 Па Па Ра (1Й, пР (1))
2 а\г6 , пр
1а ^ р а ( ,6 , пр
Дг=г6 (Л- -Т(1т (1))~ Рм (16 ('^ Ре Ц (1))-Р/ Ц(1)). /38л
& п + р] (1,6
П +р/ (1г 6 (1))+ РЕ (1 16 ия ()
В связи с предполагаемым отсутствием или малой величиной ускорений источника как гравитационных взаимодействий, так и электромагнитных можно считать, что
Хе (1,6 (1))- ХЕ (1) •
И тогда выражение (36) можно переписать в виде
ХЕ (1,6,пр (1 ))-
Хе (1)-т(1,6(1 ))• хХе (1 )+г(1,6(1 ))• ^^ХЕ (1)
Т(1„пр (1 ))• ХХе (1 ) + Га(1,6,пр (1 ))• Щп^
(39)
Хе (1 )-
-1
Хе (1 )-,(1,6 (1 ))• Хе (1 )-Га(1,6,пр (1 ))• Хе (1) + [А,1 +АГ2]• ХХе (1)
Таким образом, для источников электромагнитных сигналов, на которые оказывает существенное влияние гравитационное поле, общая задержка от эффектов ТЗС будет равна
Ат - А, + Ат2 -
Ре1 (16 (1)) Ре Ц (1))- Р1 Ц (1)) | Ра(1,6, пр(1 ^ РаЦ пр(1)) (40)
п + Р] (1)) п +Р] (1))+ Ре (1)) па па + Рапр(1))
Преобразование же координат нужно осуществлять в два этапа. Сначала учитывая запаздывание электромагнитного сигнала, можно получить
Л /л\ Ре! ({)- т({о ({))• Ре ({) (л л л
РШ ^ (Г))% Ре ({« *))-Р, ^ ({)) Ре () (41) п + р! (/„, О))+Ре {о ({)) ятрЛ^^
А затем
Ре! -({))=^тР^ЖРШ, (42)
{ - --тТ.-ттЛ--~-}
Ги + А (7)) п~
а га\ю, прх // а
Здесь
Ра ({¡6, пр ({ ))=ет ({ 6п пр ({ ))• Хе ({) радиальная скорость источника сигнала в момент предшествия предшествию, т(, м) \Хе (<бппр({))-X ({)]Т
е пр ({))=-Г-ТТЛ-
р Ра {6, пр ({ ))
На первый взгляд, полученные выражения в рамках ТЗС, не существенно отличаются от известных соответствующих выражений в ОТО и СТО, но позволяют объяснить многие эффекты, до настоящего времени не объясненные физикой или объясненные в рамках гипотез. Один из эффектов упомянут в начале данной статьи. Математическому описанию таких эффектов будут посвящены последующие публикации. Выводы
1. ТЗС, основываясь на классическом представлении о пространстве и времени, с учетом запаздывания сигнала от источника до приемника, приводит к преобразованиям координат и времени, похожим на преобразования Лоренца.
2. ТЗС позволяет учесть, как запаздывание электромагнитного сигнала, излучаемого источником и принимаемого приемником, так и запаздывание гравитационного взаимодействия между его источником и центральным полем тяготения.
3. ТЗС объясняет появление неучтенных временных задержек и истинное значение радиальной дальности в момент предшествия с точки зрения физики. В СТО эти явления проявляются как уход шкал времени
движущихся источников сигналов и преобразование координат, и выводятся только с помощью математического анализа.
Библиографический список
1. Жданов Л.С., Учебник по физике для средних специальных учебных заведений. Изд. 2-е,стереотипное, М: Наука, 1977, 592 стр. с ил.
2. Вовасов В.Е. - Доплеровское смещение частоты для навигационных приемников, установленных на высокодинамичных объектах.
3. Журнал "Телекоммуникации ", №11, 2012 г., стр.24-31.
4. ГЛОНАСС принципы построения и функционирования/ Под ред. Перова А.И. Харисова В.Н. Изд. 4-е, перераб. - М: Радиотехника, 2010, 800 c., ил.
5. Поваляев А.А., Спутниковые радионавигационные системы: время, показания часов, формирование измерений и определение относительных координат. - М.: Радиотехника, 2008, 328 c., ил.
6. Иродов И.Е., Основные законы механики, Учебное пособие, М. "Высшая школа", 1985, 251 с.
7. Ацюковский В.А., Логические и экспериментальные основы теории относительности, М., Издательство МПИ, 1990, 55 стр.
8. Relativity in the Global Positioning System, Neil Ashby, Dept. of Physics, University of Colorado Boulder, CO 80309-0390 U.S.A., 2003-01-28