Некоторые методы решения некорректно поставленных задач математической физики
На правах рукописи
Джураев Хайрулло Шарофович НЕКОТОРЫЕ МЕТОДЫ РЕШЕНИЯ НЕКОРРЕКТНО ПОСТАВЛЕННЫХ ЗАДАЧ МАТЕМАТИЧЕСКОЙ ФИЗИКИ 01.01.02 -Дифференциальные уравнения, динамические системы и оптимальное управление
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора физико-математических наук
ДУШАНБЕ - 2012
Работа выполнена в Таджикском национальном университете Официальные Танана Виталий Павлович оппоненты: доктор физико-математических наук, профессор Южно-Уральский государственный университет заведующий кафедрой вычислительной математика Мухамадиев Эргашбой доктор физико-математических наук чл.-корр. АН РТ, профессор Вологодский государственный технический университет профессор кафедры информационных систем и технологий Сафаров Джумабой доктор физико-математических наук, доцент Курган-Тюбинский государственный университет им. Н.Хусрава заведующий кафедрой математического анализа Ведущая Московский государственный университет им. М.В. Ломоносова, организация: факультет вычислительной математики и кибирнетики
Защита состоится 28 декабря 2012 г. в 11 часов на заседании диссертационного совета ДМ 047.007.01 при Институте математики Академии наук Республики Таджикистан, 734063, г.Душанбе, ул. Айни 299/4.
С диссертацией можно ознакомиться в библиотеке Института математики Академии наук Республики Таджикистан.
Автореферат разослан " " 2012 г.
Ученый секретарь Диссертационного совета Халилов Ш.Б.
I.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность проблемы. Многие задачи прикладного характера, такие как гео- и биофизические, электродинамические, газо-, гидро- и аэродинамические, задачи физики плазмы и др. сводятся к уравнениям математической физики. Фактически само построение уравнения математической физики, адекватно описывающего те или иные физические закономерности окружающего нас мира, представляет собой решение некоторой задачи, которую естественно назвать "обратной".
Исследователь наблюдает явление и пытается построить такое уравнение, решение которого обладает наблюдаемыми свойствами. Обычно в основе получаемых уравнений лежат физические законы, которые позволяют сформулировать общий вид дифференциальных соотношений. Как правило, в них присутствует некоторое число произвольных функций, определяющих свойства физической среды. Если свойства среды известны, то уравнение математической физики в сочетании с краевыми и начальными условиями позволяет предсказать развитие физического явления в пространственно-временной области. Это классическая задача для уравнений математической физики. В теории обратных задач подобные задачи называются "прямыми". В современном естествознании очень часто возникают следующие обратные задачи: известен общий вид уравнения математической физики, но характеристические свойства среды не известны, их требуется определить по наблюдаемым решениям уравнения. Типичной является ситуация, когда непосредственные измерения внутри некоторой области невозможны по тем или иным причинам, однако возможно косвенное наблюдение и качественно и количественно измерять физические поля на границе или вне этой области. В математическом отношении такие задачи должны обеспечить корректность постановки задачи.
Понятие корректности постановки задачи математической физики было сформулировано в начале ХХ века известным французским математиком Ж.
Адамаром1. Задача математической физики называется поставленной корректно, если выполняются следующие условия: 1) решение задачи существует;
2) решение задачи Адамар Ж. Задачи Коши для линейных уравнений с частными производными гиперболического типа.-М.: Наука, 1978.–357 с. Первое издание этой книги на английском языке вышло в 1932 г.
единственно;
3) решение задачи непрерывно зависит от данных задачи. Сформулировав понятие корректности, Ж. Адамар привел пример некорректной задачи для уравнения математической физики, которая, по его мнению, не соответствовала никакой реальной физической постановке. Ж. Адамар показал это на примере задачи Коши для уравнения Лапласа, который стал классическим примером некорректной задачи.
Необходимость рассмотрения некорректных в классическом смысле (по Адамару) задач математической физики в связи с проблемами интерпретации данных геофизических наблюдений была впервые указана в 1943 г. дважды Героем Социалистического труда, академиком АН СССР А.Н. Тихоновым2. Он показал, что если сузить класс возможных решений до компакта, то из существования и единственности решения следует его устойчивость.
Пути развития теории и методы решения некорректных задач связаны с именами видных математиков А.Н. Тихонова, М.М. Лаврентьева, В.К. Иванова, а также с созданными ими научными математическими школами, во многом определившими пути развития теорий и приложений некорректных задач.
Большое число задач математической физики, не удовлетворяющих условиям корректности по Адамару, сводится к операторному уравнению первого рода, то есть к уравнению вида Az = u, () где z – искомая характеристика модели, u (следствие z) – величина, сопоставляемая с реально наблюдаемой, а A – оператор связи между этими величинами (z и u элементы некоторых пространств U, F соответственно, а A – непрерывное отображение U в F ).
Такая постановка задачи исследована многими авторами. В частности, развитию теории решения некорректных задач значительно способствовал метод регуляризации, предложенный А.Н. Тихоновым3. Метод основан на радикальной идее о стабилизации минимума уклонения значений Az от заданной правой части u при помощи некоторого вспомогательного неотрицательного сглаживающего функционала (z), определенного Тихонов А.Н. Об устойчивости обратных задач. //ДАН СССР.–1943.–Т.39,№5.
Тихонов А.Н. 1) О решении некорректно поставленных задач и методе регуляризации. //ДАН СССР/–1963.–Т.151,№3. 2) О регуляризации некорректно поставленных задач. //ДАН СССР.–1963.– Т.153,№1.
на некоторой части U0 DA, которая сама является метрическим пространством.
Требуется, чтобы множества MC {z U0 ;
(z)C} были компактны в U при любом C0. Относительно решения (*) предполагается, что оно содержится в MC при некотором значении C+. Тогда решается задача минимизации по zU параметрического функционала Тихонова M [z, u]2 (Az, u) + [(z)]2, 0.
F Решение этой задачи z при определенном выборе параметра регуляризации =() принимается за приближение к искомому решению z. Доказывается, что при элементы z z в пространстве U при равномерной по ограниченности отношения 0. В дальнейшем В.Я. Арсенин, А.Б. Бакушинский, А.В. Гончарский, А.С. Леонов, И.К. Лифанов, А.Г. Ягола и другие рассматривали разные уравнения (*), и доказывали сходимость z z в основном пространстве U0, если 0 и U0 – банахово (гильбертово и другие).
М.М. Лаврентьев4 показал, что при определенных условиях, наложенных на оператор A, можно заменить задачу (*) на близкую, но уже разрешимую при любых uF. При этом существенным моментом приближенного решения (*) являлось необходимое знание как точности задания элемента u, то есть оценка уклонения F (u, u), так и функции корректности () (модуля непрерывности оператора A1 на образе множества U ). Это позволило М.М. Лаврентьеву указать алгоритм построения таких приближений zU, для которых F (u, u)0, когда U (z, z )0 (где z точное решение (*), принадлежащее компакту M ), для достаточно широкого класса задач. В этом направлении вели исследования В.Г. Васильев, С.И. Кабанихин, В.Г.
Романов, С.И. Шишатский и многие другие представители этой научной школы.
В.К. Иванов5, используя некоторые идеи математического программирования, избавился при приближенном решении (*) от задания функции корректности ().
Вместе с тем не требовалось и знания, характеризующего точность задания правой части (*). Однако метод В.К. Иванова требовал задания компактного множества M, то Лаврентьев М.М. Об интегральных уравнения первого рода. //1) ДАН СССР.–1959.–Т.127,№1. 2) ДАН СССР.–1960.–Т.133,№2.
Иванов В.К. О линейных некорректных задачах. //ДАН СССР.–1962.–Т.145,№2.
есть множества корректности задачи (*). Полученные В.К. Ивановым приближенные решения определялись как элементы zM U, на которых U (Az, u) = min U (Av, u).
vM Такое решение называется квазирешением, а сам метод – методом квазирешения.
Метод квазирешения имеет наглядную геометрическую интерпретацию, что послужило отправным пунктом для ряда исследований. В данном направлении исследования вели В.В. Васин, Ж. Лионс, В.П. Танана и другие.
Если уравнение (*) является обыкновенным дифференциальным уравнением, то операторное уравнение содержит один “главный” член, который и определяет корректную разрешимость задачи. Если же (*) является уравнением в частных производных, то соответствующие операторные уравнения содержат, как правило, уже несколько неподчиненных друг другу неограниченных операторов.
Такая ситуация возникает, например, при изучении задачи определения гравитационного поля, задачи обтекания тел сверхзвуковым потоком, ряда гео и биофизических, гидро-, аэро- и газодинамических задач. Задачи из области естествознания, в целом, приводят к необходимости решения краевых и начальных задач для уравнений эллиптических, гиперболических и параболических типов. Для этих задач не имеет место непрерывная зависимость решения от исходных данных, то есть решения этих задач неустойчивы по отношению к возмущению исходных данных.
Актуальность проблемы обусловлена потребностями в разработке новых подходов для построения регуляризирующих алгоритмов (РА), обладающих свойством устойчивости к малым изменениям исходных данных, легко реализуемых на компьютере.
Настоящая диссертационная работа посвящена исследованию этих проблем, то есть разработке алгоритмов приближенных решений ряда некорректных задач математической физики, легко реализуемых на ЭВМ.
Отметим, что В.К.Ивановым6 было построено регуляризованное решение задачи Коши для уравнения Лапласа в полосе, используя стабилизирующий множитель Иванов В.К. Задача Коши для уравнения Лапласа в бесконечной полосе. //Дифференциальные уравнения.–1965.–Т.1,№1.
вида exp(2 s2 ). М.М. Лаврентьевым7 было предложено использование классической формулы Грина для построения регуляризованного решения задачи Коши для уравнения Лапласа. Используя идеи М.М. Лаврентьева, Ш.Ярмухамедовым8 было построено регуляризованное решение задачи Коши для уравнения Лапласа. И.О.
Исломовым9, на основе работы Ш. Ярмухамедова, было построено регуляризованное решение задачи Коши для уравнения Гельмгольца. В работе И.В. Мельниковой, А.Ю.
Фрейберга10 построено регуляризованное решение первой краевой задачи для волнового уравнения. В качестве регуляризующего алгоритма взят конечный отрезок ряда.
В.Я.Арсениным11 был предложен способ построения класса регуляризирующих операторов для интегральных уравнений первого рода типа свертки с помощью интегральных преобразований, легко реализуемых на ЭВМ.
В данной диссертационной работе, изучается задача Коши для уравнения эллиптического типа, когда решение не имеет непрерывную зависимость от исходных данных. Для обеспечения непрерывной зависимости решения от исходных данных, задачи рассматриваются на классе данных Коши, состоящих из быстро убывающих вместе со всеми производными функций при |x|. Используя метод В.Я.Арсенина, построено семейство таких приближенных решений u (x, y), что при 0 равномерно по x для каждого фиксированного y сходятся к решению u(x, y). Другой проблемой, изучаемой в данной работе, являются краевые задачи для гиперболических уравнений. Характерной особенностью таких задач является некорректность, то есть не имеет место непрерывная зависимость решения от исходных данных. На Лаврентьев М.М. О некоторых некорректных задачах математической физики.–Новосибирск: СО АН СССР, 1962.–93с.
Ярмухамедов Ш. О задаче Коши для уравнения Лапласа //Автореф. дис. доктора физ.-мат.
наук.–Новосибирск, 1983.–19 с.
Исломов И.О. О задаче Коши для уравнения Гельмгольца //Автореф. дис. канд. физ.-мат. наук.– Новосибирск, 1986.–16 с.
Мельникова И.В., Фрейберг А.Ю. О регуляризации краевой задачи для уравнения колебаний.
//Журнал вычислительной матем. и матем. физики.–1985.–Т.25,№5.
Арсенин В.Я. Об одном способе приближенных решений интегральных уравнений первого рода типа сверток. //Математическая логика, теория алгоритмов и теория множеств.–Труды МИ АН СССР: Наука.–1973.–Т.133.–С.33-51.
основе регуляризирующих методов строится класс приближенных решений краевых задач для гиперболического уравнения. Также, используя эти методы, строятся семейства приближенных решений смешанных задач для уравнения эллиптического и гиперболического типов.
М.М.Лаврентьеву12 принадлежит идея замены исходного интегрального уравнения первого рода близким ему, в некотором смысле, уравнением, для которого задача нахождения решения устойчива к малым изменениям правой части и разрешима для любой правой части. Используя идею М.М.Лаврентьева, в работе построен метод приближенных решений задачи Коши для уравнения Эйлера-Пуассона-Дарбу.
Целью работы является построение семейства (класса) регуляризирующих алгоритмов (РА), обладающих свойством устойчивости к малым изменениям исходных данных.
Основные задачи
, решаемые для достижения поставленной цели, следующие:
- построение класса приближенных решений задачи Коши для уравнения эллиптического типа, обладающих свойством устойчивости к малым изменениям исходных данных, с помощью интегральных преобразований Фурье;
- построение класса РА решения краевых задач для гиперболического уравнения, обладающих свойством устойчивости к малым изменениям исходных данных, основанных на идее регуляризации задачи суммирования рядов Фурье;
- построение приближенных решений смешанной задачи для эллиптического и гиперболического уравнения, обладающих свойством устойчивости к малым изменениям исходных данных, основанных на идее регуляризации задачи суммирования рядов Фурье;
-построение приближенных решений задачи Коши для уравнения Эйлера-Пуассона Дарбу (ЭПД), исходя из возмущенной задачи;
- построение приближенных решений сингулярного интегрального уравнения первого рода, исходя из возмущенной задачи;
- исследование условий стабилизации и согласования параметра регуляризации с Лаврентьев М.М. Об интегральных уравнения первого рода. //1) ДАН СССР.–1959.–Т.127,№1. 2) ДАН СССР.–1960.–Т.133,№2.
погрешностью;
- выбор сглаживающей функции и последовательности множителей.
Методы исследования. Основные результаты диссертации получены с применением теории дифференциальных уравнений, функционального анализа, методов вычислительной математики, теории обратных и некорректных задач. При построении регуляризованных решений задач, изучаемых в работе, используется метод регуляризации для интегральных уравнений первого рода типа свертки, разработанный В.Я. Арсениным. Также используется методика М.М. Лаврентьева, разработанная для операторного уравнения первого рода, исходящая из возмущенной задачи.
Научная новизна. Впервые показана возможность аналитического метода построения семейства РА для: задачи Коши для уравнения эллиптического типа на основе интегрального преобразования Фурье;
задачи Коши для уравнения эллиптического типа на основе метода разделенных переменных (метода Фурье);
краевой задачи для гиперболического уравнения на основе суммирования рядов Фурье;
смешанной задачи для эллиптических и гиперболических уравнения, на основе суммирования рядов Фурье;
задачи Коши для уравнения Эйлера-Пуассона-Дарбу исходя из возмущенной задачи;
сингулярного интегрального уравнения первого рода исходя из возмущенной задачи.
Кроме того, в работе впервые исследованы условия стабилизации и согласования параметра регуляризации для этих задач. Показаны способы выборов сглаживающих функций и последовательности сглаживающего множителя. Установлена зависимость параметра регуляризации от погрешности.
Научная ценность работы состоит в том, что полученные в ней результаты расширяют и углубляют представления о построении приближенных решении задач математической физики, имеющих прикладное значение.
Практическая значимость результатов исследования связана с возможностью их использования для создания принципиально новых, более эффективных технологий в различных областях: геофизике, радиофизике, биофизике, гидравлике, аэродинамике, газодинамике и др.
Задачи, поставленные в ходе диссертационного исследования, решались в рамках фундаментальных и поисковых НИР, проводимых в Таджикском национальном университете (Президентского фонда фундаментальных исследований научных проектов в 2005-2007 и 2008-2010, проект 0108ТД745).
Достоверность результатов работы определяется использованием обоснованных методов построения семейств РА, способа выборов сглаживающих функций, способа выбора последовательности сглаживающего множителя и установления зависимости параметра регуляризации от погрешности.
Личный вклад автора. Выбор темы, цели и задач для диссертационного исследования осуществлялись автором лично. В совместных публикациях [11 13, 20, 23, 24, 29, 35, 46] автору в основном принадлежат: математическая постановка задач, идеи решений, общее руководство, а соавторам–конкретные выкладки и расчеты.
В диссертацию не включены результаты работ [20, 23, 24, 29, 46].
Апробация работы. Основные положения диссертационной работы докладывались и обсуждались на международных, республиканских конференциях и семинарах, в том числе: международной конференции по “Теории приближения и вычислительной математики” (Днепропетровск, 1993);
международной конференции “Математическое моделирование и вычислительный эксперимент” (Ташкент, 2002);
3-ей международной конференции по “Математическому моделированию и вычислительному эксперименту” (Душанбе, 2002);
международной конференции по “Дифференциальным и интегральным уравнениям с сингулярными коэффициентами” (Душанбе, 2003);
международной конференции “Дифференциальные и интегральные уравнения и смежные вопросы анализа” (Душанбе, 2005);
научно-теоретической конференции “Проблемы современной физики” (Душанбе, 2006, 2008);
международной научно-практической конференции “Перспективы развития науки и образования в ХХI веке” (Душанбе, март, 2007, 2009, 2010, 2011);
научно-теоретическая конференция “Проблемы современной физики конденсированных сред” (Душанбе, май, 2007, 2010, 2011);
семинаре совещания “Наука - производству” (Душанбе, июнь, 2007);
международной научно-теоретической конференции “Математические проблемы технической гидромеханики, теории фильтрации и орошаемого земледелия” (Душанбе, май, 2008);
Республиканской научной конференции “Дифференциальные и интегральные уравнения” (Душанбе, сентябрь, 2008);
международной конференции по “Математической физике и ее приложениям” (Самара, сентябрь, 2008, 2010, 2012);
Республиканской научно-методической конференции (Душанбе, 2009);
7-ой Всероссийской научной конференции с международным участием “Математическое моделирование и краевые задачи” (Самара, 2010);
городском научном семинаре ТНУ (под руководством академика АН РТ, профессора Н. Раджабова);
в Институте математики АН РТ (под руководством член-корр. АН РТ, профессора З.Х. Рахмонова) и на семинаре кафедры вычислительной математики ЮУрГУ (под руководством профессора В.П. Танана).
Публикации. По теме диссертации опубликовано 47 научных работ, в том числе 26 работ в журналах, рекомендованных ВАК для защиты докторских диссертаций.
Структура и объем работы. Диссертация состоит из введения, четырех глав, разбитых на разделы, заключения и списка использованных источников из наименований, объемом 208 страниц текста.
Основные положения, выносимые на защиту:
а) Построен класс регуляризирующих операторов решения задач математической физики, обладающих свойством устойчивости к малым изменениям исходных данных с помощью интегральных преобразований Фурье.
б) Построено семейство регуляризирующих операторов решения задач математической физики, обладающих свойством устойчивости к малым изменениям исходных данных, основанных на идее регуляризации задачи суммирования рядов Фурье.
в) Построено приближенное решение задач математической физики, исходя из возмущенной задачи. При этом важную роль играют условия стабилизации и согласования параметра регуляризации.
II. КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
Приведем основные определения и полученные результаты в том порядке, в каком они расположены в представляемой работе.
Введение содержит обзор основных публикаций и монографий, посвященных некорректным задачам и методам их решения. В этом пункте обоснована актуальность темы диссертации, сформулированы цель и основные результаты исследований, показана научная новизна и дана общая характеристика работы.
Первая глава диссертации посвящена изложению нового аппарата построения РА для одномерной задачи Коши в верхней полуплоскости, то есть рассматривается задача Коши для уравнения u(x, y) u(x, y) + cu(x, y) = f (x, y) ( x, y0) u(x, y) + a +b (0.1.1) x y с начальными условиями u(x, y) u(x, 0) = (x), lim = (x). (0.1.2) y y 2 Здесь = + - оператор Лапласа, a, b, c - заданные числа, f (x, y), (x), (x) x2 y заданные функции.
Задача (0.1.1)-(0.1.2) является некорректно поставленной по Адамару, то есть не выполняются условия непрерывной зависимости решения от исходных данных. В подразделах 1.0.1 и 1.0.2 приведены примеры, в которых задача (0.1.1)-(0.1.2) будет некорректной, то есть решение неустойчиво относительно малых возмущений исходных данных.
Пусть K класс данных Коши, f (x, y) (при любом фиксированном y0), (x) и (x) такие, что:
f (x, y) y0), (x) (x) 1) (при любом фиксированном и бесконечно дифференцируемые функции;
2) производные f (k) (x, y) (при любом фиксированном y0), (k) (x) и (k) (x) (k = 0, 1, 2,...) стремятся к нулю при |x| быстрее любой степени.
|x| Обозначим через U совокупность решений u(x, y), удовлетворяющих условиям:
а) функции u(x, y), ux (x, y) и uxx (x, y) абсолютно интегрируемы на всей оси x при любом фиксированном y0;
б) функции uy (x, y), uyy (x, y) имеют в каждом конечном интервале 0yy интегрируемую мажоранту n(x).
В общем случае, когда f (x, y) (для любого фиксированного y0), (x) и (x) только ограничена, мы будем их рассматривать как функционал над пространством K бесконечно дифференцируемых финитных функций.
Эти предположения дают возможность построения таких алгоритмов приближенных решений задачи (0.1.1)-(0.1.2), которые дают решения, обладающие свойством устойчивости к малым изменениям исходных данных, то есть к начальным условиям при любом фиксированном y0.
Пусть U класс функций одной или двух переменных бесконечно дифференцируемых по одной переменной. Рост последовательных производных этих функций будем характеризовать некоторым числом q0.
Пусть R=(, )действительная ось, D={(x, y):xR, y0}полупространство, C(D)пространство непрерывных функций u(x, y):(x, y)D с нормой u(x, y) = sup |u(x, y)|, D CH (D)={u(x, y)C(D): u(x, y) H}, C,2 (D, q)пространство функций u(x, y)C(D), бесконечно дифференцируемых по x на R и дважды дифференцируемых по y на [0, y0 ], обращающаяся в нуль вне некоторого интервала (a, a)R, причем u x u(x, y) C,2 (D,) = sup.
x 0n;
(x,y)D Пусть в (0.1.2) функции (x), (x) и f (x, y) (для каждого y[0, y0 ]) из U (C(R)) равномерно ограничены, то есть существует такое число M, что (x) M, (x) M для всех xR и всех (x), (x) и f (x, y) из L1 (R) (L2 (R)). Тогда существует решение 2u u задачи (0.1.1)-(0.1.2) u(x, y) такое, что u(x, y),, принадлежит пространству x x C,2 (D, q);
u(x, y) удовлетворяет уравнению (0.1.1) и условиям (0.1.2). Однако можно указать начальные условия, то есть функции (x), (x) и правой части (0.1.1) f (x, y), которые не обладают свойствами бесконечной дифференцируемости, например, зададим графики исходной функций в виде ломаной линии. Она не обладает свойством бесконечной дифференцируемости. Тогда вместо нахождения u(x, y) можно ставить лишь задачу о нахождении приближенного решения. В качестве приближенного решения будем брать регуляризованное решение вида u (x, y) = R(f,,, x, y, ) зависящее от параметра, где, зависит от.
Для построения РА вводятся следующие определения.
Определение 1.1. R(f,,, x, y, ) называется приближенным решением задачи Коши для уравнения эллиптического типа в окрестности заданной точности условия Коши в L1 (R)(L2 (R)), если существует функционал s(, x, y) из C(D), удовлетворяющий условию R(f,,, x, y, ) u(x, y) s(, x, y) для любого (0, 0 ] и s(, x, y)0 при 0 для всякого xR и 0yy0.
Здесь - параметр регуляризации, u(x, y) - точное решение задачи Коши для уравнения эллиптического типа, 0 и y0 - заданные числа (0 0, y0 0).
Определение 1.2. R(f,,, x, y, ) называется вполне приближенным решением задачи Коши для уравнения эллиптического типа в -окрестности заданной точности условия Коши в L1 (R)(L2 (R)), если существует функционал c(,, x, y, ) из C(D), удовлетворяющий условию R(f,,, x, y, ) R(f,,, x, y, ) c(,, x, y) для любого (0, 0 ] и c(,, x, y)0, ()0 при 0 для всякого xR и 0yy0.
Здесь f (x, y) (для любого фиксированного y0) - -приближения правой части уравнения, а (x) и (x) есть -приближения начальных данных соответственно.
Задачу (0.1.1)-(0.1.2) разбиваем на две задачи:
1) Задача Коши для однородного уравнения v(x, y) v(x, y) v(x, y) + a +b + cv(x, y) = 0 (0.1.3) x y с заданными начальными условиями v(x, y) v(x, 0) = (x), lim = (x). (0.1.4) y y 2) Задача Коши для исходного уравнения w(x, y) w(x, y) w(x, y) + a +b + cw(x, y) = f (x, y) (0.1.5) x y с нулевыми начальными условиями w w(x, 0) = 0, lim = 0. (0.1.6) y0 y Тогда решение задачи (0.1.1)-(0.1.2) запишется в виде u(x, y) = v(x, y) + w(x, y).
В разделе 1.1 рассматривается задача (0.1.3)-(0.1.4). В зависимости от коэффициентов a, b, c уравнение (0.1.3) определяется как уравнение Лапласа или как уравнение Гельмгольца. Если 4c=a2 + b2, то уравнение (0.1.3) является уравнением Лапласа, а задача (0.1.3)-(0.1.4) задачей Коши для уравнения Лапласа. Если 4c=a2 + b2, то уравнение (0.1.3) является уравнением Гельмгольца. Точнее, задача (0.1.3)-(0.1.4) есть задача Коши для уравнения Гельмгольца.
В подразделе 1.1.1 рассматривается уравнение z(x, y) = 0 (0.1.7) в области D с начальными условиями z(x, y) z(x, 0) = p(x), = q(x), (0.1.8) y y= где -оператор Лапласа, а p(x) и q(x) -заданные функции.
Тогда z(x, y) = (s, y) exp(ixs)ds, (0.1.9) 2 P (s)Ch(sy) + Q(s) Sh(sy), а P (s) и Q(s) -преобразование Фуръе где (s, y) = s соответственно функций p(x) и q(x).
Пусть r(s, ) -некоторая заданная стабилизирующая функция, определенная для всех неотрицательных значений параметра, любых sR и обладающая следующими свойствами: 1) 0r(s, )1 для всех значений 0 и sR;
2) r(s, 0)1;
3) для всякого 0 r(s, ) четная по s и r(s, )L2 (R);
4) для каждого 0 r(s, ) при |s|;
5) при 0 r(s, )1, не убывая, причем на всяком отрезке |s|s эта сходимость равномерная;
6) для любого s=0 r(s, ) -монотонно убывающая по функция и r(s, )0 при и эта сходимость равномерная на каждом отрезке [s1, s2 ];
7) для всякого 0 r(s, )·Ch(ys) и r(s, )· Sh(ys) принадлежат L2 (R) при s любом фиксированном y (y 0).
Ch(ys) Этим условиям отвечает, например, стабилизирующий фактор r(s, )= Ch(ys)+(s), где (s) заданная неотрицательная четная функция, кусочно-непрерывная на любом конеяном отрезке оси R, причем (s)0, (s=0)0 и при достаточном больших значениях |s|: (s)const0.
r(s, ) Далее, в качестве примера функции обладающей свойствами 1)-7) можно рассматривать стабилизирующий множитель вида r(s, )= exp(|s|), r(s, )= exp(2 s2n ), r(s, )= exp((s)2n ) n1.
r(s, ) Для каждой заданной сглаживающей функции рассмотрим однопараметрическое семейство R(p, q, x, y, ) вида:
R(p, q, x, y, ) = r(s, )(s, y) exp(ixs)ds. (0.1.10) 2 Обозначим:
1 f1 (x, y) = Ch(sy) exp(ixs)ds, f1 (x, y) = r(s, )Ch(sy) exp(ixs)ds, 2 1 Sh(sy) 1 Sh(sy) f2 (x, y) = exp(ixs)ds, f2 (x, y) = r(s, ) exp(ixs)ds.
2 s 2 s Доказана следующая теорема.
Теорема 1.1. Пусть функция z(x, y) вида (0.1.9) есть точное решение задачи (0.1.7)-(0.1.8). Тогда справедлива оценка R(p, q, x, y, ) z(x, y) s(, x, y), где R(p, q, x, y, ) -регуляризирующее решение вида (0.1.10) для задачи (0.1.7)-(0.1.8), s(, x, y) = n(, y) p(x) + m(, y) q(x) L1, L |p(t)|dt, |q(t)|dt, p(x) = q(x) = l1 L n(, x) и m(, y) -модули непрерывности соответственно функций f1 (x, y) и f2 (x, y), - параметр регуляризации.
L1, L Теорема 1.2. Пусть выполняются неравенства p(x)p(x) q (x)q(x) и y - заданное число (y0). Тогда для регуляризирующего решения задачи (0.1.7) (0.1.8) имеет место оценка R(, q, x, y, ) R(p, q, x, y, ) c(,, x, y), p где R(, q, x, y, ) - функция вида (0.1.10) с заменой p(x) и q(x) на p(x) и q (x), а p c(,, x, y, ) = ·( f1 (x, y) + f2 (x, y) L1 ), L |f1 (t x, y)|dt, |f2 (t x, y)|dt.
f1 (x, y) L1 = f2 (x, y) L1 = Теорема 1.3. Параметр = () является корнем уравнения =.
f1 (x, y) + f2 (x, y) L1 L Причем, если lim () = 0, то lim R(, q, x, y, ) = z(x, y).
p 0 При p(x) p(x) L2 и q (x) q(x) L имеет место:
Теорема 1.4. Если функция r(s, ) удовлетворяет условиям 1)-7), то функция R(p, q, x, y, ) вида (0.1.10) является регуляризирующим решением для задачи (0.1.7) (0.1.8).
В подразделе 1.1.2 методика, разработанная в подразделе 1.1.1, применяется для задачи: найти u(x, y) из уравнения z + 2 z = 0, = const (0.1.11) в области D с начальными условиями z(x, y, ) z(x, 0, ) = p(x), = q(x), (0.1.12) y y= где - оператор Лапласа, а p(x) и q(x) - заданные функции.
Подраздел 1.1.3 посвящен изучению задачи (0.1.3)-(0.1.4). Сделав в (0.1.3) (0.1.4) замену v(x, y) = exp ax+by ·z(x, y), получим задачи (0.1.11)-(0.1.12). Здесь 1 ax b ax = 4 (4c a2 b2 ), а p(x)=(x) · exp, q(x)=((x) + 2 (x)) · exp.
2 Итак, решение задачи Коши для однородного уравнения эллиптического типа есть решение задачи (0.1.3)-(0.1.4), которое имеет вид:
ax + by at b v(x, y) = exp {f1 (t x, y)(t) + f2 (t x, y)((t) + (t))}dt.
exp 2 2 В качестве приближенного решения принимаем значение однопараметрического семейства решения вида ax + by R(,, x, y, ) = exp ·R(, q, x, y, ).
p R(, q, x, y, ) p Здесь -приближенное решение задачи Коши для уравнения Лапласа или Гельмгольца в зависимости от коэффициентов ax b ax a, b, c, а p(x)=(x)· exp q (x)=((x) + (x))· exp, -приближенные 2 2 значения начальных данных.
В разделе 1.2 рассматривается задача Коши для неоднородного уравнения эллиптического типа, то есть уравнения (0.1.5), в полупространстве D с начальными условиями (0.1.6). Уравнение (0.1.5) определяется свойством коэффициентов a, b, c Если 4c=a2 + b2, то уравнение (0.1.5) как уравнение Пуассона или Гельмгольца.
является уравнением Пуассона. В этом случае задача (0.1.5)-(0.1.6) является задачей Если 4c=a2 + b2, то уравнение (0.1.5) является Коши для уравнения Пуассона.
неоднородным уравнением Гельмгольца. Точнее, задача (0.1.5)-(0.1.6) есть задача Коши для неоднородного уравнения Гельмгольца.
Решение задачи (0.1.5)-(0.1.6) можно представить в виде:
y w(x, y) = (x, y, )d, где (x, y, ) -решение задачи Коши для однородного уравнения эллиптического типа с начальными данными (x, y, ) (x, y, ) = 0, = f (x, ).
y y= y= Методика, разработанная в разделе 1.1 применяется в разделе 1.2 для изучения задачи Коши неоднородного уравнения эллиптического типа, то есть задачи (0.1.5) (0.1.6). В подразделах 1.2.1 и 1.2.2 применяется методика подраздела 1.1.1. Показано, что построенное приближенное решение является регуляризирующим для задачи Коши уравнения Пуассона и неоднородного уравнения Гельмгольца соответственно.
В разделе 1.3 рассматривается задача Коши для уравнения µ u(x, y) · u(x, y) + = 0, µ 0 (0.1.13) y y в области D с начальными условиями u(x, y) u(x, 0, µ) = (x), lim y µ · = (x), (0.1.14) y y где -оператор Лапласа, а (x) и (x) -заданные функции.
Получено решение (s, y, µ)·eixs ds, u(x, y, µ) = (0.1.15) 2 где (s, y, µ) = a(s, y, µ)(s) + y 1µ ·b(s, y, µ)(s), y 2k s2k (1 + µ)(k + 1 + µ ) y 2k s2k (1 µ)(k + 1 µ ) 2 a(s, y, µ) =, b(s, y, µ) =, µ µ k!(2k + 2 µ)(1 2 ) k!(2k + 1 + µ)(1 + 2 ) k=0 k= а (s) и (s) -преобразования Фурье функций (x) и (x) соответственно.
Методика подраздела 1.1.1 применяется для изучения задачи (0.1.13)-(0.1.14).
Показано, что класс однопараметрических решений вида r(s, )(s, y, µ)·eixs ds R(,, x, y, µ, ) = 2 является регуляризирующим решением для задачи (0.1.13)-(0.1.14).
Все результаты настоящей главы, то есть методы построения устойчивого к малым изменениям исходных данных решения одномерных задач Коши, можно распространять на многомерный случай.
Вторая глава диссертации посвящена изложению краевых задач для гиперболического уравнения второго порядка, то есть решение уравнения 2 u(x, t) u(x, t) q(x)u(x, t) + f (x, t) p(x) = k(x) (0.2.1) t x x в области G = {(x, t) : 0 x h, 0 t T }, удовлетворяющее граничным условиям при x = 0 и x = h, t [0, T ], u(x, t) 1 + 2 u(x, t) = (x) при t = 0, x [0, h], (0.2.2) n (x) при t = T, x [0, h].
Здесь коэффициенты p(x), k(x) и q(x) определяются свойствами среды, где происходит колебательный процесс, свободный член f (x, t) выражает интенсивность внешнего воздействия, -граница прямоугольника G, 1 и 2 -числа, а (x) и (x) -граничные скорости смещения. Если 1 =0, то имеем первую краевую задачу, если 2 =0, то вторую, а при 1 =0 и 2 =0 -третью.
Краевые задачи (0.2.1)-(0.2.2) являются некорректно поставленными по Адамару, то есть не выполняются условия непрерывной зависимости решения от исходных данных.
Некорректность задачи (0.2.1)-(0.2.2) изучена в подразделе 2.0.3. Приведены примеры, в которых решения задачи (0.2.1)-(0.2.2) являются неустойчивыми относительно малых возмущений исходных данных.
Величины p(x), k(x), q(x), f (x, t) (для каждого t[0, T ]), (x) и (x) будем [0, h].
считать неотрицательными кусочно-непрерывными функциями на k(x)C (1) (0, h), p(x), q(x), f (x, t)(для каждого Причем p(x)0 k(x)c0, q(x)0, t[0, T ]) C(0, h). При таких условиях можно доказать, что решение задачи (0.2.1)–(0.2.2) существует и единственно. Нашей задачей является обсужденние вопросов, связанных с построением приближенных решений обладающими свойствами устойчивости к малым возмущениям исходных данных.
Примем F =(0, h), Ft =(0, T ). Будем предполагать, что решение u(x, t) задачи (0.2.1) (0.2.2) для каждого t принадлежит множеству D(A) (оператор A= x2 + x – действует на функции u(x, t)) с областью определения F. Каждый элемент множества D(A) непрерывен, имеет непрерывно дифференцируемую производную q(x)u(x, t), k(x) u x такую, что квадратично суммируема на F Ft. Кроме того, все функции множества D(A) удовлетворяют однородным граничным условиям (0.2.2). Предполагается также, что решение u(x, t) задачи (0.2.1)-(0.2.2) обладает производными p(x) 2, u t квадратично суммируемыми на F Ft. Множество D(A) есть класс возможных решений. Однако нельзя брать в качестве приближенного решения задачи (0.2.1)-(0.2.2) с приближенными исходными данными произвольный элемент u(x, t) из D(A), так как такое приближенное решение не будет устойчивым относительно малых возмущений исходных данных. Поэтому в настоящей главе построен регуляризующий алгоритм решения рассматриваемой задачи (0.2.1)-(0.2.2) в обычном для некорректных задач предположении, что решение при точно заданных исходных данных существует. В качестве приближенного решения будем брать регуляризованное решение вида u (x, t) = R(f,,, x, t, ) зависящее от параметра, где, зависит от. При этом используются условия стабилизации и согласования параметра регуляризации для задачи (0.2.1)-(0.2.2). Эти понятия можно ввести следующим образом:
Определение 2.1. R(f,,, x, t, ) называется регуляризуемым решением для краевой задачи гиперболического уравнения второго порядка в окрестности заданной точности краевых данных в L2 (0, l), если существует функционал s(, x, t) из C(G), удовлетворяющий условию R(f,,, x, t, ) u(x, t) s(, x, t) для любого, принадлежащего (0, 0 ] и s(, x, t)0 при 0 для всякого x[0, l] и 0tt1.
Определение 2.2. R(f,,, x, t, ) называется вполне регуляризуемым решением для краевой задачи гиперболического уравнения второго порядка в -окрестности заданной точности краевых данных в L2 (0, l), если существует функционал c(,, x, t) из C(G), удовлетворяющий условию R(f,,, x, t, ) R(f,,, x, t, ) c(,, x, t) для любого, принадлежащего (0, 0 ] и c(,, x, t)0, ()0 при 0 для всех x[0, l], 0tt1.
В разделе 2.1 рассматривается уравнение 2 u(x, t) u(x, t) q(x)u(x, t) = k(x) (0.2.3) t x x в области G, удовлетворяющее краевым условиям при x = 0 и x = h, t [0, T ], u(x, t) 1 + 2 u(x, t) = (x) при t = 0, x [0, h], (0.2.4) n (x) при t = T, x [0, h].
Здесь –граница прямоугольника G, 1, 2 - числа, а (x) и (x) - заданные непрерывные функции на [0, l]. Предполагается, что k(x)c0, q(x)0, k(x)C (1) (0, l), q(x)C(0, l). Если 1 =0, то имеем первую краевую задачу, если 2 =0 - вторую, а при 1 =0 и 2 =0 - третью.
Решение задачи (0.2.3)-(0.2.4) можно представить в виде u(x, t) = vn (t, n )n (x, n ), (0.2.5) n= где vn (t, n )=An Fn (t) + Bn Gn (t), причем l 2 sin(n (T t)) + n 1 cos(n (T t)) Fn (t) =, An = (x)n (x, n )dx, [2 + (n 1 )2 ] sin(n T ) l 2 sin(n t) n 1 cos(n t) Gn (t) =, Bn = (x)n (x, n )dx, [2 + (n 1 )2 ] sin(n T ) а числа n и функции n (x, n ) являются собственными значениями и собственными функциями задачи Штурма-Лиувилля:
d d d d (q(x) 2 ) = 0, k(x) + 2 = 1 + 2 = 0.
dx dx dx dx x=0 x=l Здесь An, Bn - коэффициенты Фурье функций (x) и (x) по системе {n (x, n )}, n T - иррациональные числа.
Регуляризирующий алгоритм можно представить в виде R(,, x, t, ) = r(n, )vn (t, n )n (x, n ). (0.2.6) n= Здесь r(n, ) стабилизирующий множитель, обладающий следующими свойствами: 1) r(n, ) определено для всех 0 и любого натурального n;
2) 0r(n, )1 для любых 0 и n = 1, ;
3) r(n, 0)1;
4) {r(n, )}l2, для всякого ;
5) limn r(n, )= равномерно для любого (0, 0 ], где 0 – любое фиксированное положительное число;
6) lim r(n, )= для всякого n не убывая;
7) для любого n r(n, ) – монотонно убывающая по функция и lim r(n, )=0;
8) для каждого 0 {r(n, )·Fn (t)} и {r(n, )·Gn (t)} принадлежат l2 при всяком фиксированном t [0, T ].
[ 2 +( )2 ] sin( T ) r(n, )=en r(n, )= [ 2 +(n 1 )21 sin(n T n 2, n Если положить, например, или ] )+ n r(n, )=(1 + n )s, где {n } -последовательность положительных чисел, порядок роста которых при n не ниже, чем ns, где s 1, то условия 1)-8) для таких r(n, ) выполняются.
- постоянная равномерной ограниченности {n (x, n )}, Пусть M то есть |n (x, n )|M, тогда справедлива следующая теорема:
Теорема 2.1. Пусть функция u(x, t) является точным решением задачи (0.2.3) (0.2.4), а (x) и (x) непрерывные -приближения к (x) и (x) в L2 (0, l) и t заданное положительное число из [0, T ]. Тогда для каждой последовательности {r(n, )}, удовлетворяющей условиям 1)-8), однопараметрическое семейство вида (0.2.6) является регуляризирующим алгоритмом задачи (0.2.3)-(0.2.4) и для любого t[0, T ] lim R(,, x, t, ()) = u(x, t), причем lim () = 0.
Параметр =() является корнем уравнения (, t) + l(, t) =, 3M где |r(n, )Gn (t)|2.
|r(n, )Fn (t)|, l(, t) = (, t) = n=1 n= В разделе 2.2 рассмотрено уравнение с частными производными, объяснившее многие физические явления:
2 u(x, t) 2 u(x, t) u(x, t) + a(x) + b(x)u(x, t) = c(x), (x, t)G. (0.2.7) x2 t x Оно соответствует свободным колебаниям материальных тел (скорости сигналов в линии при одновременном способе управления). Считая, что a(x)a0 0, b(x)0, c(x)c0 0 - непрерывно заданные на [0, h] функции и связаные определенными соотношениями с упругими характеристиками и плотностью неоднородной среды, будем искать решение уравнения (0.2.7) при предельных условиях при x = 0 и x = h, t [0, T ], u(x, t) 1 + 2 u(x, t) = (x) при t = 0, x [0, h], (0.2.8) n (x) при t = T, x [0, h].
Здесь –граница прямоугольника G, 1, 2 -числа, а (x) и (x) - заданные непрерывные функции на [0, h]. Если 1 =0, то имеем первую краевую задачу, если 2 =0 - вторую, а при 1 =0 и 2 =0 - третью.
Применяя методики раздела 2.1, построен класс РА решения (0.2.7)-(0.2.8) с приближенными исходными данными.
Методика раздела 2.1 применяется в разделе 2.3 для изучения уравнения с частными производными, также объяснящего многие физические явления:
2 u(x, t) 2 u(x, t) u(x, t) u(x, t) A(x) + B(x) + C(x)u(x, t) = D(t) + E(t), (x, t)G.
2 x x t t (0.2.9) Оно соответствует свободным колебаниям электромагнитного поля (скорости сигналов Считая, что A(x)a0 0, в линии при последовательном способе управления).
D(t)d0 0, B(x)0, C(x)0, E(t)0 - непрерывные заданные на [0, T ] и [0, h] функции соответственно, будем искать решение уравнения (0.2.9) при предельных условиях:
при x = 0 и x = h, t [0, T ], u(x, t) 1 + 2 u(x, t) = (x) при t = 0, x [0, h], (0.2.10) n (x) при t = T, x [0, h].
Здесь - граница прямоугольника G, 1, 2 -числа, а (x) и (x) – заданные непрерывные функции на [0, l]. Если 1 =0, то имеем первую краевую задачу, если 2 =0 - вторую, а при 1 =0 и 2 =0 - третью.
Показано, что класс однопараметрических решения вида u (x, t)=R(,, x, t, ) является регуляризирующим алгоритмом для задачи (0.2.9)-(0.2.10).
Третья глава состоит из трех разделов и в ней приведен РА смешанной задачи для уравнения Лапласа и для гиперболического уравнения.
В разделе 3.1 рассматривается следующая смешанная задача для уравнения Лапласа u(x, y) = 0, G = {(x, y) : 0 x l, y 0}, (0.3.1) u |y=0 = (x), x [0, l], u(x, 0) = (x), (0.3.2) y 0 при x = 0, u(x, y) y 0.
1 u(x, y) + 2 = (0.3.3) x 0 при x = l, Здесь постоянные 1, 2 таковы, что 1 2 + 2 2 =0, а (x) и (x) непрерывные на [0, l] функции (данные задачи).
Проведя рассуждения, аналогичные разделу 2.1, получим решение задачи (0.3.1)– (0.3.3) в виде Sh(k y) u(x, y) = k Ch(k y) + k vk (x, k ), (0.3.4) k k= где k = (, vk )L2, k = (, vk )L2 - коэффициенты Фурье функций (x), (x) по системе {vk (x, k )}.
В качестве приближенного решения (0.3.1)-(0.3.3) с приближенными исходными данными будем брать значение однопараметрического семейства операторов Sh(k y) vk (x, k ), R(,, x, y, ) = r(k, ) k Ch(k y) + k (0.3.5) k k= где r(k, ) - стабилизирующие множители, определенные для всех 0 и k = 1,.
Кроме этого, предположим, что они удовлетворяют еще условиям 1)-7) раздела 2.1 и свойству:
8) для всякого 0 r(k, )·Ch(k y) и r(k, ) Sh(k y) принадлежат l2 при любом k фиксированном y0.
Теорема 3.1. Пусть функция u(x, y) является точным решением задачи (0.3.1)-(0.3.3) вида (0.3.4), а (x), (x) непрерывные -приближения к (x), (x) в пространстве L2 (0, l) и y - заданное положительное число. Тогда для каждой последовательности {r(k, )}, удовлетворяющей вышеприведенным условиям 1)-8), однопараметрическое семейство вида (0.3.5) является регуляризирующим решением задачи (0.3.1)-(0.3.3) и для любого y0 выполняется равенство lim R(,, x, y, ) = u(x, y), причем lim () = 0.
Раздел 3.2 посвящен изучению смешанной задачи для уравнения гиперболического типа второго порядка, то есть 2u 2u u u A(t) + C(x) 2 + D(t) + E(x) + [F1 (t) + F2 (x)]u = 0 (0.3.6) t x t x при x [0, l], t [0, t0 ], в предположении достаточной гладкости коэффициентов (A(t)a0 0, C(x)c0 0), D(t)0, E(x)0, F1 (t)0, F2 (x)0. Рассмотрим задачу нахождения его решения, удовлетворяющего краевым условиям u(x, 0) = (x), u(x, t0 ) = (x), 0xl;
(0.3.7) u(0, t) u(l, t) 0 u(0, t) + 0 = 0, 1 u(l, t) + 1 = 0, 0tt0, (0.3.8) x x где постоянные i, i таковы, что i 2 + i 2 =0, i = 0;
1 а (x) и (x) - непрерывные на [0, l] функции (данные задачи).
Методика, разработанная в разделах 2.1 и 3.1 применяется для изучения задачи (0.3.6)-(0.3.8). Показано, что построенное приближенное решение является регуляризирующим для задачи (0.3.6)-(0.3.8).
В разделе 3.3 методика, разработанная в разделе 3.1, применяется к задаче Коши для уравнения Лапласа в области G = {(r, ) : r0 rR, 0}, то есть для неустойчивых задач о продолжении потенциального поля.
В четвертой главе построено приближенное решение задачи Коши для уравнения эллиптического типа исходя из возмущенной задачи. А также этим методом исследовано сингулярное интегральное уравнение первого рода. При этом важную роль играют условия стабилизации и согласования параметра регуляризации.
В раздела 4.1 рассмотрено уравнение (0.1.13) в полупространстве D с начальными условиями (0.1.14).
Заметим, что решение задачи (0.1.13)-(0.1.14) можно представить в виде (0.1.15). В практических задачах начальные значения получаются в результате измерений, то есть u u(x, 0) = (x), lim y µ = (x).
y y Здесь (x) и (x) - суть -приближения (x) и (x) соответственно, которые понимаются в следующем смысле:
(x) (x), (x) (x).
(0.4.1) Введем вспомогательные задачи µ u (x, y) · u (x, y) + = 0, µ 0, (0.4.2) y+ y u u (x, 0) = (x), lim (y + )µ = (x), (0.4.3) y y µ u (x, y) · (x, y) + u = 0, µ 0, (0.4.4) y+ y u u (x, 0) = (x), lim (y + )µ = (x), (0.4.5) y y где (0, 0 ] (0 0).
Введем понятие условия стабилизации и согласования параметра регуляризации для задачи (0.1.13)-(0.1.14), то есть замены задачи (0.1.13)-(0.1.14) возмущенной (0.4.2) (0.4.5) следующим образом:
Определение 4.1. Задача (0.4.2)-(0.4.3) называется стабилизуемой для задачи (0.1.13)-(0.1.14), а соответственно u (x, y) - регуляризуемым решением, если существует функционал s(, x, y), удовлетворяющий условию:
u (x, y) u(x, y) s(, x, y) для любого (0, 0 ] и s(, x, y)0, при 0 для каждого xR и 0yy0.
Определение 4.2. Задача (0.4.4)-(0.4.5) называется вполне стабилизуемой для задачи (0.1.13)-(0.1.14), а соответственно u (x, y) - вполне регуляризуемым решением, если существует функционал (,, x, y), удовлетворяющий условию u (x, y) u (x, y) c(,, x, y) для каждого (0, 0 ] и c(,, x, y)0, ()0 при 0 для любого xR и 0yy0.
Проведя рассуждения, аналогичные рассуждениям в разделе 1.3, получим решение задачи (0.4.2)-(0.4.3) в виде f1 (x t, y, µ)(t)dt + f2 (x t, y, µ)(t)dt, u (x, y) = (0.4.6) где (y + )1µ a(s, y+, µ)eixs ds, f2 (x, y, µ) = b(s, y +, µ)eixs ds.
f1 (x, y, µ) = 2 Функцию u (x, y), определяемую (0.4.6), можно рассматривать как приближенное решение задачи (0.1.13)-(0.1.14).
Теорема 4.1. Пусть u(x, y) есть точное решение задачи (0.1.13)-(0.1.14). Тогда справедлива оценка u (x, y) u(x, y) s(, x, y), (0.4.7) где u (x, y) - решение задачи (0.4.2)-(0.4.3), s(, x, y) = (, y) (x) + l(, y) (x) L1, L |(t)|dt, (x) |(t)|dt, (x) = = L1 L (, y), l(, y) - модули непрерывности функций f1 (x, y, µ) и f2 (x, y, µ) соответственно, - параметр регуляризации.
Следствие. Для того, чтобы решение задачи (0.4.2)-(0.4.3) u (x, y)u(x, y) при 0, где u(x, y) - решение задачи (0.1.13)-(0.1.14), необходимо и достаточно, чтобы s(, x, y)0 при 0.
Теорема 4.2. Пусть u (x, y) -решение задачи (0.4.2)-(0.4.3). Тогда имеет место оценка u (x, y) u (x, y) c(,, x, y), (0.4.8) где c(,, x, y) = |f1 (x, y, µ)|L1 + |f2 (x, y, µ)|L1, |f1 (x, y, µ)|L1 = |f1 (x t, y, µ)|dt, |f2 (x, y, µ)|L1 = |f2 (x t, y, µ)|dt.
Теорема 4.3. Параметр = () является корнем уравнения =, f1 (x, y, µ) + f2 (x, y, µ) L1 L причем, если lim () = 0, то lim u (x, y) = u(x, y).
0 В разделе 4.2 рассмотривается сингулярное интегральное уравнение первого рода l u(x) dx = v(t), l t l, (0.4.9) xt l где функция v(t) удовлетворяет условию Гельдера с показателем (01;
v(t)H ) в области интегрирования.
Введем вспомогательные сингулярные интегральные уравнения l (t l)(l + t ) u (x) u (t)ln + dx = v(t), (0.4.10) (l + t)(t l ) xt+ l l (t l)(l + t ) u (x) u (t)ln + dx = v (t), (0.4.11) (l + t)(t l ) xt+ l где (0, 0 ]) (0 0).
Понятие условия стабилизации и согласования параметра регуляризации для уравнения (0.4.9), то есть понятие замены уравнения (0.4.9) возмущенным уравнением (0.4.10) ((0.4.11)), введем следующим образом:
Определение 4.3. Уравнение (0.4.10) называется стабилизуемым для уравнения (0.4.9), а соответственно u (t) - регуляризуемым решением, если существует такой функционал s(, t), что u (t) u(t) s(, t) для каждого (0, 0 ] и s(, t)0 при 0, для всех t[l, l].
Определение 4.4. Уравнение (0.4.11) называется вполне стабилизуемым для уравнения (0.4.9), а соответственно u (t) - вполне регуляризуемым решением, если существует такой функционал c(,, t), что u (t) u (t) c(,, t) для любого (0, 0 ] и c(,, t)0, ()0 при 0 для каждого t[l, l].
Теорема 4.4. Пусть u(t) - некоторое решение уравнения (0.4.9) из пространства С[-l,l]. Тогда справедлива оценка u (t) u(t) s(, t), (0.4.12) где u (t) - решение уравнения (0.4.10), ln tl tl ln l+t l+t s(, t) = ·M · или s(, t) = ·M ·, ln tl ln (tl)(l+t) + ln tl l+t (l+t)(tl) l+t а M 0 - постоянная Гельдера, - параметр регуляризации.
Теорема 4.5. Пусть u (t) - решение уравнения (0.4.11) при заданной максимальной абсолютной погрешности для v(t). Тогда имеет место оценка u (t) u (t) c(,, t), где c(,, t) = или c(,, t) =.
ln tl (tl)(l+t) tl ln + ln l+t l+t (l+t)(tl) Теорема 4.6. Параметр = () является корнем уравнения tl ·M · ln =, l+t причем, если lim () = 0, то lim u (t) = u(t).
0 В заключение сформулируем основные результаты, полученные в работе.
III. ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 1. С помощью интегральных преобразований Фурье построены классы приближенных решений задач Коши для уравнений эллиптического типа в полупространстве, обладающих свойством устойчивости к малым изменениям исходных данных.
2. Построен класс приближенных решений краевых задач для гиперболического уравнения в прямоугольной области, то есть РА, обладающих свойством устойчивости к малым изменениям исходных данных, основанных на идее регуляризации задачи суммирования рядов Фурье.
3. Построено семейство приближенных решений смешанной задачи для уравнения Лапласа, а также для уравнения гиперболического типа второго порядка в прямоугольной области, обладающих свойством устойчивости к малым изменениям исходных данных, основанных на идее регуляризации задачи суммирования рядов Фурье.
4. Построено приближенное решение задачи продолжения потенциального поля, то есть задачи Коши для уравнения Лапласа в ограниченной области.
5. Построено приближенное решение задачи Коши для уравнения Эйлера-Пуассона Дарбу исходя из возмущенной задачи в полупространстве.
6. Построено приближенное решение сингулярного интегрального уравнения исходя из возмущенного уравнения.
7. Представлен способ выбора сглаживающей функции и последовательности множителей, метод построения РА и стабилизации задач, метод определения и согласования параметра регуляризации с погрешностью.
Полученные результаты являются основой для построения приближенного решения задач, имеющих важное прикладное значение, реализуемого на компьютере.
IV. ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ Содержание диссертации опубликовано в работах:
ПУБЛИКАЦИИ В ВЕДУЩИХ РЕЦЕНЗИРУЕМЫХ ИЗДАНИЯХ, РЕКОМЕНДОВАННЫХ ДЕЙСТВУЮЩИМ ПЕРЕЧНЕМ ВАК РФ 1. Джураев Х.Ш. О решении задачи Коши для уравнения Лапласа /Х.Ш.Джураев //ДАН Тадж. ССР.–1986.–Т.29, №9.–С.506-509.
2. Джураев Х.Ш. О решении задачи Коши для уравнения Эйлера - Пуассона - Дарбу /Х.Ш.Джураев //ДАН Тадж. ССР.–1988.–Т.31,№7.–С.432-436.
3. Джураев Х.Ш. О решении задачи Коши для уравнения Эйлера-Пуассона-Дарбу (ЭПД) в многомерном случае /Х.Ш.Джураев //ДАН Тадж. ССР.–1989.–Т.32,№5.– С.292-295.
4. Джураев Х.Ш. О решении задачи Коши для уравнения Пуассона в многомерном случае /Х.Ш.Джураев //ДАН Тадж. ССР.–1989.–Т.32,№6.–С.292-295.
Джураев Х.Ш. О решении задачи Коши для системы уравнений Эйлера 5.
Пуассона-Дарбу (ЭПД) /Х.Ш.Джураев //Изв. АН Тадж. ССР. Отд. физ.-мат. хим.-геолог. наук. –1990.–№1.–С.9-12.
Джураев Х.Ш. Об одном регуляризирующем алгоритме построения 6.
приближенных решений задачи Коши для уравнения Гельмгольца /Х.Ш.Джураев //ДАН Тадж. ССР.–1991.–Т.34,№2.–С.292-295.
7. Джураев Х.Ш. Об одном регуляризирующем алгоритме решений задачи Коши для уравнения Эйлера-Пуассона-Дарбу (ЭПД) /Х.Ш.Джураев //Изв. АН Тадж.
ССР. Отд. физ.-мат.-хим.-геолог. наук.–1991.–№4.–С.16-20.
8. Джураев Х.Ш. О регуляризации решения задачи Коши для уравнения Эйлера Пуассона -Дарбу (ЭПД) /Х.Ш.Джураев //ДАН РТ.–1993.–Т.36,№6-7.–С.65-70.
9. Джураев Х.Ш. Об одном решения задачи Коши для уравнения Эйлера-Пуассона -Дарбу /Х.Ш.Джураев //ДАН РТ. –2002.–Т.45,№5-6.–С.55-60.
10. Джураев Х.Ш. Об одном регуляризуещем алгоритме решения сингулярного интегрального уравнения 1-го рода /Х.Ш.Джураев //ДАН РТ.–2003.–Т.46,№3-4.– С.53-58.
Джураев Х.Ш. О регуляризации задачи Дирихле для гиперболического 11.
уравнения /Б.Алиев, Х.Ш.Джураев //ДАН РТ.–2004.–Т.47,№3.–С.51-55.
Джураев Х.Ш. О регуляризации граничных задач для гиперболического 12.
уравнения /Б.Алиев, Х.Ш.Джураев //Вестник национального университета, ТГНУ.–2005.–№2.–С.13-16.
13. Джураев Х.Ш. О регуляризации решения краевой задачи для гиперболического уравнения /Б.Алиев, Х.Ш.Джураев //ДАН РТ.–2006.–Т.49,№4.–С.301-305.
14. Джураев Х.Ш. Об одном регуляризирующем алгоритме неустойчивой задачи о продолжении потенциального поля /Х.Ш.Джураев //Вестник национального университета, ТГНУ.–2006.–№5.–С.40-43.
Джураев Х.Ш. Об одном подходе к проблеме регуляризации задачи 15.
Коши для уравнения Эйлера-Пуассона-Дарбу /Х.Ш.Джураев //Дифференциальные уравнения.–2007.–Т.43,№5.–С.701-705.
Джураев Х.Ш. Об одном подходе к проблеме регуляризации краевых 16.
задач для гиперболического уравнения /Х.Ш.Джураев //Вестник Самарского госуниверситета, Естественнонаучная серия.–2008.–№8/2(67).–С.
17. Джураев Х.Ш. О регуляризации скорости сигналов в линии при одновременном управлении /Х.Ш.Джураев //ДАН РТ.–2009.–Т.52,№1.–С.23-29.
Джураев Х.Ш. Об одном подходе к проблеме регуляризации задачи Коши 18.
для уравнения эллиптического типа с постоянными коэффициентами/Х.Ш.Джураев //Вестник Таджикского национального университета, серия естественных наук.–2010.–3(59).–С.16-27.
19. Джураев Х.Ш. О регуляризации задачи распространения волн в анизотропной неоднородной среде /Х.Ш.Джураев //ДАН РТ.–2010.–Т.53,№2.–С.104-109.
Джураев Х.Ш. Об одном устойчивом методе определения приближенного 20.
решения уравнения тепло- и массопереноса /Х.Ш.Джураев, З.С.Норматов //Вестник Таджикского национального университета, спецвыпуск, посвящен. Году образования и технических знаний.–2010.–С.50-57.
21. Джураев Х.Ш. О регуляризации стационарного поля температуры в среде /Х.Ш.Джураев //Известия АН РТ. Отделение физико-математических, химических, геологических и технических наук.–2010.–№2(139).–С.27-33.
22. Джураев Х.Ш. О решениях смешанной задачи для уравнения Лапласа /Х.Ш.Джураев //ДАН РТ.–2011.–Т.54,№4.–С.27-33.
23. Джураев Х.Ш. Искусственная гиперболизация уравнения теплопроводности Таджикского /Х.Ш.Джураев, З.С.Норматов, Г.К.Собирова //Вестник национального университета.–2011.–4(68).–С.3-7.
Джураев Х.Ш. Метод искуственной гиперболизации для обратной 24.
задачи теплопроводности /Х.Ш.Джураев, З.С.Норматов //Вестник Таджикского технического университета.–2011.–№4.–С.3-7.
25. Джураев Х.Ш.Регуляризация краевых задач для гиперболического уравнения /Х.Ш.Джураев //Научно - технический вестник Поволжья.–2011.–№6.–C.33-36.
Джураев Х.Ш.О решениях краевых задач для волнового уравнения 26.
/Х.Ш.Джураев //В мире научных открытий. Математика. Механика.
Информатика.–2012.–№1.1(25). –C.129-142.
ПРОЧИЕ ПУБЛИКАЦИИ, СОДЕРЖАЩИЕ ОСНОВНЫЕ НАУЧНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ Джураев Х.Ш. О приближенных решениях краевых задач для уравнения 27.
диффузии /Х.Ш.Джураев //Теории приближения и вычислительной математики:
тезисы докл. междун. конф. –Днепропетровск, 26-28 мая 1993 г.–С.74.
28. Джураев Х.Ш. Об одном регуляризирующем алгоритме решения интегрального уравнения Абеля /Х.Ш.Джураев //Дифференциальные и интегральные уравнения и их приложения: сб.науч.тр., выпуск 3.–Душанбе: Изд-во ДГПУ им. К. Джураева, 1995.– С.15-17.
29. Джураев Х.Ш. Модель и решения обратной задачи определения мощности утеплителя из условия обеспечения требуемых температурных режимов в пахотном слое почвы /Х.Ш.Джураев, Б.Хамидов //Математическое моделирование и вычислительный эксперимент: тезисы докл. междун. конф. –Ташкент, 25-27 марта 2002 г.–С. 30. Джураев Х.Ш. Об одном регуляризирующем алгоритме построения решений задачи Коши для эллиптического уравнения /Х.Ш.Джураев //Математическое моделированию и вычислительный эксперимент: материалы 3-ей междун. конф. – Душанбе, 10-12 декабря 2002 г.–С.35-42.
31. Джураев Х.Ш. О построении регуляризующего алгоритма решения общего линейного гиперболического уравнения второго порядка с регулярными и сингулярными коэффициентами /Х.Ш.Джураев //Дифференциальные и интегральные уравнения с сингулярными коэффициентами: труды междун. конф. –Душанбе, 25-28 октября г.–С.66-70.
Джураев Х.Ш. О решениях задачи Коши для сингулярно-возмущенных 32.
обыкновенных дифференциальных уравнений второго порядка с постоянными коэффициентами /Х.Ш.Джураев //Перспективы развития науки и образования в ХХI веке. Ч.2: материалы II междун. научно-практической конф. –Душанбе: Изд-во ТТУ им. М.С. Осими, 2006.–С.21-23.
33. Джураев Х.Ш. Об одном регуляризующем алгоритме неустойчивой задачи о продолжении потенциального поля /Х.Ш.Джураев //Проблемы современной физики:
прог. и тезисы докл. научно-теоретической конф. посвященной 65-летию д.ф.м.н., профессора Д.С.Саидова. –Душанбе, 9 декабря 2006 г.–С.40-41.
34. Джураев Х.Ш. Регуляризация задачи Коши для уравнения Эйлера-Пуассона Дарбу в многомерном случае /Х.Ш.Джураев //Современные проблемы физики конденсированных сред: прог. и тезисы докл. научно-теоретической конф.
посвященной памяти заслуженного деятеля науки Таджикистана, профессора Б.Н.
Нарзуллоева. –Душанбе, 4-5 мая 2007 г. –С.51-56.
Джураев Х.Ш. Построения типа метода регуляризации для некоторых 35.
теоретических и прикладных задач /Б.Алиев, Х.Ш.Джураев //Наука – производству:
материалы семинара-совещания. –Душанбе: Изд-во ТГНУ, -2007.–С.34-39.
36. Джураев Х.Ш. О регуляризации решения краевой задачи для гиперболического уравнения с сингулярной точкой /Х.Ш.Джураев //Математические проблемы технической гидромеханики, теории фильтрации и орошаемого земледелия: материалы междун. научно-практической конф. посвященной 70-летию д.ф.м.н., профессора М.А.
Саттарова. –Душанбе, 27-28 мая 2008 г.–С.26.
37. Джураев Х.Ш. Об одном подходе к проблеме регуляризации краевой задачи для гиперболического уравнения /Х.Ш.Джураев //Математическая физика и ее приложения: тезисы докл. междун. конф.–Самара, 8-13 сентября 2008 г.–С.65.
Джураев Х.Ш. Регуляризации задачи Коши для эллиптического уравнения 38.
на основе теории возмущений /Х.Ш.Джураев //Дифференциальные и интегральные уравнения: материалы респуб. научной конф. посвященной 60-летию образования ТГНУ и 70-летию академика АН РТ Н.Раджабова. –Душанбе, 25-26 сентября г.–С.28-30.
39. Джураев Х.Ш. О решении некоторых некорректных задач математической физики /Х.Ш.Джураев //Наука и современное образование: проблемы и переспективы:
материалы междун. научно-практической конф. посвященной 60-летию ТГНУ. – Душанбе, 24-25 октября 2008 г.–С.106-108.
Джураев Х.Ш. О регуляризации решения граничных задач для уравнения 40.
колебания сигналов в линии /Х.Ш.Джураев //Проблемы современной физики конденсированных сред: прог. и тезисы докл. научно-теоретической конф.
посвященной 80-летию академика АН РТ А.А. Адхамова. –Душанбе, 15-ноября г.–С.59-62.
Джураев Х.Ш. О регуляризации решения граничных задач для уравнения 41.
колебания /Х.Ш.Джураев //Современные проблемы физики конденсированных сред и астрофизики: материалы междун. конф. посвященной 70-летию профессора Н.С.
Султонова. –Душанбе, 21-22 мая 2010 г.–С.76-77.
Джураев Х.Ш. Об одном подходе к проблеме регуляризации задачи Коши 42.
для уравнения Лапласа /Х.Ш.Джураев //Перспективы развития науки и образования:
материалы IV междун. научно-практической конф. –Душанбе, 20-22 мая 2010 г. – С.182-183.
43. Джураев Х.Ш. О регуляризирующем алгоритме задачи распространения волн в анизотропной неоднородной среде /Х.Ш.Джураев //Математическое моделирование и краевые задачи: труды седьмой Всероссийской научной конф. с международными участием. Ч.3. –Самара: Изд-во СамГТУ, 3-6 июня 2010 г. –С.88-92.
44. Джураев Х.Ш. О регуляризирующем алгоритме решений задачи Коши для уравнения эллиптического типа с постоянными коэффициентами /Х.Ш.Джураев //Современные проблемы математического анализа и их приложений: материалы междун. научной конф. посвященной 60-летию академика К.Х. Бойматова. –Душанбе, 23-24 июня 2010 г.–С.33-34.
Джураев Х.Ш. О регуляризирующем алгоритме решений многомерных 45.
задач Коши для эллиптических дифференциальных уравнений /Х.Ш.Джураев //Математическая физика и ее приложения: материалы второй междун. конф. – Самара, 29 августа - 4 сентября 2010 г.–С.109-110.
46. Джураев Х.Ш. Искусственная гиперболизация уравнения теплопроводности /Х.Ш.Джураев, З.С.Норматов //Современный научный вестник.–2011.–№15 (111).– С.90-94. –№16 (112).–С.92-98.
47. Джураев Х.Ш. Задачи Коши для аналога уравнения Лапласа в бесконечной полосе /Х.Ш.Джураев //Перспективы применения инновационных технологий и усовершенствования технического образования в высших учебных заведениях стран СНГ. Ч.2. : материалы V междун. научно-практической конф. –Душанбе, 13- октября 2011 г.–С.8-9.