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

Согласование экономических интересов в корпоративных структурах - Гераськин М.И

..pdf
Скачиваний:
47
Добавлен:
24.05.2014
Размер:
2.24 Mб
Скачать

121

и коэффициенты a,b1 ,b2 вычисляются по формулам

 

b1

= b / , b2 =

b / , a = ξ13 (ξ11)b1 (ξ12 )b2 ,

 

 

1

2

где

 

 

 

 

= (lnξ11 lnξ12 )(lnξ12 lnξ23 )(lnξ11 lnξ13 )(lnξ13 lnξ22 ),

b

= (lnξ

32 lnξ13 )(lnξ12 lnξ23 )(lnξ33 lnξ13 )(lnξ12 lnξ22 ),

1

 

 

b2 = (lnξ11 lnξ12 )(lnξ33 lnξ13 )(lnξ11 lnξ13 )(lnξ23 lnξ13 ).

Рис. 3.1 – Формирование гипербол, аппроксимирующих множество Парето В общем случае К критериев уравнение гиперболической поверхности,

проходящей через К точек аппроксимации Аk (ξ1k ,ξ2k ,...,ξkK ),k K , имеет вид

ξК = a(ξ1)b1 (ξ2 )b2 ...(ξК 1)bК 1

(3.14)

с коэффициентами a,b1 ,b2 ,...,bК 1, получаемыми в результате

решения

системы уравнений

 

ξkK = a(ξ1k )b1 (ξ2k )b2 ...(ξkK 1)bК 1 ,k = 1,2,...,K .

(3.15)

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

122

первого ортанта или, в двумерном случае, первого координатного угла.

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

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

определяются коэффициенты аппроксимирующей поверхности путем решения системы (3.13), после чего вычисляются координаты центра

аппроксимирующей поверхности по формуле

1

ξ1C = ξC2 = ... = ξCK = ξC = (a) b1+b2 +...+bK 1+1 .

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

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

Алгоритм решения двухкритериальной задачи. В случае К=2

аппроксимирующая зависимость является гиперболой. С учетом нормализации вершина гиперболы принадлежит началу координат, так как асимптотами являются координатные оси. Линия ϑ( П ) π -оптимальных сочетаний критериев имеет вид, показанный на рис. 3.1.

Предлагается следующий алгоритм формирования минимаксно- оптимального управления:

123

н~

1)выбирается начальный закон управления u U , которому

соответствует сочетание критериев в точке Ан на рис. 3.1 (индексом "н" обозначено начальное значение);

2) определяется опорное управление, тождественное начальному ui = uн

при i = 0 (i - номер итерации) или полученному на предыдущей итерации

ui = uiC1 при i > 0 ;

3) определяется критерий с наибольшим нормализованным значением

ξi

= max ξi

[ u ]

 

 

k'

 

k K

k

i

 

 

 

 

 

 

 

 

 

 

и фиксируется значение другого

критерия

ˆ

= Rk [ui ],k ¹ k' ; область

~

Rk

U

дополняется ограничением

 

 

 

 

 

 

 

~

 

 

 

ˆ

,k ¹ k'};

 

U' = {u ÎU ,Rk [u ] = Rk

 

4)формируется управление uik' , удовлетворяющее условию

минимальности

ξi

[ ui

] = min ξi

[ u ],k' K,

k'

k'

~

k'

 

 

 

u U'

 

 

и вычисляются координаты точки Ai' ( ξ'1i ,ξ'2i ), принадлежащей множеству

ϑ( П );

5) определяется критерий с наибольшим нормализованным значением

ξi = max ξi [ ui ] ;

k'' k k' k K

6) задается приращение Δξi и вычисляется значение критерия с номером k ¹ k'' , соответствующее этому приращению, по формуле, обратной формуле

(3.7):

Rk

= (ξk [ uk' ] + Dξ

 

)(Rk

- Rk

)+ Rk

;

ˆ

i i

i

max

min

min

~

область U дополняется ограничением

~ = { Î = ˆ ¹ }

U'' u U ,Rk [ u ] Rk ,k k'' ;

124

7)формируется управление uik'' , удовлетворяющее условию

минимальности

xi

[ ui

] = min xi

[ u ],k'' ÎK,

k''

k''

~

k''

 

 

 

u U''

 

 

ивычисляются координаты точки Ai'' (ξ''1i ,ξ''2i );

8)вычисляются координаты центра гиперболы

1

ξC1i = ξC2i = ... = ξCKi = ξCi = (a)b+c+...+ z+1 ;

9) формируется управление uiC , соответствующее точке Ci (x1Ci ,xC2i ) или

ближайшей к ней точке,

если Ci Ф ;

для этого по координатам точки Ci

определяются значения исходных критериев

 

 

 

 

 

RC = ξC

( Rmax

Rmin ) + Rmin

,k K

k

ki

k

 

 

 

 

k

k

 

 

 

 

 

и находится управление u из условия принадлежности области

~C

 

~

 

 

k

 

 

C

ÎK},

U

= {u ÎU , R

 

[ u ] = Rk ,k

если Ci Ф , или, если Ci Ф , из условия

 

 

 

 

 

 

 

 

 

minmax

 

Rk

[ u ] - RC

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

k

 

 

 

 

 

 

 

 

u U k K

 

 

 

 

 

 

 

 

 

 

xiC1 - xiC

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10) проверяется условие

окончания

итераций

 

 

£ e. Если оно

 

 

выполнено, то точка Ci

считается приближенным минимаксно-оптимальным

сочетанием

критериев,

а ее прообраз uC

- минимаксно-оптимальным

 

 

i

 

управлением

u0 ; в

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

приращение уменьшается

Δξi+1 = 0,5Δξi и вычисления повторяются, начиная с шага 2.

Алгоритм решения многокритериальной задачи в общем случае. В

общем случае K критериев алгоритм имеет вид:

1)

выбирается начальный закон управления u

н

~

 

ÎU ;

2)

определяется опорное управление по правилу

 

 

ì

u

при i = 0,

 

 

 

u = í

н

при i >0;

 

 

 

i

uC

 

 

 

 

î

i 1

 

 

 

125

3) формируется

 

K

 

управляющих

зависимостей

 

 

uk ,k K путем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

последовательного решения K задач минимизации

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

 

 

 

ξki

[ uki ] = min maxξki [ u ],k = 1,...K ,

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

~

k K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

ì

 

 

 

 

 

 

 

 

u U ki

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ü

 

í ~

 

~

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

k

 

 

 

k

для всех k ¹ arg maxξ

i

[ u

i

 

],k Î K

ï

 

 

 

ïU

ki

= íu ÎU ,R

0

[ u ] = R

0

k

k

1

ý,

 

ï

ï

 

 

 

 

 

 

 

 

 

k K

 

 

 

 

 

ï

 

 

 

î

 

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

þ

 

в каждой из которых

 

k

=

 

i

( Rmax

Rmin ) + Rmin

,

 

i

 

=

 

i

 

+ Δξ

 

, где

 

i

-

 

ξ

ξ

 

ξ

 

 

ξ

R

 

 

i

 

 

 

 

 

 

 

 

 

 

k

k

 

k

k

 

 

k

 

 

 

 

k 1

 

 

 

k 1

 

значение k-го критерия, полученное в результате предыдущей задачи. В каждой из K задач начальным приближением служит управление ui при k=1,

uki 1

при k=2,3,...,K;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4)

 

вычисляются

координаты

точек

 

Ak

( ξk

,ξk

,...,ξk

),k K,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

1i

2i

Ki

принадлежащих множеству ϑ( П );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5)

вычисляются

координаты

 

 

центра гиперболической

поверхности

C = ( ξC

,ξC

,...,ξC )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

1i

2i

Ki

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξC

= ξC

=...= ξC

= ξC

= (a)

 

 

 

 

1

 

 

 

,

 

 

 

 

 

 

b +b

2

+...+b

К 1

+1

 

 

 

 

 

 

1i

2i

 

Ki

 

 

 

 

i

1

 

 

 

 

 

 

 

 

где коэффициенты b1 ,b2 ,...,bК 1 определяются из решения системы

 

 

 

 

ξkK = a(ξ1k )b1 (ξ2k )b2 ...(ξkK 1)bК 1 ,k = 1,2,...,K ;

 

 

6)

формируется

управление

 

 

uC ,

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

точке

C или

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

i

ближайшей к ней точке,

если

Ci

Ф ;

для этого по координатам точки Ci

определяются значения исходных критериев

 

 

 

 

 

 

 

 

 

 

 

 

 

RC = ξC ( Rmax Rmin ) + Rmin

,k K

 

 

 

 

 

 

 

 

 

 

k

ki

k

 

 

 

 

k

 

k

 

 

 

 

 

 

 

 

 

и находится управление u из условия принадлежности области

 

 

 

 

 

 

~C

 

~

 

 

 

k

 

 

C

 

K},

 

 

 

 

 

 

 

 

 

 

U

= {u U , R

 

[ u ] = Rk ,k

 

 

 

 

 

 

если Ci Ф , или, если Ci Ф , из условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

minmax

 

Rk [ u ] RC

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

k K

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

126

 

 

 

 

 

7) проверяется условие окончания итераций

 

ξiC1 − ξiC

 

≤ ε . Если оно

 

 

выполнено, то точка Ci

считается приближенным минимаксно-оптимальным

сочетанием

критериев,

а ее прообраз uC

-

минимаксно-оптимальным

 

 

i

 

 

 

 

 

управлением

u0 ; в

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

приращение

уменьшается

Δξi+1 = 0,5Δξi и вычисления повторяются, начиная с шага 2.

Условия сходимости алгоритма. Алгоритм позволяет определить минимаксно-оптимальное сочетание критериев ξ0 за конечное число итераций. Докажем необходимое условие сходимости алгоритма в виде следующей теоремы.

Теорема 3.1. Для заданной точности решения ε > 0 всегда найдется такой номер итерации i, что различие между минимаксно-оптимальным сочетанием критериев и решением, полученным с помощью алгоритма, не превысит этой точности, то есть

ξ0 − ξiC

≤ ε .

(3.16)

Доказательство. Для случая К=2 гиперболы, соответствующие смежным итерациям, определяются уравнениями (рис. 3.1)

Γi 1: ξ2 = fi 1(ξ1 ), Γi : ξ2 = fi (ξ1 ).

Пусть приращение

Δξi

 

подбирается таким образом что на каждой

итерации точки

А'

и А''

лежат на гиперболе по разные стороны от точки

 

 

 

i

i

 

 

 

 

 

 

 

 

 

Сi . В этом случае можно подобрать такое δ [0,1], что

 

 

 

 

ξСk = δξ'k + (1 − δ)ξ''k , 0 ≤ δ ≤ 1, k = 1,2 .

Поскольку гипербола является выпуклой кривой, то из условия

выпуклости

(

) k )

 

(

k )

(

)

 

(

k )

[ ]

 

(

k

≤ δf

 

f

 

δξ' +

1

− δ ξ''

 

ξ'

+ 1 − δ

 

f

 

ξ''

для любого δ 0,1

вытекает, что

 

 

 

 

 

 

 

 

 

 

 

ξ2iC = fi (ξ1Ci )= fi (δξ'1 +(1− δ)ξ''1)≤δf (ξ'1)+(1− δ)f (ξ''1)

=δξ2iC 1 +(1− δ)f (ξ''1). (3.17)

127

Так как по построению точки Аi'' верно неравенство

f (ξ''1)= fi (ξ1iC1)≤ ξ2iC1 ,

(3.18)

то при подстановке ξC2i1 вместо f (ξ''1) неравенство (3.17) не изменит знака:

ξ2iC = δξ2iC1 + (1− δ)ξ2iC1 = ξ2iC1.

(3.19)

С учетом того, что по свойству симметричности

 

ξ1C = ξ2C = ξC ,

(3.20)

из (3.19) следует невозрастание последовательности точек {ξiC}:

 

ξCi ≤ ξCi1.

С другой стороны, по свойству симметричности (3.20) и

свойству

минимальности ξi

[ ui

] = minξi

[ u ],k K сочетание критериев

ξО

при

k

k

~ k

 

 

 

 

 

u U

 

 

 

минимаксно-оптимальном управлении uО ограничивает последовательность

{

i

}

снизу ξО

~

k

, k = arg max ξ

k

. Таким

образом, существует

точек

ξC

 

= minξ

 

 

 

 

 

 

 

u U

 

k K

 

 

 

предел

 

lim ξiC = ξО ,

а это

означает, что,

 

начиная

с некоторого номера i,

 

 

i→∞

 

 

 

 

 

 

выполнится условие (3.16). Теорема доказана.

Особенности применения метода аппроксимации. Решение многокритериальной задачи на основе метода аппроксимации множества Парето сводится к последовательности скалярных оптимизационных задач и предусматривает: а) формирование К Парето-оптимальных управлений; б) построение в соответствии со значениями критериев при этих управлениях гиперболических поверхностей (кривые Гi ,Гi 1 на рис. 3.1),

аппроксимирующих поверхность Парето в пределах малой окрестности опорного управления; в) нахождение точки сочетания критериев,

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

128

Предложенный метод позволяет определять минимаксно-оптимальное сочетание критериев ξ0 как в случае выпуклого к началу координат множества ϑ( П ), так и в невыпуклом случае, поскольку на предпоследнем шаге в невыпуклом случае ищется точка, ближайшая к C в смысле

u = argminmax

 

ξk [ u ] − ξC

.

~

k K

 

 

i

i

 

u U

 

 

 

 

 

 

 

 

 

Метод, кроме того, позволяет

 

учесть

приоритеты критериев,

 

 

 

 

 

K

 

задаваемые коэффициентами важности βk : åβk

= 1, βk > 0, k K . В этом

k =1

случае алгоритм применяется в неизменном виде, но нормализованные критерии подвергаются преобразованию: ξk = βk ξk , k K .

Общие результаты предложенного метода заключаются в следующем:

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

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

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

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

129

3.3. Метод многокритериальной оптимизации механизма корпоративного управления на основе теории графов

Рассмотрим метод, позволяющий осуществить многокритериальный

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

Граф Парето-оптимальных управлений. Управление, оптимальное по критерию (3.10), может быть выбрано путем сопоставления Парето-

оптимальных управлений u*k , сформированных по критерию (3.8).

Введем в рассмотрение параметр

hnm =

R

[u*

]R

[u*

]

, n,m K ,

(3.21)

k

m

k

n

 

 

 

 

 

 

k

 

 

Rk*

 

 

 

 

 

 

 

 

 

 

 

отражающий долю прироста (потерь) k-го критерия относительно его

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

u*n к управлению u*m (рис. 3.2). В случае hknm > 0 управление u*m является более предпочтительным по критерию Rk по сравнению с управлением u*n , в

противном случае более предпочтительным является управление u*n .

Сформируем граф управлений [101], вершинам которого поставим в соответствие Парето-оптимальные управления u*k ,k K , а дугам процессы

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

Определим веса дуг графа как алгебраическую сумму относительных приростов (потерь) критериев системы при переходе от управления u*n к

управлению u*m :

130

K

 

S nm = åhknm , n,m K .

(3.23)

k =1

Вес Snm представляет собой векторную характеристику дуги (перехода) от управления u*n к управлению u*m : при Snm>0 управление u*m является более предпочтительным по векторному критерию (3.3), чем управление u*n .

Выражение

для

параметров

S nm через нормализованные

значения

критериев получим, подставив (3.22) в (3.23):

 

K

K

*

*

K

 

S nm = åhknm = å

Rk [um ]Rk [un ]

= å(

 

k [u*m ]

 

k [u*n ]), n,m K .

(3.24)

R

R

 

k =1

k=1

 

Rk*

k=1

 

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

в которой первая и последняя вершины совпадают, например, u1* ,u*2 ,...,u*K ,u*1 .

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

S 12 + S 23 + ...+ S( K 1 )K + S K 1 = 0 .

В самом деле,

S 12 + S 23 + ...+ S( K 1 )K + S K 1 = K (R [u* ]R [u* ])+ K (R [u* ]R [u* ])+ ...+

å i 2 i 1 å i 3 i 2

 

 

 

 

 

 

 

 

i=1

i=1

+ åK (

Ri [u*K ]

Ri [u*K 1 ])+ åK (

Ri [u1* ]

Ri [u*K ])= −åK

Ri [u1* ]+ åK (

Ri [u*2 ]

Ri [u*2 ])+

i=1

 

 

i=1

 

 

i=1

i=1

+ ...+ åK (

Ri [u*K ]

Ri [u*K ])+ åK

Ri [u1* ]= 0.

 

 

 

 

i=1

 

 

i=1

 

 

 

 

 

 

Следовательно, сумма приростов (потерь) критериев при переходе от

управления u*n

к управлению u*m

равна взятой с противоположным знаком

сумме приростов (потерь) критериев при последовательном сравнении всех

Парето-оптимальных управлений, кроме u*n и u*m :

S nm = −[S12 + S 23 + ...+ S( n1 )n + S m( m+1 ) + ...+ S( K 1 )K + S K 1 ] .

(3.25)

Вершины графа управлений u*k ,k K характеризуются

значениями

параметров

Соседние файлы в предмете Экономика