Текст произведения
(PDF):
Читать
Скачать
Введение
Оценка пространственных координат источников радиосигналов (ИРС) представляет большой исследовательский интерес, т. к. позволяет существенно повысить скорость передачи информации телекоммуникационных систем, улучшить точность устройств специального назначения и т. п. [1, 2]. Для этого происходит оценка углов сигналов по азимуту с последующим формированием диаграммы направленности, максимум которой направлен на полезный источник, а нули – в направлениях помех. Таким образом, отношение мощности сиг-нала к мощности помех и шума повышается, в результате чего, согласно теореме Шеннона, повышается скорость передачи информации. В основе такого рода конструкций на базовом уровне – многоканальные антенные решетки (АР), сигналы
с выходов которых подвергаются обработке, как правило, в цифровом вычислительном модуле. На заре пространственной цифровой обработки сигналов были популярны однородные линейные АР, основный недостаток которых заключается в возможности оценки координат только по азимуту
[3–9]. Из проведенных исследований, в частности [10], стало известно, что геометрия АР (т. е. координаты антенных элементов (АЭ) в декартовой системе координат) оказывает значительное влияние на точность оценок угловых координат ИРС. Необходимо упомянуть, что в последние несколько лет приобрели некоторую потребность формы АР, благодаря которым появилась возможность определения пространственных координат сигналов, как по азимуту, так и углу места, такие как кольцевые
и концентрические [11–14]. Таким образом, понятно, что на величину ошибок определяющее влияние оказывает набор АЭ в декартовой системе координат. И тогда, задавая значения антенн вдоль осей
X, Y и Z, можно значительно повысить точность оценок координат ИРС.
В работах [15, 16] рассматривается проблема получения нижней границы Крамера – Рао (ГКР) для оценки пространственных координат по азимуту
и углу места с использованием АР с всенаправленными элементами. Они используют полярные координаты элементов АР и не дают точного выражения, которое учитывало бы как положение источника, так и местоположение в пространстве каждого АЭ. Этот подход не позволяет точно оценить влияние размещения АЭ на характеристики данных оценок и построить решетку с наилучшими характери-стиками в определенных секторах сканирования.
Далее в работе описывается подход по снижению дисперсии оценок угловых координат радиосигналов за счет оптимального размещения АЭ
в пространстве. Предложенный подход основан на выражении, описывающем зависимость точности пеленгации от координат АЭ в декартовой системе координат, на основании ГКР. Вычисление новой формы АР осуществляется путем оптимизации по заданным критериям, т. е. положения антенн вдоль осей X и Y. Кроме того, полученное уравнение будет очень полезно для анализа важных факторов, определяющих точность радиопеленгации при использовании АР такого рода. Влияние этих факторов позволит создать такую конфигурацию, которая будет обладать лучшими характеристиками точности и разрешающей способности методов радиопеленгации.
Описание антенных решеток
Прежде чем перейти к изложению информации, касающейся подходов к оценке дисперсии при определении угловых координат радиосигналов, необходимо сделать предварительные замечания. Пусть имеется АР (рис. 1), которая состоит из
N ненаправленных элементов, произвольным образом распределенным на плоскости XY. Кроме того, допустим, что в раскрыв АР падает электромагнитная волна источника сигнала s(t) на несущей частоте ω0 с пространственными координатами θ (азимут)
и φ (угол места) относительно осей X и Z.
Рис. 1. Плоская антенная решетка произвольной формы
Fig. 1. Planar antenna array of arbitrary shape
Выражение узкополосной волны имеет вид [7]
где u(t) – амплитуда; ω0 – несущая частота; v(t) – фаза в зависимости от времени. Так как сигнал узкополосный, задержка τi вызывает сдвиг фазы :
где c – скорость света; i – индекс сигнала; j – мнимая единица. И теперь, если сигналы на АЭ описываются как x1, x2, …, xN, они выглядят в векторной записи следующим образом:
где а – направляющий вектор; λ – длина волны; – волновое число; – радиус-вектор, указывающий на n-й АЭ; Т – транспонирование.
Выражение дисперсии оценок радиопеленгации
Нижняя ГКР, связанная с концепцией информации Фишера, представляет собой теоретический предел средней точности оценки угловых координат радиосигналов для конкретной АР. При условии, что это объективная оценка вектора параметров, основанная на выборках, нижняя ГКР для произвольного числа сигналов и их параметров
(в первую очередь, определение азимута и угла места) выведена в работе [17]. Таким образом, матрица дисперсии ошибок оценок пространственных координат по азимуту и углу места может быть записана [18]
где – дисперсия шума; Tr – след матрицы;
– действительная часть; – поэлементное умножение; K – число отсчетов корреляционной матрицы; ; ; ; ; ; H – эрмитово сопряжение, S – корреляционная матрица сигналов;
A – матрица направляющих векторов;
R – корреляционная матрица сигналов и шума.
Покажем далее матрицы более подробно:
(1)
Далее необходимо раскрыть выражение (1) путем определения производных направляющих по азимуту и углу места:
здесь RT – матрица радиус-векторов антенн.
Для дальнейшего упрощения (1) на первом этапе определяем производную направляющего вектора по азимуту для k-го сигнала на n-м АЭ [12]:
(2)
где ψ – фаза сигнала.
На следующем этапе, подобно выражению (2), определим производную направляющего вектора относительно координаты по углу места для k-го сигнала на n-м АЭ:
(3)
После объединения выражений (2) и (3) получим дисперсию ошибок для одного сигнала:
(4)
где I – единичная матрица.
В более упрощенном виде формула (4) будет выглядеть следующим образом:
(5)
Оценка дисперсии для решеток с ненаправленными элементами
Подводя итог вышесказанному и на основании выражений (4) и (5), можно сделать вывод, что точность пеленгации по азимуту и углу места
в основном определяется координатами АЭ, а именно коэффициентом AR (5). В частности, для АР из двух элементов и единственного ИРС поступающего сигнала получается [18]:
(6)
Если имеется плоская АР, состоящая из трех элементов в плоскости XY и одного источника сигнала с произвольными координатами θ, φ, то [18]:
(7)
Согласно формулам (6), (7), точность радиопеленгации определяется суммой квадратов разностей координат антенн.
Два источника сигнала – три антенных
элемента
Далее рассмотрим сценарий, в котором есть
2 ИРС, а также 3 всенаправленных АЭ. Из линей-ной алгебры хорошо известно, что нельзя брать менее трех АЭ, поскольку в этом случае, как видно из выражения (4), обратная матрица внутри квадратных скобок становится сингулярной. Поэтому для совместной оценки координат двух источников сигнала необходимо иметь 3 или более антенн. Тогда можно записать:
Здесь матрица A имеет размерность 3 × 2, т. е.
3 антенны и 2 источника. Рассмотрим получение PN более детально:
Рассмотрим среднюю часть внутри квадратных скобок матрицы PN подробнее:
где коэффициент C равен
(8)
Чтобы получить уравнение (8) без мнимых
частей, применим следующее свойство:
После учета коэффициента C матрица PN
принимает следующие значения:
1.1
2.1
3.1
1.2
2.2 (9)
3.2
1.3
2.3
3.3 .
Элементы матрицы PN (9) занимают следующие позиции:
Принимаем следующие обозначения:
и
Для получения оценки нижней ГКР первой координаты сигнала необходимо выражение для матрицы PN производных D. В данном конкретном случае нас в первую очередь интересует следующее выражение:
Полученное выражение для AR для первого сигнала равно:
(10)
Предположим, что выражение коэффициента AR для второго сигнала равно
(11)
Тогда получается, что окончательное выражение для матрицы D1, которое в основном определяет точность оценки пеленгов, будет выглядеть так:
Обратная матрица D1 будет выглядеть следующим образом:
(12)
Поскольку выше было сказано, что нужны только значения первого и второго сигналов, то необходимо учитывать только элементы главной диагонали матрицы. Значения за пределами глав-ной диагонали отвечают за взаимное влияние между координатами двух сигналов, следовательно
и (13)
поэтому последние выражения после подстановки формул (10)–(12) в (13) будут определять точность оценок координат с помощью трехэлементной планарной АР в случае совместного пеленгационного измерения двух сигналов. Из двух последних формул можно сделать вывод, что для уменьшения погрешностей пеленгования как можно меньше, при условии, что известны приблизительные сек-торы обработки, AR11 и/или AR22 должны быть как можно большими:
и/или (14)
Рассмотрим ситуацию с двумя источниками сигналов, координаты которых имеют следующие значения: θ1 = 25°, φ1 = 45° и θ2 = 85°, φ2 = 45°.
В качестве эталона используем трехэлементную стандартную кольцевую АР с радиусом r = 0,0442 м. Далее, применив любой алгоритм оптимизации, будем одновременно максимизировать полученные выражения (14) относительно радиус-векторов r1, r2 и r3. Кроме того, предполагается, что приблизительные сигнальные секторы известны, а координаты x и y ограничены радиусом r. После оптимизации (14) получаем следующие координаты новой АР для приема и обработки сигналов из описанных выше зон θ1, φ1 и θ2, φ2: x1 = 0,0405, x2 = –0,0301,
x3 = –0,0246, y1 = –0,0182, y2 = 0,0328, y3 = –0,0369 (рис. 2).
Рис. 2. Схема стандартной кольцевой «○» и оптимизированной «●» решеток
Fig. 2. Scheme of standard ring “○” and optimized “●” arrays
Ниже приведены оценки ГКР полученной АР
в сравнении со стандартной традиционной кольцевой АР для различных ситуаций, т. е. для одного (рис. 3) и двух сигналов в зависимости от отношения сигнал/шум (ОСШ) (рис. 4), кроме того, приведен график разницы ошибок между оптимизированной и стандартной АР (рис. 5).
Рис. 3. Граница Крамера – Рао антенных решеток: сплошная кривая «○» – обычная решетка;
пунктирная кривая «◊» – оптимизированная решетка (см. рис. 2) с одним сигналом: θ1 = 0°–180°, φ1 = 45°
Fig. 3. Cramer-Rao boundary of antenna arrays: solid curve “○” - ordinary array;
dotted curve “◊” - optimized array (see Fig. 2) with one signal: θ1 = 0°–180°, φ1 = 45°
Рис. 4. Граница Крамера – Рао антенных решеток: сплошная кривая «○» – обычная решетка;
пунктирная кривая «◊» – оптимизированная решетка (см. рис. 2)
в зависимости от ОСШ; два сигнала: θ1 = 25°, φ1 = 45° и θ2 = 85°, φ2 = 45°
Fig. 4. Cramer-Rao boundary of antenna arrays: solid curve “○” - ordinary array;
dotted curve “◊” - optimized array (see Fig. 2) depending on SNR; two signals: θ1 = 25°, φ1 = 45° and θ2 = 85°, φ2 = 45°
Рис. 5. Разница между совместной оценкой границы Крамера – Рао двух координат сигналов
оптимизированной АР и стандартной кольцевой АР
Fig. 5. Difference between joint estimation of Cramer-Rao boundary of two signal coordinates
of an optimized antenna array and a standard ring antenna array
Два сигнала поступают в новую АР (см. рис 3), которая построена таким образом, что минимум ошибок приходится на целевой диапазон предполагаемых азимутальных местоположений сигналов,
т. е. несовместная оценка здесь выполняется. Кроме того, на графиках рис. 4, 5 проиллюстрировано, что сумма ошибок совместной оценки координат по азимуту через новую решетку немного ниже, чем при использовании стандартной кольцевой решетки.
Рассмотрим другую ситуацию с двумя источниками сигналов: θ1 = 50°, φ1 = 45° и θ2 = 140°,
φ2 = 45°. Также за эталонную принимаем трехэлементную стандартную кольцевую АР радиусом
r = 0,0442 м. После оптимизации (14) получаем следующие координаты новой АР: x1 = 0,0419,
x2 = –0,0079, x3 = –0,0042, y1 = 0,0142, y2 = 0,0373,
y3 = –0,0440 (рис. 6).
Рис. 6. Схема стандартной кольцевой «○» и оптимизированной «●» решеток
Fig. 6. Scheme of standard ring “○” and optimized “●” arrays
Ниже приведены рис. 7–9, иллюстрирующие результаты исследования случая с двумя источниками сигналов.
Рис. 7. Граница Крамера – Рао антенных решеток: сплошная кривая «○» – обычная решетка;
пунктирная кривая «◊» – оптимизированная решетка (см. рис. 6); один сигнал: θ1 = 0°–180°, φ1 = 45°
Fig. 7. Cramer-Rao boundary of antenna arrays: solid curve “○” - ordinary array;
dotted curve “◊” - optimized array (see Fig. 6); one signal: θ1 = 0°–180°, φ1 = 45°
Рис. 8. Граница Крамера – Рао антенных решеток: сплошная кривая «○» – обычная решетка;
пунктирная кривая «◊» – оптимизированная решетка (см. рис. 6) в зависимости
от ОСШ; два сигнала: θ1 = 50°, φ1 = 45° и θ2 = 140°, φ2 = 45°
Fig. 8. Cramer-Rao boundary of antenna arrays: solid curve “○” - ordinary array;
the dotted curve “◊” is the optimized array (see Fig. 6) depending on
the signal-to-noise ratio; two signals: θ1 = 50°, φ1 = 45° and θ2 = 140°, φ2 = 45°
Рис. 9. Разница между совместной оценкой границ Крамера – Рао двух сигналов координат
оптимизированной АР и стандартной кольцевой АР
Fig. 9. Difference between joint estimation of Cramer-Rao boundary of two coordinate signals
of an optimized antenna array and a standard ring antenna array
Из кривых на рис. 7 оценки точности пеленгации одного сигнала видно, что элементы новой АР расположены для оптимального приема и обработки сигналов, приходящих из секторов ≈ 20° и 140° по азимуту с помощью пеленгаторных методов. Общие ошибки радиопеленгации через новую решетку немного ниже по сравнению с обычной кольцевой решеткой после использования совместной оценки координат сигналов по азимуту, как показано на рис. 8. Кроме того, на рис. 9 показано, что разница между ошибками через новые АР меньше по сравнению со стандартными кольцевыми, поскольку графики в целевых секторах имеют отрицательные значения.
Таким образом, можно сделать вывод, что точность пеленгации может быть повышена только за счет оптимального расположения АЭ в пространстве без привлечения дополнительных вычислительных средств и методов.
Заключение
В работе рассмотрены планарные АР, состоящие из всенаправленных элементов и размещенные на плоскости XY. Основное назначение подобных решеток в современной связи – оценка координат пеленгации по азимуту. Однако в настоящее время очень актуальна проблема выбора наилучшей конфигурации решетки для радиопеленгации. Это уменьшает погрешности пеленга и количество АЭ и, таким образом, снижает вычислительные затраты. Нижняя ГКР является критерием оптимальности способности конкретной решетки оценивать координаты сигналов. Получены точные общие выражения ГКР для планарных АР. Формулы описывают зависимость точности пеленгования по азимуту или углу места от положения элементов антенны. Выявлено, что при наличии источника сигнала выражение в общем случае квадратичное. С другой стороны, две волны с произвольными координатами поступают на решетку, и зави-симость между точностью оценки пеленгации
и расположением АЭ выглядит как сумма косинусов разностей волновых чисел. Полученные выражения позволяют расположить АЭ так, чтобы ошибки пеленга были значительно уменьшены в секторе сканирования цели, превышая при этом измерения круглых АР.
Представленный подход использован для создания двух- и трехэлементных АР путем минимизации полученных точных выражений ГКР для оценки пеленгации по азимуту в случае одного источника и в случае двух одновременных волн. Было показано и доказано, что дисперсия оценок радиопеленгации снижается на 0,3–0,5° по сравнению с круговым массивом. Более того, рассматриваемый подход может быть расширен как на произвольное количество антенн, так и на угол места.