способ оперативной синхронизации мер времени и частоты, размещенных на подвижных объектах

Формула изобретения

Способ оперативной синхронизации мер времени и частоты, размещенных на подвижных объектах, основанный на определении запросной и беззапросной дальности до подвижного объекта, отличающийся тем, что в беззапросном канале радиотехнической системы наземного объекта определяют момент времени достижения минимального расхождения бортовой и наземной шкал времени, в запросном канале радиотехнической системы наземного объекта определяют момент времени прохождения подвижного объекта через точку траверза, а относительную разность частот бортового (подвижного) и наземного (неподвижного) генераторов определяют по формуле



где f_н-б разность частот бортового (подвижного) и наземного (неподвижного) генераторов;

f_н - частота наземного генератора;

v_с - скорость подвижного объекта;

с - скорость света;

R_о - высота орбиты подвижного объекта;

t_min - момент времени достижения минимального расхождения бортовой и наземной шкал времени;

t_о - момент времени прохождения подвижного объекта через точку траверза.

Описание изобретения к патенту

Изобретение относится к измерительной технике и может быть использовано для определения частотно-временных параметров высокостабильных генераторов (мер времени и частоты), размещенных на подвижных объектах, например самолетах или навигационных космических аппаратах (НКА) системы "ГЛОНАСС".

Известен наиболее близкий к изобретению по техническому решению способ определения частотно-временных параметров высокостабильных генераторов, который по существу сводится к определению запросной и беззапросной дальности до НКА.

Рассмотрим подробнее способ определения шкал времени бортовой и наземной аппаратуры [1]. По запросу с Земли бортовая аппаратура излучает сигнал, несущий информацию о бортовой шкале времени. Излучается импульс, временное положение которого совпадает с секундной меткой бортовой шкалы. В наземной аппаратуре осуществляется измерение временного интервала между моментом приема сигнала с борта и ближней секундной отметкой наземной шкалы.

Очевидно, что поправка U к шкале бортового времени связана с измеренным интервалом t и временем распространения сигнала _p соотношением

U= t-_p. (1)

Учет непрерывно изменяющегося времени распространения сигнала может выполняться при обработке результатов измерений расхождения шкал времени в соответствии с соотношением (1). При этом обычно имеется возможность достаточно точного учета времени распространения сигнала, например по результатам измерения запросной дальности _p=D/C, где D - запросная (наклонная) дальность до НКА; C - скорость света.

Так как разработанные способы непосредственного сравнения частот при больших расстояниях обладают большой погрешностью [1], то для сравнения частот высокостабильных генераторов в настоящее время применяются способы, связанные с использованием с результатов измерений расхождений шкал времени, сформированных исследуемыми генераторами и реализуется двухкратным измерением расхождений шкал времени. Если U₁ и U₂ - суть расхождения шкал времени соответственно при первой и второй сверках, то искомая величина расхождения частот определяется соотношением



где

f_э - номинальное значение частоты исследуемых генераторов;

t - уход одной шкалы относительно другой за время T_н, обусловленный различием действительных значений частот.

Соотношение (2) используется при сравнении частот наземных неподвижных генераторов.

При сравнении этим способом частот бортового и наземного генераторов необходимо дополнительно учитывать эффекты Доплера первого и второго порядка и гравитационное смещение частоты. Порядок учета этих эффектов описан в [3], поэтому подробно на этом останавливаться не будем.

Из известных этот способ сравнения частот является наиболее высокоточным, но он обладает существенными недостатками: требует большого интервала времени измерения (накопления информации) и обладает малой оперативностью.

Анализ выражения (2) показывает, что средняя квадратическая погрешность определения относительной разности частот бортового и наземного генераторов обратно пропорциональна времени измерения (наблюдения НКА). Так, если требуется определить указанным способом относительную разность частот бортового и наземного генераторов со средней квадратической погрешностью (f/f)=110^-12,, при (t) 100нс,, то время наблюдения НКА должно составлять T_н = 27 ч, т.е. более суток, что снижает оперативность закладки частотно-временной информации на борт НКА. На коротких интервалах времени этот способ не обладает достаточной точностью. Так, при T_н = 1 ч и (t) 100нс величина (f/f)~110^-10, что не обеспечивает потребности современных радионавигационных систем [4].

Цель изобретения - повышение точности частотно-временной синхронизации мер времени и частоты, размещенных на подвижных объектах, на коротких интервалах времени.

Поставленная цель достигается тем, что на основе полученных аналитических соотношений, связывающих частотно-временные измерения при синхронизации бортовых мер времени и частоты в одной зоне радиовидимости НКА, разработан новый способ синхронизации этих мер.

Для большей наглядности основных выводов рассмотрим средневысотный НКА, движущийся по круговой орбите.

При этом примем следующие допущения [5]:

1. Земля является сферой радиуса R_з = 6371 км.

2. Движение спутника происходит под действием только силы притяжения Земли.

3. Ускорение силы земного притяжения всегда направлено к центру Земли, и его абсолютная величина g определяется по формуле Ньютона



где

R - расстояние от центра Земли до спутника;

- коэффициент, равный произведению гравитационной постоянной на массу Земли.

С точностью, достаточной для решения ряда практических задач, можно считать, что

=398610⁵ км³/с².

4. Величина и направление начальной скорости спутника (скорости вывода его на орбиту) подобраны из условия получения круговой орбиты.

При принятых выше допущениях, с учетом релятивистских эффектов, измеренная на неподвижном объекте разность частот f_изм между частотой f_н генератора неподвижного объекта и частотой f_б бортового генератора НКА определяется соотношением [3-6]

(3)

где

- релятивистское отношение скорости движения объекта V_c к скорости света C;

- угол между радиус-вектором и вектором скорости КА, измеренный в системе отсчета, связанной с подвижным объектом (фиг. 1).

Разложив выражение (3) в ряд по степеням и пренебрегая членами выше второй степени, получаем [7]

f_изм=f_н-f_б(1-cos-²/2+²cos²). (4)

Если ограничиться двумя первыми членами разложения (4), то получим приближенное выражение для измеренной на неподвижном объекте разности частот f_изм:

f_изм=f_н-f_б(1-cos). (5)

Будем считать, что частота наземного генератора f_н равна частоте бортового генератора f_б (f_н = f_б), тогда выражение (5) можно переписать в виде

f_изм=f_бcos . (6)

Соответствующая измеренная разность шкал времени наземного и бортового генераторов, связанных с системой отсчета наземного генератора, с учетом (6) будет иметь вид [2]:

(7)

где U_н = U_н-б(t_н) - начальное (на момент времени t_н) значение измеренной разности шкал времени.

Первая частная производная функции U_н-б(t) по времени равна

(8)

Момент времени t_мин достижения U_н-б(t) своего минимального значения (dU_н-б(t)/dt = 0) будет при cos (t_мин) = 0, тогда

(9)

что соответствует моменту времени t₀(t_мин = t₀) прохождения КА непосредственно над радиоприемником наземного генератора (точка траверза).

Рассмотрим случай, когда частота наземного генератора не равна частоте бортового генератора и отсутствует систематическое изменение частоты (f_н-б=f_н-f_б=const), что соответствует принятой в КНС "ГЛОНАСС" модели функционирования средств частотно-временного обеспечения [2]. Тогда выражение (5) будет иметь вид

. (10)

Соответствующая измеренная разность шкал времени наземного и бортового генераторов, связанных с системой отсчета наземного генератора, с учетом (10) запишется как:

(11)

Первую частотную производную функции U_н-б(t) по времени из (11) запишем

(12)

Обозначив

получим

(13)

Соответствующий момент времени t_мин достижения U_н-б(t) своего минимального значения (dU_н-б(t)/dt = 0) будет иметь вид

(14)

Сопоставление выражений (9) и (14) показывает, что, если разность частот бортового и наземного генераторов не равна нулю (f_н-б0), т.е. частота бортового генератора не равна частоте наземного генератора (f_нf_б), то момент времени t_мин нахождения минимальной разности бортовой и наземной шкал времени U_н-б(t) не соответствует времени t₀ прохождения источника излучения (НКА) непосредственно над радиоприемником неподвижного объекта (точкой траверза). Причем разность моментов времени t_мин - t₀ несет информацию о разности частот бортового и наземного генераторов.

Действительно, накопленная разность бортовой и наземной шкал времени за определенный интервал времени T_н обусловлена разностью частот соответствующих генераторов за этот же интервал времени.

Так как измеренная разность частот бортового и наземного генераторов в системе отсчета наземного генератора в точке траверза не равна нулю, то и соответствующий момент времени t_мин минимальной разности бортовой и наземной шкал времени не соответствует моменту времени t₀(t_мин-t₀), как показано на фиг. 2 и 3.

Приведем выражение (14) к виду

(15)

Предварительный анализ выражения (15) показывает, что существует принципиальная возможность определения относительной разности частот бортового и наземного генераторов в одной зоне радиовидимости КА, если измерить скорость V_с движения КА и угол между плоскостью распространения радиоволн, излученных бортовым передатчиком и направлением на наземный генератор в момент t_мин, когда U_н-б(t), т.е. измеренная разность шкал времени бортового и наземного генераторов, минимальна.

Величина (t_мин) = arccos(af_н/(1-af_н)) при условии использования высокостабильного генератора af_н=110^-10 - 110^-12) на борту средневысотного КА (=210^-5) близка к 90^o, поэтому можно считать [8], что на участке траектории полета КА, соответствующему интервалу времени от момента t₀ до момента t_мин, движение является прямолинейным и равномерным (фиг. 4).

В этом случае выражение для cos(t) будет иметь вид

(16)

где

t₀ - момент времени прохождения движущегося объекта непосредственно над радиоприемником неподвижного объекта (точка траверза);

R₀ - высота орбиты КА.

Подставляя (16) в (11) и используя первую подставновку Эйлера для нахождения интеграла, получим

(17)

Введем обозначение a_р=V_с²/СR₀ и перепишем выражение (17) в виде

U_н-б(t)=U_н+af_н(t-t_н)+1/2a_р х (1-af_н)(t-t_н) х (t-2t₀+t_н). (18)

Первая частная производная функция U_н-б(t) по времени с учетом выражения (18) имеет вид

(19)

Соответствующий момент времени t_мин достижения U_н-б(t) своего минимального значения (dU_н-б(t)/dt=0) будет

(20)

В таблице и на фиг. 5 приведены значения разностей, полученные в соответствии с выражением (20) для различных типов бортовых мер частоты и времени.

Приведем выражение (20) к виду

(20)

Таким образом, относительная разность частот бортового и наземного генераторов на одном витке КА может быть определена в соответствии с (21).

Исходными данными для расчета являются разность моментов времени прохождения движущегося объекта непосредственно над радиоприемником неподвижного объекта (точкой траверза) t₀ и момент времени t_мин нахождения минимальной разности бортовой и наземной шкал времени, а также параметры орбиты НКА - орбитальная скорость движения V_с и высота орбиты R₀ известные априорно или полученные в процессе сеанса измерений. При этом следует отметить, что разность моментов времени t_мин и t₀ составляет для КНС "ГЛОНАСС" величину от 1 с до 50 мкс, в зависимости от типа бортового стандарта частоты (см. таблицу и фиг. 5), что может быть измерено современными техническими средствами, применяемыми в наземном комплексе управления КНС.

Рассмотрим возможности технической реализации предлагаемого способа. Материальным объектом, на котором может быть реализован предлагаемый способ, является совмещенная радиотехническая система (РТС), в которой имеется запросный и беззапросный фазовые каналы измерения дальности (например, СУИК - система управления измерения и контроля КНС "ГЛОНАСС"). При этом момент времени t_мин нахождения минимальной разности наземной и бортовой шкал времени может быть определен в беззапросном канале измерения дальности ОТС, а момент времени t₀ нахождения точки траверза в запросном канале измерения дальности РТС (в этом случае получаемые измерения дальности не зависят от частоты бортового генератора). Остальные исходные данные, необходимые для расчета относительной разности частот бортового и наземного генераторов в соответствии с выражением (21), а именно орбитальная скорость V_c движения КА и высота орбиты R₀ измеряются известными методами [8] радиотехническими средствами наземного комплекса управления (НКУ) или передаются потребителю в составе эфемеридной информации.

Порядок определения разности моментов времени t_мин и t₀ следующий. Момент времени t₀ находится в запросном канале измерения дальности РТС. В этом случае на наземном пункте и на борту НКА устанавливается приемопередающая аппаратура в соответствии с тем, как показано на фиг. 6. Высокочастотные колебания частоты f₁, излучаемые передатчиком наземного пункта, принимаются на борту и после частотного преобразования излучаются бортовым передатчиком, работающем на частоте f₂. Такое преобразование частоты необходимо для частотного разделения прямого канала и канала ретрансляции. Измерение фазовых сдвигов проводится на наземном пункте) осуществляется на одной и той же частоте. Для этого частота ретрансляции f₂ выбирается такой, чтобы в месте сравнения (на наземном пункте) колебания частот f₁ и f₂ можно было бы легко путем преобразования привести к одной и той же частоте. Такое преобразование удобно выполнить, если частоты f₁ и f₂ удовлетворяют условиям



где

m и n - целые числа.

Ретранслированные сигналы на частоте f₂ принимаются наземным приемником, преобразуются к частоте f₁ и вместе с колебаниями передатчика подаются на измеритель разности фаз (измеритель интервалов времени), который регистрирует фазовый сдвиг между двумя этими колебаниями. После устранения фазовой неоднозначности (используя значение примерного, заданного априорно, расстояния от НКА до наземного пункта) на наземном пункте выполняются измерения в соответствии с выражением

U(t_н)=U(t_н+_ин)-U(t_н), (22)

где

Ut_н - измеренные значения разности фаз для времени t_н;

U(t_н+_ин ) - измеренные значения разности фаз для величины t_н+_ин;

_ин - период следования импульсов сигналов времени неподвижного объекта.

По результатам сравнения последовательных измерений значений U(t) определяется минимальное значение U(t_н) и соответствующий ему момент времени t₀ прохождения НКА через точку траверза.

Момент времени t_min находится в беззапросном канале РТС. В этом случае на наземном пункте и на борту НКА устанавливается аппаратура так, как это показано на фиг. 7. Колебания частоты f₃, излучаемые передатчиком НКА, наземного генератора подаются на измеритель интервалов времени (измеритель разности фаз). В соответствии с (22) производятся измерения разности фаз U_н-б(t) сигналов принятого с борта НКА и сформированного наземным генератором. По результатам сравнения последовательных измерений значений U_н-б(t) определяется минимальное значение U_н-б(t) и соответствующий ему момент времени нахождения минимальной разности бортовой и наземной шкал времени.

Следует отметить, что момент времени t₀ и t_мин в запросном и беззапросном каналах РТС определяются одновременно, в связи с этим частоты в этих каналах должны быть разнесены (f₁f₂f₃).

Проведенный анализ показывает практическую возможность определения моментов времени t₀ и t_мин с помощью технических средств современных РТС.

Оценим точность оперативного определения относительной разности частот бортового и наземного генераторов разработанным способом.

Поскольку значение относительной разности частот боритового и наземного генераторов определяется на основании полученной зависимости (21), связывающей ее с другими величинами, подвергаемыми прямым измерениям, то процесс определения этого значения относится к косвенным измерениям [9].

В общем случае можно считать, что требуется оценить истинное значение величины Q, которая связана с величинами Q_j (j = 1, 2,...m), измеряемыми прямыми способами, некоторым нелинейным уравнением

Q = F(Q₁, Q₂,...Q_m). (23)

В [9] показано, что в качестве наиболее достоверного значения косвенно измеряемой величины Q следует понимать значение, получаемое подстановкой в (23) косвенного измерения средних арифметических X_j рядов измерений исходных величин

(24)

Так как исходные величины в выражениях (21) измеряют с помощью различных средств измерения в разное время, то можно с большой степенью достоверности считать результаты измерений некоррелированными, поэтому дисперсия оценки (24) определяется из выражения

(25)

где

- дисперсия измерения средних арифметических рядов измерений исходных величин.

В соответствии с выражениями (24) и (25), средняя квадратичная погрешность определения разности частот бортового и наземного генераторов может быть оценена с помощью равенства

(26)

где

= (t_мин-t₀);

- средняя квадратичная погрешность измерения высоты орбиты КА;

- средняя квадратическая погрешность измерения орбитальной скорости движения КА;

S - средняя квадратическая погрешность измерения разности моментов времени t_мин и t₀ (t_мин-t₀).

Переменными в (26) являются:

1) относительная разность частот бортового и наземного генераторов;

2) средняя квадратическая погрешность измерения высоты орбиты КА;

3) средняя квадратическая погрешность измерения орбитальной скорости движения КА;

4) средняя квадратическая погрешность измерения разности моментов времени прохождения движущегося объекта непосредственно на радиоприемником неподвижного объекта (точкой траверза) (t₀) и момента времени нахождения минимальной разности бортовой и наземной шкал времени (t_мин).

Расчеты, выполненные в соответствии с (26) показали, что средняя квадратическая погрешность определения относительной разности частот бортового и наземного генераторов на ограниченном интервале (в зоне радиовидимости КА) с помощью разработанного способа, при относительной погрешности по частоте бортового генератора 210^-11, средней квадратической погрешности измерения высоты орбиты КА 1 м и средней квадратической погрешности измерения скорости КА 0,001 м/с составляет величину 1-210^-13, что примерно в 2-3 раза меньше погрешности оперативного определения разности частот бортового и наземного генераторов известным способом.

Таким образом, разработан новый способ оперативной синхронизации мер времени и частоты, размещенных на подвижных объектах, основанный на выявленной функциональной зависимости разности частот бортового (подвижного) и наземного (неподвижного) генераторов от разности моментов времени достижения минимального расхождения бортовой и наземной шкал времени и прохождения НКА над точкой траверза и параметров движения космического аппарата (движущегося объекта), позволяет повысить в 2-3 раза точность частотно-временной синхронизации генераторов, расположенных на подвижных объектах на ограниченных интервалах времени.

Источники информации

1. Шебшаевич В.С., Дмитриев П.П., Иванцевич Н.В. и др. Сетевые спутниковые радионавигационные системы / Под ред. П.П. Дмитриева и В.С. Шебшаевича. - М.: Радио и связь, 1982.

2. Тищенко В.А., Рассадин А.А. Центральный синхронизатор системы ГЛОНАСС//Радионавигация и время/РИРВ. - СПб. - 1994, - N 1. - с. 15-16.

3. Фатеев В.Ф. Введение в релятивистскую теорию радионавигации и синхронизации. Л.: ВИКИ им. А.Ф. Можайского, 1988.

4. Дейл С.А., Дейли П. Результаты наблюдения советских спутников глобальной орбитальной навигационной системы / Пер. с англ. IEEE PLANS, Las-Vegas.

5. Тихонравов М.К. Элементы теории полета и элементы проектирования искусственных спутников Земли. М.: Машиностроение, 1974.

6. Тоннелла М.-А. Основы электромагнетизма и теории относительности. - М.: Иностр. лит., 1962.

7. Корн Г. , Корн Т. Справочник по математике для научных работников и инженеров. - М.: Наука, Гл. ред. физ.-мат. литер., 1984.

8. Космические радиотехнические комплексы /Под ред. С.И. Бычкова. М.: Сов. радио, 1967.

9. Бурдун Г.Д., Марков Б.Н. Основы метрологии. - М.: Изд-во стандартов, 1975.