книги / Линейная алгебра
..pdfanXi+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пТп) |
Хг |
= |
0т+ <*г,г+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.Определители
Определителем |Л|п-го порядка, соответствующим квадратной