Будь умным!


У вас вопросы?
У нас ответы:) SamZan.ru

Лабораторная работа 4 по курсу Численные методы Выполнил- студент группы Э42 Воробченко Е

Работа добавлена на сайт samzan.ru: 2016-03-30


Инженерно-технологическая Академия

Южный Федеральный Университет

Лабораторная работа №4 по курсу

«Численные методы»

 

Выполнил: студент группы Э-42

Воробченко Е.А

Вариант №3

Проверил: Фоменко Наталья Алексеевна

Таганрог 2013

1.Формулировка поставленной задачи.

Решите систему уравнений методом Якоби, Зейделя, наименьших невязок и методом скорейшего спуска

3)

2. Описание методов.

Итерационные методы Якоби и Зейделя.

Рассмотрим простейшие итерационные методы решения СЛАУ – методы Якоби (простой итерации) и Зейделя. Будем рассматривать систему

,

(4.1)

где

(4.2)

- матрица, имеющая обратную матрицу, , – соответственно векторы-столбцы неизвестных и правых частей.

Начнем рассмотрение с метода Якоби. Будем считать, что все , , где – элементы матрицы А, лежащие на главной диагонали. Преобразуем систему (4.1) к виду

(4.3)

В равенстве (4.3) значение суммы считается, как это принято, равным нулю, если значение верхнего предела суммирования меньше значения нижнего предела суммирования.

Далее верхний индекс будет указывать номер итерации (приближения), например,

Зададим произвольным образом начальное приближение номера , например, так

Исходя из равенства (4.3), итерационное равенство (формулу) метода Якоби запишем так

,

(4.4)

Окончание итераций определяется:

1) по числу итераций, заданному заранее (максимальному числу итераций);

2) по величине так называемой стабилизации

– заданное число, близкое к нулю;

3) по убыванию нормы вектора невязки на -ой итерации, где

,

(4.5)

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

Напомним, что в качестве нормы вектора можно взять величины: или  

С вычислительной точки зрения способ 3) хотя и более громоздкий по сравнению со способами 1) и 2), однако, он дает наиболее объективную информацию о погрешности полученного решения.

Идея метода Зейделя состоит в том, чтобы найденные значения использовать для вычисления , усовершенствовав формулу (4.4) так:

.

(4.6)

4.2. Каноническая форма итерационных методов.

Для исследования сходимости итерационных методов, т.е. установления справедливости равенства где – точное решение системы (2.1), удобнее записывать эти методы в матричной, а не в координатной форме.

Представим матрицу в виде суммы трех матриц

, где

- диагональная часть матрицы А.

Очевидно, метод Якоби с использованием введенных обозначений в векторной форме принимает вид

,

где – матрица, обратная к матрице :

 

Метод Якоби также можно записать следующим образом:

.

(4.7)

Аналогичным образом из соотношений (4.6) можно получить представление метода Зейделя в векторной форме:

.

(4.8)

Далее мы увидим, что векторные равенства (4.7) и (4.8) являются частными случаями так называемой канонической формы одношаговых (двухслойных) итерационных схем вида

(4.9)

где – квадратная невырожденная матрица размера , называемая стабилизатором, – число, называемое итерационным параметром.

Матрица называется положительно определенной, если скалярное произведение для всех ненулевых векторов , или, что то же самое, .

Сформулируем теорему, принадлежащую А.А. Самарскому.

Теорема. Пусть – симметричная положительно определенная матрица, , и пусть выполнено неравенство для любого ненулевого вектора из -мерного пространства.

.

Тогда итерационный метод (4.9) сходится, т.е.

Покажем, как использовать данную теорему для доказательства сходимости, например, метода Зейделя.

Сравнивая (4.8) и (4.9) приходим к равенствам

, .

Таким образом, если и – положительно определенная матрица, то при условии выполнения неравенства , что является краткой формой записи неравенства относительно скалярных произведений

, ,

метод Зейделя сходится.

Заметим, что

.

(4.10)

Нетрудно проверить, что для любого -мерного вектора

.

(4.11)

C другой стороны, из неравенства   вытекает неравенство

.

(4.12)

В самом деле, выберем , где . Тогда

Поскольку – любое, то все , . Значит, справедливо неравенство (4.12). В силу (4.11) и (4.12) из (4.10) имеем

.

4.3. Вариационно-итерационные методы решения СЛАУ.

Преимущество данных методов – они не используют никакой дополнительной информации об операторе , т.е. и , входящие в оценку и необходимые для выбора , здесь не требуются. Рассмотрим методы минимальных невязок и скорейшего спуска.

1. Метод минимальных невязок.

(4.13)

Для получим равенство, умножив обе части равенства (4.13) на матрицу :

.

Меняя знаки и группируя слагаемые соответствующим образом, получаем:

или

Параметр , будем выбирать из условия минимума невязки по норме

.

Продифференцируем по , получим

,

.

(4.14)

2. Метод скорейшего спуска.

Получается из условия минимума энергетической нормы погрешности где , – точное решение исходной системы. Поскольку , и учитывая, что

, получим

Дифференцируя по , получим

, откуда

 

(4.15)

3.  Программная реализация данных методов.




1. Ключевые компетенции для Европы Берн 1996 г
2. Властивості і якості особистості
3. политический режим
4. 1 Понятие предпринимательства и его роль Понятия предпринимательства и предпринимательской деятельности
5. за технологических или эксплуатационных погрешностей зазоров погрешностей формы контактирующих поверхно
6. Утверждаю.
7. КРЕПЫШ МУНИЦИПАЛЬНОГО ОБРАЗОВАНИЯ ГОРОД НОЯБРЬСК 629800 г
8. История русской культуры VXVI вв
9. Сущность роль и место дисциплины в процессе подготовки к профессиональной деятельности
10. Содержание калькулирования и его роль в управлении производством
11. Мотив карточной игры в рассказе Шаламова На представку
12. Лабораторная работа 4 Задание 1
13. i igi i ii1 i1 элемент памятиблок умножения
14. Міжнародні економічні відносини
15. мышь Проводник Работа с окнами
16. Нові слова як особливість розвитку сучасної англійської мови У статті розглянуто діалектичний розвиток.
17. Введение Межбанковские расчеты являются одним из элементом платежной системы государства
18. 1 Технико ~ эксплуатационная характеристика пригородного участка Пригородные пассажирские пер
19. Производство чугуна и стали
20. Электрооборудование промышленных и гражданских зданий Выполнил- ст