Материал: лекция 2 ММОСУ

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Лекция 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 24. Для исключения x1 в

третьем уравнения из него вычтем первое, умноженное на l31 14. В результате после первого шага исключения получим

 

 

9

 

 

 

 

 

 

4

2

6

 

 

 

 

1

 

 

 

 

 

 

 

0

 

3

3

.

2

 

 

 

 

 

 

 

 

 

 

0

1

 

5

 

5

 

 

 

 

 

 

 

 

 

 

4

 

2

 

2

 

 

 

 

 

На втором шаге алгоритма исключим x2 в третьем уравнении, для чего из третьего уравнения вычтем первое, умноженное на l32 12. После второго шага исключения получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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 a21a11 ;

из n-ого уравнения первого, умноженного на ln1 an1a11 . После первого шага система уравнений преобразуется к виду

 

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 - номер столбца.