CC BY
60
УДК 621.396.96
ПРИМЕНЕНИЕ
РАДИОМЕТЕОРНОЙ АППАРАТУРЫ ДЛЯ ОПЕРАТИВНОЙ СИНХРОНИЗАЦИИ СТАНДАРТОВ ВРЕМЕНИ И ЧАСТОТЫ В ПРЕДЕЛАХ ПРЯМОЙ ВИДИМОСТИ
КОВАЛЬ Ю.А., КУНДЮКОВ С.Г.,
БАВЫКИНА В.В, ТКАЧУКА.А.,
ЛУЦЕНКО А.Б, НЕСТЕРЕНКО Г.В.
Показывается возможность применения сигналов радиометеорной аппаратуры «МЕТКА-6» для оперативной синхронизации стандартов времени и частоты в зоне прямой видимости. Рассматриваются особенности структуры и принципов работы этой аппаратуры и имитаторов метеорного радиоканала. Приводятся экспериментальные результаты, полученные на трассе ХГНИИМ - ХНУРЭ.
Для сличения эталонов времени и частоты в службах Госстандартов Украины и России применяются разработанные в Харькове экспериментальные образцы радиометеорной аппаратуры типа «МЕТКА». В 1993 г. разработан и внедрен в Киеве (УкрЦСМ), НАКУ КА и Харькове (ХГНИИМ) комплекс «МЕТКА-6М». В эту аппаратуру включена ЭВМ типа IBM для оперативной статистической обработки результатов измерений, их хранения и вывода на печать [1, 2].
С 1985 г. действует многопозиционная радиометеорная система сравнения (МРМСС) эталонов Украины и России [1]. МРМСС оснащена разработанными в ХНУРЭ образцами аппаратуры типа «МЕТ -КА-6» и «МЕТКА-6М». Высокие метрологические показатели аппаратуры «МЕТКА-6» и «МЕТ-КА-6М» подтверждены при аттестации службы времени СССР в 1988 г. и службы времени Украины в 1995 г.
времени в пунктах A и B соответственно, а AT -сдвиг между ними. Аппаратура синхронизации в обоих пунктах излучает зондирующие сигналы, привязанные к собственным шкалам времени. При наличии метеорного следа, обеспечивающего распространение радиоволн между пунктами A и B, в обоих пунктах принимаются сигналы друг друга, задержанные на величину тр . При этом измеряется временное положение принятых сигналов tA и tB :
tA — ДТ + тB +тBA +тА Ч +тПРД +тp +ТПРМ ,
tB ——AT + тА +tAB +тВ t1 - AL +т ПРД +Tp + ТПРМ ,
где гпрд и тпрд — аппаратурные задержки в устройствах формирования и излучения сигналов в пунктах A и B; т^м и тПРм - задержки в устройствах приема сигналов и измерения их временного положения.
HA(t)
< > < >
0
ДТ
HB(t)
X
^ t0 \
/ t А \
/ 1 n(to -tf)'
,(to -tB) \
/1А
ж-
N ✓
N ✓
Л
✓ N
/ ч
tB
t0
tB t
t2
0
>•
t
*■
Рис. 1. Опорно-ретрансляционный алгоритм
Свойство обратимости метеорного распространения радиоволн предполагает равенство времени распространения сигналов в обоих направлениях:
_AB _ НА _
ТР _ ТР _ ТР .
Далее в каждом пункте выполняется переизлучение сигналов в моменты времени, зеркально обращенные относительно фиксированных опорных меток
t0 : 2t0 - tA и 2t0 - tB (для n = 1). При приеме
До настоящего времени не рассматривалась очевидная возможность оперативной синхронизации стандартов времени и частоты, размещенных в пределах прямой видимости аппаратуры «МЕТКА-6М». Такое использование комплексов «МЕТКА-6М» позволит без существенных затрат расширить круг потребителей эталонной частотно-временной информации.
1. Краткая характеристика и принцип работы аппаратуры “МЕТКА — 6М”
В основе функционирования комплекса аппаратуры “ МЕТКА - 6М” лежит высокая кратковременная стабильность и обратимость задержки сигналов в метеорном канале связи. Для измерения сдвига шкал времени используется опорно-ретрансляционный алгоритм (ОРА) обмена сигналами в сличаемых пунктах (рис. 1). На(1) и Hb(1) - шкалы
переизлученных сигналов измеряется их временное положение:
1а = 2t0 _tB +ТПРД +тp +тITPM =
= 2t0 + АТ + (тПрм _ тАрд ) _ (тПрм + тПрд ) =
= ДТ + 2t0 +тдп ;
t2 - 2t0 _ t1 + ТПРД + тp + тПРМ -
= 2t0 - ДТ - (тПрм _ тПрд ) + (тПрм + тПрд ) =
= -ДТ + 2t0 -Тдп ,
где ТАП = (тАРм _ тПрД ) _ (тПрм + тПрД ) - аППаратурная составляющая, для оценки которой периодически выполняется измерение аппаратурных задержек (пилотирование).
14
РИ, 2004, № 1
В результате в обоих пунктах может быть получена оценка сдвига шкал:
ДТА = 2t0 - tf + тдп ; А1!6 = 2t0 - tB + тап •
Сверяемые пункты оснащаются идентичными полукомплектами аппаратуры “МЕТКА-6М”, каждый из которых содержит антенно-фидерное, передающее, приемно-измерительное устройство и ЭВМ. Основные параметры аппаратуры приведены в табл. 1.
Таблица 1
Алгоритм Вид сигнала Рабочая частота, МГц
ОРА, n=1 Код из 16-ти импульсов 57,3
Импульсная Ширина Длительность Частота
мощность, спектра, сигнала, повторения,
кВт МГц мкс Гц
20...40 1 16х2 50
В силу специфики метеорного распространения радиоволн необходимо, чтобы один из пунктов сличения постоянно изучал сигналы для поиска метеорного следа, пригодного для связи с корреспондентом (активный режим). При этом аппаратура второго пункта может работать как в активном, так и в пассивном режиме, когда сигналы излучаются только после обнаружения полезного метеорного следа.
Для формирования излучаемого радиосигнала используются опорные сигналы «5 МГц» и «1 с», поступающие от эталона. Радиосигнал представляет собой кодовую посылку из 8 пар импульсов, имеющих длительность 2 мкс. Длительность кодового интервала в паре импульсов равна 4 мкс, а кодовый интервал между парами составляет 16 мкс, за исключением одного интервала, который имеет длительность 32 мкс. Этот интервал может находиться в кодовой посылке радиоимпульсов на любом месте между второй и седьмой парой импульсов. Таким образом, база кода составляет 132 мкс и возможны пять различных кодовых комбинаций сигнала. Кодирование в системе необходимо для различения сигналов «свой — чужой» во избежание ложных измерений при локационных отражениях своих сигналов «назад». Следует отметить, что две пары кодов (№1 и 5; №2 и 4) являются зеркальными, а код №3 — симметричным во времени.
Период зондирования составляет 20 мс. В каждом периоде на приемник подаются стробы С1 и С2 для приема основного и переизлученного сигналов. Длительность и временное положение стробов фиксировано. Параметры строба С1 о беспечивают возможность приема зондов любого корреспондента на расстоянии до 2400 км и собственных сигналов в режиме локации. Временное положение и длительность строба С2 позволяют определять сдвиг шкал сличаемых эталонов в пределах ±200 мкс. При появлении полезного сигнала и его захвате
устройствами обнаружителя (декодером и схемами проверки по периоду) формируются узкие селектирующие стробы под принимаемые радиоимпульсы для измерителя интервалов времени.
После измерения временного положения сигнала Г , принятого под стробом С1, запускается устройство формирования момента переизлучения, а затем — формирователь радиосигнала для переизлучения. Результат измерения временного сдвига сигнала, принятого под стробом С2 относительно
опорного момента 2to , является искомой оценкой сдвига шкал. Последовательность результатов измерения сдвига шкал записывается в оперативную память и индицируется на цифровом табло, а также вводится в компьютер для статистической обработки и документирования измерительной информации.
Метрологические характеристики аппаратуры типа «МЕТКА-6» и «МЕТКА-6М» оценивались: в режимах «локация сам на себя» и «локация друг на друга»; по независимым измерениям с применением перевозимых квантовых часов (ПКЧ) и спутниковых систем; по рассеянию результатов наблюдений от среднего значения в каждом сеансе и от сеанса к сеансу на интервалах до года; при сравнении между собой результатов измерений в обоих пунктах. Диапазон изменений случайной и аппаратурной систематической погрешности (АСП), оцениваемый с помощью ПКЧ, не превышал 2030 нс.
2. Анализ возможностей имитаторов метеорного радиоканала для сличения шкал времени в пределах прямой видимости
Для диагностики радиометеорной аппаратуры сличения и оценки ее АСП разработано несколько типов аналоговых и цифровых имитаторов метеорного радиоканала (ИМР) [3]. В ИМР, предназначенных для измерения АСП, предложен и реализован алгоритм зеркальной ретрансляции сигналов (ЗРС) [4]. Этот алгоритм имеет сходство с ОРА при n = 1 и отличается от него тем, что ретранслированный сигнал зеркально обращен относительно опорного импульса (рис. 2). Это обстоятельство позволяет совместно использовать ОРА и ЗРС при установке зеркальной пары кодов.
tx
s(t)
S2 (t)
Шкала А
,* t
2t ОП
t s|(l)\ t r?'
1(2tОП _ t)
Шкала ИМР --------►
ОП
0
0
Рис. 2. Зеркальная ретрансляция сигналов
При подключении комплекса радиометеорной аппаратуры и ИМР с ЗРС к общему хранителю времени (рис. 3), а также при условии обратимости
РИ, 2004, № 1
15
задержки в АФУ ИМР и
s(t _ 2t ОП ~ т РИ
среде распространения между основной и измерительной антеннами (тр) сдвиг шкал, измеренный аппаратурой, работающей с опорно-ретрансляционным алгоритмом при n =
1, составит:
ДТ = (ТПРМ _тпрд)/2 + ТРИ /2 ,(1)
+ тПРД _тПРМ ) s(t _ ч s(t _ _ ч
s(t -т ПРД ч s(t -т ПРД -Т Р ч
Рис. 3. Измерение АСП с применением ИМР
где тПРМ , ТПРД — соответственно полные задержки трактов приема и передачи полукомплекта аппаратуры, включая АФУ;
три — задержка в ИМР. Отсюда следует, что можно определить АСП аппаратуры
(тпрм _тпрд)/2 = ДТ _ТРИ/2 ,
если известны и стабильны параметры ИМР (t0n
и ТРИ )-
Реализовать ЗРС можно, используя параметрическое (спиновое) эхо в диамагнитных или ферромагнитных средах. При воздействии на входы устройства, основанного на этом эффекте, сигнала Sj(t) и считывающего импульса в момент времени ton на выходе появляется эхо-сигнал s^^ = s^2t on _ 0 •
Можно реализовать алгоритм ЗРС на основе цифровой обработки сигналов. На рис. 4 приведена функциональная схема цифрового устройства ЗРС. Особенностью устройства является применение реверсивного счетчика адреса, который изменяет направление счета на обратное после воздействия
считывающего импульса в момент времени ton .
держки при подключении входа и выхода к одной измерительной антенне исключается из-за малого коэффициента передачи (амплитуда принимаемого сигнала составляет десятки вольт, переизлучаемого —десятки милливольт). Счетчик адреса непрерывно считает в одном направлении. После каждой смены адреса ОЗУ сначала выполняется считывание выборки, записанной при предыдущем обращении к этому адресу, а затем запись новых данных. Величина задержки радиосигнала в ИМР
составляет 13 = N / F + (tдС +1вс) и может регулироваться с дискретностью 1/F посредством изменения коэффициента пересчета счетчика адреса N. Максимальная величина задержки составляет 3,4 мс. Этого достаточно для работы через ИМР двух полукомплектов аппаратуры сличения, со сдвигом шкал, обеспечивающим отсутствие перекрытия сигналов во времени.
Режим задержки может быть использован для решения следующих задач:
— контроль работоспособности аппаратуры сличения;
Чтение из ОЗУ начинается с того адреса, на котором закончилась запись. Вся сохраненная с момента t3n последовательность цифровых отсчетов выдается в обратном порядке на ретрансляцию. Обратный счет имитирует реверсирование времени. Для работы устройства необходимо обеспечить кратность периода частоты дискретизации F и интервала ton _ tзп.
На рис. 5 приведена структурная схема цифрового ИМР с ЗРС, разработанного для аппаратуры «МЕТ-КА-6М». Данный ИМР может работать как в режиме ЗРС, так и в режиме задержки сигналов. В режиме задержки ИМР пе-реизлучает сигналы с управляемой задержкой и коэффициентом передачи.
Таким образом имитируется отражение от метеорного следа. В режиме задержки непрерывно происходит и прием, и переиз-лучение. Циркуляция сигнала в ИМР в режиме за-
— оценка случайной составляющей погрешности одного полукомплекта аппаратуры и ее зависимости от соотношения сигнал/помеха;
— оценка суммарной погрешности измерения сдвига шкал двумя полукомплектами при их размещении в одном пункте;
— решение вспомогательных задач по отладке аппаратуры, ремонту, обучению персонала.
Рис. 4. Цифровое устройство ЗРС
16
РИ, 2004, № 1
В режиме ЗРС имитатор через измерительную антенну принимает сигналы, излучаемые аппаратурой сличения, и переизлучает их, выполнив зеркальное обращение принятой реализации относительно момента времени tоп • Интервалы приема и переизлучения сигналов разделены во времени, что исключает циркуляцию сигналов в устройстве. Длительность этих интервалов определяется частотой дискретизации F и емкостью ОЗУ N: tз/С = N/F , а величина tоп зависит от 13/с , времени запуска tзп и разности задержек в устройствах дискретизации и восстановления t Вс сигнала: tоп = t3n +13/с _ (tдс _ tвс)/ 2 •
Используемая в формуле (1) величина трИ представляет собой АСП И МР и определяется отличием ton имитатора от опорного интервала t0 аппаратуры «МЕТКА-6М»: Три = ton — to •
Если шкалы времени аппаратуры «МЕТКА-6» и ИМР формируются от разных хранителей времени, то измеренный аппаратурой временной интервал с
точностью до АСП и три / 2 будет равен сдвигу шкал эталона, подключенного к аппаратуре «МЕТКА-6», и хранителя времени, подключенного к ИМР. Это обстоятельство и положено в основу предлагаемой методики синхронизации стандартов времени и частоты, расположенных в пределах прямой видимости от эталона, сопряженного с аппаратурой «МЕТКА-6М».
Суммарная величина АСП используемого полукомплекта аппаратуры «МЕТКА-6» и имитатора может быть непосредственно измерена при их подключении к общему хранителю времени. Затем ИМР размещается у потребителя, находящегося в зоне прямой радиовидимости аппаратуры «МЕТКА-6», и подключается к его хранителю времени. Благодаря этому достигается высокая точность предлагаемой методики измерения сдвига шкал.
3. Экспериментальная проверка канала синхронизации с применением ЗРС
Эксперименты по проверке предлагаемой методики проводились на трассе Харьковский государственный научно-исследовательский институт метрологии (ХЕНИИМ)—ХНУРЭ. Состав и структура экспериментального стенда приведены на рис. 6. В ХЕНИИМ комплекс «МЕТКА-6М», сопряженный с государственным эталоном времени и частоты, работал в штатном режиме с установкой «зеркальных» во времени кодов: код передачи №1; код приема № 5. Для излучения использовалась антенна («Ужгородская») на крыше корпуса, в котором установлен эталон, ориентированная на время эксперимента в направлении ХНУРЭ.
Для приема и переизлучения сигналов на крыше корпуса «И» ХНУРЭ установлена антенна в виде пятиэлементного волнового канала. Антенна проверена и ориентирована на здание ХЕНИИМ. В качестве ретранслятора использован ИМР с ЗРС, разработанный для аппаратуры «МЕТКА-6М» (см. рис. 5). В качестве стандарта использован рубидиевый генератор СЧВ-74. Контроль временного положения и оценка уровня радиосигналов проводились с помощью осциллографа С1-81.
Непосредственному проведению измерений предшествовал визуальный контроль уровня принимаемого сигнала и его ослабление до уровня, обеспечивающего отсутствие ограничения в ИМР. При этом также выполнялось грубое сведение шкалы сверяемого стандарта со шкалой эталона для того, чтобы ретранслируемые ИМР сигналы попадали во
второй строб приема вблизи момента 2ф . Затем в ИМР включался режим ЗРС и аппаратура «МЕТ-КА-6М» регистрировала сдвиг шкал между эталоном (ХЕНИИМ) и стандартом СЧВ-74 (ХНУРЭ).
Вх.
УВЧ 1
57,5 МГц
Вых.
УВЧ 2
1o МГц 1 с
От эталона
АЦП
ф Буф. регистр
ЦАП $=
АЛУ
60 МГц
Устройство управления
1 МГц
120 МГц
Синтезатор частоты
Буф. регистр
а' Ч Коэф. і к
передачи
Счетчик адреса
Рис. 5. ИМР для аппарпатуры «МЕТКА-6»
РИ, 2004, № 1
17
Рис. 6. Структура экспериментального стенда
Результаты эксперимента показали, что на данной трассе уровень принимаемых радиосигналов достаточный для работы ИМР без дополнительных входного и выходного усилителей. Для более протяженных трасс такие усилители могут потребоваться. В этом случае необходимо выполнить оценку уровня сигналов у потребителя и сформулировать требования к дополнительным усилителям.
Из-за особенностей программы обработки результатов в аппаратуре «МЕТКА-6М», ориентированной на прием сигналов через прерывистый метеор -ный радиоканал, единичные результаты получаются в виде пакетов размером до 254, между которыми имеются перерывы по 1020 с. Интервал между соседними результатами в пакете равен периоду зондирования и составляет 20 мс. В ходе экспериментального сеанса сличения, проведенного 16.09.2002 г., получено 1166 результатов единичных измерений (5 пакетов).
На рис. 7 приведена последовательность единичных результатов измерений сдвига шкал AT , входящих в пакет № 2, и линия регрессии, наклон которой характеризует погрешность установки частоты сличаемого хранителя (а), а также последовательность остатков (б). На рис. 8 приведена гистограмма распределения остатков линейной регрессии для того же пакета.
Результаты обработки всех 5 пакетов сведены в табл. 2. Разброс значений оценки AT/t по пакету не превышает 1-10-9 , а СКО остатков составляет приблизительно 15 нс, что согласуется с СКО единичных измерений аппаратуры «МЕТКА-6М».
Очевидно, что такая погрешность оценки AT/t по одному пакету слишком велика для оценки погрешности установки частоты рубидиевых стандартов, которая обычно имеет величину 10_10.. .10_12. Для этого необходимо использовать результаты измерения сдвига шкал, получаемые за более продолжительный период наблюдений.
Результаты эксперимента не позволяют оценить АСП измерения сдвига шкал, которая определяется нестабильностью аппаратурных задержек на интервале времени между измерениями АСП и сдвига шкал. Однако очевидно, что при использовании
предлагаемой методики величина АСП будет не больше, чем при применении аппаратуры «МЕТКА-6» в штатном режиме.
Рис. 7. Последовательность значений AT , линия регрессии (а) и остатки (б) для оного пакета
Рис. 8. Гистограмма распределения остатков
Таблица 2
Номер Число измерений AT/t • 10“8 СКО, нс
сеанса всего с ATj > 3ст
1 150 4 5,67 15,1
2 254 8 5,62 14,7
3 254 9 5,53 15,7
4 254 5 5,66 16,0
5 254 4 5,57 15,5
В ходе эксперимента было обнаружено, что результаты измерений, принадлежащие нескольким соседним пакетам, не удается объединить в одну последовательность в целях повышения точности оценки AT/t.
На рис. 9 показаны все результаты единичных измерений для пяти пакетов и линия регрессии,
18
РИ, 2004, № 1
построенная по всем 1166 точкам, наклон которой существенно отличается от наклона регрессии для одного пакета. Это обстоятельство можно объяснить тем, что время начала пакета фиксируется с дискретностью 1 с и с большой погрешностью (иногда более 10 с), обусловленной использованием системного времени компьютера и особенностями программы. Для того чтобы обеспечить сшивку результатов в непродолжительной последовательности пакетов, необходима доработка программного обеспечения ЭВМ в аппаратуре «МЕТКА-6М».
Рис. 9. Все значения AT для пяти пакетов
Для сличения частот высокостабильных генерато -ров можно использовать совместную обработку двух или более пакетов, полученных во время разных сеансов, с большим интервалом между сеансами, например, сутки, неделя или месяц. При этом относительная погрешность фиксации времени начала пакета окажется небольшой и не повлияет на оценку AT/t. Следовательно, такие измерения могут выполняться без изменения программы.
Выводы
Проведенный теоретический анализ и выполненные экспериментальные исследования показали, что имеется принципиальная возможность проводить высокоточные измерения сдвигов шкал и нестабильностей по частоте стандартов, размещенных в пределах прямой видимости от аппаратуры «МЕТ-КА-6М», излучающей эталонные сигналы. Для проведения измерений пригоден ИМР с ЗРС, разработанный для аппаратуры «МЕТКА-6М». Такие измерения оперативны (50 единичных измерений в секунду) и имеют случайную погрешность единичного измерения около 15 нс. Для конкретного места установки сверяемого стандарта необходимо оценить уровень принимаемых сигналов и включить в состав ИМР соответствующие усилители.
Для повышения точности измерений целесообразно разработать методику оценки АСП данного канала. Точность и оперативность предлагаемой методики существенно повысится при использовании фазового рад иометеорного метода синхронизации. Фазовые измерения позволяют исключить погрешность алгоритма ЗРС, обусловленную несимметричностью импульсных сигналов, которая неизбежно возникает при их прохождении через реальные тракты аппаратуры
Литература: 1. Кащеев Б.Л., Бондарь Б.Г., ГорбачВ.И., КовальЮ.А. Метеоры сегодня. К.: Техніка, 1996. 196 с. 2. Кащеев Б.Л., Коваль Ю.А., Кундюков С.Г. Высокоточные радиометеорные системы сравнения эталонов времени и частоты / / Радиоэлектроника и информатика. 1997. N 1. С.9-18. 3. КундюковСГ, Луценко А.Б. Измерение аппаратурных задержек в радиометеорных системах синхронизации // Радиотехника. 2002. Вып. 130. С. 121-126. 4. Патент України 20380 А, G0401/02. Спосіб обробки сигналів при звірянні шкал часу та пристрій для його реалізації / Б. С. Дудник, Б. Л .Кащеев Ю.О.Коваль, С.Г. Кундюков, Г.В. Нестеренко, Б.1.Макаренко, Ю.О. Леман, С.Ф. Семенов. № 95094289; Заявл. 26.09.95; Опубл. 15.07.97.
Поступила в редколегию 28.04.2003
Рецензент: д-р. техн. наук Величко А.Ф.
Коваль Юрий Александрович, д-р техн. наук, профессор кафедры основ радиотехники ХНУРЭ, специалист в области радиометеорного метода синхронизации шкал времени и обработки сложных сигналов. Адрес: Украина, 61166, Харьков, пр. Ленина, 15, кв. 32, тел. (0572)40-95-55.
Кундюков Сергей Григорьевич, канд. техн. наук, старший научный сотрудник Научно-учебного центра кафедры ОРТ ХНУРЭ. Научные интересы: синхронизация шкал времени, цифровая обработка сигналов. Адрес: Украина, 61022, Харьков, ул. Сумская, 63, кв. 4, тел: (0572)45-58-84.
Бавыкина Вероника Всеволодовна, канд. техн. наук, доцент кафедры ОРТ ХНУРЭ. Научные интересы: синхронизация шкал времени, цифровая обработка сигналов Адрес: Украина, 61118, Харьков, пр. Тракторостроителей, 118-А, кв. 28, тел. (0572)40-91-13.
Ткачук Александр Александрович, канд. техн. наук, старший научный сотрудник, зав. лабораторией ХГНИ-ИМ. Научные интересы: метрология времени и частоты, обработка сигналов. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. (0572) 40-30-96.
Луценко Александр Борисович, аспирант кафедры ОРТ ХНУРЭ. Научные интересы: синхронизация шкал времени, цифровая обработка сигналов. Увлечения: разработка прикладного программного обеспечения. Адрес: Украина, 61068, Харьков, ул. Тарасовская, 17, кв. 41, тел. (0572) 21-86-08.
Нестеренко Георгий Викторович, канд. техн. наук, старший научный сотрудник, директор научно-учебного центра кафедры ОРТ ХНУРЭ. Научные интересы: радиолокация, радиосвязь, медицинская электроника, автоматизация технологических процессов в животноводстве и птицеводстве, техническая диагностика. Увлечения: точная механика, приборостроение, специальное оборудование. Адрес: Украина, 61000, Харьков, ул. Лагерная, 71/1, кв. 13, тел. (0572)43-17-58.
РИ, 2004, № 1
19
