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

книги / Математические методы принятия решений

..pdf
Скачиваний:
2
Добавлен:
13.11.2023
Размер:
22.94 Mб
Скачать

2 ед. Из табл. 4.1 для w = 46; 24; 2 ед. находим f\(w) = 96; 96; 0 ед. Суммируем соответствующие значения: 0 + 96 = 96 (ед.), 85 + 96 = = 181 (ед.), 170 + 0 = 170 (ед.), и выбираем наибольшее из них: 181 ед. Оно получено для хг = 1. Для w = 46 заносим в табл. 4.2

£2 = 1, Д М =181. Оптимальные

результаты

вычисления / 2(ги)

и Х2 приведены в табл. 4.2.

 

 

Таблица 4.2

 

 

 

 

 

Второй шаг оптимизации

 

W

hbo)

Х2

W

Л(«0

Х2

0-21

0

0

48-65

192

0

22-23

85

1

66-67

255

3

24-43

96

0

68-69

266

2

44-45

170

2

70-71

277

1

46-47

181

1

72-87

288

0

Из табл. 4.2 следует, что при грузоподъемности транспортно­ го средства до 21 ед. ничего в него погрузить нельзя, при грузо­ подъемности 2 2 ... 23 ед. можно погрузить только один предмет второго вида, при грузоподъемности 24 ... 43 ед. — либо один пред­ мет первого вида, либо один предмет второго вида. Максимальной стоимость груза будет при погрузке предмета первого вида. При грузоподъемности 44 ... 45 ед. можно погрузить либо один предмет первого вида, либо два предмета второго вида. Большая стоимость груза будет в последнем случае / 2(11;)= 170. При грузоподъемно­ сти 4 6 ... 47 ед. можно погрузить один предмет первого вида, два предмета второго вида или по одному предмету первого и второ­ го видов. Стоимость груза в последнем варианте — максимальная. Далее аналогично.

Ш аг 3. Проведем оптимизацию целевой функции при условии, что самолет загружают предметами первых трех видов. Требуется максимизировать по £3 функцию

/ 1,2,3 = яз^з + /г(и> — хгРгУ, /з(го) = max / 1,2,3, 0 ^ £3 ^

.

Задаем значение w и для каждого такого значения получаем максимум / 1,2,3 по £ 3. Значения fi(w) берут в табл. 4.2. Например, пусть ш = 38 ед. Возможное значение £ 3 = 0, 1,2 ед., и соответ­ ствующая стоимость предметов третьего вида равна 0; 50; 100 ед.

На предметы первого и второго видов остается соответственно вес 38; 22; 6 ед. По табл. 4.2 находим fi{w) = 96; 85; 0 ед. Суммар­ ная стоимость будет равна 96; 135; 100 ед. Максимальное значение стоимости при w = 38 ед. равно 135 ед. при хз = 1. Далее анало­ гично. Результаты расчетов помещены в табл. 4.3.

Таблица 4.3

 

Третий шаг оптимизации

 

W

h(w)

Хъ

W

h(w)

XZ

0-15

0

0

44-45

170

0

16-21

50

1

46-47

181

0

22-23

85

0

48-63

192

0

24-31

96

0

64-69

242

1

32-37

100

2

70-71

277

0

38-39

135

1

72-87

288

0

40-43

146

1

 

 

 

Ш аг 4. Получим оптимальный вес груза w *9 который может принять самолет.

Зададим значение w , которое может быть занято четвертым гру­ зом, и вычислим Х4, / 4(w) (табл. 4.4).

Таблица 4.4

Четвертый шаг оптимизации

w

M w ) £4

w

M w )

X4

0-9

0

0

46-47

181

0

10-15

20

1

48-57

192

0

16-21

50

0

58-63

212

1

22-23

85

0

64-69

242

0

24-33

96

0

70-71

277

0

34-37

116

1

72-81

288

0

38-39

135

0

82-87

308

1

40-45

146

0

 

 

 

Второй этап решения задачи состоит в нахождении оптимально­ го решения. Максимальное значение f^(w) определяет соответству­ ющее значение аргумента Х4 количество груза четвертого вида, который берет самолет (xj). При w = 83 имеем

U(w*) = max f 4(w) = 308, х*Л= 1.

7 — 4077

Вес, занятый предметами четвертого вида, будет х^Рц = 10 ед. На остальные три вида груза остается го —10 = 73 (ед.). Входим с этим значением в табл. 4.3 и получаем х* = 0, а затем х* = 0,

х \ = Ъ.

Строго говоря, нет необходимости заполнять всю табл. 4.4, до­ статочно было заполнить только одну строку для w = 83.

Из табл. 4.1-4.4 видим, что они содержат «заготовки» для опти­ мальной загрузки самолета любой грузоподъемности до w = 87 ед. указанными предметами. Таким образом, мы решили не только по­ ставленную задачу, но и ряд родственных задач. При этом в памя­ ти компьютера необходимо держать табл. 4.1-4.4: столбцы f\(w ), h ( w\ hiw ), / 4(w) на протяжении следующего шага, а столбцы X], Х2, 23, Х4 и соответствующие им w — на протяжении всего решения. Для сложных задач последнее обстоятельство имеет существенное значение и ограничивает возможности решения задач рассмотрен­ ным методом.

В экономических приложениях это задача о распределении вы­ деленных средств S между предприятиями для обеспечения макси­ мального количества выпускаемой продукции.

§ 4.3. Задача о вкладе средств в производство

Задача. Планируется организация деятельности двух произ­ водств (например, швейного цеха и цеха по ремонту обуви) сроком на 5 лет = 5). Функции вклада средств в соответствующее про­ изводство имеют вид ср(ж) = 0,75а:, ф(у) = 0,3у, т. е. если средства z вкладывать только в первую отрасль, то по годам средств будет вкладываться z • 0,75; z • 0,752; z • 0,753; ... Соответственно, вклады­ вая средства z только во вторую отрасль, получаем z • 0,3; z • 0,32; г - 0,33; ...

Функции дохода как функции от объема вкладываемых средств х и у имеют вид

f ( x ) = l - e - x, g(y) = 1 е~2у.

Требуется распределить ресурсы z = 2 ед. между отраслями (це­ хами) по годам так, чтобы получить максимальный доход.

Р е ш е н и е . Решение задачи не удается получить в аналити­ ческом виде, решим ее графически. Пусть к началу пятого года

количество средств равно Z4 и на пятый год выделено средств пер­

вому производству

тогда второму производству будет выделено

У5 = z 4 ^ 5- Чтобы найти условное оптимальное управление на пя­ том шаге xs(z4), нужно для каждого 24 найти максимум функции

/ 5 = 1 - e~Xs + 1 - e~2(z4~X5) = 2 - (e~Xi + e~2(Z4~Xs)).

При фиксированном 24 максимум функции /5 достигается либо при х$ = 0, либо внутри отрезка [0 ,24]. Максимум внутри отрезка будет в точке х *, которую находим из условия

д / 5 _ Л

р - 1 5 _ /) f, - 2 ( Z 4 - X 5 ) _ Л

dxs ~

 

“ ’

=

J (224 - Ь

2),

х 5 ^ 0.

Из условия x s ^ O следует 24 ^

(In 2)/2, т. е. при 24 ^ (In 2)/2 я»

« 0,347 максимум достигается в точке

 

£5(24) = \ (224 -

In 2).

При 24 ^ (In 2)/2 максимум будет в точке 2:5(24) = 0. Итак, услов­ ное оптимальное управление на пятом шаге можно записать в сле­ дующем виде:

 

при

1п 2

 

 

 

(4.2)

(^4)— 4 1 ..

. . .

In 2

- (224 -

In 2) при

24 > — .

 

Построим графики (рис. 4.1) зависимостей х$ = 2:5(24) и f* = = / 5*(24). При построении зависимости f*(z^) пользуемся соотно­

шением (4.2).

 

 

 

 

Перейдем

к оптимизации на

 

четвертом

шаге.

Запас

средств

 

после третьего

шага

обозна­

 

чим 23. Определим, в каких

 

пределах

может

изменяться 23.

 

Наибольшее

значение

23 будет

 

достигнуто,

если

все

средства

 

вложить в первую отрасль:

4.1. Графики зависимостей

 

 

 

 

Рис.

23тах = 2 - 0,753 = 2 • 0,753 = 0,844.

/* и х* от аргумента Z4

Наименьший запас средств соответствует случаю, когда все средства на трех шагах вложены во вторую отрасль:

23 min = 2 • 0,33 = 0,054.

Будем рассматривать теперь набор значений 2:3 от 0,1 до 0,8 с шагом 0,1. Для каждого значения 2:3 найдем условное оптималь­ ное управление £4(2:3) на четвертом шаге к условный оптимальный доход /4 (2:3) на двух последних шагах.

Для этого построим серию кривых, показывающих выигрыш на двух последних шагах при любом управлении на четвертом

и при

оптимальном на пятом:

 

 

/ 4,5

= / 4(23, Х 4) + / 5* 0 4) = / 4( 23, Х 4)

+ / *(0,75æ4 + 0,3(23 -

Х 4» ,

 

/ 4(23, х4) = 2 - (е~х<+

 

значение /* (24) берем из рис. 4.1.

 

 

Значения / 4,5(^4) приведены на

рис. 4.2. Для каждого

значе­

ния 23 максимальная ордината определяет условный максимальный выигрыш / 4*5(23) на двух последних шагах, а абсцисса —условное оптимальное управление 2:4(23).

На рис. 4.3 построены графики зависимостей / 4*5(23) и 2:4(23).

Рис. 4.2. Графики завиРис. 4.3. Графики зависимостей симости / 4,5 от аргумен- / 4*5 и х* от аргумента 23

та Х4 и параметра 23

УздзЫ

О 0,1

0,3 0,5 0,7 0,9 *3

О 0,2 0,4 0,6 0,8 1,0 1,2 г2

Рис. 4.4. Графики зави-

Рис. 4.5. Графики зависимостей / *4 5

симости

/з д 5 от аргу-

и х* от аргумента гг

мента хз и параметра гг

 

Аналогично решают задачу для третьего и второго шагов. Ре­ зультаты приведены на рис. 4.4-4.7.

/г,3,4,5

0,5

1,0

1,5

х2

0

0,2

0,4

0,6

0,8

1,0

1,2

1,4

г,

Рис. 4.6.

Графики

зави-

Рис. 4.7. Графики зависимостей

f

*3,4,5

симости

/ 2,з,4,5

от

аргу-

 

и х*

от аргумента г\

 

 

 

мента хг

и параметра г\

 

 

 

 

 

 

 

 

 

 

 

Спланируем

первый

шаг. Доход

на

 

пяти

шагах

(при

любом

управлении на первом шаге и оптимальном на последующих) опре­ деляем по формулам

/ 1 , 2 , 3 , 4 , 5 ( x i ) = / l ( 2 , Х\) + Я з , 4 , 5 ( 2 1 ) =

= 2 - (е“ *' + e - 2(2- * ■>) + ЯздзСх,), zi = 0,75xi + 0,3(г - х,).

На первом шаге имеющиеся ресурсы известны: z = 2. Поэтому функция / 1,2,3,4,5(21) опишется одной линией (рис. 4.8). Из рис. 4.8

/l,2,3,4,5

 

 

 

видно, что

максимальный

доход

будет

 

 

 

 

/*2 з 4 5 = 4,35; оптимальное

управление

 

 

 

 

на первом шаге —х* = 1,60. По известно­

 

 

 

 

му оптимальному управлению на первом

 

 

 

 

шаге можно

определить

запас

средств

 

 

 

 

к концу этого шага, что, в свою оче­

 

 

 

 

редь, определит оптимальное управление

 

 

 

 

на следующем шаге.

 

 

 

 

 

 

 

 

Вычислим запас средств к концу пер­

 

 

 

 

вого шага:

 

 

 

 

 

 

Рис. 4.8.

График

зави­

z \ = 0,75х* + 0,3(z - х \)

= 1,32.

симости

/ 1,2,з,4,5

от ар­

 

 

 

 

 

 

 

гумента Х\

 

Оптимальное

управление

на

 

втором

шаге тогда (см. рис. 4.7) будет иметь значение х \ =

1,02.

Остаток

средств к концу второго шага будет

 

 

 

 

 

 

 

 

 

z\ = 0,75^2 + 0,3(z* — Xj) = 0,86;

 

 

 

 

и, далее, Х3 = 0,62, z* = 0,54, х£ = 0,30, z \ = 0,30, x |

= 0.

 

 

Получим оптимальное управление, показывающее, сколько

средств надо вкладывать в первую отрасль:

 

 

 

 

 

 

х * = ( х * ; х \ \ х * ; x j ; х | ) = (1,60; 1,02; 0,62; 0,30; 0).

 

 

Автоматически получаем количество средств,

вкладываемых

во вторую отрасль:

 

 

 

 

 

 

 

У* = z - x * =

0,40,

у \ = z* - X * =

0,30,

Уз = z 2

-

Х 3 =

0,24,

У4 = *3 - *4

= 0,24,

у5* = zX - 2 ? = 0,30.

 

 

 

 

 

Определим остаток средств: 0,30 • 0,3 = 0,09.

§ 4.4. Задача о распределении средств поражения

Рассмотрим задачу о распределении средств поражения по обо­ роняющимся целям, которая аналогична задаче о резервировании средств.

Задача. Планируются боевые действия самолетами по оборо­ няющимся целям. Цели эшелонированы по глубине территории

N2 = 12, N3 = 15,

на четырех параллельных рубежах. Перед тем как выйти на данный рубеж, самолеты проходят зону действия огневых средств этого ру­ бежа, где подвергаются обстрелу со стороны последних. Огневые средства каждого рубежа могут вести огонь не только по средствам поражения, направляющимся непосредственно на цели данного ру­ бежа, но и по тем самолетам, которые проходят через зону действия, направляясь на следующие рубежи.

Вероятность поражения одного самолета, пролетающего зону действия орудий г-го рубежа, определяется формулой

где Ni — среднее число орудий, сохранивших боеспособность на г-м рубеже; оц — эффективность стрельбы орудий по самолетам.

Среднее число орудий г-ro рубежа, пораженных j -й волной вы­ летов самолетов, направленных на цели этого рубежа, вычисляется по формуле

где Ni — число орудий на г-м рубеже; Vj — среднее число самолетов в j -й волне вылетов, сохранивших боевые свойства до г-ro рубе­ жа; Pi — средняя вероятность поражения одного орудия г-го рубежа атакующим его самолетом.

Имеется N орудий на всех рубежах:

4

где Ni — число орудий на г-м рубеже: N\ = 10, ЛГ4 = 10, отсюда N = 47.

Имеется ZQ = 80 самолетов, которые необходимо распределить на четыре волны вылетов так, чтобы сделать максимальным среднее число поражений на всех рубежах:

4

я ®)= 2 /<(*)» t=i

где fi(x) — среднее число поражений на г-м рубеже.

Волны вылетов распределены по времени так, чтобы к момен­ ту подлета следующей волны предшествующая ей волна самолетов успела выполнить боевое задание.

Пусть

Р\ = 0,4,

Рг = 0,5,

Рз = 0,4,

Р* = 1,0,

ai = 0,05,

а.2 = 0,04,

аз = 0,04,

04 = 0,05.

Р е ш е н и е . Будем пользоваться общей схемой решения задач методом динамического программирования. Проведем оптимиза­ цию на четвертом (последнем) шаге. К зоне действия орудий 4-го рубежа подойдет z3 самолетов. Очевидно, что все они должны быть направлены на поражение орудий четвертого рубежа. Условное оп­ тимальное управление на 4-м шаге определяется по формуле

x j(z 3) = z3.

На четвертом рубеже установлено N 4 орудий, этот рубеж не под­ вергался бомбежке, т. е.

N 4 = N4.

Подлетевшие z3 самолетов поразят на четвертом рубеже / 4*(z3) орудий:

Q4(x*4) = m z 3) = N 4 ( 1 - е " й Р4) ,

где V4(z3) = z3(l —V4 ); V4 — вероятность поражения одного самолета на четвертом рубеже; V4 = 1 —е-а4ЛГ4, т. е. V4(z3) = г3е-а4Л/4.

Задавая различные значения z3, получим график функции / 4*(z3) (рис. 4.9). Вообще-то следовало бы брать значения z3 в интервале (0, 80), но в данной задаче в этом нет не­

 

обходимости.

 

 

 

Проведем

оптимизацию на

третьем

 

шаге. Зададим число Z2 самолетов, пре­

 

одолевших первый и второй

рубежи:

 

гг е (10,40). Для каждого значения Z2 вы­

 

числим суммарный выигрыш: на третьем

 

шаге при любом управлении, на четвер­

Рис. 4.9. График зависи­

том — при оптимальном, т. е.

 

 

 

 

мости /* от аргумента z3

/з,4 =

Яз(хз) + f t ( z з)-

 

Из z2 самолетов надо выделить жз на подавление орудий тре­ тьего рубежа, а 22 —жз самолетов направить на четвертый рубеж через зону огня третьего рубежа. Условное оптимальное управле­ ние хЦ г2) найдем из условия максимального выигрыша на двух последних шагах.

Здесь <2з(®з) — среднее число орудий, пораженных на третьем рубеже выделенными для этой цели жз самолетами. При таком управлении до четвертого рубежа дойдут 23 самолетов. Имеем

где v3 = жз(1 — vj) = жз(1 — 1 + е~агМг) = жзе- ®3^ 3.

Вычислим среднее число самолетов из оставшихся 22 —жз са­ молетов, прошедших через огонь третьего рубежа и дошедших до четвертого рубежа. На третьем рубеже осталось в действии Щ = Щ д 3(жз) орудий. Тогда среднее число дошедших само­ летов до четвертого рубежа из полного числа 22 —жз будет

23 =(22 - Ж3)е “3^ 3, где величина N 3 зависит от Ж3. Следовательно,

Задав значения 22 и взяв значения /4 (23) согласно графику, при­ веденному на рис. 4.9, для каждого из них получим зависимости /з,4(жз) (рис. 4.10). Отмечая точки максимума функции /з,4(жз), для каждого 22 строим по значениям, приведенным на рис. 4.10, графи­ ки зависимостей / 3*4 и ж* от 22 (рис. 4.11).

Оптимизируем решение на втором шаге. Процедура аналогична оптимизации, проведенной на третьем шаге.

Для разных значений z\ (число самолетов, преодолевших пер­ вый рубеж) вычислим (рис. 4.12) целевую функцию

h ,за = Q2OE2) + /*,4(^2)-

Здесь

Q2(X 2)

= N2(1 - е“ Я Рз),

-

04N2

v2 = х2е - “^ ,

22

= (21 - х 2)е~а2Й2,

N 2 = N 2 -

Q2(X 2).

Соседние файлы в папке книги