МОДЕЛИРОВАНИЕ ЭФФЕКТА ВЗАИМНОГО ВЛИЯНИЯ МИКРОПОЛОСКОВЫХ ЭЛЕМЕНТОВ УКВ-ДИПАЗОНА
И.В. Пешков, канд. физ.-мат. наук, доцент Н.А. Фортунова, канд. техн. наук, доцент В.А. Жигулин, магистрант
Елецкий государственный университет им. И.А. Бунина (Россия, г. Елец)
DOI:10.24412/2500-1000-2023-11-4-97-103
Исследование выполнено за счёт гранта Российского научного фонда № 23-2100125, https://rscf.ru/project/23-21-00125/
Аннотация. Эффект взаимного влияния антенных элементов способен оказывать существенное воздействие на характеристики сложных, чувствительных систем. Данный эффект изменяет форму диаграмм направленности (ДН) антенных элементов, меняет их активное сопротивление, что усложняет согласование импедансов антенных решеток (АР). В областях со строгими требованиями взаимное влияние может оказать очень существенное влияние на результаты работы системы. В частности, они способны влиять на результаты работы цифровых алгоритмов оценки угловых координат источников радиоизлучения со сверхразрешением. Проведено исследование взаимного влияния микро-полосковых антенны диапазона УКВ.
Ключевые слова: взаимное влияние, антенные решётки, цифровые антенные решётки, микополосковые антенны, моделирование, диаграмма направленности.
Большая часть работ в области изучения цифровых антенных решёток выполнена с использованием классической теории. Фундаментальным для этого подхода является предположение, что взаимное влияние (ВВ) между антенными элементами (АЭ) пренебрежимо мало. Во многих случаях это предположение оправдано. С другой стороны, в областях со строгими требованиями, ВВ должно учитываться при проектировании. Области, представляющие особый интерес, включают высокочувствительные антенны для медицинской техники, спутниковой связи, дистанционного зондирования и радиоастрономии [1].
Ключевыми причинами возникновения взаимного влияния можно считать следующие: недостаточная изоляция; снижение эффективности антенны; уменьшение коэффициента усиления антенны; высокая корреляция между излучателями [2].
Рассмотрим воздействие взаимного влияния относительно цифровых антенных решёток [3]. Решаемые задачи, направленные на достижение цели: провести моде-
лирование антенных решёток с разными видами антенных элементов. Провести моделирование алгоритмов для антенных решёток с разными видами антенных элементов: диполями и патч-антеннами. Провести моделирование алгоритмов с учётом таких факторов, таких как S-параметры решётки и искажённые диаграммы направленности элементов.
Принцип работы цифровых антенных решёток
Существует достаточно большое количество алгоритмов оценки угловых координат радиосигналов для цифровых антенных решёток. Их целью является определение мощности сигнала на выходе АР в зависимости от угла прихода сигнала. По максимумам этой функции, называемой псевдоспектром, можно оценить, с какого направления поступает радиосигнал. Достаточно широко известны алгоритмы Барлетта, Кейпона, MUSIC и т.д.
На рисунке 1 представлена обобщённая схема процесса радиопеленгации цифровой антенной решётки. Её можно разделить на две части: аналоговую и цифро-
вую. В аналоговой части происходит при- реляционной матрицы, а затем и псевдо-ём сигналов антенной решёткой. В цифро- спектра. вой части происходит расчёт сначала кор-
Аналоговая часть Цифровая часть
Рис. 1. Общая схема системы цифровой пеленгации
Далее рассматривается модель работы алгоритмов пеленгации в соответствии с данной схемой на рисунке 1. Пусть имеется линейная антенная решётка, состоящая из М антенных элементов. Она принимает Б независимых сигналов из Б разных направлений с углами 0б, как показано на
рисунке 1. Здесь &>(£) - падающий сигнал, хм(к) - принятый сигнал с М-го антенного элемента, у(к) - итоговый выходной сигнал. В итоге формула для матрицы выходного сигнала решётки имеет следующий вид:
х = А ■ s + п (1)
где п - матрица шума; А - матрица сканирующих векторов, рассчитываемая по следующей формуле:
A = О (в) q (в2 )-q (Gd )]
(2)
где q (бв ) - это сканирующий вектор из Мэлементов для направления 0:
q (в)=[1
g jn sin в gj (М—1) 71 sin в
]
(3)
Здесь представлена формула для линейной решётки. Учитывается сдвиг фаз в антенных элементах, вызванный их расположением (с расстоянием между элементами в 0,5Х). Корреляционную матрицу определяется по следующей формуле:
Rx =
1 N
— УХ^к) • x,(k)*
1 N
ТгУхм(k) • х1(к)*
N к=1
1 N
N k=1
1 N
ттТхм(к)-хм(кУ
N k=1
(4)
где N - число отсчётов цифрового сигнала, хм (к)* - эрмитово-сопряжённая матрица.
Далее на основе корреляционной матрицы рассчитывается псевдоспектр. Ниже представлены формула для методов Барлетта [6], Кейпона [7, 8] и MUSIC [9] соответственно:
pb (0) = a(0) Rja (0)
pc (0) = - 1
Pmu(0)=p
(0)H R-1a (0)
1_
aH (0) ENEHNa (0)|
(5)
(6)
(7)
где a (0) - матрица сканирующих векторов для всего диапазона пеленгации, EN -собственный вектор шумов. Выражение (5) соответствует методу Барлетта, формула (6) - методу Кейпона. Метод MUSIC (7) несколько отличается от представленных выражений. В его формуле используется не корреляционная матрица (4), а вычисляемый на её основе собственный вектор шумов Еж. Подробнее математические преобразования рассмотрены в [1, 2].
Взаимное влияние и метод встроенных элементов (embedded elements)
Одним из условий в классическом анализе антенных решёток является исключе-
ние взаимной связи - предполагается, что ДН всех антенных элементов идентичны. Обычно в расчётах принимается ДН изолированного элемента. На рисунке 2а представлен пример подобной ДН. Такого рода приближение справедливо для антенных решёток с большим числом элементов. Однако для решётки с относительно небольшим количеством элементов взаимное влияние оказывает более сильное влияние и, как следствие, ВВ существенно искажает диаграммы направленности отдельных элементов (рис. 2б). Следовательно, общая ДН решётки больше не является просто произведением коэффициента решётки и ДН изолированного элемента.
а) б)
Рис. 2. Диаграммы направленности патч-антенны а) изолированного антенного элемента, б) с учётом взаимного влияния
Для проведения более приближённых к реальности расчётов требуется вычисление диаграмм направленности для каждого отдельного элемента антенной решётки. Диаграмма направленности отдельного элемента с учётом взаимного влияния представляет собой ДН антенной решётки, в которой исследуемый элемент считается активным, а все остальные элементы являются пассивными с отключёнными выходами. Такой подход принято называть вычисление по методу встроенного эле-
мента (embedded element). Ключевым моментом является то, что вычисленная таким образом ДН отдельно АЭ - это диаграмма направленности, учитывающая всю антенную решётку, включая рассеяние вблизи активного элемента.
Моделирование антенных решёток
Было проведено моделирование антенных решёток согласно методу моментов по методу встроенного элемента (рис. 2б). После расчёта каждого типа антенны вы-
полнялась оптимизация для того, чтобы максимизировать коэффициент усиления (КУ) и минимизировать параметр S11 при импедансе, равном 50 Ом. Моделирование проводилось согласно следующим параметрам:
- количество элементов решётки M = 6;
- частота принимаемого сигнала -5 ГГц.
Ниже представлены полученные модели и их характеристики.
Таблица 1. Параметры моделей антенных решёток
Номер модели
Внешний вид
Описание
Материал антенн
Материал диэлектрика
ЛАР, элементы - прямоугольные патч-антенны. элементы оптимизированы, КУ = 10 дБ, S11 = -21 дБ на f = 5 ГГц
идеальный проводник
воздух
ЛАР, элементы - прямоугольные патч-антенны, элементы оптимизирова- медь ны, КУ= 1,8 дБ, S11 = -11 дБ на f = 5 ГГц
FR4
ЛАР, элементы - прямоугольные патч-антенны, элементы оптимизирова- медь ны, КУ= 4,8 дБ, S11 = -17 дБ на f = 5 ГГц
RO4730JXR
1
2
3
По результатам проведённых экспериментов было проведено сравнение нескольких антенных решёток с различными видами антенных элементов. Наилучшие результаты в значениях отклонений и минимально различимых углах показали линейные антенные решётки с элементами -
полосковыми антеннами с воздушным диэлектриком (№1, Таблица 1). Значения параметров составили 4 дБ и 10-14° соответственно. Данный вид АЭ обладает самыми высокими значениями коэффициента усиления (рис. 3).
а) 16,23 дБ б) 11,74 дБ в) 9,18 дБ
Рис. 3. Сравнение коэффициента усиления главного лепестка ДН для АР с разными видами диэлектрика а) воздух, б) FR4, в) RO4730JXR
Более низкие результаты показали элементы с диэлектриком марки RO4730JXR (мощность появления отклонений 10 дБ,
минимально различимые углы 20°). Ещё более низкие результаты показала АР с диэлектриком FR4 (20 дБ и 38°).
Разница в результатах между антенными решётками вызвана различиями в параметрах материалов диэлектриков. При моделировании использовались следующие материалы подложек:
- Воздух
- FR4 (простой стеклотекстолит)
- RO4730JXR (ламинат антенного класса)
Основными характеристиками материала подложки являются относительная диэлектрическая проницаемость 8г и тангенс угла диэлектрических потерь tg8 [3]. Значения данных параметров, использовавшиеся при моделировании, представлены в таблице 4.
Таблица 4. Параметры материалов диэлектриков АЭ
Название Диэлектрическая проницаемость, £г Тангенс угла диэлектрических потерь, tgб
Воздух 1 0
РЯ4 4.8 0.026
Я047301ХЯ 3 0.0023
Согласно [2], коэффициент усиления антенны уменьшается с увеличением тангенса угла потерь для всех частот. Потери мощности в диэлектрике увеличиваются с ростом частоты, а также при использовании материалов с высокими tg8 и 8г (подробно данная зависимость рассмотрена в [4]). Согласно данным, представленным в таблице 4, значения tg8 и 8г у БЯ4 выше, чем у RO4730JXR и воздуха, что является причиной более низких результатов ан-
тенных решёток с данным типом диэлектрика.
В связи с этим, материалы, применимые в ВЧ-УВЧ диапазонах, при использовании в УКВ диапазоне показывают заметно меньшие результаты по сравнению с идеальным случаем (воздух) и специализированными материалами (RO4730JXR). Согласно [5], предпочтительными для использования в СВЧ технике являются материалы с 8г=2-10 и tg8=10-3-10-4.
Рис. 4. Пример распределения поверхностного тока в АР
Рассмотрим, как распределяется ток в микрополосковом антенном элементе. Наибольшее значение наблюдается в точке крепления линии питания (точка Т1 на рис. 4). Высокие значения тока располагаются между правой и левой гранями излучателя (грани, параллельные оси У на рис. 4), достигая максимума на самих гранях (точки Т2 и Т3). Меньшие значения
тока наблюдаются на гранях, параллельных оси X (точки Т4 и Т5).
Рассмотрим, как эффект взаимного влияния проявляется в распределении поверхностных токов в антенных решётках. На рисунке 5 изображена антенная решётка №3 из таблицы 1 и график распределения тока.
а) Значения тока в разных точках АР
Current distribution
Y(m)o
A1 0.18
б) Расположение точек на АР Рис. 5. Распределение токов в антенной решётке
Рассмотрим полученные данные. Значения тока в точках А1, А2 и А3 (места крепления линии питания) у 1, 2 и 3 элемента различаются между собой. Также отличаются токи в точках В1 и В2 (правый край излучателя) у 1 и 3 элементов. В целом видно, что распределение тока в каждом антенном элементе отличается. Важно отметить, что картина изменений в распределении токов симметрична между 1 -3
проявлением воздействия эффекта взаимного влияния между элементами антенной решётки.
Самый низкий результат продемонстрировали решётки 2 и 3 из таблицы 1, где АЭ расположены плотно друг к другу и, следовательно, имеют наибольшее взаимное влияние. Они имеют одинаковую форму антенных элементов и отличаются только материалом изготовления.
и 4-6 антенными элементами. Это является
Библиографический список
1. Оценка вероятности появления ложных пиков кольцевых и концентрических антенных решеток при радиопеленгации со сверхразрешением / Ю.Б. Нечаев, И.В. Пешков, Ю.О. Аальмуттар Атхеер, К.Д. Аль Хафаджи Сармад // Вестник Воронежского государ-
ственного университета. Серия: Системный анализ и информационные технологии. -2016. - № 2. - С. 16-23.
2. Frank B. G. Smart Antennas for Wireless Communications: With MATLAB. - New York: McGraw-Hill Professional, 2005. - 288 c.
3. Анпилогов, В.Р. Диссипативные потери в микрополосковых линиях и микрополоско-вых антеннах / В.Р. Анпилогов, И.В. Зимин, Ю.Н. Чекушкин // Ракетно-космическое приборостроение и информационные системы. - 2018. - Т. 5, № 3. - С. 60-69.
4. Analysis of loss tangent effect on Microstrip antenna gain / Prasanna Puttaswamy, Pramod Srivatsa Krishna Murthy, Bindu Anil Thomas // International Journal of Applied Science and Engineering Research. - 2014. - № 3.
5. Богданов, Ю. Фольгированные диэлектрики - как выбрать оптимальный вариант для печатных плат ВЧ/СВЧ-диапазонов / Ю. Богданов, В. Кочемасов, Е. Хасьянова // Печатный монтаж. - 2013. - № 2. - С. 156-168.
MODELING THE EFFECT OF MUTUAL INFLUENCE OF MICROSTRIP ELEMENTS
OF THE VHF BAND
I.V. Peshkov, Candidate of Physical and Mathematical Sciences, Associate Professor
N.A. Fortunova, Candidate of Technical Sciences, Associate Professor
V.A. Zhigulin, Graduate Student
Yelets State University named after I.A. Bunin
(Russia, Yelets)
Abstract. The effect of mutual influence of antenna elements can have a significant impact on the characteristics of complex, sensitive systems. This effect changes the shape of the radiation patterns (DP) of the antenna elements, changes their active resistance, which complicates the matching of the impedances of the antenna arrays (AR). In areas with stringent requirements, mutual influence can have a very significant impact on system performance. In particular, they are capable of influencing the results of digital algorithms for estimating the angular coordinates of radio sources with super-resolution. A study of the mutual influence of microstrip antennas in the VHF range was carried out.
Keywords: mutual influence, antenna arrays, digital antenna arrays, mycostrip antennas, modeling, radiation pattern.