CC BY
95
РАДИОТЕХНИКА
УДК 621.396.96
ВЫСОКОТОЧНЫЕ РАДИОМЕТЕОРНЫЕ СИСТЕМЫ СРАВНЕНИЯ ЭТАЛОНОВ ВРЕМЕНИ И ЧАСТОТЫ
шкалы эталона А. При этом выражения для шкал примут вид:
Ж А
АТ
8(1)
0
АТ
------»
АТ
0
*
Шкала А
Шкала В
Шкала С
АТ <--
а
«■
8(1-Т)
Т C A А АТ
АТ
Т
t
-►
t
0
Т
■>
t
Рис.1. Шкалы времени эталонов А, В, С
КАЩЕЕВ Б.Л., КОВАЛЬ Ю.А., КУНДЮКОВ С.Е.
Рассмотрены современное состояние радиометеорного метода сравнения эталонов времени и частоты, а также новые результаты исследований, полученные в ХТУРЭ в этом направлении. Описаны основные технические решения разработанных и внедренных аппаратурных комплексов. Намечены перспективы использования радиометеорного метода сравнения эталонов времени и частоты в Украине.
1. Задача сравнения эталонов времени и частоты
Синхронность или привязка территориально раз -несенных эталонов времени и частоты — основа для реализации современных частотно-временных методов, позволяющих решать актуальные научно-технические и оборонные задачи в системах координатно-временного обеспечения, метрологии времени и частоты, радиоастрономии, радиолокации, цифровой синхронной связи. Успешное решение таких задач определяется как стабильностью эталонов времени и частоты, так и точностью их сравнения. Высокие темпы совершенствования эталонов (относительная нестабильность современных эталонов составляет 10-13—10-14, а в перспективе — 10-15—10-16) определяют актуальность задачи повышения точности методов сравнения.
Основные методы сравнения реализуются в виде многопозиционных радиотехнических систем. По функциональному назначению, составу и принципам работы им присущи черты информационных, локационных и измерительных систем.
Идеальную шкалу времени аналитически можно записать с применением дельта-функций:
ад
HH) = X+ nT), (1)
п=-ад
где T — единичный интервал шкалы, размерность H(t) - с-1.
Начало отсчета времени для шкал территориально разнесенных эталонов, размещенных в одной инерциальной системе, выбирается произвольным, как единым для системы, так и индивидуальным для каждой из шкал. В последнем случае каждая из шкал описывается выражением, аналогичным (1). При выборе единого начала отсчета времени вводится сдвиг шкал времени DT, как алгебраическая величина, характеризующая взаимное опережение или запаздывание шкал. Например, на рис.1 изображены три шкалы времени, у которых единое начало отсчета времени принято совпадающим с началом отсчета
ад ад / \
HaН) = X8 + nT); HB(t) = X 8 + ( + АTBA),
п=-ад п=-ад
ад
He t)=X8 + nT + А TCA), (2)
п=-ад
где а:tca > 0 ; А TBA < 0 — сдвиги шкал эталонов C и B относительно эталона A. Соответственно можно ввести следующие сдвиги шкал а Tab =-А Tba ,
АTAC = -АTCA , АTCB = -АTBC .
Для вычисления сдвига шкал в случае
|АТ| < T / 2 можно предложить соотношение:
T/2 X
А TAB = j \{Ha +)- Hb(t))dtdx , (3)
-T/2-T/2
где HA(t) и HB(t) приведены в выражениях (2).
Интегрирование в формуле (3) выполняется по временной шкале, для которой определяется сдвиг (в данном случае по времени шкалы пунктов А или В). Для определения сдвига шкал в n-м периоде необходимо на nT изменить пределы интегрирования в выражении (3).
Результаты измерений сдвигов шкал эталонов системы можно представить квадратной матрицей АТ. По главной диагонали матрицы расположены результаты диагностики соответствующего пункта (в идеальном случае они равны нулю), а симметрично относительно главной диагонали — измеренные значения взаимных сдвигов шкал пунктов (теоретически должны отличаться лишь знаком).
Разлагая функции (1) и (2), описывающие шкалы времени, в ряд Фурье, можно показать связь шкал времени с эталонными частотами и дать фазовую трактовку сдвигов шкал:
ад
Ha (t) = 1/ T X exp(0 ),
k=-ад ад
Hb (t) = 1/T X exp[e A +Tf ), (4)
k = -ад ад
Hc (t) = 1/T X exp[[A +TT ),
k=-ад
где 0 A =kOt — полная фаза k-й гармоники спектра шкалы А; Q=2p/T — угловая частота первой гармо-
РИ, 1997, № 1
9
ники ^В =кОДТВА; =кОДТСА — начальные
фазы к-х гармоник спектров шкал В и С соответственно.
Согласно с выражением (4) получаются фазовые определения сдвигов шкал:
ДТВА = ( Q В - QA )/Ю; Д^А =( 0 C - QA )/Ю, (5)
где 0 В =kOt + кОДТВА, 0 В =kOt + кОДТСА — полные фазы к-х гармоник спектра шкал В и С соответственно.
Соотношения (5) справедливы и для случая различных частот эталонов. При этом сдвиг шкал ДТ является функцией времени. Анализ ДТД) позволяет оценить отличие частот эталонных генераторов. Так, если частота эталона В отличается от частоты эталона А на величину ДО<<О, то сдвиг шкал в n-м периоде
составит ДТВ =(кДОпТ + +кОДТВА)/кО = ДОпТ/ О + ДТВ , где ДТВ — начальный сдвиг шкал.
В рассмотренном случае ДО = О( ДТВ - ДТ$А ),
а в общем случае ДО = Od[AT(t)]/dt, где ДТД) — линейная аппроксимация дискретных значений Tn. Таким образом, сравнение шкал времени позволяет производить и сравнение (или синхронизацию) частот эталонов.
Первые методы сравнения были предложены А.Эйнштейном в период развития релятивистской теории (задача синхронизации часов). В одном из этих методов в синхронизируемых пунктах принимаются радиосигналы от общего равноудаленного излучателя, в другом — сигнал часов одного из пунктов принимается во втором пункте и ретранслируется обратно (по современной терминологии это дуплексный или ретрансляционный метод). В работах по частной теории относительности упоминается и вариант использования перевозимых часов.
К настоящему времени арсенал методов сравнения и варианты их технической реализации существенно расширились. В основе всех методов сравнения лежит передача информации о шкалах времени. В зависимости от алгоритмов передачи информации методы сравнения можно разделить на три группы.
К первой относятся методы односторонней передачи информации из одного пункта (A) в другой (B). Сдвиг БТ можно определить, если известно время задержки сигналов между пунктами. Примеры реализации методов первой группы — перевозимые квантовые часы (ПКЧ); СДВ, ДВ и спутниковые навигационные системы; система, использующая отражение сигналов от Луны; телевизионное сравнение; односторонний радиометеорный канал; кабельные, лазерные и волоконно-оптические линии. Поскольку в этих методах (кроме ПКЧ) излучаются сигналы только из одного пункта, методы первой группы принято называть пассивными.
Вторая группа методов основана на двусторонней передаче информации. При этом время задержки сигналов между пунктами может быть неизвестным, требуется только его высокая стабильность и обратимость. Методы второй группы принято называть активными. Они реализуются в радиометеорном, спутниковых и оптических каналах связи.
В третьей группе используется дополнительная высокостабильная шкала, сигналы которой принимаются в сверяемых пунктах. Сдвиг определяется с учетом известной разности задержек сигналов между пунктом дополнительной шкалы и сличаемыми пунктами. Методы этой группы реализованы в спутниковых радионавигационных системах (СРНС) Глонасс и GPS (Навстар) в дифференциальном режиме, а также в радиоинтерферометрах со сверхдлинными базами (РСДБ). В СРНС дополнительная шкала формируется эталоном на борту спутника. В РСДБ роль дополнительной шкалы играют высокостабильные сигналы космических источников (пульсары, квазары).
Основным критерием качества систем сравнения является погрешность измерения сдвига шкал. Погрешности измерений сдвигов шкал эталонов системы, подобно результатам измерений, можно представить в матричной форме. Сопоставление систем сравнения по точности возможно только при идентичных методиках расчета погрешностей (интервалы накопления результатов измерений, алгоритмы статистической обработки).
Общими для всех методов сравнения, за исключением ПКЧ, источниками погрешностей являются: нестабильность и невзаимность (для активных методов) задержки сигналов в канале; нестабильности аппаратурных задержек; внутренние шумы и помехи канала связи; динамика измеряемой величины из-за нестабильности эталонов. Каждый из методов имеет и специфические источники погрешностей, обусловленные особенностями канала, алгоритмами сравнения, аппаратурной реализацией, релятивистскими эффектами.
Опуская подробный сопоставительный анализ методов сравнения (рис.2, где имеют место сокращения: 1 — звуковое радиовещание; 2 — ДВ радиостанции; 3 — ИСЗ (низкоорбитальные); 4 — радионавигационные системы ДВ (“Чайка”, “Лоран-С”); 5 — ТВ; 6 - ПКЧ (НР-50601А, Ч1-76, 17н646); 7 -среднеорбитальные спутниковые навигационные системы ГЛОНАСС и GPS (7.1 — штатный режим, 7.2 —дифференциальный режим); 8 — радиометеорный метод (МЕТКА, КАМА, 17н91, 17н830) и связные ИСЗ; 9 — лазерные и волоконно-оптичесие линии), можно выделить три метода, обеспечивающие в настоящее время высшую точность и широко применяемые для решения научно-технических задач: ПКЧ, радиометеорный метод сравнения (РМС) и СРНС. Перспективным методом является РСДБ, поскольку он имеет не только высокую потенциальную точность, но в будущем может стать основой нового типа эталонов времени относительная нестабильность сигналов квазаров и пульсаров не хуже 1014 с.
2. Основы радиометеорного метода сравнения шкал времени и обобщенная модель метода
В основе РМС лежат высокая стабильность (оцениваемая скоростью изменения задержки — V) и обратимость задержки сигналов при распространении радиоволн метрового диапазона на дальность до 2000 км из-за рассеяния от ионизированных метеорных следов на высоте 80-105 км от Земли. РМС, наряду с методами, использующими ПКЧ и СРНС (GPS, ГЛОНАСС), в настоящее время обеспечивает
10
РИ, 1997, № 1
1 нс.
10 нс.
100 нс.
1 мкс.
10 мкс
0.1 мс.
1 мс.
10 мс.
100 мс.
1 с.
Погрешность
1 2 3 4 5 6 7.1 7.2 8 9
Рис.2. Современные средства и системы сличения шкал времени
погрешности сравнения порядка десятков наносекунд. При этом РМС превосходит указанные методы по таким показателям, как производительность измерений, автономность, оперативность, экономичность, скрытность, устойчивость к ионосферным возмущениям.
Работы по РМС были начаты в Харькове вскоре после первых публикаций на эту тему [1] и уже в начале 70-х годов разработаны первые экспериментальные образцы
ности задержки сигналов и отличные от идеальных характеристики; 3) конечное отношение сигнал/ помеха; 4) нестабильности сигналов эталонов и устройств привязки к этим сигналам; 5) непостоянство измеряемой величины из-за взаимного хода шкал сверяемых эталонов; 6) алгоритмы обмена измерительной информацией между пунктами, определяющие потенциальные погрешности методов (если реализация метода нуждается в отдельном канале связи, то ошибки при передаче информации — источник погрешностей РМС); 7) алгоритмы обработки результатов измерений.
Классификацию погрешностей РМС можно произвести по двум признакам — исходя из однотипных по характеру источников и по интервалам корреляции (см. рис.3.).
По характеру источников можно выделить погрешности: сигнальные, которые еще называют помеховыми (по причине появления) или потенциальными (из-за их неизбежности); трассовые, зависящие от характеристик МРРВ; аппаратурные; динамические, вызванные непостоянством измеряемой величины; методические. Особо в этом перечне следует отметить методические погрешности, которые зависят не только от вида метода, но и от
аппаратуры (Метка-1) для службы времени Госстандарта [2] и МО. При этом были получены погрешности сравнения 0,2...0,3 мкс [2]. Позднее к данному направлению подключился Казанский университет [3]. Головными организациями, представляющими интересы потребителей, явились: от Госстандарта — ВНИИФТРИ; от МО - ЛНИРТИ.
В дальнейшем работы проводились в направлениях: теоретические и экспериментальные исследования возможностей повышения точности РМС (исследования точностных характеристик метеорного радиоканала, анализ и синтез алгоритмов сравнения, поиск путей снижения сигнальных и аппаратурных погрешностей, совершенствование методик обработки результатов измерений) [4-7]; разработка, изготовление и внедрение образцов аппаратуры для сравнения эталонов времени и частоты Г осстандарта [8 -10]; разработка, изготовление и испытания макетов аппаратуры в рамках ОКР, проводимых в ЛНИРТИ (ныне РИРВ) [11].
Основными источниками и причинами возникновения погрешностей МС являются (рис. 3.): 1) амплитудно- и фазо-временные характеристики метеорного канала сравнения (МКС), а также некоторая нестабильность и невзаимность задержки сигналов при метеорном распространении радиоволн (МРРВ); 2) тракты аппаратуры, имеющие нестабиль-
Рис.3. Классификация погрешностей РМС
большинства иных источников погрешностей (см. рис.3.).
К основным методам сравнения принято относить встречный метод (ВМ), ретрансляционный и опорноретрансляционный методы (РМ и ОРМ). В первых образцах радиометеорной аппаратуры сравнения [2] применялся РМ; в одной из первых публикаций по РМС [1] описан ВМ, сущность которого иллюстрирует рис. 4. Более совершенным оказался ОРМ [12], временная диаграмма которого приведена на рис. 5. ОРМ запатентован и в США, причем материалы этого патента представлены после опубликования [12].
Ha(1)
0
дт
Ню
1
Шкала A ----->
А
<
tlA—Tp+ДТаВ
А
Шкала B ----->
t
0 tiB—Tp+ДТва
Рис.4. Встречный метод сличения
По интервалам корреляции погрешности РМС можно разделить на три группы. В первой группе интервал корреляции соизмерим с длительностью
РИ, 1997, № 1
11
4 HA(t) t«T«» 2АТаб Шкала A
0 АТ нБ© *• ,'^tlA ton \ t. гт; ;1 2ton<-rT t2A t 2АТба '-•-.jV-J Шкала B
1 АП Я' >
0 tlB ton І2Б 2ton t
Рис. 5. Метод опорной ретрансляции (n=1)
используемых сигналов; во второй — со временем существования метеорных следов; в третьей — превышает время ожидания метеорных связей.
Погрешности первой и второй групп являются случайными, учитывая их неизбежность и неустранимость. Погрешности третьей группы можно отнести к систематическим, разделив на безусловно систематические (постоянная величина и закономерно изменяющиеся функции времени) и условно систематические (медленно меняющиеся случайные функции). Грубые погрешности (промахи) в РМС возникают из-за аномальных измерений при малых отношениях сигнал/помеха, многолучевости при МРРВ и ошибок в канале связи.
Перечисленные источники погрешностей положены в основу обобщенной модели (ОМ) РМС. В ОМ РМС (рис. 6.) включены: методы (алгоритмы) сравнения; сличаемые эталоны с их параметрами; модель МКС; модель сигналов и помех; характеристики нестабильности и невзаимности аппаратурных задержек; алгоритмы обработки результатов измерений. В ОМ РМС не включены такие факторы, как
условия размещения аппаратуры в пунктах, выбор основных технических решений, качество реализации аппаратуры (показаны на рис. 6. пунктиром). Эти факторы целесообразно учитывать по экспериментальным результатам.
Величина единичного измерения (наблюдения) DT в соответствии с ОМ РМС может быть получена в результате решения системы уравнений, записанной в матричном виде:
и] = [А Т] + [тр] + [тс/п] + [тап], (6)
где [іи]; [АТ]; [тр]; [т/п]; [тап] — матрицы измеряемых интервалов времени, истинного сдвига шкал, задержек сигналов в РМКС, помеховых погрешностей, аппаратурных задержек соответственно.
Состав матриц уравнения (6) для основных методов (ВМ, РМ, ОРМ) приведен в табл. 1, где приняты
следующие обозначения: tf; tB; tf ; tB — измеренное время задержки принимаемых транслированного (ин-
декс 1) и ретранслированного (индекс 2) сигналов в пунктах А и В; DTAB — истинный сдвиг шкал пункта
А относительно пункта В; тПрд ; т^д; тПрм 'тПірм — величины задержек в передающих и приемных
Таблица 1
Ма- три- цы Методы
ВМ РМ (n < -l; ten = 0) ОРМ (n > -l; ton * 0)
[ти ] "tA" tb J rtA 1 t2 tb L L2 _
[ат] " АТаб " - АТаб " атАб " (n + 1>(Топ +АТАб )
[тр ] тр + АВЗ + АРВ L тр J т P + АВЗ + АРВ (l - n)lp + АВЗ + АРЛ + Х;т12 _
[тАП ] "в + т а " Т ПРД ^ТПРМ т А +тБ т ПРД ",_ТПРМ " тВ + т A " ТПРД +ТПРМ тБ + тА - п(тА +тБ ) 1прд ^ 1прм 11 у прд ^ 1прм )
[асп ] аА аБ ■ аА ■ АА2 - n • аЧ
трактах пунктов А и В соответственно; тр — обратимая часть времени задержки сигналов между пунктами при МРРВ; АВЗ; АРЛ — параметры, характеризующие невзаимность задержки сигналов из-за физических особенностей МРРВ и релятивистских эффектов соответственно;
(к . * BA AB \.
( А ВЗ + А РЛ = тр тр );
VT = dzP/dt — скорость изменения задержки сигналов в течение метеорного радиоэхо (характеризует нестабиль -ность задержки сигналов при МРРВ); т12 — временной интервал между циклами трансляции и ретрансляции сигналов (для
РМ и ОРМ); аА ; АВ; аА ; аВ — случайные абсолютные погрешности при измерениях временных положений сигналов для циклов трансляции и ретрансляции в пунктах А и В; n — коэффициент ретрансляции (для ОРМ и РМ); ton — время задержки опорного импульса (для ОРМ).
Решение уравнений (6) позволяет определить абсолютную суммарную погрешность единичного измерения в виде
А 5 =АТ-АТ =А ррв +А с/П +А АП + А РЛ , (7)
где АРРВ = АВЗ + АСТ — погрешности МРРВ; D^, Dcr — погрешности невзаимности и нестабильности задержки сигналов при МРРВ соответственно; А С / П —
Рис.6. Структура обобщенной модели РМС
12
РИ, 1997, № 1
помеховые погрешности; Д АП — аппаратурные погрешности.
Значения дт, а также слагаемые погрешности для основных методов сравнения приведены в табл. 2. Учитывая разнородный характер слагаемых погрешностей уравнения (7) и различия в интервалах корреляции, в ОМ РМС включены алгоритмы обработки результатов измерений.
Таблица 2
Результат измерения; абсолютные погрешности Методы сличения
ВМ РМ (n <_1; Ап = 0) ОРМ (n > _1; ton ^ 0)
Результат измерения - дтA tA _ tB L1 L1 tA(n_ 1)+tA _ Ап(п + 0.
2n ; (0.5tA _ Ап )_ дляп = 1
2
Погрешности МРРВ -дА ^РРВ Двз д ВЗ +д РЛ 2
Дст VtTi2; ± VtTi2 дляп = ±1 2n n
Помеховые погрешности _ ДА АСП дА _ дВ д2 _ пдЧ + A]1 (n _ l) 2n ; д2 _дВ з — L _ для n = 1 2
2
Аппаратурные _ ДА ПАП /Га тА \Св тВ Ч уПРМ ТПРД / уПРМ ТПРД/ 2
Предложенная ОМ РМС позволяет решать основные задачи анализа и синтеза рад иометеорных систем сравнения (РМСС) эталонов времени и частоты.
3. Результаты теоретических исследований РМС
При исследовании метеорного радиоканала основное внимание было обращено на параметры, определяющие потенциал РМС: полоса пропускания; нестабильность и невзаимность задержки сигналов. Экспериментальные исследования метеорного канала показали, что сигналы практически не искажаются, если их ширина спектра не превышает (10 ... 15) % от несущей частоты. Для оценки V использовались результаты измерений задержки сигналов по фазе несущих частот и по фазе разностной частоты с использованием разработанного комплекса ФАЗА, параметры которого приведены в табл.3. Выборочные средние значения для |Vt| составили: в локационном режиме — (2,5—3) 10-7; на трассе — (0,5—1) 10-7. Максимальные значения — (1— 2)10-6 и (5—7)10-7 соответственно. Закон распределения |Vtj| близок к экспоненциальному. В 5—10% случаев имели место существенные отличия V для различных частот и комплексов, что можно объяснить необратимостью задержки сигналов на метеорных следах.
Результаты экспериментального исследования кратковременной фазовой нестабильности задержки сигналов в метеорном радиоканале позволили обосновать выбор фазового принципа построения второго поколения промышленной радиометеорной
аппаратуры [11] и основных технических решений образца аппаратуры МЕТКА-11.
Таблица3
Характери- стики МЕТКА-5, 6, 6М, 7 "Фаза" Метка-11
Метод сравнения ОРМ (п=1) СРМ ОРМ (п=1) +ВМ ВМ
Вид сигнала Код (16 импульсов) ФМ импульсы + ДЧС ФМ импульсы + ДЧС
Ти, мкс 2 26 - ФМ, 500 - ДЧС 13*10 - ФМ 256*5 - ДЧС
2д&, МГц 1 1 0.25
Импульсная мощность, Ри ,кВт 20 - 40 10 2 - 5
Погрешность, нс 15 - 30 50 - 100 ФМ, 10 - ДЧС, 1 -ф О О
Использо- вание Службы времени Украины и России Прототип промышленной аппаратуры 17-Н-830 Госстандарт, пункты НАКУ КА, радиоинтерферометры
Анализ теоретических и экспериментальных исследований источников невзаимности задержки канала (поляризационные явления при распространении в ионосфере и отражении от метеорного следа, а также релятивистские эффекты) показал, что измерения с использованием современной метеорной аппаратуры сравнения пока не позволяют подтвердить модельные оценки невзаимности. При этом можно согласиться с верхней теоретической оценкой (менее 1 нс).
Результаты исследования метеорного радиоканала были положены в основу линейной параметрической модели канала сравнения. Линейная параметрическая модель метеорного канала сравнения (ЛПМ МКС) состоит из каскадного соединения линейного четырехполюсника с постоянными параметрами, идеальной линии с управляемой величиной задержки, идеального амплитудного модулятора и генераторов модулирующих функций (рис. 7). Представление метеорного канала сравнения линейной параметрической системой не противоречит физической теории метеоров, современным моделям МРРВ и одновременно позволяет использовать основные положения теории измерения параметров сигналов в рамках ОМ РМС.
Анализ прохождения сигналов в ЛПМ МКС показывает, что для периодических импульсных сигналов, временные параметры которых удовлетво -
ряют неравенству ти < тp < T < тм (ти — длительность импульса; т p — среднее время задержки сигналов между пунктами; T — период повторения импульсов; т м — средняя продолжительность метеорных связей), искажения формы импульса определяются только линейным четырехполюсником модели, а высота и задержка импульса зависят от его временного положения в пределах длительности метеорного следа. При этом выходной сигнал имеет вид суммы пачек идентичных по форме импульсов. Огибающая
РИ, 1997, № 1
13
пачек определяется амплитудно-временной характеристикой следа A(t). Длительность i-й пачки зависит
SBx(t)
Tp(t)
Г 1
SBbix(t) ------>
A(t)
N(t)
Li—lJ
N(t)
1
2
3
4
5
Рис. 7. Структурная схема линейной параметрической модели радиометеорного канала сличения:
1-линейный четырехполюсник; 2-идеальный фазовый модулятор (линия задержки); 3-идеаль-ный амплитудный модулятор; 4-генератор модулирующей функции tp(t); 5-генератор модулирующей функции 4(t); 6-генератор функции “численности” — N(t)
от вида (насыщенный или ненасыщенный) и параметров i-го метеорного следа. Следовательно, основной вклад МРРВ в погрешность сличения связан с вариациями амплитуд импульсов и изменениями задержек на интервалах, соизмеримых с периодом повторения и более.
Проведен комплексный анализ основных методов сравнения (ВМ, РМ и ОРМ) с учетом влияния конечного отношения сигнал/помеха, нестабильности задержки сигналов в канале и корреляции погрешностей в пунктах. Показаны преимущества ВМ и выполнен синтез более совершенных комбинированных и модифицированных методов — совместного ретрансляционного метода (СРМ), ВМ с обменом результатами измерений по метеорному каналу “ВМ + ОРМ”.
Для выбора вида и параметров сигналов для РМС решены теоретические задачи: анализ пороговых отношений сигнал/помеха; оценка потенциальных точностных характеристик основных типов простых и сложных сигналов; анализ влияния коррелированных помех. Анализ показал преимущества сложных сигналов и позволил наметить основные направления уменьшения сигнальных погрешностей: увеличение базы сигналов; усреднение результатов измерений по всем импульсам кода; переход к фазовым измерениям.
Анализ потенциальной точности для фазовых измерений проводился применительно к универсальному случаю когерентного дискретного составного частотного сигнала (КД СЧС), частным случаем которого является двухчастотный сигнал (ДЧС). Получены соотношения для потенциальной точности КДСЧС. Показано, что потенциальную точность КДСЧС можно реализовать при оптимальной обработке и измерениях временного положения каждого из элементов по огибающим и фазам с последующим весовым усреднением результатов этих измерений. Предложено техническое решение, реализующее такую обработку [13].
В настоящее время наибольший удельный вес в результирующей погрешности измерения сдвига шкал
в РМС имеют аппаратурные систематические погрешности (АСП).
АСП возникают вследствие нестабильности задержек в трактах формирования, излучения, приема и измерения времени прихода сигналов, а также невзаимности задержек в общих устройствах приема и передачи. Они не зависят от методов (алгоритмов) сравнения и составляют для результатов измерения сдвига шкал в пунктах:
где 4РМ , тПВрМ — задержки в приемо-измерительных каналах РМСС (АФУ и антенные коммутаторы в режиме приема, приемные устройства, измерители
временного положения сигналов); тПРд, тПРд — задержки в трактах формирования и излучения сигналов РМСС (устройства привязки к сигналам эталонов, формирователи сигналов, передающие устройства, антенные коммутаторы и АФУ в режиме
излучения); Аллп =[тПрм ~?Прд)/2 ,
Авап = (тТіРМ - тПщ )/2 — составляющие результирующих АСП, определяемые аппаратурными задержками в пунктах.
Из выражения (7) следует, что АСП становятся равными нулю, если комплекты аппаратуры в пунктах имеют идентичные параметры. Однако это невозможно обеспечить по причинам: 1) изменение питающих напряжений и температуры, старение и замена элементов приводят к нестабильностям задержек сигналов в трактах аппаратуры; 2) общие
устройства приема и передачи (АФУ и антенные коммутаторы) обладают невзаимными задержками сигналов; 3) вариация амплитуд принимаемых сигналов приводит к нелинейным искажениям и связанным с этим изменениям задержек.
На первом этапе развития РМС основное внимание уделялось первой из приведенных причин АСП. Для снижения АСП, вызванных нестабильностями задержек сигналов в трактах аппаратуры, предложены способы и устройства для контроля задержек сигналов в трактах (так называемое пилотирование задержек), а также способ компенсации АСП, получивший название “сквозной канал”. Известные способы пилотирования задержек и метод “сквозного канала” не позволяют контролировать и учитывать второй и третий источники АСП (невзаимности задержек сигналов в общих устройствах и изменения задержек сигналов при вариации амплитуд принимаемых сигналов). Более того, эти источники АСП экспериментально были обнаружены и исследованы сравнительно недавно [11].
В активных системах сравнения для исключения погрешностей, вызванных различными угловыми координатами траекторий прохождения сигналов между пунктами, используются общие АФУ для излучения и приема. При этом, на основании известного свойства обратимости антенн, задержки сигналов в АФУ при излучении и приеме полагаются одинаковыми, а следовательно, согласно выражению (7) — не влияющими на АСП.
14
РИ, 1997, № 1
Невзаимность задержек АФУ РМСС была впервые обнаружена в Харькове в 1991 г. при “локации друг на друга” двух полукомплектов аппаратуры МЕТКА-8. Изменение длин кабеля, смена АФУ и замена элементов в коммутаторе приводили к АСП до 20—30 нс.
Физически невзаимность задержек в АФУ можно объяснить изменениями формы сигнала из-за отличий в режимах согласования кабеля при приеме (антенна кабель ^ коммутатор ^ входное сопротивление приемника) и при излучении (передатчик ^ коммутатор ^ кабель ^ антенна). Поскольку идеальное согласование АФУ, что явилось бы наиболее простым решением этой задачи, невозможно, были продолжены экспериментальные исследования, а также проведен теоретический анализ искажений формы сигнала при рассогласованиях в АФУ и связанных с этим аномалий в задержках [10].
Отдельные результаты моделирования удалось подтвердить экспериментально. Для реальных длин кабеля и качества согласования АФУ РМСС аномальные задержки могут достигать единиц процентов от длительности импульсов, что составляет десятки наносекунд. Это может приводить к такого же порядка АСП. В наибольшей степени аномальные задержки проявляются в цепи передатчик ^ кабель ^ антенна, поскольку для достижения приемлемого КПД выходное сопротивление передатчика Zj>>ZB, а сопротивление антенны в полосе частот сигнала отличается от волнового не менее чем в 1,1—1,2 раза.
Для учета аномальных задержек в АФУ предложены и реализованы схемы пилотирования с использованием дополнительной измерительной антенны (ИА). Экспериментальные исследования схемы пилотирования с ИА, проведенные в режиме “локации друг на друга” двух полукомплектов аппаратуры МЕТКА-8 (полевая лаборатория ХИРЭ) и МЕТКА-6М — МЕТКА-5 (НПО “Метрология”), показали влияние на результаты пилотирования места размещения и высоты ИА. В случае с комплексами МЕТКА-6М — МЕТКА-5, высота подъема основных антенн которых составляла 2—2,5 м, пилотирование с ИА, размещенной на расстоянии 15—20 м от основных антенн, и оценка АСП по локационному режиму дали совпадающие результаты. У полукомплектов аппаратуры МЕТКА-8 антенны имели высоту 7 — 8 м. При перемещении ИА (на расстояние до 40 м) и изменении ее высоты отличия результатов пилотирования и локационных измерений составляли + 10 нс.
Таким образом, для повышения точности пилотирования АСП с применением ИА необходимо ее размещать на расстояниях, превышающих размеры несколькихзон Френеля.
Возможность измерять полную величину АСП аппаратуры в каждом из пунктов для сигналов с симметричной формой обеспечивает запатентованный способ зеркальной ретрансляции сигналов (ЗРС) [17]. При подключении аппаратуры и цифрового радиоимитатора (ЦРИ), осуществляющего ЗРС, к одному эталону, а также обратимости задержек ИА измеренная полукомплектом величина t составит:
0 = (ТПРМ - ТПРД) + 2tOn + АРИ , (8)
где тПРМ, тПРд — соответственно полные задержки трактов приема и передачи полукомплекта аппарату-
ры, включая АФУ; ton — момент времени, относительно которого ЦРИ производит ЗРС; А РИ — АСП ЦРИ.
Из выражения (8) следует, что АСП аппаратуры
ААП = (ТПРМ - ТПРД )/2 = А /2 - tOn - А-Ри/2
можно определить, если известны и стабильны параметры ЦРИ — tОП и АРИ .
Для оценки погрешностей в первых образцах аппаратуры использовались основные метрологические стандарты с учетом особенностей РМС. Однако не было принято единое понятие наблюдения, которое в одних случаях определялось как результат измерения в течение периода зондирования, в других — как усредненный результат измерений “на метеоре”. В алгоритмах обработки результатов не учитывалась их неравноточность и некоторые специфические особенности РМС.
В процессе разработки и испытаний в Харькове экспериментальных образцов аппаратуры предложены единые критерии оценки погрешностей и обобщенный показатель качества [6]:
в=ОІP2AfnBNk , (9)
где а — одна из разновидностей оценки случайной погрешности; р — излучаемая импульсная мощность; 2 А^ — ширина спектра сигнала; в — база сигнала; N — число импульсов в коде; k — число периодов, соответствующее одному измерению.
При подстановке в выражение (9) вместо а
энтропийной погрешности — АЭ либо вместо 2А^
и в соответственно — аЭФ и ВЭФ получается ряд вариантов обобщенного показателя качества РМСС.
Коэффициенты качества можно дополнить такими существенными показателями, как оперативность измерений (среднее число измерений в единицу времени либо время ожидания) и стоимость,
например, в виде Ps/C ( 5 — стоимость системы, С — численность измерений).
Предельные значения коэффициента в можно определить, используя результаты информационной теории измерений, где рассматривается близкий по физическому смыслу энергетический порог чувствительности измерительных устройств. С учетом особенностей метеорного радиоканала и временных
измерений emin И 10~2 (дж)1/2.
Подобные коэффициенту b показатели качества можно применить для сравнения не только различных типов РМСС, но и систем, основанных на других принципах. Достоинство предложенного показателя заключается в том, что с его помощью можно для определенных классов РМСС рассчитать диаграммы обмена энергии и (или) полосы частот на погрешность сличения.
Дальнейшим развитием обобщенного коэффициента (9) явилась обобщенная статистическая модель РМСС с обучением, позволяющая учесть условия размещения аппаратуры в пунктах, реальную помеховую обстановку, сезонные и суточные изменения численности, протяженность и ориентацию трассы, качество исполнения аппаратуры. В такой модели по
РИ, 1997, № 1
15
обучающей выборке, состоящей из результатов сличений, оцениваются статистические характеристики для числа сличений “на метеорном следе”, времени ожидания, погрешностей единичных измерений, а также выявляются параметры суточного (при соответствующей продолжительности обучающей выборки — сезонного) хода этих величин.
Обобщенная статистическая модель РМСС позволяет генерировать результаты сличений для измененных параметров трассы и аппаратуры (энергетика, вид и параметры сигнала, метод сравнения, техническая реализация отдельных блоков). Совер -шенствованию методов обработки результатов измерений посвящены работы [16,17].
Измерение сдвига шкал эталонов в течение сеанса, продолжающегося обычно несколько часов, дает массив результатов, который можно записать в матричной форме:
[ti, AT ij], (10)
где i = 1,2, ...., m — номер метеорного следа; m — количество метеоров за сеанс; j = 1,2, ..., n номер измерения в пределах i-го метеорного следа; n — количество измерений на i-м метеоре; t — текущее время i-го метеорного следа. В матрице (10) отдельные значения AT ^ отсутствуют, поскольку n у различных метеорных следов в общем случае отличаются.
Массив AT ij представляет собой случайную выборку из генеральной совокупности, в которой закон распределения в первом приближении можно считать нормальным, так как погрешность результата измерений представляет собой сумму большого числа независимых случайных величин. В подавляющем большинстве случаев, после амплитудной коррекции и исключения грубых промахов, гипотеза о нормальности закона распределения реальных массивов из -мерений подтверждается. Поэтому для таких массивов в качестве оценки результата измерений может выступать среднее значение (либо среднее взвешенное для неравноточных измерений), а в качестве погрешности — СКО — s. Причинами неравноточно-сти результатов измерений являются: флуктуации амплитуд принимаемых сигналов как “на метеорном следе”, так и при переходе от следа к следу; изменение уровня помех; кратковременная нестабильность и необратимость задержек сигналов “на метеорном следе”; случайные вариации времени задержки в аппаратуре. Количественно неравноточность подтверждается при сопоставлении s для групп единичных измерений AT, отличающихся амплитудами [16]. Предложены алгоритмы весовой обработки, различающиеся способами оценок весовых коэффициентов и этапами обработки [18,19].
Если гипотеза о нормальности закона распределения не подтверждается, рекомендуется применять робастные методы обработки результатов [17].
4. Разработка радиометеорной аппаратуры и многопозиционных систем сравнения эталонов времени и частоты
К настоящему времени разработано несколько поколений радиометеорной аппаратуры сравнения шкал времени — экспериментальные образцы и специальные макеты типа МЕТКА, созданные в
Харькове; макеты аппаратуры типа КАМА, разработанные в Казани (КГУ); два поколения промышленной аппаратуры 17-Н-91 и 17-Н-830, сконструированные в С.-Петербурге (РИРВ) при участии ХТУ-РЭ.
Основными потребителями аппаратуры являются службы времени Госстандарта СССР (ныне Госстандарта Украины и Госстандарта России) и других ведомств. Радиометеорная аппаратура сравнения используется для синхронизации пунктов Наземного автоматизированного комплекса управления космическими аппаратами (НАКУ КА) на территории Украины и решения ряда научных задач.
В службах времени Госстандартов Украины и России применяются не промышленные, а разработанные в Харькове экспериментальные образцы аппаратуры типа МЕТКА-6, отличающиеся более высокой точностью и меньшими габаритами. Основные характеристики комплексов МЕТКА-6 приведены в табл.3. Отличительными особенностями комплексов МЕТКА-6 являются: 1) аналоговое суммирование двух соседних импульсов кода и измерение временного положения путем аппаратурного усреднения по всем парам импульсов кода; 2) применение пяти кодовых комбинаций для отличия излучаемых и принимаемых сигналов; стробирование принимаемых сигналов в соответствии с кодом ожидаемых для приема сигналов; 3) высокоточный дискриминатор для оценки временного положения по нулевому переходу разности прямого и задержанного сигналов; 4) наличие специального АРУ, позволившего расширить динамический диапазон приемного устройства до 50 дБ; 5) ручное пилотирование аппаратурных задержек и возможность их компенсации в пределах ±90 нс с дискретностью 10 нс; 6) внутренний цифровой имитатор для диагностики измерительной части аппаратуры.
Метрологические характеристики аппаратуры типа МЕТКА-6 оценивались: в локационных режимах “сам на себя” и “друг на друга”; по независимым измерениям с применением ПКЧ и системы “Глонасс”; по рассеянию результатов наблюдений от среднего значения в каждом сеансе и от сеанса к сеансу на интервалах до месяца; при сравнении между собой результатов измерений в обоих пунктах. Диапазон изменений случайной погрешности и АСП, оцениваемой с помощью ПКЧ, не превышал 20 нс.
В 1993 г. разработан и внедрен в Киеве (УкрЦСМ), НАКУ КА и Харькове (ГНПО “Метрология”) комплекс МЕТКА-6М, сопряженный с ЭВМ. В модернизированный вариант аппаратуры МЕТКА-6М включена ЭВМ типа IBM РС/АТ 386 для оперативной статистической обработки результатов измерений, их хранения и вывода на печать. Наиболее совершенными точностными показателями обладают образцы фазовой радиометеорной аппаратуры сравнения (ФАЗА и МЕТКА-11). Аппаратура ФАЗА была разработана и испытана в Харькове в период 1985-90 гг. в ходе выполнения ОКР совместно с РИРВ. В аппаратуре ФАЗА использованы: составной сигнал, содержащий код из 15 фазоманипулирован-ных (ФМ) по Баркеру импульсов и когерентный ДЧС; комбинированный метод сравнения — ОРМ для ФМ сигнала и ВМ для ДЧС; формирование сигналов и гетеродинов по принципу ФАПЧ по сигналам эталона; цифровой фильтр для оптималь-
16
РИ, 1997, № 1
ной обработки ФМ сигнала на промежуточной частоте; фазовые принципы измерения временного положения на разностной и несущих частотах ДЧС; измерение амплитуд принимаемых сигналов и определение периода повторения в пределах секунды; передача информации о результатах измерений по метеорному каналу с помощью относительной фазовой манипуляции несущей ФМ импульсов; автоматическое пилотирование аппаратурных задержек и их компенсация в заданных пределах.
Типичные гистограммы результатов измерений макета ФАЗА приведены на рис. 8. Погрешность измерений по ФМ сигналу («60 нс) позволяет устранить неоднозначность фазовых измерений по огибающей ДЧС, имеющего ширину спектра 1 МГц. В свою очередь погрешность измерений по огибающей ДЧС (< 10 нс) позволяет, при соответствующем накоплении результатов измерений и стабильности современных эталонов, устранить неоднозначность фазовых измерений по несущим ДЧС. Таким образом, результирующая погрешность измерений может составить менее 1нс.
Аппаратура ФАЗА послужила прототипом не только для промышленной аппаратуры 17-Н-830, но и последнего, разработанного в Харькове, экспериментального образца аппаратуры МЕТКА-11 (табл .3). Основным отличием комплекса МЕТКА-11 является наличие двух режимов — сравнение шкал времени и передача информации. Отличительные особенности
Фаза элементов ДЧС Огибающая ДЧС
Рис. 8. Г истограммы результатов измерений
аппаратуры МЕТКА-11: несущая частота в диапазоне 40—50 МГц; относительно небольшая ширина спектра 0,2 МГц; применение цифровой обработки и измерения временного положения сигналов; измерение полной АСП с помощью ИА и устройства с ЗРС. Указанные особенности РМКС МЕТКА-11 позволяют уверенно достичь погрешности измерений менее 1 нс, облегчают решение задачи электромагнитной совместимости комплекса с действующими системами метрового диапазона и расширяют возможности функционального применения аппаратуры.
В аппаратуре МЕТКА-11, как и в макете ФАЗА, для вхождения в связь и устранения неоднозначности фазовых измерений используются ФМ по коду Баркера импульсы, а для высокоточных фазовых измерений по огибающим и когерентным несущим —ДЧС. При этом по-иному реализован ДЧС — в виде бинарной по закону меандра ФМ несущей частоты.
В состав аппаратуры МЕТКА-11 включена ЭВМ IBM РС/АТ-486, которая выполняет следующие функции: формирование грубой части шкалы времени с дискретностью 1 с; управление работой устройств цифрового блока; вычисление точной части значения временного положения принятого сигнала
по его цифровым отсчетам; накопление и статистическая обработка результатов сравнения, приема и передачи информации; помехоустойчивое кодирование и декодирование информации; обмен информацией с оператором, включающий отображение данных на экране дисплея и ввод данных с клавиатуры. Основные технические решения аппаратуры МЕТКА- 11 защищены изобретениями и патентом [ 13-15]. С 1985 г. действует многопозиционная радиометеорная система сравнения (МРМСС) эталонов Украины и России [9]. Структура МРМСС приведена на рис .9. МРМСС оснащена разработанными в ХТУРЭ образцами аппаратуры типа Метка-6,-6М. Аттестация служб времени в 1988 г. подтвердила высокие метрологические характеристики МРМСС.
В настоящее время, в соответствии с Соглашением о сотрудничестве стран СНГ по обеспечению единого времени и частоты (Бишкек, 9.10.92 г.), проводятся регулярные сравнения на трассе Харьков — Москва. Регулярно проводятся сеансы сравнения первичного (Харьков) и вторичных (Киев, Ужгород) эталонов Украины. Проведены экспериментальные сравнения в “треугольнике” Харьков — Киев — Москва (октябрь 1993 г.) [10].
5. Перспективы РМС
Анализ современного состояния РМС показывает, что возможности метода далеко не исчерпаны. Это относится не только к точности сличения, но и к таким характеристикам, как помехозащищенность, скрытность, дополнительные функциональные воз-
Рис.9. Многопозиционная радиометеорная система сличения эталонов времени
можности [7]. Наиболее перспективным направлением уменьшения случайных погрешностей РМС является применение сложных сигналов. Наивысшую точность при этом обеспечат КДСЧС, которые позволят реализовать высокоточные фазовые методы измерений не только на разностных, но и на несущих частотах элементов. Из алгоритмов (методов) сравнения наиболее перспективен ВМ в сочетании с обменом результатами измерений между пунктами по радиометеорному каналу связи. Наиболее актуальной и сложной остается задача снижения АСП. Прогресс в этом направлении возможен при использовании разработанных алгоритмов [15].
В новых поколениях аппаратуры расширятся функции ЭВМ. Помимо обработки результатов измерений, ЭВМ позволит осуществлять диагностику и управление режимами работы аппаратуры, адаптацию к параметрам метеорного следа и помеховой
РИ, 1997, № 1
17
обстановке, расширение функциональных возможностей аппаратуры.
Реальными являются следующие дополнительные функциональные возможности РМС: измерения разности частот и кратковременных нестабильностей эталонов частоты; измерения сдвигов шкал в пределах прямой видимости; сопряжение со связными спутниковыми системами для проведения сличений; передача информации; определение скорости ветра в метеорной зоне. Наибольший интерес представляет разработка аппаратуры для сравнения и передачи информации.
Особую актуальность имеет РМС для Украины, не располагающей собственными СРНС и предприятиями по выпуску ПКЧ. Поэтому вопросы совершенствования РМС включены в космическую программу Украины и программу создания Государственной службы единого времени и эталонных частот.
Разработка новых поколений аппаратуры может не только повысить эффективность существующей МРМСС, но и расширить структуру системы за счет ведущих метрологических центров стран СНГ и Европы.
Литература: 1. Lattorre V., Jonson G. Time synchronisation techniques//”IEE INT. Conv. Rec.”.— 1964.— Part 6.— Р. 422-428. 2. Дудник Б. С., Кащеев Б.Л., Лейкин А.Я. Использование метеорного распространения радиоволн для привязки часов службы времени и частоты // Измерительная техника.— 1971.— №12.— С.38-42. 3. Сидоров В.В., Кардоник Г.С., Плеухов А.П. Высокоточная синхронизация разнесенных радиотехнических систем.— Казань: Изд-во КГУ.— 1978.—224 с. 4. Дудник Б.С., Кащеев Б.Л., Коваль Ю.А., Пушкин С.Б. и др. Снижение погрешностей сличения шкал хранителей времени с помощью метеорного распространения радиоволн // Измерительная техника.— 1981.— №7.— С.41-44. 5. Дудник Б.С., Кащеев Б.Л., Коваль Ю.А., Пушкин С.Б. и др. Высокоточное сличение шкал Государственного и вторичного эталонов времени // Измерительная техника.— 1982.— №1.-С. 30-32. 6. Дудник Б. С., Коваль Ю.А., Пушкин С.Б. Критерии оценки точностных показателей радиометеорных комплексов для синхронизации шкал времени // Измерительная техника.— 1986.— №1.— С. 16-18. 7. Кащеев Б.Л., Коваль Ю.А., Дудник Б. С., Сидоров В.В. Перспективы совершенствования радиометеорного метода синхронизации эталонов времени и частоты // Измерительная техника.— 1992.— № 10.— С.26-28. 8. Дудник Б. С., Кащеев Б.Л., Коваль Ю.А., Семенов С.Ф. и др. Новый комплекс аппаратуры для сличения эталонов времени и частоты по радиометеорному каналу // Измерительная техника.— 1986.— №4.— С.15-16. 9. Пушкин С.Б., Кащеев Б.Л., КовальЮ.А., Семенов С.Ф. и др. Высокоточная радиометеорная система синхронизации эталонов времени Госстандарта СССР // Измерительная техника.— 1989.— №4.— С.22-24. 10. Кащеев Б.Л., Дудник Б. С., Коваль Ю.А., и др. Исследование радиометеорного канала сравнения на трассах Киев — Харьков- Москва// Метрология времени и пространства. Тез. 5-го Российского симпозиума.— Москва: 1994.— С. 143. 11. Кащеев Б.Л., Дудник Б.С., Коваль Ю.А. Семенов С.Ф. и др. Радиометеорные комплексы синхронизации эталонов времени и частоты с
погрешностями менее 10 нс//Измерительная техника.-1992.- №12.- С.31-32. 12. А.С 712807СССР,MКПG04С 11/02. Способ сравнения шкал времени/ Б.С.Дудник, Б.Л. Кащеев, В.П. Моисеев.- Опубл. 7.06.80.- Бюл. №4.- 13. А.С. 1741096СССР, ME^HG 04С11/02. Устройство для сравнения эталонов времени / Б.Л. Кащеев, Ю.А. Коваль, Б.С. Дудник и др. Опубл. 15.06.92.- Бюл. №22.- 14. А. С. 1385118 СССР, MКПG04 С11/02. Устройство для сверки шкал времени / Б.Л. Кащеев, Ю.А. Коваль, Г.В. Нестеренко.- Опубл. 30.08.88.- Бюл. 12. 15. Патент Украины 20380А, G04С11/02. Способ обработки сигналов при сличении шкал времени и устройство для его реализации /Б.С.Дудник, Б.Л.Кащеев, С.Г.Кундю-ков, Ю.А.Коваль, Г.В.Нестеренко, Б.И.Макаренко, Ю.А.Леман, С.Ф. Семенов. 16. Коваль Ю.А., Бавыкина В.В., Попович А.В. Исследование возможностей применения весовой обработки измерительной информации при радиометеорной синхронизации шкал эталонов времени // Радиотехника. Харьков: Изд-во ХГУ.- 1991.-Вып.94.- С. 38-44. 17. Бавыкина В.В., Коваль Ю.А., Попович А.В., Пугач А.А. Робастные методы обработки результатов измерения сдвига шкал времени // Методы представления и обработки случайных сигналов и полей. Тез. докл. IV Междунар. науч.-техн. конф.- Харьков: 1993.- С. 13.
Поступила в редколлегию 29.12.97
Кащеев Борис Леонидович, д-р техн. наук, профессор кафедры основ радиотехники ХТУРЭ, специалист в области радиометеорных исследований, инициатор и неизменный руководитель метеорных исследований в Харькове. Служебный адрес: 310726, Харьков, пр. Ленина, 14, тел. 43-17-58; Домашний адрес: 310024, Харьков, ул. Пушкинская, 74, кв. 14, тел. 43-33-11.
Коваль Юрий Александрович, д-р техн. наук, профессор кафедры основ радиотехники ХТУРЭ, специалист в области радиометеорного метода сличения шкал времени и обработки сложных сигналов. Домашний адрес: 310100, Харьков, пр. Ленина, 15, кв. 32, тел 40-95-55.
Кундюков Сергей Григорьевич, старший научный сотрудник ПНИЛ радиотехники ХТУРЭ, специалист в области радиометеорного метода сличения шкал времени и цифровой обработки сигналов. Домашний адрес: 310022, Харьков, ул. Сумская, 63, кв. 4, тел. 45-58-84.
18
РИ, 1997, № 1
