где аij и bi - случайные величины, di - заданные уровни вероятности. Обозначим:
n |
|
Yi = bi − ∑aij ×d j, |
(4.58) |
j=1 |
|
где Yi - случайная величина. |
|
Обычно принимают, что случайные величины cj,aij,bi,yi |
подчиняются |
закону нормального распределения с известным математическим ожиданием и дисперсией.
Подставив (4) в неравенство (3), получим:
P[уi ≥ 0] ≥ di , |
(4.59) |
P[уi ≥ 0] ≤ di , |
(4.60) |
P[уi ≤ 0] ≥ di , |
(4.61) |
P[уi ≤ 0] ≤ di . |
(4.62) |
Случайная величина yi при независимых di и Qi будет иметь математическое ожидание и дисперсию:
n
Yi = bi − ∑aij ×xi , (4.63)
j=1
δуi |
2 |
= Qi |
2 |
n |
2 ×x j |
2 , |
|
+ ∑δij |
(4.64) |
||||||
|
|
|
|
j=1 |
|
|
|
где Qi 2 - дисперсия случайной величины bi. Для первого варианта ограничения (4.54) можно записать:
|
≥ tdi ×δуi, |
(4.65) |
Yi |
где tdi - коэффициент, учитывающий закон распределения случайной величины, определяемый аналитически или таблично в зависимости от значения вероятности di;
δyi - среднеквадратическое отклонение случайной величины yi. Подставив в (4.63) значения Yi и δi , получим
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
n |
2 ×x j |
2 +Qi |
2 . |
|
||
|
bi |
− ∑ |
aij |
×x j ≥ tdi × |
∑δij |
(4.66) |
||||
|
|
j=1 |
j=1 |
|
|
|
|
|
||
После преобразований:
n |
n |
|
|
|
|
|
||
∑ |
|
×x j + tdi × |
∑δij ×x j2 +Qi |
2 ≤ |
|
. |
(4.67) |
|
aij |
bi |
|||||||
j=1 |
j=1 |
|
|
|
|
|
||
Если сравним выражение (4.67) с аналогичным ограничением в детерминированной постановке
n
∑Qij × x j ≤ bi , (4.68)
j=1
66
то увидим, что ограничение в стохастической постановке отличается двумя признаками:
1) выполнен переход от детерминированных значений к математическим ожиданиям случайных величин Qij и bi ;
2) появился дополнительный член:
ξ = t |
|
× |
n |
2 |
×x 2 |
+Q 2 |
, |
(4.69) |
di |
∑δ |
|||||||
1 |
|
j=1 |
ij |
j |
i |
|
|
|
|
|
|
|
|
|
|
|
который учитывает все вероятностные характеристики задачи:
-закон распределения с помощью tdi;
-заданный уровень вероятности di, дисперсию случайных величин Qij,
равную δij2 , и дисперсию случайных величин bi, равную Qi 2 . Таким образом, получим
n
∑Qij ×x j +ξi ≤ bi . (4.70)
j=1
Решение задач стохастического программирования в такой постановке возможно методами сепарабельного программирования, потому что ограничения задачи не являются линейными функциями.
Поскольку в ограничениях появился дополнительный положительный член ξ1 , это приведет к тому, что потребуется большая величина ресурса bi
по сравнению с детерминированной постановкой.
4.6. Модели оптимального планирования транспортного типа
Если требуется решение вопросов о выборе схемы прикрепления поставщиков и потребителей продукции, используются модели транспортного типа. Классическая транспортная задача заключается в планировании прикрепления поставщиков к потребителям продукции и формулируется следующим образом: однородный продукт, находящийся в m пунктах производства в количестве Р1,Р2,...Рm, требуется доставить в n пунктов потребления. Потребность продукции в этих пунктах равна S1,S2,...Sn.
Экономико-математическая модель задач транспортного типа: целевая функция - затраты на перевозку продукта должны быть мини-
мальными:
m n |
|
F = ∑∑cij ×xij → max. |
(4.71) |
j=1 j=1 |
|
Ограничения:
1. Вся продукция от предприятий-поставщиков отправляется потребителям:
67
n
∑xij = Pi ; i =1,m. (4.72)
j=1
2.Все потребители обеспечены продукцией:
m
∑xij =Sj; j =1,n. (4.73)
j=1
3. Мощность поставщиков равна потребности в продукции (условие закрытости):
m |
n |
|
∑Pi = ∑Sj. |
(4.74) |
|
j=1 |
j=1 |
|
Модификации транспортной задачи позволяют учитывать особенности различных хозяйственных условий, а именно:
1. Запрет каких-либо перевозок.
Если между поставщиками и потребителями продукции не существует маршрутов (связей) или ими нельзя пользоваться, можно задать стоимость перевозки сij, намного превышающую стоимость остальных перевозок (на-
пример, 99999).
2. Ограниченность пропускных способностей коммуникаций.
Это условие учитывается при введении ограничений, лимитирующих наибольшее значение объема перевозки по конкретному маршруту:
0 ≤ xij ≤ dij, |
(4.75) |
где dij - пропускная способность транспортной линии.
3. Нарушение условия равенства производства и потребления (откры-
тая транспортная задача).
Если не вся продукция нужна потребителям, т.е. ∑Pi > ∑Sj ,то огра- |
|
i |
j |
ничение на продукцию, отправляемую из пунктов производства, принимает вид
n
∑xij ≤ Pi ; i =1,m. (4.76)
j=1
Транспортная задача сводится к классическому виду путем введения фиктивного потребителя S n+1 с потребностью
Sm+1 |
= ∑Pi −∑Sj. |
(4.77) |
|
|
i |
j |
|
В целевой функции должны учитываться затраты, связанные с хранением и с потерей излишней продукции в каждом пункте производства.
Если суммарный объем производства меньше суммарного объема потребления, необходимо учитывать не только транспортные расходы, но и ущерб от недопоставок. В этой задаче
∑Pi ≤ ∑Sj |
(4.78) |
|
i |
j |
|
|
68 |
|
и ограничения на продукцию, поступающую в каждый пункт потребления, будут
∑αij ≤ Sj; j = |
|
|
(4.79) |
1, n. |
|||
i |
|
||
Этот случай также сводится к классической транспортной задаче путем введения фиктивного поставщика с объемом производства:
Pm = ∑Sj −∑Pi . |
(4.80) |
|
j |
i |
|
4.7. Решение задач по планированию перевозок
Задача по планированию перевозок формулируется следующим образом: необходимо составить план транспортирования строительных материалов, минимизирующий затраты на перевозки и издержки, связанные с тем, что часть продукции остается у поставщиков.
Исходные данные представлены в табл. 4.12.
Таблица 4.12
Исходные данные для расчета
Потребители |
|
Поставщики |
|
Потребность в |
|
|
|
|
|||
База №1 |
База №2 |
База №3 |
материалах |
||
|
|||||
|
|
||||
|
|
|
|
|
|
Объект № 1 |
10 |
12 |
8 |
100 |
|
|
|
|
|
|
|
Объект № 2 |
11 |
7 |
13 |
120 |
|
|
|
|
|
|
|
Объект № 3 |
15 |
13 |
9 |
200 |
|
|
|
|
|
|
|
Объект № 4 |
8 |
11 |
8 |
160 |
|
|
|
|
|
|
|
Мощности по- |
210 |
340 |
200 |
|
|
ставщиков, т |
|
||||
|
|
|
|
В соответствующих клетках таблицы задана стоимость перевозок 1 тонны груза от поставщиков к потребителям - сij, тыс.р. за 1 тонну. Потери, связанные с хранением продукции у поставщиков, составляют: 5, 7, и 4 тыс.р. за 1 тонну для базы N 1, 2 и 3 соответственно.
Решим задачу с использованием программы симплекс-метода. С учетом переменных, определяющих объем грузов, остающихся на базах, размерность задачи составляет 3×5 =15 переменных.
Симплекс-матрица представлена в табл. 4.13.
69
Таблица 4.13 Симплекс-матрица для решения транспортной задачи
Номер |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вид |
Правая |
Х1 |
Х2 |
Х3 |
Х4 |
Х5 |
Х6 |
Х7 |
X8 |
Х9 |
Х10 |
Х11 |
Х12 |
Х13 |
Х14 |
Х15 |
часть огра- |
||
строки |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
связи |
ничения |
F |
10 |
12 |
8 |
11 |
7 |
13 |
15 |
13 |
9 |
8 |
11 |
8 |
5 |
7 |
4 |
→ |
min |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
1 |
|
|
1 |
|
|
1 |
|
|
1 |
|
|
= |
210 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
|
1 |
|
|
1 |
|
|
1 |
|
|
1 |
|
= |
340 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
1 |
|
|
1 |
|
|
1 |
|
|
1 |
|
|
1 |
= |
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
1 |
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
= |
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
1 |
1 |
1 |
|
|
|
|
|
|
|
|
|
= |
120 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
1 |
1 |
1 |
|
|
|
|
|
|
= |
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
|
|
1 |
1 |
1 |
|
|
|
= |
160 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение задачи в системе электронных таблиц EXCEL осуществляется с помощью пункта меню «Сервис», «Анализ данных». Необходимо представить симплекс-матрицу в виде формул.
Результаты решения представлены в табл. 4.14.
|
Таблица 4.14 |
|
Результаты решения |
||
Индексы базисных переменных |
Оптимальные значения базисных пере- |
|
|
меных |
|
3 |
100 |
|
|
|
|
5 |
120 |
|
|
|
|
8 |
100 |
|
|
|
|
9 |
100 |
|
|
|
|
10 |
160 |
|
|
|
|
13 |
50 |
|
14 |
120 |
|
|
|
|
Оптимальное значение функции цели 6210,0.
По результатам расчета можно сделать следующие выводы. Поскольку x15=0, вся продукция базы N 3 направляется потребителям. На базе N 1 остается 50 тонн продукции, на базе N 2 - 120 тонн. С учетом того, что потери, связанные с хранением нереализованной продукции учтены в целевой функции, оптимальная стоимость перевозок составит
6210 - 50 × 5 - 120×7 = 5120 тыс.р.
70