Особенности итеративной регуляризации в обратных задачах рассеяния

УДК 621.396.6

ОСОБЕННОСТИ ИТЕРАТИВНОЙ РЕГУЛЯРИЗАЦИИ В ОБРАТНЫХ ЗАДАЧАХ РАССЕЯНИЯ

А.Б. Хашимов

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

Ключевые слова: задачи рассеяния, математическая модель, функциональные уравнения, итеративные регуляризирующие схемы.

Обратные задачи рассеяния формулируются следующим образом [1]: по заданной диаграмме направленности (ДН) определить или форму рассеивателя, или распределение поверхностного тока для идеально проводящего рассеивателя известной формы, или поверхностный импеданс рассеивателя. Очень близки по постановке к обратным задачам рассеяния методы восстановления амплитудно-фазового распределения (АФР) излучающих систем по результатам измерений в ближней зоне [2]. Этот класс задач относится к некорректно поставленным, что вызывает значительные трудности при решении практических задач антенной техники, радиолокации и радионавигации. При этом для произвольно заданной ДН или углового распределения эффективного поперечника рассеяния (ЭПР) решение обратной задачи существует не всегда, а если существует, то не зависит непрерывно от исходных данных [1]. Кроме того, для обратных задач характерна статистическая постановка, связанная с тем, что результаты практических измерений сопровождаются случайными погрешностями, вызванными соответствующими точностными характеристиками измерительной аппаратуры зондирования электромагнитных полей, неточным определением координат при позиционировании зонда. Это означает, что численные методы решения обратных задач должны сопровождаться статистическими оценками полученных результатов.

Рассмотрим следующую задачу: пусть известна форма идеально проводящего рассеивателя, заданная как поверхность £, удовлетворяющая условиям непрерывности Ляпунова - в любых точках р, q , принадлежащих £, существуют векторы нормалей пр, п , для которых выполняется неравенство

\пр - < с|р - Ч\ “,

где с, 0 < а < 1 - произвольные числа. Будем считать, что на произвольной вспомогательной поверхности задано множество точек с координатами (хп, уп, zn ), п = 1: N , координаты точек определены со случайными ошибками (Дхп, Дуп, &п ) : (хп, уп, zn ) = (хеп, уеп, ) + (Дхп, Дуп, ) ,

где (хеп,уеп,) - точные координаты. В точках (хп,уп,zn ) заданы АФР составляющих вектора напряженности электрического поля Ехп, Еуп, Еп , также определяемые со случайными ошибками: (Ехп, Еуп, Еп ) = ( Еуп, Еуп, Кп ) + (ДЕхп, ДЕуп, ДЕ,п ) , где ( Еуп, Еуп, Кп ) - точное значение

поля. Введем координаты точек на поверхности рассеивателя £: (х'п,у'п,zn) = (х'^у,у'^у,z'У) + + (Дх'п, Дуп, ^п), п = 1: N, где необходимо восстановить значения АФР поверхностного тока (]хп, ]уп, Лп ) . Составляющие поверхностного тока связаны с заданными составляющими вектора

напряженности электрического поля строгими электродинамическими соотношениями следующего вида:

Z

Ех = — Г

х Иг •!

Г г2 + г2 ^

2БГ + ^-^Вг г

V У

гхгу г г

- і -^LD - і -xzD

7у г ^2 г

Г Г

Е = ^ Г ik {

У

гхгУ П

- ]х --------^ + ]у

Г

( г2 +г2 ^

2Бг + ^--------z-Dr

г

V У

г г

- h—D)

2 г

г

Е2 = — Г 2 ik {

гг ГГ

- і rxrLD - і -y-zD + і

Ух г Уу ^г ' J 2

Г Г

( г2 + г2 ^

2 Бг +^---------------У-Dr

г

V

У

где k = 2л/Я - волновое число; 2С = ^/ца/ва - волновое сопротивление среды; ва, ца - абсолютные диэлектрическая и магнитная проницаемости среды; гх = х - х '; гу = у - у'; гг = г - г ';

2 1--------------------------------------2 о. п 1 Г Іk 1 Ї „ 1

г = Л г + гу + г2 ; функции Бг = —!— + — I ф; Dг = —

4л V г

( іЛ

4л

3ік

г 2

ф; ф = ехр {-ікг)/г .

Таким образом, задача восстановления АФР поверхностного тока по заданному АФР вектора напряженности электрического поля является типичной некорректной задачей, сформулированной в виде системы интегральных уравнений I рода. Строго говоря, система (2) относится к классу функциональных уравнений, так как области определения искомых и заданных функций не совпадают. Кроме того, численное исследование системы (2) затруднено влиянием случайных ошибок при определении функциональных зависимостей. Для формирования полной математической модели (ММ) рассматриваемой задачи необходим выбор метода дискретизации системы (2), приводящего к системе линейных алгебраических уравнений (СЛАУ). Этот метод должен учитывать еще одну специфическую особенность ММ, связанную с необходимостью использования квадратурных формул высокой точности, например, Гаусса-Лежандра с числом узлов не менее 16 и выше. Это требуется для обеспечения необходимой устойчивости численного решения, в значительной степени зависящей от случайных факторов задачи. Поэтому любое повышение точности вычислительных процедур в конечном итоге приводит к стабилизации математического моделирования. Без ограничения общности будем считать, что размерности множеств точек задания

АФР (Ехп, Еуп, Егп ) и определяемых АФР токов (]хп, ]уп, ]2п ) совпадают. Обсуждение вопросов

построения ММ с недоопределенными или переопределенными СЛАУ выходит за рамки настоящей статьи. Широко используемый метод моментов для решения интегральных уравнений приводит к следующей формулировке ММ задачи восстановления (выбор базисных и весовых функций соответствуют методу коллокации):

N /

I ( і:

Е =

хп

т=1

N /

Е = ^ і і ктп

уп хт ху *

1^тп  і^тп

' + ]2тК2х'

ут ух

+і

тп тп

' + І2тК’

ут уу

т=1

N

Еп = 1(

і Ктп + і

^ хт х2' ^

ут у 2

т=1

кт.п + ітКт); п = 1: N,

где коэффициенты СЛАУ определяются численным интегрированием соответствующих функциональных зависимостей ядра системы уравнений (2) по поверхностям конечных элементов. Методы дискретизации поверхности £ для построения оптимальной структуры конечных элементов рассмотрены в [3]. Необходимо отметить, что для трехмерных (3D) задач большой электрической размерности СЛАУ может иметь очень высокий порядок. В этом случае компьютерные ресурсы являются естественным ограничением возможностей математического моделирования обратных задач рассеяния в строгой электродинамической постановке. Для преодоления та-

х

3

3

г

г

г

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

Вместе с тем на практике достаточно часто встречаются задачи, в которых асимптотические представления полей излучения в дальней зоне для 3D рассеивателей соответствуют полям излучения в плоскости, перпендикулярной оси эквивалентных двумерных (2D) рассеивателей - случай Е -поляризации [4]. Например, система вибраторных антенн вертикальной поляризации, расположенных вблизи протяженных цилиндрических объектов (антенны радио- и телевещания, базовых станций мобильной связи, установленные на мачтах и башнях). Тогда система (2) значительно упрощается и может быть представлена одним уравнением [5]. Будем считать, что ось г совпадает с продольной осью рассеивателей, тогда

где L - контур, образованный сечением рассеивателей плоскостью, перпендикулярной их оси;

уравнения (4) дает возможность прогнозировать поведение АФР токов для тех рассеивателей, влияние которых может приводить к искажениям поля излучения. Уменьшить такое влияние можно, например, изменением расположения антенной системы относительно рассеивателей. В частности, учет влияния поясов и опор башни (мачты) в общей структуре поля излучения антенной системы дает возможность синтеза АФР токов, формирующих ДН в наибольшей степени приближенной к заданной.

В идеализированной постановке, когда отсутствуют погрешности определения поля излучения, координат точек, в которых находятся АФР токов, численное решение уравнения (4) может быть получено как прямыми, так и итерационными методами, которые имеют свои особенности и области применения. Для задач большой размерности широкое применение получили итерационные методы, например, минимизации обобщенных невязок, сопряженных градиентов. Используемый программный комплекс MATLAB предоставляет разработчикам ММ широкие возможности для составления программ с использованием эффективного представления матричных форм Function Handle, а рациональная организация блочной структуры СЛАУ с использованием неявных вложенных циклов дает возможность значительно повысить быстродействие программ. Важным достоинством итерационных методов является возможность пошагового контроля скорости сходимости и нормы невязки решения. Кроме того, предварительные численные эксперименты для однотипных задач могут дать практический интервал значений нормы невязки, при которых достигаются как удовлетворительная точность решения, так и приемлемое время расчетов.

Значительно расширяет возможности математического моделирования обратных задач рассеяния использование суперкомпьютеров. Большое количество вычислительных экспериментов, проведенных на кластере «СКИФ-Урал» Южно-Уральского государственного университета, показывают, что полный цикл моделирования электродинамических задач, включающий в себя

формирование СЛАУ размерностью порядка 1,3 -109 элементов и решение прямым методом исключения Гаусса не превышает 3 минут. Использование итерационных методов решения СЛАУ позволяет увеличить размерность задачи в десятки раз, но при этом быстро нарастает время решения.

Особенностью уравнений (2), (4) является ярко выраженная неустойчивость численного решения, характерная для некорректных задач. Это можно подтвердить следующим численным экспериментом. Пусть контур L представляет собой окружность радиуса a, ka = 12. Идеально проводящий рассеиватель возбуждается плоской волной. На первом этапе будем считать, что все исходные данные определяются точно, то есть все погрешности отсутствуют. Численное решение ИУ (4) методом коллокации определяет АФР тока на контуре L, при этом точки на-

функция Ганкеля второго рода нулевого порядка. Численное решение

блюдения и точки источников принадлежат L . С помощью этого тока в заданных точках вне L можно определить АФР поля рассеяния Ez . Будем считать, что это поле является исходным для решения обратной задачи рассеяния. В отсутствии погрешностей относительное изменение нормы токов прямой и обратной задачи не превышает величины 10_6, что позволяет считать эту задачу устойчивой. На втором этапе введем случайные изменения амплитуды и фазы заданного поля Ez . Без ограничения общности будем считать, что случайные изменения имеют равномерный закон распределения. С помощью генератора случайных чисел rand и процедуры масштабирования введем аддитивные случайные изменения амплитуды и фазы с относительной погрешностью 10_8. Норма тока в этом случае имеет порядок 106 относительно той же величины для устойчивой задачи. Если относительные случайные изменения значений поля имеют

_3

величину порядка 10 , что характерно для большинства практических методов повышенной точности измерений в сверхвысокочастотном диапазоне, то норма восстановленного тока имеет

порядок 1011. Следовательно, рассматриваемая задача восстановления АФР тока может быть решена только с использованием методов регуляризации.

Широкое применение для решения некорректных задач получили итеративные методы регуляризации [6]. Для детализации постановки задачи запишем уравнение (4) в операторной форме:

Az = u, z е Z, u eU.

Введем метрические пространства X , Y и отображение G: X ^ Y, заданное на подмножестве Dg с Y . Требуется по элементу x е DCj найти его образ G (x) е Y . В терминах операторного уравнения (5) G = A-1, X = U, Y = Z и задача состоит в определении A_1. В этом случае Dg = AZ с U. Для корректной задачи отображение G определено на всех X и непрерывно. Если вместо элемента x е Dg задано его неточное значение x5 е X , то вводится норма разности элементов x5,x в метрике X pX (x5,x)<5 , где 5 - число, характеризующее погрешность задания x5 . Будем требовать, чтобы численный метод решения уравнения (5) удовлетворял следующему основному требованию: чем точнее задано x5 (чем меньше 5), тем точнее можно определить x. Основываясь на результатах работы [6], сформулируем регуляризирующий алгоритм следующим образом: определить оператор R, который ставит в соответствие паре (x5,5) элемент Z5 е Z , причем Z5 ^ z в метрике Z при 5 ^ 0. Это означает, что R (x5, 5) имеет смысл приближенного решения задачи при заданном наборе входных данных. Введем следующие итеративные схемы, отвечающие свойствам оператора R (x5, 5) [6]:

ный оператор. Как показано в [6], схема (6) применяется для ограниченного оператора А, а схема (7) - в случае неограниченного замкнутого оператора А. В частности, если А - интегральный

пользование схем (6) и (7) при I = 0,1,... позволяет получить схемы регуляризации I = 0,1,... порядков.

Итеративные регуляризирующие схемы должны быть дополнены правилами останова, предоставляющие широкие возможности для коррекции параметра ц в зависимости от скорости сходимости итерационного процесса и погрешности задания исходных данных. Правило остано-

где Ц - параметр итеративной регуляризации; A - оператор, сопряженный к A; E - единич-

оператор Фредгольма из W2(l) (Q) и (Q) (что соответствует рассматриваемой задаче).

, то ис-

ва регламентирует число итераций n (g) так, чтобы

Характерной особенностью схем (6) и (7) является отсутствие четких рекомендаций по выбору оптимального параметра ц. Большое количество вычислительных экспериментов численного решения задачи (4) с помощью схем (6), (7) позволяют предложить следующую методику определения оптимального параметра ц (необходимо отметить, что этот параметр в значительной степени зависит от величины 5, поэтому эта методика не может считаться универсальной, но для задач с незначительно отличающейся геометрией и величинами 5 получаемые оценки ц могут быть использованы как стартовые для задач с другими входными значениями).

1. Вводится серия параметров ц1, ц2,..., цк . Для рассматриваемой задачи восстановления АФР тока по неточно заданным значениям поля в некоторой области D эксперименты показывают, что эти значения могут находиться в интервале 103 .4 • 104.

2. Для этой серии параметров ц последовательно применяется одна из итеративных схем для одинакового числа итераций т .

3. Для каждого параметра ц из серии находится невязка квазирешения ||гт - zm_Ц .

4. Из серии вычислительных экспериментов находится оптимальный параметр цор ^ тт||zm _ zm_1\\, например, методом интерполяции кубическими сплайнами.

5. Для цор{ проводится полная прогонка итеративной регуляризирующей процедуры с вве-

денным правилом останова, если до его выполнения не произошло прерывание процесса при условии достижения заданной точности восстановления - zг_11|/zmax < Д .

Рассмотрим практическую реализацию итеративной регуляризирующей схемы (6) для уравнения (4). Пусть D - двумерная область, в которой заданы неточные значения Ег (х, у), (х, у)е D . Тогда операторы в (6) принимают вид:

А*А = |Л(х;,У;)К(г,г;)^; К(г,г;) = 161я01)(кг)н02)(кг..)<*; (х;,у;)еL;

L 16 Б

А*и = 4 | Ez

4 Б

г = >/(х_х')2 +(У_у')2; г; = ^(х_х;)2 +(У _У;)2.

В соответствии с выбранным правилом дискретизации Б итеративная регуляризирующая схема записывается следующим образом:

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

Рассмотрим пример, реализующий предложенный подход для решения обратной задачи восстановления. Для контура L в виде окружности радиуса a, ka = 12 исходным для решения обратной задачи рассеяния является поле рассеяния плоской волны на идеально проводящем цилиндре. Область D представляет собой окружность радиуса R = 4а . Введем следующие параметры аддитивного случайного процесса, интерпретирующие неточно заданные исходные данные: амплитуда amp = 0,1Ez max, фаза phase = 5°. На рис. 1 представлены результаты восстановления АФР тока с помощью итеративной регуляризирующей схемы (6). Число итераций для раз-

_3

личных статистических выборок при относительной точности восстановления 10 не превышает 200, оптимальный параметр ц = 4540 . Для сравнения приведены результаты строгого решения задачи определения АФР тока при отсутствии погрешностей исходных данных.

Рис. 1. Амплитудные распределения токов рассеивателя: сплошная линия - строгое решение; • • • - восстановленное АФР

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

Большой практический интерес представляет вопрос о границах предложенного метода восстановления, то есть, при каких ошибках в исходных данных можно получить численные результаты, отвечающие достоверным АФР токов и ДН полей рассеяния в рамках выбранных критериев точности. Этот вопрос можно решить с помощью серии вычислительных экспериментов, в частности, для рассмотренного примера срыв достоверного восстановления наступает при amp = 0,34Ez max , phase = 25°.

Рис. 2. Амплитудные ДН в дальней зоне: сплошная линия - строгое решение; • • • - восстановленное поле

Сравнение этих результатов иллюстрирует высокую эффективность итеративной регуляри-зирующей процедуры для неточных (зашумленных) исходных данных.

Выводы

1. Для решения обратных задач восстановления АФР токов по неточно заданным исходным данным для напряженности электрического поля принципиально необходимо использование итеративных регуляризирующих схем [6]. Применение прямых методов в этом случае приводит к недопустимо высоким ошибкам восстановления. Переход к 2D задачам особенно эффективен для случая E -поляризации, для этих задач значительно расширяются возможности математического моделирования сложных антенных систем.

2. Сравнение полученных результатов с известными тестовыми решениями показывает высокую эффективность и универсальность предложенных математических моделей некорректных задач.

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

Литература

1. Колтон, Д. Методы интегральных уравнений в теории рассеяния / Д. Колтон, Р. Кресс. -М.: Мир, 1987. - 312 с.

2. Бахрах, Л.Д. Радиоголография в микроволновой технике / Л.Д. Бахрах, А.П. Курочкин. -М. : Советское радио, 1979. - 320 с.

3. Салихов, Р.Р. Суперкомпьютерное моделирование полей рассеяния на объектах сложной формы /Р.Р. Салихов, А.Б. Хашимов // Вестн. Юж.-Урал. гос. ун-та. Сер. «Компьютерные технологии, управление, радиоэлектроника». - 2013. - Т. 13, № 1. - С. 55-60.

4. Войтович Н.И., Хашимов А.Б. О соответствии асимптотических решений двумерных и трехмерных задач в антенной технике / Н.И. Войтович, А.Б. Хашимов // Радиотехника и электроника. - 2010. - Т. 55, № 12. - С. 1471-1476.

5. Галишникова, Т.Н. Численные методы в задачах дифракции / Т.Н. Галишникова, А. С. Ильинский. - М. : Изд-во Московского университета, 1987. - 208 с.

6. Бакушинский, А.Б. Итеративные методы решения некорректных задач / А.Б. Бакушин-ский, А.В. Гончарский. - М. : Наука, 1989. - 128 с.

Хашимов Амур Бариевич, канд. физ.-мат. наук, доцент кафедры конструирования и производства радиоаппаратуры, Южно-Уральский государственный университет (г. Челябинск); xab@kipr.susu.ac.ru.

Bulletin of the South Ural State University Series “Computer Technologies, Automatic Control, Radio Electronics”

2013, vol. 13, no. 4, pp. 78-85

FEATURES OF ITERATIVE REGULARIZATION FOR INVERSE SCATTERING PROBLEMS

A.B. Khashimov, South Ural State University, Chelyabinsk, Russian Federation, xab@kipr. susu. ac. ru

This paper proposes the mathematical models of restoring the current distribution on surface of perfect conducting scatterer in case inaccurate input data of the electromagnetic field. The rigorous electrodynamics formulations which lead to the ill-posed functional

equations of I kind are used. It is shown that the use of iterative regularization schemes provides a stable numerical solution under the certain criteria.

Keywords: scattering problems, mathematical model, functional equations, iterative regularization schemes.

References

1. Colton D., Cress R. Metody integralnikh uravneny v teorii rasseyaniya [Integral Equations Methods in Scattering Theory]. Moscow, World, 1987. 312 p.

2. Bakhrah L.D., Kurochkin A.P. Radiogolografiya v mikrovolnovoii tekhnike [Radiohologram Method in Microwave Technique]. Moscow, Soviet Radio, 1979. 320 p.

3. Salikhov R.R., Khashimov A.B. Supercomputer Simulation of the Scattered Fields on Complex Shape Objects [Superkomp’yuternoye modelirovaniye v zadachakh rasseyaniya na ob’ektakh slozhnoy formy]. Bulletin of the South Ural University, Series “Computer Technologies, Automatic Control, Radio Electronics”, 2013, vol. 13, no. 1, pp. 55-61.

4. Voitivich N.I., Khashimov A.B. On the Correspondence of Asymptotic Solutions to 2D and 3D Problems in Antenna Engineering [O sootvetstvii asimptoticheskikh resheniy dvumernykh i trekhmer-nykh zadach v antennoy tekhnike]. Journal of Communications Technology and Electronics, 2010, vol. 55, no. 12, pp. 1374-1379.

5. Galishnikova T.N., Il’insky A.S. Chislennye metody v zadachakh difraktsyi [Numerical Method for Diffraction Problems]. Moscow, Publishing Center MSU, 1987. 208 p.

6. Bakushinsky A.B., Goncharsky A.V. Iterativnyi metody resheniya necorrectnikh zadach [Iterative Methods for Ill-Posed Problems]. Moscow, Science, 1989. 128 p.

Поступила в редакцию 25 июля 2013 г.