
- •Оглавление
- •Предисловие
- •Лекция 1. Понятие погрешности решения систем алгебраических уравнений
- •1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
- •Источники погрешностей
- •1.2. Абсолютные и относительные погрешности
- •1.3. Понятие о системе линейных уравнений
- •1.4. Матричные уравнения
- •Лекция 2. Точные методы решения систем линейных алгебраических уравнений (слау)
- •2.1. Формулы Крамера для решения системы линейных уравнений
- •2.2. Решение системы линейных уравнений методом последовательного исключения неизвестных (метод Гаусса)
- •Пусть дана система
- •Обозначим
- •Лекция 3. Приближенные методы решения слау
- •3.1. Понятие предела для векторов и матриц Пусть дана последовательность векторов
- •Аналогично, если имеется последовательность квадратных матриц
- •3.2. Метод простой итерации Пусть дана система линейных уравнений:
- •В матричном виде:
- •Условия сходимости итерационного процесса
- •3.4. Оценка погрешности приближенного процесса методом итерации
- •Лекция 4. Метод зейделя
- •4.1. Условия сходимости процесса Зейделя
- •Пусть дана линейная система
- •4.2. Оценка погрешности процесса Зейделя
- •4.3. Приведение системы линейных уравнений к виду, удобному для итераций
- •Лекция 5. Методы решения нелинейных уравнений
- •5.1. Итерационные методы решения нелинейных уравнений а). Метод простой итерации.
- •Число итераций при использовании этого метода
- •Если корень нас не устраивает, то мы находим
- •Г). Метод Ньютона.
- •5.2. Метод Рыбакова
- •5.3. Метод наискорейшего спуска
- •Лекция 6. Интерполирование и экстраполирование
- •6.1. Интерполяционный многочлен Лагранжа
- •Перепишем этот многочлен в другой форме
- •Окончательно получаем
- •Лекция 7. Конечные разности
- •Рассмотрим некоторые свойства конечных разностей
- •7.1. Первая интерполяционная формула Ньютона для равноотстоящих узлов интерполяции
- •7.2. Вторая интерполяционная формула Ньютона
- •7.3. Оценка погрешностей интерполяционных формул Ньютона
- •Лекция 8. Линейное интерполирование
- •8.1. Интерполирование по Эйткину
- •8.2. Интерполяция и приближение сплайном
- •Приближение линейными сплайнами
- •Лекция 9. Численное интегрирование функций
- •9.1. Постановка задачи
- •9.2. Формула прямоугольников
- •Лекция 10. Численные методы решения задачи коши для обыкновенных дифференциальных уравнений
- •10.1. Задача Коши. Общие замечания. Постановка задачи
- •10.2. Метод Эйлера
- •10.4. Метод Рунге-Кутта
- •Лекция 11. Разностный метод решения уравнений математической физики (метод сеток)
- •11.1. Метод сеток для уравнения параболического типа
- •Лекция 12. Метод сеток для уравнения гиперболического типа
- •12.2. Метод сеток для уравнений Пуассона и Лапласа
- •Образцы выполнения контрольных работ Контрольная работа № 5
- •Контрольная работа № 6
- •Решение уравнения методом Эйлера
- •Решение уравнения модифицированным методом Эйлера
- •Решение уравнения методом Рунге-Кутта
- •Аналитическое решение заданного уравнения
- •Сравнение точного решения и приближенных решений исходного дифференциального уравнения
- •Контрольная работа №7
- •Контрольная работа №8
- •Решение дифференциального уравнения методом Рунге-Кутта
- •Аналитические решения заданного уравнения
- •Примеры заданий для самостоятельной работы Задание 1
- •Задание 2
- •Задание 3
- •Задание 4
- •Задание 5
- •Варианты задания
- •Задание 6
- •Варианты задания
- •Задание 7
- •Варианты задания
- •Задание 8
- •Варианты задания
- •Задача для самостоятельного решения
- •Вопросы для самостоятельной подготовки
- •Библиографический список
Лекция 4. Метод зейделя
4.1. Условия сходимости процесса Зейделя
Метод Зейделя представляет собой некоторую модификацию метода последовательных приближений. В методе Зейделя при вычислении (k + 1) приближения неизвестного xi учитываются уже найденные ранее (k+1) приближение неизвестных.
Пусть дана линейная система
(4.1)
Выбираем
произвольно начальное приближение
корней
и подставляем в первое уравнение системы
(4.1)
,
полученное первое приближение подставляем во второе:
.
Полученные первые приближения х1(1) и х2(1) подставляем в третье уравнение системы (4.1)
и т. д.
.
Аналогично строим вторые и третьи итерации.
Таким образом, предполагая, что k приближение корней хki известно, по методу Зейделя строим (k + 1) приближение
где k = 0, 1, 2, …, n.
Пример. Методом Зейделя решим систему:
Приведем систему к нормальному виду:
За нулевые значения возьмем соответствующие значения свободных членов
.
Строим итерации по методу Зейделя
Второе приближение
И т.д.
-
№ итерации
0
0,19
0,97
-0,14
1
0,2207
1,0703
-0,1915
2
0,2354
1,0988
-0,2118
3
0,2424
1,1088
-0,2196
4
0,2454
1,1124
-0,2226
5
0,2467
1,1138
-0,2237
6
0,2472
1,1143
-0,2241
7
0,2474
1,1145
-0,2243
8
0,2475
1,1145
-0,2243
Построенный процесс заканчивается, когда с заданной степенью точности получаем одинаковые значения в двух итерациях подряд.
;
Процесс Зейделя для линейной системы Х = + Х также, так и процесс последовательных приближений, сходится к единственному решению при любом выборе начального приближения, если какая-нибудь из норм матрицы меньше единицы. То есть
либо
,
либо
.
Процесс Зейделя сходится к единственному решению быстрее процесса простой итерации.
4.2. Оценка погрешности процесса Зейделя
Пусть дана линейная система
x = + х ,
где xi – точное значение корней линейной системы, xi(k) – k приближение, вычисленное по методу Зейделя. Тогда оценка погрешности этого метода делается по формуле
;
.
Пример. Подсчитать, сколько итераций по методу Зейделя необходимо выполнить, чтобы с точностью до 10-4 найти корни системы.
Решение:
;
;
.
Значит
Или
Имеем :
,
т.е.
,
по формуле (3.8) определяем
,
,
.
4.3. Приведение системы линейных уравнений к виду, удобному для итераций
Процессы последовательных приближений и метод Зейделя для линейных систем х = + сходятся к единому решению, независимо от выбора начального вектора, если
или
.
Таким образом, для сходимости вышеуказанных итерационных процессов достаточно, чтобы значения элементов ij при i j были небольшими по абсолютной величине. Это равносильно тому, что если для линейной системы АХ = В модули диагональных коэффициентов каждого уравнения системы больше суммы модулей всех остальных коэффициентов (не считая свободных членов), то итерационные процессы для этой системы сходятся, т.е. мы имеем систему
.
Причем,
если
то процессы последовательных приближений
и Зейделя для данной системы сходятся.
Применяя элементарные преобразования,
линейную систему АХ
= В
можно заменить такой эквивалентной
системой Х
= р
+ Х,
для которой условия сходимости будут
выполнены.