Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика / Шилкин Лекции.PDF
Скачиваний:
33
Добавлен:
11.05.2015
Размер:
7.25 Mб
Скачать

ЛЕКЦИЯ 8

Рассмотрим однородную систему

α11x1

α x

21 1

. . .

αm1x1

+ a12 x2 +... + a1n xn = 0

 

+α22 x2 +... + a2n xn = 0

(6)

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

 

+ am2 x2 +... + amn xn = 0

 

Или в матричной форме: Αx = 0 .

Эта система всегда совместна, т.к. она имеет нулевое решение x1 = x2 = ... = xn = 0 . Возникает вопрос: при каких условиях система (6) имеет и

ненулевые (нетривиальные) решения.

По теореме Кронкера-Капелли, это будет тогда и только тогда, когда r(Α) < n . Если r(Α) = n и матрица Α квадратная, то система имеет единственное нулевое решение.

Если в квадратной системе r(Α) < n , то система имеет ненулевое решение. То, что r(Α) < n , означает, что определитель системы равен нулю, т.е. для

того, чтобы однородная система линейных уравнений Αx = 0 , состоящая из n уравнений с n неизвестными, имела ненулевое решение, необходимо и достаточно, чтобы Α = 0 .

Докажем некоторые свойства решений однородной системы (6):

1°

Если x - решение однородной системы, то и cx , c R - тоже решение

этой системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В самом деле, если Αx = 0 Α(cx) = c(Αx) = c 0 = 0 .

 

 

 

 

 

2° Если x(1) и x(2)

- решения системы Αx = 0 , то и x(1) + x(2)

- тоже реше-

ние этой системы. Действительно, т.к. Αx(1)

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

Αx

(2)

 

(1)

+ x

(2)

 

 

(1)

+ Αx

(2)

= 0 + 0

 

= 0 .

 

 

 

 

 

 

 

 

 

= 0 Α x

 

 

= Αx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из этих свойств вытекает и более общее утверждение: если x(1) , x(2) ,..., x(k )

- решение

 

однородной системы

Αx = 0 ,

то

и их

линейная

комбинация

x = c x(1)

+ c

 

x(2) +... + c

 

x

(k ) является решением этой системы, c

 

R,i =

 

.

2

k

i

1, k

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Линейно независимая система решений x(1) , x(2) ,..., x(n) называется фундаментальной, если каждое решение системы (6) является линейной комбинацией решений x(1) , x(2) ,..., x(k ) .

Приведём без доказательства следующую теорему:

1

- если r(Α) = r < n , то система (6) имеет (n r) линейно независимых решений, т.е. обладает фундаментальными системами решений.

Найдём фундаментальные системы решений для системы (6).

Пусть r(Α) = r < n и пусть для определённости базисный минор D r -го порядка расположен в верхнем углу матрицы Α:

 

 

 

a11

a12

...

a1r

 

D

 

=

a21

a22

...

a2r

0

 

 

 

 

ar1

ar 2

...

arr

 

Перенесём слагаемые с неизвестными xr+1,..., xn первых r темы (6) в правые части. Получим систему:

α

x

+ a x

2

+... + a

x

r

= −a

x

r+1

... a

x

n

 

11 1

12

1r

 

 

1,r+1

 

1n

 

α

21x1 +α22 x2 +... + a2r xr

= −a2,r+1xr+1 ... a2n xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

+ ar 2 x2

+... + arr xr

= −ar,r+1xr+1

... arn xn

αr1x1

уравнений сис-

(7)

Придадим свободным неизвестным xr+1,..., xn последовательно следующие (n r) наборов значений:

1)xr+1 =1, xr+2 = 0,..., xn = 0

2)xr+1 = 0, xr+2 =1,..., xn = 0

(n r) ) xr+1 = 0, xr+2 = 0,..., xn =1

При этом получим (n r) решений системы (6), которые находятся лю-

бым способом (например, по формулам Крамера, т.к. x(1) = (α1(1) ,α2(1) ,...,αr(1) ,1,0,...,0)

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

x(n2) = (α1(nr) ,α2(nr) ,...,αr(nr)

D 0 ).

. .

,0,0,...,1)

Эти решения образуют матрицу:

 

α1(1) ,α2(1) ,...,αr(1)

1 0 ... 0

 

 

 

α1(2) ,α2(2) ,...,αr(2)

0 1 ... 0

 

,

Χ =

 

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

 

 

 

α(nr ) ,α(nr) ,...,α(nr)

0 0 ... 1

 

 

 

 

 

 

1

2

 

r

 

 

 

причём r(Χ)= (n r) (справа стоит

Εnr ,

 

Εnr

 

=1)

и, следовательно, решения

 

 

x(1) ,..., x(nr) линейно зависимы. Это и есть ФНСР.

Можно было бы придавать свободным неизвестным xr+1,..., xn произ-

вольные наборы значений, лишь бы определитель, составленный из этих наборов, был отличен от нуля. Но если в качестве наборов свободных неизвестных выбираются строки матрицы Εnr , то это удобнее всего, и такая система фун-

даментальных решений называется ФНСР. Общее решение системы (6) имеет вид:

2

 

 

 

x = c x

(1) + c

2

x(2) +... + c

nr

x(nr) ,

где

x(1) ,..., x(nr) - какая угодно

ФСР, а

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c1,...,cnr

- произвольные постоянные.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Структура решений неоднородных систем.

 

 

 

 

Αx(0)

 

Пусть

Αx = b -

неоднородная

система (4). Вектор

x(0) ,

для

которого

= b , называется частным решением этой системы.

 

 

 

 

 

 

 

 

Теорема:

 

 

 

 

 

 

 

 

 

 

 

 

~

(0)

 

 

(0)

 

 

Всякое решение неоднородной системы (4) имеет вид:

+ x

, где x

 

 

x = x

 

 

- некоторое частное решение этой системы,

а x - общее решение соответст-

вующей (отбрасывается вектор b правых частей) однородной системы.

 

 

 

 

 

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть x(0)

- решение неоднородной системы (4), а x - решение (6). Тогда

 

(0)

 

 

 

(0)

+ Αx

= b + 0 = b , т.е.

~

 

(0)

+ x есть решение

неоднородной

Α x

 

 

+ x = Αx

 

 

 

x = x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системы (4).

 

Пусть теперь

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0)

-

 

x - произвольное решение неоднородной системы, а x

 

некоторое её частное решение. Тогда

 

~

= b , Αx

(0)

 

~

− Αx

(0)

= b b

= 0 ,

Αx

 

 

 

= b и Αx

 

 

~

= x

(0)

+ x

есть решение однородной системы Αx = 0 . Отсюда,

~

= x

(0)

+ x .

т.е. x

 

x

 

Теорема доказана.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, всякое решение неоднородной системы состоит из суммы некото-

рого её частного решения и любого решения x

соответствующей однородной

 

 

 

~

= x

(0)

+ x = x

(0)

+ c1 x

(1)

+ c2 x

(2)

+... + cnr x

(nr)

.

 

 

 

 

 

 

 

 

системы, т.е. x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти общее решение и ФСР.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 + 2x2 +3x3 + 4x4 +5x5 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

+3x

2

+ 4x

3

+5x

4

+ x

5

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 4x2

+5x3

+ x4 + 2x5 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

+3x

2

+5x

3

+12x

4

+9x

5

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4x

+5x

2

+ 6x

3

3x

4

+3x

5

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение:

3

 

 

 

 

 

1 2 3 4 5 1 2

 

 

3

4

 

 

 

5

1 2

3

4 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 3 4 5 1 0

2 3 9

0

1 2 3 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

0

0

 

 

~

 

 

 

 

 

3 4 5 1 2 ~

0

4 11 13

~ 0

5 5

 

 

 

 

 

1 3 5 12 9 0 1

 

 

2

8

 

 

 

4

0 0

0 5 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

0 10 10

 

 

 

 

 

 

 

4 5 6

 

3 3 0

6 19 17

0 0

 

 

 

 

 

 

 

 

1

2

 

3

 

 

 

4

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ 0

 

 

 

3 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть свободные неизвестные

 

x3

и

x5 .

Тогда

x4 = x5 ,

x2 = −2x3 12x5 ,

 

 

 

 

 

 

 

 

 

 

 

x3

+15x5

 

 

 

1

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = x

 

+15x

 

 

 

 

 

2x3 12x5

 

2

 

 

12

 

 

 

 

 

 

 

, т.е. x =

x

3

+ 0 x

5

= x

 

1

+ x

 

 

0

 

 

 

 

 

 

1

3

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

3

0

 

 

5

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 x3 + x5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 x3 + x5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 + x2 + x3 + x4 + x5 = 7

 

 

 

 

 

1 1 1 1 1

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2x2

+ x3 + x4

3x5 = −2

 

 

0 -1 - 2 - 2 -1

 

 

 

 

 

 

 

 

 

3x1

~

 

 

 

23

 

 

 

 

 

 

 

 

 

x2

+ 2x3

+ 2x4

+ 6x5 = 23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0 0 0

 

 

0

 

 

 

 

 

 

5x

+ 4x

 

+3x

 

+3x

 

+ −x

 

 

=12

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0 0 0

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

3

 

 

 

4

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Линейные неравенства и их применение в экономике.

Множество точек называется выпуклым, если наряду с любыми его точ-

ками x(1) и x(2) оно содержит и все точки отрезка, соединяющего эти точки. Например:

 

 

 

 

 

 

 

 

 

, но

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отрезок, соединяющий точки x(1)

и x(2) описывается соотношениями:

x = x(1) + t(x(2) x(1) )= tx(2) + (1t)x(1) , 0 t 1.

(8)

 

 

 

В самом деле, видим, что x

(1)

x = t x

(1)

x

(2)

 

(2)

r

(1)

t 1, но

 

 

 

= r

 

 

t , 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r (1) + x(1) x = r x(1) x = r r (1) .

4

Подставляя, получим:

r

r

(1)

 

(2)

r

(1)

 

 

= r

 

 

t ,

 

 

 

 

 

 

 

 

откуда r = r (1) +tr (2) tr (1) = tr (2)

+ (1t)r (1)

- параметрическое уравнение отрез-

ка, соединяющего точки x(1) и x(2) . Учитывая, что координаты точки в пространстве и её радиус-вектора совпадают, получаем (8).

Если точки x(1) и x(2) находятся на плоскости, то в координатной форме

(8) примет вид:

x

 

= tx(2)

+ (1t)x(1)

 

1

1

1

(9)

 

 

= tx(2)

+ (1t)x(1)

x

2

, 0 t 1

 

2

2

 

В пространстве R3 добавится ещё одно такое неравенство на третью ко-

ординату x3 , в пространстве Rn

 

таких равенств будет n.

Если в (8) или (9) взять t > 0 , то получим луч, исходящий из точки x(1) в точку x(2) .

Утверждение:

Пересечение любого числа выпуклых множеств также есть выпуклое множество.

Доказательство:

Рассмотрим любые две точки x(1) и x(2) пересечения. Эти точки принадлежат каждому из пересекающихся выпуклых множеств. Отсюда следует, что и весь отрезок, соединяющий эти точки, принадлежит пересечению, т.е. и пересечение также выпуклое множество.

Утверждение доказано.

Выпуклое множество, все границы которого линейны, называется выпуклым многогранным телом или выпуклым многогранником.

Множество точек n-мерного пространства, удовлетворяющих уравнению a1x1 + a2 x2 +... + an xn = c , (10)

называется гиперплоскостью.

5

 

 

В случае,

когда

n = 2 , имеем

a1x1 + a2 x2 = c

- прямая. При n = 3 -

a1x1 + a2 x2 + a3 x3 = c - плоскость. При n > 3 имеем гиперплоскость.

 

 

Если

рассмотреть

прямую

a1x1 + a2 x2 = c ,

то

она разбивает плоскость

x o x

2

на

две

полуплоскости

p+

и p,

описываемые неравенствами

1

 

 

 

 

 

 

 

 

a1x1 + a2 x2 > c и a1x1 + a2 x2 < c соответственно. Общей границей этих полуплоскостей является множество точек прямой.

Аналогично,

всякая

плоскость в R3 ,

заданная общим уравнением

Ax + By + Cz + D = 0

разбивает это пространство на два полупространства: Η+ ,

описываемое неравенством

Ax + By + Cz + D > 0 ,

и Η, описываемое неравенст-

вом Ax + By + Cz + D < 0 , причём эти полупространства есть выпуклые множества. Для доказательства достаточно взять две точки, провести отрезок между ними и подставить в соответствующее неравенство.

По аналогии определяются и полупространства в R4 .

Η+ : a1x1 + a2 x2 +... + an xn > c ,

Η: a1x1 + a2 x2 +... + an xn < c .

Гиперплоскость a1x1 + a2 x2 +... + an xn = c , являясь общей границей полу-

пространств Η+ и Η, отделяет одно от другого. Системой линейных неравенств называется система:

α

x

+ a

 

x

2

+... + a

 

 

x

n

 

b

 

 

11

 

1

12

 

 

 

1n

 

 

 

 

1

 

α

21

x

+α

22

x

2

+... + a

2r

x

n

b

(11)

 

 

1

 

 

 

 

 

 

 

 

2

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

α

m1

x

+ a

m2

x

2

+... + a

mn

x

n

b

 

1

 

 

 

 

 

 

 

 

m

Каждое неравенство этой системы определяет полупространство, которое является выпуклым множеством. Решением же системы (11) является пере-

сечение этих полупространств, т.е. выпуклым множеством в пространстве Rn , ограниченном гиперплоскостями αi1x1 + ai2 x2 +... + ain xn = bi , i =1, m .

Вектор a = (a1, a2 ,..., an ) Rn называется градиентом гиперплоскости (10) a1x1 + a2 x2 +... + an xn = c .

Покажем, что в направлении вектора a (вектора нормали плоскости), значение выражения a1x1 + a2 x2 +... + an xn возрастает.

В самом деле, уравнение луча, исходящее из произвольной точки x0 ги-

перплоскости (10)

 

(a, x)= c имеет вид

 

 

x = x0 + ta

 

0

 

=

 

 

0

 

+t(a, a)= c +t

 

a

 

2

. Так как

 

 

 

(a, x)= a, x

 

+ta

a, x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дует, что с ростом t

значение (a, x)растёт.

, t 0 . Так как

 

0

= c , то

a, x

 

 

 

 

 

a 0 , то a > 0 . Отсюда и сле-

Системы линейных неравенств используются при построении линейной модели производства, которая описывает производство продукции и расход ресурсов в единицу времени. Система неравенств (11) описывает модель произ-

6

водства и видов продукции x1, x2 ,..., xn с использованием m видов производственных ресурсов, имеющих запасы ресурса i в количестве bi , причём требуется aij ресурса i для производства единицы продукции x j .

Пример:

Предприятие имеет возможность выпускать два вида продукции Π1 и Π2 , имея в расположении три вида ресурсов Ρ1 , Ρ2 и Ρ3 . Нормы расхода ресурсов и

запасы сырья заданы в таблице. Составить оптимальный план выпуска продукции.

Π

1

Π2

 

Запасы

 

 

 

 

сырья

Ρ

 

1

3

21

Ρ1

 

Ρ32

 

3

4

33

 

 

3

1

24

Доход

5

2

 

Составим математическую модель задачи. Пусть выпускается x1 единиц продукции Π1 и x2 единиц продукции Π2 . Тогда затраты ресурсов не должны превосходить запасов сырья.

x

+3x

 

21

 

 

1

 

 

2

 

 

 

 

3x1

+ 4x2 33

 

 

 

+ x2

24

 

3x1

 

x

0,

 

x

2

0

 

1

 

 

 

 

 

 

При этом мы должны получить максимальный доход: f = 5x1 + 2x2 max

Для решения этой задачи:

1.строим многоугольник членов;

2.строим вектор-градиент максимизируемой функции;

3.перпендикулярно вектору-градиенту через начало координат проводим прямую (гиперплоскость);

4.параллельно перемещая её (гиперплоскость) в направлении градиента, получим оптимальный план на множестве дополнительных планов.

Строим прямые:

x1 +3x2 = 21

(I.)

3x1 + 4x2 = 33

(I.I.)

3x1 + x2 = 24

(Ι.Ι.Ι.)

x1 = 0, x2 = 0

 

Определяем полуплоскости, задаваемые неравенствами (просто подставляем точку (0,0) в каждое).

7

Видим, что максимальное значение функции f = 5x1 + 2x2 достигается в крайней точке множества допустимых планов x1 = 7 , x2 = 3.

fmax = 5 7 + 2 3 = 35 + 6 = 41, т.е. необходимо выпускать 7 единиц продукции Π1 и 3 единицы продукции Π2 . При этом остаются неиспользованными 5

единиц ресурса Ρ1 . Ресурсы Ρ2 и Ρ3 используются целиком.

В моделях с большим количеством продукции и ресурсов используется

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

8