Анализ устойчивости


Рассмотрим пару двойственных задач. Задачу об оптимальном плане выпуска продукции (1) и задачу определения оптимальных оценок ресурсов (2).

З.1.clip_image002 (1) clip_image004 (1`)

clip_image006 (2) clip_image008 (2`)

clip_image010 (3) clip_image012 (3`)

Решим задачу (1) симплекс методом и найдем ее оптимальное решение и максимальное значение целевой функции clip_image014 и найдем решение двойственной задачи и минимальное значение целевой функции clip_image016. clip_image018 clip_image020

Под анализом устойчивости оптимального плана производства будем понимать решение трех вопросов:

1.Как влияет на максимальное значение целевой функции «малое» изменение какого-либо i-того ресурса.

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

3.Как оценить целесообразность выпуска какого-либо вида продукции.

На 1 и 3 вопросы можно дать ответы решив задачу заново с измененной информацией, но этого делать не надо. Ответ на эти вопросы можно получить используя имеющуюся информацию решения задачи

1.Как влияет на максимальное значение целевой функции «малое» изменение какого-либо i-того ресурса.

Обозначим через clip_image022 вектор ресурсов исходной задачи (1):

clip_image024, а максимальное значение целевой функции полученное при решении этой задачи через clip_image026. Изменим объем i-го ресурса на clip_image028. Вектор изменения ресурса обозначим clip_image030= clip_image032. Новый вектор ресурсов clip_image034. Пусть при этом новом векторе ресурсов план clip_image036 будет допустимым планом. При этом новом векторе ресурсов значение целевой функции изменяется на величину clip_image038. Найдем отношение clip_image040 и перейдем к пределу clip_image042 clip_image044. clip_image046

Если clip_image048, clip_image050, то при увеличении i-того ресурса на clip_image052 и при неизменных всех остальных ресурсах максимальная прибыль увеличится на clip_image054.

Если clip_image056, clip_image058, то уменьшится на clip_image059.

Если i-ый ресурс недефицитный, то его оценка clip_image061, следовательно «малое» изменение недефицитного ресурса на максимальное значение целевой функции не влияет.

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

Обозначим исходный вектор ресурсов через clip_image024[1]. Изменим объем ресурса вида Р на величину clip_image063, тогда вектор изменения ресурсов clip_image065=clip_image067. Новый вектор ресурсов будет равен clip_image069 или clip_image071 clip_image073 =clip_image075. Решим исходную задачу (1) симплексным методом и в последней оптимальной симплексной таблице, в последнем столбце, получим преобразованный по формулам полного исключения неизвестных вектор, который обозначим через clip_image077. Если одновременно в симплексных таблицах преобразовывать clip_image079, то в последней симплексной таблице получим преобразованный вектор clip_image081. Можно показать clip_image083 (4).

Вывод формулы (4):

Известно, что при решении ЗЛП симплексным методом правые части ограничений на каждом шаге преобразуются по формулам:

clip_image085; clip_image087 ;clip_image089 (*)

где r – разрешающая строка, а s – разрешающий столбец. Если правая часть состоит из суммы двух векторов:

clip_image065[1]=clip_image091 и clip_image093,

т.е. clip_image095 =clip_image097, то применяя формулы (*), получим:

clip_image099clip_image101 clip_image103;clip_image087[1]

Т.о., чтобы преобразовать вектор clip_image105 по формулам (*), достаточно преобразовать векторы clip_image107 в векторы clip_image109 и найти выражение clip_image111. В нашем случае clip_image113.

Решив (1) симплексным методом в последней таблице мы найдем вектор clip_image114 и clip_image116. Поэтому, не решая задачу с новым вектором clip_image118, мы можем найти его величину относительно оптимального базиса задачи (1) используя формулу (4) решение задачи (1). Чтобы решить задачу (1) симплексным методом, приведем ее к стандартному виду и запишем условие в исходную симплексную таблицу.

Век

A1

An

E1

Ep

Em

B

БП

X1

Xn

Xn+1

Xn+p

Xn+m

bi

Хn+1

a11

a1n

1

0

0

b1

Хn+2

a22

a2n

0

0

0

b2

Хn+p

ap1

apn

0

1

0

bp

Хn+m

am1

amn

0

0

1

bm

Z

clip_image1201

clip_image120[1]n

0

0

0

0

Решив задачу симплексным методом мы получим:

Век

E`1

E`p

E`m

B`

БП

X1

Xn

Xn+1

Xn+p

Xn+m

С.Ч.

Xi1

clip_image122

b`1

Xi2

b`2

Xim

b`m

Z

D`1

D`n

D`n+1

D`n+p

D`n+m

Q

В этой таблице все оценки D`j>=0; j=1…m+n; b`i>=0; i=1…m. БП Хi1…Хim. Решение задачи X*i1=b`1,…, X*im =b`m и Zmax=Q. clip_image124 и clip_image126. Все остальные свободные переменные равны нулю. И тогда в оптимальном базисе задачи (1): clip_image128+clip_image130=clip_image132. Для того, чтобы оптимальный базис задачи (1) оставался оптимальным и для clip_image134, необходимо чтобы значения БП: Хi1,…,Xim были неотрицательны, то есть Xi1=b`1+Dbp*q1>=0,…,Xim=b`m+Dbp*qm>=0, то есть необходимо решить следующую систему неравенств b`i+Dbp*qi>=0, где i=1,…,m. Решив эту систему неравенств мы найдем: clip_image136

Т.о. номенклатура выпускаемой продукции и состав БП в оптимальном плане задачи (1) не изменится, если . При изменении i-того ресурса на величину Dbi: clip_image138. Интервал (clip_image140) наз. Интервалом устойчивости оценки y*i. Если объем i-того ресурса менять в этом интервале, а объемы всех остальных ресурсов не менять, то в оптимальном плане задачи (1) сохраняется состав БП, т.е. номенклатура выпускаемой продукции, т.е. оценки всех ресурсов не меняются. Если объем i-того ресурса изменяется за пределами интервала устойчивости оценок, то состав bp меняется, а из-за этого изменяется оценки Z строки и следовательно изменяется значение оценок y*i. Изменение максимальной прибыли DZmax при изменении i-того ресурса на малую величину Dbi можно определять по формуле: DZmax=Dbi· y*i, если Dbi принадлежит интервалу clip_image142.

3.Как оценить целесообразность выпуска какого-либо вида продукции.

Согласно второй теореме двойственности оптимальные решения clip_image144 и clip_image146 должны удовлетворять условиям: clip_image148, тогда для тех видов продукции, которые выпускаются, стоимость ресурсов идущих на изготовление ед. продукции вида j должна быть в точности равна прибыли, от ед. этой продукции, т.е. clip_image150.

Если clip_image151, то такой вид продукции выпускать не выгодно. Пусть теперь имеется новый s-тый вид продукции с нормами затрат ресурсов на ед. продукции ais ,i=1,…,m и прибылью от ед. продукции сs. Тогда стоимость всех ресурсов в оптимальных оценках, идущих на изготовление ед. продукции вида s будет равна: clip_image153

И если эта величина будет < прибыли, то тогда данный вид продукции выпускать целесообразно, а если больше прибыли — не целесообразно. Т.о. двойственные оценки y*i формируют показатель clip_image155. Если Dj>0 то j–тый вид продукции выпускать невыгодно. А если Dj<0 то j–тый вид продукции выпускать выгодно. А для тех видов продукции, которые выпускаются Dj=0.

Ответ на третий вопрос.

Согласно второй теореме двойственности оптимальное решение clip_image157и clip_image159должны удовлетворять условиям:

clip_image161

тогда для всех видов продукции, которые выпускаются, т.е. clip_image163 стоимость ресурсов идущих на изготовление единицы продукции вида j должна быть равна прибыли от реализации единицы продукции т.к. clip_image165

если стоимость всех ресурсов идущих на изготовление единицы продукции вида j больше прибыли от реализации единицы этой продукции, то такой вид продукции выпускать невыгодно, то clip_image167

Пусть имеем новый S-й вид продукции с нормами затрат ресурсов на единицу продукции clip_image169 и прибылью clip_image171, тогда стоимость всех ресурсов в оптимальных оценках идущих на изготовление единицы продукции вида S будет выражаться следующей величиной:clip_image173

и если эта величина будет меньше прибыли, то тогда данный вид продукции выпускать целесообразно, а если эта величина больше или равна прибыли, то нецелесообразно. Таким образом, двойственные оценки clip_image175 формируют показатель:clip_image177

Если clip_image179, то j-й вид продукции выпускать невыгодно, а если clip_image181 — то выгодно.

Загрузка...