Материал: Методические указания для организации выполнения курсовой работы по курсу Высшая математика. Пантелеев И.Н

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

Базис

СБ

P0

2

3

0

0

0

 

P1

P2

P3

P4

P5

P6

 

 

 

 

1

P2

3

1

0

1

2/5

0

-1/5

1/5

2

P4

0

5

0

0

0

1

1

-1

3

P1

2

1

1

0

-1/5

0

-2/5

2/5

4

 

 

5

0

0

4/5

0

-7/5

7/5

Последняя строка снова содержит отрицательное число. В базис вводим вектор P5 , ис-

ключаем P4 .

 

Базис

СБ

P0

2

3

0

0

0

 

P1

P2

P3

P4

P5

P6

 

 

 

 

1

P2

3

2

0

1

2/5

1/5

0

0

2

P5

0

5

0

0

0

1

1

-1

3

P1

2

3

1

0

-1/5

2/5

0

0

4

 

 

12

0

0

4/5

7/5

0

0

В 4-й строке последней симплексной таблицы нет отрицательных чисел. Значит найденный опорный план X * = (3, 2, 0, 0, 5, 0) является оптимальным. Значение целевой функции zmax = 12 .

Составим двойственную задачу.

Умножим третье ограничение на-1, тогда все неравенства будут содержать знак«≤». Задача примет вид исходной задачи симметричной пары 1:

z = 2x1 + 3x2 ® max

ì - x1 + 2x2 £ 1, y1

ï

+ x2

£ 8, y2

í 2x1

ï

- x2

£ -3, y3

î- 2x1

x1 , x2 ³ 0.

Число переменных в двойственной задаче равно числу ограничений в исходной задаче, т.е. трём: y1 , y2 , y3 .

Умножим правые части ограничений на соответствующие переменные двойственной задачи и сложим их, получим целевую функцию: g = y1 + 8y2 - 3y3 .

Целевая функция исходной задачи исследуется на максимум, следовательно, целевая функция двойственной задачи исследуется на минимум.

Матрица системы ограничений исходной задачи имеет вид:

 

æ -1

2

ö

 

A = ç

2

1

÷ . Транспони-

 

 

ç

- 2

 

÷

 

 

 

è

-1ø

 

руем её и получим аналогичную матрицу двойственной задачи - A

T

æ-1

2

- 2

ö

 

= ç

2

1

-1

÷ . Правыми

 

 

è

ø

частями в ограничениях двойственной задачи являются коэффициенты при

неизвестных в

целевой функции исходной задачи.

 

 

 

 

 

 

11

Окончательно двойственная задача имеет следующий вид:

g = y1 + 8 y2

- 3y3 ® min .

ì- y1 + 2 y

2 - 2 y3 ³ 2,

ï

2 y1 + y2 - y3 ³ 3,

í

ï

y1 , y2

, y3 ³ 0.

î

Найдём её решение, используя теоремы двойственности. По первой теореме двойст-

венности

оптимальные решения исходной и двойственной задач равны, следовательно,

g min = zmax

= 12 .

Из соотношений второй теоремы двойственности следует, что если какое-то ограничение исходной задачи выполняется в виде строгого неравенства, то соответствующая двойственная оценка равна нулю. Подставим найденное оптимальное решение в систему ограничений исходной задачи:

-1* 3 + 2 * 2 = 1 = 1, 2 * 3 +1* 2 = 8 = 8, 2 * 3 +1* 2 = 7 > 3.

Третье ограничение выполняется в виде строгого неравенства, следовательно, y3* = 0 .

Если некоторая компонента xi* оптимального плана исходной задачи отлична от нуля, то соответствующее ограничение двойственной задачи выполняется в виде равенства. В нашем примере и x1* = 3 ¹ 0 , и x2* = 2 ¹ 0 , следовательно, оба ограничения двойственной задачи выполняются в виде равенства.

ì- y1 + 2 y2 - 2 y3 = 2,

Учитывая, что

*

= 0 , получим:

ì- y1 + 2 y2 = 2,

í

2 y + y

2

- y

3

= 3.

y3

í

2 y

1

+ y

2

= 3.

î

1

 

 

 

 

 

î

 

 

 

Решив систему, получим y1 = 4 / 5, y2 = 7 / 5.

Окончательно Y * = (4 / 5, 7 / 5, 0), g min = 12.

Решение двойственной задачи можно получить другим способом, используя формулу

Y * = CБ P-1 .

Матрица P-1 находится в последней симплексной таблице. Ее столбцы расположены под столбцами единичной матрицы, образующими базис начального опорного решения, т. е. под векторами P3 , P4 , P6 (именно для этой цели мы продолжали вычислять вектор P6 ):

æ 2 / 5 1/ 5 0

ö

æ3* 2 / 5 + 0 * 0 + 2 * (-1/ 5)

ö

ç

 

÷

ç

 

÷

Y * = (3 0 2)ç 0

1 -1÷ = ç 3 *1/ 5 + 0 *1 + 2 * 2 / 5

÷ = (4 / 5 7 / 5 0).

ç

 

÷

ç

3 * 0 + 0 * (-1) + 2 * 0

÷

è-1/ 5 2 / 5 0

ø è

ø

Ответ: X * = (3, 2) , zmax

= 12 ; Y *

= (4 / 5, 7 / 5, 0), g min = 12.

 

12

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

Пусть имеется m поставщиков А1, А2, ..., Аm однородного груза в количествах соответственно а1, а2, .., .аm единиц и n потребителей В1, В2, ..., Вn этого груза, потребность которых составляет соответственно b1, b2 ..., bn единиц.

Известны стоимости перевозок (тариф) единицы груза от i-го поставщика к j-му потре-

бителю - сij (i=1,m; j=1,n).

Требуется составить такой план перевозок, который обеспечит минимальные транспортные расходы.

Возможны три ситуации:

1) количество груза у всех поставщиков равно потребности в данном грузе всех потребителей:

m

n

a1 + a2 + ... + am = b1 + b2 + ... + bn или åai

= åb j .

i=1

j =1

2) количество груза у всех поставщиков больше потребности в данном грузе всех -по требителей:

 

 

 

 

 

m

n

a1 + a2

+ ... + am

> b1

+ b2

+ ... + bn

или åai

> åb j .

 

 

 

 

 

i=1

j =1

3) количество груза у всех поставщиков меньше потребности в данном грузе всех по-

требителей:

 

 

 

 

 

 

 

 

 

 

 

m

n

a1 + a2

+ ... + am

< b1

+ b2

+ ... + bn

или åai

< åb j .

 

 

 

 

 

i=1

j=1

В первом случае модель задачи называется закрытой, во втором и третьем – открытой. Теорема. Для разрешимости транспортной задачи необходимо и достаточно, чтобы за-

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

m n

е. чтобы выполнялось равенство åai = åb j .

i=1 j =1

В случае превышения запаса над потребностью вводится фиктивный (n + 1)-й пункт на-

m n

значения с потребностью bn+1 = åai - åb j и соответствующие тарифы считаются равными

i =1 j =1

нулю.

 

 

Аналогично

вводится фиктивный(m + 1)-й пункт отправления с запасом груза

n

m

 

am+1 = åb j

- åai

и тарифы полагаются равными нулю. Этим задача сводится к закрытой

j =1

i =1

 

транспортной задаче, из оптимального плана которой получается оптимальный план исходной задачи.

Решение транспортной задачи включает следующие этапы:

1. Нахождение первоначального опорного плана(метод северо-западного угла, метод минимальной стоимости). При этом число заполненных клеток должно быть равно m+n-1.

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

Заполнение начинается с левой верхней клетки (северо-западный угол). На каждом шаге, исходя из запасов очередного поставщика и запросов очередного потребителя, заполняется только одна клетка и исключается из рассмотрения один поставщик или потребитель. Если в очередную клетку таблицы требуется поставить перевозку, а поставщик или потреби-

13

тель имеют нулевые запасы или запросы, то в клетку ставится перевозка, равная нулю (базисный нуль), и после этого исключается из рассмотрения поставщик или потребитель.

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

2. Проверка опорного плана на оптимальность, например, методом потенциалов.

Пример. Четыре предприятия используют три вида сырья. Потребности в сырье каждого из предприятий соответственно равны100, 90, 170 и 30 ед. Сырьё сосредоточено в трёх пунктах, а запасы соответственно равны 200, 160, и 140 ед. Тарифы перевозок заданы матрицей

æ12

15

21

14

ö

С = ç13

8

15

10

÷ .

ç

16

26

12

÷

è19

ø

Составить такой план перевозок, при котором общая стоимость перевозок является минимальной?

Данная задача является открытой, так как потребности в сырье100+90+170+30=390 меньше запасов 200+160+140=500. Введём 5-го фиктивного потребителя с потребностью b5 = 500 - 390 = 110 .

Теперь исходные данные задачи запишем в виде таблицы, а опорный план получим методом северо-западного угла.

Заполнение таблицы начинаем с клетки(1,1). х11 = min(a1=200,b1=100)=b1=100 – запасы А1 позволяют полностью удовлетворить потребности пункта В1, значит исключаем этого потребителя из рассмотрения. Теперь запасы пункта A1 считаем равными a1=200-100=100 ед. В оставшейся части таблицы левой верхней клеткой является (1,2): х12 = min(a1,b2)=b2=90

– снова запасы удовлетворяют потребность полностью. Внесём значение в соответствующую клетку и исключим из рассмотрения столбец В2. Запасы пункта А1 считаем равными a1=10090=10 ед. Теперь «северо-западным углом» является клетка (1,3). х13 = min(a1,b3)=a1=10 – запасы могут удовлетворить потребность пунктаB3 частично. Заполняем клетку (1,3) и исключаем из рассмотрения строкуA1. Потребности пункта B3 считаем равными b3=17010=160. х23 = min(a2,b3)=a2=b3=160 – запасы A2 исчерпаны, потребность B3 удовлетворена. Но по правилам мы не можем вычеркнуть и строку и столбец одновременно. Поэтому исключим из рассмотрения сначала столбец B3, а в клетку (2,4) запишем х24 =0 (так как запасы А2 уже исчерпаны) и только теперь вычеркнем строку А2. И так далее. Получим следующую таблицу.

 

Потребности

b1=100

b2=90

b3=170

b4=30

 

b5=110

Запасы

 

β1=12

β2=15

β3=21

β4=16

 

β 5=4

a1=200

α1=0

12

100

15

90

21

10

14

0

0

 

a2=160

α2=-6

13

 

8

 

15

160

10

0

 

a3=140

α3=-4

19

 

16

 

26

 

12

30

0

110

Число заполненных клеток равно7 и m+n-1=3+5-1=7 – план невырожденный. Оптимальный план найдём методом потенциалов.

14

Теорема. В оптимальном плане транспортной задачи заполненным клеткам отвечают

равенства αi + β j = cij , а пустым неравенства αi

+ β j

£ cij .

Расставим потенциалы:

 

 

 

ì α1

 

 

 

ì α1 + β1 = с11 = 12

 

 

 

= 0

 

 

 

 

 

ïα2 = -6

 

ïα + β = с = 15

 

 

 

 

ïα + β = с = 21

 

 

 

ïα3 = -4

 

 

 

1

 

2

12

 

 

 

 

ï

 

 

 

 

 

ï

 

1

 

3

13

= 15 , положим α1 = 0 , тогда

ï

β

= 12

.

íα2 + β3

= с23

í

 

1

 

 

ïα

 

+ β

 

= с

 

= 10

 

 

 

ï

β2 = 15

 

ï

α

2

+ β

4

= с

24

= 12

 

 

 

ï

β3

= 21

ï

3

4

34

 

 

 

ï

β

 

= 16

 

 

 

 

 

 

 

 

 

 

îα3 + β5 = с35 = 0

 

 

 

ï

 

4

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

î

β5

 

 

(Обычно равным нулю принимают потенциал строки или столбца с наибольшим чис-

лом заполненных клеток.)

 

 

 

 

 

 

 

 

 

Теперь проверим пустые клетки на выполнение неравенства αi + β j £ cij .

 

 

 

 

 

 

 

ìα1 + β4 = 0 +16 = 16 > 14;

 

14 = 16 -14 = 2

 

 

 

 

 

 

 

ïα + β

5

= 0 + 4 = 4 > 0;

 

15

= 4 - 0 = 4

 

 

 

 

 

 

 

ï 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

+ β = -6 +12 = 6 < 13

 

 

 

 

 

 

 

 

 

ï 2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

+ β2 = -6 +15 = 9 > 8;

22 = 9 - 8 = 1 .

 

 

 

 

 

 

 

íα2

 

 

 

 

 

 

 

ïα2 + β5 = -6 + 4 = -2 < 0

 

 

 

 

 

 

 

 

 

ïα3 + β1 = -4 +12 = 8 < 19

 

 

 

 

 

 

 

 

 

 

ïα3

+ β2

= -4 +15 = 11 < 16

 

 

 

 

 

 

 

 

 

ïα

+ β

3

= -4 + 21 = 17 < 26

 

 

 

 

 

 

 

 

 

î 3

 

 

 

 

 

 

 

 

Для клеток (1,4), (1,5), (2,2) неравенство не выполняется, значит опорный план не является оптимальным. В одну из этих клеток нужно "ввезти" груз. Выбираем ту, для которой разница Dij максимальна, т. е. в (1,5). Строим цикл.

Цикл перерасчёта таблицы - это последовательность ячеек, начинающаяся и заканчивающаяся в одной и той же клетке, с вершинами, лежащими в занятых клетках, кроме одной.

Вершина цикла – клетка, в которой происходит поворот под прямым углом . "Перемещаем" груз по следующим правилам:

1.каждой из клеток, связанных циклом присваивается знак: пустой ячейке "+", остальным - поочерёдно знаки "-" и "+" .

2.среди минусовых клеток находим числоx = min( xij ) и прибавляем его к числам,

стоящим в плюсовых клетках, и вычитаем из чисел, стоящих в минусовых клетках; остальные клетки вне цикла остаются без изменения.

В нашем примере цикл образуют шесть ячеек: (1,5) – пустая, для которой не выполня-

ется неравенство, и (3,5), (3,4), (2,4), (2,3), (1,3) – заполненные.

х = min(10,0,110)=0. Значит в плюсовые клетки"завозим" 0 ед. груза, из минусовых "вывозим". Получим новый опорный план:

 

Потребности

b1=100

b2=90

b3=170

b4=30

 

b5=110

Запасы

 

β1=12

β2=15

β3=21

β4=12

 

β 5=0

a1=200

α1=0

12

100

15

90

21

10

14

 

0

0

a2=160

α2=-6

13

 

8

 

15

160

10

30

0

110

a3=140

α3=0

19

 

16

 

26

 

12

0

 

 

 

 

 

15

 

 

 

 

 

 

Источник: https://studfile.net/preview/16565055/