Документ взят из кэша поисковой машины. Адрес оригинального документа : http://mph.cmc.msu.ru/Home/Opus/a65.doc Дата изменения: Tue Apr 6 14:56:18 2010 Дата индексирования: Mon Oct 1 19:32:51 2012 Кодировка: koi8-r

Доклады Академии наук СССР 1960. Том 131, ? 6

А. Н. Тихонов, А. А. Самарский

О коэффициенто-устойчивости разностных схем



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

п. 1. Рассмотрим на отрезке [pic] класс краевых задач
[pic] (1)
Коэффициенты уравнения принадлежат классу [pic] кусочно-непрерывных
функций и удовлетворяют условиям:
[pic], (2)
где [pic] и [pic]- положительные постоянные.
Пусть [pic]= { [pic]} - равномерная разностная сетка с шагом [pic], а
[pic] - однородная трехточечная разностная схема, соответствующая оператору

[pic]; (3)
[pic]
Функционалы [pic], [pic], [pic] и [pic] удовлетворяют условиям [pic]
работы [1], т. е. мы рассматриваем тот же исходный класс разностных схем,
что и в работе [1]. При этом предполагается, что [pic] и [pic] - линейные
функционалы.

п. 2. Если [pic], то разностный оператор [pic] называется
консервативным. Консервативный оператор можно записать в самосопряженной
форме
[pic][pic],
где [pic] (3()
Заметим, что разностная схема (3) может быть консерватизирована путем
умножения на множитель
[pic]. (4)
В результате получим консервативную, вообще говоря, неоднородную
схему.

п. 3. Разностная функция Грина [pic] определяется условиями
[pic] (5)
Решение краевой задачи
[pic] (6)
дается формулой
[pic]. (7)
Функция Грина [pic] удовлетворяет следующему условию «симметрии» [pic],
где [pic] дается формулой (4). Для консервативного оператора [pic], и мы
получаем условие симметрии [pic].

Лемма 1. Если коэффициенты [pic] из класса [pic]удовлетворяют условиям
(2) , а [pic] - исходная, разностная схема вида (3), то разностная функция
Грина [pic], определяемая условиями (5), и ее первые разностные отношения
[pic]
ограничены по абсолютной величине постоянной, зависящей только от
[pic],[pic].

п. 4. При решении разностных краевых задач может оказаться, что
коэффициенты разностных уравнений по тем или иным причинам определяются
неточно. Однако желательно, чтобы при малом искажении коэффициентов решение
задачи менялось мало.
Пусть [pic] и [pic] - решения разностных краевых задач
[pic] (8)
[pic] . (9)
При этом коэффициенты уравнения искажаются либо за счет искажения
коэффициентов дифференциального уравнения, либо за счет неточности
вычисления функционалов [pic], [pic], [pic] и [pic], либо, наконец, в
результате обеих указанных причин.
Будем говорить, что разностная схема (3) удовлетворяет принципу ко-
устойчивости, если из условий
[pic] (10)
где [pic] при [pic], следует сходимость решения [pic], разностной краевой
задачи (9) к решению [pic] задачи (1), т. е.
[pic] при [pic] . (11)
Отсюда, в частности, следует, что для ко-устойчивой схемы
[pic]
где [pic] при [pic].
Если в условиях (10) и (11) заменить
[pic], (12)
то мы получим принцип ко-устойчивости n-го ранга.

п. 5. Формулируем необходимое условие ко-устойчивости.
Пусть [pic] имеет разрыв в точке [pic], [pic], так что [pic]. Введем
функцию [pic], совпадающую с [pic] всюду, кроме интервалов [pic] и [pic].
Тогда для ко-устойчивости схемы [pic] вида (3) необходимо выполнение
условия
[pic] при [pic]. (13)
Нетрудно заметить, что необходимое условие сходимости в классе разрывных
коэффициентов, полученное ранее в работах [3,4], является следствием
необходимого условия ко-устойчивости (13) (при [pic]).

п. 6. Лемма 2. Любая консервативная схема [pic] из исходного семейства
схем удовлетворяет необходимому условию ко-устойчивости (13).

Лемма 3. Пусть[pic],[pic] - решения краевых задач
[pic]
где [pic]- консервативные разностные операторы вида (3(), коэффициенты
которых удовлетворяют условиям
[pic] . (2()
Тогда имеет место неравенство
[pic] (14)
где С - постоянная, зависящая только от [pic]и[pic].

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

Теорема 1. Если консервативная схема [pic] из исходного семейства имеет
в некотором классе [pic] n-й интегральный порядок точности, то она имеет
этот же n-й порядок точности для коэффициентов из класса[pic]*, т. е.
для [pic] [pic].


п. 7. Рассмотрим теперь каноническую схему (см. [1])
[pic] , (15)
функционалы которой не зависят от h3, и потребуем, чтобы она удовлетворяла
необходимому условию (13).

Теорема 2. Если каноническая разностная схема (15) из исходного
семейства схем удовлетворяет необходимому условию ко-устойчивости (13), то
она консервативна, т. е. [pic] или [pic].

Опираясь на теорему 2 и леммы 2 и 3, можно убедиться в том, что:

Теорема 3. Любая однородная консервативная схема из исходного семейства
схем удовлетворяет принципу ко-устойчивости.

В результате приходим к следующей основной теореме.
Теорема 4. Необходимым и достаточным условием ко-устойчивости
канонической схемы [pic] является ее консервативность.

Теорема 5. Любая консервативная схема [pic] имеет 1-й интегральный
порядок точности в классе[pic]*.

п. 8. Потребуем теперь, чтобы разностная схема [pic] удовлетворяла
необходимым условиям ко-устойчивости 2-го ранга.

Теорема 6. Существует единственная каноническая схема («наилучшая
консервативная схема»), имеющая 2-й интегральный порядок точности в [pic] и
удовлетворяющая принципу ко-устойчивости 2-го ранга; эта схема [pic]
консервативна и определяется при помощи функционалов
[pic]. (16)
Заметим, что при доказательстве этой теоремы используется, в частности,
лемма 1.

п. 9. Заменяя интеграл, определяющий [pic] в формуле (16), схемой по
какой-нибудь квадратурной формуле, мы получим вместо наилучшей канонической
схемы [pic] неканоническую схему
[pic], где [pic] ,
[pic].

Теорема 7. Для того чтобы, определенная выше неканоническая схема [pic]
имела второй интегральный порядок точности в классе [pic], необходимо и
достаточно, чтобы [pic] при [pic] ([pic]).
Аналогичная теорема имеет место и для схемы [pic], функционалы [pic] и
[pic] которой вычисляются по аналогии с функционалом [pic].

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

Литература

1. А.Н.Тихонов, А.А.Самарский, ДАН, т. 131, ?4, 1960.
2. А.Н.Тихонов, А.А.Самарский, ДАН, т. 131, ?3, 1960.
3. А.Н.Тихонов, А.А.Самарский, ДАН, т. 108, ?3, 1956.
4. А.Н.Тихонов, А.А.Самарский, ДАН, т. 122, ?4, 1959.
* [pic] - класс функций, имеющих на отрезке [0,1] непрерывную производную
= -го порядка, удовлетворяющую условию Гельдера порядка [pic].
* [pic] - класс функций, имеющих на отрезке [0,1] m кусочно-непрерывных
производных, причем m -я производная удовлетворяет в интервалах ее
непрерывности условию Гельдера порядка [pic].