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

Методы линейной и сетевой оптимизации

.pdf
Скачиваний:
0
Добавлен:
24.11.2025
Размер:
1.45 Mб
Скачать

 

Вариант 18

 

 

 

Продукция

А

В

С

Запасы

Сырье

 

 

 

сырья, ед.

I

-

-

2

12

II

4

1

2

30

III

4

1

4

42

Прибыль, ден. ед.

14

2

8

 

 

Вариант 19

 

 

 

Продукция

А

 

В

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

2

 

-

-

8

II

2

 

3

1

18

III

4

 

3

-

24

Прибыль, ден. ед.

6

 

9

1

 

 

 

Вариант 20

 

 

Продукция

А

В

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

 

3

3

1

22

II

 

5

3

-

28

III

 

6

6

1

42

Прибыль, ден. ед.

 

12

18

2

 

 

 

Вариант 21

 

 

Продукция

А

В

 

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

2

1

 

4

20

II

-

-

 

1

4

III

3

-

 

2

18

Прибыль, ден. ед.

3

1

 

6

 

 

 

Вариант 22

 

 

 

Продукция

А

В

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

 

5

1

6

38

II

 

2

1

5

24

III

 

3

-

3

22

Прибыль, ден. ед.

 

6

2

12

 

 

 

Вариант 23

 

 

 

Продукция

А

В

 

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

-

2

 

2

16

II

1

1

 

-

4

III

-

1

 

2

14

Прибыль, ден. ед.

1

3

 

2

 

 

 

Вариант 24

 

 

 

Продукция

А

В

С

Запасы

Сырье

 

 

 

 

сырья, ед.

I

 

1

3

2

20

II

 

1

2

2

18

III

 

-

3

4

30

Прибыль, ден. ед.

 

3

9

6

 

20

III. Транспортная задача

Общая постановка транспортной задачи состоит в определении оптималь-

ного плана перевозок некоторого однородного груза из m пунктов отправления

А1, А2, ... , Аm в n пунктов назначения B1, B2, ... , Bn . При этом в качестве критерия оптимальности берется либо минимальная стоимость перевозок всего груза, либо минимальное время его доставки.

Математическая модель транспортной задачи имеет вид:

m

n

 

 

 

Z cij xij

min ,

(3.1)

i 1

j 1

 

 

 

m

 

 

 

 

 

xij

bj

( j

1,..., n)

(3.2)

i 1

 

 

 

 

 

n

 

 

 

 

 

xij

ai

(i

1,..., m)

(3.3)

j 1

 

 

 

 

 

xij 0,

(i 1,..., m; j 1, ..., n)

(3.4)

Здесь cij – тарифы перевозки единицы груза из i - го пункта отправления

Аi в j - ый пункт назначения Bj , bj – потребность в грузе в j - ом пункте на-

значения, ai – запасы груза в i - ом пункте отправления.

Наличие линейных уравнений (3.2) и (3.3) обеспечивает доставку необхо-

димого количества груза в каждый из пунктов назначения и вывоз всего имею-

щегося груза из всех пунктов отправления. При этом, ввиду (3.4), исключаются обратные перевозки. Задача (3.1) – (3.4) является частным случаем задачи линейного программирования, однако, в силу своей специфики решается

специальным методом.

Если выполняется так называемое условие баланса

m

n

 

ai b j ,

(3.5)

i 1

j 1

 

21

то такая транспортная задача называется закрытой. Если условие баланса (3.5)

не выполняется, то задача называется открытой.

Теорема 1. Для разрешимости транспортной задачи необходимо и доста-

точно, чтобы выполнялось условие баланса (3.5).

Для решения задачи (3.1) – (3.4) применяется метод потенциалов, который по существу является другой формой симплекс-метода.

 

 

 

 

 

 

m

n

Замечание. Если условие баланса нарушено,

т.е. ai bj , то вводят

 

 

 

 

 

 

i 1

j 1

фиктивного

 

при

m

n

 

или

потребителя

поставщика

ai bj

 

 

 

i 1

j 1

 

 

 

 

при

m

n

 

с поставкой

m

n

или потребностями

 

ai bj

am 1 ai bj

 

 

i 1

j 1

 

 

i 1

j 1

 

 

n

m

bn 1

b j ai соответственно. Стоимости соответствующих перевозок

 

j 1

i 1

полагаются равными нулю. При построении начального опорного плана в этом случае фиктивные клетки заполняются в последнюю очередь. Опишем алгоритм метода на примере.

Пример. Компания контролирует 3 фабрики, производительность кото-

рых в неделю (в тыс. изделий) задается вектором a (19, 12, 11) . Компания заключила договоры с четырьмя заказчиками, еженедельные потребности кото-

рых (в тыс.изделий)

задаются вектором b (9, 10, 8, 11) Стоимость транс-

портировки 1

тыс. изделий с i-ой фабрики j-му заказчику задается матрицей

 

 

 

8

5

3

7

 

тарифов

C

 

2

4

6

3

 

 

.

 

 

 

5

2

3

7

 

 

 

 

 

Требуется определить оптимальный план перевозок с целью минимиза-

ции суммарных затрат на транспортировку.

22

Р е ш е н и е. 1. Построение математической модели. Обозначим через

xij

(i

1,3;

 

j

1,4)

количество продукта, перевозимого от i-го поставщика j-му

потребителю. Тогда общая стоимость перевозок

 

 

 

 

 

 

 

 

 

 

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f cij xij 8x11 5x12 3x13

7x14 2x21 4x22 6x23 3x24

 

 

 

i 1

j 1

 

 

 

 

 

 

 

 

 

 

 

(3.6)

 

 

 

 

 

 

5x31 2x32 3x33 7x34.

 

 

 

 

 

 

 

 

 

 

Поскольку суммарный объем вывезенных изделий не может превышать

 

 

 

 

 

 

 

 

 

 

 

их

объемов производства,

то

переменные

xij (i

1,3;

j

1,4)

должны

удовлетворять следующим ограничениям:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

x

x

19,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

12

13

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x21 x22

x23

x24

12,

 

 

 

 

 

 

 

(3.7)

 

 

 

 

 

 

x

x

x

x

11.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

31

32

33

34

 

 

 

 

 

 

 

 

 

 

 

Объем суммарных поставок каждому заказчику должен удовлетворять

 

 

 

 

 

 

его потребность, т.е. для

переменных xij

(i

1,3;

j

1,4)

 

выполняются

следующие равенства:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

x

9,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

21

31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x12 x22

x32

10,

 

 

 

 

 

 

 

 

(3.8)

 

 

 

 

 

 

 

x23

x33 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

x

11.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

24

34

 

 

 

 

 

 

 

 

 

 

 

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

xij 0, (i

1,3;

j

1,4)

.

 

(3.9)

Таким образом, сформулированная выше задача свелась к задаче нахо-

 

 

 

 

ждения таких значений переменных xij (i

1,3;

j

1,4)

, которые удовлетворяли

бы условиям (3.7) – (3.9) и обеспечивали бы минимальное значение целевой функции (3.6). Такая задача (3.6) – (3.9) называется открытой моделью транспортной задачи (ТЗ).

2. Сведение полученной модели к стандартной транспортной задаче.

Стандартная ТЗ разрешима только в том случае, если выполняется условие баланса:

3

4

 

ai bj .

(3.10)

i 1

j 1

 

23

3

4

В нашей задаче ai 42;

bj 38 , т.е. условие баланса (3.10)

i 1

j 1

нарушено. В связи с этим необходимо привести ТЗ (3.6) – (3.9) к стандартной (закрытой модели), для которой будет выполняться условие баланса (3.10).

3

4

 

 

 

Поскольку ai bj (42 > 38), введем фиктивного потребителя под

i 1

j 1

 

 

 

 

3

4

 

 

номером 5 с объемом потребления b5 ai bj 42 38 4 (тыс. изд.).

 

i 1

j 1

Стоимости перевозок от каждого поставщика этому фиктивному

 

 

потребителю полагаем равными нулю, т.е.

ci5 0; i

1,3

. Получим следующую

стандартную ТЗ:

 

 

 

 

3

4

f cij xij

i 1

j 1

8x11 5x12 3x13 7x14 0x15 2x21 4x22 6x23 3x24

0x25 5x31 2x32 3x33 7x34 0x35 min.

x

x

x

x

x

19,

 

11

12

13

14

15

 

 

x21 x22

x23

x24

x25

12,

x

x

x

x

x

 

11.

 

31

32

33

34

35

 

 

x11 x21

x31

9,

 

 

x

x

x

10,

 

 

 

 

12

22

32

 

 

 

 

 

 

x23

x33

8,

 

 

 

x13

 

 

 

x

x

x

11,

 

 

 

 

14

24

34

 

 

 

 

x

x

x

4.

 

 

 

 

15

25

35

 

 

 

 

(3.11)

(3.12)

(3.13)

xij 0

(i

1,3;

j

1,5)

(3.14)

Решение ТЗ (3.11) – (3.14) будет решением исходной открытой задачи

(3.6) – (3.9).

3. Решение задачи методом потенциалов. 1. Построение начального до-

пустимого опорного плана. Занесем данные задачи (3.11) – (3.13) в табл. 3.1, в которой в верхнем правом углу каждой клетки (i, j), расположенной в i-строке и

j-м столбце, помещены соответствующие стоимости cij (i 1,3; j 1,4) и

найдем начальный опорный план методом минимальной стоимости (можно построить его и другими методами, например, методом минимального элемента).

24

Таблица 3.1

bj

 

 

 

9

 

 

 

 

10

8

 

11

 

4

ai

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

5

3

 

7

 

0

19

 

 

 

 

 

 

 

 

 

8

 

8

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

4

6

 

3

 

0

12

 

 

 

9

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

2

3

 

7

 

0

11

 

 

 

 

 

 

 

 

10

 

 

 

1

 

 

 

 

 

 

 

 

Выпишем начальный опорный план в виде матрицы

 

 

 

 

0

0

8

8

3

 

 

 

 

 

 

X

0

 

9

0

0

3

0

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

0

10

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Стоимость перевозки при этом плане составляет

 

 

 

z

z(X 0 ) 3 8 7 8 0 3 2 9 3 3 2 10 0 1 127 (ден. ед.).

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как число занятых клеток в таблице равно 7,

m n 1 3 5 1 7 , то

этот опорный план является

невырожденным. Пустые клетки начального

опорного плана в табл. 3.1

соответствуют

небазисным переменным xij .

Множество индексов базисных переменных (занятых клеток

табл. 3.1)

Uб 1,1 , 1,2 , 2,2 , 2,3 , 2,4 , 3,4 , 3,5

образует базисное

множество.

Тогда множество индексов небазисных переменных (пустых клеток табл. 3.1) образует небазисное множество

Uн U \ Uб ,

где U – множество индексов всех клеток табл. 3.1.

2. Улучшение допустимого опорного плана. Проверим выполнение критерия оптимальности для построенного опорного плана. Строке 1 (с наибольшим числом базисных клеток) припишем нулевой потенциал u1 0 и, используя

уравнения

ui vj cij , i, j Uб ,

найдем потенциалы всех строк и столбцов для базисных (заполненных) клеток табл. 3.1. Имеем

u1 v3 3, u2 v1 2, u3 v2 2 , u1 v4 7, u2 v4 3, u3 v5 0 , u1 v5 0.

25

Полагая u1 0 , получаем v3 3, v4 7,v5 0, u3 0, v2 2, u2 4, v1 6 . Нулевой потенциал можно выбирать для произвольной строки (столбца). Далее по формулам

ij cij (ui vj ), i, j Uн

подсчитаем оценки небазисных (пустых) клеток и занесем их в левые нижние углы пустых клеток табл. 3.2.

Имеем

11 8 (6 0) 2 0, 22 4 ( 4 2) 6 0, 31 5 (0 6) 1 0 ,12 5 (2 0) 3 0, 23 6 ( 4 3) 7 0, 33 3 (0 3) 0 ,

25 0 ( 4 0) 4 0, 34 7 (0 7) 0 .

Таблица 3.2

 

bj

9

 

 

10

 

8

11

 

 

4

 

ui

 

 

ai

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

5

 

3

┌ ─ 7

 

+

┐0

 

 

 

 

19

2

 

3

 

 

8

ا 8

 

 

ا

3

 

u 0

 

 

 

 

 

 

 

ا

 

 

ا

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– ┌ ─ 2

─ ─ 4

─ ─ 6

+ ┘ 3

 

 

ا

0

 

u2 4

 

 

12

ا 9

 

6

 

 

7

3

 

 

ا

 

 

 

 

 

ا

 

 

 

 

 

4

ا

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ └ ─ 5

 

─ ─ 2

 

─ ─ 3

─ ─ 7

– ┘ 0

 

u3 0

 

 

11

-1

 

 

10

 

 

 

 

 

1

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v j

v1 = 6

 

 

v2 = 2

 

v3 = 3

v4 = 7

v5 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Анализ

оценок

ij ,

i, j Uн показывает, что

данный

опорный план

неоптимален, т.к. среди оценок ij

есть отрицательные. Поэтому переходим к

итерации построения нового плана, лучшего в том смысле, что значение целевой функции (3.6) будет на нем меньше, чем на начальном.

Найдем такую клетку в табл. 3.2, в которой находится минимальное

отрицательное значение оценки ij ,

i, j Uн .

Это будет

клетка

(3, 1).

Следовательно, переменная x31 будет вводиться в базис.

 

 

С помощью

этой клетки и

базисных

клеток

построим

цикл

{ 3,1 , 2,1 , 2,4 , 1,4

, 1,5 , 3,5 , 3,3 } ,

считая

первым

звеном

цикла,

например, вертикальное звено из клетки (3, 1). В табл. 3.2 цикл изображен пунктирами. Помечаем клетки цикла поочередно знаками « + »,« – », причем клетка (3, 1), вводимая в базис, помечается знаком « + ». Среди клеток, помеченных знаком « – », выбираем ту, в которой значение xij минимально. В

26

нашем случае это клетка (3, 5), где x35 1. Обозначим x35 1. Число = 1

добавим к перевозкам в клетках, помеченных знаком « + », и вычтем из величин xij в клетках, помеченных знаком « – ». Объемы xij остальных

перевозок не изменяем. Клетка (3, 5) удаляется из базисного множества, а клетка (3, 1) вводится в базисное множество. В результате получим новый опорный план (табл. 3.3).

Таблица 3.3

bj

9

 

10

 

8

11

4

ui

ai

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

5

 

3

7

0

 

19

2

2

 

 

8

7

4

u1 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4

 

6

3

0

u2 4

12

8

5

 

7

 

4

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

2

 

3

7

0

u3 1

11

1

 

10

 

 

 

 

 

 

 

 

1

 

1

1

 

v j

v1 = 6

 

v2 = 3

 

v3 = 3

v4 = 7

v5 = 0

 

Новый опорный план имеет вид:

 

 

 

0

0

8

7

4

 

X

1

 

8

0

0

4

0

 

 

 

.

 

 

 

1

10

0

0

0

 

 

 

 

 

Общая стоимость перевозок

z1 z(X 1) 3 8 7 7 0 4 2 8 3 4 5 1 2 10 126 (ден. ед.).

Проверим полученное решение на оптимальность. Для этого найдем

потенциалы занятых и оценки свободных клеток по ранее указанному правилу.

Поскольку среди оценок нет отрицательных, то найденный план является

оптимальным.

Выписываем матрицу Х* (без последнего столбца):

27

 

 

 

0

0

8

7

 

X

 

 

8

0

0

4

 

 

 

.

 

 

 

1

10

0

0

 

 

 

 

 

Минимальные суммарные затраты по оптимальному плану составляют zmin 126 (ден. ед.).

Из таблицы 3.3 видно, что избыточная продукция в количестве 4 тыс.

изделий остается на первой фабрике.

Контрольные задания для самостоятельного решения

Задание 3

Постановка задачи. Компания контролирует 4 фабрики,

производительность которых на неделю (в тыс. изделий) задается вектором a = =(а1, а2 , а3 , а4 ). Компания заключила договоры с пятью заказчиками, потребность которых еженедельно (в тыс.изделий) задается вектором b = (b1 , b2 , b3 , b4 , b5). Стоимость транспортировки 1 тысячи изделий j - му заказчику с

i -ой фабрики-изготовителя задается матрицей C Cij 4 5 .

Требуется:

1)составить математическую модель задачи;

2)привести ее к стандартной задаче ( с балансом);

3)построить начальный опорный план;

4)решить задачу методом потенциалов;

5)проанализировать результаты решения.

 

Вариант

 

 

 

 

С

 

 

 

 

 

 

 

 

 

a

b

 

 

 

 

 

 

 

 

3 9

2

1 6

 

1

(50,25,25,15)

(15,20,15,20,30)

 

6 2

4 1 7

 

 

 

 

 

 

 

 

1 4

8 3 2

 

 

 

 

 

 

 

6 2 4

 

 

 

 

 

1 4

 

 

 

 

 

3 5

1 4 2

 

2

(25,20,15,12)

(20,15,20,15,10)

 

5 3

2 1 4

 

 

 

 

 

 

 

 

 

6 3

1 4 2

 

 

 

 

 

 

 

5 7 4

 

 

 

 

 

1 3

 

28

 

 

 

2

6

5

1 4

3

(12,14,13,9)

(8,10,10,7,7)

 

3

6

4 5 6

 

 

 

 

 

 

 

6

7

4

5 6

 

 

 

 

 

6

8

7

2 5

 

 

 

 

 

 

 

 

 

4 6

2

5 3

4

(24,18,114,10)

(18,15,18,13,8)

 

5 3

2 1 4

 

 

 

 

 

 

 

7 4

2

5 3

 

 

 

 

 

 

 

5

7

4

 

 

 

 

1 3

 

 

 

 

5 2

1

5

3

5

(20,13,6,9)

(8,11,7,5,12)

 

8 7

8 1 2

 

 

 

 

 

 

 

6

5

4

3 2

 

 

 

 

 

3 6

1 5

4

 

 

 

 

 

 

 

 

 

5 7

3

6

4

6

(20,14,10,6)

(14,11,14,9,4)

 

6 4

3 2

5

 

 

 

 

 

 

 

9

6

4

7

7

 

 

 

 

 

2

4

6 8 5

 

 

 

 

 

 

 

 

 

4 3

2 1 4

7

(12,14,8,10)

(9,11,7,6,9)

 

6 3

6 2 4

 

 

 

 

 

 

1

5

3

4

6

 

 

 

 

 

4 1

3

5 2

 

 

 

 

 

 

 

 

 

6 8

4

7 5

8

(21,15,11,7)

(15,12,15,10,5)

 

5 3

2 1 4

 

 

 

 

 

 

 

8 5

3 6 4

 

 

 

 

 

 

 

5

7

4

 

 

 

 

1 3

 

 

 

 

4 3

2 1

4

9

(16,11,12,13)

(10,9,8,15,7)

 

7 3

6 2 4

 

 

 

 

 

 

1

4

2 4 6

 

 

 

 

 

4

1

3 5 2

 

 

 

 

 

 

 

 

 

7 9

5

8

6

10

(31,25,21,17)

(25,22,25,20,15)

 

 

 

 

 

 

 

 

8 6

5 4

7

 

 

 

 

6

3

1 4

2

 

 

 

 

 

3 5

7

9

6

 

 

 

 

 

 

 

 

 

3 6

4

5

6

11

(11,13,12,8)

(8,9,9,6,5)

 

2 6

5 1 4

 

 

 

 

 

 

 

6

7

4 5

6

 

 

 

 

 

6

8

7 2

5

 

 

 

 

 

 

29

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]