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

Все в одном

.pdf
Скачиваний:
39
Добавлен:
13.02.2015
Размер:
23.06 Mб
Скачать

 

 

 

 

 

 

 

(t) t (t)

 

 

t 2

 

 

 

 

(0) 1

(t) e

 

2

 

 

 

 

 

 

 

 

 

 

1. Теорема Карлина

fi:X R, где X Rn, i=1..m, т.е. есть n-мерная вектор-функция f(x)=(f1(x)..fm(x)). Поставим задачу f(x)= (f1(x)..fm(x)) maxx X (1)

Опр: доп. решение (1) наз. альтернативой. X - множество всех альтернатив.

Rn – пространство альтернатив. f(x) Rm – критериальное пространство (множество образов множества альтернатив).

 

 

1

 

2

i

fi

fi

Введем на Rm отношение порядка “ ”. Пишем f1 f2, если

 

f 1

f 2

i

0

 

i0

 

i0

На X вв. отношение доминирования x1 x2 f(x1) f(x2)

Опр: x X называется эффективной точкой Парето, если не x X: x x .

Множество всех точек Парето обозначим П.

F= {f(x)} Rm f(П) F (f(П) – образ множества Парето)

Свертывание критериев

n

 

Рассмотрим (1). Вв. (x)= i fi (x) (2),

I – произв.

i 1

 

Есть 2 способа выбора I: a) I=1, I 0;

б) 1=1, остальные – произв.

Сформулируем задачу (x) maxx X (3) – задача с 1 критерием.

Опр: (x) – линейная свертка

После введения (x) отношение выглядит: x1 x2 (x) (x) (4)

Теорема: (Карлина) 1) X – выпукл. мн-во, все fi(X) – вогн., x0 П задачи (1) вектор 0: x0 – оптимальное решение задачи (3);

2) 0 – фикс. вектор, x0 –решение задачи (3) при этом x0 П задачи (1)

▲ x0 П. Построим H={y Rm|y=f(x)-f(x0),x X}. Обозначим Em Rm – неотрицательный октант критериального пространства. Em ={y|y 0} – с вершиной в начале координат. Докажем, что H Em может содержать только

0 . это не так, y 0, y H x:f(x) f(x0) – это противоречие x0 П H Em ={ 0 }. Обозначим K – L(H) (лин-ая

p

 

 

 

оболочка, натянутая на H), H K. y K y= i y i H , yi H, i

1, i 0 . Покажем, что K Em ={

 

}:

0

i 1

 

 

 

берем произвольную y K y= y1+(1- )y2 (для простоты). y1,y2 H x1,x2 X 0 y= (f(x1) -f(x0))+(1- )(f(x2)-f(x0)) = f(x1)+(1- )f(x2)-f(x0) т.к. вогн f( x1+(1- )x2) -f(x0) (*), ( x1+(1- )x2 X, т.к. X – выпукл.) Предп., что в K есть 0 и 0

из (*) получаем противоречие x0 – эфф.

Есть множества K – выпук. и замкн. мн-во, Em -выпук., внутр. точки их не перес. K Em содержит лишь гран.

точки можно отделить, т.е. c Rm( 0) и -скаляр, определяющий гиперплоскость (c,x)= : (c,y) (c,y1)

y Rm,y1 Em (5), т.к. 0 H K, 0 Em из (5) =0 0 (c,y1) c 0

m

Строим : ci 0 i=ci/ 0, i=1 по первому неравенству из (5) (c,y) 0 y K, но K H (c,y) 0

i 1

y H (c,f(x)-f(x0)) 0. Делим на (x) (x0) x X x0 – опт (3). Доп., что x0 – опт. для (3) (x) (x0)

x X. Доп., что x0 П x : x x0

f (x ) f (x0 ) (x ) (x0 ) , а это противоречит x0 – опт. для (3).

2. Теорема Форда-Фолкерсона о величине max потока и min пропускной способности разреза.

Опр: дан ориентированный граф G, на дугах которого веса. Этот граф называется сетью, а веса – пропускной способностью дуги (i,j). G=(X,U),(i,j) U, cij 0.

Опр: Стац. потоком сети G из источника s в сток t наз. n-мер-я ф-я f:U R+ , уд. след. усл.:

fij

f ji 0, i s, t;0 fij cij (i, j) U .

 

 

j i

j i

Смысл неравенства: на дуге значение fij не может превосходить пропускной способности дуги. fij – поток по дуге.

f(i,X)=f(i, +i) –общий поток из i, f(i,X)-f(X,i) – чистый поток через i.

Опр: Стац. потоком (2-е опр) называется такая потоковая функция, которая удовлетворяет условию сохранения потока во всех промежуточных вершинах и на дуге ее значение не превосходит пропуск способности дуги. {f(i,X)-f(X,i)=0, i s,t; 0 f c}

Опр: Величина потока: f: U R+-ст поток, через v=f(s,X)-f(X,s) (где v – чистый поток из источника s) будем обозначать величину стац. потока.

Таким образом, задача о нахождении max величины потока сводится к

v max, f (s, X ) f ( X , s) v; f (i, X ) f ( X ,i) 0,i 1..t 1f (t, X ) f ( X ,t) v,0 f c

Опр: Поток – max, если он имеет наибольшую величину из всех возможных.

G – сеть, s,t – источник и сток. Разобьем X на 2 подмножества X Y Y ,Y X ,Y X \ Y . s Y,t Y множество дуг <Y, Y > , направленных из Y в Y .

Опр: <Y,Y >={(i,j) U|i Y&j Y } называется разрезом, отделяющим источник s от стока t. Одновременно с <Y,Y > может <Y ,Y>. Если удалить все дуги разреза, то не будет дуг, соед. s и t. Пропуск. cпос-ть разреза –

 

 

 

cij

проп. cпос-ть его дуг C( Y,Y

)

 

 

 

 

 

 

 

 

 

(i, j) Y ,Y

Лемма: G-сеть, f – стационарный поток величины v разреза, отделяющий источник от стока, выполняется следующее соотношение: v f ( Y ,Y ) f ( Y ,Y ) C( Y ,Y )

f (s, X ) f (X , s) v,

f (i, X ) f (X , i) 0, i 1..t 1

▲ В системе

 

 

 

 

 

 

просуммируем только те уравнения,

f (t, X ) f (X , t) v

 

 

 

 

 

 

 

которые соответствуют вершинам i мн-ва Y: v f ( Y ,Y

) f (Y ,Y )

 

 

 

)

0

 

C ( Y ,Y

Теорема: (Форда-Фолкерсона)

В сети величина max потока = (*) пропускной способности min разреза (т.е. min пропускной способности)

f,v – max поток и его величина. Рассм. процедуру, строящую мн-во вершин Y X:

1.s Y (включаем s в Y);

2.i Y, j Y: (i,j) U&fij cij j включаем в Y;

3. i Y, j Y: (j,i) U&fji 0 j включаем в Y.

Покажем, что t Y. это не так на 2 или 3 шаге t включаем в Y

a) i t б) i t

f c

f 0

Опр: Цепь наз увеличивающей, если на прямых дугах f c, а на обратных f 0.

Продолжая эти рассуждения относительно i, мы построим цепь из s в t, на прямых дугах которой f c , на обратных f 0 эта цепь – увел-ая (но по Л: чтобы f-поток имел max величину в G не увеличивающих цепей из s в t) противоречие с f - max

t Y Y 0 Y ,Y образуют разрез, отделяющий s от t.

По 2) на дугах этого разреза f=c (если нет Y мы можем расширить).

Y ,Y : потоковая функция=0, иначе мы можем расширить по Л:

vf (Y ,Y ) f (Y ,Y ) c(Y ,Y ) построили разрез – min.

0

Для всех ост разрезов v c(Y ,Y )

Алгоритм (Форда-Фолкерсона): Вход – сеть G=(X,U), функция пропускной способности c; s, t.

1)s : , zs - не определено, s- помечена, но не просмотрена, f 0

2)Выбираем i: помечена, но не просмотрена, если выбор - невозможен goto 7

3) j i : fij cij , j непомечена j : min{ i ,cij fij }; zi : i

4) j i : f ji 0, j непомечена j : min{ i , c f ij }; zi : i

5)Если t-помечена найдена увеличивающая цепь goto 6, else i-просмотрена, goto 2

6)Изменить поток f, используя найденную увеличивающую цепь и t .

Новый поток: v=v+ , снять все пометки, s : , goto 2

7) end.

4. Теорема о разрешимости матричных игр в смешанных стратегиях.

Опр: Игра – формализованная модель, описывающая конфликтную ситуацию.

Игра с нулевой суммой – участвуют 2 игрока с противоположными интересами.

] из игроков владеет набором стратегий (набор правил, опред. действия игрока).

I1 : {s1 ,..., sn1 } 1 , I 2 : {Q1 ,..., Qn21 } 2

Опр: Розыгрыш в этой игре - из игроков независимо выбирает стратегию (не имея информации о противнике). <i,j> - розыгрыш. Выигрыш 1-го= aij . Выигрыш 2-го= - aij .

Возникает платежная матрица n1 xn2 , A (aij ) .

Опр: 1 , 2 , A - матричная игра 2х лиц с нулевой суммой. ] 1 {1,.., n1}, 2 {1,.., n2 }

Опр: Гарантированный выигрыш 1-го – выигрыш независимо от стратегии 2-го ( i это min aij ). Чистая

1 j n2

стратегия – игрок, выбрав ее, не меняет при любом розыгрыше.

Опр: Тройка p0 , q0 , v называется opt реш. матр. игры в чистых стратегиях, если aij0 ai0 j0 v ai0 j ij . opt решение i0 , j0 - седловая точка платежной матрицы А.

Опр: i* - максиминная стратегия 1-го, если гарантир. выигрыш –max.

Опр: j * - минимаксная стратегия 2-го, если проигрыш – min.

v = min ai* j

max min aij - нижняя цена игры. v max aij* min max aij - верхняя.

1 j n2

1 i n1 1 j n2

1 i n1

1 i n1 1 j n2

Утверждение: v a *

* .

Тройка < i* , j * ,v> - оптимальное решение v v

 

i j

 

 

 

 

Рассмотрим вектор распределения вероятности p ( p1 ,.., pn

) . pi

- вероятность, с которой 1-й выберет i-ю

 

 

 

1

 

n

 

 

 

 

 

1

 

n1

 

 

 

 

стратегию,

p pi 1, pi 0, i 1..n1 , для 2-го игрока аналогично получим вектор q – это множество

 

i1

 

 

 

 

смешанных стратегий. Если игроки исп. наборы стратегий p и q соответственно, то

 

 

n1

n2

 

 

выигрыш= M ( p, q) aij pi q j (6).

 

 

 

i 1

j 1

 

 

Опр: p* - максиминная стратегия 1-го, если гарант. средний выигрыш min M ( p, q) максимален.

q Q

Опр: Тройка p0 , q0 , v называется opt реш. матр. игры в смеш. стратегиях, если

M ( p, q0 ) M ( p0 , q0 ) v M ( p0 , q) p P, q Q .

Теорема: матричная игра 2х лиц с нулевой разрешима в смешанных стратегиях, при этом p* , q* , v ,

где v M ( p* , q* ) явл. opt решением и max min M ( p, q) min max M ( p, q)

P Q Q P

▲] i, j aij 0 . Если это не так, прибавим к aij одно и тоже число 0 . Седловая точка А не изм., а изм.

лишь знач-е выигрыша ( на то же самое число). Рассм. 2 з-чи ЛП:

n2

 

 

 

 

n1

 

 

 

 

 

 

xi

max

 

и

yi

min

- пара взаимно двойственных задач.

 

i 1

 

 

 

i 1

 

 

Ax 1, x 0, x R n2

 

 

AT y 1, y 0, y R n1

 

 

 

 

 

 

 

 

 

 

 

(c, x) max

(b, y) min

 

 

c : (1,...,1) R n2 , b : (1,...,1) R n1 Ax b

и AT y c

обе задачи одновр. разр. или неразр.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

y 0

 

 

 

] X-огр. (I) – разр. (II) – разр. и w= (c, x* ) (b, y* ) -opt реш. I и II

 

 

n2

n1

 

 

 

 

 

 

 

 

 

 

xi* y*j

w .

Т.к.

А

сост. из aij 0 w 0 .

Введем

v=1/w. Обозн.

p* : vy* 0,

i 1

j 1

 

 

 

 

 

 

 

 

 

 

q* : vx* 0

p* P, q* Q, т.к. p* 1, q*

1.

Рассм.

pAq* M ( p, q* ) pA(vx* )

v( p, Ax* ) v( p, b) v M ( p* , q) ,

т.к.

 

p P ( p,b) 1.

Аналогично

p* Aq M ( p* , q) vy* Aq v( AT y* , q) v(c, q) v v p* Aq* p* , q* , v - opt.

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Симплекс-метод.

Решает задачи ЛП (оптимизации лин. функции в области, задаваемой линейными ограничениями), запис. в канонической форме. НР1 –если ограничения противоречат.

 

 

 

 

n

 

 

f (x1 ,...,xn ) c j x j

max

 

 

 

 

 

j 1

f (x) (c, x)

 

 

 

 

 

 

 

n

 

bi , i 1..m,

 

 

aij x j

(I) Ax b

(II)

j 1

 

 

 

 

x

 

0, j 1..n

x 0

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Опр: Реш. СЛУ Ax=b наз. базисным, если ненулевым координатам этого решения соответствуют линейно независимые (ЛНЗ) столбцы матрицы A.

Опр: Базисное решение СЛУ Ax=b называется опорным, если все его координаты 0.

1 n

~

 

1

,..., x

k

,0,...0), k n .

 

 

 

 

 

] A=(A ,…,A ). Возьмем x

(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

k

 

~

 

~

 

 

 

1

k

– ЛНЗ, т.е. j

A

j

0

2

0

x - базисное решение Ax

b и A ,…,A

 

 

j

 

 

 

 

 

 

 

 

j 1

 

 

 

j 1

 

Для базисного решения k m=rang(A).

Опр: Базисное решение вырожденное k m. Невырожденное k=m.

Опр: Выпук. комбинацией 2х точек x1 , x 2 X наз-ся точка x( ) x1 (1 )x 2 , где [0,1] .

Опр: Множество X – выпуклое, если x1 , x 2 X x( ) x1 (1 )x2 X , [0,1]

Опр: X – выпук. x X - крайняя точка (вершина), если она не представима в виде выпук. комбинации двух других точек этого мн-ва.

Схема: сиплекс-метода. ] есть задача II. Считаем, что Amxn, r(A)=m,m n.

~

 

 

 

 

 

 

~

1

,..., x

m

,0,...0) .

Предп. x - опорное решение, все опорные реш. невырождены (задача невырождена).

x

(x

 

Базис, соотв. этому опорному реш., образуют столбцы B ( A1 ,.., Am ) M mxm . A1,..,Am - ЛНЗ B-1.

 

 

n

 

 

 

 

 

n

 

 

 

 

 

Bx B A j x j b (7), xB=(x1,…,xm). (7) умножим на B-1:

x B (B 1 A j )x j B 1b мы фактически

 

j m 1

 

 

 

 

 

j m 1

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

решили СЛУ относительно x базисного. Обозначим zj=B-1Aj, z0=B-1b x B z j x j z 0

(8) – СЛУ, эквив. Ax=b

 

 

 

 

 

 

j m 1

 

 

 

 

 

- исходная.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

~ B

 

 

 

 

 

 

Опр: Она наз. Жордановой формой. x

B

 

x j z

j

значение целевой функции

 

 

 

 

 

 

 

x

 

 

 

 

 

j m 1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

~ B

 

n

 

 

n

 

 

 

~ B

 

 

 

 

n

 

 

 

 

n

 

 

(C, x) (C

B

, x

B

)

C j

x j (C

B

 

x j z

j

) C j x j

(C

B

)

( C

B

z

j

x j

C j

x j

)

 

 

 

, x

 

 

 

, x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j m 1

 

 

 

 

 

 

 

 

j m 1

 

 

j m 1

 

 

 

 

 

 

 

j m 1

 

 

 

 

j m 1

 

 

 

 

~

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(C

B

z

j

C j

)z

j

(C, x)

 

 

j x j

(9) , где j C

B

z

j

C j

, j=m+1..n

 

 

(C, x )

 

 

 

(C, x )

 

 

 

 

 

 

 

 

 

j m 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j m 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Опр: Введенные j

называются относительными оценками замещения (они сравнивают значение целевой

функции в произвольной точке со значением в опорной точке).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8,9) запишем в виде симплекс-таблицы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xm 1

 

...

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

z

 

 

...

 

z

 

 

z

 

 

 

 

 

 

X N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1,m 1

 

 

 

 

 

 

1,n

 

 

1,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

... ... ...

 

...

или

X B

Z

 

 

Z 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

m

z

m,m 1

 

...

 

z

m,n

 

z

m,0

 

 

 

 

f

 

 

 

(C B , z0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

m 1

 

...

 

 

n

(C B , z 0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема: Все j , j 1..n , отвечающие некоторому опорному решению, 0 данное опорное решение явл. opt решением задачи максимизации.

Теорема: В симплекс-таблице столбец с оценкой 0 , в котором нет элементов 0 целевая функция неограниченна сверху в ОДР (это воспринимается как неразрешимость – НР2).

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

Теорема: Невырожденное опорное решение доставляет max целевой функции все оценки 0 .

При жордановых преобразованиях С-Т с ведущим элементом zlr переменная xl выводится из базиса, на ее место приходит переменная xr, т.е. xr xl .

Для новой С-Т надо выразить из l-того уравнения xr через остальные переменные, а в других уравнениях исключить.

Правила преобразования С-Т:

1)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Главный элемент заменяем на обратный zlr

zlr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

Остальные элементы в главной строке делим на главный z

 

zlj

 

 

 

 

 

 

 

 

 

 

 

 

 

lj

 

zlr

 

 

 

 

 

 

 

 

 

 

 

 

3)

 

 

 

 

 

 

 

zir

 

 

 

r

 

 

 

 

 

 

 

 

Аналогично в столбце дел. на главный с минусом zir

zlr

, r

zlr

 

 

 

 

 

 

 

 

 

 

 

 

4)

Для не из главных строки и столбца: z

 

 

zij zlr

zlj zir

i l, j r

ij

 

zlr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Правило выбора главного столбца: с наибольшей по модулю оценкой 0 .

Правило выбора главной строки: составляем отношение правых частей ограничений к элементам0 главного столбца. Выбираем ту строку, где отношение min.

Построение начальной С-Т:

f

(c, x) max

Ax Eu b, x 0,u 0 (*). Перешли к канон. форме. A содержит m ортов

1)

, где b 0 .

Ax b, x 0

 

строим С-Т. x=0,u=b – опор. решение, соотв. новым ограничениям.

f (c, x) max

2)

, где b 0, A ( A , E)

Ax b, x 0

 

x1,..., xm - базисные.

x1 xm

x m 1

 

 

 

A

 

b

 

 

 

x n

 

 

 

 

 

 

 

 

3) А не содержит полного набора ортов (столько же, сколько уравнений) используем метод искусственного базиса. Строим вспомогательную задачу.

опорное решение исх. задачи будет построена С-Т, определяющая опорное решение.

 

Решение не исходная задача имеет противоречащие ограничения ОДР – 0

НР1.

 

 

m

 

 

 

 

 

Вспомогательная задача:

F wi max

(2)

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

Ax w b, x 0, w 0

 

 

Вспомогательная задача всегда разрешима (ОДР 0 , т.к. x=0,w=b – допустимое нет НР1,

F 0 , т.е.

ограничение сверху нет НР2).

 

 

 

 

 

Теорема: У исх-й з-чи ограничения непротиворечивы opt целевой ф-ии вспом-й задачи =0.

 

~

~

~

~ ~

~

~ ~

▲ непротиворечивы x , удов-ее

Ax

b, x

0 можем рассм. (x, w) , где w 0 F(x, w) 0 . С

др. стороны, F 0 F *

0 max

 

 

 

 

 

F * 0 opt 2 (x* , w* ) , где w* 0 x* ОДР исх.

 

 

Следствия:

1)opt F 0 исходная задача имеет НР1.

2)opt значение F 0 по opt С-Т можно построить начальную С-Т исходной задачи.

Возможны 3 случая при переходе от вспомогательной задачи к исходной: ] A=Mmxn

1)r(A)=m, задача невырождена в opt С-Т вспомогательной задачи все искусств-ые переменные находятся среди небазисных элементов для перехода к начальной С-Т исходной задачи надо вычеркнуть столбцы, отвеч. искусств. переменным, строку оценок и найти новую строку оценок, исп. целевую функцию исх. задачи.

2)r( A) m , задача невырождена в opt С-Т вспом. задачи окажутся строки с искусств. переменными,

все числовые элементы которых = 0. Необходимо вычеркнуть эти строки, затем goto 1).

3)r(A)=m, задача вырождена в opt С-Т вспом. задачи могут оказаться строки с искусств. переменными

с правой частью=0 и элементами 0 выбираем в строке элемент 0 как главный, выполняем жорданово преобр-ние С-Т все искусств. переменные в таких строках будут преобразованы в небазисные goto 1)