Разработка строгих математических моделей для плоских щелевых антенн бегущей волны с использованием тензорной функции Грина и их экспериментальное обоснование
- Авторы: Нефёдов Е.И.1, Заярный В.П.2, Кирьянова Н.И.3, Пономарев И.Н.3
-
Учреждения:
- Фрязинский филиал Института радиотехники и электроники имени В.А. Котельникова РАН
- Волгоградский государственный технический университет
- Волгоградский государственный университет
- Выпуск: Том 24, № 1 (2021)
- Страницы: 22-31
- Раздел: Статьи
- URL: https://journals.ssau.ru/pwp/article/view/8693
- DOI: https://doi.org/10.18469/1810-3189.2021.24.1.22-31
- ID: 8693
Цитировать
Полный текст
Аннотация
Произведена разработка строгих математических моделей для антенн бегущей волны с постоянной шириной щели на основе тензорной функции Грина с учетом как основной, так и кросс-поляризационной составляющей поля в дальней зоне излучения. Обоснована возможность использования этих моделей в качестве базовых моделей для исследования электродинамических свойств антенн бегущей волны с другой конфигурацией щели (линейно-, экспоненциально расширяющиеся и др.). С применением полученных математических моделей изучалось влияние геометрических параметров исследовавшихся антенн бегущей волны на их характеристики. Установлено, что для электродинамической плоскости Е при увеличении длины и ширины щели главный луч диаграммы направленности сужается, что приводит к повышению направленных свойств антенн. В электродинамической плоскости Н ширина главного луча диаграммы направленности уменьшается только при возрастании длины щели антенн бегущей волны и практически не зависит от ее ширины. Показано, что полученные теоретические результаты исследования хорошо согласуются с результатами экспериментов, которые также позволили выявить ограничения применимости разработанных моделей. Исследование характеристик антенн бегущей волны производилось на средней частоте 10 ГГц в диапазоне частот 8–12 ГГц.
Полный текст
Введение
В современной радиоэлектронной аппаратуре, в частности в составе систем сверхбыстрой обработки информации с применением объемных интегральных схем (ОИС), являются востребованными плоские щелевые антенны (ЩА) различных конфигураций [1; 2]. Известны различные виды подобных ЩА, работающих в режиме продольного осевого излучения и поэтому принадлежащих к классу антенн бегущей волны (АБВ). К ним относятся антенны с постоянной шириной щели, с линейно расширяющейся щелью, с экспоненциально расширяющейся щелью (антенна Вивальди) и щели со специальным расширением в виде плавного согласующего перехода на входе АБВ со свободным пространством на ее выходе. Конструкция таких антенн обеспечивает высокий коэффициент усиления, низкий уровень боковых лепестков и достаточно широкую полосу частот [3; 4]. Подобные антенны особенно востребованы при разработке ОИС микроволнового и оптического диапазонов, поэтому дальнейшее изучение их электродинамических свойств представляется важным и актуальным.
В данной статье внимание уделено приоритетному исследованию электродинамических свойств щелевых антенн с постоянной шириной щели (ПЩА) по следующим причинам. Во-первых, именно для этих АБВ оказалась возможной разработка строгих математических моделей для расчета поля на основе тензорной функции Грина, учитывающих как основную, так и кросс-поляризационную составляющую поля излучения. Во-вторых, модели для ПЩА применимы для разработки математических моделей АБВ, имеющих другую конфигурацию щели, для которых возможно их представление в виде набора регулярных участков постоянной ширины, скачкообразно изменяющих ширину через определенные интервалы их длины в соответствии с законом изменения щели данной АБВ [5–7]. В-третьих, ПЩА относительно просты и технологичны в изготовлении, что способствует качественному производству требуемого количества образцов АБВ различной конфигурации для проведения экспериментальных исследований.
1. Математическая модель поля излучения щелевой антенны с постоянной шириной щели на основе тензорной функции Грина
Одна из первых электродинамических моделей для расчета поля излучения ПЩА была предложена в [8; 9]. Однако она не учитывает кроссполяризационную составляющую поля излучения, и, как следствие, в ней отсутствует соответствующая компонента тензорной функции Грина. В данной статье разработана строгая математическая модель поля излучения ПЩА, учитывающая обе компоненты функции Грина, что позволяет определить полное электромагнитное поле в дальней зоне излучения. При этом в ней физическая конструкция ПЩА заменялась идеально проводящей бесконечной полуплоскостью, возбуждаемой внешним сторонним электрическим полем в щели постоянной ширины, врезанной перпендикулярно краю полуплоскости (рис. 1). В такой модели первоначально устанавливается поле в апертуре регулярной щели, а затем решается задача возбуждения идеальной полуплоскости этим полем.
Рис. 1. Идеальная полуплоскость с постоянной по ширине щелью и сторонним полем
Fig. 1. An ideal half-plane with constant width slot and its external electric field
В общем случае обе компоненты электрического поля излучения для таких щелевых АБВ были получены в [10] и в сферической системе координат (рис. 1) имеют вид:
(1)
(2)
где – внешнее стороннее электрическое поле возбуждения в апертуре регулярной щели ; и – соответственно, поперечная и продольная координаты источников стороннего поля в пределах щели; и – элементы тензорной (диадной) функции Грина; – коэффициент фазы (волновое число). При этом является основной составляющей поля излучения, а – составляющей кросс-поляризации.
Функции и связывают поле в точке источника с полем излучения в точке наблюдения (рис. 1) и определяются следующими соотношениями:
(3)
(4)
где
– комплексный интеграл Френеля; – аргумент комплексного интеграла Френеля. В элементе в первом слагаемом знак «+» берется для случая , а знак «–» – для случая . Интегрирование в (1), (2) производилось по поверхности регулярной щели (рис. 1), при этом согласующий переход Lf и питающая щелевая линия не учитывались.
Из выражений (3), (4) следует, что элементы тензорной функции Грина и состоят из двух слагаемых, которые разным образом определяют действие внешних сторонних источников поля в точке и формируют поле излучения в (1), (2). Первые слагаемые в и при обращаются в нуль. Поэтому они определяют поле излучения от внешних источников , расположенных практически во всех точках щели исключая краевые (при и некоторую область возле края щели шириной в пределах которой для в интервале значением интеграла Френеля можно пренебречь и считать . Следовательно, для короткой щели длиной и для всех точек наблюдения поля , расположенных на краю идеальной полуплоскости при и , а также в Е-плоскости при и это слагаемое обращается в нуль и на поле излучения не влияет. Вторые слагаемые в и наоборот, при достигают бесконечного значения и вносят наибольший вклад в поле излучения для источников , расположенных на краю излучающей щели . По мере удаления источников от края щели вдоль оси вклад этого слагаемого в поле дальней зоны монотонно уменьшается по закону . Второе слагаемое обеспечивает выполнение закона сохранения энергии для точек наблюдения , приближающихся к краю полуплоскости при и .
Таким образом, из выражений (3), (4) следует, что в формировании основной компоненты поля излучения в плоскости Н при и в элементе участвуют оба слагаемых, а в плоскости Е при и – только второе слагаемое. При этом как в плоскости Н, так и в плоскости Е элемент обращается в нуль. Такое поведение основной и кросс-поляризационной составляющих является характерным свойством излучения для идеальной полуплоскости, возбуждаемой щелью с любым законом расширения.
В качестве внешнего стороннего поля возбуждения для постоянной щели в (1), (2) использовалось квазистатическое приближение. Согласно [6; 7], это внешнее поле определяется выражением
(5)
Выражение (5) представляет собой статическое решение по координате для электрического поля в поперечном сечении щели с учетом особенностей поведения поля на краях этой щели. При этом внешнее поле (5) имеет характер бегущей волны, распространяющейся вдоль щели по координате к краю полуплоскости.
Приведенные модельные представления для ПЩА уже частично апробировались и проверялись в соответствующих разработках [5; 6; 10–12].
2. Расчет диаграмм направленности исследовавшихся щелевых антенн с применением полученных строгих математических моделей
Полученные в разделе 2 соотношения явились основой для расчета диаграмм направленности исследуемых ПЩА. Расчет обеих компонент поля излучения в дальней зоне в электродинамических плоскостях Е и Н для них производился с использованием выражений (1)–(5) в микроволновом диапазоне, на средней частоте излучения 10 ГГц (соответствующая длина волны = 3 см).
На рис. 2–4 приведены семейства расчетных ДН для ПЩА, у которых длина L принимала значения 3 см, 9 см и 15 см. Для каждого из этих значений длины щели, ее ширина W имела значения 3 см, 6 см и 9 см. На всех рисунках диаграммы направленности нормированы по мощности и показаны половины их изображений в плоскостях Е (справа) и Н (слева), исходя из условия их симметрии.
Рис. 2. Расчетные ДН исследовавшихся ПЩА при постоянной длине L = 3 cм для значений ширины W = 3 cм (кривая 1), W = 6 cм (кривая 2), W = 9 cм (кривая 3)
Fig. 2. Calculated RP of investigated CWSA on constant length L = 3 cm for width values W = 3 cm (graph 1), W = 6 cm (graph 2), W = 9 cm (graph 3)
Рис. 3. Расчетные ДН исследовавшихся ПЩА при постоянной длине L = 9 cм для значений ширины W = 3 cм (кривая 1), W = 6 cм (кривая 2), W = 9 cм (кривая 3)
Fig. 3. Calculated RP of investigated CWSA on constant length L = 9 cm for width values W = 3 cm (graph 1), W = 6 cm (graph 2), W = 9 cm (graph 3)
Рис. 4. Расчетные ДН исследовавшихся ПЩА при постоянной длине L = 15 cм для значений ширины W = 3 cм (кривая 1), W = 6 cм (кривая 2), W = 9 cм (кривая 3)
Fig. 4. Calculated RP of investigated CWSA on constant length L = 15 cm for width values W = 3 cm (graph 1), W = 6 cm (graph 2), W = 9 cm (graph 3)
Из приведенных графиков ДН видно, что с увеличением длины щели L и ее ширины W ширина ДН уменьшается, что свидетельствует об увеличении коэффициента направленного действия. При этом уровень боковых лепестков (УБЛ) ДН снижается, что говорит о повышении усилительных свойства антенны.
Для теоретических (расчетных) ДН бесконечная полуплоскость считается идеально ровной, а щель регулярная и абсолютно симметричная, поэтому влияние кросс-поляризационной составляющей поля в главных электродинамических плоскостях Е и Н отсутствует (показано в разделе 1). Расчет кросс-поляризационной составляющей поля для произвольных значений углов в пределах показал, что она хаотически изменялась относительно среднего уровня порядка –15 дБ, а ее максимальные значения не превышали –8 дБ (в системе координат рис. 2–4 соответственно, уровни 0,032 и 0,158), что для данных ПЩА ниже уровня половинной мощности. Установление более точных характеристик распределения кросс-поляризационной составляющей для ПЩА, подобных исследовавшимся, является задачей достаточно трудоемкой и требует отдельного изучения.
3. Экспериментальное измерение диаграмм направленности исследовавшихся плоских щелевых антенн и их анализ
Для проведения натурного эксперимента был изготовлен ряд опытных образцов плоских симметричных щелевых антенн с постоянной шириной щели, аналогичных показанным на рис. 1. При этом для них длина L и ширина W щели имели те же значения, что и в случае расчетных ДН. Для изготовления опытных ПЩА использовалась металлическая поверхность, имеющая высокую проводимость (медь, алюминий), нанесенная на тонкую диэлектрическую подложку, диэлектрическая проницаемость которой . Поэтому ее наличие практически не оказывает влияния на фазовую скорость волны и позволяет разработку точных математических моделей для данных ПЩА.
Измерение ДН опытных ПЩА производилось также на частоте f0 = 10 ГГц ( = 3 см) на физической установке, описанной в [13; 14], модернизированной для проведения данных исследований. Питание ПЩА производилось аналогично питанию щелевых линий, описанному в [5].
Измеренные диаграммы направленности для исследовавшихся ПЩА приведены на рис. 5–7. Они также нормированы по мощности и совмещены для удобства сравнения.
Рис. 5. Экспериментально измеренные ДН щелевых антенн при постоянной длине L = 3 cм для значений ширины W = 3 cм (кривая 1); W = 6 cм (кривая 2) и W = 9 cм (кривая 3)
Fig. 5. Experimentally measured RP of slot antennas on constant length L = 3 cm for width values W = 3 cm (graph 1); W = 6 cm (graph 2) and W = 9 cm (graph 3)
Рис. 6. Экспериментально измеренные ДН щелевых антенн при постоянной длине L = 9 cм для значений ширины W = 3 cм (кривая 1); W = 6 cм (кривая 2) и W = 9 cм (кривая 3)
Fig. 6. Experimentally measured RP of slot antennas on constant length L = 9 cm for width values W = 3 cm (graph 1); W = 6 cm (graph 2) and W = 9 cm (graph 3)
Рис. 7. Экспериментально измеренные ДН щелевых антенн при постоянной длине L = 15 cм для значений ширины W = 3 cм (кривая 1); W = 6 cм (кривая 2) и W = 9 cм (кривая 3)
Fig. 7. Experimentally measured RP of slot antennas on constant length L = 15 cm for width values W = 3 cm (graph 1); W = 6 cm (graph 2) and W = 9 cm (graph 3)
Из анализа приведенных на рис. 5–7 результатов следует, что для Е-плоскости экспериментально измеренные ДН по уровню половинной мощности очень хорошо согласуются с теорией для значений L = 9 см и 15 cм (когда ). Для более коротких антенн при приближении значений к 1 рассогласование между измеренными и расчетными ДН нарастает и принимает значения, превышающие значение статистической погрешности (по углу ее значение порядка 4°). Экспериментально установлено, что при дальнейшем уменьшении длины ПЩА форма их ДН искажалась, рассогласование с расчетными ДН существенно увеличивалось. Таким образом, для исследовавшихся ПЩА условие определяет границу применимости разработанных модельных представлений.
Из приведенных графиков видно, что для Н-плоскости ширина измеренных ДН по уровню половинной мощности существенно зависит от длины щели и слабо зависит от ее ширины. С увеличением длины щели антенны главный лепесток ДН сужается, что также хорошо согласуется с расчетными данными.
Во всех случаях ширина ДН в плоскости Н получалась несколько шире, чем в плоскости Е, что характерно и для расчетных ДН. Также видно, что у экспериментально полученных ДН вершина главного лепестка более протяженная и пологая, чем в случае расчетных ДН, что многократно подтверждалось при проведении экспериментов.
Эксперимент также показал, что форма ДН исследовавшихся ПЩА в диапазоне частот 8–12 ГГц существенно не изменялась. Уровень боковых лепестков экспериментально измеренных ДН во всех случаях был выше, чем для расчетных ДН, но не превышал значения половинной мощности. Высокий УБЛ, вероятнее всего, объясняется конечной поверхностью исследовавшихся антенн, их неидеальной формой и наличием отражения от стенок измерительной камеры.
Возможность задания формы ДН изменением геометрических размеров щели и удобство размещения подобных ЩА на проводящих поверхностях также позволяют их апробацию и использование в составе более сложных антенных систем, например в составе дисковых антенных решеток [6; 7].
Заключение
Для плоских симметричных щелевых антенн микроволнового диапазона с постоянной шириной щели разработаны математические модели поля излучения в дальней зоне с применением тензорной функции Грина, позволяющие учесть как основную, так и кросс-поляризационную составляющую поля излучения. При этом физическая конструкция ПЩА заменялась идеально проводящей бесконечной полуплоскостью, возбуждаемой внешним сторонним электрическим полем частотой f0 = 10 ГГц (длина волны = 3 см).
С использованием этих моделей теоретически рассчитаны семейства диаграмм направленности исследовавшихся ПЩА в Е- и Н-плоскостях для значений длины щели 3 см, 9 см, 15 см и ее ширины 3 см, 6 см, 9 см. Анализ этих ДН показал, что в Е-плоскости при увеличении длины щели и ее ширины ДН сужаются, что повышает их направленные свойства. При этом уровень боковых лепестков ДН уменьшается, что приводит к увеличению усилительных свойств антенн. В Н-плоскости при изменении длины и ширины щели в тех же пределах ширина главного луча ее ДН по половинной мощности остается практически неизменной и несколько больше, чем для всех ДН в Е-плоскости. Расчет также показал, что для идеализированного случая (бесконечная полуплоскость идеально ровная, щель регулярная и абсолютно симметричная) кроссполяризационная составляющая поля в главных плоскостях отсутствовала. В других плоскостях она имела хаотически изменяющиеся значения относительно среднего уровня –15 дБ, ее максимальные значения не превышали –8 дБ и не достигали уровня половинной мощности излучения.
Результаты эксперимента для аналогичных опытных образцов ПЩА, у которых длина и ширина щели имели те же значения, что и модельные (рис. 2–4), показали, что их ДН имеют те же закономерности, что и расчетные ДН. Сравнительный анализ выявил, что для Е-плоскости экспериментально измеренные ДН по уровню половинной мощности хорошо согласуются с теорией для значений длины щели L, когда отношение . Для более коротких антенн при приближении значений к 1 рассогласование между расчетными и измеренными значениями ширины ДН увеличивается и превышает значение статистической погрешности по углу (порядка 4°). Экспериментально установлено, что при дальнейшем уменьшении длины ЩА форма их ДН искажалась, рассогласование с модельными представлениями существенно увеличивалось, что может свидетельствовать об ограничении возможностей использования разработанных математических моделей. Также установлено, что для Н-плоскости ширина измеренных ДН по уровню половинной мощности хорошо согласуется с расчетными значениями, она тоже определяется длиной щели и слабо (в пределах статистической погрешности) зависит от ее ширины.
Эксперимент также показал, что для исследовавшихся ПЩА форма ДН существенно не изменялась в диапазоне частот 8–12 ГГц. Уровень боковых лепестков экспериментально полученных ДН во всех случаях был выше, чем для расчетных ДН, но не превышал значения половинной мощности. Высокий УБЛ у измеренных ДН, вероятнее всего, объясняется конечной поверхностью исследовавшихся антенн, их неидеальной формой и наличием отражения от стенок измерительной камеры.
Об авторах
Евгений Иванович Нефёдов
Фрязинский филиал Института радиотехники и электроники имени В.А. Котельникова РАН
Email: nefyodov32@mail.ru
1932 г. р., доктор физико-математических наук, профессор, главный научный сотрудник Фрязинского филиала Института радиотехники и электроники имени В.А. Котельникова РАН, г. Фрязино, Московская обл., Россия. Окончил радиотехнический факультет Военной краснознаменной инженерной академии связи имени С.М. Буденного (г. Ленинград) в 1956 г. С 1961 г. по настоящее время работает во Фрязинском филиале Института радиотехники и электроники имени В.А. Котельникова РАН. Академик Академии инженерных наук имени А.М. Прохорова, Академии медико-технических наук, Академии военных наук и Международной академии информатизации. Автор 30 книг, 50 изобретений и 250 статей. Подготовил к защите более 20 кандидатов и 10 докторов наук. Награжден знаками «Почетный радист СССР» и «Почетный радист РФ». Создатель ряда приоритетных научных школ.
Область научных интересов: асимптотическая теория дифракции, системы сверхбыстрой обработки информации на объемных интегральных схемах СВЧ-, КВЧ- и оптических частот, теория единого информационного поля ноосферы.
Россия, 141190, Московская обл., г. Фрязино, пл. Введенского, 1Вячеслав Петрович Заярный
Волгоградский государственный технический университет
Email: zvp2000@mail.ru
1954 г. р., доктор технических наук (2002 г.), профессор, заведующий кафедрой экспериментальной физики Волгоградского государственного технического университета (ВолгГТУ), г. Волгоград, Россия. Окончил радиотехнический факультет Уральского политехнического института (УПИ, г. Свердловск) в 1978 г. по специальности «Радиоэлектронные устройства». До 1984 г. работал в одном из оборонных НИИ, затем поступил в очную аспирантуру (УПИ, физико-технический факультет), которую окончил с защитой кандидатской диссертации (кандидат физико-математических наук, 1988 г.). Академик метрологической академии РФ. Имеет более 200 научных публикаций. Член диссертационного совета Д 212.028.05 ВолгГТУ.
Область научных интересов: электроника (в т. ч. физическая электроника, электроника СВЧ и КВЧ), радиоэлектроника, измерительные преобразователи, датчики, измерительные системы (в т. ч. информационно-измерительные системы), физический эксперимент.
РоссияНина Ивановна Кирьянова
Волгоградский государственный университет
Email: nikiryanova@mail.ru
1949 г. р., кандидат экономических наук, доцент кафедры телекоммуникационных систем института приоритетных технологий Волгоградского государственного университета, г. Волгоград, Россия. Окончила Куйбышевский электротехнический институт связи в 1976 г. и имеет квалификацию инженера автоматической электросвязи. Работала инспектором в Государственном надзоре за связью по Волгоградской области. С 2002 г. преподает в Волгоградском государственном университете. В 2008 г. защитила кандидатскую диссертацию по теме «Развитие инфокоммуникаций в России в условиях перехода к информационному обществу».
Область научных интересов: сети и системы электросвязи, построение телекоммуникационных систем и сетей с применением беспроводных технологий.
Россия, 400062, г. Волгоград, Университетский пр., 100Игорь Николаевич Пономарев
Волгоградский государственный университет
Автор, ответственный за переписку.
Email: ponomarev.igor@volsu.ru
1963 г. р., старший преподаватель Волгоградского государственного университета, г. Волгоград, Россия. Окончил Таганрогский радиотехнический институт имени В.Д. Калмыкова в 1986 г. по специальности «Радиотехника». Окончил заочную аспирантуру под руководством доктора физико-математических наук, профессора Е.И. Нефёдова (ФИРЭ РАН, г. Фрязино). Тема диссертации связана с исследованием плоских антенн бегущей волны СВЧ-диапазона. Имеет более 20 научных публикаций.
Область научных интересов: решение граничных задач электродинамики и математическое моделирование антенн и устройств согласования микроволнового диапазона.
Россия, 400062, г. Волгоград, Университетский пр., 100Список литературы
- Nefyodov E.I., Smolsky S.M. Electromagnetic Fields and Waves. Microwaves and mmWave Engineering with Generalized Macroscopic Electrodynamics. New York: Springer, 2019. 360 p. DOI: https://doi.org/10.1007/978-3-319-90847-2
- Nefyodov E.I., Smolsky S.M. Understanding of Electrodynamics, Radio Wave Propagation and Antennas: Lecture Course for Students and Engineers. Wuhan: Scientific Research Publishing, 2012. 449 p.
- Gibson P.J. The Vivaldi Aerial // 1979 9th European Microwave Conference. 1979. Р. 101–105. DOI: https://doi.org/10.1109/euma.1979.332681
- Широкополосный полосковый излучающий элемент пазового типа для ФАР / Ю.А. Вайнер [и др.] // Антенны. 1980. Вып. 28. С. 95–100.
- Заярный В.П., Парпула С.А., Гирич В.С. Излучение плоских антенн СВЧ-диапазона с линейно расширяющимся раскрывом для дисковых антенных решеток // Журнал технической физики. 2014. Т. 84, № 11. С. 106–111. URL: http://journals.ioffe.ru/articles/41073
- Изучение влияния угла раскрыва плоских коротких щелевых антенн микроволнового диапазона на их электродинамические характеристики / В.П. Заярный [и др.] // Известия вузов. Радиофизика. 2016. Т. 59, № 6. С. 529-534. URL: https://radiophysics.unn.ru/ issues/2016/6/529
- Фролов А.А., Гирич С.В., Заярный В.П. Антенна кругового обзора сверхвысокочастотного диапазона // Известия вузов. Радиофизика. 2012. Т. 55. № 10-11. С. 698–703. URL: https://radiophysics.unn.ru/issues/2012/10/698
- Janaswamy R., Schaubert D. Analysis of the tapered slot antenna // IEEE Transactions on Antennas and Propagation. 1987. Vol. 35, no. 9. P. 1058–1065. DOI: https://doi.org/10.1109/tap.1987.1144218
- Janaswamy R., Shaubert D.H., Pozar D.M. Analysis of the transverse electromagnetic mode linearly tapered slot antenna // Radio Science. 1986. Vol. 21, no. 5. P. 797–804. DOI: https://doi.org/10.1029/rs021i005p00797
- Нефёдов Е.И., Пономарев И.Н. Возбуждение идеально проводящей полуплоскости расширяющейся щелью, прорезанной перпендикулярно ее краю // Электродинамика и техника СВЧ, КВЧ и оптических частот. 2014. Т. 18, № 1. С. 58–84.
- Заярный В.П., Нефедов Е.И., Пономарев И.Н. Особенности электродинамических характеристик симметричных щелевых линий на проводящей бесконечной полуплоскости в микроволновом диапазоне // Известия вузов. Радиофизика. 2018. Т. 61, № 4. С. 315–320. URL: https://radiophysics.unn.ru/issues/2018/4/315
- Виленский А.Р. Метод анализа пространственно-временных характеристик излучения печатных щелевых антенн бегущей волны // Наука и образование. 2014. № 5. С. 139–151. URL: http://technomag.edu.ru/doc/710740.html
- Заярный В.П. Радиофизические свойства твердотельных слоистых структур с зарядовой связью: методы и информационные возможности для их определения. М.: Радио и связь, 2001. 220 c.
- Фролов А.А., Гирич С.В., Заярный В.П. Изучение электродинамических характеристик антенн и антенных систем СВЧ-диапазона // Известия вузов. Радиофизика. 2009. Т. 52, № 4. С. 328–336.