Сетевое издание
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,940

ЧИСЛЕННЫЙ МЕТОД РЕШЕНИЯ КРАЕВОЙ ЗАДАЧИ ВТОРОГО РОДА ДЛЯ НАГРУЖЕННОГО ОПЕРАТОРА ШТУРМА-ЛИУВИЛЛЯ

Абрегов М.Х. 1 Бечелова А.Р. 1 Нахушева Ф.М. 1
1 ФГБОУ ВПО «Кабардино-Балкарский государственный университет им. Х.М. Бербекова»
Работа посвящена численному методу решения краевой задачи второго рода для нагруженного обыкновенного дифференциального уравнения. В работе также получены необходимые и достаточные условия существования и единственности решения рассматриваемой задачи. Нагруженные дифференциальные уравнения возникают при изучении движения почвенной влаги в задачах управления качеством водных ресурсов, когда в водоем поступает из точечных источников загрязняющее вещество определенной интенсивности. В классе достаточно гладких коэффициентов доказана сходимость решения разностной задачи к решению дифференциальной задачи в равномерной метрике со вторым порядком точности по шагу сетки. Основным методом исследования задачи является принцип максимума. С помощью принципа максимума получены априорные оценки в равномерной метрике, откуда следует устойчивость и сходимость схемы.
нагруженное линейное дифференциальное уравнение; однозначная разрешимость; численный метод решения.
1. Абрегов М.Х., Бечелова А.Р. Вторая краевая задача для нагруженного линейного дифференциального уравнения второго порядка // Известия КБНЦ РАН. – 2010. – № 3 (35). – C. 120-126.
2. Алиханов А.А., Березгов А.М., Шхануков-Лафишев М.Х. Краевые задачи для некоторых классов нагруженных дифференциальных уравнений и разностные методы их численной реализации // ЖВМ и МФ РАН. – 2008. – Т. 48, № 9. – С. 1-10.
3. Ильин В.А., Моисеев Е.И. Нелокальная краевая задача для оператора Штурма-Лиувилля в дифференциальной и в разностной трактовке // Доклады АН СССР. – 1986. – Т. 291, № 3. – C. 534–539.
4. Ильин В.А., Позняк Э.Г. Основы математического анализа. – М. : Наука. – Ч. 1. – 1982. – 684 с.
5. Самарский А.А. Теория разностных схем. – М. : Наука, 1983. – 616 с.
6. Тихонов А.Н., Васильева А.Б., Свешников А.Г. Дифференциальные уравнения. – М. : Наука, 1980. – 231 с.
Математические модели, возникающие при изучении ряда прикладных задач, приводят к необходимости решения краевых задач для нагруженного дифференциального уравнения второго рода. Такие примеры  можно найти в математической физике, математической биологии и других областях.

В настоящей работе  будем изучать численный метод решения задачи

                                      (1)

                              (2)

где    - оператор Штурма-Лиувилля, - фиксированная точка интервала (0,1).

            В начале установим  условия однозначной разрешимости задачи (1), (2).

Теорема 1. Пусть        и функция   такова, что для всех

            .                                                          (3)

Тогда решение задачи (1),(2) существует единственно и принадлежит классу .

Пусть - решение задачи

                                                                          (4) 

а  - решение задачи

         .                                          (5)

            В работе  методом функции Грина  показано, что необходимым и достаточным условием однозначной разрешимости  задачи  (1), (2) является условие  

                                                                                                                         (6)

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

                                                            (7)

             Покажем, что условие (3) обеспечивает выполнение (6), т.е. является достаточным условием однозначной разрешимости задачи (1), (2).  В силу непрерывности  и , (3)  означает, что либо      , либо       для всех .  Введём обозначение

                                                         (8)

и получим нижнюю оценку выражения .

Рассмотрим два случая.

1.   Пусть . Получим верхнюю оценку наибольшего значения  функции  на отрезке. Пусть  - точка максимума . Из равенства

                                   (9)

в силу      следует:

.                                                          (10)

            Если наибольшее значение  достигается на одном из концов отрезка, то переходя к пределу  при      в (9), получим оценку (10). Тогда

                   (11)

2. Пусть .  В этом случае получаем нижнюю оценку наименьшего значения  функции  на отрезке :

откуда находим

                                         (12)

Далее будем считать,  что выполнены условия А:

     .         

           

            Теорема 2. Если      удовлетворяют условию А  и выполнено (3),  то решение задачи   (1), (2) принадлежит классу

            Для численного решения задачи (1), (2) на отрезке введём равномерную сетку  .  Шаг сетки   выберем  меньше половины длины меньшего из отрезков . Номер    выберем из  условия  .

            Проекции решений задач (1)-(2), (4), (5) на сетку обозначим  соответственно.

Пусть сеточная функция   - решение конечно-разностной задачи

            ,                                                                   (13)

            

а сеточная функция  - решение конечно-разностной задачи

  ,                                                                    (14)

       

где 

      (15)

Введём обозначения:

                                          (16)

                                             (17)

 и в качестве приближённого решения задачи (1), (2) на сетке  выберем сеточную функцию :

.                               (18)

 

            Теорема 3. Пусть  выполнены условия А и (3).  Тогда сеточная функция ,  определённая по формуле (18), сходится при    к решению  задачи (1), (2) со вторым порядком точности по шагу  в равномерной метрике.

            Доказательство теоремы 3  будет заключаться в получении априорной оценки

,                                                         (19)

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

=

                                                      .                                     (20)

            Оценим слагаемые в правой части (20). Как известно , конечно-разностные схемы (13) и (14) сходятся к решениям, соответственно, дифференциальных задач (4) и (5) с  порядком O(h2). Следовательно, существуют положительные постоянные  и ,  что

,                                                     (21)   

            Значения  и  аппроксимируются   и  с точностью   , следовательно, существуют  положительные числа  и ,  не зависящие от , что

              .                       (22)

            Для задачи (4) известна априорная оценка

                                                         (23)

Из априорной оценки    , с учётом (10),   находим:     

                                                    (24)

            Получим нижнюю оценку выражения  

            Пусть  . В этом случае, в силу (15),    Оценим сверху максимальное значение   функции . Если  где , то в силу  из уравнения (14) получаем оценку .  Если  то из левого краевого условия (14) следует оценка  .  Если  то из правого краевого условия (14) следует оценка  .  Тогда

.

            Пусть  . В этом случае минимальное значение  функции   на сетке  имеет оценку: , где  и, следовательно:   

                                                 

            Таким образом, при условии (3)

                                                     (25)

Применяя оценки  (11), (12), (21)-(25), из (20) получаем:

,       (26)

где

            Из (26) следует утверждение теоремы 3.

При  выражения ,  в силу принципа максимума для задач (5) и (14)  соответственно. В этом случае остаются в силе все изложенные выше результаты, и, как следует из (26),  . Эти результаты аналогичны результатам, полученным в работе .

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

Рецензенты:  

Шхануков-Лафишев М.Х.,  д.ф.-м.н., профессор, ФГБУН «Институт информатики и проблем регионального управления   Кабардино-Балкарского научного центра РАН», г. Нальчик;

Ашабоков Б.А., д.ф.-м.н., профессор, ФГБОУ ВПО «Кабардино-Балкарский государственный аграрный университет им. В.М. Кокова», г. Нальчик.


Библиографическая ссылка

Абрегов М.Х., Бечелова А.Р., Нахушева Ф.М. ЧИСЛЕННЫЙ МЕТОД РЕШЕНИЯ КРАЕВОЙ ЗАДАЧИ ВТОРОГО РОДА ДЛЯ НАГРУЖЕННОГО ОПЕРАТОРА ШТУРМА-ЛИУВИЛЛЯ // Современные проблемы науки и образования. – 2015. – № 1-1. ;
URL: https://science-education.ru/ru/article/view?id=18383 (дата обращения: 24.09.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074