Линейная алгебра. Векторная алгебра. Аналитическая геометрия
.pdf1) r(A + B) 6 r(A) + r(B); 2) r(A + B) > jr(A) r(B)j;
3)r(AB) 6 minfr(A); r(B)g;
4)r(AB) = r(A), åñëè B невырожденная квадратная матрица.
Åñëè ðàíã r(A) = r, то минор порядка r, отличный от 0, называется
базисным, а входящие в него строки и столбцы базисными. Заметим, что о базисных строках и столбцах можно говорить только после выбора базисного минора.
Строки матрицы можно рассматривать как n-мерные вектора. Поэто-
му понятия линейной комбинации è линейной зависимости строк
вводятся также как в x10 эти понятия вводились для векторов.
Теорема 12.2. (теорема о базисном миноре) Любая строка (любой столбец) матрицы является линейной комбинацией ее базисных строк (столбцов).
Доказательство. Пусть r(A) = r. Не умаляя общности можно считать, что отлич- ный от 0 минор M стоит в правом верхнем углу (при перестановке строк и столбцов по теореме 12.1 ранг не меняется). Пусть s и t целые числа такие, что 1 6 s 6 m,
1 6 t 6 n.
Рассмотрим определитель порядка r + 1
|
|
a21 |
a22 |
: : : a2r |
a2t |
|
|
|
a11 |
a12 |
: : : a1r |
a1t |
|
D = |
|
: : : |
: : : : : : : : : : : : |
|
||
|
|
|
|
|
|
|
|
|
ar1 |
ar2 |
: : : arr art |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
as1 |
as2 |
: : : asr ast |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Он равен 0. В самом деле, при s < r в определителе две одинаковых строки, а при
t < r два одинаковых столбца. Значит, D = 0. При s > r, t > r определитель D = 0
как минор порядка r + 1.
Разложим D по последнему столбцу:
D = a1tA1t + a2tA2t + : : : + artArt + astAst = 0:
Алгебраическое дополнение Ast = M 6= 0 и не зависит от выбора чисел s и t. Обо-
значим через si = Ait (i = 1; r). Тогда ast = 1sa1t + 2sa2t + : : : + rsart (äëÿ
Ast
всех t = 1; n). Таким образом, строка с номером s является линейной комбинацией базисных строк.
Сформулируем несколько следствий из теоремы о базисном миноре.
49
Следствие 12.3. Если число строк матрицы больше ее ранга r, то
строки матрицы линейно зависимы. Если число строк совпадает с рангом матрицы, то строки матрицы линейно независимы.
Следствие 12.4. Определитель квадратной матрицы A равен 0 тогда
и только тогда, когда одна из строк является линейной комбинацией остальных.
Доказательство. 1. Необходимость. Пусть A матрица порядка n, определитель которой jAj = 0. Тогда r(A) = r < n. После выделения базисного минора в матрице A
найдется строка, не вошедшая в этот минор. По теореме 12.2 она является линейной комбинацией базисных строк.
2. Достаточность. Пусть k-тая строка матрицы A является линейной комбинацией
строк с номерами i1; i2; : : : ; ir: aki = 1ai1j + 21ai2j + : : : + rairj. Вычтем из k-той строки i1-ю строку, умноженную на 1, i2-ю, умноженную на 2; : : : ; ir-þ, умножен- íóþ íà r. Получим akj 1ai1j 21ai2j : : : rairj = 0. Значит, определитель jAj = 0.
Следствие 12.5. Определитель квадратной матрицы A равен 0 тогда
и только тогда, когда ее строки линейно зависимы.
Так как строки матрицы, входящие в базисный минор, линейно независимы (следствие 12.3), и любая совокупность строк, число которых больше ранга матрицы, линейно зависима, то справедлива теорема
Теорема 12.6. Ранг матрицы равен числу ее линейно независимых строк.
Эта теорема может быть взята за определение ранга матрицы. Теорема о ранге матрицы играет принципиально важную роль в мат-
ричном анализе, в частности, при исследовании систем линейных уравнений.
Задания для самостоятельного решения
Задание 12.1. В матрице A укажите базисный минор, если
|
0 |
2 |
1 |
1, |
|
|
0 |
1 |
2 |
1, |
|
0 3 |
0 |
0 |
1 |
1 |
||
|
|
|
|
|
|
|
|
A = B |
0 |
0 |
1 |
0 |
C. |
|||||
à) A = |
4 |
2 |
á) |
A = |
3 |
6 |
â) |
5 |
4 |
7 |
7 |
|||||||
|
B |
3 |
1 |
C |
|
|
B |
5 |
10 |
C |
|
B |
0 |
0 |
|
|
C |
|
|
@ |
|
|
A |
|
|
@ |
|
|
A |
|
B |
0 |
0 |
C |
|||
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
A |
50
Задание 12.2. Найдите ранг матрицы A, åñëè
0 1 0 1
|
B |
1 |
2 |
3 |
|
1 |
C |
|
à) A = |
B |
3 |
1 |
2 |
|
4 |
C |
; |
|
5 |
4 |
7 |
|
7 |
|
||
|
B |
|
|
|
|
|
C |
|
|
B |
|
|
|
|
|
C |
|
|
@ |
1 |
5 |
8 |
|
2 |
A |
|
|
0 |
3 |
4 |
11 |
|
9 |
1 |
|
|
B |
1 |
3 |
5 |
|
1 |
C; |
|
â) A = |
5 |
11 |
27 |
|
17 |
|||
|
B |
5 |
|
|
1 |
|
11 |
C |
|
B |
2 |
|
|
C |
|||
|
B |
|
|
|
|
|
C |
|
|
@ |
|
|
|
|
|
|
A |
|
B |
2 |
|
1 1 2 |
C |
|
|
|
|
|||
á) A = |
B |
3 |
|
2 |
5 |
1 |
C; |
|
|
|
|
|
5 |
|
3 |
13 |
5 |
|
|
|
|
||||
|
B |
|
|
|
|
|
|
C |
|
|
|
|
|
B |
|
|
|
|
|
|
C |
|
|
|
|
|
@ |
6 |
|
4 |
19 |
4 |
A |
|
|
|
|
|
|
0 |
2 |
|
1 |
1 |
2 |
0 |
1 |
||||
|
B |
1 |
|
1 |
2 |
3 |
4 |
C |
||||
ã) A = |
|
1 |
2 |
1 |
1 |
3 |
||||||
|
B |
|
|
|
5 |
|
8 |
|
5 |
|
12 |
C |
|
B |
1 |
|
|
|
|
C |
|||||
|
B |
3 |
|
7 |
|
|
|
C |
||||
|
B |
|
8 |
9 |
13 |
C |
||||||
|
B |
|
|
|
|
|
|
|
|
|
|
C |
|
@ |
|
|
|
|
|
|
|
|
|
|
A |
Задание 12.3. Найдите ранги матриц
|
|
|
0 3 |
5 |
6 |
4 |
1 |
|
|
|
0 3 |
5 |
6 |
7 |
1 |
||||
|
|
|
B |
1 |
2 |
4 |
|
4 |
C |
|
|
|
B |
1 |
2 |
4 |
3 |
C. |
|
|
1 |
|
4 5 |
|
2 8 |
|
2 |
|
4 5 |
4 9 |
|||||||||
A |
|
= |
B |
|
|
|
|
|
C |
è A |
|
= |
B |
|
|
|
|
C |
|
|
|
|
B |
3 |
8 |
24 |
|
28 |
C |
|
|
|
B |
3 |
8 |
26 |
10 |
C |
|
|
|
|
B |
|
|
|
|
|
C |
|
|
|
B |
|
|
|
|
C |
|
|
|
|
@ |
|
|
|
|
|
|
A |
|
|
|
@ |
|
|
|
|
A |
Задание 12.4. При каких значениях m è n ранг матрицы A равен 2,
åñëè |
0 |
2 |
3 |
4 |
1 |
1 |
0 |
3 |
|
|
|||||||||
à) A = B |
1 |
2 |
1 |
3 |
C; |
á) A = B |
2 |
||
5 |
7 |
m |
6 |
4 |
|||||
|
B |
1 |
3 |
|
7 |
|
C |
B |
|
|
B |
|
n |
C |
B n |
||||
|
B |
|
|
|
|
C |
B |
|
|
|
@ |
|
|
|
|
|
A |
@ |
|
Задание 12.5. При каких значениях p è q
0 1
2 |
4 |
1 |
1 |
3 |
1 |
2 |
C. |
m |
7 |
4 |
|
4 |
6 |
|
C |
3 |
C |
||
|
|
|
C |
|
|
|
A |
ранг матрицы
1 2 1 3 4
|
B |
1 |
|
3 |
7 |
4 |
p |
C равен 2? |
|
A = |
B |
3 |
1 |
5 |
2 |
8 |
C |
||
|
7 |
|
|
q |
|
|
12 |
||
|
B |
4 |
7 |
|
C |
||||
|
B |
|
|
|
|
|
|
C |
|
|
@ |
|
|
|
|
|
|
|
A |
Задание 12.6. При каких значениях a, b è c ранг матрицы
|
0 |
4 |
5 |
3 |
4 |
|
|
3 |
1 |
|
|
B |
3 |
2 |
1 |
1 |
5 |
C равен 2? |
|||
|
2 |
|
1 |
1 |
2 |
a |
||||
A = |
B |
|
|
|
|
|
|
|
C |
|
|
B b |
8 |
5 |
c |
|
11 |
C |
|||
|
B |
|
|
|
|
|
|
|
C |
|
|
@ |
|
|
|
|
|
|
|
|
A |
Задание 12.7. Найдите ранг матрицы A в зависимости от значений k,
51
|
0 |
4 |
|
|
|
1 |
|
3 |
0 |
1 |
|
|
B |
1 |
|
|
2 |
|
1 |
1 |
C. |
||
|
5 |
|
|
1 |
|
k + 1 |
1 |
||||
åñëè A = |
B |
|
|
|
|
|
|
|
|
C |
|
|
3 |
k |
2 |
|
|
4 |
4 |
|
1 |
||
|
B |
|
|
|
C |
||||||
|
B |
|
|
|
|
|
|
C |
|||
|
@ |
|
|
|
|
|
|
|
|
|
A |
13. Теорема Кронекера Капелли.
Рассмотрим систему m линейных уравнений с n неизвестными
8
>
>
>
>
<
a11x1 + a12x2 + : : : + a1nxn a21x1 + a22x2 + : : : + a2nxn
>: : : : : : : : : : : :
>
>
>
: am1x1 + am2x2 + : : : + amnxn
= b1;
= b2;
(13:1)
:: :
=bm:
Обозначим через A = (aij) основную матрицу системы, через Ae = (aijjbi) расширенную матрицу системы.
Ответ на вопрос о разрешимости системы (13:1) дает теорема
Теорема 13.1. (Кронекера-Капелли) Система (13:1) совместна тогда и только тогда, когда ранг основной матрицы системы совпадает с рангом расширенной матрицы, то есть r(A) = r(Ae).
Доказательство. а) Необходимость. Пусть система (13:1) совместна и вектор (x1; x2; : : : ; xn) ее решение. Тогда, подставив его в систему, получим
8 a21x1 |
+ a22x2 |
+ : : : + a2nxn |
|
> |
a11x1 |
+ a12x2 |
+ : : : + a1nxn |
> |
: : : : : : |
: : : : : : |
|
> |
|
|
|
< |
|
|
|
> am1x1 |
+ am2x2 |
+ : : : + amnxn |
|
> |
|
|
|
>
:
=b1;
=b2;
:: :
=bm:
Это означает, что последний столбец матрицы Ae является линейной комбинацией остальных столбцов и его вычеркивание не меняет ранга матрицы. Следовательно, r(A) = r(Ae).
б) Достаточность. Пусть r(A) = r(Ae). Это значит, что базисный минор матрицы A
является базисным минором матрицы Ae. Столбец свободных членов, не вошедший в базисный минор, является линейной комбинацией остальных столбцов, то есть
0 b2 |
1 |
|
0 a21 |
1 |
|
0 a22 |
||||||
B |
b1 |
C |
|
B |
a11 |
C |
|
B |
a12 |
|||
: : : |
= 1 |
: : : |
+ 2 |
: : : |
||||||||
B b |
m |
C |
|
B a |
m1 |
C |
|
B a |
m2 |
|||
B |
|
C |
|
B |
|
C |
|
B |
|
|||
@ |
|
|
A |
|
@ |
|
|
A |
|
@ |
|
|
1 |
0 a2n |
1 |
|
||
C |
+ : : : + n B |
a1n |
C |
|
|
: : : |
: |
||||
C |
B a |
|
C |
|
|
C |
B |
|
mn |
C |
|
A |
@ |
|
|
A |
|
52
Таким образом вектор ( 1; 2; : : : ; n) является решением системы (13:1), то есть система (13:1) совместна.
14. Исследование систем линейных уравнений.
Рассмотрим систему m линейных уравнений с n неизвестными
8
>
>
>
>
<
a11x1 + a12x2 + : : : + a1nxn a21x1 + a22x2 + : : : + a2nxn
>: : : : : : : : : : : :
>
>
>
: am1x1 + am2x2 + : : : + amnxn
= b1;
= b2;
(14:1)
:: :
=bm:
Найдем ранги r(A) и r(Ae) основной и расширенной матриц системы. Если r(A) 6= r(Ae), то система несовместна. Пусть r(A) = r(Ae) = r, òî
есть система совместна. Будем считать, что базисный минор M расположен в левом верхнем углу. Тогда последние m r уравнений являются линейной комбинацией первых r уравнений (являются следствиями первых уравнений) и система (14:1) примет вид:
8
> a11x1 + a12x2 + : : : + a1nxn
>
>
>
< a21x1 + a22x2 + : : : + a2nxn
>: : : : : : : : : : : :
>
>
>
: ar1x1 + ar2x2 + : : : + arnxn
= b1
= b2
(14:2)
:: :
=br
Если n = r, то по теореме Крамера система имеет единственное решение, так как = M 6= 0.
Пусть n > r. Оставим в левой части системы первые r неизвестных, остальные перенесем в правую часть.
8
> a11x1 + a12x2 + : : : + a1rxr
>
>
>
< a21x1 + a22x2 + : : : + a2rxr
>: : : : : : : : : : : :
>
>
>
: ar1x1 + ar2x2 + : : : + arnxr
= b1 a1 r+1xr+1 : : : a1nxn;
= b2 a2 r+1xr+1 : : : a2nxn;
(14:3)
:: :
=bm ar r+1xr+1 : : : arnxn:
Неизвестные, коэффициенты при которых не входят в базисный минор, называются свободными. Неизвестные, коэффициенты при которых
53
входят в базисный минор, называются базисными èëè зависимыми. Очевидно, число свободных неизвестных равно n r.
Решая систему (14:3) любым известным нам способом, найдем зависимые неизвестные
|
8 |
:x: |
1: |
= |
:f:1(:xr+1;:::::: ; xn); |
(14:4) |
|
> |
|
|
|
|
|
ãäå fi линейные |
< xr = |
fr(xr+1; : : : ; xn); |
|
|||
> |
|
|
|
xr+1; xr+2; : : : ; xn . |
|
|
|
: |
|
|
|
|
|
функции от переменных По формулам (14:4) находят общее решение системы (14:3), à, çíà-
чит, и системы (14:1). Придавая свободным неизвестным произвольные
значения, будем получать частные решения системы. Сформулируем основные теоремы о числе решений системы линейных
уравнений:
Теорема 14.1. Система линейных уравнений (14:1) имеет единствен-
ное решение тогда и только тогда, когда ранги основной и расширенной матриц совпадают с числом неизвестных, то есть r(A) = r(Ae) = n (система является определенной).
Теорема 14.2. Система линейных уравнений (14:1) имеет бесконечно
много решение тогда и только тогда, когда ранги основной и расширенной матриц совпадают, но меньше числа неизвестных, то есть r(A) = r(Ae) < n (система является неопределенной).
При решении системы (14:1) наиболее часто применяют метод Гаусса,
причем преобразования выполняют не с самими уравнениями, а с расширенной матрицей системы. Достоинствами метода Гаусса являются: а) меньшая по сравнению с другими методами трудоемкость, б) возможность одновременно исследовать систему на совместность, и если система является совместной, получить ее общее решение.
|
8 4x11 |
|
2x22 |
2x33 |
+ 2x44 |
= 2; |
||
Пример 14.1. Решите систему |
> |
x |
+ x |
x |
+ x |
= 3; |
||
|
|
|
|
|
|
|
||
|
> |
2x1 + 5x2 |
|
|
3x4 = 1; |
|||
|
> |
|
|
|||||
|
> |
|
|
|
|
|
|
|
|
< |
|
|
|
x3 |
|
|
|
|
> 3x1 + 3x2 |
x4 = 1: |
>
>
>
:
54
Решение . Хотя число уравнений равно числу неизвестных m = n, но определитель
системы = 0 и по методу Крамера систему решать нельзя. |
|
|
Ae, найдем ее ранг |
|||||||||||||||||||||||
и сравним с рангом матрицы A. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Будем решать систему методом Гаусса. Преобразуем матрицу |
|
|
||||||||||||||||||||||||
A = |
0 4 |
1 |
2 1 j 3 |
1 |
|
0 0 |
9 |
2 7 |
j 5 1 |
|
||||||||||||||||
|
|
B |
1 |
2 |
1 2 |
|
j |
2 |
C |
B |
|
1 |
2 |
1 2 |
j |
2 |
C |
|
||||||||
e |
|
2 |
5 |
0 3 j 1 |
|
0 |
9 |
2 7 j 5 |
|
|||||||||||||||||
|
B |
3 |
3 |
|
1 |
|
1 |
|
|
1 |
C B |
|
0 |
9 |
2 7 |
|
|
5 |
C |
|
||||||
|
|
B |
|
|
|
|
|
|
|
|
|
C |
|
B |
|
|
|
|
|
|
|
|
|
C |
|
|
0 0 |
@ |
|
2 |
|
|
|
|
|
j |
|
A |
1 |
@ |
2 |
|
1 |
2 |
|
j |
|
|
A |
|
|||
9 |
|
7 |
j |
|
5 1 |
|
|
|
|
2 |
|
|
|
|
|
|||||||||||
|
1 |
2 |
1 |
2 |
|
2 |
|
C |
|
0 9 |
|
2 7 j 5 |
! |
|
|
|
|
|||||||||
B 0 0 |
|
0 |
0 |
|
j 0 |
|
|
|
|
|
|
|
||||||||||||||
B |
0 |
0 |
|
0 |
|
|
j |
|
|
|
C |
|
|
|
|
|
|
|
|
j |
|
: |
|
|
|
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
B |
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
|
j |
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Значит, r(A) = r(A) = 2 и по теореме Кронекера-Капелли система совместна. За |
||||||||||
свободные |
|
e |
|
x2 |
|
x4. Получим систему |
|
|||
( |
2x3 |
неизвестные примем |
|
è |
|
|
|
|
||
= 5 9x2 + 7x4: |
|
|
|
|
|
|
||||
x1 |
x3 |
= 2 + 2x2 2x4; |
|
|
|
|
|
|
|
|
|
|
|
|
( x3 |
= 2; 5 |
4; 5x2 |
|
|
||
Ее общее решение имеет вид |
+ 3; 5x4 |
: |
||||||||
|
|
|
|
|
x1 |
= 0; 5 |
2; 5x2 |
+ 1; 5x4 |
; |
|
|
|
|
|
|
|
|
|
|
|
|
Пример 14.2. Корм для птицы, составленный из четырех видов зерна, должен содержать 17 единиц вещества А, 35 единиц вещества В, 44 единицы вещества С и 26 единиц вещества D. Известно содержание единиц полезных веществ в 1 кг. зерна каждого вида и цена 1 кг. зерна каждого вида:
Âèä |
Содержание полезных веществ (ед. в 1 кг.) |
Öåíà 1 êã. (ä.å.) |
|||
|
|
|
|
|
|
|
A |
B |
C |
D |
|
|
|
|
|
|
|
1 |
1 |
3 |
4 |
2 |
5 |
2 |
2 |
4 |
5 |
3 |
7 |
3 |
1 |
1 |
1 |
1 |
2 |
4 |
1 |
3 |
4 |
2 |
5 |
|
|
|
|
|
|
Итого |
17 |
35 |
44 |
26 |
|
|
|
|
|
|
|
Найдите все возможные составы корма, обеспечивающие необходимое количество полезных веществ. Существует ли состав корма, стоимость которого равна 60 д. е.?
Решение . Обозначим через xi вес зерна i-того вида (i = 1; 2; 3; 4) и составим систему уравнений для нахождения состава корма на каждый день
8
> |
x + 2x2 + x3 + x4 |
= 17; |
3x1 1+ 4x2 + x3 + 3x4 |
= 35; |
|
> |
4x1 + 5x2 + x3 + 4x4 |
= 44; |
> |
|
|
< |
|
|
> |
2x1 + 3x2 + x3 + 2x4 |
= 26: |
> |
|
|
>
:
Решим эту систему методом Гаусса. Будем последовательно исключать неизвестные x1 è x2.
55
0 3 |
4 |
1 |
3 |
j |
35 |
1 |
0 0 |
2 |
|
2 0 |
j |
|
16 1 0 |
0 |
1 |
1 |
0 |
j |
8 |
1 |
|
||||||
B |
1 |
2 |
1 |
1 |
|
17 |
C |
B |
1 |
2 |
1 |
1 |
|
17 |
C |
B |
1 |
2 |
1 |
1 |
|
17 |
C |
|
|||
|
|
|
|
jj |
|
|
|
|
|
|
j |
|
|
|
|
|
|
|
j |
|
|
||||||
4 |
5 |
1 |
4 |
44 |
0 |
3 |
3 0 |
j 24 |
0 |
0 |
0 |
0 |
j |
0 |
|||||||||||||
B |
2 |
3 |
1 |
2 |
|
26 |
C B |
0 |
1 |
|
1 0 |
|
|
|
8 |
C B |
0 |
0 |
0 |
0 |
|
0 |
C |
|
|||
B |
|
|
|
|
|
|
C |
B |
|
|
|
|
|
|
|
|
|
C |
B |
|
|
|
|
|
|
C |
|
@ |
|
|
|
|
j |
j |
A @ |
|
|
|
|
j |
|
A @ |
|
|
|
|
j |
|
A |
|
|||||
|
|
0 |
1 |
1 |
0 |
8 |
!: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
1 |
1 |
|
17 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j
Ранги основной и расширенной матриц равны. Значит, система имеет решение.
(Решение системы
x1 = 1 + x3 x4; x2 = 8 x3:
Это решение имеет смысл, если все переменные неотрицательны. Тогда необходимо, чтобы x3 2 [0; 8], x4 2 [0; x3 + 1].
Ответим на второй вопрос примера. Для этого в систему нужно добавить еще одно уравнение:
|
8 3x1 |
+ 4x2 |
+ x3 |
3+ 3x4 |
= 35; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
> |
|
x1 + 2x2 + x + x4 |
= 17; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
4x |
|
+ 5x |
+ x |
+ 4x |
4 |
= 44; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
> |
|
1 |
|
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
2x1 + 3x2 + x3 + 2x4 = 26; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
> |
5x1 + 7x2 + 2x3 + 5x4 = 60: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
> |
|
|
эту систему, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Решая> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
|
|
|
|
|
35 1 0 |
|
|
|
|
|
|
|
16 1 |
|
|
|
|
|
|
|
|||||||
0 3 |
4 1 |
|
3 |
j |
0 |
2 |
|
2 0 |
j |
|
|
|
1 2 1 1 |
|
17 |
|
|||||||||||||
|
1 |
|
2 |
1 |
|
1 |
j |
17 |
|
|
|
|
1 |
2 |
1 |
1 |
|
17 |
|
|
0 |
|
j |
|
|
1 |
|||
B |
4 |
5 1 |
|
4 |
44 |
C |
B |
0 |
3 |
3 0 |
j 24 |
C |
|
0 1 1 0 |
8 |
||||||||||||||
2 |
3 1 |
|
2 |
j |
26 |
0 |
|
1 0 |
j |
|
8 |
|
0 0 0 0 |
j |
|
1 |
|
||||||||||||
B |
|
j |
C B |
1 |
|
j |
C |
|
|
j |
C |
||||||||||||||||||
B |
|
|
|
|
|
|
|
C B |
|
|
|
|
|
C B |
|
|
|||||||||||||
B |
|
|
|
|
|
|
|
|
C |
B |
|
|
|
|
|
|
|
|
|
C |
|
@ |
|
|
|
|
A |
||
B |
5 |
7 2 |
|
4 |
|
60 |
C B |
0 |
3 |
|
3 0 |
|
|
25 |
C |
|
|
|
|
|
|
|
|||||||
@ |
|
|
|
|
|
|
j |
|
A |
@ |
|
|
|
|
|
j |
|
|
|
A |
|
|
|
|
|
|
|
получим, что ранги основной и расширенной матриц различны. Значит, эта система решений не имеет. Нельзя составить рацион, стоимость которого равнялась бы 60 д. е. Проверьте, можно ли составить рацион, стоимость которого равнялась бы 61 д. е.
Метод Гаусса применяют при решении компьтерном систем линейных уравнений. Программистский вариант метода Гаусса имеет три отличия от математического:
1.индексы строк и столбцов матрицы начинаются с нуля;
2.недостаточно найти просто ненулевой элемент в столбце. В программировании все действия с вещественными числами производятся приближенно, поэтому можно считать, что точного равенства вещественных чи- сел вообще не бывает. Поэтому вместо проверки на равенство нулю числа
aij следует сравнивать его абсолютную величину с достаточно малень-
ким числом ". Åñëè aij < ", то следует считать элемент aij нулевым.
56
3. при обнулении элементов j-го столбца, начиная со строки i + 1,
ìû ê s-й строке, где s > i, прибавляем i-ю строку, умноженную на ко-
эффициент k = akj
aij . Такая схема работает хорошо только тогда, когда коэффициент k по абсолютной величине не превосходит единицы. В
противном случае, ошибки округления умножаются на большой коэффициент и, таким образом, растут. Математики называют это явление неустойчивостью вычислительной схемы. Если вычислительная схема неустойчива, то полученные с ее помощью результаты не имеют никакого
отношения к исходной задаче. Схема Гаусса устойчива, когда коэффициент k = akj 6 1. Поэтому при поиске разрешающего элемента в j-ì
aij
столбце необходимо найти не первый попавшийся ненулевой элемент, а
максимальный по абсолютной величине .
Задания для самостоятельного решения
Решите методом Гаусса системы уравнений
Задание 14.1. |
8 |
2x1 |
+ |
3x2 |
+ |
x3 |
|
= |
1; |
|
|
|
|
||||
|
> |
3x1 |
+ 2x2 |
+ x3 |
|
= 5; |
|
|
|
|
|||||||
|
< |
2x |
1 |
+ x |
2 |
+ 3x |
3 |
|
= 11: |
|
|
|
|
||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
||||
Задание 14.2. |
: |
2x1 + |
x2 |
|
x |
|
+ 4x4 = 1; |
||||||||||
8 |
3x1 |
+ |
2x2 |
4x33 |
+ |
7x4 |
= |
|
1; |
||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
4x |
1 |
|
4x |
2 |
+ 10x |
3 |
+ 11x |
4 |
= 4: |
||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 3x1 |
+ |
2x2 |
+ 3x3 |
|
3x4 |
= 1; |
||||||||||
Задание 14.3. |
> |
2x1 |
+ x2 |
|
3x3 |
+ 4x4 |
= 1; |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
> |
|
|
|
3x2 |
|
9x3 + 5x4 = 3; |
||||||||||
|
> x1 |
|
|
||||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> x1 + 13x2 + 6x3 + 17x4 = 2: |
||||||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: x1 + 2x2 |
+ 3x3 |
|
|
2x4 = 6; |
|
|||||||||||
Задание 14.4. |
8 |
2x1 |
|
x2 |
|
2x3 |
|
|
3x4 = 8; |
|
|||||||
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
> |
3x1 + 2x2 |
|
x3 + 2x4 = 4; |
|
||||||||||||
|
> |
|
|
||||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> 2x1 3x2 |
+ 2x3 + x4 = 8: |
|
>
>
>
:
57
|
8 3x1 |
+ 4x2 |
|
x3 |
|
|
x4 |
= 0; |
|
|||||
Задание 14.5. |
> |
2x1 |
+ x2 |
+ x3 |
+ x4 |
= 1; |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
> |
|
|
|
|
|
x3 |
|
|
x4 = 1; |
|
|||
|
> x1 + 3x2 |
|
|
|
|
|||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
5x1 3x2 + 6x3 + 3x4 = 3: |
|
|||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
2x1 |
|
4x2 + 3x3 |
|
|
6x4 = 1; |
|
||||||
|
8 x1 |
2x2 |
|
x3 |
+ x4 = 2; |
|||||||||
Задание 14.6. |
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
3x1 |
|
6x2 + 2x3 |
|
|
4x4 = 1; |
|||||||
|
> |
|
|
|
||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
4x1 8x2 2x3 4x4 = 6: |
||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
4x1 + 3x2 |
|
3x3 |
|
|
x4 = 4; |
|
||||||
|
8 |
3x1 |
|
x2 |
+ 3x3 |
|
|
2x4 = 1; |
|
|||||
Задание 14.7. |
> |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||||
|
> |
3x1 + x2 |
|
|
|
|
x4 = 0; |
|
||||||
|
> |
|
|
|
|
|
||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
2x1 3x2 + 6x3 2x4 = 3: |
||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: x1 + x2 |
+ x3 |
|
+ 2x4 = 2; |
||||||||||
Задание 14.8. |
8 |
2x1 + 4x2 + 5x3 + x4 = 5; |
||||||||||||
> |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
> |
|
|
|
x2 |
|
2x3 + 8x4 = 1; |
|||||||
|
> x1 |
|
|
|||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
6x1 + 13x2 + 16x3 + 6x4 = 17: |
||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: x1 + 2x2 + 2x3 |
|
|
|
5x4 = 3; |
|||||||||
|
8 |
2x1 + 5x2 |
|
x3 |
+ 3x4 |
= 7; |
||||||||
Задание 14.9. |
> |
|
|
|
|
|
|
|
|
|
|
|
||
|
> |
5x1 + 12x2 |
|
|
+ x4 |
= 11; |
||||||||
|
> |
|
|
|||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> x1 + x2 |
+ 7x3 18x4 = 16: |
||||||||||||
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
|
x1 |
|
x2 |
+ 2x3 |
|
|
3x4 = 5; |
|||||
|
|
8 |
2x1 |
3x2 |
+ 5x3 |
+ x4 = 1; |
||||||||
Задание 14.10. |
> |
|
|
|
|
|
|
|
|
|
||||
|
|
> |
3x1 |
|
x2 + 10x3 + 6x4 = 3; |
|||||||||
|
|
> |
|
|||||||||||
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
6x1 5x2 + 17x3 + 4x4 = 7: |
|||||||||||
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
: x1 + 2x2 |
|
x3 |
|
|
|
3x4 = 5; |
||||||
|
|
8 |
2x1 + 3x2 |
5x3 |
+ 2x4 = 3; |
|||||||||
Задание 14.11. |
> |
|
|
|
|
|
|
|
|
|
|
|||
|
|
> |
3x1 + 4x2 |
|
6x3 + 2x4 = 11; |
|||||||||
|
|
> |
|
|||||||||||
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> 4x1 + 7x2 7x3 4x4 = 13: |
>
>
>
:
58