Лекция 2. Системы линейных алгебраических уравнений
2.1 Решение систем линейных алгебраических уравнений
2.1.1 Определения и основные свойства матриц
Для начала напомним некоторые сведения из линейной алгебры. Рассмотрим систему из n линейных алгебраических уравнений (СЛАУ) с
n неизвестными x1, ,xn :
a11x1 a12x2 a1nxn b1, |
|
|
|
|
|
||||||||||||||
|
|
|
|
a22x2 a2nxn b2, |
|
|
|
|
|
||||||||||
a21x2 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
n1 |
a |
n2 |
x |
2 |
a |
nn |
x |
n |
b . |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|||||
Запишем СЛАУ в матричной форме, для этого введем векторы |
|||||||||||||||||||
|
x |
|
|
|
b |
|
|
|
|
|
a |
|
a |
|
|||||
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
11 |
|
1n |
|||
x |
,b |
|
и матрицу A |
|
|
|
. |
||||||||||||
|
x |
n |
|
|
|
b |
|
|
|
|
|
a |
n1 |
a |
|
|
|||
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
nn |
|||||
Тогда СЛАУ примет вид |
|
|
|
|
|
||||||||||||||
Ax b. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Определение: A 1 |
– обратная матрица для квадратной матрицы А, если |
||||||||||||||||||
AA 1 I , где I |
– единичная матрица. |
|
|
|
|
|
|||||||||||||
Определение: А — невырожденная, если для А существует A 1. Теорема: следующие утверждения эквивалентны:
1.А — невырождена;
2.det A 0;
3.линейная однородная система Ax 0 имеет единственное решение
x0;
4.для любого b система Ax b имеет единственное решение;
5.столбцы (строки) матрицы А линейно независимы, то есть для любой
комбинации чисел 1, 2, n , не все из которых равны нулю, линейная комбинация
|
a11 |
|
|
|
a12 |
|
|
|
a1n |
|
|
|
0 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a21 |
|
|
|
a22 |
|
|
|
a2n |
|
|
|
0 |
|
|
|||
|
|
|
2 |
|
|
|
n |
|
|
|
|
|
|
|
||||
1 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
an1 |
|
|
|
an2 |
|
|
|
ann |
|
|
|
|
|
||||
Определение: ранг матрицы – |
|
rank(A) – число линейно независимых |
||||||||||||||||
столбцов (строк). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Следствие: А — невырожденная, если и только если А — полного ранга |
||||||||||||||||||
(rank(A)=n). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Определение: комплексное или вещественное число и вектор |
x 0, |
|||||||||||||||||
x Rn |
|
называются |
собственным |
значением и собственным вектором |
||||||||||||||
матрицы А, если они удовлетворяют алгебраическому уравнению
Ax x
Собственный вектор — это такой вектор, который, будучи умножен на матрицу А, изменяет лишь свою длину.
Из определения собственного значения следует:
Ax x 0 A I x 0.
Так как x — собственный вектор и по определению x 0, то в соответствии с теоремой A I — вырожденная матрица, а det A I 0.
Левая часть последнего выражения — полином степени n относительно комплексной переменной , называемый характеристическим полиномом матрицы, следовательно собственные значения матрицы есть корни характеристического полинома матрицы A, их ровно n с учетом кратности.
Определение: совокупность всех собственных значений i in 1 матрицы
А называется спектром матрицы А, а величина
A max i
1 i n
называется спектральным радиусом. Определение: нормой матрицы А называется
A |
|
|
|
sup |
|
|
Ax |
|
|
AT A1 2. |
||
|
|
|
|
|
|
|||||||
|
|
|||||||||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
x |
|
|
|
|||
|
|
|
|
|
|
|
||||||
|
|
|
|
x Rn |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
Решение СЛАУ методом исключения
Одной из наиболее распространенных задач является задача решения СЛАУ
Ax b,
для заданных квадратной матрицы А и вектора b.
Если А — невырожденная матрица, то для нее существует обратная A 1.
Умножая СЛАУ на A 1 слева, получим решение: x A 1b.
Полученный метод решения СЛАУ неэффективен в силу большого времени вычислений.
Рассмотрим более эффективный алгоритм решения СЛАУ – метод исключения Гаусса. Этот алгоритм с некоторыми модификациями используют для решения СЛАУ до двухсотого порядка.
Вначале продемонстрируем работу алгоритма на примере следующей СЛАУ:
4x1 9x2 2x3 6, |
|
|
|
|
|||||||
|
|
|
4x3 6, |
|
|
|
|
||||
2x1 4x2 |
|
|
|
|
|||||||
x 2x |
2 |
2x 1, |
|
|
|
||||||
|
1 |
|
3 |
|
|
|
|
|
|||
которую запишем в матричной форме |
|
||||||||||
4 |
9 |
|
2 x1 |
|
|
6 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
4 |
|
4 |
x2 |
|
|
6 |
. |
|
|
|
1 2 |
|
2 |
x |
3 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Введем расширенную матрицу системы: |
|
||||||||||
|
4 |
9 |
|
2 |
6 |
|
|
|
|
|
|
|
2 |
4 |
|
4 |
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
1 |
2 |
|
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
На |
первом |
шаге |
алгоритма исключим x1 из второго |
и третьего |
|||||||
уравнений системы. |
Для исключения x1 во втором уравнении, |
из второго |
|||||||||
уравнения вычтем первое, умноженное на l21 2
4. Для исключения x1 в
третьем уравнения из него вычтем первое, умноженное на l31 1
4. В результате после первого шага исключения получим
|
|
9 |
|
|
|
|
|
||
|
4 |
2 |
6 |
|
|||||
|
|
|
1 |
|
|
|
|
|
|
|
0 |
|
3 |
3 |
. |
||||
2 |
|
||||||||
|
|
|
|
|
|
|
|
||
|
0 |
1 |
|
5 |
|
5 |
|
||
|
|
|
|
|
|
|
|
|
|
4 |
|
2 |
|
2 |
|||||
|
|
|
|
|
|||||
На втором шаге алгоритма исключим x2 в третьем уравнении, для чего из третьего уравнения вычтем первое, умноженное на l32 1
2. После второго шага исключения получим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
9 |
2 |
6 |
|
|
|
0 |
|
0 |
0 |
|
|||||
|
|
|
|
||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
||||||||
|
0 |
|
3 3 |
|
, |
L |
|
1 |
|
0 |
0 |
. |
|||||
2 |
|
2 |
|
||||||||||||||
|
0 |
|
4 |
4 |
|
|
|
|
|
1 |
|
|
|||||
|
0 |
|
|
|
|
|
|
1 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
4 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|||
За два шага закончилась первая часть алгоритма — приведение системы к треугольной форме. Решение полученной треугольной системы осуществляется снизу вверх
x3 44 1;
x2 2 3 3x3 0;
x1 14 6 2x3 9x2 1.
Перейдем к построению алгоритма исключения для СЛАУ общего вида
a11x1 a12x2 a1nxn b1, |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
a21x2 a22x2 a2nxn b2, |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
n1 |
a |
n2 |
x |
2 |
a |
nn |
x |
n |
b . |
|
|
|
|
|
n |
|||||
Шаг 1: исключим x1 из уравнений с номерами 2,…, n путем вычитания:
из второго уравнения первого, умноженного на l21 a21
a11 ;
…
из n-ого уравнения первого, умноженного на ln1 an1
a11 . После первого шага система уравнений преобразуется к виду
|
a11x1 a12x2 a1nxn b1, |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
a 1 x |
|
|
a |
1 x |
|
|
|
b 1 , |
|
|
|
|
|
||||||
|
|
|
|
2 |
n |
|
|
|
|
|
|
||||||||||||
|
|
|
22 |
|
|
|
|
2n |
|
2 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a 1 x |
2 |
|
a 1 x |
n |
b 1 . |
|
|
|
|
|
|||||||||
|
|
|
|
n2 |
|
|
|
|
|
nn |
|
n |
|
|
|
|
|
|
|||||
где |
a 1 |
a |
ij |
l |
a |
|
|
,b 1 |
b l |
b ,l |
i1 |
a |
a ,i, j 2, ,n |
||||||||||
ij |
|
|
i1 |
1j |
i |
|
|
i |
|
|
i1 1 |
|
|
|
i1 |
11 |
|||||||
|
Шаг 2: исключим |
x2 |
из уравнений с номерами 3,…, n путем вычитания: |
||||||||||||||||||||
|
Из третьего уравнения второго, умноженного на l32 a32 a22 |
||||||||||||||||||||||
|
…. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из n-ого уравнения второго‚ умноженного на ln2 an2 a22 . |
||||||||||||||||||||||
|
После второго шага получим |
|
|
|
|
|
|
|
|||||||||||||||
|
a11x1 a12x2 a13x3 a1nxn b1, |
|
|||||||||||||||||||||
|
|
|
|
a 1 x |
|
|
a |
1 x a 1 x |
|
|
b 1 , |
|
|||||||||||
|
|
|
|
2 |
n |
|
|||||||||||||||||
|
|
|
22 |
|
|
23 |
3 |
|
|
|
|
2n |
|
2 |
|
||||||||
|
|
|
|
|
|
|
|
|
a |
2 x a 2 x |
n |
b |
2 |
, |
|||||||||
|
|
|
|
|
|
|
|
|
33 |
3 |
|
|
|
|
3n |
|
|
3 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xn |
|
|
|||||||
|
|
|
|
|
|
|
|
|
an3 x3 |
ann |
bn |
. |
|||||||||||
После n-1 шага исключения исходная система придет к треугольной
форме
a11 |
a12 |
|
|
0 |
a221 |
|
||
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
||
где
a1n |
x1 |
|
|
|
|
b1 |
|
||
1 |
|
|
|
|
|
|
|
1 |
|
a2n |
x2 |
|
|
|
b2 |
|
|||
|
|
|
|
|
|
|
|
||
|
|
|
|||||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
||||
a n 1 x |
n |
|
|
b |
n 1 |
||||
nn |
|
|
|
|
|
n |
|
||
a k 1
lik aikk 1 ,bik bik 1 likbkk 1 ,aijk aijk 1 likaijk 1 , kk
k=1,…n–1– номер шага исключения, i=k+1, …, n – номер строки, j=k+1, …, n - номер столбца.