Научная статья на тему 'Двойственная регуляризация в одномерных обратных задачах рассеяния'

Двойственная регуляризация в одномерных обратных задачах рассеяния Текст научной статьи по специальности «Математика»

CC BY
231
49
Читать
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕОДНОРОДНЫЕ СРЕДЫ / МНОГОСЛОЙНЫЕ НАНОСТРУКТУРЫ / ДИЭЛЕКТРИЧЕСКАЯ ПРОНИЦАЕМОСТЬ / ПРОВОДИМОСТЬ ЗЕМНОЙ КОРЫ / РАССЕЯННОЕ ПОЛЕ / НЕКОРРЕКТНЫЕ ОБРАТНЫЕ ЗАДАЧИ / ДВОЙСТВЕННАЯ РЕГУЛЯРИЗАЦИЯ

Аннотация научной статьи по математике, автор научной работы — Гайкович Константин Павлович, Гайкович Петр Константинович, Галкин Олег Евгеньевич, Смирнов Александр Ильич, Сумин Михаил Иосифович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Гайкович Константин Павлович, Гайкович Петр Константинович, Галкин Олег Евгеньевич, Смирнов Александр Ильич, Сумин Михаил Иосифович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
Предварительный просмотр
i Надоели баннеры? Вы всегда можете отключить рекламу.

DUAL REGULARIZATION IN ONE-DIMENSIONAL INVERSE SCATTERING PROBLEMS

A method is proposed to restore subsurface permittivity profiles of one-dimensional inhomogeneous media by scattered-field measurements based on the Lagrange formalism in the solution of nonlinear ill-posed problems. This method enables the solution of problems in their natural statements for initial Maxwell equations. The method efficiency is demonstrated in numerical simulation of the retrieval of Earth’s crust conductivity profiles and the profiles of permittivity inhomogeneities in multilayer periodic structures by spectral measurement data.

Текст научной работы на тему «Двойственная регуляризация в одномерных обратных задачах рассеяния»

РАДИОФИЗИКА

УДК 621.371+537.87

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

© 2013 г. К.П. Гайкович1, П.К. Гайкович1, О.Е. Галкин2, А.И. Смирнов3, М.И. Сумин2

Институт физики микроструктур РАН, Н. Новгород 2Нижегородский госуниверситет им. Н.И. Лобачевского 3Институт прикладной физики РАН, Н. Новгород

gai@ipm. sci-nnov.ru

Пт1тупило кр1докцию 18.07.2012

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

Ключевые слова: неоднородные среды, многослойные наноструктуры, диэлектрическая проницаемость, проводимость земной коры, рассеянное поле, некорректные обратные задачи, двойственная регуляризация.

Введение

Необходимость решения одномерной обратной задачи рассеяния возникает в самых разных методах электромагнитной и акустической диагностики одномерно-неоднородных сред. Впервые такая задача была сформулирована и решалась А.Н.Тихоновым [1] для магнитотеллурического зондирования земной коры со слоистой структурой глубинного профиля проводимости по данным многочастотных измерений электромагнитных полей на ультранизких частотах, способных проникать на несколько километров в толщу земли [2]. На классе непрерывных функций эта задача решалась в [3-5] путем ее сведения к нелинейному интегральному уравнению (типа Фредгольма или Вольтерра) 1-го рода и использования итерационной схемы, начиная с борновского приближения, методом обобщенной невязки Тихонова [6]. Исследования показали, что такой подход неэффективен не только для сильных возмущений проводимости среды, но и в случае слабых неоднородностей теряет свою сходимость после 2-3 итераций.

Более сложные проблемы возникают при диагностике неоднородного полупространства по угловой и спектральной зависимости коэффициента отражения [7-10], в частности разнообразных многослойных структур [8-10]. В данном случае даже, когда борновское приближение обладает достаточной точностью [10], возможны сильные искажения, связанные с ошиб-

кой линеаризации при выделении вклада возмущения в измеряемый коэффициент отражения по мощности.

Таким образом, возникает необходимость иных подходов к решению нелинейных обратных задач. Применительно к геомагнитному зондированию в [9] был предложен метод двойственной регуляризации, алгоритм которого в самой простой форме (градиентной минимизации невязки в функции Лагранжа) исследовался в [10] для сильных возмущений профиля проводимости. Этот метод, основанный на лагранже-вом подходе в общей теории оптимизации и оптимального управления [11-16], позволяет решать сложные обратные задачи непосредственно в рамках их естественных дифференциальных постановок, без перехода к нелинейному интегральному уравнению.

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

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

При решении большинства обратных задач физической диагностики [4] часто переходят от постановки в терминах дифференциального

уравнения (обыкновенного или в частных производных) к постановке, в которой основным объектом является интегральное уравнение. Как правило, для линейных задач это линейные интегральные уравнения первого рода типа Фред-гольма или Вольтерра, для решения которых имеется широкий арсенал эффективных вычислительных алгоритмов, например, таких, как метод Тихонова, метод регуляризации на компактах и ряд других [6].

В то же время для нелинейных обратных задач процедура перехода от дифференциальной постановки к ее «интегральному аналогу» является сама по себе самостоятельной и подчас весьма сложной проблемой. Для перехода к интегральным уравнениям часто приходится использовать те или иные приближения (например, борновское). При этом неизбежно теряется точность и искажается первоначальная, соответствующая сути дела, физическая модель.

Существуют методы, являющиеся составной частью общей теории методов оптимизации и оптимального управления, которые позволяют решать обратные задачи непосредственно в рамках их естественных дифференциальных постановок. К ним относится рассматриваемый в данной работе метод, основанный на теории двойственности, или, другими словами, метод двойственной регуляризации [14,16]. Он является устойчивым к ошибкам исходных данных, ориентирован на решение нелинейных оптимизационных задач и в своей основе опирается на одноименный метод [12,13,15,17], разработанный для решения линейно-выпуклых оптимизационных задач. Его особенностью является то, что он, во-первых, развивает классический в теории оптимизации лагранжев подход; во-вторых, служит удобной универсальной основой для решения различных обратных задач физической диагностики и, в-третьих, позволяет решать такие обратные задачи в их естественных дифференциальных постановках. Отметим, что обратные задачи непосредственно в рамках их естественных дифференциальных постановок позволяют решать и методы, основанные на минимизации нелинейного функционала невязки (см., например, [18]), которые с формальной точки зрения являются частными случаями двойственных методов.

В качестве общей характеристики рассматриваемого подхода следует отметить, что двойственные алгоритмы являются одними из самых эффективных при решении оптимизационных задач с ограничениями [19, 20], в которых непосредственным образом используется классическая идея «снятия» ограничений, заложенная в принципе Лагранжа. Отметим также, что ха-

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

По-видимому, первым двойственным алгоритмом в теории оптимизации явился предложенный в 1958 г. так называемый алгоритм Удзавы [21] (см. также [19]), заключающийся в непосредственном решении на основе градиентного метода задачи, двойственной к исходной оптимизационной задаче, и быстро завоевавший популярность. Вопросы его обоснования рассматривались в [22, 23]; при этом предполагалось, что 1) исходные данные известны точно; 2) функция Лагранжа задачи должна иметь седловую точку. Вместе с тем оба эти предположения являются весьма ограничительными, так как, во-первых, для практических задач характерно именно наличие ошибок в исходных данных и, во-вторых, для таких задач, как правило, само доказательство существования или несуществования седловой точки представляет собой нелегкую математическую задачу. Формальное же применение алгоритма в общей ситуации может привести и приводит к стандартным эффектам неустойчивости приближенного решения [13,15,17].

В 1968 г. работы М. Хестенса и М. Пауэлла [24, 25] дали толчок дальнейшему интенсивному развитию основанных на двойственности методов решения нелинейных оптимизационных задач с ограничениями. Многие подробности, связанные с двойственными методами решения нелинейных оптимизационных задач математического программирования в конечномерных пространствах на основе идеологии [24, 25], которые называют также методами множителей Лагранжа в условной оптимизации, а также соответствующую библиографию можно найти в [18, 26]. Основной конструкцией этих методов служит так называемая модифицированная функция Лагранжа, представляющая собою сумму обычной классической функции Лагранжа и взятого с некоторым положительным весом штрафного слагаемого, которое может иметь различный вид. В силу нелинейности оптимизационной задачи двойственные методы, основанные на конструкциях модифицированной функции Лагранжа, носят локальный характер, и об их сходимости речь можно вести только при условии, что начальное приближе-

ние как по основной, так и по двойственной переменной в этом процессе находится достаточно близко к искомой оптимальной паре этих переменных [18, 26]. Важнейшими условиями этой сходимости являются некоторые априорные предположения об исходных данных оптимизационных задач, выражающиеся в терминах

V.» и * и

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

Дх) ^ тт, g(x) = С?!(х),...,gm (х)) = ° х е Ят (1) с дважды непрерывно дифференцируемыми функциями /, gi, /=1..., т, одним из основных априорных условий локальной сходимости метода множителей Лагранжа в варианте монографии [26, стр. 111-120] является линейная независимость градиентов g1,(x*),..., ¿т (х*). Наряду с другими важнейшими априорными предположениями, на основе хорошо известной для этой конечномерной задачи теории необходимых и достаточных условий оптимальности второго порядка и классического свойства компактности ограниченного замкнутого множества в конечномерном пространстве (подробности см. на стр. 111-120 в [26]), это предположение приводит к следующему неравенству при некотором е > 0 и при некотором достаточно большом коэффициенте штрафа с > 0:

/(х*) < /(х) + (Я*, g(х)) + c|g(х)|2 = (2) = Lc (х, Я*), х е Sе (х*), х ^ х*, где (у) - скалярное произведение, Я* е Ят -некоторый вектор, 5”6 (х*) = {х е Ят : |х - х* | < е},

L6 (х, Я) - модифицированная функция Лагранжа. Последнее неравенство говорит о том, что вектор Я является обобщенным вектором Куна-Таккера в смысле работы [16] и что модифицированная функция Лагранжа, взятая в точке Я , совпадающей с обобщенным вектором Куна-Таккера, принимает строго минимальное значение в искомой оптимальной точке х . Это обстоятельство является определяющим, и именно оно и лежит в основе метода множителей Лагранжа в теории конечномерной оптимизации, предложенного впервые в [24, 25] и развитого в дальнейшем для конечномерных оптимизационных задач в огромном числе различных публикаций (см., например, библиографию в [18, 26, 27]).

К сожалению, в задачах бесконечномерной оптимизации, к классу которых принадлежат и рассматриваемые здесь нелинейные обратные задачи поиска минимальных по норме решений,

описанная выше ситуация принципиально меняется. Связано это прежде всего с тем, что в бесконечномерных пространствах, и в частности в бесконечномерных гильбертовых пространствах, несправедливы многие важнейшие факты, лежащие в основе схемы, описанной выше для случая конечномерных задач математического программирования. Например, в бесконечномерном гильбертовом пространстве ограниченное замкнутое множество не обязательно является компактом. Другая неприятность, приводящая к принципиальным затруднениям в обосновании приведенной выше «конечномерной» схемы, связана с отсутствием адекватной теории условий оптимальности второго порядка для задач бесконечномерной оптимизации, аналогичных (1):

f (x) ^ min, g (x) = 0, x e X, (3)

где X, H - гильбертовы пространства, f: X ^ R

- некоторый функционал, g : X ^ H - некоторый оператор. Оказывается, однако, что, несмотря на отсутствие этих важных обстоятельств в ситуации бесконечномерных задач (3), их свойство, выражаемое неравенством

f (x*) < f (x) + (^\g(x)) + c||g(x)f = (4) = Lc (x, Я* ),x e Ss (x*), x ^ x*

(где Я* e H - некоторый элемент, Ss (x ) = = {x e X: |x - x* | < s}), аналогичным неравенству (2) для конечномерного случая, является достаточно характерным. И по этой причине именно оно и может служить тем априорным предположением, на котором основывается «бесконечномерная» двойственная схема [14, 16]. Характерность этого свойства обусловлена тем фактом (подробности см. в [16]), что для параметрической задачи

f (x) ^ min, g (x) = h, x e X (5)

(очевидно, задача (3) входит в семейство задач (5) при h=0) в любой окрестности любой точки h, для которой эта задача разрешима, обязательно имеются другие точки h', такие, что для соответствующих задач (5) при h = h' выполняется неравенство

f (xh'*) < f (x) + (Я>(x) - h') +

2 (6)

+ c||g(x) - h'|| , x e Ss (xh ), x Ф xh ,

где xh * - решение задачи (5) при h = h'. Можно

также утверждать, что таким же характерным свойством (подробности см. в [16]) для задачи (3) является выполнение неравенства

f (x*) < f (x) + (Я*, g(x)) + (7)

+ c||g(x)||, x e Ss (x* ), x Ф x ,

отличающееся от (4) только типом штрафного слагаемого. Именно по этой причине ниже в статье, как и в [14, 16], используется комбинированное штрафное слагаемое вида с(||g(х)|| +

II и 2

+ ^(х) ). Важнейшими особенностями «нелинейного» двойственного алгоритма [14, 16] конструирования минимизирующей последовательности в оптимизационной бесконечномерной задаче являются учет неточного задания исходных данных и применение для решения двойственной задачи процедуры регуляризации по Тихонову. Регуляризация при максимизации двойственной вогнутой функции является необходимой составляющей двойственного алгоритма [14, 16], так как неустойчивость при решении бесконечномерных оптимизационных задач, в том числе и выпуклых, на основе стандартных градиентных процедур является хорошо известным фактом оптимизационной теории [19]. Подчеркнем, что в известных нам публикациях [24, 25] (см. также библиографию этих работ), посвященных «конечномерному» методу множителей Лагранжа, рассматривались лишь задачи с точными исходными данными и процедуры какой-либо регуляризации для них не применялись.

Поле в одномерно-неоднородной среде.

Сведение уравнений Максвелла к интегральным уравнениям

В среде с неоднородным распределением комплексной диэлектрической проницаемости в (г) = вг +81(г) в весьма общем случае неоднородности 81(г) в многослойной среде со значе-

ниями вг в г-м слое распределение комплексных амплитуд векторов гармонического [ ~ ехр(-/ш t) ] электрического E и магнитного H полей определяется комплексной амплитудой плотности электрического тока источника j из уравнений Максвелла:

Ух E = — H, (8)

с

_ тт . ш _ 4л .

Ух H + г—sE = — ] , (9)

с с

где с - скорость света, ш - циклическая частота. Вклад неоднородности может быть выделен в качестве эффективного источника

гш

} ^ =----е1 (z)E путем перенесения этого вклада

11 4л

в правую часть (9):

^ тт . ш . ш 4л. 4л.. . ,.т

Ух H + г—в г E = -г — SlE + — j = — 0е// + j). (10) с с с с

Далее задача может быть сведена к интегральному уравнению [28] с использованием функций Грина для многослойной среды, которые получаются путем разложения поля источника по плоским волнам с использованием формализма входных импедансов:

E(r) = E0 (г) + ВД = E0 (г) -

- О (г ', г Н(г' )E(r ' )^', (11)

V

Ео (г) = | О (г ', г)](г ' )Ог'.

V

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

Для случая одномерно-неоднородной среды в^) = вг +в1( z) решение аналогичных прямой и обратной задач существенно упрощается. Рассмотрим такую многослойную среду в полупространстве z > 0, зондируемую плоской электромагнитной волной

Е0 (г) = Е0еХР(^хХ + ЯуУ + ,

К = >/ kl2 - ^2 - ^ = 7 kl2 - ^,

падающей на границу полупространства из

вакуума (kl2 = (—)2). Пусть для определенности

с

волновой вектор этой волны лежит в плоскости ху, так что поле падающей волны на границе может быть представлено как

^) = Е; (Zl) + El0(Zl) =

= ^[сова уо + вта (^-Хо -^Zо)Г (12)

kl kl

где а - угол, определяющий соотношение между компонентами горизонтальной и вертикальной

22

поляризации, так что < Е 0(^)| >а= Е0 . Тогда

зондирующее поле в 1-м слое многослойной среды Е0(z')и отраженное структурой (без учета непрерывной неоднородности) поле Ед (zl)на ее границе могут быть вычислены с

использованием формализма метода входных импедансов [28]:

E0 (z') = E¿ + E0 = E0 [cosa G1 (Zi, z')y0 +

+ sina G"(zi,Z’)(kzLx0 -kLz0)] = (i3)

kl kl

= E01 (z')cosayo + E0l(z')sina(kix0 - z0),

kl kl

ER (Zi) = Ef1 (Zi) + E„'"(Zi) = Eo [Recosa y 0 +

к k (i4)

+ sim^+i (-~ZLx0--xz0)], ki ki

z ) - z

— Tll-líaikz-1 (Z'-Z‘ ^ in 1 a-2ikz,ldl -ikz,l (z'-zl К

G'u (zx, z') _ ЭД1 (e z-zt +

z 1° + z,

z (2) - z 7-

p||,l _ in 1 T 1 _ 1 T i

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

RN+1 - _(2) + _ > Ли~ 11 7(i Zi« ^ L\ i

T¡" _J^ r¿ — ДГЇ

le 1 с A A

^ 7 (i) + 7

Z. tz. ,k .d.

in_i - ikz,i di

z (i) + z i—2 ¿in ^i-1

01 T1 H _ |01 1 T z in + zieikz idi

0/ Vo; i—2 zi>j) + zi-1

z(k) _ zink +1) - izk tg(kz,kdk ) z z(N+1) _ z

zin . (k+1) , Й \ k , zin zN+1 :

zk - izin )tg(kz,kdk )

k, i. 1 k„ , ю

c

zk _,—, z1_^4 k0 _-, dw+1 _ 0.

z _

k0 0 k

k

2o

J01( z,)

k.

(z')

7 2—11 11 0, (z)

kz x0 + kxkzz0 — kxkz z"

o + o1(z)

+ (15)

-Ey (z')[k2 yo ]-

+ Ez (z')[+kxk2 Xo + klz0 - k2^% ,o

o + o1(z)

zo Uaz'-

o1 (z)

E0z (z)z0 •

+ smaR1,n+1 (-kzLXo - ^zo)] + k1 k1

(16)

r i k2

+ E0cosaJo1 (z')--j^-G1 (z1,z')G1 (z',z1)dz'yo -

(kiz kx2 2o k„

- En sinaJ o1 (z ') —— (—--)G!I (z1, z ') x

x G"(z ',z^dz '(-kz1 Xo -^zo), k1 k1

G11,1 (z ', z1) _

тнд Jl ,1

1 dII,i d II,i ^>2ikz,idi

1 Rl,1 Rl, N+1 Є

x (e-ikz,1 (zM-z) + r\

N1 z(l-1) - z

pll,1 in •

r1,1 “ T(l-1)

,1 e_ikz’l (z-zl-1 -2dl))

N +1 Є ),

l-1 (i) +

l T'1 I T ^in ^i Jkz,idi

zl + z

T1_ ГЇ

in

z ^i) + z

i_1 ^in ^^i+1

Т!^?ГЇ

0 J-1 z (i) + z

o; 1 T z. t z. ik .d-

711 in i _e kz>iai

0 A A z ^i) + z &1 i_1 ^/+1

(l) _ zin~1) - izl tg(kz,lal) (1)

zK¿ _

zl - izí-1)tg(kz,ldl )

z , z;' _ z.

k

k

Полное поле от среды создается суммой поля, отраженного основной структурой, и рассеянного поля, которое представляется одномерным аналогом формулы (11):

E(z) = ЕД (z) + Е1 ^) = ЕД (z) +

, г ^ехр{±(z-

zil =-^^—, г,1 = —^—, dl = 0.

1 ^ в/ 1 kг,/ 1

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

<|Е(г1)|2 >= |[|Е0111 (г) + )|2 +

(17

22 + |Е1 (г1)\ + |Е1( г1)| +

+ 2Яе(Е0К ;* (г1 )Е; (г1)) + 2Яе(Е0К|1* (г, )Е1 (гх))],

)

где < > - усреднение по углу а . Из (17) получаем полный коэффициент отражения структуры:

в + в,(г)

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

Е(г,) = Ед(г,) + Е,(г,) = Е0[со8аД^Уо +

R _-

2 _ |R0| + \R1\ +

2

En

(18)

x < Re[(ER1* (z1 )E1 (z)) + (ER"* (z1 )E1 (z1))] >

'2 Ь ,1 |2,|„|| I2-

где = -(|^1,N+1 + Д1,N+1 ), а звездочкой

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

x

x

2

Ar = \R\2 - |R0|2 = — < Re[(ERx* (z, )EX (z,)) +

E>| (19)

+ (ER11* (Zi)El( Zi))] >+| RJ

для которого имеем

(20)

Аг _ |Rj|2 + Re{Jo¡(z')-0[—R¡N+1G1 (z¡,z')

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

z z

k2 - k2

x G1 (z', z1)-Z—-- R¡'*n +1G" (z¡, z')G" (z', z¡)]dz'}.

kz

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

Задача геомагнитного зондирования• Для УНЧ-СНЧ-диапазонов радиоволн на земной поверхности выполняются граничные условия Леонтовича, в силу чего можно считать волну распространяющейся вглубь по нормали к земной поверхности z = 0, а земную кору можно считать проводником с чисто мнимой комплексной диэлектрической проницаемостью o_o' + io"« i4гост / ю , которую определяет профиль проводимости ст = ct(z). Для определенности будем считать, что поле имеет компоненты E , Hx, и в дальнейшем опускаем их индексы.

В этом случае уравнения Максвелла предельно упрощаются:

d2E . 4roj(z)ra

■ + i -

-E = 0,

dz2

TT c dE

H = —i-------, zn < z < 0,

ш dz

Е(г = 0, ш) = Е0(ш), Н(г = 0,ш) = Н0(ш),

ш е [ш,,ш2]. (22)

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

Из (21) легко получить эквивалентное интегральное уравнение [4]:

4лгш

H(ш,z) = JH(ш,z”)dz"

a(z')dz'. (23)

Если в (23) положить г = 0, получаем интегральное уравнение, связывающее измеренную на поверхности компоненту поля Н0 = Н(ш,0) и профиль ст(г):

H0(ш) = ^J JH(ш,z ")dz"

ст(z")dz". (24)

Аналогичные соотношения можно получить и для электрического поля.

Ядром уравнения (24) служит интеграл от неизвестного распределения поля H(ra,z), и задача определения профиля ct(z) является нелинейной. Анализ начинается с естественного выбора однородного профиля проводимости <j(z) = Сто = const в качестве первого приближения. В этом предположении, если на поверхности такой однородной среды источником создается электрическое поле с распределением по частоте Е°(ш) = Е0°(ю), то ему соответствует вполне определенная частотная зависимость магнитного поля H°(ш) = H°>(ю), а соответствующие распределения полей в среде выражаются так:

z

. zN

E (ш, z) = E0 (ш )exp^ — i -) ,

Ö Ö

H °(ш, z) = C(-—1 Е0(ш)( z — i^) = ш0 Ö Ö

= H 00(ш)ехр^^ — iz)

ÖÖ

(25)

(26)

(21)

где ö =

толщина скин-слоя. Тогда от-

іша„

где гп <0 - нижняя граница интервала анализа, в котором заведомо находится область искомой неоднородности профиля проводимости, а измеренными считаются электрические и магнитные поля на поверхности г = 0 в диапазоне частот [ш,,ш2]:

клонение АН0 (ш) = Н0(ш) - Н00 (ш) от реально измеренной зависимости свидетельствует о неоднородности подповерхностного профиля проводимости. Таким образом, в первом приближении, которое по существу является борнов-ским, в ядре интегрального уравнения (24) используется Н 0(ш,г). В качестве ст0 естественно

2

С

0 01 ООН 0.012 0.013 ООН 0.015

Рис. 1. Слева: исходные «данные измерений» при среднеквадратичном значении случайной погрешности 5Н0 = 1 в зависимости от частоты f= ш/2л (в условных единицах); справа: жирная сплошная кривая - исходный профиль проводимости, точки - результат восстановления (1-я итерация), сплошная - результат восстановления (3-я итерация)

выбрать приповерхностное значение проводимости, которое полагается известным. Такой подход позволяет свести задачу к решению интегрального уравнения Фредгольма 1-го рода и решать его итерационным методом относительно ст(г):

АН

0

(ю) = ^ К' (ю, г ')Аст1+1(г ' )dz', (27)

Кг(ш,г ') = -Н (г ",ш)dz ", с *

-да

где ядро К пересчитывается на (г +1)-м шаге итерационного процесса по восстановленному профилю стг(г). Регуляризированное решение задачи (27) получается на каждом шаге с использованием алгоритма метода обобщённой невязки, разработанного в [28] для решения на классе комплекснозначных функций в гильбертовом пространстве W21 (пространство Соболева). В этом методе невязка минимизируется до уровня погрешности данных 8Н, вычисляемой в интегральной метрике ¿2[шш 2 ].

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

Интервал частот «измерений» выбирался из естественных физических соображений таким образом, чтобы толщина скин-слоя на минимальной частоте была больше слоя, в котором ищутся неоднородности проводимости, а на максимальной частоте зондировался бы тонкий поверхностный слой, в котором проводимость не отличается от поверхностной. В примере, представленном на рис.1, где восстанавливается относительно слабая (20% от невозмущенного уровня) неоднородность проводимости на глубине 1 км, при изменении частоты от 3 до 200 Гц толщина скин-слоя изменяется от 3 до 0.35 км.

«Данные измерений» АН 0(ш) рассчитывались для неоднородного профиля, показанного рис.1, при значении ст0 = 0.01 См/м и использовались в уравнении (9) для решения задачи при среднеквадратичном уровне моделируемой погрешности измерения мнимой части комплексной амплитуды 8Н 0 = 1. Качество решения возрастает в процессе итераций, однако при данном уровне погрешности результат практически не улучшается после 3-й итерации. Выполненный анализ показал, что для сильных неоднородностей проводимости, которые чаще всего встречаются в реальных условиях, рассмотренный итерационный алгоритм приводит к существенным погрешностям восстановления. Это объясняется тем, что ядро, используемое в первом приближении, очень сильно отклоняется от реального, и, как следствие, результат решения на значительной части интервала восстановления может не приближаться, а удаляться от точного решения, что приводит к раскачке неустойчивости в последующих итерациях.

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

Математическая постановка задачи осуществлялась в следующих предположениях:

1) Область искомой неоднородности профиля проводимости ст(г) заведомо находится в некотором конечном промежутке zn < z < 0 , где zn < 0 - нижняя граница интервала анализа. Поэтому при z < zn проводимость ст^) считается величиной постоянной.

2) Используется естественное предположение, что при z < zn величины полей Е и Н ограничены при любой фиксированной частоте ш.

Это позволяет использовать классы функций, задаваемых лишь на конечном промежутке [zn, 0], что облегчает обоснование сходимости метода регуляризации. Отметим также, что изучение исходной обратной задачи УНЧ-зонди-рования на бесконечном полуинтервале (-«,0] приводит формально к иной постановке оптимизационной задачи, в которой все равно приходится определять некий конечный интервал, за пределами которого измерения становятся неинформативными. В силу названных причин такая постановка в данной работе не рассматривается.

В принципе, возможны два варианта постановки обратной задачи УНЧ-зондирования.

В первом варианте мы считаем значения полей Е^ = zn,ш), Н^ = zn,ш) на нижней границе zn дополнительными вспомогательными неизвестными, подлежащими нахождению наряду с проводимостью ст^). Существенно, что при введении вспомогательных дополнительных неизвестных рассматриваемая задача оказывается линейной по этим вспомогательным переменным. Алгоритм, реализующий такой подход, описан нами в [7].

Во втором варианте постановки задачи мы считаем, что все неоднородности искомого профиля проводимости локализуются внутри промежутка восстановления [zn ,0], то есть за пределами этого промежутка проводимость является постоянной и равной поверхностной проводимости стп = ст^ = zn) = ст0, которая доступна измерению. Поскольку численный эксперимент проводился именно для такого варианта постановки задачи, то все дальнейшее изложение следует данному варианту. Итак, на интервале z < zn предполагаем, что ст^) = стп. Решение уравнения (21) на этом интервале и граничное условие при z = zn дают формулы (25), (26), поэтому обратную задачу геомагнитного зондирования в комплексной форме на промежутке

[гя,0] можно поставить следующим образом: по-

добрать профиль проводимости ст^), z е [zn, 0], так, чтобы при любой частоте ше[ш1,ш2] решения Е[ст,ш](z) и Н[ст,ш](z) системы (21) с краевыми условиями

Н ^ = z„) =0' - 1)Л/2лстп/ шЕ ^ = Zn ), (28)

Е[ст, ш](z = 0) = Е0(ш) (29)

удовлетворяли условию

Н [ст, ш]( z = 0) = Н0(ш). (30)

Переформулируем эту задачу в вещественной форме, введя новые обозначения:

х1 = Re Е, х2 = 1т Е, х3 = Re Н, х4 = 1т Н , х10 = Re Е0, х20 = 1т Е0, х30 = Re Н0, х40 = 1т Н0. Тогда система уравнений (21) примет вид:

dx, ю

—L =---------------x4, —

dz c dz

dx3 4roj(z)

dx2

dz

- = — x,

dx4

dz

c

4roj(z)

-x,,

x2,

или

= Л(ю,а( z)) x,

dz

(31)

(32)

x =

v x4 у

' 0 0 0 -ю

1 0 0 ю 0

с - 4тсст( z) 0 0 0

v 0 - 4тсст( z) 0 0

Соответственно, обратная задача геомагнитного зондирования в вещественной форме на промежутке [zn, 0] примет следующий вид: подобрать профиль проводимости ст^), z е [zn, 0], так, чтобы при любой частоте ш е [ш 1, ш 2 ] решения х1[ст, ш](z), х2[ст, ш](z), х3[ст, ш](z) и х4[ст, ш](z) системы (31) с краевыми условиями

,[ст,ш]^ = zn ) = у12%стп

X (-x1[ст,ш](z = zn)= zn)),

Х4[ст,ш](z = zn) = ^12гостп/ш X

х (х1[ст,ш]^ = zn) - ш](z = zn)),

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

\ х1 [ст, ш]( z = 0) = х10 (ш),

[х2[ст, ш](z = 0) = х20 (ш) удовлетворяли условиям

[х3[ст, ш](z = 0) = х30 (ш),

Х4[ст,ш](z = 0) = Х40(ш) ■

(33)

ю

c

c

x

2

x

3

Эту задачу можно сформулировать в виде эквивалентной задачи минимизации: минимизировать функционал с о

10 (ст) = || ст || 2= I ст2(z)dz ^ min (35)

J z„ ст

по всем тем функциям ct(z) из множества D = {ст е L2 [zn ,0]: 0 < ст(z) < M} (где M - достаточно большое положительное число), которые для любых ю е [ю,, ю 2 Удовлетворяют условиям (34), где x3[ct,ю](z) и x4[ct,ю](z) - решения системы (31) с краевыми условиями (33).

Можно доказать, что сформулированная задача минимизации имеет решение ct,(z) . В то же время нельзя гарантировать единственность оптимального элемента из-за нелинейности в ограничении этой задачи (нелинейность типа произведения управления на траекторию).

В случае линейных обратных задач, которые сводятся к эквивалентным линейно-выпуклым оптимизационным задачам, для их решения на основе двойственной регуляризации достаточно конструкции классической функции Лагранжа. В нелинейных случаях необходимо использовать конструкции модифицированных функций Лагранжа, то есть обычных функций Лагранжа с добавленными к ним штрафными слагаемыми того или иного вида. Рассмотрим здесь один из возможных случаев модифицированной функции Лагранжа сформулированной задачи:

L^ (ст, ю) = ||ст||2 + I [X,(ro)(x3[CT, ra](z = 0) - x30(ra)) +

+ Х ^ю^Дст, ю]( z = 0) - x40 (ro))]dro +

f ш V7 2

+ H-{iJ ^[ст, ю]( z = 0) - x30(ra)|2 draj +

m2 V 7 2

+ 1 | |x4[<j, ra](z = 0) - x40(ra)|2 draj + (36)

i“2 I 12

х4[ст,ю](z = 0) - х40(ю) d“}

%<“i

i“2 I 12

х3[ст,ro](z = 0) - х30(ю) d“ +

“i

или

Lц (ст,ш) = ||ст|| +

+ (^1(ш),(Х^^ ш ](z = 0) - х30(ш ))) +

+ (я 2 (ш), (Х4 [ст, ш](z = 0) - Х40 (ш))) +

+ ц{|| х3[ст, ш]( г = 0) - х30(ш)|| +

+ | |х4[ст,ш]( г = 0) - Х40(ш)|| +

II 1|2

+ х3[ст, ш](г = 0) - х30(ш)|| +

2

+ |Х4[ст,ш](г = 0) - Х40(ш)|| },

где Я = (Я1, Я2), ц > 0 . Смысл параметра ц со-

стоит в том, что при достаточно большом его

значении для любых Я(ш) = (Я1(ш),Я2(ш)) гарантированно существует минимум модифицированной функции Лагранжа Ьц (ст, Я) по ст(г).

Правда, элемент, на котором достигается минимум, вообще говоря, не единственный.

В соответствии с результатами [14, 15] запишем регуляризированную по Тихонову с параметром регуляризации а модифицированную двойственную задачу, представляющую собой задачу максимизации вогнутого функционала на гильбертовом пространстве

L^(ra1,ш2) = L2(ra1,ш2)хL2(ra1,ш2):

К (Я) = К (Я)-а||Я II2 -

- min L (ст, Я) - а || Я |

CTeD ^

(37)

где Я = (Я1,Я2) еL^(ю1,ю2), максимум берётся по Я из множества Лц = {Яе L2(ш1,ш2): || Я11 < ц} .

Алгоритм, реализующий решение этой задачи, описан нами в [7] и использован здесь в программной реализации для решения поставленной обратной задачи. Его исследования в численном моделировании по замкнутой схеме выполнялись аналогично описанной выше процедуре для метода итерационного решения интегрального уравнения. На рис. 2 представлен пример восстановления сильной неоднородности при том же уровне невозмущенной проводимости ст0 = 0.01 См/м, что и на рис. 1.

Можно видеть, что результат восстановления практически безупречен.

На рис. 3 представлен пример неоднородности в форме меандра, характерной для слоев земной коры.

Результат показывает, что восстановление профиля со скачком восстанавливаемой функции является более сложной задачей и качество восстановления оказывается хуже, чем для гладкого профиля. Тем не менее результат позволяет верно судить о форме проводимости и дает хорошую количественную оценку амплитуды неоднородности. На практике, если есть основания полагать, что мы имеем дело с однородным слоем, можно аппроксимировать результат решения ступенчатой функцией и повторить решение, используя эту функцию в качестве 1-го приближения.

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

Рис. 2. Слева: исходные «данные измерений» при среднеквадратичном значении случайной погрешности 5Н0 =

= 1 в зависимости от частоты; справа: 1 - исходный профиль проводимости, 2 - результат восстановления

Рис. 3. Слева: исходные «данные измерений» при среднеквадратичном значении случайной погрешности 5Н0 = = 1 в зависимости от частоты; справа: 1 - исходный профиль проводимости, 2 - результат восстановления

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

Задача восстановления профиля диэлектрической проницаемости периодически неоднородной среды по коэффициенту отражения

Рассмотрим существенно более сложную задачу восстановления неоднородности в периодической эпитаксиальной структуре - такой, которая используется в рентгеновских зеркалах [8], -по частотной зависимости коэффициента отражения по мощности в одном из брэгговских максимумов отражения. Неоднородность правильной периодической структуры с периодом d из двух веществ с параметрами е01и е02, связанная с процессом диффузии материалов при ее изготовле-

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

Ф) = е01 + (е02 - ^01);[(^, ;[(г) =;[(г + d). Будем считать, что за периодической структурой находится подложка большой толщины, в которой диэлектрическая проницаемость является постоянной. Используя уравнение (13) для Ы-слойной периодической структуры, такую задачу, как показано в [8, 9], можно свести к решению интегрального уравнения. Из (20) имеем:

I |2 і |2 І 12

Ar = R -\R0\ = Ш2 +

Г i k2

+ j f (z')Re{As—[ —Rt*+fi1 (zi, z')G1 (z', z,) -

Рис. 4. Слева: исходные «данные измерений» при среднеквадратичном значении случайной погрешности 5Н0 = 1 в зависимости от частоты; справа: 1 - исходный профиль проводимости, 2 - результат восстановления

- 2 _ - 2

■ ^(г', z1)] }dz'.

—_

Или, в компактной форме,

Аг (ю) = Щ2 +| f(z’)g (0, ю,г'^' . (39)

г'

Используя условие периодичности, получаем:

і 12 Аг (ю) = Щі| +

N _1

■ | f (г)]Г g[0, ю,г + (/ _ 1)а:

(40)

N _1

2

| f0 (г)Хg[0, ю, г + (і _ 1^]±

что позволяет свести задачу к решению интегрального уравнения

Аг(ю,0) = Ї f (г)К(0,ю,-

0

а

= | f \г)К (0, ю,

(41)

и применить итерационный алгоритм

Аг (0, ю) = | /(к) (7)^ (0, ю, ^ +1^1 (/к-1 (7))|2. (42)

а

Постановка обратной задачи для решения (42) возможна на основе как угловой, так и спектральной зависимости коэффициента отражения. В данной работе мы исследуем спектральную зависимость. Решение начинается с начального приближения /(0) = 0, либо выбранного с учетом имеющейся априорной информации или модели распределения неоднородности по периоду структуры. На каждом

шаге итераций уравнение (42) решается как интегральное уравнение Фредгольма 1-го рода методом обобщенной невязки [28].

Численное моделирование этого метода было выполнено для 50-периодной структуры Мо-Si, используемой в рентгеновских зеркалах [8, 9]. При решении полагалось, что период и толщины слоев основной структуры известны, поскольку они хорошо определяются по положению брэгговских пиков в угловой структуре коэффициента отражения в жестком рентгеновском диапазоне [8], тогда как влияние неоднородности на периоде структуры проявляется в частотной зависимости брэгговского пика в мягком рентгеновском диапазоне. Для рассматриваемой структуры этот пик расположен в спектральной области X =12.5^14.5 пт при угле падения 0 = 850 , на котором поляризация падающего излучения близка к горизонтальной.

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

На рис. 6, 7 представлены примеры численного моделирования по замкнутой схеме решения обратной задачи при моделируемом уровне случайной погрешности со стандартным отклонением 8г = 0.012. Значения диэлектрической проницаемости здесь и далее представлены для средней длины волны X = 13.3 нм из диапазона измерений.

Рис. 5. Функция ядра |к (X, z)| уравнения (42) в зависимости от длины волны на одном периоде структуры Мс-^

пт А, пт

Рис. 6. Слева: 1 - профиль основной структуры, 2 - исходный профиль проводимости, 3 - результат восстановления. Справа: 1 - коэффициент отражения идеальной периодической структуры без неоднородностей |^0 (Х)|2,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2 - исходные «данные измерений» |^(Х)|2 + 5г при среднеквадратичном значении случайной погрешности 5R2 = = 0.012 в зависимости от частоты, 3 - коэффициент отражения, вычисленный по восстановленному профилю, 4

I 12

— значения квадратичного слагаемого Щ1 (е1 (2)) , рассчитанные для первого приближения

z, пт Я, пт

Рис. 7. Обозначения аналогичны рис. 6

Результаты моделирования показывают, что характерные масштабы моделируемых неоднородностей 0.5 - 2 нм восстанавливаются с точностью 20-40%. При этом коэффициент отражения, вычисленный по восстановленному профилю, удовлетворяет условию невязки с заданной точностью. Вместе с тем, как видно из

рис. 7, значения ^(е^))^ первом приближении могут существенно превышать уровень по-

грешности измерений, что в ряде случаев приводит к расходимости алгоритма, причем в получаемом решении условие невязки может быть как выполнено, так и не достигаться. К такой же расходимости могут приводить и погрешности борновского приближения, в рамках которого получено решаемое уравнение. Для примера, представленного на рис. 6, первая поправка к борновскому приближению оказывается достаточно мала, что можно видеть на рис. 8.

z,nm\

12.5

Л,пт

Рис. 8. Поправка к борновскому приближению в

ядре уравнения (42) для примера на рис. 6

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

Тогда обратная задача формулируется следующим образом: подобрать профиль неоднородности структуры f (z) так, чтобы при любой

частоте 0e[Qj,ю2] решения системы (8), (9) удовлетворяли условию

Мf ](га) = 'Ч(га). (43)

Эту задачу можно сформулировать в виде эквивалентной задачи условной минимизации

функционала

il и 2

Io( f) = | |f|| ^ min МУ ](q) = Лг0(ю), f е L2(0, d) = D, ш£ [Qj, ю2].

Модифицированный функционал Лагранжа для этой задачи имеет вид:

L,(УД) - ||f|| 2+

(44)

+

Г Х(ш)(Дг[/](и) - Дг0 (ra))dra ■

J Wj

(45)

+^Ч .L(Дг[ / ](И) -ДоИ)2 d® +

+ f (М/](“) - А^0 (га))2 dra),

J Wj

а регуляризированная модифицированная двойственная задача -

К= V(X)-а||ХII2 -

- min L (/,Х)-а || X ||2 ^ max, (46)

(/,xo)eR4xL2(0,z„) Ц

ХеЛ, = {Xe L2(гаj,га2): ||Х11 ^,}-Далее схема решения повторяет алгоритм [7]. При измерениях отражения в зависимости от угла падения зондирующей волны в алгоритме нужно сделать соответствующие очевидные замены аргументов и переменных интегрирования. Все утверждения, сделанные выше, по поводу решения задачи УНЧ-зондирования справедливы и в данном случае.

Численный эксперимент

Была создана численная реализация описанного алгоритма метода двойственной регуляризации решения задачи. Его исследование в численном моделировании показало, что отсутствие информации о фазе отраженного сигнала приводит к серьезным трудностям в решении. Зачастую при градиентной минимизации модифицированного функционала Лагранжа (45) невязка легко минимизируется до заданного уровня погрешности, но, в отличие от задачи УНЧ-зондирования, соответствующий профиль может сильно отличаться от моделируемого, т.е. решение попадает в локальный минимум. В процессе максимизации функционала (46) двойственной задачи решение может быть выведено из таких локальных минимумов, но чтобы реализовать это преимущество двойственного подхода в полной мере, необходимо в численном моделировании подобрать оптимальную последовательность параметров итерационной схемы в алгоритме, пользуясь большой свободой их выбора [7], или использовать разумно выбранное первое приближение. Численное моделирование показало, что сходимость алгоритма существенно улучшается, если в качестве первого приближения используется, например, решение линеаризованной задачи, полученное с применением метода Тихонова. Этот способ введения дополнительной априорной информации, как правило, приводит к радикальному улучшению качества решения некорректных задач (см. примеры в [4]). Вместе с тем вряд ли возможно создать какой-либо автоматически работающий алгоритм для решения столь сложной нелинейной обратной задачи для произвольного профиля неоднородности - досто-

И

г, пт х, пт

Рис. 9. Метод минимизации невязки. Слева: 1 - профиль основной структуры, 2 - исходный профиль проводимости, 3 - результат восстановления. Справа: 1 - коэффициент отражения идеальной периодической структуры

без неоднородностей |Л0 (X)2, 2 - исходные «данные измерений» при 5г = 0, 3 - коэффициент отражения, вычисленный по восстановленному профилю

Рис. 10. Метод двойственной регуляризации. Обозначения аналогичны рис. 9, но при 5г = 0.012

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

На рис. 9 представлен пример моделирования решения задачи стандартным методом минимизации функционала невязки (4-е слагаемое в (45)) градиентным методом для заданного профиля неоднородности по рассчитанным значениям коэффициента отражения без набрасывания случайной погрешности 8г = 0. Результат демонстрирует, как решение в процессе минимизации попадает в локальный минимум (при уровне невязки 0.0040), где оно сильно отличается от исходного профиля, в то время как вычисленный по этому решению коэффициент

отражения с большой точностью (0.004) совпадает с заданным, так что эту разницу трудно различить на рисунке. На рис. 10 можно видеть результат восстановления профиля неоднородности методом двойственной регуляризации по вычисленным значениям коэффициента отражения, на которые набрасывалась случайная погрешность. Эта погрешность моделировалась случайными числами, имеющими нормальное распределение с таким же значением стандартного отклонения 8г = 0.012, как и на рис. 6, 7. Максимум функционала регуляризированной модифицированной двойственной задачи (46) достигнут при значении невязки 0.0118, весьма близком к уровню моделируемой погрешности. При градиентной минимизации функционала Лагранжа (45) использовалось естественное априорное ограничение 0 < / < 1.

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

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

Сравнивая результаты моделирования рассмотренных задач низкочастотного зондирования и диагностики неоднородностей многослойной структуры, можно видеть, что в последнем случае - в отличие от первого - метод Тихонова позволяет получить неплохие результаты, что является следствием лучшей применимости борновского приближения. Однако потеря информации о фазе рассеянного сигнала при измерениях коэффициента отражения по мощности создает трудно устранимый дефицит информации, что дает для решения много больше свободы попасть в ложный минимум. Вместе с тем результаты метода двойственной регуляризации и в этом случае надежнее, поскольку не зависят от непредсказуемых ошибок линеаризации задачи и точности борновского приближения.

Заключение

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

Представленные результаты получены по программе фундаментальных исследований Отделения физических наук РАН и поддержаны грантами РФФИ (№11-02-97060_р_поволжье, №12-01-00199, №13-07-97028_р_по-волжье), а также грантом Минобрнауки РФ в рамках государственного задания на оказание услуг в 2012-2014 гг. подведомственными высшими учебными заведениями (шифр заявки 1.1907.2011).

Список литературы

1. Тихонов А.Н. Об определении электрических характеристик глубоких слоев земной коры // Доклады АН СССР. Нов. сер. 1950. Т. 73. № 2. С. 295-297.

2. Велихов Е.П., Жамалетдинов А.А., Собчаков Л.А. и др. Опыт частотного электромагнитного зондирования земной коры с применением мощной антенны СНЧ-диапазона // Доклады РАН. 1994. Т. 338. № 1. С. 106-109.

3. Поляков С.В., Ермакова Е.Н., Поляков А.С., Якунин М.Н. Формирование спектров и поляризации фонового ультранизкочастотного электромагнитного шума на поверхности земли// Г еомагнетизм и аэрономия. 2003. Т. 42. № 2. С. 240-248.

4. Gaikovich K.P. Inverse problems in physical diagnostics. New York: Nova Science Publishers Inc., 2004.

5. Гайкович К.П., Смирнов А.И. Обратные задачи низкочастотного когерентного зондирования земной коры // Вестник Нижегородского университета им. Н.И. Лобачевского. 2007. № 2. С. 73-80.

6. Тихонов А.Н., Гончарский А.В., Степанов В.В., Ягола А.Г. Регуляризирующие алгоритмы и априорная информация. М.: Наука, 1983.

7. Гайкович К.П., Кутерин Ф.А., Смирнов А.И., Сумин М.И. Двойственная регуляризация в обратной задаче УНЧ-зондирования земной коры // Вестник Нижегородского университета им. Н.И. Лобачевского. 2009. № 1. С. 47-52.

8. Barisheva M.M., Gaikovich K.P., Gaikovich P.K. et al. Reflectometry sounding of inhomogeneities in periodic multilayer structures // Proceedings of 12th International Conference on Transparent Optical Networks (IC-TON 2010, Munich, Germany, June 27-July 1, 2010). P. Tu.P5 (1-4).

9. Gaikovich P.K., Sumin M.I., Gaikovich K.P. Onedimensional inverse scattering problem // Proceedings of 2011 13th international Conference on Transparent Optical Networks (ICTON 2011, Stockholm, Sweden, June 26 - 30, 2011). P. We.A2.4 (1-4).

10. Gaikovich K.P., Gaikovich P.K., Galkin O.E., Sumin M.I. Dual regularization in one-dimensional inverse scattering problem // Proceedings of 5th International Conference «Ultrawideband and Ultrashort Impulse Signals» (6-10 September, 2010, Sevastopol, Ukraine), Sevastopol: IEEE. Р. 90-92.

11. Бакушинский А.Б., Гончарский А.В. Некорректные задачи. Численные методы и приложения. М.: Изд-во Моск. ун-та, 1989.

12. Сумин М.И. Регуляризованный градиентный двойственный метод решения обратной задачи финального наблюдения для параболического уравнения // Ж. вычисл. матем. и матем. физ. 2004. Т. 44. №

11. С. 2001-2019.

13. Сумин М.И. Регуляризация в линейновыпуклой задаче математического программирования на основе теории двойственности // Ж. вычисл. матем. и матем. физ. 2007. Т. 47. № 4. С. 602-625.

14. Сумин М.И. Регуляризованный двойственный метод решения нелинейной задачи математического программирования// Ж. вычисл. матем. и матем. физ. 2007. Т. 47. № 5. С. 796-816.

15. Sumin M.I. Parametric dual regularization in a linear-convex mathematical programming // Computational Optimization: New Research Developments. Chapter 10. New York: Nova Science Publishers Inc., 2010. Р. 265-311.

16. Sumin M.I. Parametric dual regularization in a nonlinear mathematical programming // Advances in Mathematics Research. Volume 11. Chapter 5. New-York: Nova Science Publishers Inc., 2010. Р. 103-134.

17. Сумин М.И. Некорректные задачи и методы их решения. Материалы к лекциям для студентов старших курсов: Учебное пособие. Нижний Новгород: Изд-во Нижегородского госуниверситета, 2009. 289 с.

18. Короткий А.И., Цепелев И.А. Решение ретроспективной обратной задачи для одной нелинейной эволюционной модели // Труды Института математики и механики УрО РАН. Т. 9. 2. Динамика жидкости и газа. Екатеринбург: Изд-во УрО РАН, 2003. С. 73-86.

19. Мину М. Математическое программирование. Теория и алгоритмы. М.: Наука, 1990.

20. Васильев Ф.П. Методы оптимизации. М.: Факториал Пресс, 2002.

21. Эрроу К. Дж., Гурвиц Л., Удзава Х. Исследования по линейному и нелинейному программированию. М.: ИЛ, 1962. 336 c. [Англ. оригинал: Arrow K. J., Hurwicz L., Uzawa H. Studies in linear and nonlinear programming. Stanford University Press, 1958].

22. Экланд И., Темам Р. Выпуклый анализ и вариационные проблемы. М.: Мир, 1979.

23. Гловински Р., Лионс Ж.-Л., Тремольер Р. Численное исследование вариационных неравенств. М.: Мир, 1979.

24. Hestenes M.R. Multipliers and gradient methods // J. Optim. Theory Appl. 1969. V. 4. Р. 303-320.

25. Powell M. J.D. A method for nonlinear constraints in minimization problems // Optimization (R. Fletcher, ed.). New York: Academic Press, 1969. Р. 293-298.

26. Бертсекас Д. Условная оптимизация и методы множителей Лагранжа. М.: Радио и связь, 1987.

27. Гольштейн Е.Г., Третьяков Н.В. Модифицированные функции Лагранжа. Теория и методы оптимизации. М.: Наука, 1989.

28. Gaikovich K.P., Gaikovich P.K. Inverse problem of near-field scattering in multilayer media // Inverse Problems. 2010. V. 26. № 12. P. 125013 (17 pp.).

DUAL REGULARIZATION IN ONE-DIMENSIONAL INVERSE SCATTERING PROBLEMS K.P. Gaikovich, P.K. Gaikovich, O.E. Galkin, A.I. Smirnov, M.I. Sumin

A method is proposed to restore subsurface permittivity profiles of one-dimensional inhomogeneous media by scat-tered-field measurements based on the Lagrange formalism in the solution of nonlinear ill-posed problems. This method enables the solution of problems in their natural statements for initial Maxwell equations. The method efficiency is demonstrated in numerical simulation of the retrieval of Earth’s crust conductivity profiles and the profiles of permittivity inhomogeneities in multilayer periodic structures by spectral measurement data.

Keywords: inhomogeneous media, multilayer nanostructures, permittivity, Earth’s crust conductivity, scattered field, ill-posed inverse problems, dual regularization.

i Надоели баннеры? Вы всегда можете отключить рекламу.