- •Восстановление весов частных критериев оптимальности по принятым решениям
- •II. Восстановление весов частных критериев по одному наблюдению
- •III. Восстановление весов частных критериев на основании серии наблюдений
- •IV. Построение прогноза на четвертый период
- •V. Восстановление весов частных критериев при отсутствии активных ограничений
-
Восстановление весов частных критериев оптимальности по принятым решениям
В общем случае задача может быть представлена следующим образом:
, где fk – частные критерии оптимальности, х – вектор-столбец инструментальных переменных, А – технологическая матрица, b – вектор-столбец ограничений ресурсов.
Если считать, что fk – линейна, т.е. fk (х) = (ск)Тх, то исходная задача может быть представлена в следующем виде:
, где к – искомые веса частных критериев ск.
По условию Куна-Такера необходимые и достаточные условия достижения максимума функции в точке хopt имеют вид:
, где - градиент целевой функции по направлению и , - i – тая нормаль для ограничений A*x b и , уi – i-тая двойственная оценка, vj – неотрицательное число, е – единичная компонента.
Таким образом, получаем следующую систему соотношений:
В результате могут быть получены следующие ситуации:
-
Решение единственно.
-
Множество решений.
-
Пустое множество.
II. Восстановление весов частных критериев по одному наблюдению
Приложение 1.
Исходные данные. 3-й период.
Рассматриваются 2 частных линейных критерия оптимальности, которым соответствуют вектора коэффициентов:
В этих условиях был выбран план:
Множество активных ограничений : I={1,2,3}
Множество активных ограничений : J={2,5}
В данном случае система уравнений записывается следующим образом:
Так как уi могут быть любого знака, то введем уi = уi+ - уi-. Запишем полученную систему в явном виде.
Для этого определили множество I(xopt), куда войдут ограничения типа А*хopt=b. Сама система будет записана следующим образом:
60 |
* |
1 |
+ |
30 |
* |
2 |
|
-7 |
* |
у1 |
|
|
-6 |
* |
у2 |
|
-13 |
* |
у3 |
|
|
|
+ |
q1 |
= |
0 |
11 |
* |
1 |
+ |
69 |
* |
2 |
-5 |
* |
у1 |
-6 |
* |
у2 |
-20 |
* |
у3 |
+ |
2 |
+ |
q2 |
= |
0 |
|||||
25 |
* |
1 |
+ |
50 |
* |
2 |
-10 |
* |
у1 |
0 |
* |
у2 |
-17,5 |
* |
у3 |
|
|
+ |
q3 |
= |
0 |
|||||
10 |
* |
1 |
+ |
20 |
* |
2 |
0 |
* |
у1 |
-5 |
* |
у2 |
0 |
* |
у3 |
+ |
|
+ |
q4 |
= |
0 |
|||||
10 |
* |
1 |
+ |
28,8 |
* |
2 |
-2 |
* |
у1 |
-1 |
* |
у2 |
-13,8 |
* |
у3 |
|
5 |
+ |
q5 |
= |
0 |
1 + 2 = 1
к 0, к=1,2
yis 0, iI
j 0, jJ
Следует учесть, что в результате решения может быть получено множество решений для к. Условием не единственности полученного решения является наличие небазисной переменной, для которой оценочный коэффициент равен нулю, т.е. внесение этой переменной в базис не повлияет на значение целевой функции. В рассматриваемой системе предлагается использовать следующую логику рассуждений для определения иных базисных решений:
появление в решении небазисных переменных с равными нулю оценочными коэффициентами свидетельствует, как уже было описано выше, о наличии других решений системы;
признаками принадлежности исходного решения базисному (опорному плану) являются j , причем j = 0 – переменная базисная, j > 0 – переменная не базисная;
на основе выше сказанного, переход от одного базисного решения к другому предлагается осуществлять простым исключением (приравниванием нулю) какого-либо значения j из предлагаемой и системы, тем самым, включая соответствующий xj в базис;
критерием включения того или иного хj в базис будет отсутствие в базисе искусственно введенных переменных qi , так как по условию они должны быть исключены из опорного плана (j в этом случае должно на себя оттягивать любое значение с qi), и, следовательно, включать соответствующую j - тую технологию в первоначальный оптимальный план нельзя.
Таким образом, в результате проведенного анализа может быть получено множество такое, что , 0 1.
Используя предложенный способ анализа в пакете линейного программиования BLP 88, были получены следующие результаты:
-
11 = 0., 21 = 0., при 1 = 12.14, 4 = 6.43 и невошедших в базисное решение y11 и у12, т.е. опорный план (х2,х3,х5) ;
-
12 = 0.4988, 22 = 0.5012, при 2 = 6,882 и не вошедшим в базисное решение 5.
Тогда искомое решение будет представлено следующим образом: