Система, имеющая решение, называется совместной, система, не имеющая решения,-несовместной. Системы называются эквивалентными, если они имеют одни и те же решения.
Систему (1) кратко можно записать в виде матричного уравнения
. (3)
Если определитель матрицы А не равен нулю (матрица неособенная), то решение системы (3) существует и единственно.
Метод итерации можно применять к решению линейной системы (1), если преобразовать ее к специальному виду. Прежде, чем проводить преобразования, необходимо проверить выполнение условий сходимости метода итерации.
Для того, чтобы метод итерации сходился, необходимо, чтобы выполнялось условие сходимости для системы (1): метод итерации сходится, если выполнены неравенства
> (i = 1, 2, ..., n), (4)
т. е. если модули диагональных коэффициентов для каждого уравнения системы больше суммы модулей всех остальных коэффициентов, не считая свободных членов.
Преобразуем систему (1) к специальному виду: разрешим первое уравнение системы относительно x1, второе – относительно x2, третье – относительно x3 и т. д. В результате получим эквивалентную систему
(5)
Если ввести обозначения , , при и при (i, j = 1, 2, …., n), то систему (5) можно представить в виде
, (6)
или в матричной форме
, (7)
где матрица a - это матрица с нулевыми диагональными элементами. В качестве начального приближения можно принять столбец свободных членов (хотя это не является обязательным условием). Матрицу подставляем в правую часть системы (6) и получим . Далее подставим в правую часть (7), получим и т. д. Таким образом, любое (n + 1) приближение вычисляется по формуле
. (8)
Достаточное условие сходимости процесса итерации:
если для приведенной системы (6) выполнено по меньшей мере одно из условий
< 1 (i = 1, 2, 3, …, n) или < 1 (j = 1, 2, 3, …, n), (9)
то процесс итерации (8) сходится к единственному решению этой системы, независимо от выбора начального приближения.
Замечание1. Решение системы может быть получено с заданным числом точных десятичных знаков. Чтобы избежать накопления погрешностей в промежуточных вычислениях, последовательные приближения , i = 0, 1, 2,… вычисляют до совпадения всех требуемых знаков, после чего запасные знаки округляются.
Оценка погрешности для метода итерации:
≤. (10)
Пример.
Методом итерации решить систему уравнений
и оценить число необходимых для этого шагов.
Видно, что условие сходимости (4) для данной системы не выполняется, так как в уравнениях (а) и (b) нет диагонального преобладания. Для применения метода итерации необходимо преобразовать заданную систему к такому виду, для которого будет выполняться условие (4). Для этого умножим уравнение (а) на γ, (b) на δ, сложим оба уравнения и в полученном выражении выберем γ и δ так, чтобы имело место диагональное преобладание.
.
Положив , получим .
С уравнением (b) поступаем аналогично.
.
Положив , получим .
В результате исходная система примет следующий вид:
,
в которой есть диагональное преобладание, т. е. для которой выполняется условие сходимости (4).
Далее, разрешив эту систему относительно диагональных неизвестных, получим так называемую приведенную систему уравнений
к которой можно применять метод итерации (для приведенной системы выполняется условие сходимости (9)).
В качестве начального приближения принимаем столбец свободных членов:
.
Применяя формулу (8), найдем решение системы:
.
=+=,
=+= и т. д.
Решением исходной системы будет матрица .
Решение систем нелинейных уравнений методом итерации
Пусть дана система нелинейных уравнений
, (11)
или в матричной форме
, (12)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.