CC BY
164
УДК 621.396.96
А.С. Вершинин, Е.П. Ворошилин, В.П. Денисов
Зона действия многоканальных пассивных наземных систем радиотехнической разведки
Рассмотрены два метода многоканальной обработки сигналов, позволяющих повысить точность определения момента прихода. Проведено моделирование описанных методов. Рассчитана полная среднеквадратическая ошибка оценки момента прихода. Произведена оценка зоны действия пассивных систем радиотехнической разведки.
Ключевые слова: радиотехническая разведка, пассивная система, зона действия, распространение радиоволн, радиоимпульсный сигнал, многоканальная обработка, оценка момента прихода, результаты моделирования, задержка сигнала, порог обнаружения, определения координат.
Целью данной работы является оценка зоны действия пассивных систем радиотехнической разведки (РТР) при многоканальном приеме. Под зоной действия системы в данной работе понимается область пространства, в пределах которой среднеквадратическая ошибка оценки координат не превышает заданного порога. Для определения координат источника радиоизлучения в системах РТР обычно применяют угломерный и разностно-дальномерный методы либо результат их комплексиро -вания. Информативными параметрами в таких системах являются направление на источник сигнала и время задержки.
Ограничим круг рассматриваемых систем классом наземных пассивных стационарных систем радиомониторинга, использующих пеленгационную и разностно-временную информацию для оценки координат неподвижных источников радиоизлучения в УКВ-диапазоне. Особенностью функционирования подобных систем является отсутствие информации о принимаемых сигналах.
Регистрирующая аппаратура пеленгационных постов современных средств РТР, как правило, многоканальная с пространственным и поляризационным разнесением приемных каналов [1]. Сигналы с выходов этих каналов можно использовать для определения времени задержки с последующим определением координат источника радиоизлучения.
Влияние на точность пассивных наземных систем оказывают искажения принимаемого сигнала в канале распространения радиоволн (РРВ) и собственных шумов аппаратуры.
Повышение точности угломерно-разностно-дальномерных измерений и увеличение зоны действия системы возможны при совместной обработке информации от совокупности приемных каналов. Существуют две причины, позволяющие говорить о возможном увеличении точности местооп-ределения:
- увеличение отношения мощности принимаемого сигнала к мощности собственных шумов приемных устройств;
- формирование с помощью антенной решетки направленной диаграммы приема, ориентированной главным максимумом на источник излучения, и уменьшение, таким образом, влияния отраженных от местности и местных предметов радиоволн.
Модель канала РРВ. Под каналом распространения радиоволн понимается среда, в которой распространяется поле с момента излучения антенной системой источника до момента его регистрации антенной системой пункта наблюдения. При таком подходе характеристики канала существенно зависят от направленности и пространственной ориентации передающей и приемной антенн. Будем считать, что причиной искажений радиоволн являются отражения от дискретных переотраже-ний. Описание приземных каналов РРВ как линейных систем широко используется в литературе [2-6].
Используемой в работе моделью является линейный фильтровый канал с аддитивным шумом, который представляют линейной электрической цепью с постоянными параметрами [7, 8]. Для линейного канала входной и выходной сигналы связаны импульсной характеристикой h(t) или передаточной функцией K (ja ) [9, 10]:
ад
«вых() = 40*5вх () = I (-т) Т , (1)
-ад
где * - оператор свертки.
5вых(7о) = К(7о)5вх(7°) . (2)
Для канала с дискретными переотражателями и однократными отражениями импульсная характеристика представляется в виде:
т
Н() = ч05( - ¿0) + £ у ,5(/ - ¿0 - А/,), (3)
,=1
где ¿о, У 0 - время распространения прямого сигнала и его ослабление на трассе; А/,, у , - задержка /-го отраженного сигнала относительно прямого и его ослабление на трассе; т - число отражающих объектов.
Импульсной характеристике (3) соответствует передаточная функция в виде
N / ч
К(» = у0<■-10,0 +£у,е-/“,'°+А,'> . (4)
,=1
Модель сигнала. В качестве модели сигнала в работе используется простой радиоимпульсный сигнал:
5 (,) = Яе( А(,) • е] (®0/+ф0)), 0 <, <ти ,
где А(/) - огибающая; ®0 - несущая частота; ф0 - случайная фаза сигнала.
Методы многоканальной обработки. Для борьбы с замираниями сигнала и увеличения отношения сигнал/шум в современных системах радиомониторинга используют методы адаптивного управления диаграммой направленности (ДН) антенной решетки. Данные методы позволяют формировать ДН таким образом, чтобы максимум ДН был направлен на источник сигнала, а помеховые сигналы попадали в область минимума. Поскольку в данной статье рассматривается случай одного источника радиоизлучения, то можно ограничиться рассмотрением методов додетекторного (ДСС) и последетекторного суммирования сигналов (ПСС).
Графическая схема алгоритма додетекторного суммирования сигналов с выходов пространственно-разнесенных каналов приема представлена на рис. 1, где КР - квадратурный разделитель; МШУ - малошумящий усилитель; ФВ - фазовращатель; Д - детектор.
П. Л/ГТТГЧ/
Сумматор
и
Устройство определения момента прихода
Оценка момента прихода сигнала
Оценка пеленга
Расчет фазовых сдвигов
Сигнал от системы единого времени
Рис. 1. Графическая схема алгоритма совместной обработки сигналов с последующим определением оценки момента прихода
Амплитуда суммарного сигнала при ДСС на выходе антенной решетки в п-м временном дискрете:
r N-1 | (N-1
Uдсс(п) =Л Z Ai(n)cos(9¿ -ф/ф.в) + Z Ai(n)sin^¿-фгф.в)
. i=0 ) V i=0
V
N-1 2 N-1 2
Z ( (п)сюф,-ф.в + ASi (п)^пфг'ф.в) + Z ( (n)cosФгф.в - AC (фтф,-ф.в)
i=0 i=0
где Aci (n) , Asi (n) - отсчеты квадратурных составляющих сигнала в момент времени n ; N - число суммируемых каналов; фiфв - фазовый сдвиг.
Эффективность предлагаемого алгоритма зависит от ширины синтезированной диаграммы направленности, точности разностно-временных измерений, расположения на трассе переотражаю-щих объектов.
Возможным вариантом объединения сигналов с выходов пространственно-разнесенных каналов приема является схема последетекторного суммирования сигналов (ПСС), которая изображена на рис. 2. Достоинством схемы является то, что амплитуда суммарного сигнала мало зависит от направления прихода в пределах ширины диаграммы направленности каждого элемента решетки, что исключает необходимость поиска источника излучения по направлению в указанных пределах.
Сигнал от системы единого времени
Оценка момента прихода сигнала Рис. 2. Блок-схема устройства многоканальной обработки с последетекторным суммированием сигналов
Амплитуда суммарного сигнала на выходе схемы ПСС в п-м временном дискрете описывается выражением
N-11—2--------2---
иПСС (п) = Е V Ас/ (п) + Аб1 (п) ,
I=0
где АС1 (п), А$1 (п) - отсчеты квадратурных составляющих сигнала в момент времени п; N - число суммируемых каналов.
В энергетическом отношении схема ПСС подобна последетекторному накоплению пачки из N радиоимпульсов одноканальных радиолокационных станций (РЛС) [11]. В литературе этот процесс называется некогерентным интегрированием. Некогерентное интегрирование проигрывает до-детекторному суммированию с точки зрения увеличения отношения сигнал/шум, за счет подавления слабого сигнала шумом в детекторах. Потери зависят от отношения сигнал/шум на входах детекторов [12].
Метод оценки момента прихода. В современных пассивных системах местоопределения реализуются пороговые методы определения момента прихода сигнала. Данные методы характеризуются высокой оперативностью получения оценки момента прихода при малых вычислительных затратах.
Для дальнейших исследований выбран пороговый метод, реализуемый в аппаратуре измерителем с цифровым порогом, который определяет момент прихода по характерной точке переднего фронта огибающей импульса и основан на фиксации момента пересечения фронтом некоторого по-
рогового уровня. Момент пересечения порога является оценкой временного положения импульса. Алгоритм вычисления оценки момента прихода описывается условием
5 (/)> р П 5 (/ +1)> Р П-..П 5 (/ + йок )> Р, (5)
где 5(/), 5(/ +1), ..., 5(/ + йок) - отсчеты огибающей сигнала в дискретные моменты времени; р -пороговый уровень; П - логическая операция «И»; I - моменты времени; йок - длительность окна анализа.
При выполнении условия (5) за оценку момента прихода принимается значение I. Пороговый уровень р рассчитывается по шумам аппаратуры на выходе АЦП и обычно привязан к среднему квадратическому отклонению (СКО) шума. Моменты времени I и длительность окна анализа йок
измерялись числом тактов аналого-цифрового преобразователя.
Графическая иллюстрация алгоритма определения момента прихода устройством с цифровым порогом приведена на рис. 3.
4=5
Рис. 3. Иллюстрация алгоритма определения момента прихода импульсного сигнала устройством с цифровым порогом при йок = 5
cö
П
н
я
ч
G
§
сЗ
с*
cö
К
К
cö
PQ
О
Он
к
§
Он
о
я
1,2
1,0
0,8
0,6
0,4
0,2
1 Г Г 1 г ¡Г 4 і п і ■ 1 Г \ 1 Г “I
| Ог^бающа^ сигна^ : : без ш|ума kfytТ л1 1 |И ■U .У..-L-i
Реализация огибающей принятого сигнала /
- Порог 1 / 1 к
¡¡¡! г-J.i.LJ.J
* \ \ \ Момент прихода \ І і і і
Время, такты АЦП
Данный метод обеспечивает высокую точность определения момента прихода радиоимпульсно-го сигнала в условиях сильных искажений, вызванных прохождением приземного канала РРВ [13]. Это обусловлено тем, что использование низкого порога позволяет обнаружить «прямой» сигнал в точке, наиболее близкой к его началу, значит, наименее подверженной влиянию канала РРВ. Методы оценки момента прихода, использующие для оценки всю длительность сигнала (корреляционный, центра тяжести и т.д.), выдают смещенные оценки.
Параметры схемы обработки сигналов. В работе исследуются статистические характеристики оценок момента прихода следующих методов обработки сигналов:
1) среднее значение моментов прихода в одном канале (ОДИН);
2) последетекторное суммирование сигналов (ПСС);
3) додетекторное суммирование сигналов (ДСС);
При моделировании входная реализация сигнала принималась в виде суммы регулярного сигнала, рассеянной составляющей и внутреннего шума аппаратуры, причем для рассеянной составляющей была принята модель гауссова нестационарного (во времени) случайного процесса. Неста-ционарность возникает как следствие переходных процессов, имеющих место при прохождении импульсных сигналов через канал с отражениями. Задача генерации реализаций рассеянной составляющей сигнала была сведена к задаче генерации гауссова случайного процесса с заданной пространственно-временной корреляционной функцией:
К (т, х) = К (т)- К (Ах) ,
/ Л л /
где R(х) = 2<зр exp
Ах х0 у
exp
j ряу sin (а р)
комплексная огибающая функции корреляции по
пространственной координате [14]; R () = Ра р exp
т
т0.
- огибающая функции корреляции по вре-
мени; х - пространственная координата антенной решетки; X - длина волны; ар - угол прихода
относительно нормали к решетке; то - интервал временной корреляции; хо - интервал пространственной корреляции.
В литературе рассматривают два основных метода решения поставленной задачи: метод преобразования Фурье и метод формирующего фильтра [15].
Для использования, вследствие своей простоты и скорости, выбран алгоритм генерации случайного процесса методом дискретного преобразования Фурье, который может быть представлен последовательностью следующих шагов:
1) Задаются интервалы по пространственной координате х и времени Т, в пределах которых генератор обеспечивает выдачу случайного процесса с разрешением Ах и ДТ.
2) Пространственно-временная корреляционная функции Я (Ах,?) рассеянной составляющей
сигнала дискретизируется по пространственной координате и времени. В результате Я(х,?) можно
х Т
представить в матричном виде И (матрица И содержит — строк и столбцов).
Ах ДТ
3) Генерируется матрица N той же размерности, что и матрица И, элементами которой являются независимые случайные числа с нулевым математическим ожиданием, изменяющейся дисперсией по оси времени и фиксированной по оси пространственной координаты.
4) Вычисляется преобразование Фурье от матрицы И. Вычисление производится в два этапа: сначала по времени, затем по пространственной координате:
Ю ( Ю Л
(а,га) = | | %(х,?)е
—ю\—к>
5) Вычисляется преобразование Фурье от матрицы N
Ю ( Ю ^
• е—іахдх.
• е— іахдх.
6) Полученные спектральные представления корреляционной функции рассеянной составляющей сигнала (а,га) и спектр шума (а,га) перемножаются:
8п (а,га) = (а,га)-SN (а,га).
7) Вычисляется обратное преобразование Фурье от 5П (а,га), результатом которого является реализация пространственно-временного процесса ^(х,?) с заданной корреляционной функцией %(х,?).
В табл. 1 приведены параметры, значения которых изменялись при моделировании.
Таблица 1
Параметры, значения которых изменялись при моделировании
Обозначение Наименование параметра
?0 Интервал временной корреляции рассеянной составляющей сигнала
х0 Интервал пространственной корреляции рассеянной составляющей сигнала
X Длина волны
страс (? — ?пр ) Зависимость дисперсии рассеянной составляющей сигнала от времени относительно момента прихода прямого сигнала ?пр
°2б соб Дисперсия собственных шумов
1 Размеры антенной системы
N Количество элементов антенной системы
А = 1А Расстояние между элементами антенной системы
Арег (0 Огибающая регулярного сигнала
М Количество реализаций
Характеристикой, удобной для экспериментальной оценки увеличения точности определения момента прихода при многоканальном приеме, является полная среднеквадратическая ошибка (ПСКО) оценки момента прихода, которая вычисляется по формуле
1,-1
где Т, - оценка момента прихода в , -й реализации; ?ист - истинное значение момента прихода; N -количество реализаций; апско - полная среднеквадратическая ошибка оценки момента прихода. Результаты моделирования. Получены зависимости ПСКО оценки момента прихода от числа
приемных каналов, при двух различных комбинациях отношений
max
(Арег(t)) max (рег(t))
асоб
, где
а
рас
асоб - СКО собственных шумов аппаратуры, а рас - СКО рассеянных составляющих сигнала. Соответствующие результаты приведены на рис. 4, при параметрах приведенных в табл. 2.
о
X
к
«
к
CD
СТ
О
О
W
и
с
Один
— ПСС
\ —ДСС
\\ -
\ ч
\ \
а)
max
5 10
Количество каналов приема
(Арег (t)) _ 10 max (Арег(t)) _
10
асоб
5 10
Количество каналов приема
б)
maxi
(Арег(t)) = 25 max
(Арег (t)) _
2,5
Jpac
асоб
-•рас
5 10
Количество каналов приема
в)
max (
(Арег(t)) _
= 2,5,
max
(Арег(t)) _
10
асоб
Jpac
s
и
и
и
ст
о
О
U
С
16
14
12
10
- - Один
— ПСС
\ —1— ДСС
Ч\ ■
Ч'ч \ \
—>
г)
5 10
Количество каналов приема
max(Лрег(t)) _ 10 max (рег (t)) _
2,5
асоб
-фас
Рис. 4. Зависимость полной среднеквадратичной ошибки оценки момента прихода
от числа приемных каналов
Полученные результаты подтверждаются результатами экспериментальных исследований, приведенных в [13].
Из рис. 4 следует, что увеличение числа каналов приводит к увеличению точности оценки момента прихода для всех методов многоканальной обработки. Сравнение методов ПСС и ДСС пока-
зало, что метод ДСС обеспечивает более высокую точность по сравнению с ПСС и с увеличением числа каналов приема соотношение между точностями методов ПСС и ДСС изменяется слабо.
Таблица 2
Параметры модели, при которых получены представленные на рис. 4 зависимости____
Параметр to, такт АЦП Хо , отн.ед. W ^, отн.ед. А/, отн.ед. N, ед. Арег /арас , отн. ед. Арег /асоб , отн. ед. 5 ед.
Значение 4 4 10 1 измен измен. измен. 2000
Основываясь на данных об увеличении точности оценки момента прихода за счет использования многоканального приема, был проведен расчет зон действия трехпозиционной разностно-дальномерной системы (РДС). Пункты наблюдения РДС являются многоканальными пассивными РЛС с направленными антеннами. Пункты наблюдения располагаются в вершинах равностороннего треугольника со стороной 10 км. Антенные элементы пункта наблюдения разнесены на 4 интервала пространственной корреляции (х0). Зоны действия, соответствующие одноканальной и 16-канальной системам, приведены на рис. 5.
и
о
с
л
I -11
I -22 &
§ -33
-44
-55
-55
55
-33 -11 0 11 33
Координаты по х, км а) Пункты наблюдения - одноканальные
5 -33 -11 0 11 33 55
Координаты по х, км
б) Пункты наблюдения - 16-канальные
Рис. 5. Зоны действия РДС при различном числе каналов приема
Отношение зоны действия РДС при количестве антенных элементов в приемном пункте N к зоне действия системы с одним антенным элементом проиллюстрировано на рис. 6.
Зависимость, приведенная на рис. 6, получена в случае, когда искажения сигнала за счет влияния канала РРВ в пространственно-разнесенных антенных элементах сла-бокоррелированы (т.е. пространственное разнесение больше 50 длин волн). В случае уменьшения пространственного разнесения элементов системы скорость увеличения зоны действия РДС с ростом числа каналов приема будет падать из-за увеличения коэффициента корреляции искажений.
В результате работы можно сделать вывод, что мно-16 гоканальная обработка сигналов позволяет существенно повысить точность оценки момента прихода радиоим-пульсного сигнала и расширить зону действия системы радиотехнической разведки.
0
2 4 6 8 10 12 14
* Количество каналов приема
Рис. 6. Зависимость увеличения зоны действия РДС от числа приемных каналов
Статья написана в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 гг. (госконтракт № 02.740.11.0232).
Литература
1. Кондратьев В.С. Многопозиционные радиотехнические системы / В.С. Кондратьев,
А.Ф. Котов, Л.Н. Марков. - М.: Радио и связь, 1986. - 264 с.
2. Кеннеди Р. Каналы связи с замираниями и рассеянием. - М.: Сов. радио, 1973. - 302 с.
3. Кулемин Г.П. Рассеяние миллиметровых радиоволн поверхностью земли под малыми углами / Г.П. Кулемин, В.Б. Рассказовский. - Киев: Наукова думка, 1987. - 229 с.
4. Кириллов Н.Е. Помехоустойчивая передача сообщений по линейным каналам со случайно меняющимися параметрами. - М.: Связь, 1971. - 256 с.
5. Artes H. Linear time-varying channels / H. Artes, G. Matz, F. Hlawatsch. - Tech. rep., Institute of Communications and Radio-Frequency Engineering, Vienna University of Technology, 1999.
6. Bello P. Characterization of Randomly Time-Variant Linear Channels / P. Bello, I. Adcom, M.A. Cambridge // Communications, IEEE Transactions on [legacy, pre-1988]. 1963. Vol. 11, Р 360-393.
7. Прокис Дж. Цифровая связь: пер. с англ.; под ред. Д. Д. Кловского. - М.: Радио и связь, 2000. - 800 с.
8. Экспериментальное исследование структуры электромагнитного поля при распространении радиоволн сантиметрового диапазона над земной поверхностью / Г.С. Шарыгин, Ю.М. Полищук, Н.М. Лесков, В.Ф. Слюсарчук. - Томск: Изд-во ТГУ, 1970. - 127 c.
9. Гоноровский И.С. Радиотехнические цепи и сигналы. - М.: Радио и связь, 1986. - 512 с.
10. Баскаков С.И. Радиотехнические цепи и сигналы. - М.: Высшая школа, 1983. - 536 с.
11. Амиантов И.Н. Избранные вопросы статистической теории связи. - М.: Сов. радио, 1971. - 416 с.
12. Сколник М.И. Справочник по радиолокации / М.И. Сколник, К.Н. Трофимов, Я.С. Ицхок: В 4 т. Т. 1. Основы радиолокации: справочное издание. - М.: Сов. радио, 1976. - 456 с.
13. Ворошилин Е.П.. Экспериментальная оценка увеличения точности измерения задержки сигнала на трассе распространения за счет многоканального приема / Е.П. Ворошилин.,
B.П. Денисов. Радиоэлектроника // Известия высших учебных заведений России. - СПб., 2006. -
C. 18-24.
14. Денисов В.П. Фазовые радиопеленгаторы / В.П. Денисов, Д.В. Дубинин. - Томск: Томский государственный университет систем управления и радиоэлектроники, 2001. - 251 с.
15. Монаков А.А. Основы математического моделирования радиотехнических систем: учеб. пособие. - СПб.: ГУАП, 2005. - 100 с.
Вершинин Александр Сергеевич
Аспирант каф. телекоммуникации и основ радиотехники ТУСУРа
Тел.: 8-913-804-94-64
Эл. почта: vershkoff@sibmail.com
Ворошилин Евгений Павлович
Канд. техн. наук, зав. каф. телекоммуникации и основ радиотехники ТУСУРа
Тел.: (382-2) 41-36-70
Эл. почта: nvi@sibmail.com
Денисов Вадим Прокопьевич
Д-р техн. наук, профессор каф. радиотехнических систем ТУСУРа Тел.: (382-2) 41-38-69 Эл. почта: dvp@ms.tusur.ru
Vershinin A.S., Voroshilin E.P., Denisov V.P.
The service area of passive ground-based radio monitoring multichannel system
Two methods of multichannel signal processing, which allow to increase the accuracy of time-of-arrival (TOA) determining, are considered. These two methods are modeled. The TOA total standard error is estimated. The service area for passive ground-based radio monitoring multichannel system is calculated.
Keywords: radio-monitoring, passive system, range, microwave propagation, pulse signal, multi-channel signal processing, time-of arrival, TOA, estimation, modeling results, detection threshold, location, RF emitter.
