Методы решения систем линейных алгебраических уравнений
![]()
|
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «КАЗАНСКИЙ (ПРИВОЛЖСКИЙ) ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ» ИНСТИТУТ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ И ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ Специальность: «Прикладная математика и информатика» Семестровая работа по курсу «Численные методы» на тему: «Методы решения систем линейных алгебраических уравнений» Работу выполнил: студент 4 курса группы 09-405 Баграмов Р.Р. Проверила: Гнеденкова В.Л. Казань 2017 год ОглавлениеПостановка задачи 2 Решение задачи 3 Реализация метода прогонки 4 Метод Якоби. 6 Метод Зейделя. 8 Метод релаксации. 10 Метод наискорейшего спуска. 13 Другие эксперименты. 16 Экспериментальная часть для метода Якоби. 16 Экспериментальная часть для метода Зейделя. 21 Экспериментальная часть для метода верхней релаксации. 26 Сравнение методов Якоби, Зейделя, верхней релаксации. 29 Листинг. 30 Постановка задачиДана система линейных уравнений, которая имеет следующий вид: ![]() где ![]() ![]() ![]() По условию задачи дано, что ![]() ![]() ![]() Значение функции ![]() Необходимо найти решение данной системы линейных алгебраических уравнений методом прогонки, а также итерационными методами:
Решение задачиМетод прогонки. Метод прогонки для решения системы линейных уравнений в общем виде: ![]() Решением системы (3) является вектор ( ![]() Метод прогонки заключается в том, чтобы каждый ![]() ![]() Формула (4) – формула прогонки. ![]() Исключая ![]() ![]() ![]() Далее из предыдущего равенства выражается: ![]() ![]() ![]() Сравнивая (4) и (5), получаются рекуррентные формулы для вычисления коэффициентов прогонки: ![]() ![]() ![]() ![]() ![]() Выражая из первого уравнения системы (3) ![]() ![]() ![]() ![]() Все оставшиеся коэффициенты прогонки будут искаться по формуле (6). Далее вычисляются все ![]() ![]() ![]() ![]() ![]() ![]() ![]() Остальные компоненты вектора ( ![]() Метод прогонки справедлив только для систем уравнений с трёхдиагональной матрицей. Необходимо убедиться, что исходная система уравнений имеет матрицу c диагональным преобладанием. Действительно, матрица коэффициентов системы (1) в нашем случае имеет вид: ![]() |