
- •Глава 1. Интерполирование
- •Постановка задачи интерполирования
- •1.2. Единственность интерполяционного многочлена
- •. Интерполяционный многочлен Лагранжа
- •1.4. Погрешность интерполяционной формулы Лагранжа
- •. Конечные разности и их свойства
- •1.6. Первая интерполяционная формула Ньютона
- •1.7. Вторая интерполяционная формула Ньютона
- •1.8. Интерполяционная формула Гаусса
- •1.10. Численное дифференцирование
- •Глава 2. Численное интегрирование
- •2.1. Общие замечания
- •Глава 3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1. Вводные замечания
- •3.2. Отделение корней
- •3.3. Метод половинного деления
- •3.5. Метод Ньютона (метод касательных)
- •3.6. Комбинированный метод
- •3.7. Метод итерации (метод последовательных приближений)
- •Глава 4. Решение систем линейных уравнений
- •4.1. Метод Гаусса (схема единственного деления)
- •4.2. Вычисление определителей
- •4.3. Вычисление обратной матрицы
- •4.4. Некоторые сведения из линейной алгебры
- •Рассмотрим систему линейных алгебраических уравнений
- •Вычитая второе равенство из первого, получаем
- •4.6. Приведение линейной системы к виду,
- •4.7. Метод Зейделя Пусть система линейных алгебраических уравнений приведена к виду
- •Глава 5. Обработка результатов наблюдений
- •5.1. Постановка задачи
- •5.2. Приближение функции, заданной таблично,
- •Глава 6. Численные методы решения
- •6.1. Методы решения задачи Коши. Вводные замечания
- •6.2. Решение дифференциальных уравнений
- •6.3. Метод Эйлера Пусть дано дифференциальное уравнение
- •6.4. Уточненный метод Эйлера Рассмотрим дифференциальное уравнение
- •6.6. Методы Рунге – Кутта Пусть дано дифференциальное уравнение
- •6.7. Метод Милна
- •6.8. Приближенное решение систем дифференциальных
- •6.9.Краевые задачи для обыкновенных дифференциальных
- •6.10. Решение краевых задач
- •6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
- •Глава 7. Численное решение дифференциальных уравнений в частных производных
- •7.1. Классификация дифференциальных уравнений
- •7.2. Метод сеток решения краевых задач
- •7.3. Погрешность аппроксимации дифференциальных уравнений
- •7.4. Аппроксимация граничных условий
- •7.6. Метод сеток для уравнений
4.2. Вычисление определителей
при помощи метода Гаусса
Рассмотрим для простоты определитель третьего порядка
(4.15)
Пусть . Вынося из первой строки за знак определителя, получим
1
где
(j=2, 3) вычисляются по
формулам (4.3). Умножаем первую строку
определителя последовательно на
и
и вычитаем ее из второй и третьей строк
соответственно:
1
0
=
, (4.16)
0
где
(i=2, 3; j=2, 3) вычисляются по
формулам (4.5).
Этим мы снизили порядок определителя на единицу. Применим к полученному определителю те же преобразования. Выносим из первой строки за знак определителя:
1
,
где вычисляется по формуле (4.7).
Умножаем первую строку полученного определителя на и вычитаем ее из второй строки определителя. Получаем
1
0
=
.
(4.17)
Таким образом, определитель равен произведению ведущих элементов схемы Гаусса. Все ведущие элементы должны быть отличны от нуля. Если для какого-нибудь шага ведущий элемент равен нулю или близок к нулю, следует соответствующим образом изменить порядок строк и столбцов матрицы.
П ример. Вычислить при помощи метода Гаусса определитель
2 1 4
3 2 -1
-4 1 3 .
Вычисления удобно оформить в виде табл. 4.3.
Таблица 4.3
Превый столбец |
Второй столбец |
Третий столбец |
|
|
1 |
4 |
7 |
3 |
2 |
-1 |
4 |
-4 |
1 |
3 |
0 |
1 |
1/2 |
2 |
7/2 |
|
1/2 |
-7 |
-13/2 |
|
3 |
11 |
14 |
|
1 |
-14 |
-13 |
|
|
53 |
53 |
|
Последний столбец страницы служит для контроля вычислений.
4.3. Вычисление обратной матрицы
при помощи метода Гаусса
Пусть дана неособенная матрица А:
…
А =
…
… … … …
…
Требуется найти матрицу, обратную данной. Обозначим
…
Х=А-1
…
… … … …
…
По определению
обратной матрицы
,
где Е – единичная матица, т.е.
.
(4.18)
Пусть
К-й столбец матрицы
,
а
к-й
столбец единичной матрицы.
Тогда равенство (4.18) можно заменить следующими равенствами:
.
(4.19)
Таким образом,
определение элементов обратной матрицы
эквивалентно решению n
систем линейных алгебраических уравнений
вида
.
Все они имеют одну и ту же матрицу
коэффициентов А, а отличаются только
свободными членами. Поэтому решение
этих n систем можно
объединить в одну схему, рассматривая
одновременно n столбцов
свободных членов.
Пример. Найти обратную матрицу для матрицы
1 2 -1
0 1 3
2 -1 -2 .
Вычисление оформить в виде табл. 4.4.
Таблица 4.4
|
|
|
|
|
|
|
1 |
2 |
-1 |
1 |
0 |
0 |
3 |
0 |
1 |
3 |
0 |
1 |
0 |
5 |
2 |
-1 |
-2 |
0 |
0 |
1 |
0 |
|
1 |
3 |
0 |
1 |
0 |
5 |
|
-5 |
0 |
-2 |
0 |
1 |
-6 |
|
|
15 |
-2 |
5 |
1 |
19 |
|
|
1 |
-2/15 |
5/15 |
1/15 |
19/15 |
|
|
1 |
|
|
|
19/15 |
|
1 |
|
=6/15 |
=0 |
|
18/15 |
1 |
|
|
=1/15 |
=5/15 |
|
28/15 |
Ответ:
1 5 7
6 0 -3
-2 5 1 .