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

SalimovRB_matematika_gl_1-20_2011_web

.pdf
Скачиваний:
42
Добавлен:
12.05.2015
Размер:
3.88 Mб
Скачать

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

 

a11

am1

 

 

 

 

 

 

 

A

 

.

 

a

a

mn

 

 

 

1n

 

 

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

a11

 

a1n

 

A

 

 

.

a

 

a

 

n1

 

nn

Элементы a11, a22 , , ann образуют главную диагональ матрицы. Число n

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

равное

A

 

a11

a1n

 

.

 

 

 

 

 

 

 

 

 

an1

 

ann

 

 

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

1

 

0

 

 

 

 

 

 

E

.

 

0

 

1

 

 

 

Матрица, состоящая из одной строки, называется строчной и обозначается Y y1, y2 , , yn . Матрица, состоящая из одного столбца, называется столб-

цевой, например,

 

 

 

x1

 

 

 

 

x

 

 

 

 

X

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

n

 

Пусть даны

две

матрицы с

одинаковым числом строк и столбцов:

A aij , B bij ,

i 1,

2, , m, j 1,

2, ,

n. Эти матрицы называются равны-

ми друг другу (при этом пишут A B

или aij bij ), если все их соответству-

ющие элементы равны друг

другу, т. е.

aij bij

для всех

i 1, 2, , m, j 1, 2, , n.

61

5354.ru

C AB

Суммой матриц A и B называется матрица, обозначаемая C A B , элементы которой cij aij bij для всех значений i, j . Это правило можно записать

так: aij bij aij bij . Аналогично вводится понятие разности двух матриц.

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

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

Умножение матриц. Даны матрица

A aij , имеющая m строк и

k

столбцов, и матрица B bij , имеющая k

строк и n столбцов. Произведением

этих матриц называется матрица, обозначаемая C AB ( A – первая матрица),

элементы cij которой определяются формулой

 

cij ai1b1 j ai2b2 j aik bkj , i 1, 2, , m , j 1, 2, , n .

(6)

Изобразим схематично эти матрицы и их произведение:

 

a11

a1k b11

b1 j b1n

c11

 

 

 

 

 

 

 

 

 

 

ai1

aik

 

ci1

a

m1

a

b

 

b

b

c

m1

 

 

 

mk

k1

 

kj

kn

 

 

c1 j

cij cmj

c1n .

cin cmn

В формуле (6) первые индексы означают номера строки элемента матрицы, вторые – номера столбца элемента. Формула (6) показывает, что элемент cij i

-й строки и j -го столбца матрицы равен сумме произведений элементов i -й строки первой матрицы A на соответствующие элементы j -го столбца второй матрицы B . Следовательно, чтобы получить элементы ci1, ci2 , , cin i -й строки матрицы C AB , нужно элементы i -й строки A умножить на соответствующие элементы первого столбца B , и, сложив, найти ci1 . Умножив элементы i -й строки A на соответствующие элементы второго столбца B и сложив, получим ci 2 и т. д. Умножив элементы i -й строки A на соответствующие элементы n -го столбца B и сложив, получим cin .

Таким образом, элементы i -й строки матрицы С получаются с помощью i -й строки первой матрицы A . Это относится к любой строке матрицы С. Поэтому ясно, что число строк С равно числу строк A , а число столбцов C равно числу столбцов матрицы В, так как номер столбца j элемента cij совпадает

с номером столбца j матрицы B .

62

5354.ru

Аналогично найдём C1 BA , если число столбцов матрицы B равно числу

строк матрицы A . Если это не так, то произведения BA не существует. Если даже AB и BA существуют, то легко проверить на примерах, что, вообще говоря, AB BA .

Свойства умножения матриц. Пусть даны три матрицы A , B и C . Тогда

A BC AB C ; A B C AB AC .

Пусть A – квадратная матрица, а E – единичная матрица того же порядка, что и A . Нетрудно проверить, что AE EA A .

Обратная матрица. Пусть дана квадратная матрица

 

 

 

a11

 

a1n

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

a

 

a

 

 

 

 

 

 

 

n1

 

 

nn

 

 

 

Определитель этой матрицы есть число

 

 

 

 

 

 

 

 

 

 

A

 

a11

a1n

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1 ann

 

 

 

 

 

Пусть этот определитель не равен нулю и

Aij

– алгебраическое дополнение

для элемента aij .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обратной к данной матрице A называется матрица, обозначаемая A 1

и

равная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A11

 

 

 

An1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

A1n

 

 

 

 

Ann

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

Нетрудно проверить, что AA 1

A 1 A E.

 

 

 

 

 

 

 

 

Отсюда видно, что для построения обратной матрицы A 1 для матрицы

A

нужно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

элементы A заменить на их алгебраические дополнения;

 

все дополнения поделить на A

– определитель матрицы A ;

 

полученную матрицу транспонировать.

Из приведенного определения видно, что для нахождения A 1 нужно вычислить определитель матрицы A и все алгебраические дополнения для

всех ее элементов.

63

5354.ru

 

§ 4. Системы n линейных алгебраических уравнений с n неиз-

вестными. Матричный метод решения

Дана система уравнений

a

 

x

a

 

x

... a

x

b ,

 

 

 

11

1

12

2

1n

n

 

1

 

 

 

a21x1

a22 x2

... a2n xn

 

b2

,

,

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

n1

x

a

 

x

... a

x

 

b ,

 

 

 

1

 

n2 2

nn n

 

n

 

 

 

где x1, x2 , , xn – искомые неизвестные,

a11, a12 ,

,

ann

 

заданные числа,

называемые коэффициентами уравнений системы, b1, b2 , , bn – заданные числа, называемые свободными членами системы уравнений. Нужно найти x1, x2 , , xn . Введём три матрицы

a11

 

 

a1n

 

 

A

 

 

 

 

,

(8)

a

 

 

a

 

 

 

n1

 

 

 

nn

 

 

 

x1

 

 

 

 

X

 

x

 

,

 

(9)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

n

 

 

 

 

 

b1

 

 

 

 

B

b

 

 

 

(10)

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

n

 

 

 

 

A называется матрицей коэффициентов системы (7),

X

матрицей неиз-

вестных, B матрицей свободных членов. Определитель матрицы A называ-

ется определителем системы и обозначается . Итак, определитель системы

(7) равен

 

 

 

 

a11

a1n

 

.

(11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an1 ann

 

 

 

Возьмём произведение AX

матриц (8) и (9). Так как X – столбцевая мат-

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

 

a11x1 a12 x2 a1n xn

 

 

 

 

AX

a21x1 a22 x2 a2n xn .

 

 

 

 

 

 

 

 

 

 

 

a

x a

n2

x

2

a

nn

x

 

 

 

 

 

 

 

n1 1

 

 

 

 

 

n

 

 

 

 

 

 

64

 

 

 

 

 

 

 

 

 

 

 

5354.ru

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Элементы этого произведения равны согласно системе (7) свободным членам соответствующих уравнений этой системы, т. е. соответствующим элементам матрицы B . Следовательно, эти две матрицы равны друг другу. Таким образом,

AX B .

(12)

Это есть матричная запись системы (7).

Пусть определитель системы (7), т. е. определитель (11), отличен от нуля. Тогда по известной матрице (8) коэффициентов системы (7) найдём для неё обратную матрицу A 1 . На эту матрицу (все элементы которой известны) умножим обе части (12), считая матрицу A 1 первой матрицей в произведениях, и получим

A 1 AX A 1B .

(13)

Согласно первому свойству умножения матриц, левая часть формулы (13) равна AA 1 X , но так как A 1 A E , EX X , то левая часть формулы (13) равна

X . Таким образом,

 

X A 1B .

(14)

Правая часть формулы содержит известные матрицы. Найдём произведение A 1B . Это будет столбцевая матрица с известными элементами, но эта матрица по формуле (14) равна матрице неизвестных X . Поэтому их соответствующие элементы равны друг другу. Приравняв эти элементы, найдём неизвестные x1, x2 , , xn .

§ 5. Формулы Крамера

Покажем, что решение системы (7) определяется формулами Крамера

x1

,

x2 , ,

xn

n .

 

(15)

 

 

 

 

 

 

 

 

Здесь – определитель системы (7), считается,

что 0 .

1,

2 , ,

n

определители, получаемые из определителя заменой соответственно первого, второго, … n -го его столбца на столбец свободных членов системы (7), т. е.

 

b1 a12

a1n

,

 

a11 b1

a1n

,

 

 

 

 

 

bn an2

ann

 

 

an1 bn

ann

 

 

 

 

 

…….

 

 

 

65

5354.ru

n

a11

a12

b1

.

 

 

 

 

 

an1

an2

 

bn

 

Запишем разложение определителя (11) системы (7) по элементам первого столбца:

 

 

a11

a1n

 

a11 A11 a21 A21 ... an1 An1.

(16)

 

 

 

 

 

 

 

 

 

an1

 

ann

 

 

 

Заменим

элементы первого столбца a11, a21, ,

an1 соответственно

на

b1, b2 , ,

bn – свободные члены системы (7). Тогда

 

 

 

b1 A11 b2 A21

bn An1 .

(17)

Получили разложение определителя по элементам первого столбца. Ана-

логично запишем разложение определителя

по элементам второго столбца

b1 A12

b2 A22 bn An2

(18)

и т. д. Наконец, получим разложение определителя n по элементам последнего столбца:

n b1 A1n b2 A2n bn Ann .

(19)

По формуле (14) будем иметь

x1

x2xn

 

A11

 

An1

 

b

 

 

A

A

 

 

1

 

 

 

 

 

 

 

b

.

 

 

2

 

A

 

A

 

 

 

 

 

 

 

 

1n

 

nn

 

 

 

 

 

A

 

 

bn

 

 

 

A

 

 

 

В правой части матрицы перемножим и получим столбцевую матрицу. Теперь последнюю формулу запишем так:

x1

x2xn

b1 A11 b2 A21 bn An1

 

 

 

 

 

 

 

A

 

 

 

 

b1 A12 b2 A22 bn An2

 

 

 

 

 

 

 

A

 

.

 

 

 

 

 

 

b1 A1n

b2 A2n bn Ann

 

A

 

 

 

66

5354.ru

Согласно (17) – (19) в последней формуле суммы, стоящие в числителях матрицы правой части, равны соответственно 1, 2 , , n . Следовательно, эту формулу можно записать в виде

x1

x2xn

1

/ A

 

 

2

 

 

 

 

/ A

.

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

/ A

 

В этом соотношении матрицы слева и справа равны друг другу, следовательно, их соответствующие элементы равны, т. е. получаем соотношение (15). Из формул Крамера вытекает следующая

Теорема. Если определитель системы (7) не равен нулю, то эта система имеет единственное решение (которое можно найти, например, по формулам Крамера).

§ 6. Общая система линейных алгебраических уравнений. Метод Гаусса

Дана система m линейных алгебраических уравнений с n неизвестными x1, x2 , , xn

a

x

a

 

x

... a

x

 

b ,

 

 

11 1

12

 

2

 

1n n

 

1

 

 

a21 x1

a22 x2

... a2n xn

b2

,

(20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x

a

m2

x

2

... a

mn

x

n

b .

 

 

m1 1

 

 

 

 

 

m

 

Здесь коэффициенты a11, a12 , , amn

и

свободные члены b1,b2 , ,bm – заданные

числа. Будем считать, что число m уравнений не больше числа n неизвестных (случай m n требует особого рассмотрения).

Система (20) называется совместной, если она имеет решение, т. е. существуют числа x1, x2 , , xn , удовлетворяющие всем уравнениям системы. Си-

стема называется несовместной, если она не имеет решения. Две системы называются равносильными (эквивалентными), если любое решение одной из них является решением другой, и наоборот.

Следующие преобразования, называемые элементарными, переводят заданную систему в равносильную (эквивалентную) ей:

перестановка любых двух уравнений системы;

умножение любого уравнения системы на ненулевое число;

67

5354.ru

прибавление к обеим частям данного уравнения соответствующих частей другого уравнения, умноженных на любое ненулевое число.

Если в процессе элементарных преобразований системы (20) появится уравнение вида 0x1 0x2 0xn 0 , то это соотношение отбрасывается, так

как ему удовлетворяют любые значения неизвестных x1, x2 , , xn (что приводит к уменьшению числа уравнений системы). Если в процессе элементарных преобразований системы (20) появится соотношение 0x1 0x2 0xn b, b 0 ,

т. е. противоречивое соотношение, которое не может быть выполнено, то система (20) является несовместной.

Метод Гаусса заключается в следующем. Пусть a11 0 (если a11 0 , то

переставим уравнения так, чтобы в первом уравнении коэффициент при первом неизвестном не равнялся нулю, или с этой же целью перенумеруем неизвестные, что приведет к перестановке соответствующих столбцов коэффициентов). Из всех уравнений, кроме первого, в системе (20) исключим неизвест-

ную x1 ,

для этого ко второму уравнению прибавим первое, умноженное на

a21 / a11 ,

к третьему уравнению прибавим первое, умноженное на

a31 / a11 , и

т. д. Тогда придём к системе вида

 

 

 

 

 

 

a

x a

x

... a

x

b ,

 

 

11

1 12

2

1n

n

1

 

 

 

 

 

 

 

 

 

 

a22 x2

... a2n xn b2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

am2 x2

... ann xn bm .

 

Пусть a22' 0 (если a22' 0 , то снова переставим уравнения или перенумеруем неизвестные x2 , , xn ). Теперь аналогично предыдущему из всех уравнений, кроме первого и второго, исключим x2 . Если система (20) совместна, т. е. при указанных преобразованиях в ней не окажется противоречивого соотношения 0x1 0x2 0xn b, b 0 , процесс продолжим. В конечном счёте путём вышеуказанных преобразований придём к одному из следующих случаев:

к ступенчатой системе

b x

b x

... b x

... b x

B ,

 

 

 

11 1

12 2

1r r

1n n

1

 

 

 

 

b22 x2

... b2r xr ... b2n xn B2

,

(21)

 

 

 

 

 

 

 

 

 

 

 

brr xr ... brn xn Br ,

 

 

 

 

 

68

5354.ru

здесь число

уравнений

r n ,

так как система содержит неизвестные

xr 1, xr 2 , , xn

(если xr 1,

xr 2 , ,

xn не входят в систему (21), то их не будет и

висходной системе (20));

к треугольной системе

b x

b

x

... b x

B ,

 

 

 

11 1

12

2

1n n

1

 

 

 

 

b22 x2

... b2n xn B2

,

(22)

 

 

 

 

 

 

 

 

 

 

 

 

bnn xn Bn .

 

 

 

 

 

 

В системах

(21), (22) по построению все коэффициенты

b11, b22 , brr , , bnn

отличны от нуля. В случае системы (20), приведённой к

системе (22), далее поступим так: из последнего уравнения (22) найдем xn ; из предпоследнего найдем xn 1 , затем xn 2 , и, наконец, x1 , т. е. найдём все иско-

мые неизвестные.

Итак, в этом случае система (20) имеет единственное решение. Определитель преобразованной системы (22) обозначим . Он равен

b11

 

b1n

 

 

b11b22 bnn 0.

0bnn

Впоследнем легко убедиться, разложив этот определитель по элементам первого столбца, в котором только один элемент ( b11 ) отличен от 0, и разложив

аналогично оставшиеся миноры также по элементам первого столбца. Определитель исходной системы (20), когда m n , обозначим . Он равен 1 , т. е. может отличаться лишь знаком от . В самом деле, прибавлению к одному из уравнений системы (20) другого уравнения, умноженного на определённое число, отвечает соответствующая операция над строками определителя , которая не изменяет этот определитель. Перестановке уравнений в исходной системе отвечает перестановка строк в определителе системы , а перенумерации неизвестных – перестановка столбцов, каждая из которых изменит лишь знак определителя. Как видно из предыдущей формулы, 1 0 , следова-

тельно, и 0 . Итак, определитель системы (20) при m n отличен от нуля, если эта система приводится к треугольной системе (22). Таким образом, при m n система (20), приводящаяся к треугольной системе, имеет единственное решение и ее определитель отличен от нуля.

69

5354.ru

m n.

Пусть система (20) приводится к ступенчатой системе (21). Перенесем в

ее правую часть все слагаемые, содержащие неизвестные xr 1,

xr 2 , , xn ,

 

b

x b x

... b x

B b x

... b x

 

 

11 1 12 2

1r r

1

1r 1 r 1

1n n

 

 

 

b22 x2 ... b2r xr B2

b2r 1xr 1

... b2n xn

(23)

 

 

 

 

 

 

 

 

 

brr xr

Br

brr 1xr 1

... brn xn .

 

 

 

 

 

В этой системе всем неизвестным xr 1, xr 2 , , xn придадим произвольные (по нашему выбору) значения. Тогда в правых частях (23) будут извест-

ные числа, и из последнего уравнения найдём

xr , из предыдущего – xr 1 и

т. д., найдём x1 . Так как значения xr 1, xr 2 , , xn

выбраны нами произвольно,

то система (23), следовательно, и (20), имеет бесконечное множество решений. Итак, система (20), приводимая к ступенчатой системе, имеет бесконечное множество решений.

Отметим, что метод Гаусса применим и в случае

§ 7. Ранг матрицы. Теорема Кронекера – Капелли

Поставим в соответствие системе (20) две матрицы

a11

 

a1n

 

 

a11

a1n

b1

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

A

 

A

b

.

a

m1

 

a

 

 

 

 

a

a

mn

 

 

 

 

mn

 

 

m1

 

 

n

Матрица A называется основной матрицей системы,

 

 

называется расши-

 

A

ренной матрицей системы (20). Элементарным преобразованиям над (20) отвечают соответствующие преобразования над строками матриц A и A . Матрица, получаемая из данной путём элементарных преобразований над строками, а также перестановкой столбцов, называется матрицей, эквивалентной данной. Основные матрицы систем (21) и (22) называются соответственно

ступенчатой и треугольной.

Строка матрицы называется нулевой, если все её элементы равны нулю, и ненулевой, если она содержит хотя бы один отличный от нуля элемент. Например, если a11 0 , a12 0 , …, a1n 0 , b1 0 , то первая строка матрицы A будет нулевой, а первая строка матрицы A будет ненулевой.

Ранг матрицы – это такое число r , что по крайней мере один определитель порядка r , получаемый из этой матрицы при удалении некоторого числа строк и (или) столбцов, отличен от нуля, а все определители (r 1) -го порядка

70

5354.ru

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]