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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 100

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Матрицы Llt f m задаются в форме ikf-матриц путем сравнения уравнения (VI,36) с другой его формой (VI,30):

dx\1)= 2 j 2 L i7 (nW!n>dW'n>

п І

Далее невозможно последовательно решать уравнения (VI,38)

без выбора значений для

2ДѴ т , поскольку каждая такая матрица,

за исключением L ™(1), зависит от переменных, связанных со всеми

установками в замкнутом

контуре.

Сначала необходимо выбрать значения составляющих вектора ш (cL^lш уже известен) и выходного вектора Затем можно выразить х(2) через х)1’ и Wa>, используя уравнения для установки 1, и, следовательно, выразить левую часть уравнения (ѴІ,38)(1) в виде функции Wa \ После этого можно решить уравнения и определить х(2>. Не нужно выбирать значения составляющих вектора cL[U{2\ так как он находится как линейная комбинация cLQ (1) и с£(112) (1)-:

сДі) <2)= cL'i) + СІ& Ч'Л/'і) ; (VI ,30)

Уравнение (ѴІ,38)(2)

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

относительно

W<2), по которому, в свою очередь, определяется х(3). Из уравнения

(VI,35) видно, что c Z /j Y сз>

м о ж н о получить из с Ѵ ^ (2) на

основе соот­

ношения

CL<P (3)=сі'1) (2)Д/(2)

(VI,40)

 

поскольку М (2) находится при решении уравнения для установки 2. Однако, чтобы решить уравнение (ѴІ,38)(3) относительно W VJ),

необходимо выразить

через

ТѴ(3) и уже известные величины.

Поскольку TV]31 в общем случае

является функцией

обоих входов

установки 3, это нельзя сделать без выбора значений

некоторых со­

ставляющих взаимосвязанных векторов. Например,

если векторы

х(3), хш и хт имеют

одну и

ту же размерность,

достаточно вы­

брать значения составляющих одного из двух векторов хш или х(6). Зададимся значением хш , тогда уравнение для установки 3 может

быть решено по

хш относительно х(3\

W(3) и выбранного значения

хш . Используя

полученное выражение

для

N 13\ можно выразить

левую часть уравнения (ѴІ,38)(3) через Wi3)

и известные (или вы­

бранные) величины. Оно может быть поэтому решено относительно W(3), что определяет хш.

Теперь мы подошли к уравнению (ѴІ,38)(5), которое можно ре­ шить без задания значений каких-либо других величии, поскольку c L ^ (в> находится из cL[^ <3) в соответствии с

c L [u M = cL [\> (a>M{f

(VI,41)

а матрица М^І известна из решения уравнения для установки 3. Левая часть уравнения (VI,38) может быть поэтому выражена как функция Ww \ следовательно, уравнение можно решить относи­ тельно х2и . Однако можно также рассчитать х2П по уже определен­ ным значениям xw и Wn\ причем, как правило, значения векто­ ров, вычисленных этими двумя различными способами, не совпа­ дают.

'312


Наконец, приступим к решению уравнения (ѴІ,38)(4). Нет необходимости выбирать значение с£(11И4), поскольку оно находится

из

согласно

 

 

 

 

 

 

cL[v

<«> = cLft (3ШЙ>

 

(VI ,42)

а матрица

М^І известна из

решения уравнения

для

установки 3.

Так как значение

ж(4) было

выбрано ранее, уравнение (ѴІ,38)(4)

может быть решено относительно \УШ, что определяет

значение

х(6>,

которое не совпадает с заданным входным вектором х.

мат­

Решив уравнения, вычисляют все матрицы

и,

значит,

рицу Я -1, так что,

исходя из значения вектора

(2), выраженного

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

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

подустановок, как показано на рис. ѴІ-7, то

число переменных

в итерации не возрастает. Просто матрицы ІѴ<2)

и N m распадаются

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

4. МНОГОУРОВНЕВЫЙ МЕТОД ОПТИМИЗАЦИИ ХТС

Метод многоуровневой оптимизации позволяет провести декомпо­ зицию большой задачи оптимизации на последовательность более мелких задач оптимизации. В основном метод осуществляется на двух уровнях. На первом уровне подсистемы (элементы) ХТС опти­ мизируют независимо друг от друга. Второй уровень служит для согласования первых уровней оптимизации с целью достижения об­ щего оптимума системы. Если оптимизация подсистемы ХТС сама выполняется посредством двухуровневого алгоритма, полный алго­ ритм оптимизации имеет многоуровневую иерархическую дерево­ видную структуру.

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

313


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

Сущность метода многоуровневой оптимизации поясним на при­ мере оптимизации ХТС, состоящей из двух элементов, которые охва­ чены рециклом (рис. ѴІ-8, а). Предполагают, что выходы элементов У1 , У2) zi и z 2 являются непрерывными функциями входов М х и М г.

а б

Рис. ѴІ-8. Пояснение сущности метода многоуровневой оптимиза­ ции (о) и декомпозиция задачи полной оптимизации (б) простой контурной ХТС.

Задача оптимизации заключается в определении М г и М 2, так чтобы максимизировать функцию прибыли обеих подсистем, т. е. максимизировать

{fl (Mi, 4) + h(M2, x2))

 

при следующих условиях:

 

а) x!=z2 (М2, Х2 )

 

б) xi = z 1 (M1, л*)

(VI,43)

в) Ми M 2£ S

 

где fi, /2 — прибыли, связанные с эксплуатацией обеих подсистем; S — мно­ жество возможных входов.

Приведенная задача оптимизации далее называется о с н о в н о й или з а д а ч е й п о л н о й о п т и м и з а ц и и . Принимают, что

максимум существует и находится в точках М, X, где М иХ удовлет­ воряют условиям (VI,43), т. е. максимум

{fi (Ми хі)+fz (М2, **)} = /і(АГі, ч )+h (Мъ,ч)

Декомпозицию задачи полной оптимизации на две более мелкие задачи осуществляют путем разрыва потоков, соединяющих две под­ системы, как показано на рис. ѴІ-8, б. Каждому из разрываемых потоков приписывают некоторые цены Р г и Р 2 (цены некоторых входного и выходного продуктов подсистемы). Первую подсистему поэтому оптимизируют посредством выбора М х и х ѵ Формально имеем:

314


П о д з а д а ч а 1. Максимизировать

{fi (Мъ хі) + РіЧ {Мх, *і) — iV l}

Mi, хх£ .S’l

П о д з а д а ч а II. Максимизировать

{/2 (ЛГ2, *г) + Т2г2 {M2, <л 2)—P xx2]

M%, %2(z $2

где S = S i U S 2.

Ограниченные множества S x и S 2 заключают оптимизирующие переменные в реальные физические пределы и тем самым гаранти­ руют, что максимумы подзадач существуют для любого множества цен Р.

Решение подзадач I и II дает решение основной задачи, если цены Р ± и Р 2 выбраны правильно. Выбор правильных цен является зада­ чей второго уровня алгоритма. Этот уровень выбирает цены так, чтобы свести разности (хг — z2) и (х 2 — zx) к нулю; тем самым удо­ влетворяется условие равенства одноименных технологических по­ токов ХТС.

Действительный алгоритм для подбора цен разрабатывают ис­ ходя из выражения функции Лагранжа для первоначальной задачи оптимизации:

L (М, X, Р) = fi (Mi, xi) + / 2 (AT* х2)-\-Рі [zi {Mi, хі)ш^х2]-\-

Р2 [ н { М 2, x2)x i]= {/i + PiZi— P2xi] + {/2 + T2z2—Pix2) (VI,44)

где L — функция Лагранжа основной задачи; Рх и Р 2 множители Лагранжа (или цены).

Функцию Лагранжа обычно используют для того, чтобы устано­ вить условия стационарности, которые справедливы в оптимуме основной задачи. В следующем разделе, однако, показано, что при относительно умеренных ограничениях необходимым и достаточным условием для основной задачи, которая должна быть оптимизиро­ вана, является максимизация функции Лагранжа посредством под­ бора М ъ х х и М 2, х 2.

Задача оптимизации функции Лагранжа разлагается непосредст­ венно на подзадачи I и II, поскольку первый член в скобках уравне­

ния (VI,45) представляет собой целевую функцию

подзадачи I,

а второй член — целевую

функцию подзадачи II, т. е.

max {fi + Р і Ч — iV i} + max {/2-f P 2z2 — Ріх2}ш*

 

M t X 1

М г Х г

 

=

max L (M, x) = y (P)

(VI,45)

 

Mx

 

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

315


Нише выводятся следующие свойства двойственной функции:

У (Р) Sä /і {Mi, хх) +

/2 {м2, х2)

(VI ,46)

nun у (P)—fi {Mi, xi) + f2{M2, x2)

(VI ,47)

dy

dL (M*, X * )

x l) - x l

 

~dPl

dPi

 

 

 

dy

dL (M*, X * ) —22

Qi ^2) *^1

(VI,48)

~dP^ =

dP2

 

 

Здесь M* и x* — значения M и х ,

которые максимизируют под­

задачи I и II.

 

 

 

Уравнение (VI,47)

справедливо только при ограничениях, изло­

женных в следующем разделе (стр. 318)- В нем говорится, что основ­ ной оптимум может быть достигнут посредством подбора цен Р так, чтобы максимизировать прибыль, достигаемую при оптимизации подзадач. Уравнение (VI,47) констатирует далее, что градиенты двойственной функции представляют собой просто разность между количеством товара, потребляемого одной подсистемой (например, х\ для первой подсистемы), и количеством товара (например, z*), которое другая подсистема решает поставить по существующей цене. Таким образом, градиенты двойственной функции имеются при условии выполнения небольших дополнительных вычислений, и задача подбора цены является просто задачей выпуклого програм­ мирования без ограничений. Оценку оптимальных цен можно по­ лучить, решая относительно Р г и Р 2 следующие линейные уравне­ ния:

/ 1 (M i, * 1 ) +

- P i

н ( M i , xi) — P 2 О

 

Q

Q

 

 

/ 2 ( M 2 , x 2) - \ - P 2

- z2 (M2> x 2)—P l = 0

(VI,49)

в которых все члены оценивают в возможных рабочих точках.

Теоретические основы метода многоуровневой оптимизации ХТС

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

В методе многоуровневой оптимизации используется свойство функции Лагранжа для многих реальных систем распадаться на ряд независимых «подзадач». Теоретически, однако, разложение данной функции не представляет значительных трудностей. Действительно важным является то, что максимум функции Лагранжа существует и что он идентичен максимуму исходной задачи. Поэтому рассмотрим

316