CC BY
77
РАДЮЕЛЕКТРОШКА ТА ТЕЛЕКОМУШКАЦП
РАДИОЕЛЕКТРОНИКА И ТЕЛЕКОММУНИКАЦИИ
КАБЮ ЕЬЕСТКОШСБ АКБ ТЕЬЕСОММиШСАТЮЫБ
УДК 621.391.8
О. С. Антропов, С. М. Вовк, В. Ф. Борулько
МЕТОД МИНИМУМА ПРОСТРАНСТВЕННОЙ ПРОТЯЖЕННОСТИ ДЛЯ ОПРЕДЕЛЕНИЯ ДИАГРАММЫ НАПРАВЛЕННОСТИ АНТЕННЫ ПО ИЗМЕРЕНИЯМ БЛИЖНЕГО ПОЛЯ
Рассмотрена задача определения диаграммы направленности антенны путем восстановления распределения токов в плоскости раскрыва антенны по данным ближнего поля. Для решения обратной задачи восстановления распределения токов применен метод минимума пространственной протяженности. Диаграмма направленности антенны получена путем преобразования Фурье поля в раскрыве. Представлены результаты численного моделирования для случая одномерной апертуры.
ВВЕДЕНИЕ
Одной из важных и актуальных задач теории антенн является задача определения диаграммы направленности антенны (ДНА) по результатам измерения ее ближнего поля. На практике этот вид измерений наиболее распространен, поскольку позволяет минимизировать искажающее влияние окружающей среды и дает возможность оценить ДНА посредством преобразования Фурье. Однако при практическом использовании результатов измерений ближнего поля возникают трудности, обусловленные конечным размером области измерений. Так, при использовании стандарт© Антропов О. С., Вовк С. М., Борулько В. Ф., 2008
ной процедуры пересчета ближнего поля в дальнее поле, часть ближнего поля, которая не охвачена областью измерений, обычно полагается равной нулю [1]. Это приводит к неточному восстановлению распределения поля в дальней зоне вне некоторого диапазона углов, определяемого размерами антенны и области измерений, из-за появления в расчетах так называемой ошибки усечения. Для уменьшения ошибки усечения необходимо либо увеличить размер области измерений, либо ввести ограничения на физический размер исследуемых антенн. Это обуславливает необходимость поиска методов, позволяющих обойти указанные выше трудности.
Для решения указанной задачи было предложено несколько методов. Один из методов [2] использует априорную информацию о геометрии антенны и данные измерений ближнего поля для оценивания эквивалентного распределения магнитных токов на плоскости раскрыва антенны с последующим получением дальнего поля. Другой метод [3] использует дополнительную математическую обработку результатов измерений
ближнего поля на основе геометрических характеристик антенны. Идея этого метода основана на экстраполяции измеренных данных ближнего поля за пределы области измерений. Это позволяет получить эквивалентную область измерений, превышающую физически доступную. Так как данный метод основан на построении квадратичного функционала для модуля невязки измеренного поля и вычисляемого поля, то полученное решение будет неустойчивым и для получения устойчивого решения используется техника регуляризации Тихонова [4]. После восстановления поля в ближней зоне, оба метода используют прямой пересчет распределения поля в ближней зоне в распределение поля в дальней зоне с помощью стандартных алгоритмов преобразований, основанных на алгоритмах быстрого преобразования Фурье. Существенным ограничением указанных и других аналогичных методов является необходимость располагать точными априорными сведениями о геометрических характеристиках исследуемой антенны. Если же такие априорные сведения являются неточными или отсутствуют, то применение указанных методов становится проблематичным.
В [5] мы предложили новый подход к обработке результатов измерений ближнего поля. Идея этого подхода заключалась в замене стандартной техники регуляризации на новую технику, которая основана на критерии минимума пространственной протяженности решения. Применение этого подхода позволяет достичь существенного улучшения результатов восстановления поля точечных источников, расположенных в раскры-ве антенны, при отсутствии априорных сведений об их местоположении. В данной работе предлагается использовать этот подход для решения задачи определения ДНА по результатам измерений ближнего поля антенны. Ключевым моментом предлагаемого подхода является предположение о финитности пространственного распределения токов в раскрыве антенны, представленного набором точечных источников при отсутствии априорных сведений об их местоположении.
1 ОСНОВНЫЕ РЕЗУЛЬТАТЫ
Представление распределения поля в дальней зоне посредством преобразования Фурье от пространственного распределения токов в раскрыве антенны [6] является вычислительно эффективным и не содержит в себе погрешности, обусловленные ошибкой усечения. Следовательно, успешность решения задачи определения ДНА определяется результатом решения обратной задачи восстановления токов в раскрыве антенны. Для демонстрации общего подхода к решению данной обратной задачи используем простую модель, которая описывает измеренные данные 6(4) ближнего поля ан-
тенны и распределение токов а(г) в области 5 апертуры антенны в рамках скалярной задачи уравнением:
-гк\4 - г\
6(4) = | а(г)йг + N(4); 4е М, (1)
5о
где М - область измерений, N(4) - аддитивный пространственный шум, 5о - неизвестная область распределения токов в области 5, т. е. 5о с 5. Последнее означает, что если г е 5о, то а(г) ^ 0, но если г £ 5о, то а(г) = 0. Феноменологическая постановка обратной задачи восстановления токов в области апертуры антенны формулируется следующим образом: найти распределение токов а(г) по известным значениям поля 6(4) в условиях, когда эти известные данные искажены шумом и область 5о априорно неизвестна.
Отметим, что из (1) следует принципиальная невозможность точного определения а(г) в условиях ограниченной области измерений М и наличия шума. Действительно, источники тока всегда ограничены в пространстве и создают поле во всем пространстве. Но если только часть этого пространства доступна для непосредственных измерений, то, во-первых, происходит потеря части полезной информации и, во-вторых, доступная полезная информация является искаженной из-за наличия измерительного шума.
Ясно, что потерянная информация может быть восстановлена только искусственно на основе априорной информации, объем которой должен быть адекватен объему потерянной информации. Априорная информация, которая явно присутствует в модели (1), заключается в том, что искомое решение а(г) должно описываться финитной функцией, так как область 5о априорно ограничена. Теоретически такая априорная информация является достаточной для восстановления функции а(г) по любому известному сегменту данных 6(4). Однако наличие шума и отсутствие знания или неточное знание области 5о не позволяют эффективно применять соответствующие известные схемы восстановления. В этих условиях мы предлагаем решать уравнение (1) таким образом, чтобы пространственная протяженность области 5о, где сосредоточены источники, была минимальной.
Рассмотрим случай, когда источники тока расположены вдоль оси ОХ, а измерения ближнего поля выполняются вдоль линии, отстоящей от источников на расстоянии г, причем линия измерений параллельна оси ОХ и лежит в плоскости у = о. Тогда можно переписать (1) в виде:
-гЦ(4 - х)2 + г
6(4) = | а(х)е-—=== йх + N(4). (2) 5о 4(4 - х)2 + г2
10
1бо7-3274 «Радюелектрошка. 1нформатика. Управлшня» № 1, 2оо8
Поскольку задача решения уравнения (2) относительно а(х) является некорректной обратной задачей, то для получения устойчивого решения нужна процедура регуляризации. Вместо стандартной квадратичной регуляризации, основанной на добавлении члена в виде энергии решения, используем нелинейную регуляризацию, добавляя неэнергетический член в форме, описывающей пространственную протяженность решения [5]. Это приводит к формулировке следующей задачи минимизации:
F (u) = — f a (x )ej dx, (5)
2nJ
S
где u = k sin 0 = 2(n/X) sin 0; 0 - угол между направлением излучения и внешней нормалью к апертуре антенны. После этого результат можно перевести в пространство углов с учетом только тех пространственных частот, которые соответствуют диапазону [-п/2, п/2].
b(4) - f a(x)'
-ik>i(£, - x)2
V(4 - x)2 + z2
- dx
d^ +
- y2 fy[a(x);a, P, ...]dx ^min,
a (x)
(3)
где у2 - параметр регуляризации. Первый член в (3) является мерой ошибки полученного решения по области измерений; второй член определяет пространственную протяженность решения в качестве регуля-ризирующего члена. Отметим, что функция у, представленная в (3), может иметь различные варианты [7]; здесь используем следующую форму представления:
5(x);a, Р] = [\s(x)|2/Л2 + a2]в - a2p; 0 <р< 0, 5,
(4)
где параметры а, р позволяют управлять поведением функции у. Результаты численного моделирования показали [5, 7, 8], что наиболее удачным является такой выбор: а « а; р =1/16, где а - стандартное отклонение измерительного шума.
Решение задачи (3) будем называть решением по методу минимума пространственной протяженности (РММПП). Из (3) следует, что если для всех х выполняется: а2 >> |а(х)|2, то РММПП стремится к решению, получаемому при обычной регуляризации Тихонова с использованием энергетического члена. Кроме того, увеличение у2 до очень больших значений приближает РММПП к нулевому решению, а использование малых у2 приводит к решению без регуляризации. Таким образом, можно утверждать, что для каждого значения дисперсии шума существует некоторое
2
наилучшее значение параметра регуляризации у , при использовании которого абсолютная ошибка решения будет минимальной.
После решения задачи (3) с целью восстановления функции а(х), описывающей распределение токов антенны, соответствующая ДНА определяется через преобразование Фурье:
2 ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ
Цель проведения численного моделирования заключалась в сопоставлении результатов определения ДНА по измерениям ближнего поля при применении двух методов регуляризации решения обратной задачи: квадратичной регуляризации Тихонова [3, 4] и регуляризации по методу минимума пространственной протяженности. Для моделирования исходное распределение токов было задано в виде набора точечных источников, расположенных на оси ОХ в заданных узлах дискретной сетки, которая имела равномерный шаг Ах. Предполагалось, что измерение комплексной амплитуды ближнего поля было проведено вдоль отрезка прямой линии с шагом А^, отстоящей от линии источников на расстояние г. Учитывая (2), для численного моделирования была использована следующая система уравнений:
-i2n
bp =
У[ q(Ax/X ) -
p ( А£,/X ) ] 2 + (z /X ) 2
q = i
Хд/[ q(Ax/X) - p(A^/X)]2 + (z/X)2
p = i7P,
-N
p
(6)
где Ьр - известные данные, представленные комплексными значениями, Р =100; Ах/Х =0,1; А^/Х =0,1; г/Х = 2; Q = 100, причем амплитуды точечных источников ад =1 для д = 40, 50, 60, 70, 80, и ад =0 для других д, Ыр - комплексный аддитивный гауссовский шум.
Полученные результаты численного решения обратной задачи (6) с использованием квадратичной регуляризации Тихонова и регуляризации по методу минимума пространственной протяженности приведены на рис. 1. На рис. 2 приведены графики соответствующих ДНА, которые были рассчитаны посредством дискретного преобразования Фурье и приведены к диапазону углов [-п/2, п/2]. При численном моделировании, выполненном для 100 различных реализаций шума при отношении сигнал/шум +20 дБ, средний квадрат ошибки расчета ДНА по методу стандартной регуляризации Тихонова относительно истинной ДНА составил величину 6,87 %, тогда как для регуляризации по методу
2
S
S
a
q
V / i / \
V v t /
2 4 6 8 x/j
Рисунок 1 - Пространственное распределение токов:
пунктирная линия - стандартная регуляризация; непрерывная линия - предлагаемая регуляризация
90 -60 -30 0 30 60 в, град
Рисунок 2 - Восстановленные ДНА:
пунктирная линия — стандартная регуляризация; непрерывная линия - предлагаемая регуляризация
минимума пространственной протяженности эта величина составила менее 0,27 %, причем при применении предлагаемого метода наблюдалось визуальное совпадение полученной ДНА с истинной ДНА.
Таким образом, можно заключить, что предложенный метод позволяет уменьшить ошибку расчета ДНА по сравнению со стандартным подходом.
ВЫВОДЫ
Решение задачи определения ДНА целесообразно выполнять путем восстановления распределения токов в плоскости раскрыва с последующим преобразованием полученного распределения токов в значения ДНА. Данный подход обеспечивает возможность получения оценок ДНА в произвольном диапазоне углов, так как
при этом не возникает ошибка усечения данных из-за ограниченности области измерений ближнего поля. Применение метода минимума пространственной протяженности обеспечивает возможность решения обратной задачи при отсутствии априорных сведений о геометрии излучающей системы, причем применение функционала минимума пространственной протяженности в качестве регуляризирующего члена имеет преимущество по сравнению с традиционным подходом на основе квадратичной регуляризации. Результаты численного моделирования для случая одномерной апертуры подтвердили работоспособность и эффективность предложенного подхода.
ПЕРЕЧЕНЬ ССЫЛОК
1. Yaghjian A. D. An overview of near-field antenna measurements // IEEE Transactions Antennas Propagation. -1986. - Vol. AP-34. - Pp. 30-45.
2. Petre P., Sarkar T. K. Near-field to far-field transformation using an equivalent magnetic current approach // IEEE Transactions Antennas Propagation. - 1992. -Vol. 40, No. 11. - Pp. 1348-1356.
3. Bolomey J.-C, Bucci O. M, Casavola L, D'Elia G, Mig-liore M. D., Ziyyat A. Reduction of truncation error in near-field measurements of antennas of base-station mobile communication systems // IEEE Transactions Antennas Propagation. - 2004. - Vol. AP-52, No. 2. -Pp. 593-602.
4. Тихонов A. H., Арсенин В. Я. Методы решения некорректных задач. - М.: Наука, 1986. - 286 с.
5. Вовк С. М., Борулько В. Ф. Регуляризашя даних антен-них вим1рювань на основ! обмеження просторовоТ протяжное^ струмш. // Вюник Днтропетровського нацю-нального ушверситету. - 2006. - № 2/3. - C. 9-11.
6. Минкович Б. М., Яковлев В. П. Теория синтеза антенн. -М: Сов. Радио, 1969. - 269 с.
7. Вовк С. М., Борулько В. Ф. Метод минимума длительности для восстановления финитных сигналов // Известия вузов. Радиоэлектроника. - 1991. - Т. 34, № 9. - С. 66-69.
8. Вовк С. М., Борулько В. Ф. В1дновлення рефлектограм за допомогою екстраполяцп широкосмугових НВЧ-ви-м1рювань методом м1н1муму тривалост // Радюелект-рошка. ¡нформатика. Управлшня. - 2005. - № 1. -С. 5-8.
Надшшла 5.03.2008 Шеля доробки 9.04.2008
Розглянута задача визначення д1аграми спрямованост1 антени шляхом в1дновлення розпод1лу струм1в у площит розкриву антени за даними ближнъого поля. Для роз-в'язку оберненоЧ задач1 в1дновлення розпод1лу струм1в використано метод мгтмуму просторовоЧ протяжност1. Д1аграма спрямованост1 антени отримана за допомогою перетворення Фур'е поля у площиш розкриву. Представлен резулътати чиселъного моделювання для випадку одновим1рноЧ апертури.
The problem of far-field pattern determination from near-field data measurements through restoration of currents distribution on the aperture plane is considered. Extension minimum method is applied for currents estimation inverse problem. Fourier transform applied afterwards to determine far-field pattern. Simulation results are presented for one-dimensional aperture case.
12
ISSN 1607-3274 «Радюелектрошка. 1нформатика. Управлшня» № 1, 2008
