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

книги из ГПНТБ / Крылов Н.В. Организация и планирование кислородного производства [учеб. пособие]

.pdf
Скачиваний:
17
Добавлен:
25.10.2023
Размер:
16.71 Mб
Скачать

жать самое большее 1 кг вещества А. Ограничение, что готовая смесь должна содержать 2 кг вещества В, выражается равен­ ством В —2 кг и изображается линией в. Неравенство вещества В 2 кг на рис. 16 заштриховано площадью над линией в. Область (площадь затушевана и заштрихована), одновременно удовлетворяющая двум отдельным условиям, и называется об­ ластью допустимых решений. В этой области каждая точка со­ ответствует не менее 1 кг вещества Л и 2 кг вещества В. В точ­ ке пересечения линий а и в одновременно удовлетворяются оба равенства (Л = 1 кг, В = 2 кг).

 

Вещество Агкг

Рис. 15. Линейные ограничения,

Рис. 16. Линейные ограничения, вы­

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

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

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

Рассмотрим пример линейного программирования, включаю­ щий в себя два ограничения и три переменные. Ставится зада­ ча составления смеси, которая приготовляется из двух компо­ нентов (переменные величины). Ограничения задачи состоят в том, что готовая смесь должна содержать не менее 2 кг ве­ щества В и 1 кг вещества А. Требуется при соблюдении этих ограничений получить смесь, имеющую минимальную стои­ мость.

При решении задачи можно использовать следующие три утверждения (постулата), которые будут справедливыми:

1. Наибольшее число компонентов, которое может быть ис­ пользовано, равно числу ограничений. В нашем примере с дву­ мя ограничениями будет использовано не более двух компо­ нентов, хотя выбор компонентов и может быть расширен.

189

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

3.Отрицательные количества компонента недопустимы и

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

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

 

Эта

задача

может быть решена

графически. На рис. 17

по оси ординат отложены количества

компонента I, а по

оси

абсцисс

компонента — II. Ограничения

по веществу В выпол­

няются,

если в

смеси будет или 5 кг

компонента II или 3

кг

0 1 г

3 4 5 6 7

Г---1-- 1---

5 6 7

П 1

? J 4

Компонент И, нг

 

 

Компонент 11, нг

Рис. 17. Ограничения по ве­

Рис.

18.

Область

допустимых

ществу А и В .

 

 

решений.

компонента

I. Ограничения

по веществу

А

выполняются, если

в смеси будет 5 кг компонента I или 4 кг компонента II. Соеди­ нив точки с координатами (5) и (4), а также точки с координа­ тами (3) и (5), получим две прямые. Линия ch и область пра­ вее ее удовлетворяют требованию о содержании в готовой сме­ си не менее 1 кг вещества А, а линия fe и область правее ее удовлетворяют требованию о содержании в готовой смеси не менее 2 кг вещества В. В обеих этих областях лежат допусти­ мые решения. На рис. 18 рассмотрены оба ограничения. Общее решение задачи должно удовлетворять обоим ограничениям одновременно. При этом точка, соответствующая решению за­ дачи, может находиться справа от обеих линий или на одной из них и быть при этом справа от другой или на пересечении двух линий. Вся эта область (частично заштрихована на гра­ фике) является областью допустимых решений. В нашем слу­ чае она представляет собой выпуклую многоугольную фигуру, некоторые стороны которой уходят в бесконечность. В общем случае область допустимых решений является выпуклым мно­ гогранником.

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

190

ка. Применительно к рассматриваемой задаче это означает, что точка, соответствующая ее оптимальному решению (наимень­ шая стоимость смеси), лежит в вершинах с, d или е. Для слу­ чая, когда оптимальному решению отвечают две соседние вер­

шины (например, с и d),

все точки,

лежащие на прямой cd, бу­

дут соответствовать оптимальному

решению. Точка, соответ­

ствующая оптимальному

решению,

будет лежать в вершине

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

Произведем оценку трех вершин. Вершина с указывает, что в решении не используется компонент II, а принимается 5 кг компонента I. Килограмм компонента I стоит 2 руб. Следова­ тельно, общая стоимость смеси составит 5X 2=10 руб. Верши­ на е соответствует решению, в котором предусматривается ис­ пользование компонента II в количестве 5 кг. Отсюда при стоимости 1 кг компонента 1,5 руб. общая стоимость смеси со­ ставит 7,5 руб. Вершина d изображает смесь I и II компонен­ тов. Для расчета в точке d количества компонентов I и II, используемых для приготовления смеси, решим систему двух уравнений, определяющих координаты точки пересечения. Каж­

дое уравнение имеет вид: у — ах + Ъили г/= — (ограниче-

 

о fx

 

 

ние по веществу Л); у = — Qg +of (ограничение по веществу В).

Совместное решение этих уравнений дает координаты

иско­

мой точки: # = уз; •*==із •

Или 1,1539 кг компонента I, 3,0769 кг

компонента II. Общая

стоимость смеси равна

1,1539X2,0 +

+ 3,0769X1,5 = 6,923 руб.

Она меньше стоимости

в двух

других

вершинах, поэтому точка d соответствует оптимальному реше­

нию задачи.

решении

задачи мы

ограничивались

При

графическом

двумя

компонентами.

Но для

составления

смеси требуе­

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

Успешная работа предприятия, в частности ритмичность его

191

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

Составление производственных планов — это важная об­ ласть применения линейного программирования, где оптимум определяется с учетом целого ряда факторов, основным из ко­ торых является обеспечение максимального чистого дохода при минимальных затратах.

Для уяснения сущности линейного программирования, ис­ пользуемого при планировании, разберем простейший пример. Положим, что заводу по плану необходимо выпустить продук­ ции Л = 250 т, продукции ß = 150 т, продукции С= 225 т. Про­ дукция изготовляется аппаратами I и II, каждый из которых, имея разную производительность, способен изготовлять все три вида продукции. Время, полезно затрачиваемое каждым аппа­ ратом, составляет 1250 мин. В табл. 22 даны затраты времени на изготовление единицы продукции каждым аппаратом.

 

 

Таблица 22

Затраты времени на изготовление единицы продукции,

мин

Наименование

 

Продукция

 

 

 

 

аппарата

А

В

С

 

I

2

5

5

и

3

4

10

Необходимо составить план загрузки аппаратов, выполне­ ние которого требует наименьших затрат времени.

Допустим, что аппарат I выполняет программу по продук­ ции А. Ее выполнение требует 250X2 = 500 мин рабочего вре­ мени. Оставшиеся 1250—500=750 мин можно использовать для

выпуска

продукции В (150X5 = 750) мин.

Таким образом, ап­

парат I

загружен полностью, а аппарат II

целиком производит

продукцию С. Однако при норме времени 10

мин

на аппара­

те II можно изготовить продукции С только

1250:

10=125 ед.

Следовательно, при данной

загрузке аппаратов (I вариант),

для изготовления оставшихся

225—125=100 ед. продукции С

не хватает мощности.

 

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

Из таблицы видно, что не хватает производственной мощ­ ности на изготовление только 225—(150—65) = 10 ед. продук­ ции С. Значит, при новом распределении и неизменных произ­ водственных мощностях аппаратов получен более эффективный

192

Таблица 23

Загрузка аппаратов по II варианту

Наименование

Продукция

 

 

 

Сумма

потребного

в

 

с

 

аппарата

 

 

времени, мин

 

 

 

 

 

I

250-2 = 500

150

- 5 =

750

 

1250

II

150-4 = 600

65

-10 =

650

1250

план производства продукции. Но и он не является окончатель­ ным. Оказывается при данных условиях наилучшим планом будет следующий (табл. 24):

Т а б л и ц а 24

Оптимальный вариант загрузки аппаратов

Наименование

 

 

 

Продукция

 

Сумма потребного

аппарата

 

Л

 

в

С

времени,

мин

 

 

 

 

 

I

62 -2

=

124

 

225 -5 = 1125

1249

 

и

188-3

=

564

150-4 = 600

1164

 

Как видно из табл. 24, для выполнения исходного производ­

ственного

задания

требуется 1248 + 1164 = 2413 мин. Таким

об­

разом, остается еще недоиспользуемая мощность, т. е. оборудо­ вание можно дополнительно загрузить в течение 2500—2413 = = 87 мин.

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

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

Цена всех установок і-й группы равна Ць где і=1, 2,

3, ....

т,

а фонд времени работы их, например, за год— 7’фі. Завод

мо­

жет изготовлять продукцию п наименований.

продук­

Общие затраты времени на производство единицы

ции ß-ro наименования (ß = l, 2, 3, ..., п) по всем технологиче­ ским операциям, выполняемым на і-й группе установок, обозна­ чим через Гр,. Требуется найти такие значения выпуска про­ дукции каждого наименования Pp (Рр!>0), которые позволили

бы получить

минимум недогрузки

всех групп установок с уче-

 

 

т

том их цены,

т. е. достигался бы

минимум функции V

 

 

і = I

13 З ак . 420

193

п\

 

P ß Tfi

Цг при условии

недопущения перегрузки какой-

ß=l

!

установЪк, т. е.

при

либо

группы

п

Т ф і - ^ Р ѵ - т ^ о .

Учет цены установок необходим для того, чтобы при отыска­ нии оптимального варианта добиться наиболее полного исполь­ зования дорогостоящих установок.

В конечном итоге отыскивается вариант, при котором до­ стигается минимальная недогрузка установок и минимум себе­ стоимости продукции каждого наименования в течение всего периода времени.

При неравномерном выпуске продукции в отдельные проме­

жутки времени (квартал, месяц) задача решается

раздельно

для каждого промежутка времени нахождением

минимума

следующей функции:

 

где

t — число всех отрезков времени в рассматриваемом

перио­

де;

Pß6 — выпуск

продукции ß-ro наименования в б-й

отрезок

времени (6 = 1, 2,

3, ..., ^); Т$і общие затраты времени на из­

готовление единицы продукции ß-ro наименования на і-й груп­ пе установок в б-й отрезок времени.

Искомые значения выпуска по каждому виду продукции в

•отдельные промежутки времени также должны быть неотрица­ тельными величинами:

п

P ß > 0 и п , - 2 ^ Р * > 0 .

Общий выпуск продукции ß-ro наименования за весь период

времени равен Pß = ^ Р[| > 0.

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

Т\і = rain

194

при

п

п

12

р ; > о , ^ P ß = P ß -

ß = l

ß = l

5 - 1

где Рр — годовой выпуск продукции ß-ro наименования по за­ данию.

7. Математические основы линейного программирования

Для удобства, точности, одинаковости формы обычно зада­ чи линейного программирования формулируются с помощью соответствующих математических обозначений. Математиче­ скую задачу можно записать двумя способами: 1) как систему линейных алгебраических уравнений, 2) с помощью векторов

иматриц.

Вобщем виде задача линейного программирования записы­

вается так:

1. Горизонтальный способ anXi + al2x2 + a13x3+ . . .+ ßinxn= bi

а2\Х\-\-а22х2-\-а23х3-\- ... -\-а2пхп— Ь2

йт\Х\ + ат2х2+ атЗ-^3“Ь • • * &mn%7i Ьті

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

п> т задача имеет бесконечное число решений. Оптимальное решение определяется из условия:

С1Х1+ С2Х2 -f- с3х 3+ ... + спх п -> оптимум.

Все выражения записаны как равенства. Если исходная задача содержит неравенства, то они преобразуются с по­ мощью дополнительных переменных в равенства. Если пере­ менные Х\, х2, ..., хп являются решением задачи, то среди них переменных > 0 не более т. Ниже дается пример составления уравнений.

Имеются следующие данные:

 

1 +1,3 • количество

кг

компонента

11 +

1.

5 ■количество

 

кг компонента

+1,35 • количество кг

компонента III > 4 ,0 -кг вещества

В.

 

 

 

2.

20 •количество

кг компонента

I + 30 • количество

кг

компонента

11 +

+ 25 • количество кг компонента III >

2 кг вещества А.

 

 

руб.-коли­

3.

Необходимо

минимизировать

стоимость смеси С: С = 1,1

чество

компонента

1

+ 2,2

руб. • количество компонента

II

+ 1,6

руб.-коли­

чество компонента III.

I

через

компонент I I — х 2, компонент III— х 3.

Обозначим компонент

Тогда запись уравнений будет такой:

 

 

 

 

 

1,5 Хі +

1,3 дг2 +

1,35 хз»4 кг вещества В ,

 

 

 

 

20 Хі

+

30 х 2 +

25 х 3 >

2 кг вещества А ,

 

 

 

 

1,1 X,

+

2,2 х 2 +

1,6 х ъ

=

min.

 

 

 

 

 

2. По вертикальному способу уравнения записываются следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

13*

 

 

 

 

 

 

 

 

 

 

 

196

■ 1,5]

Г

і,з ]

Г 1.351

>

•*1 20,0 + * 2

L

30,0

+ *з 25,0

>

1,1

2,2

1,6

=

I

д-с ■«t

 

m in _

3. Матричная система обозначений. Вертикальная запись уравнений мо­ жет быть представлена в виде следующей компактной системы обозначений:

Г 1,6

2,2

1,6

1

[min]

1,3

1,35 і

X ,

= 4

[20

30

 

.*3.

=2.

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

ким образом:

Сх = min,

А х = В, где С — строка стоимостей;

х — столбец

переменных;

В — столбец

ограничений; А — таблица,

состоящая

из заданных

исходных данных (исключая стоимости).

строк (в

нашем случае

А представляет собой

матрицу, состоящую из иг

из двух) и п столбцов (в нашем случае из трех).

 

 

8. Решение транспортной задачи

Задача состоит в том, чтобы удовлетворить требования отгрузки про­ дукции потребителям’ с минимальными транспортными расходами (например, отгрузка жидкого кислорода с завода нескольким потребителям). Величина транспортных расходов для каждого маршрута перевозки рассчитывается по транспортным тарифам в зависимости от расстояния транспортировки и .типа применяемых транспортных средств. В задаче предусматриваются два допу­ щения:

1.Затраты на транспортировку являются линейными, т. е. затраты на единовременную транспортировку 10 ед. продукции принимаются такими, как

ина десятикратную транспортировку единицы продукции.

2.Общее количество продукции равно общей потребности.

Для облегчения решения задачи исходные данные записываются в таб­ личной форме (табл. 25).

 

 

 

 

Т а б л и ц а 2 5

 

 

Стоимостная матрица

 

 

 

Р,=4

Ра= 5

Р,=6

Жг =

3

15

8

160

=

3

24

18

36

=

4

30

8

12

М \ — 5

34

8

6

Р п — объем

производства продукции

( Р \ — на первом

предприятии,

Р г — на втором,

Р3

— на третьем).

М п —-объем потребляемой

продукции

(М і — первым потребителем,

М г— вторым,

Л43— третьим, М 4 — четвертым).

В матрице приведены затраты на транспортировку единицы

продукции.

Например, доставка

единицы

продукции

с

предприятия Р і потребителю Mt

обходится в 16 руб.,

а с предприятия

Р 2

М 4— 8 руб. и т. д.

 

Решение этой задачи основано на методе последовательного улучшения

начального решения.

Из приведенных

в таблице данных ясно,

что один из

маршрутов неприемлем и не используется, так как стоимость доставки еди­ ницы продукции чрезмерно высока (160 руб.). Для дальнейшего решения за­ дачи необходимо установить начальное распределение. Оно может быть по­ лучено двумя способами: а) если оно сделано опытным аналитиком и работниками транспортной службы, то может быть достаточно близким к оптимуму; б) использованием правила «северо-западного угла».

196

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

Таблица 2 6

Матрица первоначального распределения

 

 

Р,=4

Pä=5

Р3= 6

лк =

з

3

 

 

м 2 =

з

1

2

1

Л13 =

4

 

3

М А=

5

 

 

5

В табл. 26 представлено распределение доставки продукции потребите­ лям с соблюдением указанных равенств. Возникает вопрос, будет ли это ре­

шение

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

Оно характеризуется общими затратами в

сумме

174 руб.

(16-3+24 • 1 + 18 • 2 + 8 • 3+12 • 1+

6 •. 5 = 174 руб.). Чтобы дать оконча­

тельный ответ на этот

вопрос, проверим

неиспользованные маршруты

пере­

возок. Например, маршрут с предприятия Р 2 к потребителю М х не был использован. Если по данному маршруту отгрузить одну единицу продукции, то столбец и строка матрицы будут сбалансированы только тогда, когда пе­

ревозки

по каждому из маршрутов Р ХМ Х и Р 2М 2

уменьшатся на единицу.

Поэтому

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

по маршруту Р \ М 2 надо

увеличить на одну единицу. Транспортировка единицы продукции по марш­ руту Р ^ М ! стоит 8 руб. Изъятие перевозок единицы продукции по маршру­

там Р \ М Х и Р 2М 2 уменьшит транспортные расходы на 16 и 18

руб. Увеличе­

ние перевозок на одну единицу

по маршруту Р \ М 2

потребует

дополнительно

24 руб. Матрица распределения

для данного случая представлена в табли­

це 27. Таким образом,

конечный

результат составит: 8—16— 18+ 24 = —2 руб.

 

 

 

 

Т аблица 27

Матрица распределения при введенном

 

 

маршруте Р

М А

 

 

І м

Р п

Р; 1

Р , - 5

Р, 6

 

 

 

 

 

 

 

 

Л4, =

3

2

1

 

 

М , =

3

2

1

 

 

М 3 =

4

 

3

1

 

М А=

5

 

 

5

 

Знак (—) указывает, что использование маршрута Р 2М Х уменьшило транс­ портные расходы на 1 руб. на каждую единицу продукции. Как показывает матрица, не используются маршруты Р 3М Х, Р 3М 2, Р ХМ 3, Р 2М А) Р ХМ А. Исполь­ зование каждого из этих маршрутов приводит к потерям, т. е. по каждому из них конечный результат будет иметь знак ( + ). Например, составляем ре­ шение по маршруту Р ХМ А. Оно требует определенных усилий для сохранения баланса по строкам и столбцам, поскольку включение данного маршрута изменяет объем перевозок по маршрутам Р ХМ А> Р 3М Аі Р 3М 3, Р 2М 3, Р 2М 2, Р ХМ 2. Распределение перевозок по маршрутам при использовании маршрута Р \ М А показано в табл. 28. Конечный результат при использовании маршрута Р \ М А будет 34—6+12—8+ 18—24= 26, т. е. использование данного маршрута дает

197

Т а б л и ц а 2 8

Матрица распределения при введении маршрута РіЛІ4

 

рп

Я,=4

яа= 5

я3= 6

 

 

Л4, =

3

1

2

+

М2 = 3

2

1

+

•Л43— 4

-+

2

2

М і =

5

1

-1-

4

потери в сумме 26 руб. Маршрутом с наибольшим выигрышем в нашем слу­ чае является Р г М ь По нему в рассматриваемой задаче можно отгрузить две единицы продукции, после чего мы придем к окончательному распределению маршрутов, которое представлено в табл. 29. Конечный результат, опреде-

 

 

 

 

 

Т а б л и ц а

2 9

Матрица конечного распределения

 

 

 

Р « = 4

Я а— 5

Я 3= 6

 

Af1 =

3

1

2

+

 

уИ2 =

3

3

-+

+

 

Л43 =

4

-f

3

 

1

 

М4 =

5

Ч-

+

 

5

 

ляющий транспортные

расходы,

составит:

 

1 • 16+2 • 8 + 3 • 24+3 • 8+1 • 12+

+ 5 - 6 = 1 7 0 руб., что

на (174—170) = 4 руб.

меньше, чем

по первоначаль­

ному распределению.

Знак ( + ),

поставленный

во всех неиспользованных

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

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

9.Применение математического анализа при решении экономических задач

Ряд экономических задач в области проектирования и экс­ плуатации оборудования, используемого в газоразделении, мо­ жет с успехом решаться обычными методами классического анализа (выбор оптимальной толщины изоляции, оценка срав­ нительной экономической эффективности различных типов теп­ лообменных аппаратов, выбор оптимального размера предприя­ тий и пр.). Приведем пример использования метода математиче­ ского анализа при решении экономической задачи (выбор опти­ мальной толщины изоляции агрегатов глубокого охлаждения).

Выбор основывается на установлении математической свя­ зи капитальных, эксплуатационных и расчетных затрат с

198

Соседние файлы в папке книги из ГПНТБ