- •ЛЕКЦИЯ 1
- •СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ (СЛАУ)
- •РЕШЕНИЕ СЛАУ
- •РЕШЕНИЕ СЛАУ
- •РЕШЕНИЕ СЛАУ МЕТОДОМ ОБРАТНОЙ МАТРИЦЫ
- •РЕШЕНИЕ СЛАУ МЕТОДОМ КРАМЕРА
- •Матрицы специального вида
- •МЕТОД ГАУССА
- •МЕТОД ГАУССА
- •МЕТОД ГАУССА
- •МЕТОД ГАУССА
- •МЕТОД ГАУССА
- •МЕТОД ГАУССА
- •Lsolve (A,b)
- •Симметричные матрицы
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •МЕТОД ПРОГОНКИ
- •ИТЕРАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ СЛАУ
- •Одношаговые итерационные методы
- •МЕТОД ЯКОБИ
- •МЕТОД РЕЛАКСАЦИИ
- •Вычислительный блок Given/Find
- •Вычислительный блок Given/Find
МЕТОД ГАУССА
КОЛИЧЕСТВО ОПЕРАЦИЙ m3
КОЛИЧЕСТВО ОПЕРАЦИЙ НА ОДНУ НЕИЗВЕСТНУЮ
m2
ПОНЯТИЕ ЭКОНОМИЧНОГО МЕТОДА:
КОЛИЧЕСТВО ОПЕРАЦИЙ НА ОДНУ НЕИЗВЕСТНУЮ НЕ ЗАВИСИТ ОТ ЧИСЛА НЕИЗВЕСТНЫХ
МЕТОД ГАУССА НЕ ЯВЛЯЕТСЯ ЭКОНОМИЧНЫМ, ПОЭТОМУ НЕ ИСПОЛЬЗУЕТСЯ ПРИ m>40
НЕТ ТОЧНЫХ ЭКОНОМИЧНЫХ МЕТОДОВ РЕШЕНИЯ |
|
СЛАУ ОБЩЕГО ВИДА |
13 |
Функция MathCAD для решения СЛАУ
Реализует метод Гаусса
Lsolve (A,b)
|
|
1 |
5 |
2 |
|
|
1 |
|
|
|
0.7 |
12 |
|
|
|
2.9 |
|
A |
|
5 |
b |
|
|
|||
|
|
3 |
0 |
4 |
|
3.1 |
|
0.186 lsolve(A b) 0.129
0.915
14
МЕТОД ГАУССА ЖОРДАНО
Исходная система аналогичными преобразованиями сводится к диагональной матрице
0 |
0 |
|
ПРИМЕР 1
|
x1 x2 x3 2 |
|
|
1 |
1 |
1 |
|
2 |
||||||||
2x x |
2 |
x |
3 |
A |
|
2 |
1 |
1 |
, |
b |
|
3 |
|
|||
|
1 |
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
x x |
2 |
x 6 |
|
|
|
1 |
1 |
1 |
|
|
|
6 |
|
||
|
1 |
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
1 |
2 |
|
1 |
1 |
1 |
2 |
|
1 |
1 |
1 |
2 |
|
|
||||||||||||
2 |
1 |
1 |
3 |
|
0 |
3 |
1 |
7 |
|
0 |
1 |
1/ 3 7 / 3 |
|
1 |
1 |
1 |
6 |
|
0 |
0 |
2 |
4 |
|
0 |
0 |
2 |
4 |
1 |
0 |
2 / 3 |
1/ 3 |
|
1 |
0 |
0 |
1 |
|
|
|
1 |
|
|
|
|
|||||||||||
|
|
||||||||||||
0 |
1 |
1/ 3 7 / 3 |
|
0 |
1 |
0 |
3 |
|
|
|
3 |
|
|
|
|
|
x |
|
|||||||||
0 |
0 |
1 |
2 |
|
0 |
0 |
1 |
2 |
|
|
|
2 |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Метод Гаусса-Жордано
•Позволяет найти обратную матрицу;
•Добавить в расширенную матрицу единичную матрицу;
•Преобразования строк, приводящие исходную матрицук диагональному виду, будут переводить единичную матрицу в обратную, а вектор правых частей – в решение
|
A |
|
|
E |
|
b |
|
|
|
|
|
|
|
1 |
1 |
1 |
1 |
0 |
0 |
2 |
|
|
|
|
|
|
|
2 |
1 |
1 |
0 |
1 |
0 |
3 |
|
|
|
|
|
|
|
1 |
1 |
1 |
0 |
0 |
1 |
6 |
|
|
|
|
|
|
|
2. Прибавляем ко второй строке первую, умноженную на 2, из третьей вычитаем первую
1 |
1 |
1 |
1 |
0 |
0 |
2 |
|
|
|
|
|
|
|
0 |
3 |
1 |
2 |
1 |
0 |
7 |
|
|
|
|
|
|
|
0 |
0 |
2 |
1 |
0 |
1 |
4 |
|
|
|
|
|
|
|
3. Делим вторую строку на 3 и вычитаем ее из первой
1 |
0 |
2/3 |
1/3 |
1/3 |
0 |
1/3 |
|
|
|
|
|
|
|
0 |
1 |
1/3 |
2/3 |
1/3 |
0 |
7/3 |
|
|
|
|
|
|
|
0 |
0 |
2 |
1 |
0 |
1 |
4 |
|
|
|
|
|
|
|
4. Делим последнюю строку на 2 и прибавляем ее с коэффициентом 1/3 ко второй строке и с коэффициентом 2/3 – к первой
1 |
0 |
0 |
0 |
1/3 |
1/3 |
1 |
0 |
1 |
0 |
1/2 |
1/3 |
1/6 |
3 |
|
|
|
|
|
|
|
0 |
0 |
1 |
1/2 |
0 |
1/2 |
2 |
17
Метод Гаусса-Жордано
• Позволяет исследовать и решать недоопределенные, переопределенные, вырожденные СЛАУ
|
1 x1 0 x2 1 x3 2 |
Пример 2 |
1 x1 1 x2 0 x3 3 |
|
1 x1 1 x2 1 x3 2 |
|
2 x1 1 x2 1 x3 1 |
Число неизвестных n= 3 |
|
Число уравнений m=4 |
|
Если все уравнения независимы, то задача неразрешима.
Исследуем систему с помощью метода Гаусса-Жордано
Пример 2
Составляем таблицу:
•Выберем разрешающий элемент a11
•Выполним преобразования Гаусса –Жордано относительно этого разрешающего элемента.
•Все элементы столбца надо превратить в 0,
прибавляя или отнимая разрешающий элемент
Пример 2
Шаг 1:
|
0 |
|
|
1 |
|
|
5 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
0 |
|
|
|
|
|
|
|
|
||
|
|
|
2 |
|
4 |
|||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
1 |
|
|
|
||||
0 |
|
|
5 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
Шаг 2:
|
0 |
|
|
|
3 |
|
|
9 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|||||
|
|
0 |
|
|
|
|
||||
|
|
|
|
0 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Стр#2 - Cтр#1
Стр#3 - Cтр#1
Стр#4 -2*Cтр#1
Стр#3 + Cтр#2
Стр4 - Cтр#2
Получена нулевая строка – последнее уравнение лишнее, можно его отбросить!
Пример 2
Шаг 3
0 1
0 2
1 3
Стр#1 + Стр#3
Стр#2 - Стр#3
Стр#3 /3
Запишем СЛАУ в векторной форме:
|
1 |
|
|
|
0 |
|
|
|
0 |
|
|
|
1 |
|
Решение |
||
|
0 |
|
|
|
1 |
|
x 3 |
|
0 |
|
|
|
2 |
|
|
|
1 |
x 1 |
|
|
x 2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
|
1 |
|
3 |
X |
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
базисные вектора |
|
|
|
|
|
|
|
|
|
3 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Пример 3
Решить СЛАУ: 2 x1 0 x2 1 x3 1 x4 1
3 x1 1 x2 1 x3 0 x4 2 5 x1 1 x2 2 x3 1 x4 3 1 x1 1 x2 0 x3 1 x4 1
Число неизвестных n= 4 Число уравнений m=4
Решим систему с помощью метода Гаусса-Жордана
Пример 3
Исходная таблица
Шаг 1: переменная x2 – в базис
Шаг 2: переменная x4 – в базис
Пример 3
Шаг 2: переменная x4 – в базис
СЛАУ в векторной форме:
Две строки нулевые, отбрасываем их
x1 2 x2 0 x3 1 x4 1 13 1 1 0 2
базисные вектора
Пример 3
x1 2 x2 0 x3 1 x4 1 13 1 1 0 2
базисные вектора
Выражаем базисные переменные через свободные
x4 1 2 x1 x3 x2 2 3 x1 x3
Получили общее решение, зависящее от двух параметров – x1. x3, изменяя которое получим множество всех решений СЛАУ
Метод Гаусса Жордано
•Если в результате преобразований получены одинаковые строки, то одну из них можно отбросить, т.е. сократить число уравнений
•Если число уравнений меньше числа неизвестных, то система имеет бесконечное множество решений
•Общее решение СЛАУ зависит от параметров (свободных переменных), которые могут принимать произвольные значения
•Если в результате преобразований все коэффициенты в левой части занулились, а правые части остались ненулевыми, система несовместна (не имеет решений)