otchet_lab_4 (2) (1274776)
Текст из файла
НИУ МЭИ
ЛАБОРАТОРНАЯ РАБОТА 4
РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
УРАВНЕНИЙ ИТЕРАЦИОННЫМИ МЕТОДАМИ
Студент : Нгуен Тхань Там
Группа : Тф-13-14
Вариант : 32
Задача : 4.1, 4.2
Задача 4.2
-
постановку задачи :
Для системы уравнений
с симметричной положительно определенной матрицей найти решение методом Зейделя с точностью
, взяв нулевое начальное приближение. При программировании учесть разреженность матрицы A.
На главной диагонали элементы равны 300 , на 4-ой наддиагонали и шестой наддиагонали элементы равны 10. b=i*e^(22/i)*sin(9/i) N = 30
-
Метод Зейделя
Метод можно рассматривать как модификацию метода Якоби. Основная идея состоит в том, что при вычислении очередного (n+1)-го приближения к неизвестному xi при i >1 используют уже найденные (n+1)-е приближения к неизвестным x1, x2, ..., xi - 1, а не n-ое приближение, как в методе Якоби. Расчетная формула метода в покоординатной форме записи выглядит так:
,
i = 1, 2, ... m.. Условия сходимости и критерий окончания итераций можно взять такими же как в методе Якоби.
-
Решение задачи
-
Составить расчетные формулы покоординатной формы записи итерационного метода для индивидуального варианта.
-
На главной диагонали элементы равны 1200 , на 3-ой наддиагонали элементы равны 30, на 25-ой надддиагонали элементы равны 40.
40x24 + 30x46 1200x49 = b49
40x25 + 30x47 1200x50 = b50
Преобразуем систему к виду удобному для итерации :
X1 = -0.025x4 – 0.033x26 + 0.000833b1
X2 = -0.025x5 – 0.033x27 + 0.000833b2
X3 = -0.025x6 – 0.033x28 + 0.000833b3
X4 = -0.025x1 – 0.025x7 - 0.033x29 + 0.000833b4
.....................................................................
X25 = – 0.033x1 - 0.025x22 – 0.025x28 + 0.000833b25
......................................................................................................
X47 = – 0.033x22 - 0.025x44 – 0.025x50 + 0.000833b47
X48 = – 0.033x23 - 0.025x45 + 0.000833b48
X49 = – 0.033x24 - 0.025x46 + 0.000833b49
X50 = – 0.033x25 - 0.025x47 + 0.000833b50
В покоординатной форме записи метод Зейделя примет следующий вид:
I = 1,2,3
Xi(n+1) = -0.025xi+3(n) – 0.033xi+25(n) + 0.000833bi
I = 4,5,...,24
Xi(n+1) = -0.025xi-3(n+1) – 0.025xi+3(n) - 0.033xi+25(n) + 0.000833bi
I = 25,26 .... 47
Xi(n+1) = – 0.033xi-25(n+1) - 0.025xi-3(n+1) – 0.025xi+3(n) + 0.000833bi
I = 48,49,50
Xi(n+1) = – 0.033xi-25(n+1) - 0.025xi-3(n+1)+ 0.000833bi
Выберем вектор решения x : x(0) = b
-
С
оставить программу вычисления решения системы методом Зейделя с заданной точностью
с учетом выведенных формул п.1. В программе предусмотреть подсчет количества итераций, потребовавшихся для достижения заданной точности
-
Составить тестовый пример и отладить программу на тестовом примере.
-
Решить указанную задачу.
преобразование системы
к виду
С
оздать субматрицы
Выберем вектор решения x : x(0) = 0.000833b = c
К
ритерий окончания итераций
Используя программу решения задачи
Характеристики
Тип файла документ
Документы такого типа открываются такими программами, как Microsoft Office Word на компьютерах Windows, Apple Pages на компьютерах Mac, Open Office - бесплатная альтернатива на различных платформах, в том числе Linux. Наиболее простым и современным решением будут Google документы, так как открываются онлайн без скачивания прямо в браузере на любой платформе. Существуют российские качественные аналоги, например от Яндекса.
Будьте внимательны на мобильных устройствах, так как там используются упрощённый функционал даже в официальном приложении от Microsoft, поэтому для просмотра скачивайте PDF-версию. А если нужно редактировать файл, то используйте оригинальный файл.
Файлы такого типа обычно разбиты на страницы, а текст может быть форматированным (жирный, курсив, выбор шрифта, таблицы и т.п.), а также в него можно добавлять изображения. Формат идеально подходит для рефератов, докладов и РПЗ курсовых проектов, которые необходимо распечатать. Кстати перед печатью также сохраняйте файл в PDF, так как принтер может начудить со шрифтами.
оставить программу вычисления решения системы методом Зейделя с заданной точностью
с учетом выведенных формул п.1. В программе предусмотреть подсчет количества итераций, потребовавшихся для достижения заданной точности 














