|  | (1) | 
где  ;
;  .
.
В приведенных формулах  ,
,  — углы сферической системы координат, определяющие направление фазирования, (
 — углы сферической системы координат, определяющие направление фазирования, ( ,
,  ) — координаты центров излучателей в декартовой системе координат, ось ОХ которой ориентирована параллельно продольной оси вибраторов, а ось OZ — перпендикулярна плоскости решетки.
) — координаты центров излучателей в декартовой системе координат, ось ОХ которой ориентирована параллельно продольной оси вибраторов, а ось OZ — перпендикулярна плоскости решетки.
Возбуждение (1) позволяет рассматривать излучение элемента решетки как возбуждение электромагнитного поля в пространственном волноводе — канале Флоке [1]. Для анализа поля в пространственном волноводе используются известные соотношения
|  | (2) | 
где  , и u — векторный и скалярный потенциалы.
, и u — векторный и скалярный потенциалы.
Интегральное уравнение задачи получено из граничного условия для касательной составляющей электрического поля на поверхности идеально проводящего излучателя
|  | (3) | 
где  — стороннее электрическое поле, равное нулю всюду, за исключением области зазора между плечами вибратора, где оно принимает известное значение
 — стороннее электрическое поле, равное нулю всюду, за исключением области зазора между плечами вибратора, где оно принимает известное значение  , равное отношению возбуждающего напряжения к ширине зазора. Интегрирование граничного условия (3) с учетом (2) дает
, равное отношению возбуждающего напряжения к ширине зазора. Интегрирование граничного условия (3) с учетом (2) дает
|  | (4) | 
где С — постоянная интегрирования, a t = x, y.
Электродинамические потенциалы представляются в виде разложения по плоским волнам [2], а учет граничных условий для нормальных составляющих тока на кромках вибратора и уравнения непрерывности позволяет перейти от поверхностной плотности тока к поверхностной плотности заряда. Рассматривая волну потенциала в виде суперпозиции волн, связанных с H- и E-волнами пространственного волновода, учитываются границы раздела между диэлектрическими слоями. При этом выражение (4) сводится к интегральному уравнению вида

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

где  — базисная функция;
 — базисная функция;  — неизвестный коэффициент. В качестве базисных функций используются собственные функции области, занимаемой металлической пластинкой (для узкого прямоугольного вибратора — набор тригонометрических функций одной координаты).
 — неизвестный коэффициент. В качестве базисных функций используются собственные функции области, занимаемой металлической пластинкой (для узкого прямоугольного вибратора — набор тригонометрических функций одной координаты).
Используя базисные функции в качестве весовых функций, получаем следующую систему линейных алгебраических уравнений:

где
 
 
Решая полученную систему уравнений, можно найти распределение заряда по металлическому вибратору, через которое находятся все интегральные характеристики излучателя в решетке: парциальная диаграмма направленности, входное сопротивление и другие.
Численные эксперименты проводились при следующих параметрах решетки и излучателя (здесь и далее размеры нормированы к длине волны): шаги решетки  =
 =  = 0,6; ширина вибратора b = 0,035;
 = 0,6; ширина вибратора b = 0,035;  = 0,175;
 = 0,175;  = 2,5;
 = 2,5;  =
 =  = 0°. Предполагалось наличие диэлектрического покрытия толщиной
 = 0°. Предполагалось наличие диэлектрического покрытия толщиной  с диэлектрической проницаемостью
 с диэлектрической проницаемостью  .
.
 
  
 
Рис.1
На рис. 1,а приведены зависимости активной и реактивной составляющих входного сопротивления печатного вибратора в бесконечной решетке от его длины. Параметром является диэлектрическая проницаемость подложки  , равная: 1,5 (кривая 1); 2,0 (кривая 2) и 2,5 (кривая 3). При увеличении
, равная: 1,5 (кривая 1); 2,0 (кривая 2) и 2,5 (кривая 3). При увеличении  происходит уменьшение резонансных размеров вибратора, однако при этом сужается рабочая полоса частот. При уменьшении ширины печатного вибратора b наблюдается увеличение активной составляющей входного сопротивления в области второго резонанса и заметное сужение полосы частот. На рис. 1,б b имеет значения: 0,035 (кривая 1); 0,055 (кривая 2) и 0,075 (кривая 3). Можно также отметить некоторое смещение резонансной длины в сторону больших значений длины вибратора. Аналогичное поведение резонансной длины наблюдается при изменении высоты расположения вибратора над экраном. На рис. 1,в
 происходит уменьшение резонансных размеров вибратора, однако при этом сужается рабочая полоса частот. При уменьшении ширины печатного вибратора b наблюдается увеличение активной составляющей входного сопротивления в области второго резонанса и заметное сужение полосы частот. На рис. 1,б b имеет значения: 0,035 (кривая 1); 0,055 (кривая 2) и 0,075 (кривая 3). Можно также отметить некоторое смещение резонансной длины в сторону больших значений длины вибратора. Аналогичное поведение резонансной длины наблюдается при изменении высоты расположения вибратора над экраном. На рис. 1,в  имеет значения: 0,150 (кривая 1); 0,175 (кривая 2) и 0,2 (кривая 3).
 имеет значения: 0,150 (кривая 1); 0,175 (кривая 2) и 0,2 (кривая 3).
Входное сопротивление вибраторного излучателя в решетке существенно зависит от шага периодической структуры, что является следствием взаимодействия излучателей в ФАР. На рис. 2,а и б представлены зависимости входного сопротивления печатного вибратора шириной b = 0,075 при изменении шага решетки  и
 и  соответственно в Е- и H-плоскостях. Кривые на рис. 2, а соответствуют значениям: 1 —
 соответственно в Е- и H-плоскостях. Кривые на рис. 2, а соответствуют значениям: 1 —  = 0,7; 2 —
 = 0,7; 2 —  = 0,8 и 3 —
 = 0,8 и 3 —  = 0,9, а на рис. 2,б
 = 0,9, а на рис. 2,б  = 0,6 (кривая 1); 0,65 (кривая 2) и 0,7 (кривая 3).
 = 0,6 (кривая 1); 0,65 (кривая 2) и 0,7 (кривая 3).
 
  
 
Рис.2










