Решение систем линейных уравнений

МЕТОДОМ ПРОСТОЙ ИТЕРАЦИИ

I. ЦЕЛЬ РАБОТЫ

Приобретение навыков решения систем линейных алгебраических уравнений итерационными методами.

II. ТЕОРЕТИЧЕСКАЯ СПРАВКА

Рассмотрим метод простой итерации.

Пусть дана система линейных алгебраических уравнений

решение систем линейных уравнений - student2.ru

с неособенной матрицей ( решение систем линейных уравнений - student2.ru ). Согласно методу простой итерации ее предварительно приводят к виду

решение систем линейных уравнений - student2.ru ,

где решение систем линейных уравнений - student2.ru , т.е. первое уравнение системы разрешили относительно решение систем линейных уравнений - student2.ru , второе - относительно решение систем линейных уравнений - student2.ru и т.д.

Предположим, что известно начальное приближение

решение систем линейных уравнений - student2.ru

к точному решению решение систем линейных уравнений - student2.ru системы. Тогда все следующие приближения найдем по формуле

решение систем линейных уравнений - student2.ru

Если последовательность решение систем линейных уравнений - student2.ru сходится к некоторому предельному вектору решение систем линейных уравнений - student2.ru , то он будет решением системы. Действительно, считая решение систем линейных уравнений - student2.ru при решение систем линейных уравнений - student2.ru , получаем из выражения

решение систем линейных уравнений - student2.ru

равенство решение систем линейных уравнений - student2.ru .

Последовательность решение систем линейных уравнений - student2.ru в методе простой итерации сходится, если для матрицы решение систем линейных уравнений - student2.ru выполняется одно из неравенств

1) решение систем линейных уравнений - student2.ru ;

решение систем линейных уравнений - student2.ru .

III. ЗАДАНИЕ

Найти решение системы линейных уравнений, приведенной в лабораторной работе №1. При решении системы использовать метод простой итерации.

IV. Оформление отчета

В отчете должны быть представлены:

1. Название работы.

2. Постановка задачи.

3. Описание алгоритма (метода) решения.

4. Текст программы с описанием.

5. Результаты работы программы.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Демидович Б.П., Марон И.А. Основы вычислительной математики. - М.: Лань, 2009. 672 с.

2. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. — М.: БИНОМ. Лаборатория Знаний, 2007 636с.

3. Калиткин Н.Н. Численные методы. - М.: Наука, 1978. 512 с.

4. Крылов В.И., Бобков В.В., Монастырный П.И. Вычислительные методы. Т.1. - М.: Наука, 1976. 350 с.

Лабораторная работа № 3

РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ УРАВНЕНИЙ

МЕТОДОМ ЗЕЙДЕЛЯ

II. ЦЕЛЬ РАБОТЫ

Приобретение навыков решения систем линейных алгебраических уравнений итерационными методами.

III. ТЕОРЕТИЧЕСКАЯ СПРАВКА

Рассмотрим метод Зейделя.

Пусть система решение систем линейных уравнений - student2.ru приведена к канонической форме

решение систем линейных уравнений - student2.ru

В методе простой итерации следующее приближение решение систем линейных уравнений - student2.ru находится по предыдущему решение систем линейных уравнений - student2.ru путем подстановки решение систем линейных уравнений - student2.ru в правую часть (1). При этом порядок выбора уравнений значения не имеет.

Согласно методу Зейделя осуществляется разумный выбор порядка уравнений для подстановок и немедленный ввод в вычисления каждого из полученных приближений для неизвестных.

Предположим, что для перехода от приближения решение систем линейных уравнений - student2.ru к решение систем линейных уравнений - student2.ru выбран какой-то порядок привлечения уравнений для подстановок. Изменяя, если необходимо, нумерацию уравнений и неизвестных, можно считать, что уравнения для подстановок берутся в порядке роста их номеров. Для каждого шага порядок привлечения уравнений может быть своим. Перестановка уравнений и изменение нумераций влекут изменение матрицы решение систем линейных уравнений - student2.ru и вектора решение систем линейных уравнений - student2.ru . Чтобы отметить это, обозначим решение систем линейных уравнений - student2.ru и решение систем линейных уравнений - student2.ru для рассматриваемого шага через решение систем линейных уравнений - student2.ru и решение систем линейных уравнений - student2.ru .

Итерация в методе Зейделя выполняется в следующем порядке:

решение систем линейных уравнений - student2.ru

После нахождения вектора решение систем линейных уравнений - student2.ru устанавливается порядок подстановок в уравнения значений решение систем линейных уравнений - student2.ru и переходят к вычислению вектора решение систем линейных уравнений - student2.ru и т.д.

Приведем теперь принцип установления порядка привлечения уравнений для подстановок решение систем линейных уравнений - student2.ru . Можно пытаться улучшить ту составляющую решения, которая найдена наименее точно, чтобы при нахождении всех других составляющих употребить улучшенное ее значение.

О точности решение систем линейных уравнений - student2.ru можно судить по вектору поправки на шаге решение систем линейных уравнений - student2.ru : решение систем линейных уравнений - student2.ru , где решение систем линейных уравнений - student2.ru . Величины поправок составляющих нумеруют в порядке убывания их модулей, и в том же порядке вычисляют составляющие следующего приближения решение систем линейных уравнений - student2.ru , сначала ту составляющую, которая отвечает наибольшей по модулю поправке, и т.д.

Рассмотрим более подробно стационарный метод Зейделя, когда при итерациях порядок уравнений сохраняется, а следовательно, сохраняются решение систем линейных уравнений - student2.ru и решение систем линейных уравнений - student2.ru . Вычисления по-прежнему проводят по формуле (2).

Разложим матрицу решение систем линейных уравнений - student2.ru на сумму двух матриц решение систем линейных уравнений - student2.ru и решение систем линейных уравнений - student2.ru , где

решение систем линейных уравнений - student2.ru , решение систем линейных уравнений - student2.ru

Тогда равенства (2) можно записать в матричной форме в виде

решение систем линейных уравнений - student2.ru .

Отсюда следует, что

решение систем линейных уравнений - student2.ru ,

а так как определитель матрицы решение систем линейных уравнений - student2.ru равен единице и она имеет обратную матрицу, то равенство (2) равносильно

решение систем линейных уравнений - student2.ru

Поэтому стационарный метод Зейделя равносилен методу простой итерации, примененному к системе

решение систем линейных уравнений - student2.ru

Последовательность решение систем линейных уравнений - student2.ru в стационарном методе Зейделя сходится, если для матрицы решение систем линейных уравнений - student2.ru выполняется одно из неравенств

2) решение систем линейных уравнений - student2.ru ;

3) решение систем линейных уравнений - student2.ru .

IV. ЗАДАНИЕ

Найти решение системы линейных уравнений, приведенной в лабораторной работе №1. При решении системы использовать стационарный метод Зейделя.

IV. Оформление отчета

В отчете должны быть представлены:

1. Название работы.

2. Постановка задачи.

3. Описание алгоритма (метода) решения.

4. Текст программы с описанием.

5. Результаты работы программы.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Демидович Б.П., Марон И.А. Основы вычислительной математики. - М.: Лань, 2009. 672 с.

2. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. — М.: БИНОМ. Лаборатория Знаний, 2007 636с.

3. Калиткин Н.Н. Численные методы. - М.: Наука, 1978. 512 с.

4. Крылов В.И., Бобков В.В., Монастырный П.И. Вычислительные методы. Т.1. - М.: Наука, 1976. 350 с.

Лабораторная работа № 4

Наши рекомендации