CC BY
130
Доклады БГУИР_DoкLADY BGUIR
2017, № 5 (107) 2017, N0. 5 (107)
УДК 621. 372. 512
МЕТОД СИНТЕЗА КОНФОРМНОЙ АНТЕННОЙ РЕШЕТКИ С УЧЕТОМ ПОЛЯРИЗАЦИОННЫХ СВОЙСТВ ИЗЛУЧАТЕЛЕЙ
М О. БУСЕЛ
Военная академия Республики Беларусь, Республика Беларусь
Поступила в редакцию 20 июня 2017
Аннотация. Статья посвящена методу численного синтеза конформных антенных решеток заданной конструкции, учитывающему поляризационные свойства излучающих элементов и произвольное расположение их в раскрыве по требованиям к двум взаимно ортогональным компонентам поля. Эффективность данного метода синтеза подтверждена численным моделированием сферической антенной решетки при расчете амплитудно-фазового распределения в пределах заданной активной области решетки.
Ключевые слова: конформные антенные решетки, синтез, поляризация, кроссполяризация.
Abstract. Article is devoted a method of numerical synthesis conformal array of the set design, considering polarizing properties of radiating elements, and there any arrangement in aperture under requirements to two mutually orthogonal components of a field. Efficiency of the given method of synthesis is confirmed by numerical modeling of a spherical antenna lattice at calculation of peak-phase distribution within the set active area of a lattice.
Keywords: conformal array, synthesis, polarization, cross-polarization.
Doklady BGUIR. 2017, Vol. 107, No. 5, pp. 5-12 Synthesis method of conformal array lattice taking into account polarizing properties of radiators M.O. Busel
Введение и постановка задачи
Конформная антенная решетка (КАР) - это система излучателей, размещенных на выпуклой (конформной) поверхности. Наиболее распространенными конформными поверхностями являются: сфера, цилиндр, конус, эллипсоид и другие. Пространственная ориентация излучателей такова, что направление максимума диаграммы направленности (ДН) каждого из них совпадает с направлением нормали в месте расположения излучателя. КАР позволяют реализовать широкоугольное сканирование лучом без изменения характеристик решетки в секторе сканирования, расширить рабочую полосу частот, совместить излучающее полотно антенной решетки с поверхностью различных объектов и т. д. [1-3].
Задачей синтеза антенной решетки в широком смысле является определение размеров и формы излучающего устройства, законов распределения амплитуды, фазы и поляризации источников излучения по требованиям к векторной ДН, отражающей информацию по двум взаимно ортогональным компонентам поля. Для упрощения решения задачи синтеза по требованиям к векторной ДН считают известными: геометрию (форму) излучающей системы, поляризацию и ДН излучателей. Таким образом, в большинстве случаев задача синтеза сводится к определению амплитудно-фазового распределения (АФР) источников излучения по заданным требованиям к векторной ДН. Элементы плоских и линейных решеток имеют одинаковую ориентацию в пространстве, поэтому поляризацию поля излучения каждого из них можно не учитывать. Задача синтеза для непрерывных и дискретных излучающих линейных и плоских антенн в настоящее время достаточно хорошо проработана [4]. Такие
антенны являются частным случаем КАР. Особенности размещения излучателей, такие как различная ориентация и неэквидистантность при расположении излучателей на поверхности КАР, делают невозможным применение методов синтеза линейных и плоских к таким антеннам. Следовательно, есть необходимость в разработке метода синтеза КАР с учетом поляризационных свойств излучателей.
Цель - разработка метода синтеза КАР, отличающегося возможностью синтезировать КАР с произвольным размещением излучателей, учетом поляризационных свойств излучателей, позволяющего проводить расчет АФР в пределах заданной активной области решетки по требованиям к ДН как по основной, так и по кроссполяризационной (паразитной) составляющим.
Основная часть
В наиболее общем случае задача синтеза по требованиям к векторной ДН формулируется
следующим образом: заданы желаемая форма ДН по двум компонентам поля и конфигурация
антенной решетки, требуется найти распределение комплексных амплитуд возбуждения
излучателей антенной решетки, обеспечивающих воспроизведение требуемой векторной ДН.
Класс КАР наиболее представителен. Поэтому
для рассмотрения возьмем наиболее общую
конструкцию - сферическую поверхность.
излучатель Конструкция сферической антенной решетки (САР)
\\/ с квазигексагональным размещением излучателей
0 \ представлена на рис. 1. Она представляет собой
^— у;__\ систему из N излучателей, расположенных на
\ х поверхности радиуса R. Для управления компонентами
V^,-''' поля необходимо, чтобы излучатель имел два
^С""'^ независимо управляемых входа. В качестве такого
излучателя в общем случае рассмотрим элемент Рис. 1. Конструкция сферической антенной J J r r
решетки Гюйгенса с эллиптической поляризацией и заданными
значениями параметров поляризации [5].
Для численного решения задачи синтеза воспользуемся среднеквадратической мерой
близости требуемой FTP (9, ф) и синтезированной Fs (0, ф) ДН:
0,5л л | г г 2
а2 = J J FTP(9,ф)-F(9,ф) d9dф =аА + аФ +аР ^ min, (1)
-0,5л -л
где оА, оФ, оР - среднеквадратическая ошибка по амплитудной, фазовой и поляризационной ДН соответственно.
Требуемая векторная комплексная ДН описывается следующим выражением:
4(9,ф) = Атр(9,ф)фтр(9,ф) Ртр(9,ф), (2)
где АТР (9, ф), ФТР(9, ф), рТР(9, ф) - требуемая векторная комплексная ДН, включающая
в себя амплитудную, фазовую и поляризационную диаграммы соответственно. Выражение (2) в классической теории антенн записывается:
FTp(9, ф) = FTp(9, ф) Ртр(9, ф) , (3)
где FTP (9, ф) - требуемая комплексная ДН.
Требуемая поляризационная ДН описывается следующим выражением:
Ртр(9, ф) = ртр(9, ф)9о + ртр(9, ф)фо , (4)
где рТР (9, ф)9о, рфР (9, ф)ф0 - проекции поляризационного вектора на 9-ю и ф-ю компоненты
соответственно, причем |рТР(9, ф)| = 1 [4]; 9 о и ф0 - орты сферической системы координат.
Имея две компоненты поля, будем одну из них считать основной, например, 9-ю компоненту, а вторую, ф-ю компоненту - кроссполяризационной. Причем рТР(9, ф) = рТрн(9, ф) ,
(0, ф) = рТросс (0, ф) . Основная компонента рТрН (6, ф) = а кроссполяризационная -^кросс(9,ф). На основании этого выражение (4) примет вид:
рр(0,ф) = 00 + Ккросс(0,ф)ф0 . (5)
Выражение для векторной комплексной ДН F (0, ф) запишем в следующем виде:
_ i N
F(0, ф)=N ZAn Fin(0, фУ 0-(0,ф), (6)
N n=1
где An - АФР n-го излучателя; F1n (0, ф) - векторная ДН n-го излучателя [5], Оп (9,ф) -коэффициент, учитывающий кривизну конформной поверхности от n-го излучателя: 2 %
On (0, ф) = — (xn sin 0 cos ф+yn sin 0 sin ф + zn cos 0), (7)
X
где X - длина волны; xn, yn, zn, - декартовы координаты n-го излучателя.
Векторная комплексная ДН n-го излучателя имеет следующий вид [5]:
4 (0, ф) = 4n (0, ф)0о + ^ (0, ф)фо, (8)
где i*i0 n (0, ф) ( n (0, ф) ) - ДН n-го излучателя по 9-й (ф-й) компонентам поля соответственно:
4- (0, ф) = A Р0 (0, ф) FB- (0, ф)+ A Р0 (0, фЖГ- (0, ф) ; (9)
(0, ф) = Ain Рф (0, ф) FiBn (0, ф)+ Ain Рф (0, ф) Flip n (0, ф). (10)
Выражения (9)-(10) с учетом (3) и (5) запишем в матричном виде:
"FBn (0,ф) С (0, ф )
[1 Ккросс(0,ф)] = [Р0 (0,ф ) рф(0,ф)]
1 n
Fi n (0,ф) < n (0, ф)
(11)
10n V ' т J 1 фп
Коэффициенты Р0 (0,ф) и Рф (0,ф) находятся из решения уравнения (11):
[Р0(0,ф) Рф(0,ф)] = [1 Ккросс(0,ф)]
F1B n (0,ф) С (0,ф) n (0,ф) Fj n (0,ф)
(12)
Коэффициенты Р0 (0,ф) и Рф (0,ф) имеют следующий вид:
Ре (0ф) = С (0,ф)__^Ч0,ф)(0,ф) . (13)
- FB- (0,ф) F1in (0,ф) - Fn (0,ф) FB- (0,ф) С (0,ф) F,- (0,ф) - Fn (0,ф) FBn (0,ф) '
_„ К кросс(е,Ф)^1В (9,ф)__^ (9,ф)_ (14)
р* ^ (9,ф) (9,ф) - ^ (9,ф) (9,ф) ^ (9,ф) (9,ф) - Яфя (9,ф) ^ (9,ф) '
Синтезированная диаграмма направленности (9, ф) имеет вид:
- 1 N
Ъ (9,Ф) = N (9,ф) 0п (9,Ф), (15)
п=1
где Ап = АП + АП - АФР, включающее в себя АФР горизонтального и вертикального входов
п-го излучателя соответственно; /1п (9,ф) - векторная ДН п-го излучателя [5].
Представим АФР для горизонтального и вертикального входов на раскрыве КАР в виде разложений в бесконечные ряды Фурье:
ад ад
£ = ЁЁ йл* [9п -90' Фп - Фо ]+ п
a- = ZZ К Shs [0n -00, ф- - ф0 ]+bB q. [0n -00, ф- - ф0 ]} ; (16)
s=0 l=0 ад ад
A = ZZ К. SUs [0n -00, ф- - ф0 ] + b[s C,s [0n -00, ф- - ф0 ]} , (17)
-1
где 2 ( 2) и Ь в ( Ь в) - комплексные коэффициенты разложения (в самом общем случае в бесконечный ряд Фурье) для вертикального (горизонтального) входа; 51, 2 к -0О. фи "фо ]= sin [1 (0и -00 )]sin ифи "фо )] и С¡,s [ея-00. фи "фо ]= ^фи -00 )]с0ф(фи "фо )] -функции разложения.
Далее по тексту в целях сокращения записи вместо 5 [0и -0О, Фи -Ф0] и
С1,2 [еи -е0. Фи -Ф0 ] будет использоваться , 2 [еи , Фи ] и С1,2 [еи , Фи ] .
Подставим (16) и (17) в (15). Тогда выражение для синтезированной векторной комплексной ДН КАР можно записать следующим образом:
~ 1 N 5 Ь , .
^ (е, ф) = N £££{« + с) [0и, фи и Ъ1+¿Г) ^ [0и, фи ]} (е, фу 0(0,ф), (18)
™ и=1 5=0 1=0
где Ь = 5 - конечное число членов ряда Фурье, характеризует точность разложения двух функций [6].
(Г * \\
Ь = еМ
0,5
2л-1
2Ф0,5
( ( ( 2л ^
5 = ем 0,5 1
20 0 5
V V 0,5 JJJ
(19)
V V ' J J
Введем новое обозначение и представим выражение (18) как
Ае, Ф) = Е Е {(¿В + С)(¿¿В, (0, Ф) + К (0, Ф)) + (й£ + ¿1)(¿сВ^ (0, Ф) + ¿с?, (0, Ф))}, (20)
х=0 1=0
где ¿^ (0, ф) (¿^ (0, Ф)), ¿сс^ (0, ф) (¿с^ (0, Ф)) - базисные функции для вертикального
(горизонтального) входа, которые описываются следующими соотношениями: 1 N
К(0,Ф) = ^(0,Ф)(Ри0 (0,Ф) + рф(0,Ф)) ^ [0и,Фи] е~ю"(0,ф)} ; (21)
™ и=1
1 N
К (0, Ф) = - Е ЙВи (0, Ф) (Ри0 (0, ф) + РФ (0, Ф)) С^ [0и, Фи ] е-0и(0,ф)} ; (22)
™ и=1
1 N
К (0, Ф) = ^ Е {^1 (0, Ф) (Ри0 (0, ф) + Риф (0, ф)) 51,. [0и, Фи ] е-°"(0,ф)} ; (23)
и=1
1 N
К (0, Ф) = „ Е (0, Ф) (Ри0 (0, ф) + Риф (0, ф)) С^ [0и, Фи ] е-0и(0,ф)} . (24)
™ и=1
Как видно из выражения (16), суть предлагаемого метода заключается в представлении заданной ДН в виде разложения в ряд Фурье по системе базисных функций ¿sfs (0, ф)
(¿s¡s (0, ф)), ¿с^ (0, ф) (¿с1;. (0, ф)), нахождении комплексных коэффициентов
¿С (¿1)
и bfs (Ь^) с дальнейшей подстановкой их в формулу (15) для расчета комплексных амплитуд возбуждения излучателей КАР.
Из соотношений (21)-(24) видно, что функции ¿sfs (0, ф) (¿sJs (0, ф)), ¿с^ (0, ф)
(¿с^ (0, ф)) содержат в себе информацию как о конструктивных особенностях КАР, так
и о направленных и поляризационных свойствах ее элементов.
Коэффициенты разложения а^ (¿¡г\) и bfs (о\\) в самом общем случае находятся с использованием стандартной процедуры [7]:
~ 0,5л л ^
¿В, = — [ [ ^ТР(0, ф) ¿з^ (0, феdф ; (25)
' тт * * '
V2 ТР ^ Ч^ ^1
л Д
-0,5 л -л
/■ч 0,5л л В 2 Г Г п т
¿¡В = т [ [ ^тр(0, Ф) ¿сВ^ (0, ф)dеdф; (26)
л Д -
-0,5л -л
. г 2
аи = —
71
0,5 л л
I I /ТР(0, ф) Zsls (0, еdф ;
-0,5 л - л 0,5 л л
=— I I ^(е, ф) ^ (0, ф)^е^ф.
л 05 -
-0,5л -л
(27)
(28)
Следует отметить, что ^ТР(0, ф) в ряде случаев можно получить, задавая требуемые направленные свойства в двух главных плоскостях [8]:
4(0, ф) = 4р(0) 4Р(ф), (29)
где /¡0 (0) и /ф (ф) - требуемые векторные комплексные ДН в соответствующих главных плоскостях.
Выражения (1)-(29) описывают метод численного синтеза АФР КАР заданной конструкции по требованиям к двум взаимно ортогональным компонентам поля.
Пример
В качестве примера синтезирована сферическая антенная решетка с квазигексагональным размещением излучающих элементов (А/ ~ 0,6 X), радиус принят равным Я = 9Х с равномерным амплитудным распределением и фазовым распределением, компенсирующим набеги фазы для получения плоского фронта волны (рис. 2). Для САР задавалась активная область, ограниченная телесным углом - 20обл = 130° и 2фобл = 130°.
б
Рис. 2. Объемные ДН САР при равномерном АФР с учетом кривизны решетки: а - ДН по 0-й компоненте (основная компонента); б - ДН по ф -й компоненте (кроссполяризационная компонента)
Представленная на рис. 2 ДН по кроссполяризационной компоненте « -22 дБ с графической точностью воспроизводит кроссполяризационную ДН, представленную в [1]. Представленные на рис. 2 ДН показывают значительный уровень паразитной компоненты. Лепестки кроссполяризационной ДН с максимальным уровнем «-22 дБ находятся в главном лепестке ДН по основной компоненте. Диаграмма направленности по основной компоненте имеет уровень бокового излучения «-17 дБ, что связано с пространственным разворотом излучателей решетки.
Для проверки работоспособности предложенного метода синтеза КАР была синтезирована сферическая антенная решетка с аналогичной конструкцией. Для основной ДН (по 0-й компоненте) задавались следующие параметры: FБ =FфБ = F0Б = - 25 дБ; 20о,5 = 4°; 2ф0,5 = 4°. Кроссовая ДН не должна превышать уровень, задаваемый коэффициентом ^кросс(0,ф)= - 35 дБ.
Результаты численного синтеза САР приведены на рис. 3, 4. Из этих рисунков видно, что требования по уровню боковых лепестков FБ = - 25 дБ и ширине главного лепестка ДН в двух плоскостях (20о,5 = 4°, 2ф0,5 = 4°), синтезируемых по основной и результирующей ДН, а также требование по уровню кроссовой компоненты САР практически полностью выполнены. Это подтверждает работоспособность предложенного метода синтеза САР с учетом поляризационных свойств излучателей. Отклонения требуемых и синтезируемых ДН можно объяснить неравномерностью разбиения поверхности САР на квазигексагональную сетку при А! « 0,6 X, компенсацией кроссовой компоненты, ошибками при округлении результатов вычисления [9] и др.
б
Рис. 3. Объемные ДН САР с синтезированным АФР: а - ДН по 0-й компоненте (основная компонента); б - ДН по ф-й компоненте (кроссполяризационная компонента)
,дБк
\
Й.Л
\АЛ
м 1М
ч / 1,
Л -Ж
10 120 130 140
ш
АЛ
ТАГ
50 160 170
|\Л т
ш
АА
л л
Ш
О 10 20 3 0 40 5 0 60 70 80 90
П0,,ф)|, дБ
90 100 110 120 130 140 150 160 170
О 10 20 30 40 50 60 70 80 90
дВ
Т
Л Г
90 100 110 120 130 140 150 160 170
Л А.
А'-.А
АА
АА
Ч А
А/1
О 10 20 30 40 50 60 70
д е
Рис. 4. ДН САР: а, в, д - плоскость 0 F0(0, ф0), Fф(0, ф0), F(0, ф0) соответственно; б, г, е - плоскость ф ^(0О, ф), ^(0О, ф), ^(0О, ф) соответственно
е
ф
б
а
ф
е
в
г
е
ф
Заключение
Разработан метод синтеза КАР, отличающийся возможностью синтезировать КАР с произвольным размещением излучателей, учетом поляризационных свойств излучателей, позволяющий за счет этого обеспечить снижение уровня кроссполяризационного излучения в направлении главного лепестка ~ 25 дБ.
При моделировании равномерного амплитудного распределения и фазового распределения, учитывающего кривизну сферической поверхности на раскрыве антенной решетки, сформирована ДН на кроссполяризационной составляющей поля, которая с графической точностью воспроизводит ДН, представленную в [1].
Метод синтеза основан на представлении искомого АФР на раскрыве в виде двумерного ряда Фурье и вычислении коэффициентов для каждого из входов излучателя. На основании полученных выражений (1)-(29) рассчитывается АФР на 2^-входах КАР.
Результаты численного моделирования на примере синтеза сферической антенной решетки по расчету АФР (т. е. решения задачи синтеза) в пределах заданной активной области САР с учетом поляризационных свойств излучателей подтверждают высокую эффективность данного метода.
Список литературы
1. Josefsson L., Persson P. Conformal array antenna theory and design. New Jersey: John Wiley & Sons, Inc., Hoboken, 2006. 488 p.
2. Хансен Р.С. Фазированные антенные решетки. М.: Техносфера, 2012. 560 с.
3. Воскресенский Д.И., Пономарёв Л.И., Филиппов В.С. Выпуклые сканирующие антенны. М.: Радио и связь, 1978. 306 с.
4. Марков Л.Н. Антенные системы радиоэлектронной техники. М.: Воениздат. 1993. 368 с.
5. Бусел М.О., Калинин А.А., Романович А.Г. Математическая модель излучателя конформной антенной решетки с учетом его поляризационных свойств // Сб. научных статей Воен. акад. Респ. Бел. 2016. № 30. С. 86-95.
6. Калинин А.А., Романович А.Г. Метод численного синтеза чебышевских кольцевых антенных решеток // Доклады БГУИР. 2010. № 7 (53). С. 5-11.
7. Анго А. Математика для электро- и радиоинженеров. М.: Наука, 1964. 772 с.
8. Заневский Д.В., Носков В.В. Аппроксимация двумерной диаграммы направленности радиосредства по ее сечениям в горизонтальной и вертикальной плоскостях // Вест. Воен. акад. Респ. Бел. 2008. № 3. С. 25-28.
9. Романович А.Г., Калинин А.А. Влияние неидентичностей излучающих модулей на характеристики и параметры кольцевых антенных решеток // Вест. Воен. акад. Респ. Бел. 2010. № 4. С. 66-72.
References
1. Josefsson L., Persson P. Conformal array antenna theory and design. New Jersey: John Wiley & Sons, Inc., Hoboken, 2006. 488 p.
2. Hansen R.S. Fazirovannye antennye reshetki. M.: Tehnosfera, 2012. 560 s. (in Russ.)
3. Voskresenskij D.I., Ponomarjov L.I., Filippov V.S. Vypuklye skanirujushhie antenny. M.: Radio i svjaz', 1978. 306 s. (in Russ.)
4. Markov L.N. Antennye sistemy radiojelektronnoj tehniki. M.: Voenizdat. 1993. 368 s. (in Russ.)
5. Busel M.O., Kalinin A.A., Romanovich A.G. Matematicheskaja model' izluchatelja konformnoj antennoj reshetki s uchetom ego poljarizacionnyh svojstv // Sb. nauchnyh statej Voen. akad. Resp. Bel. 2016. № 30. S. 86-95. (in Russ.)
6. Kalinin A.A., Romanovich A.G. Metod chislennogo sinteza chebyshevskih kol'cevyh antennyh reshetok // Doklady BGUIR. 2010. № 7 (53). S. 5-11. (in Russ.)
7. Ango A. Matematika dlja jelektro- i radioinzhenerov. M.: Nauka, 1964. 772 c. (in Russ.)
8. Zanevskij D.V., Noskov V.V. Approksimacija dvumernoj diagrammy napravlennosti radiosredstva po ee sechenijam v gorizontal'noj i vertikal'noj ploskostjah // Vest. Voen. akad. Resp. Bel. 2008. № 3. S. 25-28. (in Russ.)
9. Romanovich A.G., Kalinin A.A. Vlijanie neidentichnostej izluchajushhih modulej na harakteristiki i parametry kol'cevyh antennyh reshetok // Vest. Voen. akad. Resp. Bel. 2010. № 4. S. 66-72. (in Russ.)
Сведения об авторе
Бусел М.О., м.т.н., адъюнкт кафедры радиотехники Военной академии Республики Беларусь.
Information about the author
Busel M.O., master of technical science, adjunct the engineer of radioengineering department of Military academy of Republic of Belarus.
Address for correspondence
Адрес для корреспонденции
220057, Республика Беларусь, г. Минск, пр. Независимости, д. 220, Военная академия Республики Беларусь тел. +375-29-571-76-88; e-mail: mischabox@yandex.by Бусел Михаил Олегович
220057, Republic of Belarus, Minsk, Nezavisimosty ave., 220, Military academy of Republic of Belarus tel. +375-29-571-76-88; e-mail: mischabox@yandex.by Busel Mikhail Olegovitch
