- •Вычислительная математика лабораторный практикум
- •Содержание
- •Метод исключения Гаусса
- •Введение
- •Построение алгоритма исключения Гаусса
- •3. Реализация алгоритма Гаусса в Excel
- •4. Реализация алгоритма Гаусса в пакете Mathcad
- •5. Реализация алгоритма Гаусса на языке Turbo Pascal
- •6. Вычисление определителя и обратной матрицы
- •7. Выбор ведущего элемента
- •8. Числа обусловленности
- •9. Задания для самостоятельной работы
- •Контрольные вопросы
- •1.Введение
- •Метод Якоби для решения слау
- •Метод Зейделя для решения слау
- •Задания для самостоятельной работы
- •5. Контрольные вопросы
- •Численные методы решения нелинейных уравнений
- •1. Введение
- •2. Отделение корней уравнения
- •3. Метод дихотомии для решения нелинейных уравнений
- •4. Метод Ньютона для решения нелинейных уравнений
- •5. Задания для самостоятельной работы
- •6. Контрольные вопросы
- •Полиномиальная интерполяция
- •1. Интерполяция данных каноническим полиномом
- •2. Интерполяционный полином Ньютона
- •3. Интерполяционный полином Лагранжа
- •4. Задания для самостоятельной работы
- •Контрольные вопросы
- •Метод наименьших квадратов
- •1. Введение
- •2. Линейная аппроксимация
- •3. Аппроксимация нелинейными функциями
- •4. Аппроксимация полиномом
- •Задания для самостоятельной работы
- •6. Контрольные вопросы
- •1. Введение
- •2. Постановка задачи
- •3. Численное дифференцирование с заданной точностью
- •Модификация алгоритма численного дифференцирования Использование центральной разности (6.3) для приближения производной позволяет проводить вычисления с точность порядка :
- •Результаты вычислений сведем в таблицу:
- •5. Действия над приближенными числами
- •6. Задания для самостоятельной работы
- •Контрольные вопросы
- •1. Введение
- •2. Метод прямоугольников
- •3. Метод трапеций
- •4. Метод парабол
- •5. Вычисление интегралов с заданной точностью
- •Метод Гаусса
- •7. Задания для самостоятельной работы
- •2. Провести расчеты знакомого уже нам интеграла ошибок
- •8. Контрольные вопросы
- •Список литературы
- •Учебное издание
Построение алгоритма исключения Гаусса
Метод исключения Гаусса можно разбить на два этапа, называемых прямым и обратным ходом, соответственно. На первом этапе исходная система (2) преобразуется к верхнетреугольному виду:
(1.3)
где для удобства вместо элементов в правой части системы (2) записаны Индекс в скобках вверху - индекс преобразования прямого хода, - означает, что коэффициенты изменялись на -м шаге прямого хода.
На каждом шаге первого этапа решения СЛАУ используются две формулы преобразований:
нормировка на единицу диагонального элемента (его называют ведущим или главным) и формирование в -ой строке новых коэффициентов
(1.4)
для ;
последовательное исключение всех элементов в -ом столбце, расположенных ниже диагонали и вычисление новых коэффициентов в оставшейся части матрицы:
(1.5)
Нет необходимости в памяти компьютера запоминать все полученные нули и единицы. Вычисления можно выполнять только для элементов с индексами Верхние индексы могут быть опущены, если не требуется запоминать все промежуточные вычисления.
На втором этапе решения из системы (1.3) последовательно, начиная с последнего уравнения, определяются все значения неизвестных
(1.6)
Прежде чем написать программу рекомендуется вручную решить несколько простых систем уравнений:
1. 2.
3. 4.
5. 6.
7. 8.
9. 10.
Прежде чем приступить к изучению дальнейшего материала, необходимо детально разобраться с последовательностью выполняемых операций, составить алгоритм и написать самостоятельно программу решения СЛАУ методом Гаусса.
3. Реализация алгоритма Гаусса в Excel
Рассмотрим реализацию алгоритма Гаусса в Excel на примере системы
.
Дальнейшие вычисления приведем в виде таблицы (таб.1.1).
Таблица 1.1. Вспомогательная таблица для решения системы методом Гаусса
k=0 |
3 |
-1 |
0 |
5 |
|
-2 |
1 |
1 |
0 |
|
2 |
-1 |
4 |
15 |
|
|
|
|
|
k=1 |
1 |
-0,3333 |
0 |
1,666667 |
|
0 |
0,33333 |
1 |
3,333333 |
|
0 |
-0,3333 |
4 |
11,66667 |
|
|
|
|
|
k=2 |
1 |
-0,3333 |
0 |
1,666667 |
|
0 |
1 |
3 |
10 |
|
0 |
0 |
5 |
15 |
|
|
|
|
|
k=3 |
1 |
-0,3333 |
0 |
1,666667 |
|
0 |
1 |
3 |
10 |
|
0 |
0 |
1 |
3 |
|
|
|
|
|
|
|
|
x1 = |
2 |
|
|
|
x2 = |
1 |
|
|
|
x3 = |
3 |