- •Глава 4. Численные методы алгебры. Решение систем лИнейных уравнений
- •4.1. Линейные уравнения. Теоретическое и практическое решения линейных уравнений с одним неизвестным
- •4.2. Системы линейных уравнений. Основные понятия
- •4.3. Определители (детерминанты) квадратных матриц
- •Свойства определителей
- •Расчет определителей
- •4.4. Необходимо и достаточное условие существования решения системы линейных уравнений. Методы решения
- •4.5. Прямые методы решения систем линейных уравнений
- •1. Метод с использованием обратной матрицы.
- •[Править] Метод Гаусса—Жордана
- •[Править] с помощью матрицы алгебраических дополнений
- •[Править] Использование lu/lup-разложения
- •[Править] Примеры [править] Матрица 2х2
- •2.1. Погрешности вычислений на эвм
- •1. Метод Гаусса
- •Будем рассматривать столбцы матрицы X как векторы
- •6. Итерационные методы
4.2. Системы линейных уравнений. Основные понятия
Системой m линейных уравнений с n неизвестными (или линейной системой) в линейной алгебре называют систему m линейных уравнений канонического вида с n неизвестными вида
|
(4.6) |
В (4.6) x1, x2, …, xn — неизвестные, a11, a12, …, amn — постоянные линейные коэффициенты системы (задающие матрицу системы), b1, b2, … bm — свободные члены (задающие вектор свободных коэффициентов). Индексы коэффициентов (aij) системы обозначают номера уравнения (i) и неизвестного (j), при котором стоит этот коэффициент, соответственно.
Систему (4.6) называются однородной, если все её свободные члены равны нулю (b1 = b2 = … = bm = 0), иначе — неоднородной.
Система (4.6) называется квадратной, если число m уравнений равно числу n неизвестных. Число n= m называют порядком системы.
Решение системы (4.6) — совокупность n чисел c1, c2, …, cn, таких что подстановка каждого ci вместо xi в систему (4.6) обращает все её уравнения в тождества. Найти решение системы означает выяснить, существует ли у системы решения, и если они существуют, найти их.
Система (4.6) называется совместной, если она имеет хотя бы одно решение, и несовместной, если у неё нет ни одного решения. Совместная система вида (4.6) может иметь одно или более решений.
Решения c1(1), c2(1), …, cn(1) и c1(2), c2(2), …, cn(2) совместной системы вида (4.6) называются различными, если нарушается хотя бы одно из равенств:
c1(1) = c1(2), c2(1) = c2(2), …, cn(1) = cn(2). |
Совместная система вида (4.6) называется определённой, если она имеет единственное решение; если же у неё есть хотя бы два различных решения, то она называется неопределённой. Если уравнений больше, чем неизвестных, она называется переопределённой.
В технических приложениях, как правило, рассматриваются системы с квадратными матрицами (m = n). Поэтому в дальнейшем будем рассматривать только этот случай.
Множество коэффициентов (aij) квадратной системы образуют ее матрицу А, а множество коэффициентов (bi) – вектор В свободных коэффициентов системы. Обозначив вектор (x1, x2, …, xn) неизвестных черезХ, в векторном виде систему представляют в виде:
AХ =В. (4.7)
Вопросы для проверки знаний.
1. Какие системы линейных уравнений называют однородными, а какие – неоднородными ?
2. Какие системы линейных уравнений называют квадратными ?
3. Что называют решением системы линейных уравнений ?
4. Какие системы линейных уравнений называют определённой, неопределённой и переопределёнными ?
5. Как системы линейных уравнений представляют в векторной форме ?
4.3. Определители (детерминанты) квадратных матриц
Главной характеристикой квадратной матрицы является определитель (детерминант) - численная характеристика матрицы, которая может быть определена двумя способами - 1) либо рекурсивно по порядку матрицы n либо 2) при помощи перестановок. Обозначается определитель матрицы А как (А) либо det(A).
Рекурсивное определение определителя матрицы А по ее порядку n.
n=1. Для матрицы первого порядка детерминантом является сам единственный элемент этой матрицы:
(4.8а)
n2. Для матриц второго порядка и выше детерминант задаётся рекурсивно формулой, которая называется разложением определителя по первой строке:
, (4.8б)
где — дополнительный минор к элементу a1j. - определитель матрицы порядка (n-1), которая получается из А вычеркиванием в ней первой строки и столбца с номером j .
Определение определителя при помощи перестановок. Для квадратной матрицы А порядка n справедлива формула:
, (4.9)
где α1,α2,...,αn — перестановка чисел от 1 до n, N(α1,α2,...,αn) — число инверсий (нарушений порядка) в перестановке, суммирование идёт по всем возможным перестановкам порядка n. Поскольку общее их число равно n!, то в сумму (4.9) войдут n! слагаемых, которые называют членами определителя.
Используя рекурсивное определение, найдем расчетные формулы для определителей матриц порядков n=2 и n=3.
(4.10)
(4.11)
Примеры 1. Расчет определителей: