- •Глава 1. Численные методы.
- •1.1 Структура погрешности в численном анализе.
- •Округление.
- •Погрешность вычисления функции при неточном задании аргументов. (Вопрос вычисления функции при наличии ошибки в вычислении ее аргумента)
- •1.2. Понятие близости в метрическом пространстве.
- •Примеры классов функций и соответствующих нормированных пространств.
- •1.3. Задача приближения функций. Интерполяция.
- •Задача наилучшего приближения в нормированном пространстве.
- •Задача приближения полиномами.
- •Общая задача интерполяции.
- •1.4. Интерполяционный многочлен Ньютона.
- •Пример:
- •Запись интерполяционного многочлена для равноотстоящих узлов.
- •1.5. Среднеквадратичное приближение функции.
- •1.6. Классические ортогональные многочлены и их применение в задачах приближения функций.
- •С другой стороны
- •1.7. Некоторые общие свойства ортогональных полиномов.
- •1.8. Многочлены Чебышева, их свойства и применение в задачах приближения функций.
- •Первые применения многочленов Чебышева к задаче интерполяции.
- •1.9. Равномерное приближение функций на отрезке.
- •1.10 Теоремы о равномерном приближении функций на отрезке.
- •Глава II. Численное интегрирование.
- •2.1. Квадратурные формулы на основе интерполяций.
- •Формулы Ньютона-Котеса.
- •2.2 Квадратурные формулы Гаусса-Кристоффеля.
- •Глава III. Численные методы алгебры.
- •3.1. Принцип сжатых отображений.
- •3.2. Метод простых итераций для функциональных уравнений.
- •3.3.Метод Ньютона.
- •3.4. Численные методы линейной алгебры.
- •3.4.1. Нормы матриц. Спектральные свойства матриц.
- •Для любой квадратной матрицы и любой согласованной матричной нормы имеет место неравенство:
- •Простейшие следствия из определений.
- •3.4.2. Обусловленность матриц и систем уравнений.
- •3.4.3. Итерационные методы решения систем лау.
- •3.4.4 Стационарные итерационные процедуры. Теоремы о сходимости.
- •Некоторые определения.
- •Глава IV. Численные методы решения обыкновенных дифференциальных уравнений.
- •4.1. Численное дифференцирование.
- •4.2.Численные методы решения задачи Коши.
- •4.2.1. Оценка погрешности метода Эйлера.
- •1. Ошибка растет с номером шага k.
- •4.2.2.Методы Рунге-Кутта.
- •4.3. Численные методы решения краевых задач для оду.
- •4.3.1.Задача аппроксимации на сетке.
- •4.3.2.Аппроксимация, устойчивость и сходимость разностных схем.
- •Итерационная процедура устойчива по правой части тогда и только тогда, когда
- •При этом, если выполняется условие
Итерационная процедура устойчива по правой части тогда и только тогда, когда
,
где с не зависит от h (т.е. от N)
Однако, это условие не всегда легко проверить. Поэтому для конкретных итерационных схем вырабатываются и доказываются конкретные достаточные условия. Таковы, например, условия “благонеявного преобладания” для схем прогонки.
Теорема2.(Необходимый спектральный признак устойчивости).
Пусть - собственные числа оператора Rh. Для устойчивости схемы (8) по правой части необходимо выполнение условия:
, (9)
причем константа не зависит от h (отN).
Пусть (9) не выполняется для некоторого собственного значения .То есть, не существует такой константы , для которой (9) выполнялось бы для данного .Фактически, это означает, что вместо линейного ограничения имеем:
, где 0<<1,c1-некоторая константа.
Пусть - соответствующий собственный вектор, т.е.
Оценим по сеточной норме:
.
Из последнего неравенства следует:
Заметим, что по условию на , поэтому
т.е. нарушается условие устойчивости, сформулированное ранее. Происходит экспоненциальный рост ошибки.
Рассмотрим снова сеточное уравнение вида (6)
Теорема 3.(О сходимости разностной схемы (6)).
Пусть конечно-разностная задача (6) однозначно разрешима, аппроксимирует исходную дифференциальную задачу с порядком p относительно h и устойчива. Тогда имеет место сходимость:
,
где - решение сформулированной разностной задачи;
- точное решение дифференциальной задачи, взятое на сетке.
При этом, если выполняется условие
,
то говорят, что имеет место сходимость порядка p.
Согласно условию теоремы имеет место аппроксимация порядка p:
(10)
(11)
- невязка, которая получается при подстановке точного решения в левую часть уравнения. Подставляя в (10):
. (12)
В возмущенном уравнении
в качестве возмущения выберем невязку, т.е. положим , тогда
. (13)
В силу старого определения устойчивости имеем:
. (14)
Уравнения (11) и (13) имеют одинаковые правые части. В силу однозначной разрешимости задачи (6), имеем:
,
подставим в (14)
Таким образом, мы одновременно доказали сходимость и установили, что порядок сходимости совпадает с порядком аппроксимации.