Итерационные методы

Метод простой итерации. При решении итерационными методами каждая переменная системы находится как предел последовательности приближений. Систему линейных уравнений записывают в виде:

(3)

Далее выбирают начальные приближения значений переменных x₁^{(0 )}, x₂⁽⁰⁾, …, x_n^{(0 )} (любые n чисел) и подставляя их в правую часть системы, находят следующие приближенные значения корней x₁⁽¹⁾, x₂⁽¹⁾, …, x_n⁽¹⁾ и т.д. Итерационный процесс продолжается до тех пор, пока не выполнится условие:

k = 1, …, n

Условие сходимости итерационного процесса:

(4)

Для решения систем линейных уравнений в системе MathCad существуют встроенные функции:

augment(A, B) – формирует расширенную матрицу системы, в первых столбцах которой содержится матрица А, а в последних – матрица В;

rref(A) – выполняет прямой и обратный ход метода Гаусса;

submatrix(A, i, j, l, m) – формирует столбец решения системы (то есть формирует матрицу, которая является блоком матрицы А, расположенны в строках с i по j в столбцах с l по m).

В системе MathСad также реализовано решение систем линейных и нелинейных уравнений итерационными методами. Для этого используются функции Given и Find (приложение Б).

Задание 2.1. Исследуйте систему линейных уравнений и, если решение существует, найдите его по формулам Крамера.

Вариант 1

1) 2)

Вариант 2

1) 2)

Вариант 3

1) 2)

Вариант 4

1) 2)

Вариант 5

1) 2)

Вариант 6

1) 2)

Вариант 7

1) 2)

Вариант 8

1) 2)

Вариант 9

1) 2)

Вариант 10

1) 2)

Вариант 11

1) 2)

Вариант 12

1) 2)

Порядок выполнения задания:

1. Установите режим автоматического выполнения вычислений.

2. Введите матрицу системы.

3. Вычислите определитель матрицы системы. Система имеет единственное решение, если определитель отличен от нуля.

4. Вычислите определители матриц, полученных заменой соответствующего столбца столбцом правых частей.

5. Найдите решение системы по формулам Крамера (2).

Задание 2.2. Решить систему уравнений задания 2.1 с помощью обратной матрицы.

Порядок выполнения задания:

1. Введите матрицу системы и вектор-столбец правых частей.

2. Вычислите решение системы по формуле Х=А^-1·В.

3. Проверьте правильность решения умножением матрицы системы на вектор-столбец решения.

Задание 2.3. Решитьсистемулинейных уравнений задания 2.1 методом Гаусса..

Порядок выполнения задания:

1. Присвойте переменной ORIGIN значение, равное единице.

2. Введите матрицу системы и столбец правых частей.

3. Сформируйте расширенную матрицу системы.

4. Приведите расширенную матрицу системы к ступенчатому виду.

5. Сформируйте столбец решения системы.

6. Проверьте правильность решения умножением матрицы системы на вектор-столбец решения.

Задание 2.4. Решить систему уравнений задания 2.1 с помощью встроенных функций Given и Find.

Порядок выполнения задания:

1. Записать служебное слово Given.

2. Ввести систему уравнений.

3. Вывести решение с помощью функции Find.

Задание 2.5. Решитьсистемулинейных уравнений методом простых итераций.

Вариант 1

A = , B =

Вариант 2

A = , B =

Вариант 3

A = , B =

Вариант 4

A = , B =

Вариант 5

A = , B =

Вариант 6

A = , B =

Вариант 7

A = , B =

Вариант 8

A = , B =

Вариант 9

A = , B =

Вариант 10

A = , B =

Вариант 11

A = , B =

Вариант 12

A = , B =

Порядок выполнения задания:

1. Преобразуйте исходную систему Сх=d к виду х=B+Ах.

2. Проверьте достаточное условие сходимости.

3. Введите матрицы А и B.

4. Определите нулевое (начальное) приближение решения.

5. Задайте количество итераций.

6. Введите формулу вычисления последовательных приближений решения и вычислите их.

7. Выведите на экран матрицу приближенных решений.