Алгоритм синтеза геометрической структуры антенной решетки MIMO РЛС

Доклады БГУИР

2012

№ 8

УДК 621.396.96

АЛГОРИТМ СИНТЕЗА ГЕОМЕТРИЧЕСКОЙ СТРУКТУРЫ АНТЕННОЙ

РЕШЕТКИ MIMO РЛС

В статье представлен алгоритм синтеза геометрической структуры антенной решетки MIMO (Multiple-input-multiple-output) РЛС, исходя из выполнения условия единственности главного лепестка диаграммы направленности. Исходными данными для алгоритма являются число передающих и приемных элементов, а также требуемое число формируемых каналов приема (виртуальных элементов).

Ключевые слова: MIMO (Multiple-input-multiple-output) радиолокационные системы, ортогональные сигналы, синтез геометрической структуры антенной решетки.

В последнее время в радиолокации активно развиваются MIMO (Multiple-input-multiple-output) РЛС [1,2]. В таких РЛС K различных групп, передающих элементов антенны, излучают K ортогональных сигналов, а L групп приемных элементов обеспечивают одновременный прием этих сигналов.

Одним из важных свойств MIMO РЛС является увеличением числа приемных каналов при обработке сигналов [1].

Для примера рассмотрим узкополосную MIMO РЛС с антенной решеткой (АР), у которой K=2 передающих и L=3 приемных элемента. Расположим элементы решетки так, чтобы получить V=KL=6 виртуальных приемопередающих элементов (см. рис. 1), расположенных в середине между Ti (i=1, ..., K) передающим и Rj (/=1, ..., L) приемным элементами.

На рис. 1 переменная dv обозначает расстояние между виртуальными элементами, 90 -направление на цель.

Передающие элементы синхронно излучают кодированные взаимно ортогональные

сигналы на одинаковой несущей частоте rrt •

Для упрощения будем рассматривать принятый сигнал, состоящий только из отраженного сигнала (ОС). Тогда принятый сигнал на]-м (] = 1,L) приемном элементе можно записать в виде:

С.А. ГОРШКОВ, ПИ. ОРГИШ

Военная академия Республики Беларусь Минск-57, 220057, Беларусь

Поступила в редакцию 15 июня 2012

Введение

к

где Mi (V) - комплексная огибающая ОС, 7 = 1,К ; (г(7 у) - время запаздывания ОС от 7-ого передающего элемента до у-ого приемного элемента.

О Передающий элемент X Приемный элемент О Виртуальный элемент

Рис. 1. MIMO РЛС с K=2 передающими и L=3 приемными элементами После перемножения с сигналом гетеродина (Ur(t) = exp{-](ю0 -юпр)г}) и обработки в согласованных фильтрах (СФ) получим:

F (t) = v(t - tr (h] - ]<Vr ]

i = 1, K, ] = 1, L .

Импульсные характеристики СФ1 и СФ2 согласованы с ЗМ 1!1(() и и2(7) соответственно. В таблице приведены значения времени запаздывания от Т7 (7—1,...,К) передающего до RJ (/—^...Д приемного элемента, а также время запаздывания для Vk (к— 1,.,V) виртуального элемента, принимая за начало отсчета координату элемента Т2, с учетом дальности до цели D.

Значения времени запаздывания ОС

Для реальных элементов Для виртуальных элементов

T1R1 6dv sin 60 + 2D + 5dv sin 60 11dv sin 60 + 2D V1 5.5dv sin 60 + 2D + 5.5dv sin 60 11dv sin 60 + 2D

c c cc

T1R2 6dv sin 60 + 2D + 3dv sin 60 9dv sin 60 + 2D V2 4.5dv sin 60 + 2D + 4.5dv sin 60 9dv sin 60 + 2D

c c cc

T1R3 6dv sin 60 + 2D + dv sin 60 7dv sin 60 + 2D V3 3.5dv sin 60 + 2D + 3.5dv sin 60 7dv sin 60 + 2D

cc cc

T2R1 0 + 2D + 5dv sin 60 5dv sin 60 + 2D V4 2.5dv sin 60 + 2D + 2.5dv sin 60 5dv sin 60 + 2D

cc cc

T2R2 0 + 2D + 3dv sin 60 3dv sin 60 + 2D V5 1.5dv sin 60 + 2D + 1.5dv sin 60 3dv sin 60 + 2D

cc cc

T2R3 0 + 2D + dv sin 60 dv sin 60 + 2D V6 0.5dv sin 60 + 2D + 0.5dv sin 60 dv sin 60 + 2D

cc cc

Из таблицы видно, что время запаздывания ОС для каждой пары передающих и приемных элементов равно времени запаздывания образованного ими виртуального элемента. Антенная решетка из V виртуальных элементов эквивалентна обычной приемопередающей V-элементной АР, что можно использовать для упрощения анализа MIMO АР и расчетов.

В (1) exp{-j(a0trj)} указывает на начальную фазу сигнала. Тогда, используя значения

таблицы, можно показать, что фазовые сдвиги между выходными сигналами СФ отличаются на:

2%

Am = — 2dv sin 9 .

X v

Из (2) следует, что разность фаз сигнала между двумя смежными приемными каналами всегда определяется удвоенным расстоянием между виртуальными элементами d=2dv.

Сигналы с выходов СФ перемножаются с комплексными весовыми коэффициентами и суммируются, формируя (t). Исходя из таблицы, значения комплексных весовых коэффициентов можно определить как:

.2% . . . „

j^T (Хтт + XRn )Sin 9 - -

wm,n = e X , m = 1, K; n = 1, L,

где хтт - координата передающего элемента; XRn - координата приемного элемента.

Ограничения, накладываемые на расстояние d, зависят от угла отклонения максимума диаграммы направленности (ДН) [3]. В общем случае, расстояние d чаще всего выбирают равным X / 2 . Поэтому, чтобы обеспечить единственность главного лепестка в ДН на прием, необходимо расположить приемные и передающие элементы так, чтобы обеспечить расстояние между виртуальными элементами dv = 0,25X .

Цель настоящей статьи - предложить алгоритм определения координат передающих и приемных элементов MIMO АР для получения линейной виртуальной АР, зная число передающих, приемных и виртуальных элементов, при соблюдении условия единственности максимума главного лепестка ДН.

Синтез алгоритма расчета координат элементов АР

Рассмотрим вариант, когда имеется К=2 передающих и L=3 приемных элементов, расположенных в форме равнобедренной трапеции (см. рис. 2). Пусть приемные элементы находятся на расстоянии dv = 0,5 X. Чтобы обеспечить единственность максимума ДН виртуальные элементы должны быть на расстоянии dv = 0,25Х (см. рис. 2, а). Необходимо определить расстояние между передающими элементами dt.

А.

Т/К,

ГТ/R,

d,=0.5d \

П5П

rwi

d=0SX

reí

в

а б

Рис. 2. К пояснению определения расстояния между передающими элементами: а - MIMO РЛС с K=2 передающими и L=3 приемными элементами; б - трапеция АВСD

Расстояние между первым и третьим приемными элементами - нижнее основание трапеции (отрезок DG) SL =X . Расстояние между первым и шестым виртуальными элементами -средняя линия трапеции SV = 5X/2 . Верхнее основание трапеции (отрезок АВ) - искомое расстояние между передающими элементами dt = SK. Учитывая свойство трапеции SV = (SK + SL) / 2, определим SK :

SK = 2SV -SL = 5X-X = 4X .

Рассмотрим общий случай. Пусть имеется К передающих элементов и L приемных. Желаемое число виртуальных элементов (каналов приема) - V.

Тогда (см. рис. 3) каждый ^ый передающий элемент и L приемных элементов образуют L виртуальных элементов с межэлементным интервалом dv. Тогда длина всей виртуальной решетки SV =(У -/ 2 . Если расстояние между приемными элементами d , то длина приемной решетки SL ={Ь -. Необходимо найти длину передающей решетки SK. Как видно из рис. 3, передающие, приемные и виртуальные элементы образуют равнобедренную трапецию, у которой (см. рис. 2, б) верхнее и нижнее основания равны длинам передающей и приемной решеток соответственно, а длина виртуальной решетки является средней линией этой трапеции.

T/Ri

к'

Rxi I I Rx2

!_™v j

áv=0.5d

d= 0.5).

SL

Рис. 3. К пояснению определения расстояния между передающими элементами Учитывая свойство трапеции, определим длину передающей решетки: SK = 2SV - SL = 2 • (V -/2 - (Ь - 1)d = -Ь^ .

Тогда расстояние между передающими элементами равно:

dt = V—Ь d . ' К -1

Расположим систему координат ХОУ как показано на рисунке.

с А >'

ТХ!

Тх,

dt

S,/2

Тхк

t- . .

¡_T/Rl] ¡_T/Rl] - [_T/Rl J ¡ T/Rv±] [T/Ry^J - [ T/Rv ]

rf\'=0.5rf

Rxt

Rx,

< d= 0.5X

SJ2

Rx,.

Рис. 4. MIMO антенная решетка в прямоугольной системе координат

Тогда координаты приемных и передающих элементов с учетом (3) и (4) будут рассчитываться по формулам:

yi = 0;

L -1

xi =--d + i  d, i = 0,1,...,L -1,

координаты приемных элементов;

yj = h;

V - L, .V - L, . r .

xj =--d + j-d, j = 0,1,...,K -1,

2

K-1

координаты передающих элементов.

В выражении (6) h - значение высоты трапеции. Как известно, высота трапеции не влияет на длину основания и средней линии, поэтому координату у для передающих элементов можно выбирать произвольно. Например, если задать h=0, то все элементы, приемные, передающие и виртуальные, будут располагаться на одной линии.

Геометрическую структуру решетки с расстоянием между виртуальными элементами dv = 0,25А можно получить, проводя расчет относительно передающих, а не приемных элементов (т.е. задав расстояние между передающими элементами dv = 0,5А,). Тогда уравнения (5) и (6) примут вид:

yj =0;

K-1

xj =-~

d + j • d, j = 0,1,...,K -1,

координаты передающих элементов;

y, = h;

V - к V - K x, =--d + i-d, i = 0,1,...,L -1,

2

L -1

координаты приемных элементов;

Если необходимо управлять ДН на передачу [4], или использовать принцип подреше-ток [5], или число передающих элементов больше приемных (для минимизации результирующего размера антенной решетки), то расчет необходимо проводить относительно передающих элементов (т.е. задать расстояние между передающими элементами dv = 0,5А,). В остальных случаях, нужно проводить расчет относительно приемных элементов по формулам (5) и (6).

Таким образом, для расчета геометрической структуры MIMO антенной решетки можно использовать алгоритм, блок-схема которого изображена на рисунке.

Начало

Определение значений К, L V

г ,

Расчет по формулам

Ошибка

Рис. 5. Блок-схема алгоритма синтеза геометрической структуры МГМОантенной решетки

Описание работы алгоритма. Сперва необходимо задать требуемое число передающих, приемных и виртуальных элементов (блок 1). Затем необходимо проверить, корректно ли задано число виртуальных элементов (блок 2). Максимально возможное число виртуальных элементов равно КЬ, минимально возможное - (К+Ь-1) (для случая, когда передающие и при-

2

емные элементы совмещены). Если заданные исходные данные попадают в интервал, то проверяется условие о необходимости управления диаграммой направленности антенны (ДНА) на передачу (блок 3). Если необходимо управлять ДНА на передачу, то расчет координат элементов производится по формулам (7) и (8) (блок 6). Если управлять ДНА на передачу не надо, то осуществляется сравнение числа передающих и приемных элементов, с целью минимизации результирующего размера антенной решетки (блок 4). Если число передающих элементов больше, чем приемных, то расчет проводится по формулам (7) и (8) (блок 6). Иначе - по формулам (5) и (6) (блок 7).

Результаты расчетов

Рассмотрим пример использования синтезированного алгоритма. Пусть длина волны А,=30 см, расстояние между элементами, относительно которых производится расчет d = 0,15А, = 15см, высота АР h = d, число передающих элементов К = 2, число приемных L = 4, требуемое число виртуальных элементов V = 8, требуется управлять ДН только на прием. Согласно алгоритму (рис. 5), расчет необходимо проводить по формулам (5) и (6). На рисунке представлен результат расчета координат элементов.

э с у ) О о у О V О | э с у

А А л А

X

о

Передающие элементы

Приемные элементы

Виртуальные элементы

- 0.45 - 0.36 - 0.27 - 0.18 - 0.09 0 0.09 0.18 0.27 0.36 0.45 X

Рис.6. Результат расчета по формулам

Значения координаты х (в метрах) для передающих элементов: -0,3, 0,3; приемных элементов: - 0,225, -0,075, 0,075, 0,225; виртуальных элементов: -0,262, -0,188, -0,113, -0,038, 0,038, 0,113, 0,188, 0,262.

- 90 - 80 - 70

- 50 - 40 - 30 - 20

10 20 30 40

60 70 80 90

Угол места, град

Рис. 7. Результат моделирования ДН на прием: сплошная линия - приемопередающая ДН рассчитанной АР с К=2 и Ь=4; штриховая линия - приемная ДН решетки с Ь=4 элементами; точечная линия - приемная ДН АР, состоящей из ¥=8 элементов.

Из рассчитанных значений координат элементов антенны видно, что расстояние между приемными элементами составляет 15 см = 0,5 А, расстояние между виртуальными элементами составляет 7,4 см = 0,25 А, результирующий размер антенной решетки определяется расстоянием между крайними передающими элементами и составляет 60 см = 2А,.

0.2

0.1

0

0.1

0.5

10 0

На рис. 7 представлены результаты моделирования ДН для синтезированной геометрической структуры.

Как видно из рис. 7, ширина ДН решетки с L = 4 элементами шире чем ДН рассчитанной АР с К = 2 и L = 4. При этом ширина ДН рассчитанной АР равна ширине ДН АР состоящей из V = 8 элементов.

Заключение

Имея K передающих элементов и L приемных, возможны два предельных случая, если приемные и передающие элементы совмещены (K=L), то число виртуальных элементов будет минимальным V = K + L — 1: если передающие (приемные) элементы находятся на расстоянии друг от друга Ld (Kd), то число виртуальных элементов будет максимальным V = KL . Таким образом, при заданных параметрах K, L, V, h , с использованием выражений (5) и (6) или (7) и (8) с помощью алгоритма, показанного на рис. 6, можно определить координаты передающих и приемных элементов антенной решетки для формирования ДН на передачу и прием в одной плоскости.

Также следует отметить, что одним из применений синтезированного алгоритма является возможность повышения живучести РЛС, когда передающие элементы, с менее дорогостоящими антеннами, выносятся за пределы приемной позиции. Как видно из (5)-(8), параметры виртуальной решетки не зависят от расстояния между передающей и приемной решетками как по оси Х так и по Y.

SYNTHESIS ALGORITHM OF GEOMETRICAL STRUCTURE OF MIMO RADAR

ANTENNA ARRAY

S.A. GORSHKOV, P.I. ORGISH

Abstract

In the given work the algorithm of calculation of co-ordinates of elements of MIMO radar antenna is presented. Initial data for algorithm are the number of transmitting, receiving elements, and also demanded number of channels of receiving.

Список литературы

1. Черняк В.С. // Прикладная радиоэлектроника. 2009. №4.

2. Jiane Li, Petre Stoica MIMO radar signal processing. New Jersey: A J. Wiley &sons inc., 2009.

3. Устройства СВЧ и антенны. Проектирование фазированных антенных решеток / под ред. Д.И. Воскресенского. М., 2003.

4. Горшков С.А., Оргиш П.И. // Докл. БГУИР. 2011. №6(60). С. 26-33.

5. Hongbin Li, Braham Himed// IEEE journal of selected topics in signal processing. 2010. Vol. 4, №1.