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

VM_pidr

.pdf
Скачиваний:
19
Добавлен:
19.03.2015
Размер:
2.73 Mб
Скачать

Таблиця 3. Помножимо елементи третього рядка на “-5” і “-2” і додамо до відповідних елементів другого і першого рядків табл.2. За ключовий елемент візьмемо число “-30” і поділимо на нього всі елементи другого рядка і запишемо другим рядком таблиці 4.

Таблиця 4. Інші клітинки заповнимо аналогічно. Помножимо елементи другого рядка на “-7” і “11” і додамо до відповідних елементів третього і перших рядків таблиці 3.

Такі перетворення показано справа біля кожної з таблиць (числа вказують, що даний рядок на нього треба перемножити, а стрілки – до якого рядка потрібно додати відповідні елементи).

Якщо перші два рядки табл.4 поміняти місцями, то по головній діагоналі одержаної матриці-таблиці є три ненульові елементи:

1

0

0

7

 

 

 

2

 

 

2

 

 

5

 

 

 

 

 

 

15

 

5

 

6

 

 

 

 

 

10

 

 

 

 

 

0

1

0

7

 

 

 

7

 

 

2

 

7

 

 

 

 

 

 

15

 

 

5

 

 

6

 

 

 

 

 

10

 

 

 

 

 

0

0

1

9

 

 

1

 

1

 

13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

15

 

5

 

6

 

 

 

 

 

 

 

 

 

 

 

Це означає, що ранги основної і розширеної матриць однакові і дорівнюють 3. Тобто вихідна система 3 лінійних рівнянь з 5 невідомими сумісна і має безліч розв’язків (ранг менший за кількість невідомих).

Запишемо отриману систему базисних рівнянь, виходячи із таблиці 4.

x1

x2x3 +

7

 

 

x4

+

2

 

 

x5

=

2

,

10

 

 

 

 

 

 

 

 

 

15

 

 

 

5

 

 

 

7

 

 

x4

+

7

 

 

x5

=

2

 

,

10

 

 

 

 

 

 

 

 

 

 

15

 

 

 

5

 

 

 

9

 

 

x4

+

1

 

 

x5

=

1

.

10

 

 

 

 

 

 

 

 

 

15

 

 

 

5

 

 

 

Маємо три базисні невідомі 1, х2, х3) і дві вільні (х4, х5). Друга таблиця виділяє базисну невідому х2 третя таблиця - х3 і четверта

таблиця - х1.

Звідси загальний розв’язок вихідної системи запишеться так:

x =

2

+

7

x

2

x , x

 

=

2

+

7

x

 

7

x

 

, x

 

=

1

9

x

 

1

x

 

.

1

5

10

4

15

5

2

5

10

 

4

15

 

5

 

3

5

10

 

4

15

 

5

 

71

Для перевірки правильності знайдених розв’язків, достатньо підставити ці значення х1, х2, х3 у задану систему лінійних рівнянь і отримати тотожності.

Для одержання частинних розв’язків системи, надамо вільним невідомим х4, х5 довільних значень. Наприклад, якщо х4=0, х5=0,

то: x1 =

2

, x2

=

2

, x3 =

1

, x4 = 0, x5 = 0.

 

 

5

 

5

 

5

 

 

 

 

 

 

 

Якщо х4=1, х5= −1, то розв’язок буде такий:

x1

=

37

, x2

=

47

, x3

= −

29

, x4 = 1, x5 = −1 , і т.д.

 

 

 

 

30

 

 

 

30

 

30

 

Таких частинних розв’язків можна знайти безліч.

§14. Системи лінійних однорідних алгебраїчних рівнянь

Система m лінійних алгебраїчних рівнянь з n невідомими називається однорідною, якщо вона має вигляд

a11 x1 + a12 x2

+ ... + a1n xn = 0 ,

 

+ a22 x2

+ ... + a2n xn = 0,

a21 x1

 

 

(1.8)

..................................................

 

+ am 2 x2 + ... + amn xn = 0.

am1 x1

Вона одержується із системи лінійних алгебраїчних рівнянь (1.3), якщо b1 = b2 = ... = bm = 0. Розширена матриця Ã цієї системи

одержується із основної матриці А, якщо приєднати нульовий стовпець. Значить ранг матриці А рівний рангу розширеної матриці Ã. При цьому за теоремою Кронекера-Капеллі вихідна система завжди сумісна.

Якщо ранг матриці A системи лінійних однорідних рівнянь дорівнює кількості невідомих ( r = n ), то система має єдиний нуль-

овий (тривіальний) розв’язок: x1 = x2 = ... = xn = 0. Це випливає з теореми Крамера, оскільки всі визначники j одержуються із визначника

 

a11

a12

...

a1n

 

=

a21

a22

...

a2n

(1.9)

...

... ... ...

 

 

 

an1

an2

...

ann

 

72

заміною j -го стовпця стовпцем із вільних членів b1 = 0,b2 = 0,..., bn = 0 , а тому x1 = x2 = ... = xn = 0 = 0.

Значить, система n лінійних однорідних алгебраїчних рівнянь зn невідомими має нульові розв’язки тільки тоді, коли визначник системи , складеної із коефіцієнтів, які стоять біля невідомих відмінний від нуля.

Система n лінійних однорідних алгебраїчних рівнянь з n невідомими має ненульові розв’язки тільки тоді, коли визначник системи, складеної із коефіцієнтів, які стоять біля невідомих дорівнює нулю ( = 0 ). Тобто , система має ненульові розв’язки тільки тоді,

коли r < n , де r - ранг матриці А, а n - число невідомих.

Приклад 1. Розв’язати систему лінійних однорідних рівнянь.

 

 

 

 

x1 + x2 x3 = 0,

 

 

 

 

 

= 0,

 

 

 

 

2 x1 3x2 + 2 x3

 

 

 

 

 

= 0.

 

 

 

 

3x1 + 2 x2 + x3

Розв’язування. Так як кількість рівнянь ( m = 3 ) співпадає з

кількістю невідомих( n = 3 ) , а визначник

 

 

 

1

1

1

 

 

 

 

 

 

=

 

2

3 2

 

= 3 4 + 6 9 4 2 = −16

 

 

3

2

1

 

 

 

 

 

 

 

 

 

 

 

системи відмінний від нуля, то задана система лінійних однорідних

рівнянь має тільки нульовий розв’язок:

x1 = x2 = x3 = 0.

Приклад 2. Розв’язати систему лінійних однорідних рівнянь

x1 + 3x2 x3 + 4 x4 x5 = 0 ,

 

 

 

+ x4 + x5 = 0 ,

x1 + 3x2 x3

 

+ 6 x2

2 x3

+ x4 = 0 ,

2 x1

 

+ 9 x2

3x3

+ x4 + x5 = 0.

3x1

Розв’язування. Знайдемо ранг

матриці, складеної з коефіціє-

нтів при невідомих. Для цього зведемо її до діагонального вигляду з допомогою елементарних перетворень:

73

1 3 1 4 1

1 3 1

 

4

 

1

 

 

 

 

 

 

 

 

 

 

3

 

 

 

A =

1 3

1 1 1

0 0 0

 

2

 

 

 

2 6

2 1 0

0 0 0

 

7

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 9

3 1 1

0 0 0 11 4

 

1 0 0

0

0

1 0

0

 

 

1 0 0

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

0 0 0

3 2

0 1

 

0 1 3

0 0 0

7 2

0 1

7

 

 

0 0

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0

11 2

0 1

11

 

0 0 7

 

1

0

0

 

1

0 0

1

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 1 0

 

0 1 0

0 1 0

.

 

 

0 0 4

0 0 4

 

 

 

 

 

 

 

 

0 0

4

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

0 0

 

0 0 0

 

 

 

 

 

 

 

Ранг останньої матриці, а значить, і еквівалентної їй матриці

A дорівнює

3 ( r = 3 )

і менший,

ніж число

невідомих,

а тому

вихідна система має ненульові розв’язки. Візьмемо ті рівняння заданої системи, в яких коефіцієнти при невідомих утворюють базисний мінор, наприклад:

1 4 1 М = 1 1 1 = −1 + 8 + 2 1 = 8,

2 1 0

який відмінний від нуля.

Задана система лінійних рівнянь еквівалентна такій:

x1 + 4 x4 x5 = −3x2 + x3 ,

 

+ x5 = −3x2 + x3 ,

x1 + x4

 

= −6 x2

+ 2 x3 .

2 x1 + x4

Розв’яжемо її, відносно невідомих

x1 , x4 , x5 , методом Гаусса.

Виключимо x1 в другому і третьому рівняннях

x1 + 4 x4 x5 = −3x2 + x3 ,

 

3x4 + 2 x5 = 0,

 

 

7 x4 + 2 x5 = 0.

 

Виключимо x4 в третьому рівнянні

74

 

+ 4 x4 x5 = −3x2 + x3 ,

x1

 

3x4 + 2 x5 = 0,

 

 

8

x5 = 0.

 

 

3

 

 

 

Із останнього рівняння знаходимо, що x5 = 0 , а значить, і x4 = 0 (із другого рівняння). З першого рівняння одержимо

x1 = −3x2 + x3 .

Таким чином, x1 = −3x2 + x3 ; x4 = x5 = 0.

При довільних значеннях вільних невідомих x2 та x3 одер-

жимо відповідні значення базисних невідомих. Наприклад, один із часткових розв’язків такий:

x1 = −2; x2 = 1; x3 = 1; x4 = 0; x5 = 0.

ТЕОРЕМА. Якщо визначник (1.9) системи n лінійних однорідних алгебраїчних рівнянь з n невідомими дорівнює нулю, а серед алгебраїчних доповнень Аij ( і = 1,2,...,n; j = 1,2,...,n)

елементів i го рядка є ненульові, то ця система має ненульо-

вий розв’язок: xj=Аij·t.

 

 

 

 

 

 

(1.10)

Тут t- деякий параметр.

 

 

 

 

 

 

 

Доведення. Підставивши розв’язок (1.10) в систему рівнянь

(1.8) при m = n , одержимо:

 

 

 

ak 1 x1 + ak 2 x2 + ...+ akn xn = ( ak 1 Ai1 + ak 2 Ai 2 + ...+ akn Ain )t = 0 t = 0,

при k i (за теоремою анулювання);

 

 

ak 1 x1 + ak 2 x2 + ...+ akn xn

= ( ai 1 Ai 1 + ai 2 Ai 2 + ...+ ain Ain )t = t = 0,

при k = i (за умовою теореми).

 

 

 

Приклад 3. Розв’язати систему лінійних однорідних рівнянь

 

 

 

 

x1 + x2 + x3 = 0,

 

 

 

 

 

+ 6 x2

+ 5 x3

= 0,

 

 

 

 

3x1

 

 

 

 

 

+ 4 x2

+ 3x3

= 0.

 

 

 

 

x1

Розв'язування. Визначник

 

 

 

=

 

1

1

1

 

= 18 + 12 + 5 6 9 20 = 0.

 

 

 

3 6 5

 

 

 

1

4

3

 

 

 

 

 

Знайдемо алгебраїчні доповнення до елементів, наприклад, першого рядка:

75

A = ( 1 )2

6

5

= −2, A = ( 1 )3

3

5

= −4, A = ( 1 )4

3

6

= 6.

11

4

3

12

1

3

13

1

4

 

 

 

 

 

Задана система лінійних однорідних алгебраїчних рівнянь має розв’язок

x1 = −2t ; x2 = −4t ; x3 = 6t , який залежить від параметра t .

§15. Деякі економічні задачі

Задача міжгалузевого балансу

Вдеяких задачах макроекономіки ставиться питання про ефективне ведення багатогалузевого господарства. Тут кожна галузь є і виробником , і споживачем деякої продукції (як своєї, так і продукції, виробленої іншими галузями).

Однак, з економічної точки зору, міжгалузевий баланс є більш ефективним у вартісному виразі. При цьому об’єднання окремих галузей у підгрупи полегшує складання балансів продукції.

Введемо такі позначення:

xi - загальна вартість продукції, виробленої в і-ій галузі (план валового випуску продукції) (i=1,2,…,n);

xij - вартість продукції i ої галузі, необхідної для випуску

продукції j го підрозділу ( j = 1,2,...,n ) ;

yi - вартість продукції i ої галузі, призначеної для реалізації (кін-

цевий продукт).

 

 

i ої галузі,

 

 

 

 

Прямі витрати одиниць

які використовуються

для випуску одиниці виробу продукції

j ої галузі, а також кінце-

вий продукт задані таблицею:

 

 

 

Кінцевий

 

 

Вартість

 

Прямі витрати

 

 

 

 

продукції

1

2

n

 

продукт

 

 

x1

x11

x12

X1n

 

y1

 

 

x2

x21

x22

X2n

 

y2

 

 

 

 

 

xn

xn1

xn2

xnn

 

yn

 

Зв’язок між цими величинами запишемо у вигляді системи рівнянь:

76

x1 = x11 + x12 + ...+ x1n + y1 ,

 

= x21 + x22 + + x2n + y2 ,

x2

 

 

.............................................

 

= xn1 + xn2 + + xnn + yn .

xn

Рівняння цієї системи називаються балансовими. Позначимо aij - вартість продукції i ої галузі , необхідної

для випуску одиниці продукції j ої галузі:

 

 

 

 

 

aij

=

xij

.

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

Матриця, складена із величин aij

 

a11

a12

...

a1n

 

 

 

 

 

 

називається

 

 

a22

...

 

 

A =

a21

a2n

матрицею прямих

...

...

...

...

 

 

 

 

 

 

 

витрат,

 

an1

an2

...

ann

 

 

 

 

 

 

 

 

 

а її елементи – коефіцієнтами прямих витрат.

Враховуючи, що xij = aij x j , вихідна система запишеться так:

 

x1 = a11 x1 + a12 x2 + ... + a2n xn + y1 ,

 

 

= a21 x1 + a22 x2 + + a2n xn + y2 ,

 

x2

 

 

 

 

........................................................

 

 

= an1 x1 + an2 x2 + + ann xn + yn ,

або

xn

x1 ( a11 x1 + a12 x2 + ...+ a2n xn ) = y1 ,

 

 

( a21 x1 + a22 x2 + + a2n xn ) = y2 ,

 

x2

 

 

 

 

........................................................

 

 

( an1 x1 + an2 x2 + + ann xn ) = yn .

 

xn

Позначимо через

x1

 

 

y1

 

 

 

 

 

 

 

 

X = x2

 

і назвемо вектор-планом

X , а Y = y2

 

і назвемо

...

 

 

...

 

 

 

 

 

 

 

 

xn

 

yn

 

вектором кінцевих продуктів Y .

77

Попередня система запишеться у вигляді матричного рівняння

X AX = Y , або EX AX = Y , звідси ( E A )X = Y ,

де E - одинична матриця.

Позначимо E A = B , тоді система лінійних алгебраїчних рівнянь запишеться такBX = Y .

Помножимо з лівого боку обидві частини рівняння на B1 :

B1BX = B1Y . Звідси X = B1Y .

Тобто вектор-план X можна знайти, помноживши B1 на вектор кінцевих продуктів.

Матриця B1 називається матрицею повних витрат. Елементи цієї матриці включають прямі і непрямі витрати.

Задача 1. Прямі витрати трьох галузей виробництва, а також обсяги кінцевих продуктів ( у грошових одиницях) задані у таблиці:

Продукція цехів

 

Прямі витрати

Кінцевий продукт

 

1

 

2

3

 

1

0,2

 

0,3

0,1

50

2

0,4

 

0,2

0,5

80

3

0,1

 

0,3

0,6

100

Потрібно знайти:

1)матрицю повних витрат;

2)план кожної галузі;

3)виробничу програму галузей;

4)коефіцієнти непрямих витрат.

Розв’язування. Із таблиці видно, що матриця прямих витрат

 

0

,2

0,3

0,1

буде:

 

,4

0,2

0

,5

 

A = 0

.

 

0,1

0,3

0

,6

 

Позначимо через Х - вектор - план галузей виробництва, Y - ве ктор кінцевих продуктів:

x1

 

 

50

 

X = x2

, Y =

80 .

x3

100

Зв’язок між величинами, записаних в таблиці представимо у вигляді системи лінійних рівнянь:

78

x1 x2 x3

( 0,2 x1 + 0,3x2 + 0,1x3 ) = 50, ( 0,4 x1 + 0,2 x2 + 0,5 x3 ) = 80, ( 0,1x1 + 0,3x2 + 0,6 x3 ) = 100.

В матричній формі маємо : Х −AХ=Y , або (E−A)Х=Y. Позначимо E−A=B. Система лінійних алгебраїчних рівнянь

запишеться в матричній формі:BX=Y. Звідси X= B -1Y. В нашій задачі

1 0

0 0,2 0,3 0,1 0,8

0,3 0,1

E A = 0

1 0

0,4

0,2

0,5

= − 0,4

0,8

0,5 = B.

0

0

1

0,1

0,3

0,6

0,1 0,3

0,4

Для знаходження оберненої В-1до матриці В, обчислимо визначник:

 

 

0,8

0,3

0,1

 

 

 

 

B

 

=

 

0,4

0,8

0,5

 

= 0,256 0,015 0,012 0,008

 

 

 

 

 

0,1

0,3

0,4

 

 

 

 

0,048 0,12 = 0,053.

Тому для матриці В існує обернена В-1. Знайдемо алгебраїчні доповнення до елементів матриці В :

b

= ( 1 )2

 

 

 

 

0,8

0,5

 

 

 

 

= 0,32 0,15 = 0,17 ;

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

0,3

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

= ( 1 )3

 

 

 

 

 

 

 

 

0,4

0,5

 

 

 

 

= −( 0,16 0,05 ) = 0,11;

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

0,1

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

0,4

0,8

 

 

 

 

 

 

 

 

= 0,12

+ 0,08 = 0,2;

 

 

 

 

 

 

 

 

 

b = ( 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

 

 

 

 

 

 

 

 

 

 

 

0,1

0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

= ( 1 )3

 

 

 

 

 

0,3

0,1

 

 

 

 

= −( 0,12 0,03 ) = 0,15;

 

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

 

0,3

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

0,8

0,1

 

 

= 0,32

0,01 = 0,31;

 

 

 

 

b = ( 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

 

 

 

 

 

0,1

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

0,8

0,3

 

 

= −( 0,24 0,03 ) = 0,27 ;

 

 

 

 

b = ( 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

 

 

 

 

 

 

 

 

 

0,1

0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

79

b = ( 1 )4

 

 

 

 

 

0,3

0

,1

 

= 0,15 + 0,08 = 0,23;

 

 

 

 

31

 

 

 

0,8

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

0,8

0,1

 

 

 

= −( 0,4

0,04 ) = 0,44;

 

 

 

 

 

b = ( 1 )

 

 

 

 

 

 

 

 

 

 

 

32

 

 

 

 

 

0,4

0

,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

0,8

0

,3

 

 

= 0,64

0,12

= 0,52.

 

 

 

b = ( 1 )

 

 

 

 

 

 

 

 

 

33

 

 

 

 

 

0,4

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матриця з цих алгебраїчних доповнень буде:

 

 

0 ,17

0 ,11

0 ,2

 

 

0 ,15

0 ,31

0 ,27 ,

 

 

0 ,23

0 ,44

0 ,52

 

 

0,17

0,15

0,23

а приєднана

B

П = 0,11

0,31

0,44 .

 

 

 

0,27

 

 

 

0,2

0,52

Обернена матриця має вигляд :

 

1

0 ,17

0 ,15

0 ,23

3

,21

 

2 ,83

4 ,34

B1 =

0 ,11

0 ,31

0 ,44 2

,08

 

5 ,85

8 ,3 .

 

 

0 ,053

 

 

0 ,2

0 ,27

0 ,52

3

,77

 

5 ,09

9 ,81

 

 

 

 

 

 

 

 

 

 

 

Елементи цієї матриці B-1 - це коефіцієнти повних витрат, а

сама матриця є матрицею коефіцієнтів повних витрат.

 

2) Для знаходження плану кожної галузі, помножимо B1 на

вектор кінцевих продуктів Y :

 

 

 

 

 

 

 

 

 

 

3,21

2,83

4,34

 

50

 

821

x1

 

X = B

1

 

5,85

 

 

 

 

 

 

 

Y =

2,08

8,3

 

80

 

1402

= x2

.

 

 

3,77

5,09

9,81

100

1577

x3

 

Значить: x1 = 821; x2 = 1402; x3 = 1577.

Отже, якщо відомо обсяг кінцевої продукції (у грошових одиницях) y1 = 50; y2 = 80; y3 = 100 , то потрібно запланувати такі об-

сяги виробництва для першої галузі - 821, для другої - 1402 і для третьої - 1577.

3) Для знаходження виробничої програми кожної галузі, знайдемо добуток коефіцієнтів прямих витрат і валового випуску продукції:

x11 = a11 x1 = 0,2 821 = 164 ,2; x12 = a12 x2 = 0,3 1402 = 420 ,6;

80

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