CC BY
174
УДК 621.396.969.11
В.А. Громов, Е.П. Ворошилин, М.В. Миронов
Оценка разности моментов прихода сигнала группировкой пространственно-разнесенных малых космических аппаратов
Исследован алгоритм оценки разности моментов прихода сигнала пространственно-разнесенными малыми космическими аппаратами. Проведено статистическое моделирование алгоритма оценки разности моментов прихода, предложены способы повышения точности оценки.
Ключевые слова: источник радиоизлучения, малый космический аппарат, погрешность определения координат, разность моментов прихода сигнала, моделирование.
В настоящие дни все больше разрабатываются и применяются активные радиолокационные средства радиомониторинга земной поверхности, размещаемые на космических аппаратах (КА). Однако активное однопозиционное или многопозиционное радиолокационное наблюдение с использованием космических аппаратов не позволяет решать специфические задачи обнаружения и определения координат источников радиоизлучения (ИРИ), которые могут быть замаскированы или по другим причинам недоступны для радиолокационного наблюдения.
Перспективным направлением космического радиомониторинга и контроля радиоэлектронной обстановки Земли является создание пассивной спутниковой системы радиомониторинга с применением космических аппаратов. Данная система позволит производить контроль с целью выявления несанкционированных источников излучения [1, 2], проверку соблюдения технических ограничений на параметры излучающей аппаратуры, поиск и местоопределение излучателей в случаях аварийных ситуаций, для обеспечения безопасности движения воздушных судов и т.п. Отдельной и важной областью применения подобной системы мониторинга радиоизлучений является радиоразведка в случаях военных конфликтов или для борьбы с терроризмом [3].
Космический сегмент подобной пассивной спутниковой системы радиомониторинга с применением МКА [4] может включать группировки низкоорбитальных МКА-ретрансля-торов (с высотой орбиты порядка 1500 км) для реализации разностно-дальномерного и разностно-доплеровского метода в распределенной многопозиционной системе определения координат ИРИ с высокой точностью. Также сегмент включает несколько КА на высокоэллиптических орбитах для выполнения общего обзора зоны мониторинга, решения задач обнаружения, выбора ИРИ для их дальнейшего анализа, выдачи целеуказаний группировкам МКА, обработки данных с МКА и передачи результатов в наземный центр обработки. В наземном центре решаются задачи определения координат, оптимальной фильтрации, идентификации и завязки траекторий. Координаты ИРИ вычисляются раз-ностно-дальномерным методом (РДМ) по измеренным разностям моментов прихода сигнала [5, 6]. На рис. 1 изображены примерные орбиты движения КА и диаграммы направленности их антенн. Благодаря совместной обработке информации о каждом объекте, получаемой в разнесенных точках пространства, достигается высокая точность определения координат.
Требуемая погрешность вычисления координат ИРИ в подобной системе находится в пределах нескольких сотен метров. Это накладывает следующие ограничения на точность измерения разности моментов прихода [4]:
- среднеквадратичная погрешность определения координат МКА-ретрансляторов не должна превышать 25 м;
- среднеквадратичная погрешность оценки времени задержки сигналов на МКА-ретрансляторы не должна превышать 30 нс.
Значения погрешности определения координат МКА-ретрансляторов могут быть обеспечены путем их высокоточной привязки с помощью космической радионавигационной системы ГЛОНАСС.
В статье рассматривается задача обеспечения требуемых значений среднеквадратичной погрешности оценки разности моментов прихода сигнала МКА-ретрансляторами.
Точность измерения разностей моментов прихода будет в основном зависеть от:
- отношения сигнал/шум принимаемых сигналов;
- синхронизации шкалы времени в каждом МКА;
- величины доплеровского смещения частоты сигнала, вызванного движением МКА.
Принимаемый МКА сигнал можно записать следующим образом [7]:
Ф) = «изл(0 + Л(t), (1)
где 8изл^) - излученный ИРИ сигнал; п^) - аддитивный гауссовский шум.
Такая модель сигнала является несколько упрощенной, но в то же время может быть использована при расчетах космических систем, потому что при распространении под большими углами к горизонту практически отсутствуют отраженные сигналы, а возмущения, вызванные влиянием канала распространения, можно ввести в п^) .
Определение разности моментов прихода Дтг сигналов на МКА-ретрансляторы в таком случае может быть вычислено путем их взаимно-корреляционной обработки [8]. При этом оценка Дт* поступления радиосигнала в две пространственно-разнесенные точки соответствует положению максимума модуля взаимно-корреляционной функции (ВКФ) сигналов, принятых в этих точках. Таким образом, для сигналов, зарегистрированных двумя пространственно-разнесенными МКА, можно записать:
Дт* = тах(|^ (т)|), (2)
1 Т
В (т) = T1+ т)^, (3)
T 0
где T - длительность записи сигнала; в^), - нормированные по средней мощности
сигналы ИРИ, принятые первым и вторым МКА.
Поскольку МКА перемещаются по орбитам с высокой скоростью (при высоте 1500 км радиальная скорость порядка 7 км/с), то будет наблюдаться доплеровское смещение частоты принимаемых сигналов.
Доплеровское смещение частоты, вносимое за счет движения МКА по орбите, при неподвижном ИРИ описывается выражением
Ь = Ь (1 (4)
где / - частота сигнала ИРИ; иг - скорость г-го МКА-ретранслятора; с - скорость распространения радиоволн в свободном пространстве.
Поскольку большинство современных устройств обработки сигналов являются цифровыми, то запишем формулу для вычисления взаимно-корреляционной функции дискретных сигналов, которая имеет вид
1 N-1
Е12 (/) = — Е х1 (к)Х2 (к + /), (5)
N 1=0
где N - количество отсчетов сигналов хг(й), Х2 (к); j - сдвиг одного сигнала относительно другого.
Расчет взаимно-корреляционной функции можно ускорить, используя теорему о корреляции, которая формулируется следующим образом [9]:
Д12 (У) =N рБг [х* (Щ Х2 (к)
(6)
где Х1 (Щ),Х2 (к) - быстрое преобразование Фурье (БПФ) сигналов х^к),х2 (к + j) ; -
обратное быстрое преобразование Фурье; N - количество отсчетов сигнала.
Данный подход требует выполнения двух прямых БПФ и одного обратного БПФ. Если количество отсчетов сигналов достаточно велико, метод быстрого вычисления корреляции позволяет получить результат с меньшими вычислительными затратами, чем непосредственный расчет по формуле (5).
Для исследования степени влияния различных факторов на точность оценки разности моментов прихода сигналов нами было проведено математическое моделирование в среде МА^АВ. Цель моделирования - получить статистические характеристики разностей моментов прихода для случая приема сигнала от ИРИ двумя пространственно-разнесенными МКА.
При моделировании частота сигнала ИРИ выбиралась в диапазоне от 100 МГц до 20 ГГц. Результаты, приведенные в статье, были получены при несущей частоте сигнала 1 ГГц. Моделируя супергетеродинный приемник, осуществлялся перенос сигнала на промежуточную частоту, равную 10 МГц. Затем производилась дискретизация сигнала с частотой взятия отсчетов, равной 40 МГц. ВКФ двух сигналов вычислялась по формуле (6), а разность моментов прихода сигналов оценивалась по формуле (7):
А *
ДХ; =
Д12 СО-N fd '
(7)
где Л12 (у) - взаимно-корреляционная функция сигналов; N - количество отсчетов сигналов; fd - частота дискретизации сигналов.
Истинное значение времени задержки сигнала между двумя МКА было задано равным 10 мкс. Огибающая сигнала сформирована с помощью оператора ШЩеуюЫ [10]. Алгоритм вычисления ШЩеуюЫ (Ь,г) следующий:
ю
(*) =
(
1 + cos
2-п(г -1)
L
-п
„ г - L „ т г - L
1, -+1 < г < L--
(
1 + cos
2П
г
2
(г -1)
L
-п
г - L л
0 < г <-+1 ;
2
г г - L
L--< г,
(8)
2
где L - длительность сигнала; г - коэффициент, характеризующий крутизну фронта импульса. Значению коэффициента г, равному нулю, соответствует прямоугольный импульс, равному единице - гауссовский импульс.
Модель сигнала на выходе приемника можно записать следующим образом:
s(k) = ю(к) - вт(2 - п - Дf - к - Дг) + п(к), (9)
где ю(к) - функция огибающей сигнала; Дf = ^ - fg + /д - частота сигнала на выходе приемника; ^ - несущая частота сигнала; fg - частота гетеродина; /д - доплеровское смещение частоты; п(к) - дельта-коррелированный гауссовский шум.
Доплеровский сдвиг частоты вводится в сигнал ИРИ с использованием формулы (4). Изменение скорости, приводящее к смещению частоты сигнала ИРИ, вводили в выражение (4) как случайный процесс с равномерным законом распределения, в пределах от -5 до 5 км/с. Если МКА движется по орбите по направлению в сторону ИРИ, то скорость МКА принимает положительное значение, в другом случае - отрицательное.
В работе также был исследован сигнал с линейной частотной модуляцией (ЛЧМ), частота которого задавалась следующим образом:
^ (к) = Дf + Ь - к, (10)
где Дf - начальная частота сигнала; Ь = (/1 -Дf)/к^; Д - конечная частота.
г
2
Результаты расчетов среднеквадратического отклонения (СКО) оценки разности моментов прихода сигналов приведены на рис. 2 и 3 (радиоимпульсный сигнал) и рис. 4 и 5 (ЛЧМ-сигнал).
Рис. 2. СКО оценки разности моментов Рис. 3. СКО оценки разности моментов прихода
прихода радиоимпульсного сигнала радиоимпульсного отгаата
(длительность 6,4 мкс) (длительность 12,8 мкс)
Расчеты показывают (см. рис. 2-5), что влияние доплеровского смещения частоты приводит к ухудшению точности оценки при больших отношениях сигнал/шум, где должна достигаться высокая точность оценки. Происходит это из-за смещения спектров сигналов (рис. 6), которое приводит к сглаживанию взаимно-корреляционной функции в районе максимума и уменьшению отношения сигнал/шум. В результате точность оценки разности моментов прихода уменьшается, а корреляционная обработка перестает быть оптимальной. Пример взаимно-корреляционной функции двух сигналов при доплеров-ском смещении частоты 16 кГц (при скорости МКА 5 км/с и несущей частоте сигнала 1 ГГц) приведен на рис. 7.
140—;.........................................:.........................................•• £ 180Г'зи»"."_"
\
----с доплеровским смещения частоты
- без доплеровского смещения
1,2,3- полоса си гнала 1, 2, 5 М Гц
51<;с-
----с доплеровским смещения частоты
- без доплеровского смещения
1,2,3-полоса сигнала 1, 2,5 МГц
Отношение сигнал/шум, дБ Рис. 4. СКО оценки разности моментов прихода ЛЧМ-сигнала (длительность сигнала 6,4 мкс)
Отношение сигнал/шум, дБ Рис. 5. СКО оценки разности моментов прихода ЛЧМ-сигнала (длительность сигнала 12,8 мкс)
Для повышения точности оценки разности Дт* между сигналами необходимо уменьшить влияние доплеровского смещения частоты. Устранить доплеровское смещение частоты одного сигнала относительно другого можно путем оценки величины смещения
спектра ДюД с последующим внесением этой поправки в один из сигналов:
5>(ю) = 81(ю)• 82(ю + Дюд -ДюД), (11)
где $1(ю),52(ю) - спектры сигналов; Дюд - доплеровское смещение частоты; Дюд -оценка доплеровского смещения частоты.
0.5
1 /: / : / ! .........: .....;..../... ч....... ч \ V ......V ----с
! : / / V \ \ \ Л \
.....4 А /// П / \ * \
/ / % \ \ \ ;\\ .......V: \ Д. \ \
- / /// \ \\ » л \ \ V
/ \ / \ : 4 \\ \л\ / ' \ \
R(т),
1 - Спектр первого сигнала
2 - Спектр второго сигнала (смещен на частоту Доплера)
3 - Результат умножения спектра 1 на спектр 2
4 - Идеальный случай умножения (без Доплера)
50
- 2
/ ........./.. 7 * V \
/ / / / / / Л......... \ \ \ \ \ \
// // и........... \\ \\
N
а
0. 0. 0. 0. 0. 0.
1 - Корреляционная функция сигналов с Доплером
2 - Корреляционная функция сигналов без Доплера
Рис. 6. Влияние эффекта Доплера на вычисление ВКФ
Рис. 7. Пример ВКФ двух сигналов при доплеровском смещении частоты
Т
Точность оценки частоты ДюД будет определять качество устранения влияния эффекта Доплера.
В литературе [11] активно применяют методы, основанные на вычислении центров тяжести спектров. Таким образом, алгоритм уменьшения влияния доплеровского смещения частоты на точность определения разности моментов прихода сигнала имеет вид:
1) вычисление центров тяжести спектров принятых сигналов
N-1 2 / N-1 2
юцт = I Щ IЫ , (12)
¡=0 / ¡=0
где 2. - отсчеты спектра сигнала;
2) вычисление величины оценки смещения по центрам тяжести спектров принятых сигналов
ДюД=юцт1 -юцт2; (13)
3) сдвиг спектра одного сигнала относительно второго на разностную частоту.
После этого производят вычисление оценки времени задержки между принятыми сигналами по формуле (7).
Для повышения точности оценки доплеровского смещения частоты необходимо применить прием интерполяции спектра сигнала путем добавления к массиву, содержащему спектр сигнала, нулевых отсчетов [12]. Их число выбиралось таким, чтобы длина обрабатываемого массива увеличилась в целое число раз п. На рис. 8 приведены результаты расчетов СКО оценки времени задержки для радиоимпульсного сигнала (длительность 25,6 мкс, отношение сигнал/шум 15 дБ) при увеличении разрешающей способности по частоте от 1 до 10 раз.
По результатам расчетов (см. рис. 8) следует, что увеличивать разрешающую способность по частоте более чем в 4-5 раз неэффективно, поскольку это не приносит дополнительного увеличения точности определения разности моментов прихода. На рис. 10 приведены результаты расчетов СКО оценки времени задержки при моделировании с четырехкратным увеличением количества отсчетов спектра для различных отношений сигнал/шум. По кривым, приведенным на рис. 9, видно, что при устранении эффекта Доплера при достаточно больших отношениях сигнал/шум (больше 20 дБ) точность оценки разности моментов прихода составляет 25-30 нс.
i - с доплеровским смещенем частоты --при устранении доплеровского смещения • • • без доплеровского смещения
\ , \
\ .......\
> —-
1
Коэффициент п
Рис. 8. СКО оценки разности моментов прихода радиоимпульсного сигнала при устранении эффекта Доплера
1С
и
а
Он
я и и
vV -с доплеровским смещением частоты --при устранении доплеровского смещения • • без доплеровского смещения
V \
\\
\ \
■. >
ч...... \
ч.....
— ;
10 12 14 16 18 20
Отношение сигнал/шум, дБ
Рис. 9. СКО оценки разности моментов прихода радиоимпульсного сигнала при устранении эффекта Доплера
По результатам моделирования можно сделать следующие выводы:
1. Поскольку сигналы в космических системах наблюдаются на фоне дельта коррелированного случайного процесса с нулевым средним и заданной дисперсией, то наилучшую точность оценки разности моментов прихода можно обеспечить, используя взаимно-корреляционную обработку сигналов, что приводит к необходимости размещения на борту МКА системы передачи сигналов.
2. Величина доплеровского смещения частоты в пассивных космических системах, состоящих из МКА на низких орбитах, может лежать в пределах от 0 до 330 кГц. Поэтому для обеспечения успешного устранения доплеровского смещения частоты (при сохранении точности измерения разности моментов прихода) необходимо повышать разрешающую способность спектра по частоте (т.е. увеличивать длительность реалиизации). Это приведет к необходимости увеличения объема передаваемых данных.
3. Уменьшить объем передаваемых данных с сохранением точности измерения допле-ровского смещения частоты возможно за счет использования интерполяции принимаемых сигналов.
4. Показано, что используя простой алгоритм нахождения центров тяжести спектров совместно с интерполяцией, за счет добавления нулей к реализации сигналов, можно практически полностью устранить влияние доплеровского смещения частоты, не увеличивая при этом объем передаваемых данных.
Таким образом, используя совокупность приемов, достигнута требуемая точность определения разности моментов прихода сигнала.
Статья написана в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг. (государственный контракт № 02.740.11.0183).
Литература
1. Pat. 6 417 799 US, Int. Cl. H 04 D 7/19. Method of locating an interfering transmitter for a satellite telecommunications system / A. Aubain (US), B. Lobert (US). -9.07.2002, - 5 p.
2. Pat. 5 594 452 US, Int. Cl. H 04 B 7/185. Method and System for Locating an Unknown Transmitter Using Calibrated Oscillator Phases / J.C. Webber (US), C. A. Knight (US). - 14.01.1997. - 12 p.
3. Pat. 6 839 017 US, Int. Cl. G 01 S 7/36. Method and apparatus for using anti-jam technology to determine the location of an electromagnetic radiation source / D.J. Dillman (US). - 4.01.2005. - 9 p.
4. Ворошилин Е.П. Определение координат источников радиоизлучения разностно-дальномерным методом с использованием группировки низкоорбитальных малых космических аппаратов / Е.П. Ворошилин, М.В. Миронов, В.А. Громов // Доклады Том. гос. ун-та систем управления и радиоэлектроники. - 2010. - № 1(21), ч. 2. - С. 23-28.
5. Черняк В.С. Многопозиционная радиолокация. - М.: Радио и связь, 1993. - 416 с.
6. Денисов В.П. Радиотехнические системы / В.П. Денисов, Б.П. Дудко. - Томск: Том. гос. ун-т систем управления и радиоэлектроники, 2006. - 253 с.
7. Левин Б.Р. Теоретические основы статистической радиотехники. - М.: Радио и связь, 1989. - 656 с.
8. Вознюк В.В. Определение времени задержки сигналов сверхмалыми космическими аппаратами / В.В. Вознюк, С.А. Зайцев, Д.А. Толстоусов // Известия вузов. Приборостроение. - 2005. - Т. 51, № 3. - С. 13-17.
9. Айфичер Э.С. Цифровая обработка сигналов: практический подход / Э.С. Айфичер, Б.У. Джервис. - М.: Вильямс, 2004. - 992 с.
10. Потемкин В.Г. Вычисления в среде МА^АВ. - М.: Диалог-МИФИ, 2004. -720 с.
11. Пальчик О.В. Обработка данных триангуляционного измерения / О.В. Пальчик, В.Г. Андреев // Цифровая обработка сигналов и ее применение: сб. докл. VI Междунар. конф. - М., 2004. - Т. 11. - С. 197-200.
12. Марпл С.Л. Цифровой спектральный анализ и его приложения. - М.: Мир, 1990. - 547 с.
Громов Вячеслав Александрович
Аспирант каф. радиотехнических систем ТУСУРа
Тел.: 8 (382-2) 41-36-70
Эл. почта: [email protected]
Ворошилин Евгений Павлович
Канд. техн. наук, зав. каф. телекоммуникаций и основ радиотехники ТУСУРа Тел.: 8 (382-2) 41-36-70 Эл. почта: [email protected]
Миронов Михаил Владимирович
Магистрант каф. радиотехнических систем ТУСУРа
Тел.: 8 (382-2) 41-36-70
Эл. почта: [email protected]
Gromov V.A., Voroshilin E.P., Mironov M.V.
Estimation of signal arrival time difference with the use of multi-position passive satellite system
In the paper, an algorithm of estimation of signal arrival time difference with the use of multi-position passive satellite system is investigated. Statistical modeling of the estimation algorithm is described. The ways for accuracy improvement of estimating the arrival time difference are suggested. Keywords: radio-frequency radiation source, small satellite, coordinates determination error, signal arrival time difference, simulation.
