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

книги / Типовые задачи оперативного управления непрерывным производством

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
19.07 Mб
Скачать

дартов. Эти требования можно сформулировать в виде неравенств

Vi=vi(x, ÜBx)+Vi<bi, i= 1, m, (7-23)

где bi — предельно допустимые значения показателей. Значения потоков также ограничены плановыми за­

даниями Qk

 

 

 

yk=ÿk(x)+lk>Qk, k=\,

..., /.

(7-24)

Для случайных величин

£'г эти

неравенства мо­

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

P(yh^Qh, Vi^bu i= 1, .... т,

k = l

..., /)> Р 0.

(7-25)

Критерием является прибыль

 

F= 2

ckyk — d(x),

(7-26)

ft=i

 

 

где d(x) — затраты на управление производством, вклю­ чающие энергетические, стоимость сырья и т.' п., кото­ рые предполагаются линейными относительно перемен­ ных; си— стоимость единицы продукта.

Основное отличие полученной задачи с вероятност­ ным ограничением (7-25) от задачи предыдущего пара­ графа состоит в том, что показатели качества Vi стати­ стически связаны, а следовательно, статистически свя­ заны и величины |'i.

Обычно зависимости для части показателей качест­ ва Vi(x, i>nx) и некоторых (или всех) параметров пото­ ков уи(х) определяются достаточно точно и эти огра­ ничения можно исключить из (7-25) и рассматривать как детерминированные.

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

Перейдем от вероятностных ограничений к детерми­ нированным. Совокупность показателей и* является слу­ чайным вектором V с условной плотностью распределе­

ния /?(v|x) или условной функцией распределения <F(o|x), зависящей от управляющих воздействий х. Используя p(v\x) или &~(v\x), вероятностное ограни­ чение можно представить в виде

Ь \

^тп

 

 

 

(7-27)

j ... j p(v\x)dv>P„

 

—СО

—00

 

 

 

 

ИЛИ

(Oij • •

IX)^-Po»

 

 

 

 

 

где dv означает произведение dv±, ..., dvm.

 

Переходя от случайных величин

v{ к

нормированным

погрешностям

 

V, (х,

и„х)),

получаем:

,f (и)=

”( ...

f p(V)dV»P„

(7-28)

 

•/

J

 

 

 

 

—ОС

>00

 

 

 

где щ = — (bi v{ (x, JCbx)); p (¥)—плотность

распреде-

ei

 

 

 

 

 

ления вероятностей нормированного вектора погрешно­ стей 4е; ÎF(u) —функция распределения; сГ¥—произве­ дение (P¥i, . . dWm-

Полученная эквивалентная детерминированная зада­ ча является задачей нелинейного программирования с критерием (7-26) и ограничением (7-27). Выпуклую функцию £?~(и) можно линеаризовать, например, каса­ тельными плоскостями в точках ит={иг\, ..., urm}, г— =1, ...,/? вида

R

 

(7-29)

f (« 0 + S

Аг/ («/ —urt),

г =

\

 

где

А',=* du-1, IIи = и г'

Тогда неравенство (7-27) приближенно заменим си­ стемой ограничений

Ш

f V) + 2 h'i(“/ - “r<)? P.. r = l . . „ , R (7-30)

(7-31)

в которой выбирались точки аппроксимации, где Игмин — минимальные, Мгмакс —максимальные значе­ ния щ.

Этот метод применим при небольшом числе ограни­ чений, так как необходимо вычисление m-мерных и —1)-мерных интегралов, а линеаризация приводит

кзначительному увеличению размерности.

Вслучае двух ограничений и нормальной функции

распределения ограничения (7-30) имеют следующий вид:

nr\Ui -"I- hr2«2^ Ро iï 'о>

P — коэффициент корреляции между

£2;

Ф (г)—у==г- J ехр (

dx—функция Лапласа.

—00

 

 

Детерминированная задача (7-24), (7-26), (7-27) мо­ жет быть решена другими методами, например методом аппроксимирующего программирования. В этом итера­ тивном методе (7-27) заменяется линейным ограниче­ нием в окрестности некоторой точки, которая была получена при решении предыдущей задачи. линейного программирования. В этом случае не возрастает раз­ мерность задач линейного программирования, однако основные трудности, связанные с вычислением /«-мер­ ных интегралов, сохраняются.

Пример 7-2. Сформулируем задачу оптимального смешения ко­

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

времени свойств топлив

Основная задача смешения состоит в под­

боре такого соотношения

компонентов, при котором смесь будет

иметь

показатели

качества,

удовлетворяющие

соответствующим

1

Эта задача

разработана

совместно с А.

П. Кравченко щ

В. О. Чииакалом.

 

 

 

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

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

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

Зависимости для расчета всех показателей кроме температуры застывания имеют малые погрешности, и их будем считать детер­ минированными. Ограничения, полученные на основе этих зависи­ мостей, рассматривались в § 4-5. Это ограничения на вязкость,

плотность, содержание серы, отгон при

250°С. Ограничение на от­

гон при 250°С имеет следующий вид:

 

 

п

 

 

 

 

2

(ci

Ссм) %i

О*

 

/=*

 

 

 

где Cj — отгон при 250°С,

%; / — индекс компонента;

ссм — требуе­

мое значение параметра

смеси;

Xj — количество /-го

компонента

в смеси; п — число компонентов смеси.

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

Зависимости для температур застывания t>i, v2 и расчеты опти­ мального смешения приведем для смеси трех компонентов: мазута — хи дизельного топлива АТ—х2, дизельного топлива ВТ — х2. Зависи­ мости для температур застывания имеют следующий вид:

з

Ü1= в|0'+ 2

^

^СМЗО)

1

 

 

3

 

 

V2 — Яго "Ь 2 / ^ /)

^см>

I

 

 

где tj — температура застывания /-го компонента; *См, ?смзо — пре­ дельно допустимые значения температур застывания смеси на мо­ мент приготовления и через 30 дней; a n ( tj) — линейные или квад­ ратичные зависимости; £ь £2— случайные погрешности, распреде­ ленные нормально, с дисперсиями Oj2, о22 и коэффициентом корре­ ляции р.

Численные значения всех коэффициентов и параметры распре­ деления h, !2 приведены в табл. 7-1. Здесь приняты линейные зави­ симости aa(tj).

 

Таблица

7-1

и и tz SJ

Р *с м *с м зо Лц> an ÛJ3 Г‘ 13 °20 а*î Лаа

^аз

231—23 el 3,5 |з |о,5б1 - 5 | - 5

1- 6.51—6,36| 33,811- 2,75

1

СО

| —3,93 29,831- 7.2

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

Функцию распределения F(uit и2) аппроксимируем кусочно-ли­

нейной зависимостью в области щ,

и2^ 0. В точках

1—9 (u,rj Ü2r\

r = 1 . . . 9, показанных на рис. 7-1,

аппроксимируем

F(uu ц2) ’каса­

тельными плоскостями и заменим ограничение F(uu

и2) ^ р 0 линей­

ными ограничениями вида (7-32).

 

 

2 -иг

 

 

 

 

1,5

1

2

J

 

О

О

о

 

1

4

5

6

 

0

о

о

 

0,5

7

8

3

 

0

о

о

 

_

1

. 1

1

Of

___1

0

0,5

/

1,5

Z

Рис. 7-1. Точки аппроксимации.

Рис. 7-2. Зависимость ре­

шения от

вероятности

выполнения

ограниче­

ний Ро.

 

Значения

hf, hf, F0r для численного примера приведены

в табл. 7-2. Для упрощения задачи другие ограничения кроме тем­ ператур застывания не учитываются.

 

 

 

 

 

 

 

Таблица 7-2

Строка

 

Столбцы матрицы

 

Вид

Правая часть

1

2

3

4

5

 

ограни­

ограничений

 

 

Х2

х*

 

 

чения

 

i

 

U\

Un

 

 

1

—16,36

33,81

—2,75

—3,5

—3

sa

1.6

2

—3,93

29,83

—7,2

 

=

4,3

3

1

1

1

—0,328

—0,0443

а

1

4

 

 

 

 

Л,—0.446

5

 

 

 

—0,21

—0,0749

 

PQ—0,4887

6

 

 

 

—0,1025

—0,1025

 

Ро—0*583

7

 

 

 

—0,284

—0,114

5?

Ро—0,379

8

 

 

 

—0,17

—0,17

Ро—0,41

9

 

 

 

—0,0749

-0.21

2

Р0—0,4887

10

 

 

 

—0,214

—0.214

5?

PQ—0,128

11

 

 

 

—0,114

—0,284

2

/>0—0,379

12

 

 

 

-0,0443

—0,328

 

/>„—0,446

Приме чание . Дополнительно учитывались ограничения •*,, х„

Матрица ограничений задачи линейного программирования при­ ведена в табл. 7-2. Результаты решения даны в виде графика, построенного по значениям PQ в четырех точках (рис. 7-2). Для

сравнения эта же задача решалась путем аппроксимации F(u\, ц2) секущими плоскостями по точкам щт, и2г (рис. 7-1). Полученные значения переменных отличаются от графиков на рис. 7-2 не более чем на 0,004, т. е. точность аппроксимации F (и|, ц2) девятью ли­ нейными ограничениями достаточна. Как видно из графиков (рис. 7-2), линейная аппроксимация в данном примере дает хоро­ шие результаты.

Ввиду вычислительных трудностей, возникающих при решении задачи (7-28) даже при небольшом числе огра­ ничений (т^ З ), для получения приближенного реше­ ния практических задач приходится устанавливать ве­ роятности выполнения отдельных ограничений. Таким образом, осуществляется переход к задаче с m вероят­

ностными

ограничениями Pi{vi(xt

пВх )^ 5 г}^Р,0, i=

 

 

 

m

=1,

m,

где Pio выбираются так,

чтобы п РгО—Ро-

/=1 Эта задача с вероятностными ограничениями в свою оче­

редь сводится к задаче линейного программирования путем введения запасов по каждому ограничению, опре­ деляемых дисперсией о2* и требуемой вероятностью выполнения Р»о.

б) Случайные коэффициенты матрицы

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

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

Рассмотрим такую постановку на примере задачи оптимального смешения кормов [51].

Корма получаются путем смешения различных ком­ понентов (травяная и рыбная мука, зерно и т. п.) в ко­ личествах Xj, /= 1, ..., п. Требования к содержанию пи­ тательных веществ и витаминов в кормах определены. Стоимость смеси должна быть минимальна. Содержание питательных веществ в каждом компоненте непосредст­ венно перед смешиванием не измеряется, а контроли­ руется выборочными анализами. Таким образом, содер-

жание i-го питательного вещества в /-м компоненте можно рассматривать как случайную величину aij с нормальной функцией распределения, заданную мате­ матическим ожиданием йц и дисперсией:

Предполагается, что случайные величины atj неза­ висимы. Для рассматриваемой далее задачи существен­ но, что содержание i-ro вещества в /-м продукте и со­ держание его в k-ж продукте аш независимы.

Требование к содержанию i-ro вещества в смеси в количестве не менее Ь{ записывается как неравенство

П

2 auxi f h“

(7-33)

/=5 I

 

в котором ац —случайные величины.

неравенств

Можно потребовать выполнения этих

с некоторой вероятностью

 

 

(7-34)

и затем перейти к эквивалентной детерминированной задаче, используя нормальную функцию распределе­ ния ciij.

Для^независимых и нормально распределенных случай­

ных величин atj получим дисперсию о2,, (л:) и математиче­

п

ское ожидание 8* (л) невязок (JC) = ^ &ijX; bt:

п

n

s , M = 2 auxi - bi

и вместо (7-34) получим:

Ф(Ь(х)1о1(х)) ^Poi

или

(7-35)

8/(х) — 1о/) °/ (■*) ^

Здесь Ф(0 —интеграл Лапласа,

 

 

ф « = т к —СЮ

К

/2

Ф^ЧРм)—обратная функция.

^Подставив б,- (л*) и о2,-(л) в (7-35), получим неравен­ ство

( п

\ 1/2

п __

 

2

° V 2/

< 2 a4x i ~ bi'

•••» т, (7-36)

/=i

!

/=1

 

которое является выпуклым, так как Ф-1(Рог)>0. Критерий задачи — минимальная себестоимость

смеси

F = 2 C/Jty —min,

(7-37)

/=i

 

где с,-— стоимость /-го компонента.

Полученная задача (7-36), (7-37) является задачей выпуклого программирования. В [51] эта задача при небольшой размерности решается методом штрафных

функций.

когда ац,

aih

коррелированы, получим

В случае,

квадратичные

отклонения,

подобные

(7-36):

 

( п

 

п

\ 1/2

 

2 °V2/+ 2 КфХкхА <

 

/=1

/. А=1

/

 

 

i¥=k

 

 

< 2

/■*■/'

bl’

(7-38)

 

/=i

 

 

 

где К— коэффициенты ковариации.

Вероятности Ро1 задаются так, чтобы получить тре­ буемую вероятность выполнения всех ограничений Р0 =

т

=П ^о/* Задача с одним вероятностным ограничением t=i

P{2aijXj^bi, i=l, ..., т}^Ро

не требует выбора P0i, но является более сложной.

Другая близкая по смыслу постановка задачи со­ стоит в требовании выполнения неравенств в среднем

П

2 CLijXj> bt, i =

т

(7-39)

/*=»

 

 

с ограничением на дисперсию содержания каждого ве­ щества в смеси

£ \ | 2 anxi t = l, —» /и.

Вид критерия (7-37) не меняется.

Так как 'случайные коэффициенты а*,- независимы,

получим систему неравенств:

 

2 °“*/**/< Dot» * =;1 »— » /и,

(7-40)

/=1

 

определяющих выпуклую область.

В случае, когда Щу, щъ коррелированы, получим вме­

сто (7-40)

ограничения, аналогичные (7-38):

п

п

2

°v * t/ + S KnkXjXk^ D oi, i= 1...... т.

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

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

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

сти

bi —случайные величины. При заданной вероятно­

сти

выполнения каждого ограничения Рог получим за­

дачу с ограничениями, аналогичными (7-38).

Г л а в а в о с ь м а я

СТАНДАРТНОЕ ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ДЛЯ ЗАДАЧ УПРАВЛЕНИЯ И ПЛАНИРОВАНИЯ

8-1. Основные принципы построения программного обеспечения задач АСУ

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

Рис. 8-1. Функциональные задачи.

Связи функциональных задач показаны на рис. 8-1. Системы контроля и управления работают в реальном времени и являются локальными системами, которые могут быть выполнены на ЭВМ малой мощности, в то время как задачи управления зачастую требуют боль­ ших или средних мощностей. Сбор информации и учет осуществляется в основном через системы контроля или периферийные пульты (терминалы).

Для формирования базы данных разработано спе­ циальное программное обеспечение — банки данных, по­ зволяющие организовать хранение, защиту от записи (чтения), доступ для различных программ, информа-

Соседние файлы в папке книги