При расчетах в технике высоких частот с применением зеркальных отражающих систем (параболических зеркал) всегда возникает задача поиска фазового центра антенны (ФЦА), т.к. правильная работа зеркала возможна только, если в фокусе находится антенна (именуемая облучатель, feeder, feedhorn) которая имеет фазовый фронт волны в виде сферы, и центр этой сферы находится в фокусе зеркала. При любых отклонениях, как формы фазового фронта от сферы, так и смещения ФЦА из фокуса зеркала — КПД зеркальной системы падает, потому что искажается её диаграмма направленности.
Хотя тема поиска ФЦА довольно актуальна даже в быту, потому что кроме традиционных антенн спутникового телевидения нашли распространение параболические антенны для WiFi, WiMAX и сотовой связи (UMTS/3G, LTE/4G) — тем не менее в литературе такая тема освещена слабо и пользователи часто путают фазовую ДН с обычной диаграммой направленности.
В видеороликах о программах компьютерной симуляции иногда можно встретить практические инструкции как искать ФЦА, но обычно там нет даже минимальных объяснений что ищем и что получаем.
Поэтому чтобы восполнить пробел, напишем небольшую статью с практическими примерами.
Фазовая диаграмма направленности – это зависимость фазы электромагнитного поля, излучаемого антенной, от угловых координат.
(А.П. Пудовкин, Ю.Н. Панасюк, А.А. Иванков — Основы теории антенн)
Так как в дальней зоне антенны векторы поля Е и Н синфазны, то и фазовая ДН в одинаковой степени относится к электрической и магнитной составляющей ЭМП, излучаемого антенной.
Обозначается фазовая ДН греческой буквой Пси:
Ψ = Ψ (θ, φ), при r = const.
Если Ψ(θ, φ) = const при r = const, то это означает, что антенна формирует фазовый фронт волны в виде сферы.
Центр этой сферы, в котором находится начало системы координат, называют фазовым центром антенны (ФЦА).
Фазовый центр антенны — это точка, в которую можно поместить одиночный излучатель сферической волны, эквивалентный рассматриваемой антенной системе в отношении фазы создаваемого поля.
(Драбкин А.Л., Зузенко В.Л. Антенно-фидерные устройства)
ФЦА имеют не все антенны. У антенн, имеющих фазовый центр и многолепестковую амплитудную ДН с чёткими нулями между ними, фаза поля в соседних лепестках отличается на π (180°).
Взаимосвязь между амплитудной и фазовой диаграммами направленности одной и той же антенны иллюстрируется
В реальных антеннах фазовый центр обычно рассматривается в рамках ограниченных углов главного лепестка диаграммы направленности. Положение фазового центра зависит от частоты используемого сигнала, направления излучения/приема антенны, его поляризации и других факторов. Некоторые антенны не имеют фазового центра в общепринятом понимании.
В простейших случаях, например у параболической антенны, фазовый центр совпадает с фокусом параболоида и может быть определен из геометрических соображений. В более сложных случаях, например, рупорных антенн, положение фазового центра не очевидно и требует соответствующих измерений.
Натурные измерения фазового центра очень трудоёмкие (особенно в широкой полосе частот).
В CAD-симуляторах электромагнитных полей вычисление ФЦА это очень простая задача, но она всё же требует несколько ручных манипуляций, т.к. выполняется она «брут-форсом» и требует небольшой начальной настройки функции, которую собираемся брутфорсить.
Для практических расчетов возьмем реальный облучатель парабол для диапазона Ku-band — LNB производителя Inverto, серия Black Ultra.
Этот фидер имеет такой вид (в разрезе)
Шарик размером с горошину — это и будет ФЦА, но мы ещё этого не знаем и наша задача найти его положение.
В примере будем использовать такие вводные:
— частота расчета 11538.5 МГц (длина волны 25.982 мм)
— линейная горизонтальная поляризация (в оси Y)
— сама антенна направлена по оси X, т.е. главное направление излучения θ=90, φ=0
Расчет традиционных параметров дальнего поля (Far Field) в Ansys HFSS дает такую диаграмму направленности в 3D и 2D
Мгновенные значения напряженности (Вольт/метр) электрического поля (E-field) в зависимости от фазы
Интегральная напряженность E-поля (за >1 оборот волны)
Все такие параметры дальнего поля (Far-Field) как при натурных измерениях, так и в CAD-симуляциях — рассчитывается на бесконечной сфере — Infinite Sphere. Испытуемая антенна или её компьютерная модель помещается в центр такой сферы, а измерительный зонд двигается по периметру такой сферы и измеряет амплитуду, поляризацию (амплитуду одного из компонентов) и фазу ЭМ волны. Зонд можно закрепить стационарно и поворачивать испытуемую антенну.
Главное чтобы:
— расстояние всегда было одинаковое (т.е. это была именно измерительная сфера)
— радиус сферы был достаточно большой, чтобы измерения проводились только в той области пространства где векторы электрического поля Е и магнитного H синфазны, т.е. ни одна из компонент не преобладает и не смещена по фазе (не имеет реактивности) за счет носителей заряда которые есть в металлических проводниках антенны или за счет заряженных молекул диэлектрика.
В Ansys HFSS для проведения измерений дальнего поля необходимо создать хотя бы одну бесконечную сферу: Radiation -> Insert Far Field Setup -> Infinite Sphere
φ и θ можно всегда указывать от 0 до 360, но чтобы экономить время на вычислениях иногда рационально ограничить исследуемый угол некоторым сектором. При задании шага 1 градус, полная сфера будет занимать 360*360 = 129 600 расчетных точек, а при шаге 0.1 градус почти 13 млн. Для создания 3D/2D отчетов диаграммы направленности обычно достаточно шага 2-3 градуса (14 400 расчетных точек при шаге 3 градуса). Шаг 1 градус и менее есть смысл использовать только для анализа среза
В закладке «Coordinate System» каждая сфера обязательно имеет свой центр координат. По умолчанию там всегда стоит глобальный центр координат проекта [0, 0, 0]. При желании можно добавлять любое количество других относительных координат. Как элементы геометрии модели так и пользовательскую сферу «Infinite Sphere» можно назначать относительно глобального центра координат или относительно пользовательского. Этим мы воспользуемся ниже.
Расходящийся фазовый фронт волны было видно на анимации Е-поля выше. ЭМ волна образует концентрические круги, подобные кругам на воде от брошенного камня. Фазовый центр это точка, в которую бросили такой камень. Видно что его положение находится где-то в раструбе рупора, но точное его положение не очевидно.
Метод поиска ФЦА базируется на том, что мы смотрим на направление вектора Е-поля (его фазу) по поверхности бесконечно удаленной сферы.
Для демонстрации создадим 2 анимации с векторами Е-поля на сфере с радиусом 4 лямбда (это не бесконечная сфера, но для лучшего масштаба рисунка такого радиуса вполне достаточно).
На первой анимации центр сферы размещен точно в ФЦА
На второй анимации центр размещен в точке проекта 0, 0, 0 (забегая наперед скажем что она находится на 25.06 мм позади ФЦА)
На поверхности первой сферы (она кривая, это не плоскость) видно что векторы движутся синхронно. Амплитуда (magnitude) их разная, потому что ДН антенны имеет максимум в центре (до 14.4 dBi) который плавно угасает в 2 раза (-3 dB) при углах ±20°.
Нас интересует не цвет/длина, а направление вектора. Чтобы все они двигались синхронно (синфазно).
На первой анимации все векторы двигаются синхронно, как бы мяч вращается то вправо то влево.
На второй анимации векторы несинхронны, одни уже изменили направление движения, другие ещё нет. Поверхность этой сферы постоянно претерпевает поверхностные натяжения/деформации.
Первая сфера — находится в ФЦА, вторая не находится в ФЦА.
Задача поиска ФЦА по этому методу состоит в том, чтобы с мелким шагом двигать (брутфорсить) Infinite Sphere до тех пор, пока разброс фаз на интересующем нас участке этой сферы (нас интересует только главный лепесток излучения) станет минимальным (в идеале — нулевой).
Но перед тем как перейти к брутфорсу, сначала разберемся как в HFSS можно отобразить фазовую ДН.
В отчетах дальнего поля «Results -> Create Far Field Report» мы можем вывести или традиционный прямоугольный график (Rectangular plot) или 2D круговой график (Radiation pattern) где по одной оси (например X) вывести зависимость угловой координаты (например θ), а по оси Y — значения фазы на этих углах θ.
Нужный нам отчет это rE — «излученное (radiated) E поле».
Для каждого угла [φ, θ] на бесконечной сфере рассчитывается комплексное число (вектор) электрического поля.
При построении обычных амплитудных графиков (диаграмма направленности, распределение мощности излучения по направлению) нас интересует амплитуда (mag) этого поля, которую можно получить или как mag(rE) или сразу используя более удобную переменную Gain (мощность приведена относительно мощности на порте возбуждения и относительно изотропного излучателя).
При построении фазовой ДН нас интересует мнимая часть комплексного числа (фаза вектора) в полярной нотации (в градусах). Для этого используется математическая функция ang_deg (угол_в_градусах) или cang_deg (накопленный_угол_в_градусах)
Для антенны LNA Inverto Black Ultra, фазовая ДН в плоскости XZ (φ=0) при горизонтальной поляризации возбуждения (rEY) имеет такой вид
Угол Theta=90 это излучение вперёд, Theta=0 вверх, Theta=180 вниз.
Значения ang_deg изменяются от -180 до +180, угол 181° это угол -179°, поэтому график имеет форму пилы при проходе через точки ±180°.
Значения cang_deg накапливаются если направление изменения фазы постоянно. Если фаза сделала до 3 полных оборотов (6 раз пересекла 180°) то накопленное значение достигает 1070°.
Как было написано в начале статьи, фазовая и амплитудная ДН у антенн обычно связаны одна с другой. В соседних амплитудных лепестках (beam) фазы отличаются на 180°.
Наложим один на другой графики фазовой (красный/салатовый) и амплитудной (фиолетовый) ДН
Горбы на амплитудной ДН четко следуют переломам фазы, как и написано в книгах.
Нас интересует фазовый фронт только в определенном секторе пространства, в пределах главного лепестка излучения (остальные лепестки всё равно светят мимо параболического зеркала).
Поэтому ограничим график только сектором 90 ±45° (45-135°).
Добавим на график маркеры MIN (m1) и MAX(m2) которые показывают наибольший разброс фаз в исследуемом секторе.
Кроме того добавим математическую функцию pk2pk() которая автоматически ищет на всём графике минимум и максимум и показывает разницу.
На графике выше разница m2-m1=pk2pk= 3.839 °
Задача поиска ФЦА состоит в том, чтобы двигать с мелким шагом Infinite Sphere пока значение функции pk2pk(cang_deg(rE)) не минимизируется.
Для передвижения Infinite Sphere необходимо создать ещё одну дополнительную систему координат:
Modeler -> Coordinate System -> Create -> Relative CS -> Offset
так как мы заведомо знаем, что у симметричного рупора ФЦА будет находиться на оси Х (Z=Y=0), то для Z и Y ставим 0, а двигать будет только вдоль оси X, для чего присвоим переменную Pos (с начальным значением 0 мм)
Чтобы автоматизировать процесс брут-форса, создадим задачу на оптимизацию.
Optimetrics -> Add -> Parametric
, и зададим шаг переменной Pos 1 мм, в диапазоне от 0 до 100 мм
В закладке “
Calculations -> Setup Calculation
” выберем тип отчета «Far Field» и функцию pk2pk(cang_deg(rEY)). В кнопке «Range Functions» укажем диапазон от -45 до +45 градусов (или любой другой интересующий)
Запускаем
ParametricSetup1 -> Analyze
.
Расчет выполняется достаточно быстро, т.к. все расчеты дальнего поля относятся к Post-Processing и не требуют повторного решения модели.
После завершения расчета нажимаем
ParametricSetup1 -> View analysis results
.
Видим четкий минимум при расстоянии X=25mm
Для более высокой точности редактируем параметрический анализ в диапазон 25.0-25.1 мм с шагом 0.01 мм
Получаем четкий минимум на X=25.06 мм
Чтобы визуально оценить где в модели получился ФЦА, можно нарисовать сферы (Non-model) или точки.
Вот в точку X=25.06 мм помещены 2 сферы (радиусом 2 и 4 лямбда)
Вот то же, в анимации
Вот более крупным планом нарисована плоскость и горошина в точке X=25.06
Широко распространено ошибочное мнение, что в HFSS (и других программах, например CST) при наложении графика «3D Plot» на геометрию антенны такой график автоматически помещается в ФЦА.
К сожалению это не так. График 3D всегда накладывается в центр координатной системы, которая были использована при задании «Infinite Sphere» для этого графика. Если использовалась глобальная система координат по умолчанию [0, 0, 0], то 3D Plot будет размещен в 0,0,0 (даже если сама антенна находится далеко в стороне).
Чтобы совместить графики, в настройках 3D Plot необходимо выбрать такую «Infinite Sphere» (создать ещё одну), для которой задана «Relative CS» в точке ФЦА которую мы нашли вручную.
Следует отметить, что такое наложение будет правдиво только для исследуемого сектора (например главного луча ДН), в боковых и задних лепестках ФЦ может находиться в другом месте или быть несферичным.
Также отметим, что настройки «Infinite Sphere» не имеют никакого отношения к граничному условию «Radiation Boundary». Слой Rad можно задавать как прямоугольник, конус, цилиндр, шар, элипсоид вращения и как угодно двигать его положение, форму и поворот. Положение и форма «Infinite Sphere» от этого никак не изменится. Это всегда будет сфера (шар) с бесконечным (достаточно большим) радиусом и с центром в заданной координатной системе.
Файл модели LNB_InvertoBlackUltra.aedt для изучения доступен по ссылке: https://goo.gl/RzuWxW (Google Drive). Для открытия файла требуется Ansys Electronics Desktop v19 или выше (не ниже 2018.1)
Основные теоремы направленности.
/1/
с. 49-54
Фазовый центр антенны.
ХН антенны является
комплексной функцией
Существуют антенны,
для которых можно найти такую точку на
активной поверхности, что ели поместить
в нее центр сферы большого радиуса и
записать ХН для этого начала отсчета,
то окажется, что arg
либо не меняется, либо изменяется скачком
на
при изменении направления. Такая точка
называется фазовым центром антенны.
При вращении
антенны вокруг фазового центра фаза
сигнала в точке наблюдения либо совсем
не меняется, либо претерпевает скачки
на
.
Пример: В
линейной антенне с нулевым фазовым
распределением фазовый центр совпадает
с геометрическим центром А.
Если ХН можно
представить в виде произведения некоторой
вещественной функции на
, то антенна имеет фазовый центр
расположенный в точке, положение которой
относительно начала координат
характеризуется радиусом вектором
.
Если же антенна
не имеет фазового центра, то пользуются
эквивалентным фазовым центром, называя
им центр сферы наименее отличающий ее
от фазового центра.
Теорема умножения.
Теорема
умножения формулируется следующим
образом:
ХН антенны, состоящая
из n
направленных одинаковых и одинаково
ориентированных преобразователей,
равно произведению ХН одного преобразователя
и ХН гипотетической (воображаемой)
антенны, составленной из n
ненаправленных элементов, расположенных
в фазовых центрах реальных преобразователей.
Основные
ограничения при применении теоремы
умножения:
– нельзя применять,
если преобразователи ориентированы по
разному;
– в случае
цилиндрической или сферической антенны
элементы должны быть звукопрозрачными
или размещаться в акустически жестком
экране.
Если антенна
работает вне резонанса, то теорему
умножения практически всегда можно
использовать. Для антенны, работающей
в области, близкой к резонансу, теорему
умножения можно использовать, если
взаимные сопротивления излучателей
одинаковы, но в этом случае нужно
учитывать некоторую весовую функцию,
связанную с взаимодействием элементов
в антенне.
Следствие:
ХН непрерывной
антенны, полученной путем параллельного
переноса некоторой кривой вдоль
направляющей равна произведению ХН
этой кривой и направляющей.
Например:
– для прямоугольника
со сторонами
.
Теорема умножения
может быть применена последовательно
несколько раз понимая под элементом не
один, а группу преобразователей.
Для прозрачного
цилиндра ХН равна произведению ХН
направляющей и образующей.
Теорема смещения:
ХН (антенны в бесконечном экране) в
некоторой области L
(
экрану) совпадает с ХН отрезка, являющегося
проекцией рассматриваемой антенны на
плоскость.
Пример 1. Рассмотрим
применение теоремы смещения, когда
антенна представляет собой квадрат в
плоскости XOY
(рис. 17) с равномерным амплитудным
распределением.
y
O x
x
Рис. 17 Антенна в
форме квадрата
В рассмотренном
случае ХН антенны в форме квадрата
совпадает с ХН отрезка, лежащего на
диагонали квадрата и имеющего амплитудное
распределение в виде треугольника.
Пример 2. ХН в
виде круга радиусом R
с равномерным распределением амплитуды
в плоскости XOY
совпадает с ХН отрезка, являющегося
диаметром круга и имеющего амплитудное
распределение,
Пример 3. Антенна
представляет собой кольцо с равномерным
распределением в плоскости XOY.
В этом случае ХН кольца совпадает с ХН
отрезка длиной
,
имеющего амплитудное нормированное
распределение вида
Пример 4. Пусть
имеется дискретная антенна, состоящая
из 4-х элементов. ХН антенны, состоящей
из 4-х элементов, лежащих в вершинах
квадрата в плоскости XOY
и с диагональю вдоль оси x,
совпадает с ХН антенны, состоящей из
трех элементов и обладающих амплитудным
распределением 1,2,1.
Теорема сложения:
Давление,
развиваемое антенной, имеющее распределение
колебательной скорости по всей поверхности
f(S),
равно сумме давлений, развиваемых такими
же антеннами с амплитудными распределениями
A(S)
и B(S),
при условии выполнения следующего
равенства:
f(S)=A(S)+B(S)
В простейшем
случае при равномерном распределении
Пример 1.
предположим, что антенна представлена
в виде 2-х отрезков (удобно, чтобы фазовые
центры совпадали) длиной
,
а расстояние между отрезками d.
В соответствии
с теоремой сложения ХН такой антенны
равна разности отрезков ХН длиной
и
d
с коэффициентом пропорциональности их
длинам: Имея характеристику направленности
отрезка d
как
,
в соответствии с
теоремой сложения искомое выражение
для ХН антенны с длиной d
и
запишется
в виде
Пример 2. Пусть
дана антенна в виде кольца. ХН кольца,
где R
– внешний радиус,
r
– внутренний радиус, равна разности ХН
2-х плоских антенн в виде круга с радиусом
R
и r
с учетом коэффициента пропорциональности
их площадей:
Пример 3.
Рассмотрим применение теоремы сложения
для антенны, состоящей из 3-х ненаправленных
элементов, расположенных на одной прямой
с расстоянием между элементами d.
ХН такой антенны будет равна сумме ХН
антенн, составленных из 2-х крайних
элементов, с расстоянием между ними 2d,
и среднего элемента, с коэффициентом
пропорциональности 2 и 1:
Теоремой сложения
целесообразно пользоваться в тех
случаях, когда антенна при распределении
A(S)
и B(S)
имеет фазовый центр, расположенный в
одной и той же точке, что позволяет
производить сложения алгебраических,
и не комплексных величин.
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]
- #
- #
- #
- #
- #
- #
- #
- #
- #
- #
- #
Электродинамика, микроволновая техника, антенны
УДК 621.396.67
Д. О. Хабиров, М. А. Удров ОАО “Научный центр прикладной электродинамики” (Санкт-Петербург)
Методика определения координат центра излучения антенны и практические аспекты ее применения
Представлено теоретическое обоснование и описание практического применения метода расчета координат фазового центра (центра излучения) антенн на примере синфазной антенной решетки. Метод основан на определении координат центров кривизны фазового фронта волны в зависимости от угла поворота антенны и в приведении координат всех найденных центров кривизны к единой системе координат, привязанной к оси вращения антенны.
Антенна, фазовый центр, центр излучения, эволюта
При применении излучателей электромагнитных волн в различных системах (антенные решетки (АР), зеркальные антенны) важно иметь информацию о реальных координатах их фазового центра (ФЦ) или центра излучения. Ввиду значительного разнообразия конструкций антенн необходимо разработать удобную для практического использования методику определения координат ФЦ, результат применения которой не зависит от типа антенны.
В большинстве случаев антенны не имеют ФЦ в строгом смысле [1], поскольку поверхности равных фаз не являются сферами. Однако на практике обычно важно проанализировать фазовую диаграмму направленности (ДН) в каком-либо ограниченном секторе (например, в пределах главного лепестка). В большинстве случаев в таком ограниченном телесном угле поверхности равных фаз можно аппроксимировать сегментами поверхностей сфер [2].
В этом случае следует говорить о частичном фазовом центре, т. е. центре кривизны поверхности равных фаз в направлении, заданном углами ф (азимут) и 0 (угол места). Центр кривизны поверхности – точка математически определенная; она представляет собой центр сферы, совпадающей с поверхностью равных фаз в точке, определенной направлением, заданным указанными углами.
С практической точки зрения задача определения координат центра излучения включает следующие этапы:
– поиск математических соотношений, позволяющих определять центр кривизны плоской ли-30
нии равных фаз, полученной сечением поверхности равных фаз заданной плоскостью;
– объединение решений для координат центра излучения, полученных в плоских сечениях фазовой ДН, в единой трехмерной системе координат (данный этап в настоящей статье не рассмотрен). Линия равных фаз описывается уравнением [3]
р (ф) = г+у (ф)А ,
где г – расстояние от точки пересечения оси симметрии антенны и оси ее вращения до ФЦ измерительной антенны (при этом предполагается, что ФЦ измерительной антенны расположен на продольной оси симметрии рассматриваемой антенны, т. е. обе антенны находятся на одной высоте относительно горизонтальной поверхности (пола)); у – фаза электромагнитного поля; к = 2л/Х – волновое число, причем X – длина волны излучения.
Рис. 1 иллюстрирует связь данных параметров. Начало координат 0 располагается в точке пересечения плоскостей симметрии антенны (если для какой-либо антенны плоскости симметрии отсутствуют, начало координат можно расположить в произвольной точке антенны). Оси координат п и £ означают, соответственно, продольное и поперечное смещения в плоском сечении, в котором известна фазовая ДН у(ф). Сплошной кривой показана линия равных фаз р(ф), штриховой – окружность с центром в точке (^о, “По) и радиусом-вектором р0 , совпадающая с кривой
© Хабиров Д. О., Удров М. А., 2015
р(ф) по критерию минимального среднеквадратичного отклонения в направлении, характеризуемом углом фо (на рис. 1 не показан). Таким образом, Ро = р (фо ), а точка ((, “По ) является центром кривизны линии равных фаз в направлении фо.
Координаты центра кривизны линии равных фаз в направлении 0 находятся по известным из [4] формулам для радиуса кривизны и центра кривизны кривой, заданной в полярной системе координат.
Учитывая, что r^(1/к) у(ф), и пренебрегая малыми величинами, получим:
^о (ф) = (Vк) [cos ф V (ф) -sin ф v” (ф) ;
По (ф) = (Vк) [cos ф V (ф) + sin ф V (ф)],
где у'(ф) = d ф; у”(ф) = d2y^)/ d ф2.
Приведенные формулы позволяют найти координаты частичного ФЦ (или, в общем случае, эволюты, т. е. геометрического места центров кривизны) одномерной фазовой ДН через производные функции, описывающие эту диаграмму. Эволюта представляет траекторию перемещения точки ФЦ исследуемой антенны в каком-либо плоском сечении, перпендикулярном оси вращения, при повороте антенны вокруг данной оси вращения (рис. 2). В общем случае расположение ФЦ не совпадает с местом пересечения продольной оси симметрии антенны (штрихпунктирная линия на рис. 2) и оси ее вращения на опорно-поворотном устройстве.
В рамках лабораторных исследований с помощью векторного анализатора цепей проведены измерения фазовых ДН синфазной АР, состоящей из четырех цилиндрических спиральных облучателей и экранной плоскости. Облучатели располагались в вершинах квадрата со стороной о.Тк.
На рис. 2 упрощенно показана схема проведения измерений значений фаз. Линии 1-3 представляют собой сечения экранной плоскости исследуемой АР плоскостью, перпендикулярной оси вращения, для различных углов поворота ф (ось вращения проходит через точку 0 перпендикулярно плоскости рисунка). Линия 1 соответствует исходному положению АР (ф = 0), линии 2 и 3 – повороту АР на углы ф1 =-25° и ф2 = 25° соответственно.
Положение ФЦ определялось в системе координат, связанной с экранной плоскостью АР: начало координат находилось в точке пересечения продольной оси симметрии АР и оси вращения, т. е. в центре экранной плоскости. По осям координат п и £ (на рис. 2 не показаны) происходит продольное и поперечное смещения в сечении рисунка соответственно.
В рассмотренной схеме измерения предполагается, что ФЦ АР в исходном положении (ФЦ0) расположен с поперечным смещением £0 и нулевым продольным смещением (значение ^о = 0 принято для большей наглядности объяснения). При повороте АР на углы фц и ф2 точка ФЦ перемещается, оставаясь в ее экранной плоскости (значение ^0 = 0 сохраняется) и занимая позиции ФЦ1 и ФЦ2 соответственно. Эти позиции определяются значениями поперечных смещений £1 и £,2. В процессе поворота ФЦ перемещается по некоторой кривой, обозначенной на рис. 2 как “Траектория ФЦ”. В общем случае указанная траектория не является дугой окружности с центром в точке 0, так как для различных углов поворота АР вокруг оси вращения положение точки частичного ФЦ может сместиться. Таким образом, в общем
Ц Измерительная ФЦ /Д антенна
Электродинамика, микроволновая техника, антенны
45°
90°
у …
-20 -15 -10 -5 0 ___1
0°
у = 135°
10 15 20 ф, П
-10 Зо”
-30
-40
-50
-60 Рис. 3
случае £,0 ^ £1 ^ £2, а координаты точек ФЦ необходимо определять для каждого угла поворота АР.
На рис. 3 представлены результаты измерений набега фазы у(ф) электромагнитной волны от ФЦ до измерительной антенны при углах поворота исследуемой антенны относительно продольной оси симметрии, равных 0, 45, 90 и 135°. Измерения выполнены на частоте 1246 МГц.
На рис. 4 представлены построенные описанным способом эволюты (траектории частичного ФЦ) для этих же углов попорота АР вокруг продольной оси симметрии. Система координат рис. 4 аналогична системе координат рис. 2. Круглыми маркерами на эволютах обозначены положения частичного ФЦ АР при ф = 0; маркеры в начале координат указывают проекцию оси вращения АР на опорно-поворотном устройстве.
Учитывая зависимость продольной п и поперечной £ составляющих смещения частичных ФЦ от угла поворота ф антенны вокруг оси вращения, выражения для определения координат центра излучения антенны можно записать следующим образом:
2л i-
£ц = W^2 + 0
х cos^-arctg [“(ф)/£(ф)]}|ф(ф)| dф; (1) 2 л i-
Лц = W^2 (ф) + “П2 (ф) х
0
х sin |ф – arctg [л(ф)Д(ф)]}|Ф(ф)| dФ, (2)
где |Ф(ф)| – амплитудная ДН, используемая в качестве весовой функции. Введение амплитудной ДН в формулы (1), (2) обеспечивает независимость точности расчета координат центра излучения от пределов интегрирования [3].
После преобразования [5] получим:
£ц = 1?мф)]2 +1У(ф)]2 х k 0
П, мм
у = 0°
у = 45°
/
П, мм
-5
-10
_L
-20 -15
мм
П, мм
10
мм
-5
-10
х cos
Ф – arctg
у(ф) + tg фу'(ф) У(ф)- tg фУ(ф).
|Ф(Ф)| dф;
у = 90°
П, мм
10
мм
Рис. 4
При определении положения центра излучения (координаты £ц и ^ц) необходимо учитывать, что
координаты частичного ФЦ (ЧФЦ) для угла поворота АР ф, отличного от нуля, следует относить не к исходной системе координат, а к системе, повернутой на этот угол вместе с экранной плоскостью АР (рис. 5). На рис. 5 указанные координаты обозначены как £ ц и ^ц, остальные обозначения аналогичны принятым на рис. 2. 32
-5
-10
-15
у = 135°
20
/
мм
J “i l
лц тлц чфц ^———— ▼ фт
Рис. 5
5
5
0
0
0
0
5
5
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
5
5
^ц
2 +[у(ф)]2 *
о
х sin
Ф – arctg
V Чф) + tg Ф У(Ф) у(ф)- tg фу(ф)
|ф(ф)| dф.
Приведенные выражения позволяют определить координаты центра излучения антенны, для которой известна (получена по результатам измерений либо теоретических расчетов) одномерная амплитудно-фазовая ДН. При этом координаты определяются в плоскости сечения, перпендику-
СПИСОК ЛИТЕРАТУРЫ
лярной оси вращения антенны на опорно-поворотном устройстве (рис. 2, 4, 5).
Совокупность плоских сечений ДН, проходящих через одну ось, формирует объемную ДН [3], поэтому реальное положение центра излучения может быть определено сведением решений для координат центра излучения, полученных в плоских сечениях фазовой ДН, в единую трехмерную систему координат. Указанная операция может быть проведена как аналитическим, так и графическим методами.
1. Вольперт А. Р. О фазовом центре антенны // Радиотехника. 1961. Т. 16, № 3. С. 3-12.
2. Родс Д. Р. Введение в моноимпульсную радиолокацию. М.: Сов. радио, 1960. 160 с.
3. Вендик О. Г., Парнес М. Д., Бахрах Л. Д. Антенны с электрическим сканированием (введение в теорию). М.: Science Press, 2001. 252 с.
4. Выгодский М. Я. Справочник по высшей математике. М.: Наука, 1977. 872 с.
5. Пискунов Н. С. Дифференциальное и интегральное исчисления для втузов. 13-е изд. М.: Наука. Гл. редакция физ.-мат. лит., 1985. Т. 1. 432 с.
D. O. Khabirov, M. A. Udrov Open JSC “Scientific center of applied electrodynamics” (Saint Petersburg)
Antenna radiation center coordinates finding method and its applying practical aspects
The theoretical justification and original method of antenna phase center (radiation center) coordinates calculation practical application are presented by the example of four-helix array. The method is based on wave phase front curvature center coordinates finding as function of antenna rotation angle, and on found curvature centers coordinates transformation to single coordinate system, which is related to antenna rotation axis.
Antenna, phase center, radiation center, evolute
Статья поступила в редакцию 11 июня 2015 г.
Изобретение относится к антенным измерениям с использованием сверхширокополосных (СШП) сигналов и может быть использовано при разработке, испытаниях и калибровке антенн. Измерительную и испытуемую антенны размещают в дальней зоне, а для зондирования используют СШП сигналы. При предварительном зондировании поворачивают испытуемую антенну вокруг оси вращения на выбранные углы и находят фиксированное окно приема такое, чтобы принятые сигналы попадали в него. При основном зондировании оценивают в найденном окне разность времен распространения сигналов между фазовыми центрами антенн при разных углах поворота испытуемой антенны. Для этого либо оценивают различие в положении принятых СШП сигналов во временном окне приема, либо вычисляют их фазочастотный спектр, а координаты фазового центра испытуемой антенны вычисляют для спектра частот. Находят ось фазового центра испытуемой антенны, относительно которой время распространения сигналов между фазовыми центрами антенн не зависит от угла поворота испытуемой антенны. Для определения пространственного положения фазового центра антенны изменяют ось вращения испытуемой антенны и находят другую ось фазового центра. Фазовый центр испытуемой антенны находят на пересечении осей фазового центра. Техническим результатом является обеспечение быстрого и точного определения положения фазового центра испытуемой антенны для спектра частот. 2 з.п. ф-лы, 3 ил.
Изобретение относится к измерению электрических и магнитных величин, в частности к антенным измерениям с использованием сверхширокополосных (СШП) сигналов, и может быть использовано при разработке, испытаниях и калибровке антенн.
Фазовый центр антенны – это точка, в которую можно поместить одиночный излучатель сферической волны, эквивалентный рассматриваемой антенной системе в отношении фазы создаваемого поля [1]. В реальных антеннах фазовый центр обычно рассматривается в рамках ограниченных углов главного лепестка диаграммы направленности. Положение фазового центра зависит от частоты используемого сигнала, направления излучения/приема антенны, его поляризации и других факторов. Некоторые антенны не имеют фазового центра в общепринятом понимании.
В простейших случаях, например у параболической антенны, фазовый центр совпадает с фокусом параболоида и может быть определен из геометрических соображений. Однако конструктивные и технологические погрешности приводят к смещению фазового центра даже в простейших конструкциях антенн. Для этих случаев известен способ определения фазового центра антенны [2], заключающиеся в том, что измерительную антенну помещают в область двойного фокусного расстояния исследуемой параболической антенны, перемещают антенную систему в точку максимального приема отраженного сигнала, определяют координаты указанной точки и, принимая их за центр аппроксимирующей параболоид сферы, находят фокус параболоида, который рассматривается как фазовый центр исследуемой параболической антенны. Определяя направление вектора Пойнтинга поля излучения антенны, можно определить оси фазового центра, которые пересекаются в фазовом центре антенны. Для определения точки пересечения достаточно определения двух осей.
Недостатком этого способа является ограниченная область применения – только для параболических антенн, а также значительная трудоемкость проведения измерений при необходимости определения фазового центра для спектра частот.
В более сложных случаях, например, рупорных антенн, положение фазового центра не очевидно и требует соответствующих измерений. Известен способ определения фазового центра излучающего рупора [3], заключающийся в возбуждении рупора СВЧ сигналом, приеме сигнала отраженного от специального экрана, оценке фазы принятого сигнала и определении координат фазового центра испытуемой антенны.
Недостатком этого способа является ограниченная область применения – только для рупорных антенн, а также значительная трудоемкость проведения измерений при необходимости определения фазового центра для спектра частот.
Известен способ определения фазового центра элемента антенной решетки [4], заключающийся в том, что устанавливают две антенны в дальней зоне, зондируют элемент испытуемой антенной решетки, для чего поворачивают его вокруг оси вращения, в каждом положении излучают сигналы образцовой антенной, принимают их испытуемой антенной, оценивают амплитуду и фазу принятых сигналов и находят вектор центра фазы, минимизирующий различие между измеренным и вычисленным рисунками фазы.
Недостатком способа [4] является высокая трудоемкость проведения измерений при необходимости определения фазового центра антенны для спектра частот.
Наиболее близким к заявляемому является способ определения положения фазового центра антенны [5], заключающийся в том, что устанавливают две антенны в дальней зоне, зондируют испытуемую антенну, для чего поворачивают ее вокруг оси вращения на выбранные углы, в каждом положении излучают сигналы с постоянными характеристиками одной, принимают их измерительной антенной и оценивают принятые сигналы, фазовый центр испытуемой антенны находят на пересечении осей фазового центра. При этом измерительную антенну размещают последовательно в двух точках электромагнитного поля, при каждом положении измеряют амплитуды и фазы принятого сигнала для трех ортогональных компонент вектора электрического поля, определяют две оси фазового центра, начинающиеся в фазовых центрах вспомогательной антенны (в двух положениях) и заканчивающиеся в фазовом центре испытуемой антенны.
Недостатками способа являются пригодность только при эллиптической поляризации антенны, низкая точность, а также высокая трудоемкость проведения измерений при необходимости определения фазового центра для спектра частот. Низкая точность объясняется тем, что точное определение компонент вектора электрического поля является непростой задачей, а нахождение точки пересечения двух не точно известных векторов из дальней зоны приводит к еще более существенным погрешностям.
Знание точного положения фазового центра антенны имеет большое значение для высокоточных навигационных измерений, поскольку ошибки в определении положения фазового центра непосредственно влияют на точность определения навигационных параметров объектов. Мероприятия по уточнению положения фазового центра принимаются при проектировании, изготовлении и калибровке антенн. В современных антенных системах широко применяются антенные решетки, для которых само понятие «фазовый центр» становится сложным, а тем более его измерение [6]. Например, в Глобальных Системах Позиционирования (GPS) точные эфемериды, поставляемые Международной Геодинамической Службой (IGS) и другими организациями, передают информацию о центре массы спутника, а при генерации и использовании таких эфемерид для обработки GPS данных необходимо знать точное положение фазового центра антенны по отношению центру масс спутника с учетом конструктивных погрешностей изготовления антенн и спутника, влияния самого навигационного спутника на работу антенны, углов наблюдения, поляризации и других факторов. Для повышения точности прибегают к калибровке [7]. Кроме того, современные антенные системы часто используют сложные, немонохроматические сигналы, что существенно затрудняет антенные измерения.
Задачами, решаемыми заявляемым способом, является быстрое и точное определение положения фазового центра испытуемой антенны для спектра частот.
Для решения указанной задачи способ определения положения фазового центра антенны, заключается в том, что устанавливают две антенны в дальней зоне, зондируют испытуемую антенну, для чего поворачивают ее вокруг оси вращения на выбранные углы, в каждом положении излучают сигналы с постоянными характеристиками одной, принимают их другой антенной и оценивают принятые сигналы, фазовый центр испытуемой антенны находят на пересечении осей фазового центра, для зондирования используют СШП сигналы, проводят предварительное зондирование, при котором оценивают и выбирают минимальный размер временного окна приема и его положение относительно момента излучения такие, чтобы принятые сигналы попадали в окно приема, проводят основное зондирование, при котором принимают сигналы в выбранном временном окне приема, оценивают разности времен распространения сигналов между фазовыми центрами антенн при разных углах поворота испытуемой антенны и находят параллельную оси вращения ось фазового центра испытуемой антенны, относительно которой время распространения сигналов между фазовыми центрами антенн не зависит от угла поворота испытуемой антенны, выбирают другую ось вращения испытуемой антенны, повторяют предварительное и основное зондирования и находят другую ось фазового центра.
При основном зондировании для грубого определения положения фазового центра антенны в качестве оценки разности времен распространения сигналов между фазовыми центрами антенн используют различие в их положении во временном окне приема,
Для точного определения положения фазового центра антенны для оценки разности времен распространения сигналов между фазовыми центрами антенн при основном зондировании вычисляют их фазочастотный спектр, а координаты фазового центра испытуемой антенны вычисляют для спектра частот.
Существенными отличиями заявляемого способа по сравнению с прототипом являются:
В качестве зондирующих используют СШП сигналы. Такие сигналы позволяют ускорить антенные измерения за счет их проведения одновременно в широком спектре частот.
В прототипе применяются монохроматические сигналы. Работа с такими сигналами технически проще и очевиднее, так как характеристики антенны снимаются непосредственно на каждой частоте по отдельности. Однако при необходимости определения положения фазового центра для спектра частот трудоемкость измерений многократно возрастает.
Проводят предварительное зондирование, при котором оценивают и выбирают минимальный размер временного окна приема и его положение относительно момента излучения такие, чтобы принятые сигналы попадали в окно приема. Выбор фиксированного временного окна приема позволяет учесть все информативные части принятого сигнала при всех выбранных углах поворота испытуемой антенны, дает возможность проводить дальнейшие измерения в одинаковых условиях и обеспечивает точность вычисления фазочастотной характеристики принятого сигнала при последующей обработке результатов измерений. Минимальный размер окна приема также обеспечивает повышение точности измерений, поскольку при выбранном количестве отсчетов в окне приема шаг дискретизации по времени оказывается минимальным. Кроме того, в результате выбора минимального временного окна приема из-за разницы хода в него не попадают сигналы, отраженные от посторонних объектов в зоне исследований, а за счет этого снижается стоимость антенных измерений.
В прототипе прием сигнала происходит непрерывно, не синхронизовано с моментом излучения. При этом фаза монохроматического сигнала определяется неоднозначно, что приводит к необходимости использования специальных схем, например, [3] АС СССР №1125559. Кроме того, решение проблемы с побочными отражениями в этом случае сопровождается серьезными материальными затратами на создание безэховых камер, поглощающих покрытий, работающих лишь в определенных диапазонах частот, проведение дорогостоящих облетных измерений и т.д.
Проводят основное зондирование, при котором принимают сигналы в выбранном временном окне приема, оценивают разности времен распространения сигналов между фазовыми центрами антенн при разных углах поворота испытуемой антенны и находят параллельную оси вращения ось фазового центра испытуемой антенны, относительно которой время распространения сигналов между фазовыми центрами антенн не зависит от угла поворота испытуемой антенны. Главным отличием является использование разностей во временах распространения сигналов, что позволяет перейти к измерениям положения фазового центра относительно оси вращения испытуемой антенны.
В прототипе находят оси фазового центра, исходящие из мест расположения измерительной антенны в дальней зоне, что порождает существенные погрешности измерений.
Выбирают другую ось вращения испытуемой антенны, повторяют предварительное и основное зондирования и находят другую ось фазового центра. Другая ось вращения может быть выбрана произвольно, в том числе и под углом 90° к первой оси, что повышает точность нахождения точки пересечения осей, а, следовательно, координат фазового центра.
В прототипе оси фазового центра расположены в рамках ограниченного угла основного лепестка диаграммы направленности, что снижает точность нахождения точки пересечения указанных осей.
В качестве оценки разности времен распространения сигналов при основном зондировании используют различие в их положении во временном окне приема. Такой способ позволяет лишь грубо оценить положение фазового центра без учета зависимости положения фазового центра от частоты.
В прототипе отсутствие синхронизации между моментами излучения и приема не позволяет найти отличие во временах приема сигнала при поворотах испытуемой антенны.
Более точное положение фазового центра может быть получено тогда, когда для оценки разности времен распространения сигналов при основном зондировании вычисляют их фазочастотный спектр. Разности фаз принятого сигнала на каждой частоте могут быть легко пересчитаны в разности времен распространения сигнала, а координаты фазового центра испытуемой антенны вычисляют для спектра частот. Такой подход позволяет снизить трудоемкость работ за счет того, что за один цикл измерений удается получить спектр координат фазовых центров испытуемой антенны.
В прототипе определяется фазовый центр антенны для одной частоты. При необходимости получить указанные результаты для спектра частот возникают большие временные и аппаратурные затраты на перестройку генераторов зондирующих сигналов, фазометров и других элементов устройств, реализующих такой способ.
Заявляемый способ иллюстрируют следующие графические материалы:
Фиг.1 – Схема для грубого расчета положения фазового центра антенны.
Фиг.2 – Схема для точного расчета положения фазового центра антенны.
Фиг.3 – Устройство, реализующее заявляемый способ.
Рассмотрим возможности реализации заявляемого способа.
Перед началом измерений, фиг.1, устанавливают две антенны в дальней зоне, т.е. выбирают расстояние L между предполагаемыми фазовыми центрами испытуемой 1 и измерительной 2 антенн, причем само расстояние L не имеет значения, поскольку измерения проводятся относительно оси вращения испытуемой антенны. Характеристики измерительной антенны также не имеют значения, поскольку ее амплитудно-частотная характеристика, положение ее фазового центра также не влияют на измерения. Выбирают углы поворота испытуемой антенны так, чтобы в секторе измерений не было нулей амплитудной диаграммы направленности испытуемой антенны, которым соответствуют скачки фазы. В качестве излучающей может выступать испытуемая или измерительная антенна, соответственно другая антенна оказывается принимающей.
Для измерений излучают СШП сигналы одной и принимают другой антенной. В качестве таких сигналов могут использоваться одинаковые (с постоянными характеристиками) короткие видеоимпульсы длительностью от долей до единиц пикосекунд, которые имеют спектр от нуля до десятков гигагерц. Этот импульс может быть принят стробоскопическим преобразователем [8] – устройством, которое в фиксированном по отношению к моменту излучения зондирующего сигнала окне приема из принятого сигнала выделяет один отсчет в заданный момент времени. Зондирование такими же импульсами при одном и том же положении антенн повторяется многократно (тысячи раз) с произвольной частотой, а моменты времени выделения отсчета смещают по окну приема. В результате такого способа приема происходит масштабно-временное преобразование, т.е. принятый зондирующий «быстрый» импульс воспринимается как совокупность его отсчетов, но в другом масштабе времени. При этом появляется возможность обрабатывать каждый из отсчетов «медленно», с использованием обычных аналого-цифровых преобразователей и компьютерных способов обработки.
Пусть перед началом измерений, фиг.1, ось вращения 0Y испытуемой антенны 1 перпендикулярна плоскости XOZ и не совпадает с параллельной ей осью фазового центра, проходящей через точку Zц (ФЦ). При повороте испытуемой антенны 1 на угол αi изменяется расстояние между фазовыми центрами антенн 1 и 2, а следовательно, время распространения СШП сигнала между ними на величину ΔТ и фаза принятого сигнала на ΔФ(f). Величина изменения фазы зависит от рассматриваемой частоты f.
Проводят предварительное зондирование, при котором поворачивают испытуемую антенну на выбранные углы. Излучают и принимают СШП сигналы. Подбирают время начала окна приема сигнала антенной 2 так, чтобы при любых углах поворота α испытуемой антенны 1 наблюдалось начало импульса, принятого антенной 2. Затем изменяют и находят минимальную длительность временного окна приема такую, чтобы при любых углах поворота испытуемой антенны 1 весь принятый импульс попадал в него.
Проводят основное зондирование, при котором принимают СШП сигналы в выбранном временном окне приема. Задача основного зондирования состоит в нахождении оси фазового центра, проходящей через фазовый центр (ФЦ) и параллельной оси вращения 0 антенны 1.
Для грубой оценки положения оси фазового центра антенны 1 оценивают разности времен распространения сигналов между фазовыми центрами антенн 1 и 2 при разных углах поворота испытуемой антенны 1. В частности, можно найти ось 0Z (базовую ось), проходящую через ось вращения антенны 1, фазовый центр антенны 2 и пересекающую ось фазового центра антенны 2 в точке Zц. Для этого находят такое угловое положение антенны 1, при котором принятый сигнал наиболее близок к началу окна приема (если фазовый центр смещен вперед по отношению к оси вращения) или к концу окна приема (если фазовый центр находится сзади оси вращения). Затем поворачивают антенну 2 на известный угол αi и определяют разность ДТ времен распространения СШП сигнала в указанных двух положениях. Из геометрических соображений можно определить неизвестную величину:
Zц=СΔТ/(1-cos(αi)),
где С – скорость света. Величину Zц и положение оси 0Z однозначно определяют положение оси фазового центра.
Низкая точность описанного способа определения оси фазового центра объясняется следующими обстоятельствами:
1. Положение фазового центра антенны зависит от частоты, а поэтому координата Zц найденная для СШП сигнала, является лишь неким «интегральным приближением» к фазовому центру, но может использоваться, например, при измерениях дальности.
2. Оценка временного положения СШП сигнала в окне приема неоднозначна, поскольку при поворотах антенны 1 изменяется форма принятых сигналов, так что какой момент принятого сигнала следует считать моментом приема, не очевидно.
3. Положение оси 0Z определяется не точно, поскольку в окрестности αi=0 разность ΔT времен распространения СШП сигнала меняется мало.
Для решения последней проблемы можно провести измерения по крайней мере для трех разных угловых положений антенны 1 и решить, как будет показано ниже, соответствующую систему уравнений.
Точное определение положения оси фазового центра возможно только в отдельности для каждой частоты fj требуемого спектра частот 
. Для решения этой задачи, используя дискретное преобразование Фурье (ДПФ), вычисляют фазочастотный спектр принятого сигнала при нескольких 
угловых положениях антенны 1. При n=3 методика расчета положения оси фазового центра антенны 1 для одной частоты fj состоит в следующем. Пусть ось фазового центра антенны 1 на частоте fj находится на расстоянии Zц от оси вращения Y ортогональной системы координат XYZ. По определению фазового центра расстояние Zn остается неизменным при любых выбранных углах ее поворота. Начальное угловое положение антенны 1, фиг.2, отстоящее на неизвестный от оси 0Z угол α0, будем считать базовым, при этом ось фазового центра антенны 1 проходит через точку Zц0 с неизвестными координатами. Проводят зондирование, по результатам которого с использованием ДПФ вычисляют фазу Ф0(fj) принятого сигнала. При повороте антенны 1 вокруг оси 0Y относительно базового положения на известный угол αi ось фазового центра будет проходить через точку Zц1. В результате аналогичных зондирований и вычислений находят фазу Ф1 принятого сигнала на той же частоте, но при другом положении антенны 1. Разность этих фаз позволяет оценить разность времен распространения сигналов:
ΔТ1=(Ф0-Ф1)/2πfj.
При повороте антенны 1 на известный угол α2 относительно базового положения ось фазового центра антенны проходит через точку Zц2. Аналогичным образом вычисляют разность времен распространения сигналов между базовым и текущим положениями оси фазового центра ΔT2. В результате трех зондирований может быть составлена система двух уравнений:
ΔT1=Zц(cosα0-cos(α0-α1))/C
ΔT2=Zц(cosα0-cos(α0-α2))/C,
где С – скорость света.
Эта система содержит два неизвестных α0 и Zц. и может быть решена известными методами. Полученные величины α0 и Zц являются полярными координатами оси фазового центра антенны 1 для частоты fj относительно оси ее вращения 0.
Аналогичные вычисления проводят для всех частот fj спектра частот . Для повышения точности вычислений количество углов поворота испытуемой антенны 1 выбирают большим трех, тогда система уравнений становится избыточной, а ее решение может быть получено, например, методом наименьших квадратов (МНК).
Фазовый центр антенны 1 может быть расположен не только на плоскости XOZ, но и отстоять от нее на величину Yц. Для нахождения пространственного положения фазового центра антенны 1 изменяют ось ее вращения. В качестве новой оси может быть выбрана ось 0Х. Повторяя описанные выше измерения и вычисления, находят вторую ось фазового центра. В рассматриваемом случае вторая ось вращения антенны 1 перпендикулярна первой. В идеальном случае найденные оси фазового центра пересекаются. Вычисляют точку пересечения указанных осей, которая считается фазовым центром антенны 1. В реальных условиях найденные оси оказываются перекрещивающимися. В этом случае находят фазовый центр антенны 1, минимизируя расстояние между осями, например, по МНК.
Таким образом, заявляемый способ позволяет быстро и точно определить положение фазового центра испытуемой антенны для всех выбранных частот. Полученные координаты фазовых центров позволяет использовать сложные сигналы при радиолокационных измерениях, учитывать вариации положения фазового центра, а за счет этого повысить точность измерений с использованием калиброванной антенны.
Устройство, реализующее заявляемый способ, изображено на фиг.3, где:
1 – испытуемая антенна;
2 – измерительная антенна;
3 – опорно-поворотное устройство;
4 – компьютер;
5 – линия задержки;
6 – генератор зондирующих сигналов;
7 – стробоскопический приемник;
8 – аналого-цифровой преобразователь;
9 – вход ручного управления задержкой.
Испытуемая антенна 1 является объектом измерений с неизвестным положением фазового центра. Измерительная антенна 2 предназначена для проведения измерений, ее характеристики могут быть неизвестны, поскольку не влияют на точность из-за относительного характера измерений.
Опорно-поворотное устройство 3 предназначено для вращения антенны 1 на известные углы под управлением кодовых посылок от компьютера 4.
Компьютер 4 управляет работой устройства, обрабатывает результаты измерений и вычисляет координаты фазового центра испытуемой антенны 1.
Линия задержки 5 предназначена для выбора и фиксации окна приема зондирующих сигналов, а также для изменения положения отсчетов принятого сигнала во временном окне приема. Линия задержки содержит блоки грубой задержки и точной задержки. Первый блок позволяет задержать момент излучения СШП сигнала относительно момента его приема, т.е. задает начало окна приема. Он может быть выполнен на генераторе тактовых сигналов и цифровом счетчике, число пересчета которого управляется от компьютера 4. Блок точной задержки состоит из цифроаналогового преобразователя и диода с накоплением заряда. Код, поступающий от компьютера 4, устанавливает порог срабатывания диода, который изменяет задержку запускающего сигнала.
Генератор зондирующих сигналов 6 формирует СШП сигналы.
Стробоскопический приемник 7 выделяет из принятого сигнала один отсчет по стробирующему сигналу.
Аналого-цифровой 8 преобразователь предназначен для перевода отсчетов принятого сигнала в цифровую форму.
Вход ручного управления задержкой 9 предназначен для выбора окна приема при предварительном зондировании.
Испытуемая 1 и измерительная антенны 2 устанавливаются в дальней зоне на ориентировочно известном расстоянии L. Испытуемая антенна 1 крепится на опорно-поворотном устройстве 3 с фиксированной осью вращения. Будем считать, что испытуемая антенна 1 является излучающей, а измерительная антенна 2 – приемной. Перед зондированиями компьютер 4 путем посылки управляющих кодовых сигналов в опорно-поворотное устройство 3 устанавливает требуемые углы αi поворота измеряемой антенны 1 в рамках выбранного диапазона.
Для стробоскопического приема СШП сигнала изменяют величину точной задержки в блоке 5. В результате чего момент излучения антенной 1 смещается относительно фиксированного момента приема антенной 2, а стробоскопический приемник 7 выделяет другой (по времени) отсчет принятого сигнала. Диапазон величин точной задержки определяет длительность окна приема, а шаг – точность измерения. Многократно изменяя величину точной задержки, осуществляют масштабно-временное преобразование и получают все отсчеты принятого СШП сигнала, которые отображают на экране монитора компьютера 4.
При зондированиях запускающий сигнал от компьютера 4 поступает через линию задержки 5 к генератору зондирующих СШП сигналов 8, а также по кабелю длиной, приближенно равной расстоянию между антеннами L – к стробоскопическому преобразователю 7 и аналого-цифровому преобразователю 8. СШП сигнал, излученный антенной 1, поступает к антенне 2 с задержкой на время распространения при расстоянии L между антеннами 1 и 2.
При предварительном зондировании сначала подбирают величину грубой задержки в блоке 5, а возможно, и длину упомянутого выше кабеля так, чтобы запускающий сигнал поступил на стробоскопический приемник 7 непосредственно перед поступлением в него СШП сигнала при любых углах поворота антенны 1. Для решения этой задачи компьютером 4 через опорно-поворотное устройство 3 изменяют углы поворота антенны 1, наблюдают на мониторе компьютера 4 изменение положения начала принятого сигнала в окне приема. Вручную, через вход 9 компьютера 4 изменяют величину грубой задержки в блоке 5. Затем через вход 9 компьютера 4 изменяют величину и шаг точной задержки в блоке 5, добиваясь полного приема СШП сигнала. В результате описанных действий фиксируют величину грубой задержки (момент начала приема), а также шаг и величину точной задержки (точность и длительность окна приема).
При грубом определении положения оси фазового центра антенны 1 оценивают положение принятых сигналов при известных углах поворота αi антенны 1 и по описанным выше алгоритмам вычисляют в компьютере 4 координаты точки, через которую проходит ось фазового центра.
При точном определении оси фазового центра в компьютере 4 вычисляют фазочастотный спектр принятого сигнала и находят координаты точек, через которые проходят оси фазовых центров на соответствующих частотах.
Изменяют ось вращения антенны 1 и повторяют измерения.
Фазовые центры антенны 1 определяют для каждой частоты отдельно как точки пересечения соответствующих осей фазовых центров.
Таким образом, заявляемый способ может быть реализован на современной элементной базе и позволяет быстро и точно определить координаты фазовых центров испытуемой антенны для спектра частот. Знание этих координат позволяет повысить точность антенных измерений с использованием измеряемой антенны.
ЛИТЕРАТУРА
1. Драбкин А.Л., Зузенко В.Л. Антенно-фидерные устройства. М.: Сов. радио, 1961, с.70-71.
2. АС СССР №364908.
3. АС СССР №1125559.
4. Патент JP №2000321314.
5. АС СССР №1702325.
6. Патент JP 2183172.
7. Калибровка фазовых центров антенны, GPS World, Май 2002, Издатель: Advanstar Communications Inc 859 Willamette Street, Eugene, Oregon 97401-6806, USA.
8. Рябинин Ю.А. Стробоскопическое осциллографирование. – М.: Сов. радио, 1972.
1. Способ определения положения фазового центра антенны, заключающийся в том, что устанавливают две антенны так, чтобы одна по отношению к другой находилась в дальней зоне, зондируют испытуемую антенну, для чего поворачивают ее вокруг оси вращения на выбранные углы, в каждом положении излучают сигналы с постоянными характеристиками одной, принимают их другой антенной и оценивают принятые сигналы, фазовый центр испытуемой антенны находят на пересечении осей, проходящих через ее фазовый центр, отличающийся тем, что для зондирования используют сверхширокополосные сигналы, проводят предварительное зондирование, при котором оценивают и выбирают минимальный размер временного окна приема и его положение относительно момента излучения такие, чтобы принятые сигналы попадали в окно приема, проводят основное зондирование, при котором принимают сигналы в выбранном временном окне приема, оценивают разности времен распространения сигналов между фазовыми центрами антенн при разных углах поворота испытуемой антенны и находят параллельную оси вращения ось, проходящую через фазовый центр испытуемой антенны, относительно которой время распространения сигналов между фазовыми центрами антенн не зависит от угла поворота испытуемой антенны, выбирают другую ось вращения испытуемой антенны, повторяют предварительное и основное зондирования и находят другую ось, проходящую через фазовый центр испытуемой антенны.
2. Способ по п.1, отличающийся тем, что в качестве оценки разности времен распространения сигналов при основном зондировании используют различие в их положении во временном окне приема.
3. Способ по п.1, отличающийся тем, что для оценки разности времен распространения сигналов при основном зондировании вычисляют их фазочастотный спектр, а координаты фазового центра испытуемой антенны вычисляют для спектра частот.
