Выберите количество элементов, частоту операции и азимут и угол возвышения, чтобы управлять основным лучом массива.
N = 11; fc = 28e9; az = 30; el = 20;
Спроектируйте линейный массив на желаемой частоте. Элементом по умолчанию является dipole. Найдите сдвиги фазы, которые будут применяться к каждому элементу линейного массива для управления основным лучом в азимутальном направлении. Обратите внимание на то, что расстояние разделения выбирают равным половине длины волны, чтобы гарантировать отсутствие лепестков решетки.
l = design(linearArray,fc); elem = l.Element; elem.Tilt = 90; l.NumElements = N; figure show(l)
ps_az = phaseShift(l,fc,[az;0]);
Вычислите сдвиги фазы для рулевого управления по повышению. Для этого мы изменяем геометрию нашего начального линейного массива для строения эшелона вдоль z.
elem.Tilt = 90; elem.TiltAxis = [0 1 0]; l.Tilt = 90; l.TiltAxis = [0 1 0]; l.ElementSpacing = 1.05*(elem.Length) ; figure show(l)
ps_el = phaseShift(l,fc,[0;el]);
Создайте N x N прямоугольный массив, содержащую N, 1 X N линейных массивов, сложенных вдоль положительного и отрицательного z-направлений.
l.Tilt = 0; elem.Tilt = 0; l.PhaseShift = ps_az; c = conformalArray; zposn = fliplr((-N+1)/2:1:(N-1)/2); for i = 1:N c.Element{i} = l; c.ElementPosition(i,:) = [0,0,zposn(i)*l.ElementSpacing]; end figure show(c)
Присвойте уровень подрешетки фазы сдвигов и вычислите шаблон
c.PhaseShift = ps_el; figure pattern(c,fc);
figure patternElevation(c,fc,az);
Обратите внимание, что фактическое положение пика отличается от теоретического, вычисленного из-за взаимной связи
Использование возможности установленной антенны позволяет провести первоначальный приблизительный анализ антенной решетки путем включения большой структуры в ее окрестности. В данном примере обеспечивают файл STL большого металлического отражателя, расположенного на четверть длины волны от массива. Анализ обрабатывает массив с помощью полноволнового подхода Method of Moments (MoM), и большой отражатель обрабатывается с помощью приближения Physical Optics (PO).
lambda = physconst('lightspeed')/fc; ref_offset = lambda/4; p = platform; p.FileName = 'GroundPlane.stl'; p.Units = 'm'; p.Tilt = 90; f = installedAntenna; f.Platform = p; f.Element = c.Element; f.ElementPosition = c.ElementPosition; f.ElementPosition(:,2) = ref_offset; f.FeedPhase = ps_el; figure show(f)
figure pattern(f,fc)