
- •Системы линейных уравнений Элементы теории линейных векторных пространств Понятие вектора
- •Операции над векторами Сложение векторов:
- •Умножение вектора на константу
- •Скалярное произведение двух векторов
- •Элементы матричной алгебры Основные определения
- •Виды матриц
- •Операции над матрицами
- •Числовые характеристики матрицы Норма матрицы
- •Определитель матрицы
- •Ранг матрицы, минор, алгебраические дополнения
- •Методы обращения матрицы Метод Крамера
- •Метод Гаусса-Жордана
- •Системы линейных алгебраических уравнений (слау)
- •Классификация слау
- •Условия совместности и определенности систем
- •Обусловленность слау
- •О методах решения линейных систем
- •Метод Гаусса
- •Блок-схема метода Гаусса
- •Итерационные методы решения слау
- •Метод простой итерации Приведение системы к итерационной форме
- •Оценка сходимости решения
- •Выполнение итерации
- •Проверка условия окончания решения
- •Метод Гаусса-Зейделя
- •Уточнение решения
Числовые характеристики матрицы Норма матрицы
Норма матрицы является некоторой обобщенной оценкой значений элементов матрицы. Для вычисления нормы матрицы можно использовать следующие выражения:
.
,
i
= 1, 2, … , m .
,
j = 1, 2, … , n .
Определитель матрицы
Определитель существует только для квадратных матриц. Выражения для расчета определителя в общем виде записывается как
и представляет собой сумму всевозможных произведений элементов матрицы, расположенных в разных строках и разных столбцах ее. В этих произведениях элементы обычно располагаются в порядке возрастания первых индексов. Совокупность вторых индексов образует некоторую перестановку чисел от 1 до n - j1, j2, … , jn .Так как количество таких перестановок равно n!, то сумма в выражении для определителя имеет n! слагаемых.
Степень s определяется четностью перестановки: s = 1, если перестановка нечетная и s = 2, если перестановка четная. Четность перестановки определяется числом инверсий в ней: четное число инверсий дает четную перестановку, нечетное - нечетную перестановку. Инверсией является каждое нарушение порядка следования чисел в перестановке по возрастанию.
Из этого выражения следует, что определитель равен сумме n! слагаемых, каждое из которых является произведением n элементов. Поэтому для вычисления определителя порядка n (без использования специальных приемов) требуется (n-1)n! умножений и n!-1 сложений, т.е. общее число арифметических операций равно
N = n * n! - 1 » n * n! .
Можно оценить значения N в зависимости от порядка n определителя:
n |
3 |
10 |
20 |
N |
17 |
3.6*107 |
5*1019 |
Можно подсчитать время вычисления таких определителей на ЭВМ с заданным быстродействием. Если для определенности принять среднее быстродействие равным 1 000 000 операций в секунду, то для вычисления определителя 10 порядка потребуется 36 секунд, а при n = 20 - около 1.4*1010 часов, т.е. около 5.8*108 суток » 1.6 млн. лет .
Вычисление определителя
Легко вычисляются определители невысоких порядков и некоторые специальные типы определителей. Например,
a11a22a33+a12a23a31+a21a32a13-a31a22a13-a21a12a33-a32a23a11
Свойства определителя
1. Если в матрице есть нулевая строка или нулевой столбец, то ее определитель равен нулю.
2. Если матрица содержит две одинаковые строки или два одинаковых столбца, то ее определитель равен нулю.
3. Если в матрице поменять местами две строки или два столбца, то ее определитель изменит знак на противоположный.
4. Если все элементы какой-либо строки или какого-либо столбца матрицы увеличить в a раз, то определитель матрицы также увеличится в a раз.
5. Если к какой-либо из строк матрицы прибавить другую строку, умноженную на константу, то определитель матрицы не изменится.
6. Определитель треугольной матрицы равен произведению ее диагональных элементов:
.
Ранг матрицы, минор, алгебраические дополнения
Ранг матрицы показывает число линейно независимых строк и столбцов матрицы. Ранг существует как для квадратных, так и для прямоугольных матриц. Для вычисления ранга используется понятие минора.
Минором k-го порядка для любой прямоугольной матрицы А размером m*n (k£m, k£n) называется определитель квадратной подматрицы, полученной путем выделения из исходной матрицы А элементов, стоящих на пересечении любых k строк и k столбцов ее.
Ранг матрицы численно равен порядку наивысшего отличного от нуля минора.
Примеры: Рассчитать ранги следующих матриц:
1)
,ra
= 3, т.к.
det(A)
¹
0 .
2)
,ra
= 2,
т.к.
3)
,ra
=
1, т.к.
det(A) = 0 и
все миноры второго порядка также равны
нулю.
Алгебраическое дополнение Ai j = (-1)i+j * Mi j для элемента ai j с точностью до знака равно минору порядка n-1, полученному из определителя исходной матрицы путем вычеркивания i-ой строки и j-го столбца.
Определитель матрицы равен сумме всех произведений элементов какой-либо строки или столбца на соответствующие им алгебраические дополнения, т.е.
где первая формула дает разложение определителя d по элементам i-й строки, а вторая формула - по элементам j-го столбца.
Пример. Вычислить определитель
разложив его по элементам второй строки.
Решение. Имеем d = a21A21 + a22A22 + a23A23, откуда
Этот же ответ можно было получить используя свойства определителя, т.к. его вторая строка равна сумме третьей строки и первой, умноженной на -1.