
- •Министерство образования и науки российской федерации
- •Методические указания
- •Часть 1
- •Р.Б.Салимов
- •1. Численное решение нелинейных уравнений.
- •1.1. Метод деления отрезка пополам.
- •1.2. Метод Ньютона (метод касательных).
- •1.3. Метод простой итерации.
- •2. Методы решения систем линейных алгебраических уравнений.
- •2.1. Метод Гаусса.
- •2.2. Метод обратной матрицы.
- •2.3. Метод прогонки.
- •2.4. Метод простой итерации (метод Якоби).
- •2.5. Метод Зейделя.
- •3. Численные методы решения систем нелинейных уравнений.
- •3.1. Метод простой итерации (метод Якоби) для систем нелинейных уравнений.
- •3.2. Метод Зейделя для систем нелинейных уравнений.
- •3.3. Метод Ньютона решения систем нелинейных уравнений.
- •Литература
2.4. Метод простой итерации (метод Якоби).
Суть вычислений
итерационными методами состоит в
следующем: расчет начинается с некоторого
заранее выбранного приближения
(начального приближения). Вычислительный
процесс, использующий матрицу
,
вектор
системы (2.1) и
,
приводит к новому вектору
:
,
(2.11)
Затем процесс
повторяется, только вместо
используется новое значение
.
На
-м
шаге итерационного процесса получают:
,
(2.12)
При выполнении
некоторых заранее оговоренных условий
процесс сходится при
.
Сходимость метода простой итерации
обеспечивается при выполнении условия
преобладания диагональных элементов
матрицы A:
,
(2.13)
Заданная точность достигается при выполнении условия:
(2.14)
Пример 2.5. Преобразовать систему уравнений:
(2.15)
к виду, пригодному для построения итерационного процесса методом Якоби и выполнить три итерации.
Решение. Достаточное условие сходимости (2.13) выполняется, поэтому начальное приближение может быть любым.
В
-ом
уравнении все члены, кроме
,
переносятся в правую часть:
(2.16)
Задается начальное
приближение
,
которое подставляется в правую часть
(2.16). Если
,
,
,
то результаты
первой итерации:
Результаты первой
итерации
подставляют
в правую часть (2.16) и получают результаты
второй итерации:
Результаты второй
итерации
подставляют
в правую часть (2.16) и получают результаты
третьей итерации:
Определяют достигнутую точность
|
|
A |
B |
C |
|
|
---|---|---|---|---|---|---|
|
1 |
x1 |
x2 |
x3 |
|
|
|
2 |
0,00 |
0,00 |
0,00 |
|
|
|
3 |
1,00 |
-0,33 |
1,00 |
|
|
|
4 |
1,33 |
-1,17 |
1,33 |
|
|
|
5 |
1,86 |
-1,44 |
1,63 |
|
|
|
6 |
2,06 |
-1,76 |
1,83 |
|
|
|
7 |
2,27 |
-1,93 |
1,96 |
|
|
|
… |
… |
… |
… |
|
|
|
20 |
2,66 |
-2,34 |
2,25 |
|
|
|
21 |
2,66 |
-2,35 |
2,25 |
|
|
|
22 |
2,66 |
-2,35 |
2,25 |
|
|
Рис. 2.4. Решение системы линейных алгебраических уравнений методом Якоби с помощью программы Excel |
Порядок решения.
Представить систему в виде (2.16);
Ввести в ячейки A1:C1 заголовки столбцов (рис. 2.4);
В ячейки A2:C2 – начальное приближение 0, 0, 0;
В ячейку A3 – формулу =(7-4*B2+C2)/7
В ячейку B3 – формулу =(-2-2*A2-3*C2)/6
В ячейку C3 – формулу
=(4+A2-B2)/4
Выделить столбцы A, B, C, вызвать контекстное меню Формат ячеек, установить формат числовой и указать число десятичных знаков, соответствующее необходимой точности, т.е. 2;
Выделить ячейки A3:C3 и скопировать формулы в соседние ячейки расположенных ниже строк A4:C4, A5:C5 и т.д. при помощи маркера заполнения. Каждая новая строка содержит результаты очередного приближения;
Продолжать копирование, пока результат не перестанет меняться;
Ячейки A21, B21, C21 содержат решение системы уравнений, соответствующее заданной точности.
Приближенное решение системы с точностью :
,
,