Методическая и шумовая погрешности восстановления диаграммы направленности антенны при измерениях в неидеальных условиях, представимых сходящимся пучком плоских волн
Ю.И. Чони, И.Ю. Данилов
Казанский национальный исследовательский технический
университет им. А.Н. Туполева ? КАИ
АО ”Информационные спутниковый системы
им. академика М.Ф. Решетнёва”
Аннотация
Анализируется влияние методических и технических факторов на точность восстановления диаграммы направленности (ДН) антенны, измеренной в неидеальных условиях, под которыми понимаются любые отличия облучающего поля от плоской волны, в частности, из-за относительно близкого расположения вспомогательной антенны. Процедура восстановления ДН включает реконструкцию облучающего поля в виде сходящегося пучка плоских волн путем измерения и обработки сигнал от зонда, перемещаемого по граничной поверхности рабочей зоны. На примерах двухмерной скалярной задачи исследуется устойчивость задачи восстановления ДН, потенциальная точность и влияние погрешностей измерений, демонстрируются рельефы полей пучков волн внутри и вне рабочей зоны в вариантах цилиндрической или сферической облучающей волны. Представлены зависимости, позволяющие оценить достижимую точность восстановления ДН.
Ключевые слова: диаграмма направленности, сходящийся пучок плоских волн, алгоритм восстановления, погрешности, устойчивость, точность.
Abstract
The paper discusses methodological aspects, computing peculiarities and technical factors with regard to their impact on the accuracy of antenna pattern retrieval through the data, measured in non-ideal environments. We mean any circumstances giving rise to differences between the field in the working area and the ideal plane wave in the first place, a rather short distance from the illuminating antenna. The pattern retrieval procedure is based on the illuminating field reconstruction in the form of a converging cluster of plane waves. It is obtained by processing the data from the probe that is moved along the boundary surface of the working area. With regard to the two-dimensional scalar tasks we investigate the stability of a pattern retrieval problem, potential accuracy and the influence the of measurement errors. In cases cylindrical or spherical illuminating wave, the resultant fields of converging cluster of plane waves inside and outside the working area are demonstrated. The evaluation of the achievable retrieval accuracy is presented.
Key words: antenna radiation pattern, converging cluster of plane waves, antenna measurement accuracy, reconstruction algorithm, stability, precision.
Введение
Современные векторные анализаторы цепей позволяют в автоматическом режиме проводить амплифазометрические измерения с последующей цифровой обработкой данных. Благодаря этому, находят применение процедуры реконструкции диаграммы направленности (ДН) антенн при измерениях в далеко неидеальных условиях, когда облучающее поле существенно отличается от плоской волны из-за малого расстояния до вспомогательной антенны и/или отражений от стен. Тридцатилетней давности первые работы в этом направлении носили скорее теоретический, чем практический характер [1 ? 2].
С тех пор выполнено немало работ, касающихся различных аспектов первичных измерений и алгоритмов обработки данных, например, [3]. Недавно защищенные диссертации [4, 5] подтверждают актуальность исследований в этой области и в наши дни. В работе [4] выделены три группы методов реконструкции ДН: 1) плоско параллельное сканирование ближнего поля антенны с пересчетом в дальнюю зону, 2) измерение поля антенны в зоне Френеля с пересчетом в дальнюю зону, 3) обследование облучающего поля в рабочей области с помощью антенны или зонда с известной ДН для последующей коррекции ДН испытуемой антенны. Рассматриваемый нами метод сходящегося пучка плоских волн (СППВ) относится к третьей группе, в отечественной литературе именуемой методом эталонной антенны [2], а в англоязычной ? «test zone field compensation», что полнее отражает суть дела.
Идеология сходящегося пучка плоских волн
Возможность использования СППВ для аппроксимации поля в области, лишенной источников, базируется на двух обстоятельствах [3]. Во-первых, функция Грина свободного пространства и вместе с ней произвольное поле в любой области вне источников представимы интегральным преобразованием
которое с физической точки зрения есть сплошный телесный (двумерный) пучок плоских электромагнитных волн, разбегающихся от расположенного в начале координат источника. Этот пучок включает в себя как однородные (), так и поверхностные () волны. При антенных измерениях неоднородности, способные порождать поверхностные волны, значительно удалены от испытуемой антенны. Поэтому их вкладом в облучающее поле можно пренебречь, и представить его в рабочей области как СППВ (сходящийся, т.к. источники поля расположены за ее пределами). Во-вторых, поскольку реакцию антенны на однородную плоскую волну определяет ее ДН, то сигнал на выходе антенны, облучаемой СППВ, дается интегралом типа свертки, ядром которого является ДН антенны.
Ограничимся двухмерной скалярной ситуацией (рис. 1) типичной для подавляющего большинства публикаций на эту тему, когда ДН F(ц) измеряется в плоскости, и облучающее ее поле есть СППВ в той же плоскости. Пусть облучающее поле создается источником цилиндрической волны, расположенным относительно близко в точке D. На первом этапе определяются комплексные амплитуды A(ц) волн пучка, воспроизводящего это поле. Для этого зонд, ДН которого f0(ш) известна, вращается по окружности радиуса R0 (рис. 1(а)). Здесь «радужное» кольцо изображает СППВ в виде функции A(ц)пространственного спектра волн пучка, но не порождаемое им интерференционное поле. Белая пунктирная стрелка с поперечным отрезком, символизирующим фазовый фронт, изображает отдельную плоскую волну.
Очевидно, что принимаемый зондом сигнал связан с функцией A(ц) интегралом свертки
, (1)
где ? угол поворота зонда, ? ДН зонда, смещенного из начала координат. По измеренному сигналу , из уравнения (1) находят .
(а) (б) (в)
Рис. 1. Геометрия задачи: (а) зондирование сплошного СППВ; (б) измерение ДН в сплошном СППВ; (в) дискретный СППВ, зондирование и измерение ДН
На втором этапе вращают испытуемую антенну, измеряют принимаемый ею сигнал и восстанавливают ДН антенны , решая аналогичное уравнение свертки
. (2)
Зонд следует перемещать по поверхности S, охватывающей рабочую зону: R0 > Lа/2. Действительно, усматривая в равенстве (1) совпадение касательных составляющих электрических полей[1] облучающей волны и поля СППВ, на основании теоремы о единственности решения уравнений Максвелла можно утверждать, что в случае цилиндрической облучающей волны поле СППВ точно воспроизводит облучающее поле внутри S. При этом вне S, как показывают расчеты, поле имеет едва ли предсказуемую структуру.
На практике было бы выгодно в качестве вспомогательной антенны использовать источник сферической волны, например, небольшой рупор. При этом теорема единственности побуждает переходить к замкнутой поверхности S и телесному СППВ. Но результаты расчетов свидетельствуют о том, что если речь идет об измерении ДН антенны малого поперечного (относительно плоскости измерений) размера, то и в этом случае аппроксимация облучающего поля плоским пучком дает приемлемые результаты.
Вычислительный алгоритм. Аппроксимация облучающего поля
В отсутствии аналитического решения интегральное уравнение сводят к системе линейных алгебраических уравнений (СЛАУ). Простейший прием (едва ли уступающий методу Галеркина-Ритца) ? заменить СППВ дискретным пучком из N плоских волн и осуществлять измерения в M > N точках (рис. 1(в)). При этом с учетом кардиоидной ДН зонда уравнение (1) трансформируется в систему M уравнений относительно N неизвестных An = A(цn)
, (m = 1..M), (1')
где цn = (n ? 1) Д, ц'm = (m ? 1), Д = 2р/N, д = 2р/M, и решается, стандартным приемом минимизации среднеквадратичного отклонения.
В моделируемой ситуации имеем:
,
где ? расстояние от источника в точке D до зонда (см. рис. 1(а)) и .
Интересно выяснить, что за поле создается СППВ, аппроксимирующим исходящую из точки D цилиндрическую волну? На рис. 2 для четырех значений числа волн пучка Nпредставлены нормированные амплитудные рельефы этого поля при следующей геометрии системы: R0 = 7.5 и D =25. Полоска слева на рис. 1(а) ? это единая для всех рисунков палитра близкая к цветовой гамме географической карты: от темно синего цвета (нулевой уровень) до коричневого (единичный уровень). Рабочая область S отмечена окружностью радиуса R0, что позволяет судить о масштабе изображений. Справа вверху каждого рисунка укрупненно воспроизведено поле внутри S, нормированное к локальному максимуму в S. Среднеквадратичные отклонения поля пучка E(r,) от истинного облучающего поля E0(r,) в рабочей области S приведены там же.
(а) N = 90 (б) N = 115 (в) N = 130 (г) N = 145
Рис. 2. Амплитудные рельефы полей СППВ: цилиндрическая волна из точки D
Значения N выбраны с таким расчетом, чтобы проиллюстрировать характерные ситуации. Рис. 2(б) представляет ситуацию с рациональным числом волн пучка N = 115. При этом линии уровня с визуальной точностью соответствуют цилиндрической волне, и значение свидетельствует о высокой точности аппроксимации не только амплитуды, но и фазы облучающего поля. Рис. 2(а) соответствует явно недостаточному числу волн пучка, что проявляется в изрезанности рельефа поля на периферии области S, хотя в ее центральной части явственно просматривается цилиндрическая волна. Рис. 2(в) и 2(г) соответствуют другой крайности ? излишне большому числу волн пучка. При этом устойчивость СЛАУ (1') снижается настолько, что возникает асимметрия поля относительно оси 0у, разительно проявляющаяся при N = 145. Эта асимметрия ничем другим не объяснима, кроме как ограниченной точностью вычислений и представления чисел с плавающей запятой.
Для тех же ситуаций на рис. 3 представлены спектры {A()} пучков, аппроксимирующих облучающую цилиндрическую волну из точки D. Интервал углов ограничен «передним полупространством» . Жирные линии отображают амплитуды, тонкие ? фазы волн.
Рис. 3. Комплексные амплитуды {An} волн пучка.
На рис. 3(а), 3(б) (N = 90, N = 115) ожидаемо превалируют волны лучей из окрестности D. При N = 130 (рис. 3(в)) вопреки физике растет интенсивность поперечных волн ( вблизи 0° и 180°) и появляется асимметрия, резко обостряющаяся в ситуации N = 145 (рис. 3(г)). За пределами некоторого углового сектора в окрестности направления = 90° на точку Dкомплексные амплитуды соседних волн противофазны (скачки фазы на р), что и порождает «рябь» амплитудного рельефа, ярко проявляющуюся за пределами области S, а при малом числе волн (рис. 2(а)) даже на ее периферии.
Ясно, что неумеренный рост числа N ухудшает устойчивость СЛАУ (1') из-за уменьшающегося отличия соседних волн на окружности R0 и, следовательно, малого отличия соседних столбцов матрицы СЛАУ. Но происходит это настолько катастрофично, что мы, авторы, потратили массу времени и сил на упорные поиски ошибки в разработанной ими Delphi-программе, интерфейс которой показан на (рис. 4).
Рис. 4. Интерфейс программы исследования алгоритма восстановления ДН
диаграмма направленность антенна
Для выяснения истины были выполнены тестовые расчеты в среде Mathcad, результаты которых (рис. 5) подтвердили состоятельность нашей программы.
Рис. 5. Пространственные спектры, рассчитанные в среде Mathcad
Естественно, для неустойчивых решений играют роль и формат представления чисел с плавающей запятой, и возможные отличия алгоритмов решения СЛАУ, что приводит к отличиям в деталях, но структурно сопоставляемые решения одинаковы. В области умеренных значений N < 128 решения устойчивы и симметричны. Синий цвет подчеркивает гармоничность этих решений. Их максимальные значения |An|, не подвергнутые нормировке (рис. 5), имеют значения на уровне 0,08.
При N > 130 для нашей программы и N > 135 в случае Mathcad-программы возникает алогичное асимметричное решение (черный цвет). Причем норма решения еще сохраняется на прежнем уровне, что неожиданно, и причина этого заслуживает отдельного исследования. Дальнейшее увеличение числа волн пучка вызывает резко возрастающую асимметрию решения, связанную с его неустойчивостью, закономерным признаком которой служит резкое увеличение нормы решения (на семь порядков при N = 145). Красный цвет подчеркивает катастрофичность этой ситуации, проявляющаяся и в том, что поле пучка за пределами области S на порядки интенсивнее, чем внутри S (рис. 2(г)).
Облучение цилиндрической волной связано с обременительной необходимостью использовать в качестве вспомогательной антенны линейную антенную решетку или линейный раскрыв (пилбокс, например) достаточно большого габарита, ориентированные по нормали к плоскости измерений. Соображения физического толка дают надежду на то, что, если обследуемая антенна имеет малый поперечный размер, то облучающая сферическая волна может быть представлена с помощью СППВ без большого ущерба для воспроизведения реакции испытуемой антенны на облучающую волну. В практическом же плане использование вспомогательной антенны в виде рупора небольших волновых размеров для формирования сферической облучающей волны, чрезвычайно привлекательно.
На рис. 6 приведены поля, возникающие при аппроксимации в области S сферической волны, исходящей из точки D. Геометрические соотношения прежние: R0 = 7.5, D =25. Сигнал зонда в поле сферической волны определялся очевидной формулой с учетом обозначений, поясненных на рис. 2. При изменении числа Nволн СППВ проявляются те же закономерности, что и в случае цилиндрической волны, за тем исключением, что точность аппроксимации сферической волны с помощью плоского пучка волн существенно скромнее.
(а) N = 90 (б) N = 115 (в) N = 130 (г) N = 145
Рис. 6. Амплитудные рельефы полей СППВ: облучающая сферическая волна