Лабораторная работа №2 Тема. Численные методы решения задач линейной алгебры, метод Гаусса

Задание

1. Решить систему линейных алгебраических уравнений (СЛАУ) .

2. Вычислить определитель и обратную матрицу для матрицы А методом исключения Гаусса.

3. Сделать выводы о корректности задачи (существование, единственность, устойчивость решения относительно исходных данных).

Порядок выполнения работы

1. Для расчета используйте матрицу А из приложения 2. Вектор свободных членов задайте произвольно.

2. Решите полученную СЛАУ методом Гаусса.

• прямой ход: привести СЛАУ к эквивалентной системе с треугольной матрицей системы, т.е.

• обратный ход: последовательно вычислить неизвестные х_n,..,х₂, х₁

3. Проверьте полученное решение СЛАУ, используя надстройку Excel поиск решения применительно к исходной системе .

4. Вычислите вручную определитель матрицы А методом Гаусса.

5. Вычислите вручную матрицу А^-1, обратную матрице A, методом Гаусса. Проверьте расчеты на ЭВМ, используя матричные функции Excel. Проанализируйте полученные результаты.

6. Вычислите нормы матриц А и А^-1 (можно вручную).

7. Исследуйте обусловленность матрицы, вычислив меру обусловленности (А). Сделайте заключение об обусловленности матрицы A и заданной системы.

8. Задайте небольшое возмущение исходных данных (только один элемент матрицы А, (~0.1)) и снова решите систему, используя надстройку поиск решения. Проанализируйте, как изменились результаты.

9. Проанализировав полученные данные, сделайте заключение о корректности исходной задачи (существование, единственность, устойчивость решения).

Решение слау с использованием приложения Microsoft Excel

Пример 2.1. Найти решение системы линейных алгебраических уравнений (2.1) используя алгоритм метода Гаусса

, (2.1)

Расчетная схема метода Гаусса приведена на рис (2.1).

Рис.2.1.

Пример 2.2: Сделать проверку полученного решения СЛАУ из примера 2.1, используя надстройку Поиск решения.

При решении СЛАУ с помощью надстройки приложение Excel использует итерационные (приближенные) методы. Строится последовательность приближений,i=0,1,…n. Назовем векторомневязок вектор:

(2.2)

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

Последовательность действий:

1. Заготовьте таблицу, как показано на рис.2.2.

Рис.2.2.

2. Заготовьте ячейки А7:С7, где будет сформировано решение системы (х₁, х₂, х₃). Первоначально они остаются пустыми, т.е. равными нулю. Однако для контроля правильности вводимых далее формул, удобно ввести в эти ячейки какие-либо значения, например единицы. Эти значения можно рассматривать как нулевое приближение решения системы,.

3. Введите коэффициенты системы (матрицу А) в ячейки А3:С5.

4. В столбец D введите выражения для вычисления левых частей исходной системы . Для этого можно использовать функциюСУММПРОИЗВ, принадлежащую категории Математические.

5. В столбец Е запишите значения правых частей системы (матрицу В).

6. В столбец F введите невязки в соответствии с формулой (2.2). Будет не лишним проверить правильность вычислений для случая .

7. Выберите командуменю Сервис\Поиск решения. В окне Поиск решения (рис.2.3) в поле Изменяемые ячейки укажите блок $А$7:$С$7, а в поле Ограничения – $F$3:$F$5=0. Для этого надо щелкнуть на кнопке Добавить и ввести эти ограничения.

Рис. 2.3.

Полученное решение системы (2.1.) х₁=1; х₂=-1 х₃=2 получено в ячейках А7:С7, рис.2.2.