Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
60
Добавлен:
19.03.2016
Размер:
894.43 Кб
Скачать

 

 

 

 

61

Y = DZ(

 

)

 

X

, i = 1, 2, 3.

i

DSi

 

 

 

Подставляя в данную формулу значения, приведенные в таблице выше,

получаем следующие величины для ценности каждой дополнительной единицы каждого типа сырья: Y1 = 2, Y2 = 0, Y3 = 1. Полученные результаты свидетельст- вуют, что дополнительные вложения в первую очередь следует направлять на закупку сырья S1 и лишь затем - на закупку сырья S3 .

Для ответа на четвертый вопрос опять обратимся к рис. 8. Из него видно, что при небольшом изменении коэффициентов целевой функции с1 и с2 прямая

Z(X ) =17 вращается вокруг точки C. Таким образом, точка C будет оставаться оптимальной до тех пор, пока прямая Z(X ) не выйдет за пределы, определяе- мые 1-м и 3-м ограничениями, т.е. пока вектор-градиент ÑZ (X ) будет нахо- диться между векторами-нормалями A1 и A3 к прямым, соответствующим 1-му и 3-му ограничениям.

Определим сначала диапазон изменения коэффициента с1 при фиксиро- ванном коэффициенте с2 (с2 =4). Тангенсы углов наклона для векторов ÑZ (X ), A1 , A3 соответственно равны:

tg(α ) = c2 / c1 , tg1) = 1, tg3 ) = 2 .

Поэтому диапазон изменения коэффициента с1 в целевой функции определится из соотношения 1 ≤ с2 / c1 ≤ 2 , т.е. 2 ≤ c1 ≤ 4. Диапазон изменения коэффициента с2 при фиксированном значении с1 (с1 =3) определяется аналогично, т.е.

3 ≤ c2 ≤ 6 .

Решение ЗЛП симплекс-методом

Симплекс-метод - универсальный метод решения ЗЛП, разработанный в 1949 г. американским математиком Дж. Данцигом. В его основе лежит идея по- следовательного улучшения полученного решения, для реализации которого необходимы [9]:

1)способ определения какого-либо первоначального допустимого ба- зисного решения задачи;

2)правило перехода к лучшему (по крайней мере, не худшему) реше- нию;

3)критерий проверки оптимальности найденного решения.

62

Для использования симплекс-метода ЗЛП должна быть приведена к кано- ническому (стандартному) виду, т.е. система ограничений ее модели (за исклю- чением тривиальных, или естественных ограничений) должна быть представле- на в виде уравнений (равенств). Переход к канонической ЗЛП (КЗЛП) выполня- ется следующим образом [20]:

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

ния

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

xk

(k =

 

, которые одновременно входят в целевую функцию с

n,n + m)

коэффициентом 0, т.е. не оказывают влияния на ее значение;

-переменные, на которые не наложено условие неотрицательности, представляются в виде разности двух новых неотрицательных пере-

менных, т.е. xj = xj - xj , (xj ³ 0, xj ³ 0).

Ниже рассмотрен пример использования симплекс-метода для решения задачи производственного планирования как частного случая общей задачи ли- нейного программирования (ОЗЛП).

Пример решения задачи производственного планирования симплекс-методом

Решим симплекс-методом задачу, математическая модель которой пред- ставлена соотношениями (3.12) - (3.17). Прежде всего, преобразуем с помощью дополнительных переменных неравенства (3.13) - (3.15) в уравнения. В данном случае все дополнительные переменные (их число равно 3) вводятся со знаком «плюс», так как рассматриваемые неравенства представляют собой неравенства вида «меньше или равно».

Получим систему ограничений в виде:

ìx

+ x

2

+ x

3

= 5

 

 

ï 1

 

 

 

 

 

 

 

í2x1 + x2

+ x4

=

9

(3.18)

ïx

+ 2x

2

+ x

5

=

7.

 

î 1

 

 

 

 

 

 

 

Для нахождения первоначального базисного решения разобьем перемен-

ные на две группы - основные и не основные (вспомогательные).

Для этого

можно воспользоваться следующим простым правилом для определения основ-

ных переменных [20]: в качестве основных переменных на первом шаге следует выбрать (если это возможно) такие m переменных, каждая из которых вхо-

63

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

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

I шаг. Основные переменные: x3 , x4 , x5 . Вспомогательные переменные: x1 , x2 .

В системе уравнений (3.18) выразим основные переменные через вспомо- гательные, т.е. получим следующую систему уравнений:

ì

 

 

= 5 - x1 - x2

 

 

 

 

ïx3

 

 

 

 

íx4

= 9- 2x1 - x2

 

(3.19)

ïx

5

= 7 - x

- 2x

.

 

 

 

î

 

1

2

 

 

 

 

Приравняв вспомогательные переменные к нулю, т.е. положив x1 = 0 и

x2 = 0 , получим допустимое базисное решение

 

 

1 = (0;0;5;9;7). Это решение

X

соответствует совпадающей с началом системы координат O (0; 0) вершине A пятиугольника ABCDE, который представляет ОДР задачи (см. рис. 7 выше).

Проверим, не является ли полученное допустимое базисное решение оп-

тимальным. Для

этого

рассмотрим выражение

для

целевой

функции

Z(

 

) = 3x1 + 4x2 ,

которое

для найденного решения

 

 

1

равняется

значению

X

X

функции Z1 = Z(X1) = 0. Очевидно, что значение функции Z(X ) можно увели- чить путем увеличения любой из вспомогательных переменных, входящих в выражение для Z с положительным коэффициентом. Это можно осуществить, перейдя к новому допустимому базисному решению, в котором эта переменная будет основной, т.е. сможет принимать не нулевое, а положительное значение (если новое решение окажется вырожденным, то целевая функция сохранит свое значение). При таком переходе одна из основных переменных становится вспомогательной, что в геометрической трактовке означает переход к соседней вершине пятиугольника ABCDE, в которой значение целевой функции «лучше» (по крайней мере «не хуже»). В данном случае для увеличения значения Z можно переводить в основные переменные либо x1 , либо x2 , так как обе эти пе- ременные входят в выражение для Z со знаком «плюс». Выберем переменную x2 , так как ее коэффициент больший, чем у x1 .

64

Система уравнений (3.13) - (3.15) накладывает ограничения на рост пе- ременной x2 . Поскольку все переменные должны оставаться неотрицательны- ми, выполняются следующие неравенства (при этом x1 = 0 как вспомогательная переменная):

 

 

 

ì

 

 

 

ìx3

= 5 - x2 ³ 0

 

ïx2 £ 5

 

ï

 

откуда

ï

£ 9

(3.20)

íx4 = 9- x2 ³ 0

íx2

ï

= 7 - 2x2 ³ 0

 

ï

 

7

 

îx5

 

ï

= 3,5.

 

 

 

£ 2

 

 

 

îx2

Очевидно, что сохранение условия неотрицательности всех переменных (допустимость решения) возможно, если не нарушается ни одна из границ, по- лученных для всех уравнений (3.20). В данном случае возможное наибольшее

значение для переменной x2

определяется как x2 = min{5;9; 3,5} = 3,5. При

x2 = 3,5 основная переменная x5

обращается в нуль и переходит во вспомога-

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

II шаг. Основные переменные: x2 , x3 , x4 . Вспомогательные переменные: x1 , x5 .

Выразим новые основные переменные через новые вспомогательные пе-

ременные, начиная с разрешающего уравнения:

 

 

 

 

 

ìx

2

= (7 - x - x

5

) / 2

 

 

 

 

 

 

ï

 

1

 

 

 

 

 

 

 

íx3

= 5 - x1 - (7 - x1 - x5 ) / 2

(3.21)

ïx

4

= 9 - 2x - (7 - x - x

5

) / 2

 

 

 

î

 

1

 

1

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

ì

 

 

= 3,5 - 0,5x1 - 0,5x5

 

 

 

 

 

ïx2

 

 

 

 

 

íx3

=1,5 - 0,5x1 + 0,5x5

 

 

 

 

(3.22)

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 5,5 -1,5x1 + 0,5x5 .

 

 

 

 

 

îx4

 

 

 

 

 

Второе базисное решение

 

2 = {0; 3,5;1,5; 5,5; 0} является допустимым

X

и соответствует вершине B(0; 3,5) на рис. 7. Переход от решения

 

1

к решению

X

X2 геометрически можно интерпретировать как переход от вершины А(0; 0) к вершине B(0; 3,5) пятиугольника ABCDE.

65

Проверим, не является ли оптимальным полученное решение. Запишем целевую функцию, выраженную через вспомогательные переменные, исполь-

зующиеся на данном шаге:

 

 

 

Z2 = Z(

 

2 ) = 3x1 + 4x2 = 3x1

+ 4(3,5 - 0,5x1 - 0,5x5 ) =

 

X

(3.23)

= 3x1 +14 - 2x1

- 2x5 = x1 - 2x5 +14 =14.

 

Убедимся, что полученное решение лучше,

чем предыдущее (т.е. значение Z

увеличилось): Z1 = Z2 Z1 = 14 − 0

= 14 .

 

 

Вычисленное значение Z (X 2 ) не является максимальным, поскольку в выражении для целевой функции (3.23) имеется переменная x1 с положитель- ным коэффициентом, равным 1. Система уравнений (3.22) накладывает ограни- чения на рост переменной x1 . Обеспечивая выполнение условия неотрицатель-

ности переменных, запишем следующие неравенства:

ìx

2

= 3,5 - 0,5x

- 0,5x

5

³ 0

ï

1

 

 

 

íx3 =1,5 - 0,5x1

+ 0,5x5 ³ 0

ïx

4

= 5,5 -1,5x

+ 0,5x

5

³ 0,

î

1

 

 

откуда, полагая x5 = 0 , получим

 

 

 

 

ì

 

 

 

 

 

 

 

 

ïx

£ 7

 

 

 

 

 

ï

1

 

 

 

 

 

 

 

íx1 £ 3

 

 

 

 

ï

 

11

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 3 .

 

 

 

 

îx1

 

 

 

 

(3.24)

(3.25)

Полученная система неравенств (3.25) определяет наибольшее значение для переменной x1 : x1 = min{7; 3;11/ 3} = 3. При x1 = 3 переменная x3 обращает- ся в нуль, поэтому ее предполагается перевести во вспомогательные перемен- ные. Второе уравнение системы (3.22) является разрешающим, переменная x3 переходит во вспомогательные переменные, x1 - в основные.

III шаг. Основные переменные: x1 , x2 , x4 . Вспомогательные переменные: x3 , x5 .

Как и на шаге II, выразим новые основные переменные через новые вспомогательные переменные, начиная с разрешающего уравнения.

0,5x1 = 1,5 − x3 + 0,5x5 ,

x1 = 3 − 2x3 + x5 .

(3.26)

Подставляем выражение для x1 в следующее уравнение:

 

 

 

 

 

 

 

 

 

 

66

 

 

 

 

 

 

 

x2

= 3,5 − 0,5x1 − 0,5x5 = 3,5 − 0,5(3 − 2x3 + x5 ) − 0,5x5

=

 

 

 

 

= 3,5 -1,5 + x3 - 0,5x5 - 0,5x5 = 2 + x3 - x5 ,

 

 

 

 

 

 

 

 

 

 

x2 = 2 + x3 x5 .

 

 

(3.27)

Подставляем выражение для x1 в следующее уравнение:

 

 

 

 

x4

= 5,5 −1,5x1 + 0,5x5

= 5,5 −1,5(3 − 2x3 + x5 ) + 0,5x5

=

 

 

 

 

= 5,5 - 4,5 + 3x3 -1,5x5 + 0,5x5 =1+ 3x3 - x5 ,

 

 

 

 

 

 

 

 

 

 

x4

= 1+ 3x3 x5 .

 

(3.28)

В результате преобразований (3.26) - (3.28) получаем следующую систе-

му уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ìx

= 3 - 2x

3

+ x

5

 

 

 

 

 

 

 

 

ï 1

 

 

 

 

 

 

 

 

 

 

 

 

 

íx2

= 2 + x3 - x5

 

 

(3.29)

 

 

 

 

 

 

ï

=1 + 3x3 - x5.

 

 

 

 

 

 

 

 

 

 

îx4

 

 

Базисное решение

 

 

3 = (3; 2; 0;1; 0)

 

соответствует вершине C(3; 2) пяти-

X

 

угольника ABCDE, представляющего ОДР задачи (см. рис. 7 выше). Выражаем

целевую функцию через вспомогательные переменные:

 

 

 

 

Z3 = Z (

 

3 ) = 3x1 + 4x2 = 3(3 - 2x3 + x5 ) + 4(2 + x3 - x5 ) =

 

X

(3.30)

 

 

 

 

 

 

 

= 9 - 6x3 + 3x5 + 8 + 4x3 - 4x5 =17 - 2x3 - x5.

 

 

 

 

 

 

 

 

Полагая

x3 = 0

и x5 = 0 ,

вычислим значение целевой функции

Z3 = Z (

 

3 ) =17.

Убедимся, что полученное решение лучше, чем предыдущее:

X

Z2 = Z3 Z2 = 17 −14 = 3, т.е. значение целевой функции возросло по сравне- нию со значением, вычисленным на шаге II. Кроме того, третье допустимое ба- зисное решение является оптимальным, так как коэффициенты вспомогатель- ных переменных в выражении для целевой функции (3.30) являются отрица- тельными, т.е. рост значений x3 и x5 не приведет к увеличению значения Z .

Таким образом, оптимальным решением задачи является x1 = 3 и x2 = 2 , обеспечивающие оптимальное значение целевой функции Z (x1, x2 ) = 17 . С уче-

том экономического смысла рассмотренной задачи можно сделать следующий вывод: для получения максимальной прибыли, равной 17 денежным единицам, предприятие должно выпускать и реализовать 3 единицы продукта P1 и 2 еди- ницы продукта P2 . При этом дополнительные переменные x3 , x4 и x5 показы- вают разницу между запасами сырья (ресурсов) каждого вида и их потреблени- ем, т.е. остатки ресурсов. При выполнении оптимального плана производства

67

x3 = x5 = 0 , т.е. сырье S1 и S3 расходуется полностью, а остаток сырья S2 равен

1единице.

Взаключение сформулируем критерий оптимальности решения при оты- скании максимума линейной функции, а также определим возможные подходы к отысканию минимума линейной функции.

Если при отыскании максимума линейной функции в ее выражении, со- держащем вспомогательные переменные, отсутствуют положительные коэф- фициенты, полученное решение оптимально.

При поиске минимума линейной функции F можно использовать сле- дующие два пути [9]:

1) отыскать максимум функции Z , полагая Z = −F и учитывая, что

Fmin = −Zmax ;

2)модифицировать симплекс-метод: на каждом шаге уменьшать линей- ную функцию за счет той вспомогательной переменной, которая вхо- дит в выражение линейной функции с отрицательным коэффициен- том.

68

Глава 4. Транспортная задача

К ЗЛП транспортного типа (кратко: транспортной задаче ТЗ) прихо- дят при рассмотрении различных практических ситуаций, связанных с состав- лением наиболее экономичного плана перевозок продукции, управления запа- сами, назначением персонала на рабочие места, оборотом наличного капитала и многими другими.

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

Информация, необходимая для использования модели ТЗ, включает сле- дующее [21]:

1.Список пунктов отправления (ПО) и пропускная способность каждого из них или количество поставок за определенный период.

2.Список пунктов назначения (ПН) и их показатели спроса за определен- ный период.

3.Стоимость транспортировки единицы товара из каждого пункта от-

правления в каждый пункт назначения.

Эта информация представляется в виде так называемой транспортной табли-

цы (рис. 9).

транспортировкиСтоимость единицыоднойтовара из отправленияпункта A

назначенияпунктB

в

 

1

 

 

1

ПН

B1

B2

B3

B4

B5

Запасы

 

 

ПО

пункте

 

 

 

 

 

 

 

 

A1

2

3

4

2

4

140

3

A2

8

4

1

4

1

180

в A

товараЗапас

отправления

A3

9

7

3

7

2

160

 

 

Потреб-

60

70

120

130

100

480

 

 

ности

 

 

 

 

 

 

 

 

Потребность в товаре в пункте назначения B2

Рис. 9. Пример транспортной таблицы

69

Ниже рассмотрена общая постановка ТЗ и выполнено построение ее ма- тематической модели [16].

Экономико-математическая модель ТЗ

Постановка задачи. Некоторый однородный товар (продукт, груз), нахо- дящийся у m поставщиков Ai в количестве ai единиц (i = 1,2,...,m) необходимо доставить n потребителям Bj в количестве bj единиц ( j = 1,2,...,n). Известна стоимость cij перевозки единицы товара от i -го поставщика к j -му потребите-

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

Построение математической модели. Обозначим через xij количество единиц товара, запланированных к перевозке от i -го поставщика к j -му по- требителю. Тогда математическая модель ТЗ формулируется следующим обра- зом:

 

 

 

m

n

 

Z (

 

) = ååcij xij ® min

(4.1)

X

 

 

 

i=1 j=1

 

при ограничениях

 

 

 

n

 

 

 

åxij

£ ai

, (i = 1,2,...,m) ;

(4.2)

j=1

 

 

 

m

 

 

 

åxij

³ bj

, ( j = 1,2,...,n);

(4.3)

i=1

 

 

 

xij ³ 0, (i = 1,2,...,m; j = 1,2,...,n ).

(4.4)

Ограничения (4.2) означают, что суммарный объем перевозок от i -го по- ставщика не может превышать имеющегося у него запаса товара. Ограничения (4.3) означают, что суммарные перевозки товара j -му потребителю должны полностью удовлетворить его потребности в товаре. Ограничения (4.4) исклю- чают обратные перевозки.

Из ограничений (4.2) и (4.3) следует, что

m

n

åai ³ åbj .

i=1

j=1

Если имеет место равенство

70

m

n

 

åai = åbj ,

(4.5)

i=1

j=1

 

то модель называется сбалансированной транспортной моделью. В сбаланси-

рованной модели ограничения (4.2) − (4.3) имеют вид равенств. В реальных ус- ловиях запас товара не всегда равен спросу (потребности), но транспортную модель всегда можно сбалансировать.

 

 

 

 

m

n

В случае превышения запаса над спросом, т.е. если åai >

åbj , вводится

 

 

 

 

i=1

j=1

фиктивный (n + 1) -й потребитель со спросом

 

 

 

 

m

n

 

 

 

bn+1 = åai

åbj ,

 

 

 

 

i=1

j=1

 

 

а соответствующие стоимости ci,n+1

(i = 1,2,...,m)

считаются равными нулю.

m

n

 

 

 

 

Аналогично, при åai < åbj

вводится фиктивный (m +1) -й поставщик с

i=1

j=1

 

 

 

 

запасом товара

 

 

 

 

 

 

 

n

m

 

 

 

am+1 = åbj åai ,

 

 

 

 

j=1

i=1

 

 

а соответствующие стоимости cm+1, j

( j = 1,2,...,n)

считаются равными нулю.

Таким образом, исходная задача сводится к сбалансированной ТЗ, из оп-

тимального плана которой получается оптимальный план несбалансированной ТЗ.

Примечание 1. Несбалансированную ТЗ называют также открытой ТЗ, тогда как сбалансированную ТЗ закрытой ТЗ.

Примечание 2. Стремление сбалансировать ТЗ обусловлено возможно- стью применить в этом случае эффективный вычислительный метод.

Ниже приведен пример сбалансированной ТЗ и дано ее решение с под- робными пояснениями [16].

Задача № 11. На трех базах (пунктах отправления) A1, A2 , A3 находится однородный груз в количествах, соответственно равных 140, 180 и 160 едини- цам. Этот груз требуется перевести в пять пунктов назначения B1, B2 , B3, B4 , B5 соответственно в количествах 60, 70, 120, 130 и 100 единиц. Стоимости пере- возки единицы груза из каждого пункта отправления в соответствующие пунк- ты назначения указаны в табл. 10.