
- •Модуль 1 Основы теории разностных схем
- •Комплексная цель модуля
- •Введение
- •1.3. Сетки и сеточные функции. Вводные понятия
- •1.4. Построение разностных схем
- •1.4.1. Интегро-интерполяционный метод
- •1.4.2. Консервативные и не консервативные разностные схемы
- •1.4.3. Исследование аппроксимации дифференциального уравнения
- •1.4.4 . Исследование погрешности аппроксимации граничного условия
- •1.4.5. Понятие о погрешности аппроксимации, устойчивости и сходимости разностной схемы на примере краевой задачи для оду второго порядка
- •1.4.6. Разностные операторы и их свойства
- •1. Разностные тождества
- •1.1. Формулы разностного дифференцирования произведения
- •1.2. Первая разностная формула Грина
- •1.3. Вторая разностная формула Грина
- •2. Оператор разностной производной второго порядка и его собственные значения
- •3. Обобщение оценок для скалярных произведений, содержащих оператор второй разностной производной с переменными коэффициентами
- •1.4.7. Построение разностных схем для уравнения теплопроводности
- •1.4.8. Связь аппроксимации, устойчивости, сходимости для линейных разностных задач
- •1.4.9. Исследование устойчивости разностных задач методом гармоник
- •1.4.10. Неявная схема для уравнения теплопроводности
- •1.4.11. Схема с весами
- •1.4.12. Трехслойные разностные схемы
- •Исследование устойчивости методом гармоник однородного уравнения с нулевыми граничными условиями. Необходимое условие устойчивости
- •1.4.13. Примеры трехслойных разностные схем для уравнения теплопроводности
- •1.5. Принцип максимума и следствия из него
- •1.5.1. Аппроксимация задачи Дирихле для уравнения Пуассона в прямоугольнике
- •1.5.2. Каноническая форма сеточных уравнений
- •1.5.3. Каноническая форма общего вида для сеточного уравнения
- •Разностная схема, аппроксимирующая задачу Дирихле для уравнения Пуассона в прямоугольнике:
- •Принцип максимума
- •1.5.4. Следствия из принципа максимума
- •1.5.5. Теорема сравнения . Мажоранта
- •1.5.6. Следствие из теоремы сравнения
- •1.5.7. Оценка решения неоднородного уравнения
- •1.5.8. Исследование устойчивости разностных схем на основе принципа максимума и следствий из него
- •1.5.9. Примеры решения задач на принцип максимума и следствий из него
- •1.6. Энергетический метод исследования устойчивости
- •1.6.1. Некоторые сведения из теории операторов
- •Знакоопределенные операторы
- •1.6.2. Свойства оператора симметричной разностной производной второго порядка в гильбертовом пространстве
- •Свойства разностных эллиптических операторов
- •Построение энергетической эквивалентности оператора второй разностной производной с переменными коэффициентами
- •Корректность операторно-разностных уравнений
- •1.7.1. Общие сведения
- •1.7.2. Краткие сведения об обратимости оператора в линейных нормированных пространствах
- •1.6.3. Свойства операторов разностных производных первого порядка
- •Каноническая форма двухслойных схем и ее применение к исследованию устойчивости
- •1.7.4. Необходимое и достаточное условие устойчивости в на
- •Энергетическое пространство
- •1.7.5. Исследование устойчивости в случае несамосопряженного оператора
- •Экономичные схемы для нестационарных задач математической физики
- •Схемы переменных направлений
- •Факторизованные схемы
- •Понятие о монотонных схемах
- •Понятие асимптотической устойчивости.
- •Понятие о трехслойных разностных схемах.
- •Введение в метод конечных элементов
- •Кусочно-линейные восполнения сеточных функций.
- •Понятие о методе конечных элементов.
- •10. Общее описание метода.
- •20. Пример построения и исследования конечно-элементой схемы.
- •Свойства приближенного решения мнк.
- •10. Вводные сведения.
- •20. Вариационные свойства приближенного решения.
- •30. Формулы вложения.
- •Сходимость приближенного решения к точному решению.
- •1.9. Заключение по модулю
- •Вопросы для самоконтроля
- •1.10. Проектное задание
- •Контрольная работа №1.
- •Контрольная работа №2.
- •1.11. Тест рубежного контроля
- •Бланк правильных ответов
1.4.6. Разностные операторы и их свойства
1. Разностные тождества
1.1. Формулы разностного дифференцирования произведения
Известна непрерывная
формула дифференцирования произведения:
.
Она имеет следующий разностный аналог:
. (1)
. (2)
Выведем формулу (1).
.
Выведем формулу (2).
.
Формулы суммирования по частям.
.
Далее удобно ввести обозначения трех типов скалярных произведений сеточных функций.
;
;
. (3)
______________________________________________________________
‼ Проверить аксиомы скалярного произведения самостоятельно.
______________________________________________________________
На сетке
строится разностная схема:
Таким образом, получим:
Если граничные условия 2-го и 3-го рода, то введение граничных условий 1-го рода сопровождается изменением самого оператора.
. (3)
. (4)
Докажем, например, формулу (3′).
Воспользуемся формулой:
.
.
__________________________________________________________________
‼ Доказать формулу (4) самостоятельно.
__________________________________________________________________
1.2. Первая разностная формула Грина
- непрерывный аналог.
. (5)
Пусть в формуле
(3′)
,
тогда
получаем (5).
.
Частный случай формулы Грина, когда
.
(6)
1.3. Вторая разностная формула Грина
Непрерывный аналог этой формулы:
.
Эту формулу легко получить базируясь на формуле (5) – первой формуле Грина.
. (7)
В соответствии с
формулой (5) имеем, заменив
,
. (8)
Из (5) вычтем (8).
так как
.
Из формулы (7)
следует, что оператор
является самосопряженным в пространстве
сеточных функций, определенных на сетке
и обращающихся в ноль в узлах
,
т.к. все выражения в круглых скобках в
правой части (7) обращаются в ноль, откуда
.
2. Оператор разностной производной второго порядка и его собственные значения
Далее будем
рассматривать оператор второй разностной
производной вида:
.
Этот оператор
является положительно определенным в
пространстве
,
то есть для любой функции
выполняется неравенство
.
Значит, существует
обратный оператор
и уравнение
однозначно разрешено. Этот факт следует
из соотношения (6). Положим в формуле
(6)
,
тогда имеем
.
Следовательно,
что
и требовалось доказать.
Рассмотрим классическую задачу на собственные значения для ОДУ II-го порядка вида:
. (1)
. (2)
Нужно найти
и
,
которые удовлетворяют уравнению.
Сделаем предположение
о том, что
.
Иначе среди вещественных функций не
найдется таких, которые удовлетворяют
условию (2).
При
,
любое.
.
.
Проведем нормировку:
.
,
иначе
.
Разностный аналог задачи (1), (2) имеет вид:
. (3)
. (4)
Будем искать
решение задачи (3), (4) в виде:
.
.
. (5)
. (6)
. (7)
. (8)
(9)
. (9′)
Потребуем, чтобы удовлетворялись граничные условия (4), если
.
.
.
(10)
Подставим (10) в (9′), получим
.
Собственные
функции, с точностью до коэффициента,
для дискретной задачи с учетом нормировки
совпадают с собственными функциями
непрерывной задачи. Собственные значения
- иные, нетрудно видеть, что они меньше
собственных чисел непрерывной задачи.
Даже для К=1 собственные значения
дискретной задачи меньше собственных
значений непрерывной задачи
.
Собственные функции дискретной задачи после нормировки совпадают с собственными функциями непрерывной задачи и имеют вид:
.
Если устремить
шаг
к нулю, мы получим неопределенность:
.
Возьмем к=1 – наименьшее собственное значение:
.
Найдем оценку для наибольшего собственного значения дискретной задачи, соответствующего к=N-1:
.
Для простоты
считаем, что
,
тогда имеем следующий результат:
.
______________________________________________________________
‼ Показать
самостоятельно, что
.
Найти наибольшую
______________________________________________________________