Алгебра Лекции и практика - Южный федеральный университет

advertisement
МИНИСТЕРСТВО ОБРАЗОВАНИЯ и НАУКИ РФ
ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ
РОССИЙСКОЙ ФЕДЕРАЦИИ
Федеральное государственное образовательное учреждение
высшего профессионального образования
«ЮЖНЫЙ ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ»
Факультет математики, механики и компьютерных наук
Ростовский-на-Дону государственный университет
В.Б. Дыбин
Алгебра
Лекции и практика
Методическое пособие первокурснику
Часть 1
Модуль 2
Системы линейных алгебраических уравнений
2008 г.
Гл. 2. Системы линейных алгебраических уравнений
Матричное исчисление, построенное в лекциях I-V, позволяет провести
первичное исследование систем линейных алгебраических
уравнений
(СЛАУ), опирающееся на метод, носящий имя великого немецкого математика Карла Гаусса (1777-1855).
Поскольку СЛАУ можно трактовать как матричное уравнение специального вида, обоснование метода Гаусса даётся на языке элементарных матриц и
матричных уравнений. В свою очередь, метод Гаусса позволяет построить
удобный алгоритм решения ряда матричных уравнений, в частности, алгоритм
вычисления обратной матрицы.
Лекция VI.
План
2.1 Классификация СЛАУ
2.2 Метод Гаусса решения СЛАУ
2.1 Классификация СЛАУ
Системы линейных алгебраических уравнений составляют основной аппарат линейной алгебры, а их исследование основано на алгебре матриц.
Общий вид системы m линейных уравнений с n неизвестными даётся
формулами
111  12 2    1n n  1
           
 21 1
22 2
2n n
2
(2.1)







 m11   m 2 2     mn n   m
В этих формулах 1 , ,  n называются неизвестными,  ij – коэффициентами
при неизвестных,  i – свободными членами (или правыми частями). Коэффициенты  ij порождают матрицу
 11 12

 22

A   21
 


 m1  m 2
 1n 

  2n 
,
 

  mn 
(2.2)
2
которая называется матрицей (или основной матрицей) СЛАУ (2.1),
A  M mn (R ) . Неизвестные  j , j 1, n , порождают столбец неизвестных
 1 
 
 
A 2,

 
 
 n
(2.3)
а правые части  j , j 1, m – столбец правых частей
 1 
 
 
b   2  , bRm.

 
 
 m
(2.4)
Матрица
 11 12



 A | b    21 22
 


 m1  m 2
 1n
  2n
 
  mn
1 

2 


 m 
(2.5)
называется расширенной матрицей СЛАУ (2.1).
Арифметический вектор
 1 
 
 
a   2  , a  Rn .

 
 
 n
называется решением СЛАУ (2.1), если при подстановке чисел  1 ,  2 , ,  n в
уравнения (2.1) соответственно вместо 1 ,  2 , ,  n мы получаем систему верных числовых равенств.
СЛАУ называется совместной, если она имеет хотябы одно решение, и
несовместной, если она не имеет ни одного решения. Если СЛАУ имеет только одно решение, она называется определённой, а если число решений больше
1 – неопределённой.
В результате мы получаем следующую классификацию СЛАУ. Пусть M
– множество решений СЛАУ, а M – мощность (число элементов) множества
3
M . Тогда
СЛАУ
СОВМЕСТНЫЕ
СЛАУ
НЕСОВМЕСТНЫЕ
СЛАУ
M 
M 
ОПРЕДЕЛЁННЫЕ
СЛАУ
НЕОПРЕДЕЛЁННЫЕ
СЛАУ
M 1
M 1
Если все правые части СЛАУ равны нулю,  j  0, j 1, m , она называется однородной, в противном случае СЛАУ называется неоднородной. Однородная СЛАУ всегда совместна, т.к. имеет нулевое решение
 0
 
a     , a Rn .
 0
 
Множество всех решений СЛАУ будем называть её общим решением, а
одно фиксированное решение – частным решением. СЛАУ вида (2.1) иногда
удобно представлять в виде матричного уравнения
(2.6)
Ax  b .
где A – матрица вида (2.2), x – вектор-столбец вида (2.3), а b – векторстолбец вида (2.4). Для того чтобы убедиться, что уравнения (2.6) и СЛАУ
(2.1) равносильны (т.е. одновременно неразрешимы или разрешимы, причём в
последнем случае их решения совпадают), достаточно в левой части равенства
(2.6) провести умножение матрицы A на вектор-столбец x . Применив после
этого принцип равенства матриц, получаем, что матричное равенство (2.6) эквивалентно системе равенств (2.1).
2.2 Метод Гаусса решения СЛАУ
Метод последовательного исключения неизвестных часто применяется
при решении всевозможных систем уравнений, в частности, тех, которые
встречаются в школьном курсе математики. Если выразить, например, 1 из
первого уравнения системы (2.1) и подставить это выражение в остальные
уравнения, потом выразить  2 из второго “нового” уравнения, подставив это
4
выражение в остальные “новые” уравнения и т.д., то через конечное число шагов мы получим уравнение, содержащее только  n . Найдя решение этого
уравнения, путём обратной подстановки его в предыдущие “новые” уравнения
можно найти значения всех неизвестных  n1 ,  n2 , ,  2 , 1 . Такова идеальная
схема которая молчаливо опирается на целый ряд допущений (например, на
то, что m  n , иначе для отыскания  n не хватит уравнений).
Метод Гаусса в современном изложении представляет собой, во-первых,
такую модификацию метода исключения неизвестных, которая позволяет исследовать любые системы вида (2.1). А, во-вторых, алгоритм Гаусса, используя аппарат алгебры матриц, даёт такую формализацию метода исключений,
которая позволяет существенно сократить объём выкладок, резко возрастающий вместе с ростом числа уравнений m и числа неизвестных n .
Прежде чем приступить к изложению метода Гаусса в общем случае,
продемонстрируем его на простом примере.
Пример 1. Рассмотрим систему уравнений
 1  3 2  4
.

41  5 2  1
Вычитая из второго уравнения 4 первых, получаем систему
1  3 2  4
.


17



17

2
Разделив обе части второго уравнения на (-17), приходим к системе
1  3 2  4
.

2  1

Наконец, вычитая из первого уравнения 3 вторых, получаем систему уравнений
1
1
,



1

2
которая определяет единственное решение исходной системы
 1
x    .
 1
Теперь повторим проделанные преобразования на расширенной матрице
исходной СЛАУ, заменяя каждое преобразование соответствующим строчным
5
элементарным преобразованием:

3 4  C2  4C1


 4  5  1 ~



3 4  C2 ( 17)  3 4  C1 3C2




 0  17  17  ~  0  1 ~




 0 1


 0  1 .


Заметим, что итогом проведённых преобразований явилось поученное из
основной матрицы системы матрицы приведённого вида, которая в данном
случае совпала с единичной матрицей E (2) . Кроме того, в данном случае правые части последней системы уравнений образуют единственное решение исходной системы.
Перейдём к изложению метода Гаусса в общем случае.
Элементарными преобразованиями СЛАУ называются: перемена местами двух её уравнений, умножение обеих частей какого-либо уравнения СЛАУ
на число, отличное от нуля, добавление к левой и правой частям какого-либо
уравнения соответственно левой и правой частей другого уравнения, умноженных на произвольное число.
Метод Гаусса решения СЛАУ состоит в последовательном выделении в
каждом уравнении неизвестного, которое после этого элементарными преобразованиями СЛАУ исключается из всех её остальных уравнений.
Поскольку каждой СЛАУ вида (2.1) можно поставить в соответствие её
расширенную матрицу (2.5) и это соответствие взаимооднозначно, т.е. каждая
матрица вида (2.5) является расширенной матрицей некоторой СЛАУ, вместо
элементарных преобразований СЛАУ удобно проводить строчные элементарные преобразования её расширенной матрицы. Нетрудно заметить, что после
выделения в одном уравнении некоторого неизвестного (это неизвестное
называется ведущим в данном уравнении) и последующего его исключения из
остальных уравнений СЛАУ, соответствующая строка её основной матрицы
будет иметь приведённый вид. Поэтому будем говорить, что и уравнение, отвечающее этой строке матрицы, имеет приведённый вид. Если же каждое
уравнение СЛАУ, содержащее хотя бы один ненулевой коэффициент, имеет
приведённый вид, будем говорить, что СЛАУ имеет приведённый вид.
Ясно, что в этом случае основная матрица СЛАУ тоже имеет приведённый вид.
Алгоритм метода Гаусса распадается на 3 этапа:
1) построение СЛАУ приведённого вида, равносильной исходной
СЛАУ;
2) анализ СЛАУ приведённого вида;
3) описание общего решения.
Следующее предложение является обоснованием применимости метода
Гаусса к любой СЛАУ.
Предложение 2.1. Для любой СЛАУ вида (2.1) существует равносильная
ей СЛАУ приведённого вида.
6
◄ В силу предложения 1.3 найдётся конечное число строчных элементарных преобразований, применяя которые к матрице A , мы получим матрицу
P приведённого вида. По свойству 5) элементарных преобразований существует такая матрица C , C  GM n (R ) , что P  CA . Заметим, что матрица C является произведением элементарных матриц, отвечающих указанным выше
элементарным преобразованиям. Применив те же самые элементарные преобразования к расширенной матрице  A b  системы уравнений (2.1), получаем
матрицу системы уравнений приведённого вида
C A b   CA Cb   P Cb  .
Переходя к соответствующим матричным уравнениям Ax  b и CAx  Cb , замечаем, что в силу предложения 1.9 они равносильны. Следовательно, равносильны и отвечающие им системы уравнений. ►
Лекция VII.
План
2.3 Анализ СЛАУ приведённого вида и описание общего решения
2.4 Однородные СЛАУ
2.5 Решение матричных уравнений методом Гаусса
2.6 Отыскание обратной матрицы методом Гаусса
2.3 Анализ СЛАУ приведённого вида
Исследование СЛАУ приведённого вида распадается на три случая.
1)Система уравнений приведённого вида содержит “плохое” уравнение,
т.е. уравнение вида
0  1  0   2    0   n  c ,
(2.7)
где c  0 . Так как это уравнение не имеет решений, система уравнений приведённого вида, а с ней и исходная система уравнений несовместны.
Пример 2. Решить СЛАУ методом Гаусса и найти её общее решение
21   2  3 3   4  5
21  2 2   3  3 4  1 .
41  3 2  4 3  2 4  2
◄ Переходя к матричной записи системы и применяя подходящие
строчные элементарные преобразования, получаем систему приведённого вида, равносильную исходной системе уравнений:
7
 2 1 3  5
 2 1 3  5
 2 1 3  5

 C2 3C1 
 C3  C 2 

 2  2 1  3 1 ~    5 10 0 16  ~    5 10 0 16  .

 C3  2C1 



4

3
4

2
2
0 0 0  4


 8  5 10 0 12 
 0
Мы остановили процесс получения системы приведённого вида, так как последняя система (а с ней и исходная система) уравнений несовместна. ►
2) Система уравнений приведённого вида не содержит “плохих” уравнений. Все неизвестные в системе уравнений приведённого вида делим на две
группы: неизвестные, являющиеся ведущими в своих уравнениях, называем
связанными, а остальные неизвестные – свободными. (Случай отсутствия свободных неизвестных рассмотрим отдельно). Объявляя свободные неизвестные
параметрами c1 , c2 ,  и т.д., принимающими произвольные действительные
значения, выражаем связанные неизвестные через свободные. Полученные в
результате этого формулы определяют общее решение СЛАУ. Если свободным неизвестным придать конкретные значения, а после этого вычислить по
найденным формулам значения связанных неизвестных, мы получаем некоторое частное решение рассматриваемой системы уравнений.
Таким образом, в данном случае система уравнений приведённого вида, а
с нею и исходная система уравнений совместны, неопределенны и имеют бесчисленное множество решений.
Пример 3. Следующую систему уравнений решить методом Гаусса, найти
её общее и одно частное решения
 1  3 2  2 3   4  2
21   2  3 3   4  3 .
1  2 2  5 3
5
◄ По аналогии с предыдущим примером
 1 3 2   2 
 1  3  2   2
 0  5  7   7

 C2 C1 
 C1 C2 

1 3 ~ 
2
5 0
5  ~ 
2
5 0
5 .
 2 1 3

 ( 1)C1 
 C3  C 2 

1
2
5
0
5
1
2
5
0
5
0
0 0 0




 0
Полученная система уравнений имеет приведённый вид и не имеет “плохих”
уравнений. Её свободными неизвестными являются  2 и  3 , а связанными неизвестными 1 и  4 . Полагая  2  c1 ,  3  c2 , находим 1 и  4 из уравнений
приведённой системы,
  5c1  7c2   4  7
.

 5
1  2c1  5c2
8
Из второго уравнения 1  5  2c1  5c2 , из первого уравнения 4  7  5c1  7c2 .
Поэтому общее решение рассматриваемой системы уравнений имеет вид
 5  2c1  5c2 


c1


x
 , c1 , c2  R .
c2


  7  5c  7c 

1
2
Полагая c1 и c 2 равными, например, 1, получаем частное решение
  2
 
 1
x1    . ►
1
 
 5
Замечание. В случае неопределённых СЛАУ выбор свободных и связанных неизвестных осуществляется неоднозначно и зависит от элементарных
преобразований, применённых в алгоритме Гаусса. В связи с этим и общее
решение таких систем уравнений может иметь различную форму.
3) Система уравнений приведённого вида не содержит “плохих” уравнений и свободных неизвестных. Поскольку в этом случае все неизвестные связанные, расширенная матрица приведённой СЛАУ, возможно, после перемены
местами некоторых уравнений принимает вид
 0 0  0  1 


 0  0  0 2 


 0 0   0 3 
      

.

0
0
0


n

0 0 0  0 0










0 0 0  0 0


Ясно, что данная система уравнений является определённой, а е единственное
решение имеет вид
9
 1 
 
 
x   2  .►

 
 
 n
Пример 4. Следующую систему уравнений решить методом Гаусса
 1  2 2  3 3  6
     3  3
 1
2
3
.

2




2



1
1
2
3

 31   2  3 3  1
◄ Находим расширенную матрицу системы приведённого вида, равносильной данной системе уравнений:

2
3

 1 1 3
 2 1 2

 3 1 3



6
6
2
3
 CC22CC1 
 C1 3CC2 
3 3 1  0
9 4 3  0
3
6
~
~
1 C4 3C1  0 5 8 13  ( 1)C3  0




 0 5 12 17 
 0
1



2 3 6
 C1  2C2
2 3  C3 5C2
~
5 8 13  1 4C4

0 4 4 
 0  1 0 
  0 0 1

 C1  C2 

C

2
C
2
4
3
1
2
 0 

 0  0 
~
~ 

.
C3  2 C 4

2
0
0
0

2
0
0
0




 0 0  1
 0 0  1




Полученная приведённая система уравнений является определённой, а её
единственное решение имеет вид
1
 
x  1 . ►
1
 
Итоги изучения систем линейных алгебраических уравнений методом
Гаусса сформулируем в виде ряда предложений.
Предложение 2.2. (Критерий совместности СЛАУ). Для того, чтобы
СЛАУ была совместной, необходимо и достаточно, чтобы любая равносильная ей СЛАУ приведённого вида не содержала уравнений вида (2.7).
Предложение 2.3. (Критерий определённости СЛАУ). Для того, чтобы
СЛАУ была определённой, необходимо и достаточно, чтобы любая равно10
сильная ей СЛАУ приведённого вида не содержала уравнений вида (2.7) и свободных неизвестных.
Предложение 2.4. (Критерий неопределённости СЛАУ). Для того, чтобы СЛАУ была неопределённой, необходимо и достаточно, чтобы любая равносильная ей СЛАУ приведённого вида не содержала уравнений вида (2.7) и
имела свободные неизвестные.
Если последние два условия выполнены, СЛАУ имеет бесконечное множество решений, зависящее от произвольных постоянных c1 , c2 , , ck , количество которых совпадает с числом свободных неизвестных.
2.4 Однородные СЛАУ
Общий вид однородной СЛАУ, состоящей из m уравнений с n неизвестными, даётся формулами
111  12 2    1n n  0,
           0,
 21 1
22 2
2n n

  

 
 m11   m 2 2     mn n  0.
Как уже отмечалось выше, однородная СЛАУ всегда совместна, так как имеет
нулевое решение
 0
 
 0
x   .

 
 0
Поэтому исследование такой СЛАУ сводится к выяснению существования у
неё нулевого решения. Если нулевого решения не существует, однородная
СЛАУ является определённой и подчиняется предложению 2.3. Если же ненулевое решение существует, то однородная СЛАУ является неопределённой и
подчиняется предложению 2.4.
Следующие утверждения вытекают непосредственно из предложений
2.3 и 2.4.
Предложение 2.5. (Критерий определённости однородной СЛАУ). Для
того, чтобы однородная СЛАУ была определённой, необходимо и достаточно, чтобы любая равносильная ей СЛАУ приведённого вида не содержала свободных неизвестных.
Предложение 2.6 .(Критерий существования у однородной СЛАУ ненулевого решения). Для того, чтобы однородная СЛАУ имела ненулевое решение, необходимо и достаточно, чтобы любая равносильная ей СЛАУ приведённого вида имела свободные неизвестные.
11
Если последнее условие выполнено, СЛАУ имеет бесконечное множество решений, зависящее от произвольных постоянных c1 , c2 , , ck , количество которых совпадает с числом свободных неизвестных.
Приведём ещё одно достаточное условие для существования у однородной СЛАУ ненулевого решения.
Предложение 2.7. Если число уравнений m однородной СЛАУ меньше
числа её неизвестных n , тогда СЛАУ имеет ненулевое решение.
◄ Если m  n , тогда, как следует из метода Гаусса, можно построить
СЛАУ приведённого вида, равносильную исходной СЛАУ, у которой число
уравнений не превосходит m . Следовательно, у этой СЛАУ приведённого вида число уравнений также меньше числа неизвестных. Так как каждое уравнение этой СЛАУ, имеющее хотя бы один ненулевой коэффициент, имеет точно
одну связанную неизвестную, то у приведённой СЛАУ обязательно будут свободные неизвестные. Остаётся применить предложение 2.6. ►
2.5 Решение матричных уравнений методом Гаусса
Рассмотрим матричное уравнение
(2.8)
AX  F
в предложении, что A  GM n (R ), X , F  M nm (R ) . В силу предложения 1.7
единственное решение этого уравнения имеет вид
X  A1F
В то же время матричное уравнение (2.8) ввиду правила умножения матриц эквивалентно системе m матричных уравнений
(2.9)
AX i  F i , i 1, m ,
каждое из которых, являясь уравнением типа (2.6), равносильно определённой
СЛАУ с расширенной матрицей
A F i , i 1, m .
(2.10)
Единственное решение этой системы уравнений имеет вид


i
A1  F i  A1 F , i 1, m ,
и, как следует из 2.4, получается в результате приведения основной матрицы
A системы к виду E ,
A F  ~ E A F  .
i
1
i
Однако, ничто не мешает нам решать системы уравнений с матрицами вида
(2.10) одновременно для всех значений i 1, m . Вводя расширенную матрицу
A F  и приводя строчными элементарными преобразованиями основную
матрицу A к виду E , мы получим, что
A F   A F 1  F m  ~ E A1 F 1 A1 F m  ~ E

A 1 F .
На практике обычно возникает более общая задача решения матричного
уравнения (2.8) для произвольных матриц A  M pn (R) , X  M nm (R ) ,
12
F  M pm (R) . Изложенное выше позволяет сформулировать следующий алго-
ритм решения этой задачи.
Составляем матрицу  A F  и строчными элементарными преобразованиями приводим её к виду P B  , где P – приведённая матрица,
л-эквивалентная матрице A .
1) Если P  E (в этом случае, конечно, A – квадратная матрица), уравнение (2.8) разрешимо для любой F  M mn (R ) , а B – его единственное решение.
2) Если P  E , тогда для разрешимости уравнений (2.8) необходимо и
достаточно, чтобы у матрицы P не было нулевых строк, либо при наличии
нулевой строки, например, Pi  0 выполнялось условие Bi  0 (для каждой такой строки).
3) Если P  E и уравнение (2.8) разрешимо, то для того, чтобы его решение было единственным, необходимо и достаточно, чтобы у приведённой
СЛАУ с матрицей P не было свободных неизвестных. Если последнее условие нарушено, то уравнение (2.8) имеет бесчисленное множество решений, а
его общее решение определяется способом, описанным в пункте 2.4.
Пример 5. Решить матричное уравнение (2.8), если
 1 1
 1 1 1




A   1  1 , F   1 1 1  .
 1 1
 1 1 1




◄ Применяя метод Гаусса к расширенной матрице  A F  , получаем, что

1 1 1 1 C2 C1 
1 1 1 1 C11122CC2 

 C3 C1 
 2 
 A F    1 1 1 1 1 ~  0 2 0 0 0  ~  0
 1 1 1 1 1
 0 0 0 0 0
 0





0 1 1 1

0 0 0    P B .
0 0 0 0 
Откуда следует, что матрица P  E , P3  O и B3  O , то есть уравнение (2.8)
разрешимо. Так как у приведённой СЛАУ нет свободных переменных, то его
решение X  M 23 (R ) единственно и имеет вид
1 1 1
 . ►
X  
 0 0 0
Предлагаем читателю самостоятельно убедиться, что уравнение (2.8) с
той же самой матрицей A , но с другой правой частью
13
1 1 1 


F  1 1 1 
1 2 3 


неразрешимо.
Пример 6. Решить матричное уравнение (2.8), если
1 1 1
1 2
 , F  
 .
A  
1  1 1
3 4
◄ Применяя метод Гаусса к расширенной матрице  A F  , получаем, что
1 1 1 2  C2 C1
1 1 1 2  C1 1 2 C2
 ~ 


 0  2 0 2 2  1~2 C2
3
4
1

1
1





A F   
 0 1 2 3
  P B  .
~ 


1

1
0

0


(2.11)
Откуда следует, что матрица P  E и не имеет нулевых строк, но у приведённой СЛАУ есть одна свободная неизвестная  3 . Таким образом, уравнение
(2.8) разрешимо и имеет бесчисленное множество решений. Общее решение
1
2
X об  M 32 (R ) находим из системы (2.11), определяя  X об  и  X об  соответственно из систем P B1  и P B 2  . Именно полагая  3  c1 из системы
 0 1 2 
    2

 или  1 3
 0  0 1
 1
 2


получаем, что
 X об 
1
Полагая  3  c2 из системы
 2  c1 


  1  .
 c 
 1 
 0 1 3
    3

 или  1 3
 0  0 1
 1
 2


получаем, что
14
 X об 
2
 3  c2 
 2  c1 3  c2 




   1  , то есть X об    1
 1  , c1 , c2  R .
 c 
 c
c2 
 2 
 1
Наконец, заметим, что матричное уравнение вида
XA  F
применением к нему операции транспонирования сводится к уравнению вида
(2.8)
AT X T  F T .
Оба эти уравнения разрешимы или неразрешимы одновременно, а их решения
взаимнотранспонированны. ►
2.6 Отыскание обратной матрицы методом Гаусса
Для отыскания обратной матрицы A1 методом Гаусса достаточно в
уравнение (2.8) положить F  E . Если матрица A обратима, тогда уравнение
AX  E
имеет единственное решение
X  A1  E  A1 .
В этом случае матрица A л-эквивалентна матрице E . Если же матрица
A необратима, тогда в силу предложения 1.6. л-эквивалентная ей матрица Dr
будет иметь нулевую строку.
В результате, алгоритм выяснения обратимости матрицы A и отыскания
матрицы A1 методом Гаусса принимает следующий вид.
Составляем матрицу  A E  и строчными элементарными преобразованиями приводим её к виду E B .
1) Если указанные преобразования осуществить удаётся, матрица A обратима и A1  B .
2) Если указанные преобразования осуществить не удаётся (в процессе
преобразований появляется нулевая строка), матрица A необратима.
Пример 7. Выяснить, является ли матрица A обратимой, и в случае её
обратимости найти матрицу A1 ,
 1 2  1


A   3 5  2 .
 2 1  2


◄ Применяя только что изложенный алгоритм, получаем
15
 2  1 1 0 0 

2 1 1 0 0

 C2 3C1
 1 3C3
 3 5  2 0 1 0 ~  0 1 1  3 1 0 ~

 C3 2C1

2
1

2
0
0
1


 0  3 0  2 0 1
2
 2  1 1 0
  0  1  23 0
0
3  C C

 C1 2C3 
1
2
7
~  0 1 1  3 1
0  ~  0 0   3 1  13  ~
C2 C3
C2 C3


2
1
2
1
0

0
0

0

0
0

3
3
3
3


1
1
  0 0  83 1
  83 1
3
3




1
~  0  0 23 0  13   A  GM 3 (R ). A   23 0  13  . ►

 7 1  1 
7
1
0
0


1

3
3
3
 3

Упражнения
При решении систем линейных алгебраических уравнений, в особенности, на первых парах у решающих часто появляются вычислительные ошибки.
В связи с этим полезно использовать вычислительную схему методы Гаусса с
контрольным столбцом. Контрольный столбец является суммой всех столбцов
расширенной матрицы СЛАУ и выписывается справа от неё. Элементарные
преобразования метода Гаусса применяются к контрольному столбцу наравне
с остальными столбцами. Контроль за правильностью вычислений состоит в
том, что после каждого элементарного (строчного) преобразования измененный элемент контрольного столбца должен равняться сумме всех предыдущих
элементов соответствующей строки. Если это условие нарушается ,следует
искать вычислительную ошибку!
Пример 8. Решить методом Гаусса следующую систему уравнений.
 x1  x2  x3  2
2 x1  3x2  x3  7 .
4 x1  x2  2 x3  1
◄ Выписываем расширенную матрицу системы с добавленным контрольным столбцом
  1   1  2  3


3
1 7 13  .
 2
 4 1  2 1 0


16
В качестве ведущего выбираем второй элемент первой строки и проводим элементарные преобразования C 2  3C1 , C3  C1 . Получаем матрицу
  1   1  2  3


4 13 22 
 1 0
 3 0  3  3  3


Проводим контроль вычислений второй и третьей строк:
 1  0  4  13  16  22, 3  0  3  3  3 .
Следовательно, при вычислении второй строки допущены ошибки. Действительно, в первом столбце второй строки после преобразования C2  3C1 должно быть 2  3(1)  5 . Остальные вычисления правильны: 5  4  13  22 . Проведем исправления, продолжаем вычислительный процесс:
  1   1  2  3 1C   1   1  2  3

3 3
 C1 C3
4 13 22  ~  5 0 4 13 22  ~
 5 0
 3 0  3  3  3
  0  1  1  1 C2 5C3




 0   2  3  4


9 16 27 
 0 0
  0  1  1  1


Вновь проводя вычислительный контроль для первой и второй строк, получаем:
0  1  2  3  4, 0  0  9  16  25  27 .
Вновь ошибка вычислений во второй строке, вместо 16 следует записать 18.
Далее,
 0   2  3  4  1C  0   2  3  4 

9 2
 C1  2C2
0
0
9
18
27
~
0
0

2
3



 ~
  0  1  1  1
  0  1  1  1 C3 C2




 0  0  3 2


0
0

2
3


 0 0 1 2 

.
Проводя последний вычислительный контроль замечаем, что в первой строке
допущена ошибка внимания, в четвертом столбце элемент  3 не преобразован:  3  4  1 . Наконец, матрица
 0  0 1 2


0
0

2
3


 0 0 1 2 


является приведенной
17
 x1   1 
   
x   x2    1  .
 x   2
 3  
Поскольку ошибки внимания преследуют неопытного вычислителя на
каждом шагу, следует убедиться в правильности полученного результата. Для
этого нужно сделать проверку, подставив полученные значения неизвестных в
исходную систему и убедиться, что после этого её уравнения превращаются в
числовые равенства. ►
Следующие системы уравнений решить методом Гаусса, указав общее и
одно частное решения.
1.
2.
x1  2 x2  3x3  4 x4  7
x1  2 x2  3x3  1
2 x1  x2  2 x3  3x4  6
x1  2 x2  3x3  2
3x1  2 x2  x3  2 x4  7
 x1  2 x2  3x3  3
4 x1  3x2  2 x3  x4  18
3.
.
x1  2 x2  3x3  6
x1  2 x2  3x3  4 x4  1
4.
3x1  x2  2 x3  x4  x5  1
2 x1  3x2  4 x3  5 x4  2
2 x1  x2  7 x3  3x4  5 x5  2
3x1  4 x2  5 x3  6 x4  3
x1  3x2  2 x3  5 x4  7 x5  3
4 x1  5 x2  6 x3  7 x4  4
3x1  2 x2  8 x3  5 x4  8 x5  3
5.
6.
2 x1  x2
0
2 x 2  x3
4 x1  5 x2
3
4 x1  5 x2  3x3  x4  x5  3
2 x3  x 4  1
7.
8.
2 x2  3x3  4 x4  5 x5  1
2 x1  3x2  4 x3  5 x4
 x4  6 x5  3
2
2 x3  3x4  4 x5  0
11x1  22 x2  0
 5 x1  17 x2  7
6 x2  13 x3  5
x2  4 x3  7
8 x4  3x5  31
4 x4  2 x5  10
18
9.
2 x3  6 x4
 x6  0
 2 x1  4 x2
 x5  0
10.
x1  x2 
11.
 x4  2 x5  0
 x2
 x3  x4  0
x4  0
x1  2 x2  3x3  4 x4  0
x1  x3  x6  0
2 x2
x3 
x1  3x2  6 x3  10 x4  0
x1  4 x2  10 x3  20 x4  0
4 x1  2 x2  6 x3  3x4  4 x5  0
x1  x2
2 x1  3x2  2 x3
 2 x4
0
 x5  0
Если коэффициенты системы уравнений зависят от параметров, то задача решения такой системы сводится к определению тех значений параметров,
при которых система имеет решения и нахождению этих решений.
Если параметр  входит в один коэффициент системы, целесообразно
при использовании метода Гаусса проводить такие элементарные преобразования, которые не приводят к «расползанию» параметра по матрице.
Пример 9.
x1  2 x2  x3  2
x1  3x2  4 x3  5 .
 2 x1 
x2  7 x3  6
◄
5 0  15  10  1 C
 1 2 1 2


 C1 2C3 
5 1

3
4
5
~


6
0

17

13



~
C2 3C3
  2  7 6
 2 
7
6 



  0  3  2
 0
3
 2

 C2 (  6)C1

~    6 0  17  13  ~  0 0 3  1 2  1 .

 C3 2C1 


2

7
6
1
2


 0 
На этом этапе мы замечаем, что при 3  1  0 , т.е. при    13 , система
несовместна, а при    13 система является определенной, так как проводя
преобразование 311  C2 , мы получаем, что
19
 0  3  2 
 0 0

 C1 3C2 
 0 0  32 11  ~  0 0 

 C3 2C2 
0

1
2


 0  0
5 
3 1 
2  1
3 1 
2 4 
3 1 
и, следовательно, единственное решение системы имеет вид:
5
2  4 2  1 

x  
,
,
 .
 3  1 3  1 3  1 
T
Например, при   2, x  1, 0, 1T . Предлагаем читателю для вычислительного
контроля проверить, что при   2 вектор 1, 0, 1T действительно удовлетворяет исходной системе уравнений.
Наличие параметра в СЛАУ не обязательно приводит к существованию
таких его критических значений, как в примере 9, при которых СЛАУ меняет
своё качество.
Пример 10.
x1  x2  x3  3
 2 x1  x2  x3  2 .
x1  x2  x3  
◄
 1 1  3
 1
1 
3  C1 C2  0
3  4

 C2 C1
 ( 1)C2 
 C1 3C3

2
1

1

2
~


2
0
1
~


2
0

1



 1 
 ~
C3 C1
C2  2C3

 2 C3





3 
1

1
1

0

2
0


3
0

0
2 





 0 0  321 


 0 0 2    .

3 
2 
 0  0
Таким образом, рассматриваемая СЛАУ при любых значениях  является
определенной, а её единственное решение имеет вид:
T
3   3  1 

x   2  ,
,
 .►
2
2 

Пример 11.
x1  x2  x3  3
2 x1  2 x2  2 x3  3 .
x1  x2  x3  2
20
Предлагаем читателю убедиться, что при любых значениях  данная
СЛАУ является несовместной.
Если параметр входит в несколько коэффициентов рассматриваемой системы уравнений, прежде чем применять метод Гаусса, часто целесообразно,
провести вспомогательные преобразования, учитывающие структуру матрицы
СЛАУ.
Пример 12.
x1  x2  x3  3
x1  x2  x3  3 .
x1  x2  x3  3
◄ Проведем следующее вспомогательное преобразование:
 1 1  3
   2   2   2 9

 C1 C2 C3 

1

1
3
~
1

1
3



.
  1 1 3
 
1
1 3 



Ясно, что при   2  0 (   2 ) система несовместна. Полагая   2 , проводим преобразование C1 (  2) , после этого следуем стандартному алгоритму
метода Гаусса.
1 1  9 2 
 1 1 1  9 2 
 1

 C2 C1

3(  1)
 1  1 3  ~  0   1 0  2  .
  1 1 3  C3 C1    1
0 0 3(21) 



Из вида последней матрицы следует, что мы обнаружили еще одно критическое значение параметра   1.
Если   1 , мы можем провести элементарное преобразование
C 2   1, C3   1 .
И тогда, применяя стандартный ход метода Гаусса, получаем, что:
 1 1 1

 0 1 0
 0 0

  2  C C 
 1 3
9
0
1 1
~  0  0

 0 0
3
 2 

 2 
3
 0
  2  C C 
 1 2
6
0 
~  0 


3
 2 
 0
 2 
3
3

 2 
0
 2  ,
0
 2 
3
3

то есть СЛАУ является определенной, а её единственное решение имеет вид
T
3
3 
 3
x
,
,
(2.12)

 2  2  2
Если же   1, мы получаем СЛАУ
21
 1 1 1 3


 0 0 0 0 ,
 0 0 0 0


которая совместна и неопределенна, а её общее решение имеет вид:
T
xоб  3  C 2  C3 , C 2 , C3  , где C 2 , C3  R .
(2.13)
Таким образом, при   2 исходная СЛАУ несовместна, при   2 и
  1 она определена и имеет решение вида (2.12), а при   2 и   1 она
совместная, неопределенная и имеет общее решение вида (2.13).►
Пример 13.
x1  x2  x3  x4  
x1  x2  x3  x4   1 .
x1  x2  x3  5 x4  5
◄
1  1 1 1  
 0 0 0 2   1 1 C  0 0 0 1 21 

 C1C2 
 2 1
 C1C3
1

1
1

1

1
~
1

1
1

1

1
~
1

1
1

1

1



1 
 ~
C3 C2
C3
6
1  1 1 5 5 
0 0 0 6
 0 0 0 1
 C2 C3
6
1



 

 0 0 0 0

~   1 1 0
 0 0 0 

 1


0.
1
2
Для того, чтобы эта система была совместна, необходимо, чтобы   1  0 ,
т.е.   1. Тогда общее решение имеет вид:
x1  x2  x3
x4  1
,
или
xоб  (C2  C3 , C2 , C3 , 1)T , C2 , C3  R .►
Пример 14. При каких значениях параметра система
(5   ) x1  3x2  2 x3  0
6 x1  (4   ) x2  4 x3  0
4 x1  4 x2  (5   ) x3  0
имеет ненулевые решения?
22
◄
5  
3
  C 2 C  5  
 3 

 2 1 

~
2 0 .
 6 4 4 
 2  4
 4
 2C3 (5 )C1 8  (5   ) 2 7  3 0 

4
5






Пусть   2  0 (случай   2  0 исследуем отдельно). Сократим во втором уравнений на   2 . Получим:
 5
 11  
 3 
0 

 C1 3C2 

2
 0 ~  2
 0 ,

 8  (5   ) 2 7  3 0  C3 ( 73 )C2  2(  3) 2 0 0 




так как
8  (5  λ) 2  (5  λ)(7  3 )  8  (5  λ)(5  λ  7  3 ) 
 8  (5  λ)(  2  2 )  8  2(5  λ)(λ  1)  8  2λ 2  10  2λ  10λ 
 2λ 2  12λ  18  2(λ 2  6λ  9)  2(λ  3) 2 .
Аналогично, если   3 , то имеем систему с матрицей
11   0 
 0 0 

 2C1 (11 )C3 

~
 2  0
 0  0 .
  0 0  C2 C3   0 0 




Следовательно, если (  2)(  3)  0 , система имеет лишь тривиальные
решения. При   3 имеем систему:
2 x1  3x2  2 x3  0
6 x1  7 x2  4 x3  0 .
4 x1  4 x2  2 x3  0
Используя предыдущие выкладки, получаем:
 2  3    8 0   1 C  4 0 

 
 2 1

6

7
4
~
2

0
~
2

0

 
 
.
 4  4 2  0 0 0  0 0 0

 
 

Так как x1 является свободной переменной, то общее решение СЛАУ в
этом случае имеет вид
xоб  (C1 ,  2C1 ,  4C1 )T , C1  R .
23
При   2 имеем систему
3x1  3x2  2 x3  0
6 x1  6 x2  4 x3  0 .
4 x1  4 x2  3x3  0
Используя предыдущие выкладки, получаем
 3  3   3  3 
 0 0 

 
 C1 3C3 

 6  6 4 ~  0 0 0 ~  0 0 0 .
 4  4 3   1  0 
 1  0

 



В этом случае общее решение СЛАУ имеет вид
xоб  (C1 , C1 , 0)T , C1  R .►
Исследовать систему уравнений в зависимости от параметра  . Найти
общее и одно частное решения.
12.
13.
4 x1 4 x2

2 x1  x2
 3x 4 x
2
 1
4 x3
 ,
 x3
5 x3


 x1  x2

 x2
 x1
 2x
 x2
 1
0,
2.
14.
3x3

6 x3
3x3
 6,
 3.
9,
15.
8 x2
 x1

 x1  x2
3x  2 x
2
 1
4 x3
 0,
 x3
3x3
 0,
 0.
 3x1 2 x2

7 x1 4 x2
2 x
 1  x2
 x3
 0,
 x3
 x3
 0,
 0.
16.
 x1 2 x2

 2 x1
 x
 1  x2
 x3
3x3
 0,
 0,
2 x3
 0.
17. Найти матрицу A1 , воспользовавшись методом Гаусса, если матрица
A имеет вид:
24
 2 0 3


а)  5 2 4  ,
 1 2 1


4

6
г) 
3

3
 5 2 4


б )  2 3 1 ,
 1 2 4


5 2
1

0
1 1
,
3 4 3 

2
3 2
0

3
д) 
7

2
 1 1 0 0 


0 0 1 1

в)
,,
 1 1 1 1


 1 1 1 1
1 0 1

1 0 4
.
6 2 1

2 1 1
При решении матричных уравнений вида (1.24)-(1.26) можно использовать модификации вычислительных схем, рассмотренных в примерах 5 и 6.
Пример 15. Решить матричное уравнение XA  B , где
 2 3
 1 7 
 , B  
 .
A  
1
2
6
5




 A
◄ Составим блочную матрицу   и столбцовыми элементарными
B
преобразованиями, совершаемыми над этой матрицей, приведем блок A к виду E . Тогда на месте блока B появится решение нашего уравнения,
 2


 1

 6
 2 
 0 

3
 1 2
 1 2
 C 2  2 C1 
2
0  C  2 C   0  C C  0

~ 
~
~ 
7
 1 9  ( 1)C 2  17  9 
9






5
 6 7 
 8 7 
 7
0 


.
17 

 8 
  9 17 
 является решением рассматриваеПроверим, что матрица X  
 7  8
мого уравнения,
  9 17  2 3    1 7 


  
 . ►
 7  8  1 2   6 5 
Пример 16. Решить матричное уравнение AXB  C , где
 1  1
 2 1
 5 6
 , C  
 , B  
 .
A  
1
1

1
2
7
9






◄ Применим следующую вычислительную схему.
25
1) Вводим обозначение Y  XB , тогда AXB  AY  C .
2)Составляем блочную матрицу и приводим блок A к виду E строчными элементарными преобразованиями как при решении уравнения AY  C ,
A C  ~ E Y ;
блоком B после этого наращиваем блок Y сверху и приводим блок B столбцовыми элементарными преобразованиями к виду E как при решении уравнения XB  Y ,
B  E 
  ~  .
Y   X 
После этого на месте блока Y появится решение исходного уравнения.
Действительно,
2 1 


 1 2
   1 5 6  C2 C1   1 5 6  C1 C2 
1   C C
 ~ 

A C   

 0  12 15  ~   0 17 21 ~
7
9

1
2






 0  12 15




 0
~
4

 3

1
 2 1
  C C  0
~ 
21 
4


15 
 3
0 


.
25 

18 
Таким образом,
  4 25 
 . ►
X  
  3 18 
18. Решить матричные уравнения:
7 6
 1 2 
,
4
 0 2  2
а) 
 X 
 5 4
 3
3 7
0
б) X 

,
 3 0   0 1
 1 7  1 1
в) 
X

, ,
 8 6   4 3  1 1
 1 2 4
г)  2 3 7  X
 5 8 17 


1  0 0 1
 0 2

 

 2 6 3    0 1 0  ,
 1 3 1  1 0 0 

 

1  0 0 3 
 3 5

д) X  3 9 2    0 0 2  ,
 2 5
1  0 0 1

26
 0 0  1   1 0 0   1 0 1

 
 

е)  0  1 0  X  0  2 0    0 1 0  .
  1 0 0   0 0  3   1 0 1

 
 

Метод Гаусса пригоден также для решения матричных уравнений с необратимыми, в частности, неквадратными матричными коэффициентами.
Пример 17. Решить матричное уравнение AX  B , где
 1  1 1
 3 1
 , B  
 .
A  
 2  2 2
 1 2
◄ Заметим, что в данном случае X  M 32 (R ) . Далее следуем алгоритму,
описанному в разделе 2.5,

 1 1 3 1 C2 2C1  1 1 3 1
 ~ 


 0 0 0  5 0 .
1
2
2

2
2




A B  
Матричное уравнение не имеет решений, так как СЛАУ с расширенной матрицей
   1 1 1


 0 0 0 0


совместна, а СЛАУ с расширенной матрицей
  1 1 3


 0 0 0 5


несовместна. ►
Пример 18. Решить матричное уравнение XA  B , где
 1 1 1


A   2 3  1 ,
 1  2 2


 2  3 1
 ,
а) B  
1
1

1


 2  3 17 
 .
б) B  
1
1
1


◄ Следуя вычислительной схеме примера 15, получаем
27

 0
  0 0 
1
1
0






3 1
 2
 2  3
 0  0 
2
1
1
2

2  C C   1  1
3  C 2C   1  1 0 
 A   1  2
 ~ 
 ~ 
.
а)  B   2  3
3
1
2
2 12
1
2

5

1

5
4
C

C
C

3
C





  
 1




1 1
1
0 2
1
0 2


















Откуда следует, что в случае а) уравнение не имеет решения, так как система
уравнений с матрицей
  0 1 12
1


0
 0  0 5


 0 0 0 4  2
несовместна. Следовательно, уравнение AT X T  BT не имеет решений. Поэтому равносильное ему исходное матричное уравнение тоже не имеет решений.►
б) Повторяя вычислительные операции предшествующего примера, получаем,

 0
  0 0
1
1
0






3 1
 2
 2  3
 0  0
2
1
1
2

2  C C   1  1
3  C 2C   1  1 0 
 A   1  2
 ~ 
 ~ 
.
 
3
1
2
2 12
2

3
17
2

5
15

5
0
C

C
C

3
C
B





  
 1




1 1
1
0 0
1
0 0


















Теперь система уравнений с матрицей
 0
1 12 1


 0  1  5 0


0
0
0
0
0


 
совместна. Её общее решение X T находим поочередно: X T
шение СЛАУ с матрицей
 0
1 12 


 0   1  5  , т.е. X T


 0 0 0 0
 
1
1
есть общее ре-
  C1  12 


  C1  5  ,
 C

1


28
 
а XT
2
есть общее решение СЛАУ с матрицей
 0
1 1


 0   1 0  , т.е. X T


0
0
0
0


 
2
  C 2  1


  C2  .
 C

2


Таким образом, общее решение исходного матричного уравнения
  C  12 C1  5 C1 
 , где C1 ,C 2  R . ►
X   1

C

1
C
C

2
2
2
19. Следующие матричные уравнения решить методом Гаусса:
1
 2
 1 2
 X  
 ,
а) 
  6  3
 1 0 
 4 1
 6  3
 X  
 ,
б) 
 8 2
12  6 
2   6
4
3
  
 .
в) X 
15

10
12

8

 

Историческая справка
Метод последовательного исключения неизвестных для нахождения
решения системы линейных алгебраических уравнений был впервые описан К.
Гауссом в 1849 г. Однако, уже во II веке до н.э. в Китае был известен близкий
метод под названием «фан-чен», что переводится как «выстраивание чисел
по клеткам».
Основная литература.
1. Дыбин В.Б. 12 лекций по алгебре. Пособие для первокурсника. 2006.
Электронная форма.
2. Курош А.Г. Курс высшей алгебры. М: Наука, 1973.
3. Кострикин А.И. Введение в алгебру. Часть 1. Основы алгебры. М.: Физико- математическая литература, 2000.
Задачники и дополнительные методические материалы.
4. Фаддеев Д.К., Соминский И.С. Сборник задач по высшей алгебре. М:
Наука, 1972.
29
5. Проскуряков И.В. Сборник задач по линейной алгебре. М: Лаборатория
базовых знаний, 2001.
6. Кряквин В.Д. Линейная алгебра в задачах и упражнениях. М.: Вузовская
книга, 2006.
30
Download