Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Линейная алгебра

..pdf
Скачиваний:
32
Добавлен:
12.11.2023
Размер:
13.06 Mб
Скачать

А = L U

(1.16)

невырожденной левой нижней треугольной матрицы L с единицами на главной диагонали и правой верхней треугольной матрицы U. Раз­ ложение (1.16) матрицы А называют ее LU-разложением. Оно нахо­ дит широкое применение в вычислительной практике. Такое разложе­ ние можно получить, например, с помощью метода Гаусса. Поясним это на примере.

П ример 1.

Построить LiZ-разложение матрицы

 

/

2

i

i \

 

А = I

4

1

0

 

V

- 2

2

1

Реш ение.

Первую строку матрицы А, умноженную на —2,

прибавим ко второй ее строке, затем первую строку, умноженную на 1, прибавим к третьей строке. В результате получим матрицу

2

1

1

0

-1

-2

0

3

2

Прибавим к третьей строке этой матрицы вторую ее строку, умно­ женную на 3. Тогда получим матрицу

 

2

1

1

U =

0

-1

-2

 

0

0

- 4

Проведенные элементарные операции над строками равносильны умно­ жению матрицы А слева последовательно на элементарные матрицы

/

1 0 0 \

/ 1 0 0 \

/ 1 0 0 \

Е21 = 1 - 2 1 0 ]

, £731 = [ 0

1

0

, Ез2 =

0 1 0

\ 0 0 1 /

\ 1

0

1 /

V 0 3 1 /

В итоге имеем Е32Е31Е21А = U, т.е. SA = U, гдеположено

 

/

i

0

0 \

 

S = Е32Е31Е21 = I —2 1 0 I

 

 

V -5

3

1 /

 

Умножая обе части равенства SA = U слева на S \ получим Ы7разложение

A = S~1U = LU =

где положено L = S'-1 .

Для вычисления неизвестных элементов матриц L и U разложе­ ния (1.16) можно также провести умножение матриц в правой части этого разложения, сравнить элементы полученной матрицы с соот­ ветствующими элементами матрицы А и решить полученную систему относительно неизвестных элементов матриц LnU. Решение в общем виде дает формулы для вычисления нужных элементов (см. [6], [7]).

Поясним такой подход на примере.

 

 

П ример 2.

Построить LU-разложение матрицы

 

(

I

2

—1 \

 

А - [

2

1

- 1

3 /

 

V

1

- 7

 

Решение.

Запишем подробно L/7-разложение матрицы А пока с

неизвестными элементами матриц L и U

/ 1

2 —1

\

/ 1

0

0

\

/

tin

tii2

tils\

2

1 —1

 

I = [

hi

1

 

о

I

0I

1*22

U23 )

\ 1

—7 3

/

\

/31

/32

1

/

\ 0

0

1*зз /

Проведя умножение матриц в правой части этого равенства, получим

(

u n

t*i2

^22

г*1з

 

 

hlUll

hl^l2+

/21^13 + ^*23

 

 

/з1**11

1^*12+

/з2**22

1^*13 + /32^23

+ **33

Сравнивая соответствующие элементы этих матриц, придем стеме уравнении

*11 = 1,

 

г*12 = 2,

 

 

«13 = -1 ,

/21**п =

2,

/211*12 +

t*22 = 1,

 

bl«13 + «23 = —1)

/31^*11 =

1,

/311*12 + /32^*22 =

—7,

1з1«13 +

/32«23 + «33 = 3.

Решив эту систему, найдем

 

 

 

 

 

t*n = l,

г112 =

2,

«13 =

-1 ,

 

 

hi == 2,

1*22 =

—3,

«23 =

1,

 

 

/31 = 1,

/32 =

3,

«33 =

1.

Следовательно,

 

 

 

 

 

 

 

/ 1

2

—1 \

 

/ 1

0

0 \

/ 1

2

—1 \

2

1 - 1

=

2 1 0 0 - 3

 

1 ) -

\ 1

—7

3 /

 

\ i

Э

i /

\ 0

0

1 /

- искомое £{/-раэложение матрицы А .

Выделив главную диагональ матрицы U в разложении (1.16) в диа­ гональную матрицу, матрицу U можно представить в виде произведе­ ния U = D U\ с правой верхней треугольной матрицей U\ с единицами на главной диагонали и диагональной матрицей D. Тогда разложение (1.16) превратится в разложение

A = L -D -Ih .

(1.17)

Если матрица А симметрическая, то разложение (1.17) превра­

тится в разложение

 

A = Ui ■D U\.

(1.18)

Если, кроме того, в матрице D все диагональные элементы rfi, с?2э

..., dn положительные, то можно положить V = D1! 2 Ui, где D 1! 2 - диагональная матрица с элементами \/3i, \/^2, •••, V^n по диагонали. Тогда разложение (1.18) превратится в разложение

А = V T V.

(1.19)

Разложения (1.17), (1.18), (1.19), как и разложение (1.16), широко применяются в вычислительной практике.

1.7.А лгебраические операции, группы , кольца, поля

Влюбой области деятельности приходиться рассматривать раз­ личные совокупности объектов, объединенные некоторым общим при­ знаком. Такие совокупности объектов в математике принято назы­ вать м нож ествам и, а сами объекты — элементами м нож еств. Множества будем обозначать прописными латинскими буквами А,

В, С, ..., Х } У,

..., а их элементы — малыми латинскими буквами

а, 6, с, ..., х , у,

При этом будем писать х Е А, если элемент х

принадлежит множеству А, и х £ А, если элемент х не принадлежит множеству А.

3-1307

Множество, которое не содержит ни одного элемента, будем на­ зывать пусты м м нож еством и обозначать символом 0

Пусть дано некоторое множество А, содержащее хотя бы один элемент. Будем говорить, что в множестве А определена алгебраиче­ ская операция, если указан закон, по которому любой паре элементов а и 6, взятых из этого множества в определенном порядке, ставится во взаимно однозначное соответствие некоторый элемент с, также принажлежащий этому множеству. Бели это операция сложения, то элемент с будем называть суммой элементов а и 6 и обозначать символом а + 6; если это операция умнож ения, то элемент с будем называть произведением элементов а и 6 и обозначать символом ab. В общем случае алгебраическую операцию будем обозначать символом *.

Алгебраическая операция * называется ком м утативной, если ре­ зультат ее применения не зависит от порядка выбора элементов, т.е. для любых элементов а и Ьиз рассматриваемого множества выполня­ ется равенство а*6 = Ь*а. Например, операции сложения и умножения чисел являются коммутативными операциями, а вычитание и деление - некоммутативными операциями.

Поскольку алгебраическая операция определена для пар элемен­ тов, то для трех элементов а, 6, с возможным будет рассмотрение либо элемента (а * Ь) * с, либо элемента а * (6 * с). В общем случае эти элементы могут оказаться различными.

Алгебраическую операцию * называют ассоциативной, если для любых трех элементов а, 6, с рассматриваемого множества выполня­ ется равенство (а*6)*с = а * (6*с). В этом случае результат операции *, примененный к элементам а, 6, с записывают в виде а*Ь*с. Ассоци­ ативная операция позволяет рассматривать выражения a i+аг*.. .*afc, содержащие любое конечное число элементов ai, аг, ..., а*. При этом для ассоциативной операции результат не зависит от расстановки скобок в таком выражении.

Если операция ассоциативна и коммутативна, то результат не за­ висит и от порядка расположения элементов в этом выражении.

Для алгебраической операции * часто приходится рассматривать наличие обратной операции, что равносильно решению уравнений

а* х = Ь, у* а = b

относительно элементов ж и у из А.

Решение этих уравнений приводит к правой и левой обратным операциям. В случае их существования будем говорить, что операция

* имеет обр атн ую операцию. Наличие обратной операции равно­ сильно существованию п равого и левого обратн ы х элем ентов в рассматриваемом множестве для любого элемента этого множества.

Если правый и левый обратные элементы для элемента а совпа­ дают, то этот совпадающий элемент называют обратны м элемен­ том к элементу а. В случае, когда алгебраическая операция является сложением, обратный элемент к элементу а называют п роти воп о­ лож ным элем ентом для элемента а и обозначают символом —а; в случае, когда алгебраическая операция является умножением, обрат­ ный элемент к элементу а обозначают символом а-1 .

Группой называют множество с одной ассоциативной и обраг тимой операцией. Если алгебраическая операция на группе явля­ ется сложением, то группу называют аддитивной; если алгебраиче­ ская операция на группе является умножением, то группу называют мультипликативной. Группу с коммутативной операцией назы­ вают абелевой; группу, состоящую из конечного числа элементов, называют конечной группой, а число элементов в группе - поряд­ ком группы .

Примерами групп являются:

1.Множество всех целых чисел относительно сложения.

2.Множество всех четных чисел относительно сложения.

3.Множество всех чисел, кратных данному числу п, относительно сложения.

4.Множества всех рациональных, действительных и комплексных чисел относительно сложения.

5.Множества всех рациональных, действительных и комплексных чисел, отличных от нуля, относительно умножения.

6.Множество всех корней n-й степени из единицы относительно умножения.

7.Множество всех невырожденных матриц n-го порядка относи­ тельно умножения.

8.Множества всех векторов-отрезков на прямой, на плоскости и в пространстве относительно сложения.

Если дана какая-либо группа G и подмножество Я, содержащееся в G, образует группу относительно алгебраической операции, задан­ ной в G, то группу Я называют подгруппой группы G. Например, аддитивная группа всех четных чисел является подгруппой аддитив­ ной группы всех целых чисел, которая сама является подгруппой ад­ дитивной группы всех действительных чисел.

Пусть в множестве К введены две операции - операция сложения и операция умножения. Говорят, что эти операции связаны законом дистрибутивности, если для любых элементов а, 6, с из К выпол­ няются соотношения

(а + Ь)с = ас + 6с,

а(6 + с) = аЬ+ ас.

Множество К называют кольцом, если в нем определены ассоциаг тивные операции сложения и умножения, связанные законом дистри­ бутивности, причем операция сложения коммутативная и обладает обратной операцией вычитания.

Кольцо называют ком м утативны м , если в нем операция умно­ жения коммутативная, и неком м утативны м - в противном случае. Заметим, что любое кольцо является абелевой группой по сложению.

Ненулевые элементы кольца, произведение которых равно, нулю, называют делителями нуля. Например, делителями нуля в кольце матриц второго порядка относительно сложения и умножения матриц являются ненулевые матрицы

так как их произведение

является нулевым элементом этого кольца.

Коммутативное кольцо Р, в котором есть единичный элемент и каждый ненулевой элемент имеет себе обратный элемент, называют полем.

Примерами полей являются:

1.Множества всех рациональных, всех действительных, всех ком­ плексных чисел относительно сложения и умножения.

2.Множество всех чисел вида а + 6\/2, где а и Ь- рациональные числа, относительно сложения и умножения.

3.Множество, состоящее из двух элементов 0 и 1 относительно операций сложения и умножения, заданных равенствами

0 + 0 = 0,

04-1 = 1 + 0 = 1

1 + 1 = 0,

0-0 = 0,

0 1 = 1-0 = 0,

1-1 = 1.

В любом поле отсутствуют делители нуля. Поэтому в любом поле из равенства аЬ = О следует, что либо а = 0, либо 6 = 0. Кроме того, во всяком поле сохраняются все правила действий над обык­ новенными числами. Это позволяет элементы любого поля называть числами, а само поле — числовым полем.

Из определения поля также следует, что в любом числовом поле выполняются следующие свойства его операций и элементов:

A.В поле Р определена операция сложения, которая каждой паре элементов а и 6 поля Р ставит во взаимно однозначное соот­

ветствие элемент а + 6 из Р, называемый суммой элементов а

и6. Причем:

1)Сложение коммутативно, т.е. а + 6 = 6 + а для любых эле­ ментов а и 6 из Р.

2)Сложение ассоциативно, т.е. а+ (6+ с) = (а + 6 )+ с для любых элементов а, 6, с из Р.

3)В множестве Р существует единственный нулевой элемент 0

такой, что а + 0 = а для любого элемента а из Р.

4)В множестве Р для любого элемента а из Р существует един­ ственный противоположный элемент —а такой, что а+ (—а) = 0. Это обеспечивает наличие в поле Р операции вычитания.

B.В поле Р определена операция умножения, которая каждой паре элементов а и 6 из Р ставит во взаимно однозначное соответ­ ствие элемент аб из Р, называемый произведением элементов а и 6. Причем:

1)Умножение коммутативно, т.е. аЬ= 6а для любых элементов

аи 6 из Р.

2)Умножение ассоциативно, т.е. а(6с) = (аЬ)с для любых эле­

ментов а, 6, с из Р.

3) В множестве Р существует единственный единичный эле­ мент 1 такой, что а 1 = 1 •а = а для любого элемента а из

Р.

4) В множестве Р для каждого ненулевого элемента а из Р существует единственный обратный элемент а” 1 такой, что а • а-1 = а” 1 •а = 1. Это обеспечивает наличие в поле Р операции деления на любой ненулевой элемент.

C.В множестве Р сложение и умножение связаны законом дистри­ бутивности, т.е. для любых элементов а, 6, с из Р выполняется соотношение (а + 6)с = ас + 6с.

1.8.Упражнения

1.Решить методом Гаусса системы уравнений:

1)

 

 

 

 

 

 

 

 

{

2)

 

 

 

 

 

 

 

 

 

 

r 2xi + Х 2 + Хз + Х4 = 1,

 

 

xi — 2x 2 +

Зхз — 4x 4 =

О,

 

 

 

 

3xi + 4X2 —Хз —Х4 = О,

 

 

 

2xi — 4х2 + 5хз 4* 7x 4 =

О,

 

 

 

 

XI + 3X2 —Хз +

= 2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ 5xi —3x2 + бхз -f-3X4 = 3,

 

 

 

 

 

 

 

 

 

 

 

 

 

3)

3xi + 4x2 + 5хз + 7x4 = 1,

 

 

6xi

12 x 2 4" 17хз ” 9x 4 = О,

 

 

 

 

 

7xi

14x2 4" 18хз 4” 17x4 =

0.

 

 

 

2xi + 6x2 —Зхз + 4x4 = 2,

 

4)

 

 

 

 

 

 

 

 

 

 

i

 

4xi + 2x 2 +

13хз +

10x4 =

0, {

xi 4- 2x2 — хз — 3x4 4- 4x5 =

1,

1,

к 5xi

■+■2 1хз

13x4 =

3,

 

 

 

2xi

Х2 4- Зхз 4- 2x 4 — xs =

2 ,

5)

 

 

 

 

 

 

 

 

 

6)

4” 4x2 4* 2хз — 5x 4 4" 3x 5 == 3,

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

2xi — Х2 4“ Зхз 4” 2x 4 =

1,

 

 

xi2xi4"153x2 4"- 4бхз 719x44*4"39xs5 =:=1,

 

 

 

 

3xi 4- Z2 — 5хз 4- 6x 4 4- Х5 =

О,

 

^

 

3xi 4* 3x 2 4" Зхз 4" 2x 4 = 1,

 

 

 

 

 

i

xi

— 7x2 4" 13хз — 20x 4 4" 5x 5 =

2,

 

 

3xi — Х2 хз 4- 2x4 =

—1,

 

 

 

 

 

4xi 4- 5x 2 — 14хз 4- 19x4 — хь =

—1,

к 3xi — Х2 4- Зхз —Х4 =

—1,

 

 

 

,

5xi

— 13x2 4- 21хз — 34x 4 4- llx s

= 4.

 

 

 

 

 

 

 

 

 

 

2 .

Найти произведение матриц:

 

 

 

 

 

 

 

 

 

 

 

 

1)

 

 

 

 

 

 

 

 

 

t

 

1

 

2 3

 

 

 

1 - 2

)•

/ 3 1 2 4 / 2 1

 

 

 

 

 

 

 

 

( 2

1 3 ) ( 2 2

 

 

- 2

 

1 4

 

 

 

 

2 -1

3)V

 

1

1

1 /

 

3

V

 

 

- 3

- 4

1

 

 

-3

1

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

1 2 1

 

2 3 1

 

 

 

3 5 —6 ^

 

2

 

1

8

 

 

1 3 1

 

3 1 2

 

 

 

2 4

3

 

 

-3 -1

-2

 

5)

 

2 1 3

 

1 2 1

 

 

 

3 1

1 )

 

4

 

5 - 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 1

2

 

3

 

4

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

1 2

 

 

О

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\ 2

3

1

 

4

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

Вычислить определители:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

1

2

3

4

 

2)

 

 

3

 

4

 

3)

 

 

5

6

 

 

 

 

 

 

1

 

2

 

 

3

 

4

 

 

 

 

3

6 8

11

»

 

3

 

4

5

 

6

>

1

 

■ 2

3 4

 

 

 

 

7 13 20 26

 

5

 

6 7

 

9

5

 

1 2 4

 

 

 

31

23

55

42

 

 

31

 

23

55

 

42

 

8

 

7

1

5

 

 

 

 

1

2

3

4

 

5)

 

3

4

 

 

 

в)

 

3

4

 

 

 

 

 

 

1

2

 

 

 

1

2

 

 

 

- 2

1

5

6

)

 

2

3

4

 

1

 

>

 

2

3

4

1

 

 

 

- 3

- 5

1

7

 

3 4 1

2

 

 

3 4 1

2

 

 

 

- 4

-6

■7 1

 

 

4 1

2

3

 

 

 

4 1

2

3

 

 

4.Применяя теорему Лапласа, вычислить определители:

1)

 

 

 

 

3

 

 

2)

 

0

0

0

0

 

7

 

6

5

4

2

 

1

2

 

9

 

7

8

9

4

3

 

3

4

0

0

0

0

 

7

 

4

9

7

0

0

 

7

6

5

4

0

0

J

5

 

3

6

1

0

0

 

2

3

4

5

0

0

0

 

0

5

6

0

0

 

5

1

2

6

7

3

 

0

 

0

6

8

0

0

 

2

7

5

3

4

1

 

1

5

- 5

3

 

4

 

4)

 

 

 

3

 

 

 

 

2

2

1

 

2

 

2

- 4

 

4

3

 

б

3

3 - 2

 

7

5

 

-1

 

3

 

- 1

5

 

- 9

- 5 1

3

- 1

 

- 5

- 3

 

-2

- 7

 

7

6

 

8

4

5 -6

 

4

2

 

-4

 

5

- 3

2

 

- 1

- 2

2

- 3

 

3

1

 

-2

5. а) Представить произведения определителей в виде одного определи­

теля:

 

 

 

 

 

 

 

 

 

3)

 

 

 

 

1)

 

1

2

3

2 )

3

5

 

 

 

1

2

1

3

4

1

2

1

2

3

3

1

4

2

6

1

2

3

4 .

5

б

3

5

1

4

2

3

1

4

3

1

1

2

5

 

 

 

 

 

 

б) Представить произведения определителей третьего порядка в виде опре­

делителя этого же порядка:

 

 

 

1)

3

1

1

2

5

2)

2

5

4

3

1

2

4

3

 

0 2

- 1

3

6

О

1

2

1

2

4

1

- 1

2

1

со

 

to

 

- 1

- 3

5

2

1

- 1

в) Вычислить определитель А путем умножения его на определитель 6:

1

1)

А =

-1

-1

 

 

 

 

2

 

 

-1

2)

А =

- 5

-12

 

 

 

 

9

2

 

3

 

4

0

 

со

 

ОО

1

1

 

1

 

0

- 1 3

3

<£>

5

to

15

 

 

СО

 

 

1

 

 

5

 

3

—2

-6

 

1

1

 

0

-

2

1

 

 

1

-2

- 3

11

, 6

=

0

1

0

2

0

0

1

1

 

 

 

 

0

0

0

1

 

 

1

0

0

0

6

=

-2

1

0

0

3

2

1

0

 

 

 

 

- 3

4

2

1

в. По формулам Крамера решить системы:

1)

х\ -f- 2x 2 Н” Зхз — 2x 4 =

б,

2)

xi + 2x 2 + Зхз + 4x 4

=

5,

r

'

^

2xi — Х2 2хз — 3x 4 =

8,

^

2xi + ^2 + 2хз + 3x 4 = 1 ,

 

3xi •+■2x2 —Хз + 2x4 = 4,

 

3xi

+

2x2 +

Яз + 2x 4

=

1>

k 2xi —3x2 + 2хз -f Х4 = —8,

,

4xi

+

3x 2 +

2хз + Х4

=

—5,

3)

 

 

 

2xi —Х 2 4 Зхз 4 2 x4 = 4,

Х2 Зхз 4 4X4 = —5,

Xi

— 2хз 4 3x4 = —4,

3xi 4 3x2 4 Зхз. 4 2 x4 = б,

3xi 4 2x2

— 5x4 = 12,

3xi —Х2

—Хз 4 2x4 = 6,

4xi 4 3x2 — 5хз = 5,

 

3xi —Х2

4 Зхз — Х4 = 6,

5)

 

 

 

3xi 4 5x 2 Зхз ■+•2x 4 = 12,

3xi 4 3x2 4 4хз — 5x4

=

5xi “

7x2 4 8хз 4 2x4 = 18,

4xi — 2x 2 4 5хз 4 3x 4 =

27,

4xi 4 5x2

7хз — 3x4

= ” 5,

7xi 4 8x 2 хз 4 5x4 =

40,

7xi 4 8x2 4 Зхз 4 4x4

= ” 2,

6x 1 4 4x 2 4 5хз 4 3x 4 = 41.

Найти обратные матрицы к следующим матрицам:

8. Решить системы А Х = b из упражнения б по формуле X = А ХЪ:

9. Решить матричые уравнения:

2)

3)

 

 

*

0

 

-)-(-)•

4)

(6 :)•'■(; -)■

5)(

 

 

(

 

 

 

3

X

5

-3 ^ _ ( 1 2 \

х

( 3 -5 \ _ / 9 -15 \

 

10

~6 ) ~ \ 3 4 ) '

Л

\ 9 - 1 5 ) ~ \ б —10 J

10.Для матриц из упражнения 7 построить L U -разложение А = L U и,

пользуясь соотношением A~~l = (L U )“ * =

найти обратные маг

трицы к данным.