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

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

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

anXi+ai2X2+

+

OlnXn — bl,

 

a'22X2+

+

a'2nxn = b'2,

(1.2)

 

 

 

am2*2+

 

йтпхп bm

 

В системе (1.2) следует вычеркнуть уравнения вида 0 •a?i + 0 •®2 +

... + 0 •хп = 0, если такие появились. На этом первый шаг метода Гаусса заканчивается. Элемент ац называют ведущ им элем ентом этого шага.

Следующие шаги прямого хода метода Гаусса осуществляются ана­ логично. Так, на втором шаге при а'22 Ф 0 последовательно умножают второе уравнение на a ^ /a ^ , о,\21 о!2^ •••, am2 /a22 и>соответственно, вычитаем его из 3, 4, m-го уравнений. В результате исключа­ ется неизвестное ®2 из всех уравнений, кроме 1-го и 2-го. На третьем шаге исключается неизвестное хз из всех уравнений, кроме первых трех, и т.д.

Возможно, что на некотором шаге прямого хода метода Гаусса

встретится уравнение вида

 

 

О •zi + 0 •Х2 + ... + 0 •жп —

Ь{ ф 0.

(1.3)

Тогда рассматриваемая система несовместна, и дальнейшее ее реше­ ние прекращается. Бели же при выполнении прямого хода метода Гаусса не встретятся уравнения вида (1.3), то рассматриваемая си­ стема не более чем через т шагов прямого хода преобразуется в эквивалентную систему вида

 

f flu^i + ^12^2+

+

airzr-|-

+

0>ln®n = &1,

 

i

Я22®2+

+

d2r®r+

+

fl2n®n =

&2)

(1.4)

 

 

 

 

 

 

 

 

 

 

0>rrXr+

+

arnxn =

6r.

 

Для упрощения записи в системе (1.4) штрихи над коэффициен­ тами опущены. В ней не более т уравнений, т.е. г < т , так как некоторые уравнения, возможно, были приведены к виду 0 = 0 и вы­ черкнуты.

При г = п система (1.4) имеет треугольный вид

( anzi+ai2®2+

+

ain^n =

Ьъ

 

J

<*22^2+

+

a2nSn =

&2>

(1.5)

 

 

 

 

 

Qnnxn — bn,

и в ней легко совершить обратный ход метода Гаусса. Для чего из последнего уравнения этой системы найдем значение неизвестного хп. Подставив его в предпоследнее уравнение, найдем значение xn_i. Продолжая так далее, однозначно определим все неизвестные xi, Х2,

..., хп. Следовательно, если система (1.1) при прямом ходе метода Гаусса сводится к системе треугольного вида, то такая система опре­ деленная, т.е. имеет единственное решение.

При г < п система (1.4) имеет вид трапеции. В ней неизвестные xi, Х2, ..., хг принимают за главные, а неизвестные xr+i, хг+2} •••>

хп - за свободные. Свободные неизвестные могут принимать любые фиксированные значения. Полагая

•Ег+1 — Тг+1 > ®г+2 = Тг+2)

Хп уп,

где 7г+1 , ..., 7П - произвольные постоянные, и проведя в системе обратный ход метода Гаусса, придем к системе

XI

=

01 +с*1|Г+17г+1 +

+

**1п7п >

Х2

=

02 + <*2,r+l7r+l +

+

С*2пТп)

Хг

=

+ <*г,г+17г+1+

+

агп7п)

Хг+1

=

Тг+1,

 

 

хп

=

 

 

7п,

которая является общим решением системы (1.1). Из общего реше­ ния (1.6) при конкретных значениях 7г+ъ 7г+2> •••, 7п будут полу­ чаться частные решения системы (1.1). Так как каждое свободное неизвестное может принимать бесчисленное множество значений, си­ стема (1.1) при г < п, т.е. в случае, когда она приводится к трапеце­ идальному виду, обладает бесчисленным множеством решений. Так всегда будет для совместных систем, имеющих меньше уравнений, чем неизвестных и, в частности, для однородных, имеющих меньше уравнений, чем неизвестных.

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

П ример 1. Методом Гаусса решить систему уравнений

XI

-

2х2

+

4х3

=

з,

3xi

-

я2

+

5х3

=

2,

2xi

+

S2

+

*3

=

-1 ,

2xi

-

4х2

+

Зх3

=

1.

Решение. Оставляя в расширенной матрице системы

/ 1

—2 4

3 \

3

- 1

5

2

2

1

1

- 1

2

- 4

3

1 /

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

1

- 2

4

3 \

0

5

- 7

- 7

0

5

- 7

- 7

\ °

0

- 5

- 5 /

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

1

- 2

4

3 \

0

5

- 7

- 7

 

0

0

0

0

 

\ 0

0

- 5

- 5

/

Вычеркивая здесь третью строку, придем к матрице

которая соответствует системе

Х\ — 2а?2

+

4х3

=

3,

5X2

7х3

=

-7 ,

 

 

-5 х 3

=

-5 .

Отсюда, совершая обратный ход метода Гаусса, Найдем

Xi = —1, Х2 = 0, х3 = 1.

П ример 2.Методом Гаусса решить систему уравнений

 

{ xi

+

2х2

+

Зх3 =

4,

 

2x1

+ 4х2

+

6х3 =

3,

Решение.

3xi

+

х2

х3 =

1.

Бели в расширенной матрице системы

 

/ 1

2

3

 

 

 

[

2

4

6

 

 

 

\ 3

1

—1

 

 

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

Строка (0,0,01— 5) соответствует уравнению 0 xi + 0 - x 2 + 0 x3 = = —5. Наличие такого уравнения указывает на несовместность рас­ сматриваемой системы.

Пример 3. Методом Гаусса решить систему уравнений

 

{ xi

+ 2х2

+

Зх3 =

4,

 

2xi

+ х2

х3

=

3,

Решение.

3xi

+ Зх2

+

2х3

=

7.

Элементарные преобразования прямого хода метода

Гаусса над строками расширенной матрицы системы дают следую­ щую цепочку эквивалентных матриц:

/ 1 2

 

3

4

\ / 1

 

2

 

3

4

 

1 2

1 —1

3

I ~ ( 0 - 3 - 7

-5

 

V 3 3

2

7 / \ 0 - 3 - 7

 

-5

)

( 1

2

3

 

4

(

1

2

3

4

О

- 3

- 7

 

-5

 

\ 0

- 3

- 7

-5

V о

о

О

 

О

 

 

 

 

 

 

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

Г

xi + 2х2 + Зх3

=

4,

\

—3x2 — 7хз

=

—5.

Полагая здесь х3 =

уз (уз — произвольная постоянная) и проводя

обратный ход метода Гаусса, получим общее решение:

*1 = ^(2 + 57з), ®2 = ^ (5 -7 тз), *з = 7з-

Для повышения эффективности и устойчивости метода Гаусса его модифицируют различными способами. Например, часто применяют схему, в которой на каждом шаге прямого хода ведущий коэффи­ циент выбирают наибольшим по модулю. Процесс становится бо­ лее устойчивым, если исключение неизвестных проводить с помощью вращений или отражений. О вращениях и отражениях методические указания даны в пп. 6.7 и 6.8.

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

системы

 

 

 

 

 

 

5xi

+

9х2

+

13х3

=

1,

6xi

+

5х2

-

10х3

=

з,

2xi

+

4х2

+

Зх3

=

2

целесообразно сначала из первого уравнения системы вычесть удво­ енное третье, а остальные оставить без изменения. Тогда получим систему

XI

+

*2

+

7х3

=

-3,

6xi

+

5х2

-

10х3

=

3,

2xi

+

4х2

+

Зх3

=

2,

в которой метод Гаусса проводится уже легко. Дополнительные пре­ образования совершаются также над матрицами.

В заключение отметим, что метод Гаусса и его модификации на­ ходят самое широкое применение в вычислительной практике.

1.2.Начальные сведения о матрицах. Д ействия с матрицами

Прямоугольной или ( т х п)-матрицей называется система чи­ сел, расположенных в виде таблицы

/ 0,11 Oi2

021022

\ат1 dm2

Кратко матрицу А записывают в виде А = (а1;), i = 1, 2, . . m; j = 1, 2, .. .,п. Числа ciij, составляющие данную матрицу, называют ее элементами. Первый индекс у элемента указывает номер строки, а второй — номер столбца, на пересечении которых находится этот элемент.

Матрицу называют комплексной, если хотя бы один ее элемент является комплексным числом, и действительной (вещ ественной), если все ее элементы — действительные (вещественные) числа.

Две матрицы одинакового размера т хп считают равными, если попарно равны их соответствующие элементы, т.е. элементы, стоя­ щие на одинаковых местах в этих матрицах.

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

Матрица, состоящая из нулей, называется нулевой и обознача­ ется через 0.

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

Диагональная матрица, у которой все элементы по главной диаго­ нали одинаковые, называется скалярной. Частным случаем скаляр­ ных матриц является единичная матрица

\

Матрицу, полученную из данной матрицы А заменой в ней строк соответствующими столбцами, называют транспонированной к А и обозначают через А! или Ат Если А — ( т х п)-матрица, то Ат

— (п х т)-матрица. В частности, если А — вектор-строка, то Ат— вектор-столбец, и наоборот.

Матрицу А называют сим метрической, если А = Ат ; матрицу А - комплексно сопряж енной к матрице А, если она получена из

А заменой в ней элементов на комплексно сопряженные; матрицу А*

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

Квадратная матрица А называется эрмитовой, если А* = А.

Суммой матриц А = (aij) и В = (6,; ) одинакового размера m xn называется матрица С = (cXJ) того же размера, элементы которой равны суммам aXJ- + 6XJсоответствующих элементов слагаемых ма­ триц. Таким образом,

( оц ...

а\п

( Ьц . •- bin

ЯП + 6ll •••din + bin

0>т1 •••апъп

bml •••bmn

Gml + bmi . . .

amn + bfrxn

Р азн ость матриц определяется аналогично:

 

an...

ain \

/ 6 ц . . .

bin \

/ an —6ц •.-ain —bin

\ ami ...

amn J

у 6mi ...

6mn J

у am\ bmi ...

amn bmn J

П роизведением матрицы A — (ai;) на число а называют ма­ трицу а А, все элементы которой равны произведениям соответ­ ствующих элементов исходной матрицы на это число:

( аац

otain \

a a mi

ocamn J

Вчастности, матрицу (—1) называют противоположной к матрице

Аи обозначают через —А.

Линейной комбинацией в ек тор ов -стр ок (столбцов) А \, AQ,

..., Ak называют вектор-строку (столбец)

а\А\ + OL^AI + ... + ajb-Afc.

2-1307

Сложение матриц и умножение их на числа обладают еледующими свойствами:

1)

A -f- В = В + А у

5)

1 •А = Ау

2)

Л + (Я + С) = (А + Я) + <7,

б)

(* + Р ) . А = аА + РА,

3)

А + 0 = Ау

7)

а •(А 5) OLA -f- otBу

4)

А + ( - А) = 0у

8)

а - (РА) = (ар)-А

при любых матрицах А, В, С и любых а, /?.

Умножение матриц определяется лишь для случая, когда число столбцов первого множителя равно числу строк второго множителя. Пусть

Произведением матриц А и В, заданных в указанном порядке, называется матрица С — (с,-*), элементы которой определяются по следующему правилу:

cik = anbik + ai2b2k + •••+ ainbnkl * = 1,.. .,m; к = 1,.. .,p,

т.е. элемент с,* матрицы С равен сумме произведений элементов г-й строки матрицы А на соответствующие элементы fc-ro столбца ма­ трицы В. Из этого определения следует, что матрица С будет матри­ цей размера т х р.

Например,

1-3 + 2 -2 + 3 1

1 1 + 2 •4 + 3 •3

10

18

(

 

+ 3 - 4 + 1 - 3

13

17

2 •3 + 3 •2 + 1 1 2 1

6- 3 + 1-2 + 2 1

6 •1 + 1 -4 + 2 •3

22

16

( 1 2 3 )

 

 

 

 

3 \

3 1 3 2 3-3

3 6 9

2 I ( 1 2 3 ) =

2 1

2 2 2 - 3

2 4 6

( )

1 1 1 2 1 - 3

1 2 3

Из рассмотренных примеров видно, что умножение матриц некомму­ тативно.

Умножение матриц обладает следующими свойствами:

1) А - ( В С) = {А- В) С,

2)а(А ■В) = (а ■А) ■В = А •(о •В),

3)С - ( А + В) = С - А + С- В,

4)(А + В ) С = А - С + В С .

Отметим некоторые с в о й с т в а о п е р а ц и й т р а н сп о н и р о в а н и я И к о м ­ п л е к с н о г о с о п р я ж е н и я м а т р и ц :

1)

А + В = А + В,

{аА) = а ■А,

АВ = А ■В;

2)

(А + В)г = АТ + ВТ,

(А + ВУ =А* + 5 *;

3)

(аЛ)т = а Ат,

(аЛ)* = а-А*\

4)

(АВ)Т = Вт■Ат,

(АВ)* = В*

А *;

5)

(АТ)Т = А ,

(А*У = А.

 

Умножение векторов-столбцов на числа позволяет систему линей­

ных уравнений

 

 

 

' CLil^i + <*12^2 +

+

ClinXn =

6i,

^21^1 + <*22^2 +

+

а2пХп =

62,

Q>m\X\+ ат2х2 +

+а>тпХп — Ьт

записать в векторной форме следующим образом:

ап

\

/< * 12

\

/

<*1п

 

 

XI +

 

х2 + . . . +

Хп =

<*ml

/

V <*m2

J

\ атп

а умножение матриц позволяет ту же систему записать в виде

 

А -Х

= Ь,

 

полагая

 

 

 

ain

 

®i

Ьх

А =

X

=

6 =

Под элементарными преобразованиями над матрицей по­ нимают:

1.Умножение любой г-й строки матрицы на любое число а.

2.Прибавление к любой г-й строке матрицы любой ее j строки, умноженной на число а.

3.Умножение любого j -го столбца матрицы на любое число а.

4.Прибавление к любому г-му столбцу матрицы любого ее ;-го столбца, умноженного на число а.

5.Перестановку строк (столбцов) матрицы.

Первое и второе преобразования над матрицей равносильны умно­ жению ее слева соответственно на элементарные матрицы

( 1

\

 

 

 

1

 

\

1

 

 

 

а

 

 

 

1

 

 

г

 

 

а

 

 

17

1

/

 

 

 

3

j

 

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

Применение элементарных преобразований над матрицами мы уже встречали в методе Гаусса. Применение элементарных матриц про­ иллюстрируем ниже (см. п. 1.6.) в примере на разложение матрицы на треугольные множители.

1.3.Определители

Определителем |Л|п-го порядка, соответствующим квадратной