ISSN 0236-235X (P)
ISSN 2311-2735 (E)
2

16 Июня 2024

Программный комплекс моделирования диаграммынаправленности фазированной антенной решетки и вычисления кодов фазовращателей

DOI:10.15827/0236-235X.134.301-306
Дата подачи статьи: 17.10.2020
УДК: 004.942

Калинин Т.В. (diiimaaarik@mail.ru) - Военно-космическая академия им. А.Ф. Можайского (профессор), Санкт-Петербург, Россия, кандидат технических наук, Богомаз Д.В. (vka@mil.ru) - Военно-космическая академия им. А.Ф. Можайского (старший преподаватель), Санкт-Петербург, Россия, Хрестинин Д.В. (diiimaaarik@mail.ru) - Военно-космическая академия им. А.Ф. Можайского (преподаватель), Санкт-Петербург, Россия, Вихрев А.В. (vka@mil.ru) - Военно-космическая академия им. А.Ф. Можайского (курсант), Санкт-Петербург, Россия
Ключевые слова: диаграмма направленности, фазовращатель, цифровой код, программный комплекс, среда программирования qt creator
Keywords: directivity diagram, phase shifter, digital codes, software package, qt creator programming environment


     

Сущность фазового управления диаграммой направленности (ДН) по азимуту и углу места одинакова. В качестве примера рассмотрим управление ДН по азимуту.

Пусть линейная эквидистантная решетка с шагом d имеет пять излучателей и необходимо отклонить максимум ДН на угол e от нормали к раскрыву антенны. Следует выяснить фазовые различия между сигналами, излучаемыми элементами такой фазированной антенной решетки (ФАР).

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

На рисунке 1 видно, что путь распространения электромагнитной энергии от первого излучателя до точки Д больше этого же пути от второго излучателя на величину АВ, которую обозначим через ∆r. Тогда за время распространения по пути ∆r фаза сигнала первого излучателя изменится на величину ∆φ, равную  [2–4].

Из прямоугольного треугольника ABC видно, что ∆r = dsine. Тогда

.                                      (1)

Для сложения энергий, излучаемых первым и вторым элементами, необходимо фазу сигнала, излучаемого вторым элементом, увеличить на ∆φe. Достигнуть этого можно с помощью управляемого фазовращателя, на другой вход которого подать код управления фазой ∆φe*.

Аналогичные рассуждения справедливы для любой последующей пары рядом расположенных излучателей (например, для второго и третьего, для третьего и четвертого, для четвертого и пятого). Видно, что для каждого последующего излучателя необходимо увеличить фазу сигнала на ∆φe по сравнению с каждым предыдущим излучателем. Если в решетке имеются i излучателей (i = 0, 1, 2, …), то фазовращатель с номером i должен изменить фазу сигналов [5] на величину :

.                                    (2)

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

Таким образом, для отклонения ДН на угол e от нормали необходимо сформировать коды управления фазовращателями .

Допустим, имеется плоская ФАР, содержа-щая i столбцов и j строк, на пересечении кото-рых установлены элементарные излучатели (рис. 2).

Пусть также геометрический центр этой ФАР совмещен с началом прямоугольной системы координат XYZ. При этом ось OX совпадает с нормалью к раскрыву антенны, OZ – с одной из строк, а OY – с одним из столбцов. Шаг между излучателями в строке равен dx, а в столбце – dy. Каждый из столбцов состоит из nϴ излучателей, а каждая из строк – из ne излучателей. В цепи питания каждого излучателя имеется фазовращатель, на который необходимо подать код управления положением ДН.

Положение ДН в пространстве (направление ОМ) зададим в биплоскостной системе координат с помощью двух углов – азимута e и угла места ϴ. Угол e – плоский угол двугранного угла, образованный плоскостями OXY (плоскость, перпендикулярная к оси OZ, в которой лежит нормаль к плоскости ФАР) и MOY (плоскость ϴe, в которой находится точка М, лежащая на линии ДН). Угол ϴ – плоский угол двугранного угла, образованного плоскостями XOZ (плоскость, перпендикулярная к оси OY, в которой также лежит нормаль к плоскости ФАР) и MOZ (плоскость ϴe, в которой также находится точка М). Для упрощения аппаратуры управления положением ДН в пространстве создается линейно нарастающее фазовое распределение вдоль осей OY и OZ [6, 7].

 

Рис. 2. Размещение излучателей в ФАР

Fig. 2. Positioning of emitters in a phased antenna array
Для отклонения ДН по азимуту на угол e необходимо, чтобы каждый фазовращатель, стоящий в цепи излучателя с номером i, изменял фазу входного сигнала на величину  Для этого на фазовращатель следует подать код управления .

В соответствии с выражениями (1) и (2) можно записать, что .

Тогда величину кода управления  по аналогии с последним выражением можно представить в виде ,  где  – разность величин кодов управления между последующим и предыдущим излучателями в данной строке.

Поскольку необходимо создать линейно нарастающее фазовое распределение, с учетом последнего выражения для каждого фазовращателя, стоящего в данной строке, нужно сформировать следующий код управления [8, 9]:

Максимальное значение кода управления имеет место для последнего излучателя. Обо-значим его через Ne, то есть .

Тогда для текущего фазовращателя i код управления  можно записать соотношением

.                                            (3)

Это справедливо, поскольку

.

Анализ выражений (2) и (3) показывает, что код Ne должен соответствовать коду числа, полученного от произведения .

Аналогично можно показать, что для отклонения ДН в угломестной плоскости на угол ϴ на каждый фазовращатель в столбце, стоящий в цепи излучателя с номером j, необходимо подать код управления , равный

,                                          (4)

где ;  – разность величин кодов управления между последующим и предыдущим фазовращателями в данном столбце [10, 11].

Величина кода Nϴ вырабатывается вычислительным комплексом и должна соответствовать коду числа, полученного от произведения .

Каждый из фазовращателей, расположен-ных в цепи питания вибратора, установленного на пересечении столбца i и строки j, должен участвовать в сканировании ДН как по азимуту, так и по углу места. Тогда величина подаваемого на него кода управления  должна быть равна сумме кодов  и ,  

Рис. 3. Окно диалога с пользователем

Fig. 3. Screenshot of the dialog box

 

Рис. 4. Отклонение диаграммы направленности на –35 градусов

Fig. 4. –35-degree deviation of the radiation pattern
что с учетом выражений (3) и (4) можно пред-ставить в виде  [12, 13].

Компьютерное моделирование описанной выше математической модели реализовано на языке С++ в среде Qt Creator. Программный комплекс представляет собой приложение с графическим интерфейсом (рис. 3).

Пользователю предоставлена возможность изменения следующих исходных данных:

-     количество излучателей;

-     расстояние между излучателями;

-     длина волны;

-     угол отклонения диаграммы направленности от нулевого положения;

-     количество дискрет фазы фазовращателя.

Программа рассчитывает и выводит следующие данные:

-     длина антенны;

-     ширина диаграммы направленности;

-     сдвиги фазы расчетный и дискретный (для определения ошибки установки фазы при использовании дискретных фазовращателей);

-     значение фазы на каждом фазовращателе;

-     значение кода фазы в двоичном коде на каждом фазовращателе.

На осях представлены энергетические характеристики ФАР в линейной и полярной системах координат.

На рисунке 4 представлена реализация вычисления кодов для ФАР из 20 пятиразрядных излучателей при отклонении ДН на –35 градусов.

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

Литература

1.    Гостюхин В.Л., Трусов В.Н., Гостюхин А.В. Активные фазированные антенные решетки. М.: Радиотехника, 2011. 310 с.

2.    Парахин А.С., Крайнюченко Л.В. Компьютерная модель фазированной антенной решетки // Вестн. Курганского гос. ун-та. 2013. С. 63–70.

3.    Хансен Р.С. Фазированные антенные решетки; [пер. с англ.]. М.: Техносфера, 2012. 566 с.

4.    Bhattacharyya A.K. Phased Array Antennas. NJ, Hoboken, John Wiley and Sons Publ., 2006, 516 p.

5.    Monzingo R.A., Miller T.W. Adaptive Arrays. John Wiley and Sons Publ., 1980, 541 p.

6.    Wanner S., Sekar S., Fu B., Weber R.J., Song J.M. Phased array system design and modeling. Proc. ISSSE, 2007, pp. 455–458. DOI: 10.1109/ISSSE.2007.4294511.

7.    Васильев В.С., Ивлев Д.Н. Моделирование пространственных диаграмм направленности антенных систем // Антенны. 2006. № 5. C. 39–44.

8.    Воскресенский Д.И., Канащенков А.И. Активные фазированные антенные решетки. М.: Радиотехника, 2004. 488 с.

9.    Mailloux R.J. Phased Array Antenna Handbook. MA, Norwood, Artech House Publ., 2005, 506 p.

10. Пудовкин А.П., Панасюк Ю.Н., Иванков А.А. Основы теории антенн. Тамбов: Изд-во ГОУ ВПО ТГТУ, 2011. 92 с.

11. Hansen R.C. Phased Array Antennas. John Wiley & Sons, 2009.

12. Benson K. Phased array beamforming ICs simplify antenna design. Analog Dialogue, 2019, vol. 53, no. 1.

13. Легкий Н.М., Унченко И.В. Формирование диаграммы направленности в фазированных антенных решетках // Российский технологический журнал. 2019. Т. 7. № 2. С. 29–38. DOI: 10.32362/2500-316X-2019-7-2-29-38.

References

  1. Gostyukhin V.L., Trusov V.N., Gostyukhin A.V. Active Phased Array Antennas. Moscow, 2011, 310 p. (in Russ.).
  2. Parakhin A.S., Kraynyuchenko L.V. Computer model of the phased array antenna. Vestn. of Kurgan State Univ., 2013, pp. 63–70 (in Russ.).
  3. Hansen R.S. Phased Array Antennas. 2009, 571 p. DOI: 10.1002/9780470529188 (Russ. ed.: Moscow, 2012, 566 p.).
  4. Bhattacharyya A.K. Phased Array Antennas. NJ, Hoboken, John Wiley and Sons Publ., 2006, 516 p.
  5. Monzingo R.A., Miller T.W. Adaptive Arrays. John Wiley and Sons Publ., 1980, 541 p.
  6. Wanner S., Sekar S., Fu B., Weber R.J., Song J.M. Phased array system design and modeling. Proc. ISSSE, 2007, pp. 455–458. DOI: 10.1109/ISSSE.2007.4294511.
  7. Vasilev V.S., Ivlev D.N. Modelling of antennas' spatial patterns. Antennas, 2006, no. 5, pp. 39–44 (in Russ.).
  8. Voskresensky D.I., Kanashchenkov A.I. Active Phased Array Antennas. Moscow, 2004, 488 p. (in Russ.).
  9. Mailloux R.J. Phased Array Antenna Handbook. MA, Norwood, Artech House Publ., 2005, 506 p.
  10. Pudovkin A.P., Panasyuk Yu.N., Ivankov A.A. The Basics of Antenna Theory. Tambov, 2011, 92 p. (in Russ.).
  11. Hansen R.C. Phased Array Antennas. John Wiley & Sons, 2009.
  12. Benson K. Phased array beamforming ICs simplify antenna design. Analog Dialogue, 2019, vol. 53, no. 1.
  13. Legky N.M., Unchenko I.V. Formation of the direction diagram in phased antenna array. Russian Technological Journal, 2019, vol. 7, no. 2, pp. 29–38 (in Russ.). DOI: 10.32362/2500-316X-2019-7-2-29-38.


http://swsys.ru/index.php?page=article&id=4819&lang=%E2%8C%A9=en


Perhaps, you might be interested in the following articles of similar topics: