СвёртываниеКритериев

advertisement
Методы замены векторного критерия
скалярным критерием
Одним из подходов к поиску компромиссного решения ЗВО является
сведение её к задаче параметрической оптимизации, т.е. сведение её к
однокритериальной (скалярной) оптимизации. Иначе говоря, частные
критерии Fi(X), тем или иным способом объединяются в составной
(обобщенный) критерий f(X)=Ф[F1(X), F2(X), . . . , Fm(X)], который затем
оптимизируется.
Метод взвешенных сумм (функций)
(Метод линейной свертки)
Существует несколько классических методов, относящихся к
многокритериальной оптимизации. Один из них – метод взвешенных
функций (method of objective weighting).
Идея этого метода заключается в том, что обобщённый критерий
записывается в следующем виде:
m
f ( X )    i Fi ( X ),
(1)
i 1
который называют аддитивным критерием. Здесь i0 являются весовыми
коэффициентами, которые задают предпочтение i - го критерия по
сравнению с другими критериями. Величина i определяет важность i - го
частного критерия. При этом более важному критерию приписывается
больший вес, а общая важность всех критериев принимается равной 1, т.е.
m

i
 1.
То, что решение можно получить, используя аддитивность
i 1
векторного критерия, высказал Парето. Он также ввёл понятие весовых
коэффициентов. Таким образом, мы получили однокритериальную задачу
математического программирования
min f ( X )  min
m
 F (X ) .
i 1
XD
i
i
XD
Замечание. Как правило, частные критерии имеют различную размерность.
Поэтому при образовании обобщённого критерия нужно работать не с
натуральными критериями, а с их нормированными значениями.
Нормированный критерий представляет собой отношение “натурального”
частного критерия к некоторой нормирующей величине. При этом выбор
нормирующего делителя должен быть обоснован. Возможно несколько
подходов к выбору нормирующего делителя:
 в качестве нормирующего делителя берут директивные значения
параметров, заданные заказчиком, т.е. предполагают, что в ТЗ на
проектируемый объект заданы оптимальные значения параметров:



в качестве нормирующих делителей берут максимальные
(минимальные) значения критериев, достигаемых в области
существования проектных решений (область D);
F (X )
Fi ( X )
, fi ( X ) 
;
fi ( X )  i
Fimax
Fimin
берут лучшие мировые достижения в данной области;
в качестве нормирующего берут разность между max и min значениями
критерия в области D
fi ( X ) 
Fi max  Fi ( X )
Fi ( X )  Fi min
или
.
f
(
X
)

i
Fi max  Fi min
Fi max  Fi min
Нормированные критерии будем обозначать через f i, т.е. аддитивный
m
критерий примет вид f ( X )    i f i ( X ). (2)
i 1
Какой определён принцип оптимальности?
Поскольку в области компромисса увеличение (уменьшение) одного
критерия может достигаться лишь ценой уменьшения (увеличения) другого
(или других) критериев, то справедливым является тот компромисс, при
котором абсолютный уровень снижения одного не превосходит суммарного
уровня увеличения других критериев.
Пусть имеется два решения X1 и X2. Тогда в соответствии с
изложенным принципом следует вычислить сумму абсолютных изменений
всех частных критериев, обусловленных этим переходом (переход от X 1 к
X2)
m
m
m
i 1
i 1
i 1
f   i ( f i ( X 2 )  . f i ( X 1 ))   i f i ( X 2 )   i f i ( X 1 )..
В случае f<0 решение X2 признаётся лучшим, чем X1, если f>0, то лучше
X1. Тогда оптимальным решением будет такое, для которого f0 при
переходе от него к любому другому решению, т.е.
m
m
i 1
i 1
  i f i ( X )    i f i ( X opt ), где Xopt - точка min, X любая точка из D.
Таким образом, принцип справедливой абсолютной уступки
(компенсации) приводит к утверждению, что оптимальное решение означает
минимизацию суммы нормированных частных критериев.
Иногда условия работоспособности позволяют выделить две группы
выходных параметров. В первую группу входят выходные параметры,
значения которых в процессе оптимизации нужно увеличить
Fi  ( X ) (производительность, вероятность безотказной работы), во вторую –
выходные параметры, значения которых нужно уменьшить Fi  ( X ) (расход
топлива, длительность переходного процесса). Тогда аддитивный критерий
m1
(2) примет вид f ( X )   i f i ( X ) i 1
максимизируется.

m2
 f
i 1
i

i
( X ) , где m1+m2=m. Критерий f(X) –
Замечание. Если решается задача выпуклого программирования, то
полученное решение (с использованием аддитивного критерия) является
оптимальным по Парето, т.е. оптимальное решение, полученное с
использованием метода линейной свёртки, лежит в области эффективных
решений. Доказать данное утверждение самостоятельно.
Рассмотрим пример. Переносной автомат для забивания стальных
дюбелей в бетонные стены состоит из корпуса с магазином, содержащим
запас дюбелей; подающего-спускового механизма с зарядами и ствола.
Требуется определить основные конструктивные параметры автомата –
длину ствола L и число дюбелей –N при следующих исходных данных:
число дюбелей , N12;
масса одного дюбеля с зарядом равна m=50 г;
масса ствола 1.6 кг/м;
масса корпуса 2 кг.
При фиксированной величине заряда и заданной массе дюбеля
скорость V выбрасывания связана с длиной ствола L соотношением V=k L ;
где k= 150
м 0.5
. Минимально допустимая скорость Vmin=100 м/cек. Масса
c
автомата не должна превышать 6 кг. Частными критериями являются
скорость выбрасывания и число дюбелей, помещающихся в магазине. Выбор
этих критериев объясняется тем, что чем выше V, тем надёжнее дюбеля
проникают в бетонные стены любой марки, а чем больше N, тем удобнее
работать. По мнению экспертов оба критерия V и N имеют одинаковую
важность.
Введём обозначения: F1(N,L)= k L - первый критерий (скорость);
F2(N,L)=N – второй критерий (число дюбелей).
Сформулируем задачу многокритериальной оптимизации:
maxF=max(F1, F2)
при следующих ограничениях
N12
V100
1.6L+ 0.05N+26.
Для определения оптимальных значений параметров будем использовать
m
аддитивный критерий f ( X )    i f i ( X ). Так как критерии имеют одинаковую
i 1
важность, то весовые коэффициенты можно взять равными единице. В
качестве нормирующих делителей возьмём максимальные значения
критериев, достигаемых в области существования проектных решений
(область D). Для определения нормирующих делителей будем использовать
уравнение баланса (массу автомата). Найдем Nmax из условия, что Vmin=100
м/cек. Определяем длину ствола, соответствующую минимальной скорости
2
V
L=  min  и подставляем в уравнение баланса. Получим, что Nmax=65.
 k 
Аналогично находим, т.е. берём N=12 и подставляем в уравнение баланса.
Получим Vmax=218 м/cек.
Таким образом, мы получили следующую задачу оптимизации:
найти максимум функции f ( L, N ) 
k L
N

Vmax
N max
при ограничении 1.6L+ 0.05N+26.
Так как наша функция f(L,N) монотонно возрастающая, то максимум
достигается на границе. Поэтому ограничение неравенство мы можем
заменить на ограничение равенство. Окончательно имеем
найти максимум функции f ( L, N ) 
k L N

218 65
при ограничении 1.6L+ 0.05N+2=6.
Для определения максимального значения f(L,N) с учётом ограничения
равенства мы будем использовать метод неопределённых множителей
Лагранжа. Получим следующую задачу безусловной оптимизации:
найти максимум функции ( L, N ,  ) 
k L
N

  (1.6 L  0.05 N  4).
Vmax
N max
Находим частные производные по L, N,  и приравниваем их к нулю.
Получим систему уравнений:

k
1


 1.6  0,
L Vmax 2 L

1

 0.05  0,
N N max

 1.6 L  0.05 N  4  0.

Решая
Vopt
эту
систему,
k2

,
2Vmax  1.6
N opt 
получим
следующие


k


 2Vmax  1.6 
2
значения:

1
,
0.05 N max
или =-4/13, Nopt=64, Lopt=0.49 м,
Vopt=105 м/cек.
Аддитивный критерий имеет ряд недостатков:
 Он выступает как формальный математический приём, придающий
задаче удобный для решения вид;
 В аддитивном критерии может происходить взаимная компенсация
частных критериев. Это значит, что значительное уменьшение одного
из них вплоть до нуля может быть покрыто возрастанием другого
критерия. Для ослабления этого недостатка следует вводить
ограничения на минимальные значения частных критериев и их
весовых коэффициентов.
Замечание. Хотя аддитивный критерий подвергается сильной критике, но
существуют задачи, где критерий качества должен удовлетворять
аддитивности. Например, в динамическом программировании эффект от
управления процессом F складывается из элементарных эффектов f k,
полученных на отдельных шагах процесса: F=fk.
Мультипликативный критерий преобразуется в аддитивный путём
логарифмирования целевой функции, то получим эквивалентный
аддитивный критерий, который обращается в максимум одновременно с
мультипликативным критерием. Таким образом, несмотря
на слабые
стороны, обобщённый аддитивный критерий позволят в ряде случаев
успешно решать многокритериальные задачи и получать полезные
результаты.
Мультипликативный критерий
Аддитивный критерий основан на использования принципа
справедливой компенсации значений нормированных частных критериев. Но
в ряде задач проектирования более целесообразным является оперирование
не с абсолютными, а с относительными изменениями значений частных
критериев.
Принцип справедливой относительной компенсации формулируется
следующим образом: справедливым следует считать такой компромисс,
когда суммарный уровень относительного снижения значений одного или
нескольких критерий не превышает суммарного уровня относительного
увеличения других критериев.
В математической формулировке условие оптимальности на основе
принципа справедливой относительной компенсации имеет вид
m

i 1
Fi ( X )
 0,
Fi ( X )
(3)
где Fi(X) – приращение величины i – го критерия, Fi(X) – i – го критерия.
Полагая Fi  Fi (X ) , можно представить (3) как дифференциал
натурального логарифма
m
Fi ( X )

d
(ln
F
(
X
))

d
ln
Fi ( X )  0,



i
i 1 Fi ( X )
i 1
m
(4)
Из выражения (4) следует, что принцип относительной компенсации
приводит к мультипликативному обобщённому критерию оптимальности
m
F ( X )   Fi ( X ).
(5)
i 1
Мультипликативный
критерий
образуется
путём
простого
перемножения частных критериев в том случае, когда они имеют
одинаковую важность. В случае неравноценности частных критериев
вводятся весовые коэффициенты i и мультипликативный критерий примет
вид
m
F ( X )   Fi i ( X ).
i 1
(6)
Мультипликативный критерий иногда представляется в виде отношения
произведений частных критериев (выходных параметров)
m1
F(X ) 
F

F

j
i
i 1
m2
(X )
, m1+m2=m;
(7)
(X )
j 1
где в числителе перемножаются все выходные параметры, требующие
максимизации и имеющие ограничения Fi  ( X )  TTi , а в знаменателе – все
выходные параметры, требующие минимизации и имеющие ограничения
Fi  ( X )  TTi . , где TTi – значение технического требования, предъявленного к
i– му критерию. Целевая функция (7) в дальнейшем подвергается
максимизации.
Достоинством мультипликативного критерия является то, что при его
использовании не требуется нормирование частных критериев. Недостатки
критерия: критерий компенсирует недостаточную величину одного частного
критерия избыточной величиной другого и имеет тенденцию сглаживать
уровни частных критериев за счёт неравнозначных первоначальных
значений частных критериев.
Пример 1. Производственная функция, отражающая овеществлённый
технический прогресс (модель Р. Солоу)
Yt  ALt K t1 , где Yt – выпуск продукции; L – численность рабочих; K – объём
основных производственных фондов. Здесь величины  и 1- следует
рассматривать как весовые коэффициенты.
Пример 2. Обнаружение сигналов в "белом" шуме. На вход RC – фильтра с
импульсной характеристикой h( )  e  поступает аддитивная смесь:
прямоугольный импульс s(t) плюс "белый" шум n(t). Требуется найти такое
значение , чтобы отношение сигнал/шум было максимальным, т.е. мы
желаем, чтобы значение сигнала было максимальным, а уровень шума –
минимальным ( – полоса пропускания RC – фильтра).
1.531
2
s(t)
ng
A
 1.712
T
Рис.3. Сигнал s(t)
t
0
2
0
0
100
200
g
Рис 4.Шум n(t)
300
299
2.531
4
62.152
100
2
xpg
zm
50
0
 1.712
2
0.987
0
100
0
200
g
300
299
0
0
0
Рис. 5.Сигнал плюс шум
100
200
m
250
Рис.6. Отфильтрованный сигнал
n(t)
sˆ(t)  nˆ(t)
s(t)
h()
Рис. 7.Система обнаружения сигнала
F1()=A(1-e-T)  max (уровень сигнала на выходе фильтра),
F2()=
N
 min (уровень шума на выходе фильтра).
2
где A,N,T - константы; F1 и F2 - имеют одинаковую размерность. Найти
оптимальную полосу пропускания , если справедлив принцип
относительной компенсации частных критериев. Согласно формуле (7)
мультипликативный критерий будем иметь вид:
F ( ) 
A(1  e T )
N
2
 max .
Пример. Применим мультипликативный критерий оптимальности для
определения оптимальных параметров для автомата. Мы получили
следующую задачу оптимизации:
найти максимум функции f ( L, N )  k L * N
при ограничении 1.6L+ 0.05N+2≤6.
Для определения максимального значения f(L,N) с учётом ограничения
равенства мы будем использовать метод неопределённых множителей
Лагранжа. Получим следующую задачу безусловной оптимизации:
найти максимум функции ( L, N ,  )  k L * N   (1.6 L  0.05 N  4).
Находим частные производные по L, N,  и приравниваем их к нулю.
Получим систему уравнений:
 k * N

 1.6  0,
L 2 L

 k L  0.05  0,
N

 1.6 L  0.05 N  4  0.

Решая эту систему, получим следующие значения: Nopt=53, Lopt=0.83 м,
Vopt=137 м/cек.
Использование мультипликативного критерия в задаче оптимизации
привело другим значениям параметров автомата по сравнению с решением
задачи с аддитивным критерием оптимальности. Это объясняется тем, что
диапазоны взаимной компенсации абсолютных и относительных изменений
частных критериев неодинаковы. Поэтому в каждом конкретном случае
технического проектирования следует тщательно анализировать и
обосновывать
целесообразность
учёта
либо
абсолютных,
либо
относительных изменений значений частных критериев и в зависимости от
степени важности этих отклонений выбирать либо аддитивный, либо
мультипликативный критерий оптимальности.
Заключение. Преимущества и недостатки формальных обобщённых
критериев.
Преимущества – возможность учёта в качестве fi(X) любых выходных
параметров системы, а также надёжность и стоимость.
Основные недостатки. – возможность компенсации ухудшения целевой
функции из-за ухудшения одного параметра за счёт улучшения какого-либо
другого выходного параметра,
Метод "идеальной" точки
Метод "идеальной" точки. Рассматривается m-мерное пространство
(где m число локальных критериев), в котором априори выбирается вектор,
отображающий "идеальное" решение (или, что тоже самое, "идеальная"
точка, координатами которой являются "идеальные" значения (например
минимальные или максимальные значения) локальных критериев). В этом
пространстве вводится некоторая метрика, с целью вычисления расстояния
между вектором, отображающим рассматриваемое решения, и "идеальным".
В качестве наилучшего выбирается такое решение, векторная оценка
которого наиболее близка к "идеальной" точке. Недостатками метода
являются произвол при выборе идеальной точки и введение метрики.
Определим обобщенный критерий следующим образом. Положим
ai=maxFi(X); i  1, m , т.е. ai является максимально (минимально) возможным
значением по i – му критерию. Положим a=(a1, a2, . . ., am). Точка a
называется идеальной. Смысл названия связан с тем, что такие точки
оптимальны сразу по всем критериям – получить большее (меньшее)
значение ни по одному критерию невозможно. Как правило, точка aYD.
Зададим для всех точек YYD функцию, являющуюся евклидовым
расстоянием между точками Y и a
1
m
2
 ( y , a )   ( a i  y i ) 2  .
 i 1

За целевую функцию (обобщённый критерий) берут выражение
m
f ( X )   i ai  Fi ( X ) ,
2
i 1
где i – весовые коэффициенты.
Таким образом, задача оптимизации формулируется следующим образом
m
min  i ai  Fi ( X )2
i 1
XD
С учётом нормировки
a  F (X )
min  i  i 0i

Fi
i 1


XD
m
2
Зам. Здесь принцип оптимальности выражается функцией выбора
определяемой близостью к идеальной точке.
Зам. В качестве идеальной точки берут директивные значения параметров,
заданные заказчиком, т.е. в ТЗ (техническом задании).
Какие
задачи
оптимального
проектирования
приводят
к
использованию метода идеальной точки?
Например, когда все или основные условия работоспособности имеют вид
равенств, т.е. Fi(X)=TTi, где TTi – значение технического требования,
предъявленные к i - критерию. Тогда целевая функция имеет вид:
2
 F (X)  TTi 
f (X)    i  i
  min .
Fi0

i 1


m
Пример. Пусть имеются частные критерии F1(X)=-3x1+2x2; F2(X)=4x1+3x2;
F3(X)=2x1-5x2, которые требуется максимизировать. Область D задаётся
неравенствами –x1-3x2+180; -2x1-x2+100; x10; x20. Линейная функция
F1(X) достигает максимального значения a1=12 в точке X1=(0, 6); F2(X) максимальное значение a2=24 в точке X2=(3, 4); F3(X) - максимальное
значение a3=10 в точке X3 =(5, 0). По методу идеальной точки составим
функцию
f(X)=[12-(-3x1+2x2)]2+[24-(4x1+3x2)]2+[10-(2x1-5x2)]2.
После
2
2
преобразований получим f(x1, x2)= 29 x1  38x2  8x1 x2  160 x1  92 x2  820 .
Таким образом, задача оптимизации будет такая
min f(x1, x2)
g1(X)= –x1-3x2+180
g2(X)= -2x1-x2+100
g3(X)=x1≥0; g4(X)=x2≥0
Построим область D.
y=10-2x
Xopt
y=6-2/3x
Область
D
Рис. 8. Область D
Найдем максимум функции, не учитывая ограничений. Если
полученное значение будет лежать в области D, то оно и будет решением
нашей задачи.
Находим частные производные по x1 и x2 функции f(x1, x2) и приравняем их
нулю. Получим следующую систему линейных уравнений
29x1-4x2=80
-4x1+38x2=46.
Решение этой системы: x1=2.97; x2=1.52. Эта точка находится внутри
области D. Следовательно, минимум функции достигается в точке с
координатами x1opt=2.97; x2opt=1.52.
В найденной точке Xopt, являющейся решением рассматриваемой
многокритериальной задачи линейного программирования, F1(Xopt)=-5.87;
F2(Xopt)=16.44; F3(Xopt)=-1.66.
2
2
f ( x y )  29  x  38  y  8  x  y  160  x  92  y  820
x  2
Given
y  1
—целевая функция
— начальные приближения
2  x  3  y  18  0
— ограничения
2  x  y  10  0
x 0
 2.969 

 1.523 
minimize( f  x y )  
— функция минимизации
Лекция. Проблема построения обобщённого критерия для векторных задач
оптимизации
— график целевой функции
Вопросы.
Сложности в построении обобщённого критерия; примеры.
Формальное определение обобщённого критерия. Эквивалентность
обобщённых критериев.
Локальный коэффициент замещения (ЛКЗ).
Карта безразличий. Условия постоянства ЛКЗ. Определённость обобщенного
критерия картой безразличий. Задача.
1.
2. Рассмотрим теперь в общем виде проблему построения обобщённого
критерия для многокритериальных задач принятия решений. Ограничимся
случаем двух критериев, оценки по которым будем обозначать через u и v
соответственно; тогда каждая векторная оценка может быть представлена
точкой на координатной плоскости (u, v).
Download