Научная статья на тему 'Алгоритмы диаграммообразования адаптивных антенных решеток в условиях многолучевого распространения радиоволн'

Алгоритмы диаграммообразования адаптивных антенных решеток в условиях многолучевого распространения радиоволн Текст научной статьи по специальности «Физика»

CC BY
1473
440
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
АДАПТИВНЫЕ АНТЕННЫЕ РЕШЕТКИ / МНОГОЛУЧЕВОЕ РАСПРОСТРАНЕНИЕ РАДИОВОЛН / АЛГОРИТМЫ ДИАГРАММООБРАЗОВАНИЯ

Аннотация научной статьи по физике, автор научной работы — Нечаев Ю. Б., Борисов Д. Н., Пешков И. В.

Проводится анализ различных алгоритмов определения угловых координат источников радиоизлучения в составе линейных и кольцевых антенных решеток в условиях многолучевого распространения радиоволн. Исследуются как классические методы формирования диаграммы направленности, так и с расширенными нулями при воздействии множества рассеянных в пространстве переотраженных сигналов.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Алгоритмы диаграммообразования адаптивных антенных решеток в условиях многолучевого распространения радиоволн»

УДК 621.396

АЛГОРИТМЫ ДИАГРАММООБРАЗОВАНИЯ АДАПТИВНЫХАНТЕННЫХ РЕШЕТОК В УСЛОВИЯХ МНОГОЛУЧЕВОГО РАСПРОСТРАНЕНИЯ РАДИОВОЛН

Проводится анализ различных алгоритмов определения угловых координат источников радиоизлучения в составе линейных и кольцевых антенных решеток в условиях многолучевого распространения радиоволн. Исследуются как классические методы формирования диаграммы направленности, так и с расширенными нулями при воздействии множества рассеянных в пространстве переотраженных сигналов.

Ключевые слова: адаптивные антенные решетки, многолучевое распространение радиоволн, алгоритмы диаграммо-образования.

borisov@sc.vsu.ru

peshkov_i_v@sc.vsu.ru

Ю.Б. НЕЧАЕВ1 Д.Н. БОРИСОВ2 И.В. ПЕШКОВ2

ОАО «Концерн «Созвездие», г. Воронеж

Воронежский государственный университет

e-mail:

nechaev ub@mail.ru

Введение.

Перспективным решением задач повышения помехозащищенности, емкости и пропускной способности беспроводных систем является использование адаптивных антенных решеток (ААР). Данный тип антенн, в отличие от фазированных антенных решеток (АР), позволяет осуществлять формирование не только максимума диаграммы направленности (ДН), но и нулей в произвольных направлениях, вследствие чего происходит максимальное подавление шума и помех при минимальных искажениях полезного сигнала [1]. Основной причиной снижения эффективности пространственной фильтрации ААР является множественное воздействие рассеянных в пространстве переотражен-ных сигналов [2]. Переотражения сигналов, вызванные различными рассеивателями, приводят к отклонению распределения поля в раскрыве АР от номинального и, как следствие, ведет к снижению точности алгоритмов определения угловых координат, а также установки нулей и максимумов диаграммы направленности антенной системы [3].

1. Параметрическое описание сигналов.

Пусть ААР состоит из N антенных элементов (АЭ), сигналы с которых умножаются на коэффициенты весового вектора w, после чего суммируются. Предположим, что M радиосигналов приходят на АР с произвольных направлений {[ т ¡ хп' , и описываются выражением [i]:

„/«>s(q„ П^,а'"(Р п() = ехр{/[2П ^sin(D„()]¡, где г - радиус АР, n — ^ =2G/DB; - длина волны да-го сигнала, 1) - межэле-

ментное расстояние. Для произвольной АР комплексный вектор сигналов на выходе антенных элементов описывается выражением [1]:

х(0 = ADs(O0.

где x(t) - N -мерный вектор, описывающий сигналы на выходе каждого антенного эле-s (t) •-»

мента АР, v ' - M -мерный вектор сигналов; n(t) - вектор шума пространственного канала и 11 -го канала ЦАР; А - N М матрица направляющих векторов. Далее обозначение вектора или угла с индексом 1 относится к полезному сигналу.

Если s (t) - стационарный случайный процесс, n(t) - гауссовский случайный процесс с нулевым средним и ковариационной матрицей 21 ( 2 - дисперсия шума,

I - единичная матрица), то сигналы и шум некоррелированы. Тогда пространственная корреляционная матрица может быть записана в виде [l]:

R = £[х(0хя (0] = ARm Ая , □ I, (1)

где E [...] - статистическое усреднение, R^ = E[s(t)sH (t)] - ковариационная матрица сигналов. На практике пространственная ковариационная матрица точно неизвестна и поэтому приходится ее получать из набора K временных отсчетов [l]:

- 1 к

R=—[\xH(k)x(k\

^ А- I

где X(k) - выборка сигнала с элементов АР в k -й момент времени, символ " " означает

усреднение по K выборкам.

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

min E{W Hxin}, при WH äj=l, (2)

w

где xt n - выборка сигнала с элементов АР, содержащая только помехи и шум.

Решение задачи (2) основано на использовании алгоритма обращения выборочной корреляционной матрицы (в англ. литературе - SMI (sample matrix inversion)) [4]:

W = IR1a 1.

Кроме того, для решения задачи (2) можно использовать алгоритм управляющий нулями ДН, предложенный Годара [5]:

wff=u[ ПАЯ(А □ Ая □ 21)\ где uT - декартовский базисный вектор.

2. Алгоритмы оценки угловых координат.

Основополагающим принципом диаграммообразования фазированных антенных решеток является возможность помещения максимума ДН на источник радиоизлучения (ИРИ). Поэтому при разработке подобных устройств эффективное получение информации о местоположении ИРИ является крайне важной. Далее рассматриваются алгоритмы оценки угловых координат источников радиосигналов.

2.1. Алгоритм MUSIC.

При реализации алгоритма MUSIC используется свойство ортогональности векторов шумового подпространства и н аправляющих векторов [5]:

Е„яа(П:„) = 0, □ „,!_{[ j, ..!=„}.

Координаты сигналов будут соответствовать максимумам функции [5]:

p/D ч _ äg (□ )а(Р )

ая (□ )E(;Ef »(□)

Если сигналы коррелированны, то метод MUSIC работает нестабильно [5].

2.2. Метод проецирования сигнального подпространства.

Данный метод основан на решении следующей нелинейной задачи оптимизации [5]:

{Г ,Т} = argminPExATP“,F

где P...PF - норма Фробениуса, T = (AHA)1 AH ] s . Задача минимизации решается в предположении определенного значения матрицы A [5]:

= arg (min Тг'П/Е . äEf}},

где П А = I f А)'АЯ - оператор проецирования на подпространство, образованное

столбцами матрицы A, Tr{...} - след матрицы.

Метод проецирования сигнального подпространства может использоваться, если сигналы коррелированны [5].

2.3. Метод Кейпона.

Координаты сигналов будут соответствовать максимумам функции [5]:

Р(П)=----------—----

aff(l)R1a(l)

Данный метод способен разрешать коррелированные сигналы. Причем его точность в целом ниже, чем у алгоритма MUSIC и метода проецирования подпространства [5].

3. Алгоритмы определения углов прихода радиосигналов в условиях многолучевого распространения.

Пусть базовая станция расположена на некоторой высоте, а подвижный объект на поверхности земли. В таком случае рассеиватели электромагнитных волн в основном оказываются расположенными вблизи подвижного объекта. Предположим, существует множество рассеивающих элементов, с каждым из которых связана переменная задержка распространения k (I) и переменный множитель k (I). Тогда сигнал на антенной решетке принимает вид [4]:

Л/ Ns ^

Х(0 = □ □П,яДОехр(У“„ДО)а(П„, , ад s (0,г(0, (3)

772 = 1 к=[

где тк (t) - амплитуда к -го луча сигнала т -го ИРИ, отраженного от рассеивателей, тк- (0 ~ фазовая задержка к -го луча сигнала т -го источника, Г т С тк — угол прихода к -го луча, Ns - число рассеивателей. Величины тк (I) и тк (I) носят случайный характер. Геометрия канала связи описывается моделью Асзтели (4). Тогда смещения углов прихода m -го ИРИ при отражении от Ns препятствий описываются выражением:

~ 1 N I N )

Г = L П k к = '

1 тк Л Т 1 K ’ К

Ns 1 BW ’ 2"’ '2

Проведем оценку влияния углового разброса ( BW) и числа рассеивателей Ns на

характеристики алгоритмов оценки угловых координат MUSIC, Кейпона и проецирования сигнального подпространства, которые являются представителями алгоритмов собственноструктурных, несобственноструктурных и проецирования подпространств соответственно. В (4) приводятся данные углового разброса в разных условиях распространения радиоволн: в пригороде угловой разброс принимает значения 5-80, а в городских условиях - до 200. В связи с этим в работе использовались крайние значения, которые могут быть достигнуты в реальных условиях.

3.1. Точность алгоритмов для кольцевой АР.

Моделирование алгоритмов определения углов прихода радиосигналов в условиях многолучевого распространения проводилось для кольцевой АР с межэлементным расстоянием, равным 0.5 и количеством антенных элементов (АЭ), равным 8. В каналы АР вносился шум со значением отношения мощности шума к мощности сигнала (ОСШ), равным 20 дБ, угловой разброс BW изменялся в пределах (1-26)0. При этом учитывалось разное количество рассеивателей, в частности, 10, 20 и 30 для каждого из которых проводилось 150 испытаний. В результате моделирования проводилась оценка среднеквадратического отклонения пеленгов (СКО ).

Зависимость CKO алгоритмов определения углов прихода радиосигналов (MUSIC, Кейпона, проецирования подпространства сигналов) в условиях многолучевого распространения для различных значений Ns, предполагает, что точечный сигнал излучается с азимутом 30° относительно первого антенного элемента.

СКО рассматриваемых алгоритмов пеленгования совпадают и практически не зависят от количества рассеивателей, причем в большей степени определяется угловым разбросом, а также, что зависимость СКО от углового разброса носит линейный характер. Кроме того, при воздействии одного сигнала, переотраженного несколько раз, величина СКО пеленгов имеет небольшое значение, что говорит об устойчивости кольцевой

антенной решетки к многолучевому распространению одного излучаемого сигнала.

Поскольку главной задачей адаптивной антенной решетки является подавление помеховых сигналов с минимальными искажениями полезного сигнала, то наибольший интерес представляет работа алгоритмов пеленгования при воздействии двух сигналов равной мощности. Это и является зависимостью СКО алгоритмов определения углов прихода радиосигналов, рассмотренных выше, в условиях многолучевого распространения при воздействии сигналов с азимутами их источников 30 ° и 70 ° относительно первого АЭ.

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

Таким образом, СКО алгоритмов определения углов прихода радиосигналов MUSIC, Кейпона и проецирования сигнального подпространства для кольцевой антенной решетки показывают приблизительно одинаковые значения, которые незначительно зависят от количества рассеивателей. При этом зависимость от углового разброса носит ярко выраженный характер. Суммарное значение СКО алгоритмов определения пеленгов для двух сигналов составляет 3,50.

3.2. Точность алгоритмов для линейной АР.

Моделирование алгоритмов определения углов прихода радиосигналов в условиях многолучевого распространения проводилось для линейной АР с межэлементным расстоянием, равным 0.5 и количеством АЭ, равным 8. Условия проведения численного эксперимента, значение шума, углового разброса, а также количество рассеивателей и испытаний были идентичны эксперименту с кольцевой АР.

Представленная зависимость СКО от углового разброса при воздействии одного сигнала, переотраженна разным числом рассеивателей. Исходный сигнал излучается из направления 30 ° относительно перпендикулярной оси, проведенной от центра АР.

Из представленных данных видно, что зависимость СКО обусловлена угловым разбросом, а не количеством рассеивателей. Для наихудших случаев работы алгоритмов для линейной АР СКО составляет не менее 7 °, в то время как, СКО кольцевой АР - не более 1.5 °.

Описанная зависимость СКО алгоритмов определения углов прихода радиосигналов в условиях многолучевого распространения при излучении сигналов из направле-

°°

ний 30 и 70 относительно перпендикулярна оси, проведенной от центра АР.

Анализ полученных данных показал, что алгоритм проецирования сигнального подпространства показал значения СКО превосходящие значения, полученные алгоритмами MUSIC и Кейпона. Для двух последних методов СКО практически идентично.

Как и для кольцевой АР, СКО анализируемых алгоритмов для линейной АР в большей

степени зависит от углового разброса, чем от количества рассеивателей. Значения СКО

для величины , имеющей значение не более 16°, сопоставимо со значениями для кольцевой антенной решетки для самого худшего рассматриваемого случая.

В целом можно заключить, что независимо от типа алгоритмы определения углов прихода радиосигналов (собственносруктруные, несобственносруктруные или проецирования подпространства) ведут себя схожим образом для линейной и кольцевой антенных решеток, т. е. их СКО незначительно зависят от количества рассеивателей и в большей степени определяются угловым разбросом. Причем линейная антенная решетка при определении угловых координат ИРИ не обладает устойчивостью к многолучевому распространению. Напротив кольцевая антенная решетка такой устойчивостью обладает, как при воздействии одного, так и двух исходных сигналов. Численный эксперимент показал, что суммарное значение СКО для двух сигналов составляло и ° для линейной АР и 3.5 °

для кольцевой АР, а также 8 ° и 1 ° для одного сигнала при угловом разбросе 26 ° для линейной и кольцевой АР соответственно.

4. Алгоритмы формирования ДН с расширенными нулями.

Для увеличения эффективности пространственной фильтрации ААР в условиях многолучевого распространения радиоволн применяются алгоритмы формирования ДН с расширенными нулями. Существует достаточно много алгоритмов расширения нулей ДН [2, 6, 7, 8, 9]. Однако у всех них есть один недостаток - они предназначены для работы только с линейными эквидистантными АР.

4.1. Диаграммобразование линейной АР.

4.1.1. Алгоритм формирования ДН Мэйлу-Затмана.

Для получения провала ДН шириной Ж в каждом из направлений т (для

т = 2..М 1) Мэйлу в [6] рассмотрел группу из I равномощных некогерентных источников сигнала, расположенных на прямой линии. В отличие от разбросанных дискретных фиктивных источников, которые предложил Мэйлу, Затман [7] использовал непрерывно распределенные фиктивные источники. Тогда модифицированная корреляционная матрица может быть выражена [8]:

6 = Я □ Т,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где Тп = ¡1] и б/ = Ъ„и,/2/0 = Ж/2 .

Для достаточно большого I оба алгоритма становятся идентичными [9].

4.1.2. Алгоритм формирования ДН Тафернера.

Здесь для получения весового вектора используется замена матрицы (1) на новую матрицу, которая вычисляется оконным сглаживанием направляющих векторов а( ) [8]:

) = а(Г„,) Xу(Гш для/?? = 1..М,

I

где У(1) = Ж> /15Ш( т),...,5Ш((# Гт)—=, а = (Ш/2)П С08(Ч, (для

Ф 1)а„, □ }

т = 2..М) - половина ширины нуля ЫЖ и N0 - константа, определяющая глубину ну-

_ □ 2%т(ВЖ зт((# \)В1¥^2)и

ля, --- ---—----—и , В1¥- = ВЖ~ С08(П[), ВЖ -

, 17 2П В1¥П ’"'’Ш 1)(2П ВШи)

1, V /V > ширина

□ С

главного луча ДН. Тогда новая корреляционная матрица записывается [8]:

Ы.та1'=АА ,А= 1 _а(П2),а(^),а((

4.1.3. Алгоритм формирования ДН Гершмана.

Алгоритм основан на свойстве ортогональности весового вектора подпространству помех, т. е. [ю]:

\у*П аЩ т), дляда = 2..М.

Для увеличения ширины нулей вводится ограничение р -го порядка [10]:

□ □

= 0, дляда = 2..Мир = 1...Р,

2 '

где □=— 8т(1 и, используя формулу а(П )= _1,ехр(7'й?П ехр(7'(# 1)^/

) )-•

получаем [ю]:

н □ {В'аЩ,, )},

где В = diag{1,ё,...,(N 1)ё} . Тогда корреляционная матрица для формирования широ-

ких нулей ДН записывается следующим образом [10]:

__ Р

я = П£Рвр квр,

р=о

где - определяет соглашение между глубиной и шириной нуля ДН.

4.1.4. Алгоритм формирования ДН Риба.

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

к

пространственная корреляционная матрица по углам прихода радиосигналов за время л

> 0 может быть выражена [2]:

_ м

ВД = С Й.[5(С„ (0))ая (Ц, (0))] и <ЗР „(0), (*)).□ =1, <}(□ „( (0), Г.; (*)) О

т=1

1 Зг2^-Л, ,

т

= ехр(2[Пс1(р д)уП-(к)со*2(0 „,(0))), где рт - мощность т -го радиосигнала, (0) - среднее значение угла прихода радиосигнала, которое является установленным пеленгом радиосигнала, 2 (к) - дисперсия пеленга т -го источника сигнала в момент времени к, р и ^ - номер строки и столбца матрицы, соответственно.

Для получения ДН с широкими нулями вводится условие оптимизации:

а^тт ^ н Я. (к)^ при ^ н Я. (к)^ = р1.

#

Данная задача решается следующим образом [2]:

(4)

w = е

тах

еmax Я <* (к)е„

где е тах - собственный вектор, связанный с максимальным собственным значением матриц Ял (к) и Я(к) .

Возьмем наихудший случай матрицы СК „,(0), 2 (к)), тогда [2]:

<Х^тах) = еХр(2р£/(^ ^)] ^т^К^ ),

где Ктах - интервал времени между адаптационными периодами, когда производятся

вычисления весового вектора, тах^тах - предполагаемая дисперсия углового разброса пеленгов. Тогда получаем [2]:

Й1(*) = ВД □ 0(Ктах),Щк) = Щ0) □ 0(Ктах).

Р

Серия История. Политология. Экономика. Информатика. 2012. №1 (120). Выпуск 21/1

После получения матриц Ял (к) и Я(к) из (4) вычисляется весовой вектор ^ [2].

4.1.5. Анализ алгоритмов диаграммообразования с расширенными нулями.

Оценка среднего значения отношения мощности полезного сигнала к мощности суммы помехи и шума (ОСПШ) алгоритмов диаграммообразования, устойчивых к подвижным источникам помех осуществляется для линейной эквидистантной АР с N = 8, расстояние между АЭ составляет 0.5 при воздействии равномощных некоррелированных сигналов с азимутами 20 ° и -20 ° относительно перпендикулярной оси, проведенной от центра АР. Предполагается, что помехой является сигнал, приходящий с направления -20 ° .В каналы вносится некоррелированный шум со значениями ОСШ, равными 30 дБ,

20 дБ, 10 дБ и 5 дБ. После построения матрицы Я производится ее регуляризация на 10 дБ над уровнем собственных шумов. При моделировании для каждого из алгоритмов были выбраны коэффициенты, которые, на взгляд авторов, в наибольшей степени соответствуют принятой окружающей обстановке. Для алгоритма Риба: ^ах = 4, регуляризация

отсутствует; для алгоритма Гершмана: р = 1, п = 1, регуляризация присутствует; для алгоритма Тафернера: N0 = 40, NW =20 °, БЖ = 1 °, регуляризация присутствует; для алгоритма Мэйлу-Затмана: Ж = 20 , регуляризация присутствует. Угловой разброс □ Й(Г

при этом изменяется в диапазоне 1-25 °, число испытаний составляет 200, Ns =15 . Значение ОСПШ определялось как отношение мощности суммы всех переотраженных компонентов полезного сигнала к мощности суммы всех переотраженных компонентов сигналов помех и шума [2]:

Получаем средние значения величины ОСПШ адаптивной антенной решетки в условиях многолучевого распространения радиосигналов в зависимости от углового разброса при различных ОСШ. Из приведенных данных видно, что лучшие результаты по подавлению распределенной в пространстве помехи показал алгоритм Тафернера. Алгоритм Мэйлу-Затмана показал стабильную эффективность к пространственной фильтрации в условиях многолучевого распространения во всем рассматриваемом шумовом диапазоне, но в пределе до 10 ° углового разброса.

Алгоритмы Риба и Гершмана не дали достаточной ширины провалов ДН, поскольку их работа основа на использовании матрицы Я, полученной без учета усредненных данных с выходов АР, распределение поля в раскрыве которой отличается от номинального, вследствие наличия переотраженных сигналов, например, таких как по методу Та-фернера. Поэтому для работы двух последних рассмотренных алгоритмов в условии многолучевого распространения необходимы дополнительные преобразования.

4.2. Диаграммообразование кольцевой АР.

В данной работе разработана процедура формирования расширенного нуля ДН для кольцевой АР, предполагающая случайное непрерывное распределение источников сигналов.

4.2.1. Гауссовское распределение.

Функция плотности вероятности при гауссовском распределении записывается [11]:

ОСПШ

□*-1

I

(5)

Серия История. Политология. Экономика. Информатика. 200

2012. №1 (120).mii2Sid№21/fl20). Выпуск 21/1

р(~ ) =------^=г—-------ехР((^ Г)2/2Ш2), дляП ],

erf(—)V2D

л/ С

2D

где - среднее направление прихода сигнала, - среднеквадратическое отклонение распределения. Согласно [11], точная корреляция между p -м и q -м антенным элементом кольцевой АР выражается:

к ^ ~ _ г-

Re[R{p,q)\== -f= Q expCj2 )[J0{Zc),J2k{Zc)cos{2k(U 4lUу,□ ))]ф,

лШ □ к=1

жг

к V2 ■ ■ _ □

Im[R(p,q)] = -±= □ехр(/ )[2[l/2i.1(Zc)sin((2Aril)(n , -JlUyUWy.

л/ П - ¿-=о

и

Для небольшого или среднего значения , можно принять, что [н]:

sin(2_ ) □ zD , cos(zD) □ дляг = 1,2,..., □.

1

Тогда запишем приближенные выражения для корреляции между элементами кольцевой АР [11]:

R(p,q) □ К exp(./z с sin(D □ ))exf>J

i-------- R R

где Zc = 1/Z2iZ22, Zj =2П— [cos(q,)] cos(U g)], Z2 = 2D—[sin(l^) sin(CI g)],

sin(D ) = Zj/Zc , cos(C ) = Z2/Zc .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4.2.2. Равномерное распределение.

Функция распределения плотности вероятности записывается в следующем виде [11]:

1

р( ) =— для □ [□ U],

20

где 2 - диапазон углов вокруг центральной угловой координаты .

Из [и] действительная и мнимая части точной корреляционной функции кольцевой АР может быть записана:

Re[R(p, q)} = J0 (Zc) J2k (Zc) cos(2k(□ □)) sin(2kO),

k=1

Im[R(p, q)\ = 2 □ J2kj (Zc) cos((2£ 1 )(□ С ))sinc((2k 1) _).

k=0

Для небольшого или среднего значения, можно принять, что [и]:

R(p,q) □ exp(j'Zc sin(D Г ))sinc{Zc L cos(l Г )). (7)

В [и] приближенные функции межэлементной корреляции показывают достаточно хорошее взаимное соответствие с точными функциями при <10° для равномерного распределения и при <10° для нормального распределения без ограничения на отношение R/ , а также для R! <0.25 при < 30° и < 30°. При таких ограничениях приближенная функция может быть использована вместо точной.

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

cosfD ПУ 2

(6)

Серия История. Политология. Экономика. Информатика. 201

2012. №1 (120). Выпуск 21/1

основании законов корреляции (6) и (7) для кольцевой антенной решетки можно записать аналитическое выражение получения вектора весовых коэффициентов для формирования ДН с расширенными провалами в направлениях помех [12]:

Rj, □ R», □ „f nP a{D j),

□ m=2 □

где Rj = a( I, )a( \ )H - пространственная корреляционная матрица полезного сигнала,

R m - пространственная корреляционная матрица, pq -й элемент которой получен из приближенной функции межэлементной корреляции (6) или (7).

4.2.3. Моделирование разработанной процедуры расширения нулей диаграммы направленности.

Кольцевая АР состоит из 8 антенных элементов с межэлементным расстоянием, равным половине длины волны при воздействии равномощных сигналов с угловыми координатами 20 ° и 60 ° относительно первого АЭ. В каналы вносится некоррелированный шум со значениями ОСШ, равными 30 дБ, 20 дБ, 10 дБ и 5 дБ. Формирование ДН производится из предположения как равномерного, так и гауссовского распределения источников радиоизлучения по угловой координате. Для сравнения были использованы алгоритм SMI с регуляризацией, равной 10 дБ относительно собственных шумов, а также алгоритм, управляющий нулями ДН. Количество испытаний для каждого ОСШ и угла смещения составляет 200, ОСПШ вычисляется из выражения (5). Сигнал на антенной решетке принимает вид выражения (3). Углы тк получены при использовании модели Асзтели. Для условия многолучевого распространения значения и выбраны максимальными для данных условий приближенных выражений корреляции (6) и (7), т. е. 10 °, т.к. угловой разброс ш составляет 1-25 , число испытаний составляет 200, ОСПШ получено из (5), для всех BW количество рассеивателей Ns =15. Получаем результаты работы

алгоритмов расширения нулей ДН кольцевой АР в условиях многолучевого распространения.

Кривые средних ОСПШ для ДН с расширенным нулем имеют малую крутизну в отличие от кривых ОСПШ для обычных алгоритмов диаграммообразования, что является подтверждением того, что алгоритмы, расширяющие нули ДН, имеют хорошую устойчивость к распределенным в пространстве сигналам. Однако в их работе, как отмечалось выше, имеется недостаток, заключающийся в слабом подавлении некоррелированного шума, где кривые ОСПШ для ДН с расширенным нулем практически не изменяют своего значения с увеличением угла разброса, но имеют меньшее значение ОСПШ в сравнении с традиционными алгоритмами формирования ДН.

Получаем зависимость ОСПШ при многолучевом распространении сигналов с числом рассеивателей, равным 15, от углового разброса для алгоритмов расширения нулей ДН кольцевой АР с гауссовским и равномерным распределениями для = 1°, =1°

и С =10",

□ =10°.

Как видно, результаты подобны тем, которые рассмотрены ранее в том смысле, что лучшее подавление разбросанной в пространстве помехи дают способы с равномерным распределением. Также анализ показывает, что кривые ОСПШ для распределения с = 10°, = 10° слабо изменяются с увеличением угла разброса, в то время как, кривые

ОСПШ для распределений с = 1°, =1° имеют заметное снижение с увеличением разброса. Кроме того, формирование ДН для распределения с =10°, =10°, как и прежде,

хуже справляется некоррелированным шумом.

Выводы.

Проведена оценка алгоритмов определения углов прихода радиосигналов в условиях многолучевого распространения в составе кольцевой и линейной АР. Установлено, что оценки углов прихода радиосигналов не зависят от числа рассеивателей и типа алго-

НАУЧНЫЕ ВЕДОМОСТИ Серия История. Политол?)1'йя.Э|1Ь(н1ом1и1кВ.ып1ускр1м,атика.

2012. №1 (120). Выпуск 21/1

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

Список

литературы

1. Litva J. Digital beamforming in wireless communications / J. Litva, T. K.-Y. Lo. — Boston : Artech House, 1996. — 301 p.

2. Riba J. Robust beamforming for interference rejection in mobile communications / J. Riba,

J. Goldberg, G. Vazquez // IEEE Trans. Signal Processing. — 1997 vol. 45. — 271-275 p.

3. Asztely D. On antenna arrays in mobile communication systems : fast fading and GSM base station receiver algorithms : technical report / D. Asztely. — Stockholm, 1996. — 5 p.

4. Reed I. S. Rapid convergence rate in adaptive arrays / I. S. Reed, J. D. Mallett, L. E. Brennan

IEEE Trans. Aerosp. Electron. Syst. — 1974. — Vol. AES-10. — 853-863 p.

5. Godara L. Application of Antenna Arrays to Mobile Communications, Part II: Beam-Forming and Direction-of-Arrival Considerations / L. Godara // Proceedings of the IEEE. — 1997. — Vol. 85. -No 8. — P. 11951245.

6. Mailloux R. J. Covariance matrix augmentation to produce adaptive array pattern troughs. — Electron. Lett. -- 1995. -- vol. 31 -- No. 10. -- 771-772 p.

7. Zatman M. Production of adaptive array troughs by dispersion synthesis. - Electron. Lett. -

1995. — vol. 31 — no. 25. — 2141-2142 p.

8. Taferner M. A Novel DOA-based Beamforming Algorithm with Broad Nulls / M. Taferner, A. Kuchar, M.C. Lang, M. Tangemann, C. Hock. — 10th International Symposium on Personal, Indoor and Mobile Radio Corn. PIMRC '99, Osaka, Sept., 1999. — P. 342-247.

9. Guerci J. R. Theory and application of covariance matrix tapers to robust adaptive beamform- ing. - IEEE Trans. Signal Processing. — 2000. — vol. 47 — 977-985 p.

10. Gershman A. B. Adaptive beamforming algorithms with robustness against jammer motion

/ A. B. Gershman, U. Nickel, and J. F. B o hme. — IEEE Trans. Signal Processing. - 1997. - vol. 45. - 18781885 p.

11. Performance of MIMO with circular antenna array using correlation matrix / J. Zhou [et. al.]

// Asia-Pacific Conference on Communications (APCC '06). - Busan, 2006. - P. 15.

12. Nechaev Yu.B. Beamforming Algorithm for Circular Antenna Array Immune to Multipath

Propagation and Non-Stationary Interference Sources / Yu. B. Nechaev, D. N. Borisov, and I. V. Peshkov

// Radioelectronics and Communications Systems. - 2011. - Vol. 54. - No 11. - P. 604612.

научные ведомостиВЕАЙ^ММШШ ITHMSOF SMARTANTENNAS IN 203

MULTIPATH FADINGCHANNELS

YU.B NECHAEV1 D.N. BORISOV2 I.V. PESHKOV2

1)JSCConcern"Sozvesdie", Voronezh

Voronezh State University e-mail: nechaev_ub@mail.ru

borisov@sc.vsu.ru peshkov_i_v@sc.vsu.ru

Analysis of different kinds of DOA-estimation algorithms with linear and circular antenna arrays in multipath fading channels is fulfilled. Both the clas- sical beamforming algorithms and the broad-nulls algorithms is researched when influenced by multiple spread signals.

Key words: smart antennas, multipath fading, beamforming algorithms.

i Надоели баннеры? Вы всегда можете отключить рекламу.