БлогNot. Mathcad: метод Гаусса-Зейделя и диагональное преобладание

Помощь дата->рейтинг Поиск Почта RSS канал Статистика nickolay.info Домой

Mathcad: метод Гаусса-Зейделя и диагональное преобладание

Вопрос возник по этой заметке, где реализованы прямые и итерационные методы решения систем линейных алгебраических уравнений (СЛАУ) в Mathcad.

Для сходимости как метода Якоби, так и модификации Гаусса-Зейделя, следует приводить матрицу к виду с диагональным преобладанием, это можно всегда сделать, применяя элементарные преобразования над строками и столбцами матрицы.

Приложенный документ содержит модификацию итерационного метода Якоби, известную как метод Гаусса-Зейделя и показывает, как привести матрицу к виду с диагональным преобладанием :)

Исходя из определения, достаточно, чтобы модуль элемента главной диагонали был больше остальных элементов своей строки. Поэтому для исходной матрицы, показанной на скрине, мы поменяли местами первую и вторую строки, а затем вычли из третьей строки вторую (бывшую первую). То же самое не забыли сделать для вектора правой части b, на практике лучше работать с расширенной матрицей системы, в которую последним столбцом дописан вектор правой части b.

Как видно из результата проверки, система решена с точностью не менее 15 знаков после запятой.

Метод Гаусса-Зейделя для решения СЛАУ
Метод Гаусса-Зейделя для решения СЛАУ

 Скачать метод Гаусса-Зейделя для решения СЛАУ, документ Mathcad 14/15 .xmcd в архиве .zip (22 Кб)


теги: mathcad числа

комментарии (0)

02.12.2017, 13:11; рейтинг: 51

  свежие записипоиск по блогукомментироватьстатистика

Наверх Яндекс.Метрика
© PerS
вход