Методы решения линейных некорректных задач с априорной информацией и оценка погрешностей тема автореферата и диссертации по математике, 01.01.03 ВАК РФ

Чжан Е АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2014 ГОД ЗАЩИТЫ
   
01.01.03 КОД ВАК РФ
Диссертация по математике на тему «Методы решения линейных некорректных задач с априорной информацией и оценка погрешностей»
 
Автореферат диссертации на тему "Методы решения линейных некорректных задач с априорной информацией и оценка погрешностей"

На правах рукописи

Чжан Е

Методы решения линейных некорректных задач с априорной информацией и оценка погрешностей

01.01.03 — Математическая физика АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

3 О ЯНВ ¿014

Москва — 2014

005544670

005544670

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

Научный доктор физико-математических наук,

руководитель профессор Ягола Анатолий Григорьевич

Официальные доктор физико-математических наук,

оппоненты: профессор Леонов Александр Сергеевич,

Национальный исследовательский ядерный университет «МИФИ»

доктор физико-математических наук, гл.н.с. Инна Эдуардовна Степанова, Институт физики Земли им. О.Ю. Шмидта Российской академии наук

Ведущая Российский университет

организация дружбы народов (РУДН)

Защита состоится 20 февраля 2014 г. в 16:30 часов на заседании диссертационного совета Д 501.002.10 при Московском государственном университете имени М. В. Ломоносова по адресу: 119991, г. Москва, ГСП-1, Ленинские горы, МГУ, д. 1, стр. 2, физический факультет, СФА.

С диссертацией можно ознакомиться в фундаментальной библиотеке МГУ.

Автореферат разослан « » января 2014 г.

Ученый секретарь -" Поляков П. А.

диссертационного совета, профессор

Общая характеристика работы

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

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

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

Аг = и, г е г, и в и, (1)

где Z - пространство решений, а и - пространство измерений. Физический смысл 2, А и и следующий: г - искомая физическая характеристика исследуемого объекта, А - оператор, определяющий преобразование искомого решения г в результат измерений и. Большинство обратных задач заключается в том, что даны неточные правая часть и оператор, и нужно найти приближённое решение, хорошо аппроксимирующее точное.

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

Российским математиком А. Н. Тихоновым в 60-х годах прошлого века была заложены основы теории решения некорректных задач. Им же был предложен регуляризирующий алгоритм, основанный на

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

В настоящее время двумя важными направлениями исследований в этой области являются:

1) поиск хотя бы одного регуляризирующего алгоритма, с помощью которого можно решить поставленную практическую задачу;

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

В диссертационной работе проведена работа по этим обоим направлениям.

1) Решена задача восстановления функции распределения размеров частиц аэрозоля в атмосфере.

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

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

Цель работы

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

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

Положения, выносимые на защиту

1) Метод точек перегиба для решения линейных некорректно поставленных обратных задач математической физики с априорной информацией о принадлежности решения к множеству ограниченных кусочно-выпуклых функций и способ оценки погрешности полученного приближённого решения.

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

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

4) Многошаговые алгоритмы и соответствующие программные комплексы решения обратной задачи для уравнения теплопроводности и задачи восстановления истинного изображения по дефокусиро-ванному.

Научная новизна

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

Теоретическая и практическая значимость

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

Личный вклад автора

Основные результаты, полученные в данной диссертационной работе, были впервые получены автором. Постановка математической задачи и анализ научных результатов проводились под руководством А. Г. Яго-лы и при совместном обсуждении с Д. В. Лукьяненко. Постановка задачи восстановления функции распределения размеров частиц аэрозоля в атмосфере проводилась совместно с Янфей Ваном из Института геологии и геофизики Китайской академии наук. Основное содержание диссертационной работы и её результатов полностью отражено в девяти научных публикациях автора. В материалах совместных публикаций личный вклад автора является определяющим.

Апробация работы

Основные результаты диссертационной работы были представлены: на международной конференции «The Second International Workshop on Computational Inverse Problems and Applications» (Китай, Пекин, 1215 июля 2010 года, Институт геологии и геофизики Китайской Академии Наук); на научном семинаре «Обратные задачи математической физики» под руководством А. В. Вакушинского, А. В. Тихонравова и А. Г. Яголы, проводящемся в НИВЦ МГУ (28 марта 2012 года и 11 декабря 2013 года); на международной конференции «4th International conference on «Function spaces. Differential operators. General topology. Problems of mathematical education» » (Москва, 25-29 марта 2013 года,

РУДН); на международной конференции «4th International Symposium on «Inverse problems, Design and Optimization Symposium» » (Франция, Альби, 26-28 июня 2013 года); на международной конференции «The Third International Workshop on Computational Inverse Problems and Applications» (Китай, Нанчанг, 8-12 июля 2013 года, Восточно-китайский технологический институт); на международной конференции «Методы создания, исследования и идентификации математических моделей» (Новосибирск, Академгородок, 10-13 октября 2013 года); на международной конференции «Теория и численные методы решения обратных и некорректных задач» (Новосибирск, Академгородок, 8-13 октября 2013 года); на научном семинаре кафедры математики физического факультета МГУ имени М. В. Ломоносова под руководством профессора В. Ф. Бутузова (4 декабря 2013 года).

Публикации

По теме диссертации опубликовано 9 работ, из которых 3 статьи в рецензируемых печатных научных журналах [1-3] и 6 тезисов конференций [4-9]. В журналах из списка ВАК опубликовано 3 статьи [1-3].

Структура работы

Диссертация состоит из введения, двух глав, заключения и списка литературы (74 наименования). Общий объем работы составляет 122 страниц.

Краткое содержание работы

В ПЕРВОЙ ГЛАВЕ рассматривается теория регуляризации на множестве ограниченных кусочно-выпуклых функций и построен метод точек перегиба для задач, численное решение которых сводится к необходимости решать недоопределённые системы линейных алгебраических уравнений (в диссертационной работе в качестве примера применения разработанного метода приведено решение задачи восстановления функции распределения размеров частиц аэрозоля в атмосфере) и переопределённые системы линейных алгебраических уравнений (в диссертационной работе в качестве примера приведено решение интегрального уравнения с ядром функции Грина).

Если предположить, что частицы аэрозоля в атмосфере можно описать шаровыми частицами, показатели преломления вещества которых известны, то соотношение между оптической толщиной рассеяния аэрозолем излучения с определённой длиной волны и функцией рас-

пределения размеров частиц можно описать в интегральном виде:

оо оо

Таего(А) = 11 пг2(Эехг(г, А, г?) п(г, г) йгйг. (2)

о о

Здесь п(г, г) йг - функция плотности частиц аэрозоля с радиусами частиц, лежащими в интервале между г и г+йг, на высоте г, (Зехг{г, А, т?) — эффективный фактор поглощения излучения с определенной длиной волны частицами аэрозоля по теории Ми, т) - комплексный показатель преломления, Л - длина волны излучения, таег0 - измеряемая с помощью фотометра величина оптической толщины рассеяния аэрозоля. Пусть п(г) = /0°° п(г, г) йг, уравнение (2) можно переписать в виде:

00

= / тг г2д„«(г, А, Г}) п( г) ¿г. (3)

о

В этом случае п(г) играет роль функции плотности частиц аэрозоля радиуса г в вертикальном сечении атмосферы. Она показывает, что в единичном объёме атмосферы находится п(г) ¿г частиц аэрозоля с радиусами, распределёнными в интервале между г и г + <1т. Кроме того, практический интерес представляют частицы с радиусом, заключенным в интервале [0.1,2.0] мкм.

Функцию распределения размеров частиц аэрозоля п(г) обычно определяют как произведение функций Л(г) и /(г), где Л(г) - быстро убывающая функция, /(г) - медленно меняющаяся функция. Исходя из того, что большинство измерений функции распределения размеров частиц аэрозоля над континентами показывают, что эти функции могут быть описаны юнговским распределением Л(г) = г"^'"1"1' (у* -константа формы частицы, которые обычно имеют характерный размер [0.1,4.0] мкм), имеет смысл использовать Н(г) в виде указанного юнговского распределения с весовым коэффициентом /(г).

На практике в качестве входных данных обычно используют результаты измерения оптической толщины рассеяния аэрозоля только для небольшого числа длин волн ({А;}у=1), что связано с техническими ограничениями используемых измерительных приборов.

Т.е. на практике получается система из четырёх уравнений

10\ тгг^Сг, Ау, т?) Цг) /(г) йт = тоего(Ау), ] = 1,2,3,4, (4)

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

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

В диссертации рассматриваются два случая.

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

и /"(г).

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

Для проверки разработанных методов была рассмотрена модельная задача. Дана функция распределения размеров частиц аэрозоля (Рис. 1). Положим V* = 0.15, 77 = 1.45, уровень погрешности задания оператора и правой части задачи (2) составляет 45% относительно величины точных данных. Количество узлов и уравнений составляет 100 и 4 соответственно. С помощью предложенного метода было получено значение параметра регуляризации метода точек перегиба: а = (Р,к) — (3,(25,49,74)) (где Р - число точек перегиба, к - вектор координат расположения точек перегиба). Если умножить медленно меняющуюся функцию /(г) на быстро убывающую функцию Л(г), то получится функция распределения размеров частиц аэрозоля п(г) (Рис. 1, в двойном логарифмическом масштабе).

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

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

logn(log(r))

Рис. 1. Зависимость логарифмической функции распределения размеров частиц аэрозоля к^ п(к^(г)) (единица измерения: м-2 • мкм-1) от логарифмического радиуса частиц аэрозоля к^(г) (единица измерения: мкм) при уровнях погрешностей Н = 5 = 45%. Точное решение (—). Приближённое решение по методу А.Н. Тихонова с выбором параметра регуляризации по обобщённому принципу невязки (-А-) и методу точек перегиба (— о —).

Рис. 2. Применение метода точек перегиба для решения реальной задачи. Результаты: восстановленная логарифмическая функция распределения размеров частиц аэрозоля log 7i(log(r)).

Рис. 3. Сравнение результатов, полученных с помощью различных методов при одинаковом уровне погрешностей входных данных (5 = /г = 20%). Точное решение (—); приближённое решение по методу ОПН (-о-); приближённое решение по методу МТП (— Л —).

следующей задачи: 1

I К{х, з)г(а)с1з = и(х), ф) е Ь2[О,1], и(ж) 6 Ь2[0,1], о

где ядро имеет вид

4 ' 1 й(1 — X) при Я ^ X.

Вместе с методом МТП рассмотрен также и регуляризируюгций алгоритм выбора параметра регуляризации по обобщённому принципу невязки (ОПН). На Рис. 3 представлены результаты применения обоих методов для решения поставленной задачи при одинаковых уровнях погрешности.

Во ВТОРОЙ ГЛАВЕ описываются методы построения поточечного псевдооптимального регуляризирующего алгоритма решения линейных некорректно поставленных обратных задач математической физики и оценки апостериорной погрешности решения, полученного с помощью этого алгоритма, на основе использования априорной информации о решении.

Пусть X = [ЬХ,Н.Х], У = [Ьу.Яу] и Z := С{Х), и := Ь2(У), а оператор А: Z —> I/ - линейный непрерывный инъективный оператор.

Предполагается, что априорная информация о решении Z является некоторым замкнутым ограниченным выпуклым и уравновешенным множеством в пространстве Z. Рассматривается следующее операторное уравнение:

Az-й, z€Z, й £ U. (5)

Пусть вместо точно заданных оператора А и правой части й известны лишь такие их приближения {Аьл,что Ц.А — ^ ^л и||й — гг^Цу^й. Для многих конкретных операторов А аналитическое решение z операторного уравнения (5) найти невозможно. Поэтому при решении прикладных задач математической физики соответствующее операторное уравнение (5) решается с использованием численных методов, т.е. в качестве решения задачи (5) ищется не сама функция z(x), а её конечномерная аппроксимация {z(xk)}k=i (где К - число узлов сетки, на которой ищется неизвестная функция). Далее, с помощью полученной конечномерной аппроксимации строится приближённое непрерывное решение задачи (5), и доказывается сходимость этого приближенного непрерывного решения в пространстве Z{= С(Х)).

Сначала рассматривается задача восстановления значения функции z(x) в конкретной точке Хк (узле к).

Определение 1 Методом восстановления значения функции z(x) в точке Хк (по информации Z) называется любой функционал и*: U —t R1, а погрешностью восстановления с помощью метода и* называется величина

Ao(xk,Z,£,u*) sup \z{xk) - u*(u)|, (6)

zez,

УиеСЛ||АьЛ*-и||<е

где e := 5 + Ид ■ sup2gg \\zl\z- Оптимальной погрешностью восстановления значения функции z(x) в точке Хк называется величина

Ai(®к, Z, е) := inf Д0(хк, Z, е, и'), (7)

где точная нижняя грань берётся по всем функционалам и*: U —> R1. Функционал й*, на котором эта точная нижняя грань достигается, называется оптимальным методом восстановления значения функции z(x) в точке Хк-

По теореме Рисса об общем виде линейного непрерывного функционала и теореме Смоляка предыдущую постановку (6)-(7) экстремальной задачи можно переписать в следующем виде:

A1(xk,Z,e) = inf sup \z(xk) ~ (A,u)|, (8)

где (-, •) - скалярное произведение в Ь2(У). Ассоциированная задача к задаче (8) будет ставится как

вир 2(2*), где г0 := 1п{г : \\АНлг\\ < е}. (9)

Если определить функцию Лагранжа £ : {г х £7) х £7 -> Е1 в виде £((г, и), А) —г(хк) + (А, и), то справедлива следующая теорема, которая описывает связь между исходной задачей (8) и ее ассоциированной задачей (9).

Теорема 1 (Принцип Лагранжа) Если элемент £ является допустимой тонкой в задаче (9) (т.е. г £ Еа), то

1) следующие два условия эквивалентны:

а) г является решением задачи (9);

б)3\£и*: £((г,0),А)= ш£ £((*,«),А);

2) при выполнении этих двух эквивалентных условий линейный функционал А = А является методом оптимального восстановления в задаче (8), и его погрешность равна

А!{хк,г,е) = 2{хк) = -£((г, 0),А).

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

В диссертационной работе с целью построения конкретного выпуклого и уравновешенного множества априорной информации об искомом решении рассматриваются некоторые предположения о свойствах задачи (5). В том числе, предполагается истокопредставимость решения, т.е. г = Ву, где V 6 Ьг и В — интегральный оператор.

Введём множество всей априорной информации задачи (8) П := {(г, и) 6 гПохи : \\и-АкАг\\и < е}, где ё := 5 + кА-{\\ВНв\\у^2 + кв)-па, число По находится с помощью метода расширяющихся компактов, и ВНв - приближённый оператор к В такой, что ЦБ - Вив\\v-tz ^ Лв-

Далее в диссертации рассматривается конечномерный аналог задачи (8) для всех узлов к = 1, К

Ы зир к = (10)

Кеии (г,и)еП**

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

Здесь — конечномерный аналог для множества П и zk — к-ая компонента вектора z. Связь между конечномерной задачей (10) и исходной задачей (8) также рассматривается в диссертации.

Исходя из принципа Лагранжа, вместо задачи (10) нам нужно решить неравенство

{2k-zk) + ( Af,u)^0, V(z,u)en£, k = lj<, (11)

или в векторном виде:

(z — z) + Ли У 0, V(z,u)en^, (12)

где символ z 0 означает, что все элементы вектора z неотрицательны. Л (Af.....называется матрицей множителей Лагранжа.

Множество решений неравенства (12) (или (11)) обозначим как Л. Далее в диссертации подробно рассматриваются некоторые методы нахождения одной (оптимальной в некотором смысле) матрицы множителей Лагранжа Л. Сформулируем один из них.

Среди всех решений (12) из множества Л выбираем такой элемент Л, что

л= (ATA + aiKy1AT, (13)

где Ik - единичная матрица в х иа^О - параметр регуляризации.

Для любой матрицы А верно сингулярное разложение А = EHFT, где Е, F - -ортогональные матрицы размеров М х М и К х К соответственно, Е = diagfcri,£гГ| 0,0) — диагональная матрица порядка М х К, причём cti > <72 > ... > оу > 0, число j— ранг матрицы А. Тогда с помощью разложения в ряд Тейлора можно показать, что неравенство (11) эквивалентно неравенству

0(a2) + dka + ^^dk, k = Tj<, V(Z,u)eii£, (14)

где dkt dk и dk — некоторые функции, зависящие от точки (z, и).

Очевидно, что для любого положительного числа Ео существует а0 > 0 такое, что Va G (0, а0] и Vfc: 1 ^ к ^ К верно \0(а2)\ < е0а.

В диссертации рассматривается метод выбора числа ао при заданном числе е0. Кроме того, в диссертации сформулировало достаточное условие, при котором система неравенств (14) при условии а £ (0, ао] имеет решение {а : 0 < а ^ к}.

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

(Ь)

(с) (с1)

Рис. 4. Изображения, реконструированные для уровня ошибки К = 8 = 10%. (а) неточная функция рассеяния точки К^х, у, я, ¿), (Ь) размытое изображение и(с) изображение, реконструированное с априорным выбором параметра регуляризации а = а^, (с!) изображение, реконструированное с апостериорным выбором параметра регуляризации а = а^-

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

Далее в диссертации рассматриваются два способа выбора параметра регуляризации а — априорный и апостериорный.

а) Априорный метод: а 1 где 0 < а < 2.

б) Апостериорный метод: а2 := ппп(/с, а*), где а* — параметр регуляризации, выбирающийся по обобщённому принципу невязки.

Достаточные условия существования матрицы множителей Ла-гранжа, имеющей вид (13), формулируются во второй главе диссертации. Для такой матрицы множителей Лагранжа определяется метод Яа задачи (5) по формуле: Даи := Ли. В диссертации доказывается псевдооптимальность и регуляризирующие свойства этого оператора.

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

В ЗАКЛЮЧЕНИИ диссертационной работы сформулированы основные результаты.

1) Создан метод точек перегиба решения линейных некорректно поставленных обратных задач математической физики на множестве ограниченных кусочно-выпуклых функций и вычисления оценки погрешности полученного приближённого решения.

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

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

4) Разработан алгоритм выбора параметра регуляризации для матрицы множителей Лагранжа, получены достаточные условия существования оптимального регуляризирующего алгоритма.

5) С помощью поточечного псевдооптимального регуляризирующего алгоритма решены обратная задача для уравнения теплопроводности и задача устранения размытия изображения.

Основные результаты диссертации опубликованы в следующих изданиях

Публикации в изданиях из Перечня ВАК:

[1] Ван Я., Чжан Б, Лукьяненко Д.В., Ягола А.Г. Метод решения обратной задачи восстановления функции распределения размеров частиц аэрозоля в атмосфере на множестве кусочно-выпуклых функций // Вычислительные методы и программирование,—2012,—т. 13, —с. 49-66.

[2] Wang Y. F. , Zhang Y., Lukyanenkob D. V. and Yagola A. G. Recovering aerosol particle size distribution function on the set of bounded piecewise-convex functions // Inverse Problems in Science and Engineering. —2013. —V. 21, —P. 339-354.

[3] Чжан E, Лукьяненко Д.В., Ягола А.Г. Применение принципа Лагранжа для решения линейных некорректно поставленных обратных задач с использованием априорной информации о решении ¡I Вычислительные методы и программирование. — 2013.— т. 14,—с. 468-482.

Публикации в других научных изданиях:

[4] Zhang Y. A kind of numerical methods for recovering aerosol partical size distribution function in compact set // The Second International Workshop on Computational Inverse Problems and Applications, A workshop at the Institute of Geology and Geophysics, The Chinese Academy of Sciences, Beijing, China, July 12 - July 15, 2010, —P. 49.

[5] Zhang Y. Recovering aerosol particle size distribution function on the set of bounded piecewise-convex functions // The fourth International conference on "Function spaces. Differential operators. General topology. Problems of mathematical education", Peoples' Friendship University of Russia, Moscow, Russia, March 25 - March 29, 2013. — P. 372-373.

[6] Yagola A.G., Zhang Y., Lukyanenko D.V. Regularizing algorithm for recovering solutions of ill-posed problems on the set of bounded piecewise-convex functions // The fourth International Symposium on "Inverse problem, Design and Optimization

Symposium", The Chinese Academy of Sciences, Albi, Prance, June 26 - July 28, 2013. — P. 77-79.

[7] Zhang Y. Using Lagrange Principle for solving linear inverse and ill-posed problems // The Third International Workshop on Computational Inverse Problems and Applications, East China Institute of Technology, Nanchang, China, July 8 - July 12, 2013.— P. 18-20.

[8] Yagola A.G., Zhang Y., Lukyanenko D.V. A method for solving one dimensional Fredholm integral equation of the first kind on the set of bounded piecewise-convex functions // Международная научная конференция "Методы создания, исследования и идентификации математических моделей", Novosibirsk, Akademgorodok, 10 - 13 October 2013 —P. 103.

[9] Yagola A.G., Zhang Y., Lukyanenko D.V. A method for solving one dimensional Fredholm integral equation of the first kind on the set of bounded piecewise-convex functions // The Fifth International Scientific Conference and Young Scientists School "Theory and Computational Methods for Inverse and Ill-posed Problems", Novosibirsk, Akademgorodok, 8-13 October 2013 — P. 108.

Подписано в печать: 13.12.2014 Объем: 1,0 п. л. Тираж: 100 экз. Заказ № 812 Отпечатано в типографии «Реглет» 119526, г. Москва, пр-т Вернадского, д. (495) 363-78-90; www.reglet.ru

 
Текст научной работы диссертации и автореферата по математике, кандидата физико-математических наук, Чжан Е, Москва

Московский государственный университет им. М. В. Ломоносова физический факультет, кафедра математики

На правах рукописи

04201455951

Чжан Е

Методы решения линейных некорректных задач с априорной информацией и оценка погрешностей

Специальность 01.01.03: Математическая физика

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель доктор физико-математических наук,

профессор А. Г. Ягола

Москва — 2014

Содержание

Введение............................................................................................................................................3

1 Метод решения линейных некорректно поставленных обратных задач математической физики на множестве ограниченных кусочно-выпуклых функций............................................................................20

1.1 Теория регуляризации на множестве ограниченных кусочно-выпуклых функций....................................................................................................20

1.2 Постановка физической задачи......................................................................22

1.3 Общая схема поиска приближенного решения и методы оценки погрешности..................................................................................................30

1.4 Конечномерная аппроксимация......................................................................32

1.5 Метод точек перегиба..............................................................................................37

1.5.1 Случай одной точки перегиба..........................................................39

1.5.2 Случай нескольких точек перегиба............................................48

1.5.3 Теорема о сложности метода точек перегиба....................51

1.5.4 Метод локального полного перебора........................................53

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

1.6.1 Модельная задача 1 (недоопределённая система)________56

1.6.2 Реальная задача............................................................................................63

1.6.3 Модельная задача 2 (переопределённая система)________67

2 Применение принципа Лагранжа к линейным некорректно поставленным обратным задачам..................................................................................71

2.1 Постановка задачи и метод решения..........................................................72

2.2 Выпуклость и уравновешенность множества априорной информации ..............................................................................................................................75

2.3 Новая постановка задачи и ее конечномерная аппроксимация ..............................................................................................................................................79

2.4 Принцип выбора нормального множителя Лагранжа..............86

2.4.1 Первый подход................................................................................................86

2.4.2 Второй подход................................................................................................86

2.4.3 Третий подход................................................................................................92

2.4.4 Четвертый подход........................................................................................96

2.5 Методы выбора параметра регуляризации а....................................97

2.5.1 Априорный выбор параметра регуляризации..................97

2.5.2 Апостериорный выбор параметра регуляризации________98

2.6 Регуляризирующее и оптимальное свойства алгоритма________99

2.7 Общий алгоритм решения операторного уравнения..................102

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

2.8.1 Пример 1..............................................................................................................105

2.8.2 Пример 2..............................................................................................................107

Заключение....................................................................................................................................111

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

Введение

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

В работе строится регуляризирующий алгоритм решения линейных некорректно поставленных обратных задач на множестве ограниченных кусочно-выпуклых функций — метод точек перегиба (МТП), и строится оценка погрешности приближённого решения. С помощью этого метода решаются задачи математической физики, численное решение которых сводитс!^1^еобх6лдмости решать недоопределённые и переопределённые системы линейных алгебраических уравнений.

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

Актуальность темы Многие современные задачи математической физики являются обратными задачами, которые могут быть представлены в виде операторного уравнения

Аг — и, г е г, и е и, (0.1)

где Z - пространство решений, а и - пространство измерений. Физический смысл г, А и и следующий: г - искомая физическая характеристика исследуемого объекта, А - оператор, определяющий преобразование искомого решения 2; в результат измерений и. Большинство обратных задач заключается в том, что даны неточные правая часть и оператор, и нужно найти приближённое решение, хорошо аппроксимирующее точное решение.

Большинство обратных задач, к которым сводятся прикладные задачи математической физики, являются некорректно поставленными. По определению Адамара корректной (корректно поставленной) задачей называется любая задача, у которой решение: (а) существует, (б) единственно и (в) непрерывно зависит от входных данных. А все остальные задачи Адамар называл некорректными (некорректно поставленными). Другими словами, некорректной считалась задача, у которой нарушается хотя бы одно из трёх приведённых выше свойств. Проблемы существования и единственности решения на практике обычно решаются с помощью введения понятия обобщённого (или квази-) решения. А вот проблема устойчивости решения задачи по отношению к возмущениям входных данных имеет принципиальное значение при практическом решении обратных задач на ЭВМ. Именно в связи с наличием неустойчивости решений реальных задач и появилась теория регуляризации.

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

В настоящее время двумя важными направлениями исследований в этой области являются:

1) поиск хотя бы одного регуляризирующего алгоритма, с помощью

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

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

В диссертационной работе проведена работа по этим обоим направлениям.

1) Решена задача восстановления функции распределения размеров частиц аэрозоля в атмосфере.

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

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

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

Цель работы

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

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

Положения, выносимые на защиту На защиту выносятся положения теоретического и прикладного характера.

1) Метод точек перегиба для решения линейных некорректно поставленных обратных задач математической физики с априорной информацией о принадлежности решения к множеству ограниченных

кусочно-выпуклых функций и способ оценки погрешности полученного приближённого решения.

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

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

4) Многошаговые алгоритмы и соответствующие программные комплексы решения обратной задачи для уравнения теплопроводности и задачи восстановления истинного изображения по дефокусирован-ному.

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

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

Эти методы решения конкретных задач реализованы в виде комплекса программ.

Теоретическая И практическая значимость Результаты, полученные в диссертации, могут быть применены как для решения рассмот-

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

ЛИЧНЫЙ вклад автора Основные результаты, полученные в данной диссертационной работе, были впервые получены автором. Постановка математической задачи и анализ научных результатов проводились под руководством А. Г. Яголы и при совместном обсуждении с Д. В. Лукья-ненко. Постановка задачи восстановления функции распределения размеров частиц аэрозоля в атмосфере проводилась совместно с Я. Ван из Институте геологии и геофизики Китайской академии наук. Основное содержание диссертационной работы и её результатов полностью отражено в девятых научных публикациях автора. В материалах совместных публикаций личный вклад автора является определяющим. Апробация работы Основные результаты диссертационной работы были представлены: на международной конференции «The Second International Workshop on Computational Inverse Problems and Applications» (Китай, Пекин, 12-15 июля 2010 года, Институт геологии и геофизики Китайской Академии Наук); на научном семинаре «Обратные задачи математической физики» под руководством А. Б. Бакушинского,

А. В. Тихонравова и А. Г. Яголы, проводящемся в НИВЦ МГУ (28 марта 2012 года и 11 декабря 2013 года); на международной конференции «4th International conference on « Function spaces. Differential operators. General topology. Problems of mathematical education» (Москва, 25-29 марта 2013 года, РУДН); на международной конференции «4th International Symposiurr on «Inverse problems, Design and Optimization Symposium» » (Франция, Альби, 26-28 июня 2013 года); на международной конференции «The Third International Workshop on Computational Inverse Problems and Applications» (Китай, Нанчанг, 8-12 июля 2013 года, Восточно-китайский технологический институт); на международной конференции «Методы создания, исследования и идентификации математических моделей» (Новосибирск, Академгородок, 10-13 октября 2013 года); на международной конференции «Теория и численные методы решения обратных и некорректных задач» (Новосибирск, Академгородок, 8-13 октября 2013 года) ; на научном семинаре кафедры математики физического факультета МГУ имени М. В. Ломоносова под руководством профессора В. Ф. Буту-зова (04 декабря 2013 года).

Публикации По теме диссертации опубликовано 9 работ, из которых 3 статьи в рецензируемых печатных научных журналах [1-3] и 6 тезисов конференций [4-9]. В журналах из списка ВАК опубликовано 3 статьи [1-3].

Структура работы Диссертация состоит из введения, двух глав, заключения и списка литературы (79 наименований). Общий объём работы составляет 122 страниц.

В ПЕРВОЙ ГЛАВЕ рассматривается теория регуляризации на множестве ограниченных кусочно-выпуклых функций и построен метод точек перегиба для задач, численное решение которых сводится к необходимости решать недоопределённые системы линейных алгебраических уравнений (в диссертационной работе в качестве примера приме-

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

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

ОО 00

Таего{Л) = 11 тгг2(Зех^г,\,г1)п(г,г)(1г(1г. (0.2)

о о

Здесь п(г, г) <1г - функция плотности частиц аэрозоля с радиусами частиц, лежащими в интервале между гиг + йг1, на высоте г, С^ех^т> А, т)) -эффективный фактор поглощения излучения с определенной длиной волны частицами аэрозоля по теории Ми, 77 - комплексный показатель преломления, Л - длина волны излучения, таего - измеряемая с помощью фотометра величина оптической толщины рассеяния аэрозоля. Пусть п(г) = /0°° п(г, г) (1х, уравнение (0.2) можно переписать в виде:

ОО

Таето(А) = I ТГГ2де^(г, Л, 7?) п(г) &Г. (0.2')

О

В этом случае п(г) играет роль функции плотности частиц аэрозоля радиуса г в вертикальном сечении атмосферы. Она показывает, что в единичном объёме атмосферы находится п(г) дх частиц аэрозоля с радиусами, распределёнными в интервале между г и т + йг. Кроме того, практический интерес представляют частицы с радиусом, заключенным в интервале [0.1, 2.0] мкм.

Функцию распределения размеров частиц аэрозоля п(г) обычно определяют как произведение функций Н(г) и /(г), где Н(г) - быст-

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

константа формы частицы, которые обычно имеют характерный размер [0.1,4.0] мкм), имее�