- •Операций в экономике Учебное пособие
- •1 Сибирский федеральный университет, 2007
- •Раздел 1. Исследование операций: линейные модели в экономике 5
- •Тема 1.1. Линейные модели в операционном анализе экономических систем ..5 Лекция 1.1.1. Основы методологии моделирования 5
- •Тема 1.2. Теория двойственности в операционном анализе экономических систем 37
- •Раздел 2. Нелинейные и специальные модели исследования операций 71
- •Тема 2.1. Нелинейность в экономических процессах 71
- •Тема 2.2. Специальные модели исследования операций в экономике 98
- •Раздел 1. Исследование операций: линейные модели в экономике
- •Тема 1.1. Линейные модели в операционном анализе экономических систем
- •Сфера применимости и предпосылки построения линейных моделей в экономике
- •Задача оптимального планирования производства
- •Задача о диете
- •Задача о раскрое
- •Транспортная задача
- •Задача о назначениях
- •Постановка задачи и основные определения злп
- •Формы записи общей злп
- •Метод прямого перебора
- •Метод искусственного базиса
- •Тема 1.2. Теория двойственности в операционном анализе экономических систем
- •Правила построения
- •Лекция 1.2.2. Постоптимальный анализ решения линейных моделей с использованием двойственных оценок: геометрическая интерпретация
- •Лекция 1.2.3. Постоптимальный анализ решения линейных моделей с использованием двойственных оценок: геометрическая интерпретация (продолжение)
- •Ресурса наиболее выгодно
- •Аналитический подход
- •Этап 1. Определение статуса ресурсов
- •Лекция 1.2.5. Постоптимальный анализ решения линейных моделей с использованием двойственных оценок: аналитический подход (продолжение)
- •Раздел 2. Нелинейные и специальные модели исследования операций
- •Тема 2.1. Нелинейность в экономических процессах
- •Нелинейного программирования
- •Метод множителей Лагранжа
- •Выпуклого программирования
- •Ограничения типа «неравенств»
- •Ограничения «смешанного типа»
- •Динамических задач.
- •Задача распределения капиталовложений
- •Модель динамического программирования
- •Задача о загрузке
- •Тема 2.2. Специальные модели исследования операций в экономике
- •Определение опорного плана транспортной задачи
- •Метод северо-западного угла
- •Транспортной задачи
- •Метод потенциалов
- •Лекция 2.2.3. Целочисленные задачи исследования операций. Метод Гомори для нахождения их решения. Задача о назначениях и венгерский метод
- •Задача о назначениях
- •Метод Гомори
- •Входящий поток заявок на обслуживание
- •Механизм обслуживания
- •Дисциплина очереди
- •Классификация систем массового обслуживания
- •«Первая пришла — первая обслужена»,
- •«Последняя пришла — первая обслужена»
- •Потоки событий
- •Понятие марковского случайного процесса
- •Задачи анализа поведения системы
- •Статистические задачи
- •Операционные задачи
- •Лекция 2.2.5. Многоканальные смо: многоканальные и одноканальные системы массового обслуживания
- •Одноканальные системы массового обслуживания с ограниченной длиной очереди
- •Лекция 2.2.6. Модели принятия решений в условиях неопределенности и риска
- •Свойства решений матричных игр
- •Игры порядка 2×2
- •Графический метод решения игр 2×n и m×2
- •Критерий Лапласа
- •Критерий Вальда
- •Критерий Гурвица
- •Критерий Сэвиджа
- •1. Ежегодные затраты, не зависящие от числа построенных комнат
- •2. Ежегодные затраты, пропорциональные числу построенных комнат, в долл. (табл. 2.26)
- •3. Ежегодные затраты, пропорциональные среднему числу занятых комнат r (табл. 2.27)
- •660041 Красноярск, пр. Свободный, 79.
Транспортной задачи
Для определения оптимального плана транспортной задачи разработано несколько методов. Однако наиболее часто используется метод потенциалов.
Метод потенциалов
Общий принцип определения оптимального плана транспортной задачи методом потенциалов аналогичен принципу решения задачи линейного программирования симплексным методом, а именно сначала находят опорный план транспортной задачи, а затем его последовательно улучшают до получения оптимального плана.
Для определения опорного плана транспортной задачи будем пользоваться одним из методов, изложенных ранее. Эти методы гарантируют получение занятых в исходном плане п+т-1 клеток, причем в некоторых из них могут стоять нули. Построенный план следует проверить на оптимальность.
* , *.
(i = 1,m j = 1,n) транспортной задачи существуют такие n+m чисел α1,α2,...,αm,β1,β2,...,βn,что
ai+Pj=cijnpux*j>0 (2.21)
и
αi+βj≤cij при jc*. =0 (2.22)
д ля всех i = 1,m и j = 1,n, то X = (xij*) — оптимальный план транспортной задачи.
Ч
исла
αiи
βj,
(i
= 1,m j=1,
n)
называются потенциалами,
соответственно
пунктов назначения и пунктов потребления.
109
Сформулированная теорема позволяет построить алгоритм нахождения решения транспортной задачи. Он состоит в следующем. Пусть одним из рассмотренных выше методов найден опорный план транспортной задачи. Для каждого из пунктов отправления и назначения определяют потенциалы αi и βj (i = 1,m j = 1,n). Эти числа находят из системы уравнений
αi+βj=cij, (2.23)
где Су — тарифы, стоящие в заполненных клетках таблицы условий транспортной задачи.
Так как число заполненных клеток равно п + т - 1, то система (2.23) с п+m неизвестными содержит п + т - 1 уравнений. Поскольку число неизвестных превышает на единицу число уравнений, одно из неизвестных можно положить равным произвольному числу, например 0α 1 =, и найти последовательно из уравнений (2.23) значения остальных неизвестных. Для исследования плана на оптимальность для каждой свободной клетки проверяется условие αi + βj ≤ cij. Если хотя бы одна свободная клетка не удовлетворяет
данному условию, то опорный план не будет оптимальным, его можно улучшить за счет загрузки этой клетки. Если таких клеток несколько, то наиболее перспективна для загрузки клетка, для которой разность (платеж) между тарифом клетки и суммой потенциалов наименьшая, т. е.
Экономически платеж показывает, на сколько денежных единиц уменьшатся транспортные издержки от загрузки данной клетки единицей груза.
Если для всех свободных клеток чисела 0γ > 0, то опорный план перевозок является оптимальным.
Если для опорного плана перевозок указанное условие оптимальности не выполняется, то за счет загрузки свободной клетки с отрицательным платежом план перевозок улучшается. Для наиболее перспективной свободной клетки строится замкнутый цикл с вершинами в загруженных клетках. Для
ПО
этого рассматривают все свободные клетки, для которых уг < 0 и среди данных чисел выбирают наименьшее. Клетку, которой это число соответствует, следует заполнить, она далее становится вершиной цикла. Вершинам этого цикла условно приписываются знаки: свободной клетке — плюс, следующей по часовой или против часовой стрелки занятой клетке — минус, следующей — снова плюс и т. д. Из поставок в клетках цикла с «отрицательными» вершинами выбирается наименьшее количество груза, которое и перемещается по клеткам этого цикла: прибавляется к поставкам в положительных вершинах и вычитается из поставок в отрицательных, в результате чего баланс цикла не нарушится.
Циклом в таблице условий транспортной задачи называется ломаная линия, вершины которой расположены в занятых клетках таблицы, а звенья — вдоль строк и столбцов, причем в каждой вершине цикла встречается ровно два звена, одно из которых находится в строке, а другое — в столбце.
Если ломаная линия, образующая цикл, пересекается, то точки самопересечения не являются вершинами. Примеры некоторых циклов показаны на рис. 2.5.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
—1__ |
id— |
Рис. 2.5. Примеры построения циклов в транспортной задаче
При правильном построении опорного плана для любой свободной клетки можно построить лишь один цикл.
В результате указанных выше перемещений грузов в пределах клеток, связанных циклом с данной свободной клеткой, определяют новый опорный план транспортной задачи.
Описанный выше переход от одного опорного плана транспортной задачи к другому называется сдвигом по циклу пересчета.
111
Следует отметить, что при сдвиге по циклу пересчета число занятых клеток остается неизменным, а именно остается равным п + т - 1. При этом если в минусовых клетках имеется два (или более) одинаковых числа xij, то освобождают лишь одну из таких клеток, а остальные оставляют занятыми (с нулевыми поставками).
Новый опорный план транспортной задачи проверяют на оптимальность. Для этого определяют потенциалы пунктов отправления и назначения и находят числа уу для всех свободных клеток. Если среди этих чисел не
окажется отрицательных, это свидетельствует о получении оптимального плана. Если же отрицательные числа имеются, то следует перейти к новому опорному плану. В результате итерационного процесса после конечного числа шагов получают оптимальный план задачи.
Из изложенного следует, что процесс нахождения решения транспортной задачи методом потенциалов включает такие этапы:
Находят опорный план. При этом число заполненных клеток должно быть равным п + т - 1.
Находят потенциалы βj и αi, соответственно, пунктов назначения и
отправления.
Для каждой свободной клетки определяют число у у. Если среди чи сел γ у нет отрицательных, то получен оптимальный план транс портной задачи; если же они имеются, то переходят к новому опор ному плану.
Среди отрицательных чисел у у выбирают наименьшее, строят для
свободной клетки, которой оно соответствует, цикл пересчета и производят сдвиг по циклу пересчета.
5. Полученный опорный план проверяют на оптимальность, т. е. снова повторяют все действия, начиная с этапа 2.
В заключение отметим, что при определении опорного плана или в процессе решения задачи может быть получен вырожденный опорный план. Чтобы избежать в этом случае зацикливания, следует соответствующие нулевые элементы опорного плана заменить сколь угодно малым положитель-
112
ным числом ε и решать задачу как невырожденную. В оптимальном плане такой задачи необходимо считать ε равным нулю.
Пример. В трех хранилищах горючего ежедневно хранится 50, 30 и 10т бензина соответственно. Этот бензин ежедневно получают четыре заправочные станции в количествах, соответственно равных 30, 30, 10 и 20 т. Стоимость перевозки одной тонны бензина из хранилищ к заправочным станциям задаётся матрицей С:
'[ 2 4 1
С =
3 15
2 4 4
Составить такой план перевозок бензина, при котором общая стоимость перевозок будет наименьшей.
Решение. Построим таблицу транспортной задачи.
Таблица 2.20
Пункты отправления |
Пункты назначения |
Запасы |
|||||||
|
B2 |
B3 |
B4 |
|
|||||
Ах |
30 |
1 |
90 |
2 |
|
4 |
|
1 |
50 |
|
|
|
|
|
|
|
|
||
A2 |
|
2 |
10 |
3 |
10 |
1 |
10 |
5 |
30 |
|
|
|
|
|
|
|
|
||
Аз |
|
3 |
|
2 |
|
4 |
10 |
4 |
10 |
|
|
|
|
|
|
1W |
|
||
Потребности |
30 |
30 |
10 |
20 |
90 |
||||
Сначала, используя метод северо-западного угла, находим опорный план задачи. Этот план записан в табл. 2.20. Найденный опорный план проверяем на оптимальность. В связи с этим находим потенциалы пунктов отправления и назначения.
113
Для определения потенциалов получаем систему
«!+Д=1, а1+Р2=2, a2+j32=3, a2+j33=\, a2+j34=5,
α3+β4=4,
содержащую шесть уравнений с семью неизвестными. Полагая 0α 1 =, находим:
Д =1, 2β2 = , 1α2=, ръ =0, рА =4, 0α3=.
Для каждой свободной клетки вычисляем число γij =cij −(αi + βj):
7i3=4> 7i4=-3, 72i=°> 7з1=2' 7з2=°> 7зз=4- Так как сРеДи чисел Гц имеется отрицательное {уи =-3), то построенный план перевозок не является оптимальным и надо перейти к новому опорному плану.
Таблица 2.21
Пункты |
Пункты назначения |
|
|
|||||||
отправления |
Вх |
B2 |
Вз |
ВА |
|
|||||
Ах |
30 |
1 |
20 |
2 |
|
4 |
|
1 |
50 |
0 |
|
|
i— |
|
> |
||||||
A2 |
|
2 |
104 |
3 |
10 |
1 |
|
30 |
1 |
|
|
|
lut_ _ |
|
|
||||||
Аз |
|
3 |
|
2 |
|
4 |
10 |
4 |
10 |
0 |
Потребности |
30 |
30 |
10 |
20 |
90 |
|
||||
fij |
1 |
2 |
0 |
4 |
|
|
||||
Для свободной клетки с 3γ14=− строим цикл пересчета (табл. 2.21) и производим сдвиг по этому циклу. Наименьшее из чисел в минусовых клетках равно 10 = min(20, 10). Клетка, в которой находится это число, становится свободной в новой табл. 2.22. Другие числа в табл. 2.22 получаются так: к числу 10, стоящему в плюсовой клетке табл. 2.21, добавим 10 и вычтем 10 из числа 20, находящегося в минусовой клетке табл. 2.21. Клетка на пересе-
114
чении строки А2 и столбца В4 становится свободной. После этих преобразований получаем новый опорный план (табл. 2.22).
Этот план проверяем на оптимальность. Снова находим потенциалы пунктов отправления и назначения. Для этого составляем следующую систему уравнений:
α1+β1=1, α1+β2=2, α1+β4=1, α2+β2=3, α2+β3=1, α3+β4=4. Полагаем α1 = 0, получаем
Д =β4 = 1, β2 =2, β3 =0, α3 =3, а2 =. Для каждой свободной клетки вычисляем число у у; имеем:
У\Ъ =4> 721=0> /24=3, Гп=~1> /32="3, 733 =1-
Таблица 2.22
Пункты |
Пункты назначения |
|
ОС |
|||||||
отправления |
Вх |
B2 |
B3 |
ВА |
|
|||||
Ах |
30 |
1 |
10 |
2 |
|
4 |
10 |
1 |
50 |
0 |
A2 |
|
2 |
20 |
3 |
10 |
1 |
|
5 |
30 |
1 |
Аз |
|
3 |
|
2 |
|
4 |
10 |
4 |
10 |
3 |
Потребности |
30 |
30 |
10 |
20 |
90 |
|
||||
fij |
1 |
2 |
0 |
1 |
|
|
||||
Таким образом, видим, что данный план перевозок не является оптимальным. Поэтому переходим к новому опорному плану (табл. 2.23).
Сравнивая значения αi + βj новых потенциалов, отвечающих свободным клеткам табл. 2.23, с соответствующими числами сij, видим, что указанные значения сумм потенциалов для всех свободных клеток не превосходят соответствующих чисел сij.
115
Таблица 2.23
Пункты |
Пункты назначения |
Запасы |
а |
|||||||
отправления |
Вх |
B2 |
Вз |
ВА |
|
|||||
Ах |
30 |
1 |
|
2 |
|
4 |
20 |
1 |
50 |
0 |
A2 |
|
2 |
20 |
3 |
10 |
1 |
8 |
5 |
30 |
0 |
Аз |
|
3 |
10 |
2 |
|
4 |
|
4 |
10 |
1 |
Потребности |
30 |
30 |
10 |
20 |
90 |
|
||||
βj |
1 |
3 |
1 |
1 |
|
|
||||
Следовательно, полученный план
'30 0 0 20^
0 20 10 0 v0 10 0 О,
является оптимальным.
При данном плане стоимость перевозок
F=1·30 +1·20 + 3·20+1· 10 + 2·10= 140.
116
