Учебное пособие по решению задач по курсу Экономико-математические методы и модели - Алесинская Т.В
.pdf• от существующего запаса в i-й строке отнимается величина сделанной перевозки, прежний запас зачеркивается, а вместо него записывается остаток,
т.е. (aiтек − bтекj .
Если существующий запас не позволяет перевезти всю потребность, то
•в клетку (i,j) в качестве перевозки вписывается значение запаса aiтек ;
•i-я строка вычеркивается, поскольку ее запас уже исчерпан;
•от существующей потребности в j-й строке отнимается величина сделанной перевозки, прежняя потребность зачеркивается, а вместо нее
записывается остаток, т.е. (bтекj − aiтек .
Нахождение опорного плана продолжается до тех пор, пока не будут вычеркнуты все строки и столбцы.
Метод минимального элемента
На каждом шаге метода минимального элемента из всех не вычеркнутых клеток транспортной матрицы выбирается клетка с минимальной стоимостью перевозки min cij . Заполнение выбранной клетки производится по правилам, описанным выше.
Метод Фогеля
На каждом шаге метода Фогеля для каждой i-й строки вычисляются штрафы di как разность между двумя наименьшими тарифами строки. Таким же образом вычисляются штрафы d j для каждого j-го столбца. После чего выбирается максимальный штраф из всех штрафов строк и столбцов. В строке или столбце, соответствующем выбранному штрафу, для заполнения выбирается не вычеркнутая клетка с минимальным тарифом min cij .
Если существует несколько одинаковых по величине максимальных штрафов в матрице, то в соответствующих строках или столбцах выбирается одна не вычеркнутая клетка с минимальным тарифом min cij .
71
Если клеток с минимальным тарифом также несколько, то из них выбирается клетка (i,j) с максимальным суммарным штрафом, т.е. суммой штрафов по i-й строке и j-му столбцу.
5.2. Методические рекомендации
Формально и реальные и фиктивные столбцы и строки в транспортной матрице абсолютно равноправны. Поэтому при нахождении опорных планов фиктивные строки, столбцы и тарифы необходимо анализировать и использовать точно так же как и реальные. Но при вычислении значения ЦФ фиктивные перевозки не учитываются, поскольку они реально не были выполнены и оплачены.
Если величина фиктивных тарифов превышает максимальный из реальных тарифов задачи [cф > max cij (i = 1,n; j = 1,m)], то методы
минимального элемента и Фогеля позволяют получить более дешевые планы перевозок, чем в случае с нулевыми фиктивными тарифами.
Задача № 5.01
Найти тремя методами опорный план ТЗ, в которой запасы на трех складах равны 210, 170, 65 ед. продукции, потребности четырех магазинов равны 125, 90, 130, 100 ед. продукции, тарифы перевозки в рублях за единицу продукции следующие:
|
5 |
8 |
1 |
2 |
|
|
2 |
5 |
4 |
9 |
|
|
. |
||||
|
9 |
2 |
3 |
1 |
|
|
|
Решение
Проверка сбалансированности задачи показывает, что суммарный объем запасов равен суммарному объему потребностей, т.е. введение фиктивных столбцов или строк не потребуется
72
запасы потребности
(""'""& ("""'"""&
210 +170 + 65 = 125 + 90 +130 +100 . |
|
%""$""# |
%"""$"""# |
445 ед.товара |
445 ед.товара |
Результаты нахождения опорного плана различными методами представлены в табл. 5.1, 5.2 и 5.3.
|
|
|
|
|
|
Таблица 5.1 |
Транспортная таблица с опорным планом северо-западного угла |
||||||
|
|
|
|
|
|
|
Пункты |
|
Пункты потребления, Bj |
Запасы, |
|||
отправления, Ai |
|
|
|
|
|
ед. продукции |
В1 |
|
В2 |
B3 |
B4 |
||
А1 |
125 |
|
85 |
|
|
210/85/0 |
|
5 |
8 |
1 |
2 |
||
A2 |
|
|
5 |
130 |
35 |
170/165/35/0 |
|
2 |
5 |
4 |
9 |
||
A3 |
|
|
|
|
65 |
65/0 |
|
9 |
2 |
3 |
1 |
||
Потребность, |
125/0 |
|
90/5/0 |
130/0 |
100/65/0 |
|
ед. продукции |
|
|
||||
|
|
|
|
|
|
|
Опорный план XСЗУ , найденный методом северо-западного угла |
125 |
85 |
0 |
0 |
|
|
|
|
0 |
5 |
130 |
35 |
|
[ед. товара]. |
XСЗУ = |
|
|||||
|
0 |
0 |
0 |
65 |
|
|
|
|
|
Соответствующая ЦФ (общие затраты на перевозку)
L (X СЗУ )=125 5 + 85 8 + 5 5 + 130 4 + 35 9 + 65 1= 2230 [руб.].
73
Таблица 5.2
Транспортная таблица с опорным планом минимального элемента
Пункты |
|
|
|
|
Пункты потребления, Bj |
|
|
Запасы, |
|
||||||||||||||
отправления, Ai |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ед. продукции |
|||||||
|
|
В1 |
|
|
|
В2 |
|
|
|
B3 |
|
B4 |
|
||||||||||
А1 |
|
|
|
|
|
5 |
|
45 |
|
8 |
|
130 |
|
1 |
35 |
2 |
210/80/45/0 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
A2 |
|
|
125 |
|
2 |
|
45 |
|
5 |
|
|
|
4 |
|
9 |
170/45/0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
A3 |
|
|
|
|
|
9 |
|
|
|
2 |
|
|
|
3 |
65 |
1 |
|
|
65/0 |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Потребность, |
|
125/0 |
|
|
90/45/0 |
|
130/0 |
|
100/35/0 |
|
|
|
|
|
|
||||||||
ед. продукции |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Опорный план XМЭ , найденный методом минимального элемента |
|
||||||||||||||||||||||
|
|
0 |
45 |
130 |
35 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
45 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
XМЭ = 125 |
[ед. товара], L(XМЭ )= 1100 [руб.]. |
|
|||||||||||||||||||||
|
|
0 |
0 |
0 |
65 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 5.3 |
|||
|
Транспортная таблица с опорным планом Фогеля |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
В1 |
|
В2 |
|
|
B3 |
|
|
B4 |
|
|
bi |
Штрафы строк,di |
||||||||||
А1 |
|
|
|
|
|
|
110 |
|
|
100 |
210/110/0 |
1 |
|
1 |
|
1 |
|
7 |
|||||
5 |
|
|
|
8 |
|
|
1 |
|
|
2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
A2 |
125 |
|
|
25 |
|
|
20 |
|
|
|
|
|
170/45/25/0 |
2 |
|
1 |
|
1 |
|
1 |
|||
2 |
|
|
|
5 |
|
|
4 |
|
|
9 |
|
|
|
||||||||||
A3 |
9 |
|
65 |
|
2 |
|
|
3 |
|
|
1 |
|
65/0 |
1 |
|
1 |
|
– |
|
– |
|||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
a j |
125/0 |
|
90/25/ |
130/20 |
|
100/0 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
0 |
|
/0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
3 |
|
|
|
3 |
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Штрафы |
– |
|
|
3 |
|
2 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||
столбцов, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d j |
– |
|
|
3 |
|
3 |
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|
|||
|
– |
|
|
3 |
|
3 |
|
|
|
|
– |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
74
На первом шаге нахождения опорного плана методом Фогеля возникает ситуация равенства значений максимальных штрафов транспортной матрицы
(см. табл. 5.3)
d 1столбца= d 2столбца= 3.
Минимальные тарифы в этих столбцах также совпадают
с21 = с32 = 2 .
Поэтому необходимо сравнить суммарные штрафы dij клеток (2,1) и (3,2)
d 21 = d 2строки+ d 1столбца= 2 + 3 = 5; d 32 = d 3строки+ d 2столбца =1+ 3 = 4.
Т.к. d21 > d32 , то выбираем на первом шаге для заполнения клетку (2,1).
Опорный план XФ , найденный методом Фогеля
|
0 |
0 |
110 |
100 |
|
|
|
|
25 |
20 |
0 |
|
[ед. товара], L(XФ )= 895 [руб.]. |
XФ = 125 |
|
|||||
|
0 |
65 |
0 |
0 |
|
|
|
|
|
5.3. Варианты задач для самостоятельного решения
Задача № 5.1
Найти тремя методами опорный план транспортной задачи, в которой запасы на трех складах равны 160, 140, 170 ед. продукции, потребности четырех магазинов равны 120, 50, 200, 110 ед. продукции, тарифы перевозки в рублях за единицу продукции следующие
|
7 |
8 |
1 |
2 |
|
|
4 |
5 |
9 |
8 |
|
|
. |
||||
|
9 |
2 |
3 |
6 |
|
|
|
Решите задачу для следующих случаев:
• фиктивные тарифы нулевые;
75
• фиктивные тарифы одинаковы по величине и превышают максимальный из реальных тарифов.
Сравните полученные опорные планы, соответствующие ЦФ и объясните причину их различия.
Задача № 5.2
Найти тремя методами опорный план транспортной задачи № 4.1 для случая, когда фиктивные тарифы больше максимального реального тарифа.
Задача № 5.3
Найти тремя методами опорный план транспортной задачи № 4.2 для случая, когда фиктивные тарифы больше максимального реального тарифа.
Задача № 5.4
Найти тремя методами опорный план транспортной задачи № 4.3 для случая, когда фиктивные тарифы больше максимального реального тарифа.
Задача № 5.5
Найти тремя методами опорный план транспортной задачи № 4.4 для случая, когда фиктивные тарифы больше максимального реального тарифа.
6.ОБЩАЯ РАСПРЕДЕЛИТЕЛЬНАЯ ЗАДАЧА ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
6.1. Теоретическое введение
Общая распределительная задача ЛП – это РЗ, в которой работы и ресурсы (исполнители) выражаются в различных единицах измерения. Типичным примером такой задачи является организация выпуска разнородной продукции на оборудовании различных типов.
76
Исходные параметры модели РЗ
1) n – количество исполнителей;
2)m – количество видов выполняемых работ;
3)ai – запас рабочего ресурса исполнителя Ai (i =1,n ) [ед. ресурса];
4)b j – план по выполнению работы Bj ( j =1,m ) [ед. работ];
5) |
cij – |
стоимость выполнения работы |
Bj |
исполнителем |
Ai |
[руб./ед. работ]; |
|
|
|
||
6) |
λij – |
интенсивность выполнения работы |
Bj |
исполнителем |
Ai |
[ед. работ / ед. ресурса]. |
|
|
|
||
|
|
Искомые параметры модели РЗ |
|
|
|
1) |
xij – планируемая загрузка исполнителя Ai при выполнении работ Bj |
||||
[ед. ресурса]; |
|
|
|
|
|
2) |
xijк – |
количество работ Bj, которые должен |
будет произвести |
||
исполнитель Ai [ед. работ]; |
|
|
|
||
3) |
L(X) – общие расходы на выполнение всего запланированного объема |
работ [руб.].
Этапы построения модели
I.Определение переменных.
II. Построение распределительной матрицы (см. табл. 6.1). III. Задание ЦФ.
IV. Задание ограничений.
77
Таблица 6.1
Общий вид распределительной матрицы
Исполнители, Ai |
|
|
|
|
|
Работы, |
Bj |
|
Запас ресурса, |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ед. ресурса |
|
|
В1 |
|
|
|
|
В2 |
|
… |
Bm |
|||||||
|
|
|
|
|
|
|
||||||||||
А1 |
λ11 |
|
|
|
λ12 |
|
… |
λ1m |
a1 |
|||||||
|
|
c11 |
|
|
|
|
c12 |
|
c1m |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
А2 |
λ21 |
|
|
|
λ22 |
|
… |
λ2m |
a2 |
|||||||
|
|
c21 |
|
|
|
|
c22 |
|
c2m |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
… |
|
… |
|
|
|
|
… |
|
… |
… |
… |
|||||
An |
λn1 |
|
|
|
λn2 |
|
… |
λnm |
an |
|||||||
|
|
cn1 |
|
|
|
|
cn2 |
|
cnm |
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
План, ед. работы |
|
b1 |
|
|
|
|
b2 |
|
… |
bm |
|
|||||
|
|
|
|
|
|
Модель РЗ |
|
|
|
|
||||||
L(X) |
|
n |
|
m |
|
|
|
|
||||||||
= |
∑ |
∑ сij(λijxij)→ min ; |
|
|||||||||||||
|
|
i = 1j = 1 |
|
|
|
|
||||||||||
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑xij = ai , i = 1,n, |
|
|
|
|
||||||||||
|
j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
∑λijxij = b j, j = 1,m, |
|
|
||||||||||||
|
i=1 |
|
≥ 0 (i = 1,n; j = 1,m |
), |
|
|
||||||||||
|
x |
ij |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где λijxij – это количество работ j-го вида, выполненных i-м исполнителем.
Этапы решения РЗ
I. Преобразование РЗ в ТЗ:
1) выбор базового ресурса и расчет нормированных производительностей ресурсов αi :
αi = |
λij |
; |
(6.2) |
|
λбаз j |
||||
|
|
|
78
2) |
пересчет запаса рабочего ресурса исполнителей a′i : |
|
|
|
||||||||
|
|
|
a′i = αiai [ед. ресурса]; |
|
|
|
|
(6.3) |
||||
3) |
пересчет планового задания b′j: |
|
|
|
|
|
||||||
|
b′ |
= |
b j |
|
|
ед.работ ед. ресурса |
= ед. ресурса ; |
(6.4) |
||||
|
|
|
|
|
||||||||
|
j |
|
λбаз j |
ед. работ |
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||
4) |
пересчет себестоимостей работ: |
|
|
|
|
|
||||||
|
c′ij = cijλбаз j |
|
руб. ед.работ |
= |
руб. |
|
|
|
||||
|
|
|
|
|
|
. |
(6.5) |
|||||
|
ед. работ ед. ресурса |
|
|
|||||||||
|
|
|
|
|
|
|
ед. ресурса |
|
n |
m |
II. Проверка баланса пересчитанных параметров ∑a′i = ∑b′j и |
|
i=1 |
j=1 |
построение транспортной матрицы.
III. Поиск оптимального решения ТЗX'* = (x'*ij ).
IV. Преобразование оптимального решения ТЗ X'* в оптимальное решение РЗ X*, причем переход X'* → X* выполняется по формуле (6.6)
xij = |
xij' |
(6.6) |
|
|
[ед. ресурса], |
||
|
|||
|
αi |
|
где xij и xij' – соответственно элементы решения РЗ и ТЗ.
V. Определение количества работ Xк* = (xijк*), соответствующее
оптимальному решению РЗ X*:
xijк |
|
ед. работ ед. ресурса |
|
|
|
= λijxij |
|
= ед. работ . |
(6.7) |
||
ед. ресурса |
|||||
|
|
|
|
VI. Определение ЦФ распределительной задачи L( X*) согласно (6.1).
79
6.2. Методические рекомендации
Задача № 6.01
На фабрике эксплуатируются три типа ткацких станков, которые могут выпускать четыре вида тканей. Известны следующие данные о производственном процессе:
• производительности станков по каждому виду ткани, м/ч
|
24 |
30 |
18 |
|
42 |
|
|
|
|
12 15 |
|
|
|
|
; |
||
(λij )= |
9 21 |
|||||||
|
8 |
|
10 |
6 |
14 |
|
|
|
|
|
|
|
|||||
• себестоимость тканей, руб./м |
|
|
|
|
|
|
||
|
|
2 |
1 |
3 |
1 |
|
|
|
|
|
3 2 |
4 |
1 |
|
|
|
|
(cij )= |
; |
|
|
|||||
|
|
6 |
3 |
5 |
2 |
|
|
|
|
|
|
|
|
•фонды рабочего времени станков (ai ): 90, 220, 180 ч;
•планируемый объем выпуска тканей ( b j): 1200, 900, 1800, 840 м.
Требуется распределить выпуск ткани по станкам с целью минимизации общей себестоимости производства ткани.
Решение
Пусть переменные xij – это время, в течение которого i-й станок будет
выпускать j-ю ткань. Сведем исходные данные задачи в распределительную таблицу (табл. 6.2).
80