Математическая энциклопедия - нелинейная краевая задача

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

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

Пусть вектор таков, что система уравнений

однозначно определяет вектор y=a(g). при каждом векторе напр., часто целесообразно взять

Тогда задача (1), (2) может быть сведена к системе (операторному уравнению)

где значение при х=Х решения уравнения (1) при начальном условии . При этом для получения значений необходимо численно решать (см. [1], [2]) соответствующие системы дифференциальных уравнений, что и будет означать дискретизацию исходной задачи. Для решения системы (4) можно применять различные итерационные методы решения нелинейных уравнений. Найдя решение g* уравнения (4), можно определить вектор из системы (3) и, решая уравнение (1) при начальном условии у(0)=у*, получить решение задачи (1), (2). В ряде методов наряду со значениями функции используются и значения ее производной, отыскиваемые либо путем интегрирования уравнения

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

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

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

с линейными краевыми условиями

где квадратная функциональная матрица порядка l, В п и Dnчисловые матрицы размеров и , n=0, 1, .... Начальное приближение у п (х)считается известным. В случае метода Ньютона Канторовича уравнения (6), (7) приобретают вид

Если f(x,y), В (у), D (у)дважды непрерывно дифференцируемы по у, задачи (8), (9) корректны и у о (х)достаточно близко к искомому решению у(х), то

где постоянная,(см. [1]).

Упомянутые методы обобщаются и на случай более общих условий вида

где

Многие важные теоретические и прикладные задачи приводят к необходимости решать Н. к. з. (и родственные к ним проблемы) для уравнений и систем уравнений эллиптич. типа (см., напр., [4] [8]). Для такого класса задач основными численными методами являются проекционные методы (проекционно-разностные, вариационно-разностные, конечных элементов) и разностные методы (см. [7] [14]). Их конструкции во многом аналогичны конструкциям соответствующих методов для линейных краевых задач (см. Лапласа уравнение;численные методы решения; Пуассона уравнение;численные методы решения), но получаемые в этих методах дискретные (сеточные) аналоги краевых задач, являющиеся системами нелинейных уравнений

где

часто вызывают дополнительные значительные трудности как в плане анализа самих систем и близости в том или ином смысле их решений к решению исходной задачи, так и в плане затрат вычислительной работы на формирование самих систем (10) и отыскание их решений. Специфика нелинейности краевой задачи заставляет особое внимание уделять выбору банаховых пространств и их сеточных аналогов, в к-рых и сама задача и ее конечномерные аппроксимации поддаются анализу (см., напр., [5] [12], [15]). Но наиболее подробно изучены особенно в алгоритмич. плане лишь численные методы для классов нелинейных задач, родственных в нек-ром смысле линейным задачам (слабонелинейные уравнения, уравнения с ограниченной нелинейностью) и позволяющих обойтись использованием гильбертовых пространств и их конечномерных аналогов (евклидовых пространств).

Пусть евклидово пространство вектор-функций со скалярным произведением пространство линейных самосопряженных положительно определенных операторов, отображающих H в Н;наряду с Нбудет использоваться и евклидово пространство , , отличающееся от Нлишь видом скалярного произведения:

Пусть, далее,

Если оператор непрерывен на и таков, что для любого ис справедливо неравенство

то система (10) имеет по крайней мере одно решение из .Это утверждение (см., напр., [8], [9]) является одним из следствий классич. теорем существования решений у нелинейных операторных уравнений, основанных на топологич. принципах (см., напр., [14]) и на монотонности операторов (см., напр., [6], [11], [14]). В частности, если непрерывен всюду и

для любого и, то система (10) всегда имеет решение и любое ее решение принадлежит

Если же для любых справедливы неравенства

означающие сильную монотонность и Липшиц-непрерывность оператора LN, то система (10) и ее единственное решение может быть найдено с помощью итерационного метода