АкушерствоАнатомияАнестезиологияВакцинопрофилактикаВалеологияВетеринарияГигиенаЗаболеванияИммунологияКардиологияНеврологияНефрологияОнкологияОториноларингологияОфтальмологияПаразитологияПедиатрияПервая помощьПсихиатрияПульмонологияРеанимацияРевматологияСтоматологияТерапияТоксикологияТравматологияУрологияФармакологияФармацевтикаФизиотерапияФтизиатрияХирургияЭндокринологияЭпидемиология

Выполнение теоремы об оценке

Прочитайте:
  1. E Выполнение приема Геймлица
  2. II. Выполнение процедуры
  3. II. Выполнение процедуры
  4. II. Выполнение процедуры
  5. II. Выполнение процедуры
  6. II. Выполнение процедуры
  7. II. Выполнение процедуры.
  8. II. Выполнение процедуры.
  9. Алгоритм по оценке степени тяжести вирусного гепатита
  10. Анализ некоторых конкретных проблем, подлежащих переоценке

В условиях задачи сказано, что недельный рацион должен содержать не менее 4, но не более 6 г кальция; не менее 110 г белка; не более 25 г клетчатки; и этой смеси должно быть ровно 500 г. Предположим, что эти величины могут меняться. Чтобы узнать, на сколько при этом изменится оптимальная стоимость рациона, в соответствии с третьей теоремой двойственности надо воспользоваться теневыми ценами y1-5, которые соответствуют ограничениям прямой задачи. Но для этого вначале следует провести анализ устойчивости теневых цен.

 

Обозначим изменение минимального содержания кальция Db1; максимального - Db2; минимального содержания белка - Db3; максимального содержания клетчатки - Db4, а массы смеси – Db5 (все эти величины измеряются в граммах). Теперь в таблице 17 вместо свободных членов (b1, b2; b3; b4; b5) следует подставить свободные члены (b1 + Db1, b2 + Db2; b3 + Db3; b4 + Db4; b5 + Db5), т.е. (4 + Db1, 6 + Db2; 44 + Db3; 25 + Db4; 0,5 + Db5).

Однако, если мы просто заполним диапазон ячеек электронной таблицы D3:D7 этими выражениями, программа Microsoft Excel не сможет осуществить никаких вычислений над ними, поскольку эти ячейки станут текстовыми. Чтобы избежать этого, вставим перед столбцом Е еще пять столбцов. Теперь столбец свободных членов будет занимать 6 столбцов электронной таблицы (D, E, F, G, H, I). В столбце E запишем коэффициент, который будет умножен на Db1, в столбце F – коэффициент при Db2, и т.д., в столбце I – коэффициент при Db5, а в D – то слагаемое в выражении для свободного члена, которое ни на что не умножается (т.е. прежнее значение свободного члена). Результат представлен в таблице 22.

 

Теперь новые значения свободных членов, т.е. новый столбец B, записанный в диапазоне D3:I7 электронной таблицы, необходимо подвергнуть тем же линейным преобразованиям, которым подвергались ограничения прямой задачи в таблицах 17-19. Для этого нужно выделить диапазон ячеек D8:D37 и скопировать его на диапазон Е8:I 37. В результате этого новые столбцы симплексных таблиц будут пересчитаны по тем же формулам, что и столбец D. Результаты вычислений приведены в последней строке таблицы 22 и в таблице 23. Заголовки столбцов в таблице 23 отредактированы.

 


Таблица 22 – Подготовка исходной симплексной таблицы к проведению анализа устойчивости двойственных оценок

  A B C D E F G H I J K L M N
                             
  N xб cб B Db1 Db2 Db3 Db4 Db5 x1 x2 x3 x4 x5
    у1                     -1  
    x5                        
    у2                        
    x7                        
    у3   0,5                    
  m+1     48,5                 -1  

 

Из таблицы 23 видно, что во второй симплексной таблице теперь базисная искусственная переменная у1 = 3,511 + Db1 - 0,011Db3 (см. строку 10 электронной таблицы); а не 3,511, как в таблице 17. Базисная переменная x5 =
= 5,511 + Db2 - 0,011Db3 (см. строку 11 электронной таблицы); а не 5,511, как в таблице 17, и т.д.

Оптимальный план прямой задачи примет вид Х* = (0,011 - 0,011Db3 +
+ Db5; 0,489 + 0,011Db3; 0; 0,711 - Db1 - 4,211Db3 + 380Db5; 1,289 + Db2 +
+ 4,211Db3 - 380Db5; 0; 15,222 - 0,022Db3 + Db4), оптимум будет равен 7,378 + 0,122Db3 + 4Db5.

 

Таблица 23 – Анализ устойчивости двойственных оценок

  A B C D E F G H I J K L
  N xб cб B Db1 Db2 Db3 Db4 Db5 x1 x2 x3
    у1   3,511     -0,011         1,444
    x5   5,511     -0,011         1,444
    x2   0,489     0,011         0,556
    x7   15,222     -0,222         68,889
    у3   0,011     -0,011         0,444
  m+1     3,522     -0,022         1,889
  N xб cб B Db1 Db2 Db3 Db4 Db5 x1 x2 x3
    x1   0,009 0,003             0,004
    x5     -1              
    x2   0,489     0,011         0,556
    x7   15,222     -0,222         68,889
    у3   0,002 -0,003   -0,011         0,441
  m+1     0,002 -0,003   -0,011         0,441

Продолжение таблицы 23

                         
  N xб cб B Db1 Db2 Db3 Db4 Db5 x1 x2 x3
    x1   0,009 0,003       -0,009      
    x5     -1              
    x2   0,487 0,003   0,025   -1,261      
    x7   14,93 0,41   1,51   -156,34      
    х3   0,004 -0,006   -0,025   2,269      
  m+1     7,505 -0,179   -0,629   71,830      
  N xб cб B Db1 Db2 Db3 Db4 Db5 x1 x2 x3
    x1   0,011     -0,011         0,444
    x5   1,289     4,211   -380     -167,444
    х2   0,489     0,011         0,556
    x7   15,222     -0,222         68,889
    х4   0,711 -1   -4,211         167,444
  m+1     7,378     0,122         -29,889

 

Полученный ответ будет иметь смысл лишь в том случае, если значения всех переменных будут неотрицательны (т.е. последняя таблица будет допустимой). Поэтому необходимо проверить знак базисных переменных x1, x5, х2, x7, х4; т.е. определить, при каких изменениях свободных членов они будут неотрицательны. Составим систему неравенств:

0,011 - 0,011Db3 + Db5 ³ 0

0,489 + 0,011Db3 ³ 0

0,711 - Db1 - 4,211Db3 + 380Db5 ³ 0

1,289 + Db2 + 4,211Db3 - 380Db5 ³ 0

15,222 - 0,022Db3 + Db4 ³ 0

 

Чтобы решить эту систему, вначале предположим, что меняется только минимальное содержание кальция, т.е. первый свободный член. Тогда Db2 = Db3 = Db4 = Db5 = 0, и система примет вид:

0,011 ³ 0

0,489 ³ 0

0,711 - Db1 ³ 0

1,289 ³ 0

15,222 ³ 0

Отсюда Db1 £ 0,711. В первоначальном варианте исходных данных b1 =
= 4. Следовательно, двойственная оценка y1 = 0 будет устойчивой лишь при норме минимального содержания кальция от любого значения до 4,711 г (4 +
+ 0,711 = 4,711). Можно сказать, что оптимальная стоимость корма не изменится (изменится на 0 руб.), если потребуется, чтобы в корме содержалось, например, не менее 4,5 г кальция, или не менее 4,7 г кальция, или не менее 2 г кальция, или будет допускаться, чтобы вообще не содержалось кальция. Если эта норма будет повышена, например, до 5 г кальция (5,5 г, 4,8 г, 100 г и т.д.), то неизвестно, на сколько изменится оптимум задачи.

 

Теперь предположим, что меняется только максимальное содержание кальция, т.е. второй свободный член. Тогда Db1 = Db3 = Db4 = Db5 = 0, и система примет вид:

0,011 ³ 0

0,489 ³ 0

0,711 ³ 0

1,289 + Db2 ³ 0

15,222 ³ 0

Отсюда Db2 ³ -1,289. Следовательно, двойственная оценка y2 = 0 будет устойчивой при норме максимального содержания кальция от 4,711 г
(6 - 1,289 = 4,711) до любого значения. Можно сказать, что оптимальная стоимость корма не изменится, если потребуется, чтобы в корме содержалось, например, не более 4,8 г кальция, или не более 10 г кальция, или будет допускаться, чтобы кальция было сколь угодно много. Если эта норма будет снижена, например, до 4,5 г кальция (4,7 г, 1 г, 0 г и т.д.), то неизвестно, на сколько изменится оптимум задачи.

 

Если предположить, что меняется только минимальное содержание белка, то Db1 = Db2 = Db4 = Db5 = 0, и система примет вид*:

Следовательно, двойственная оценка y3 = 0,122 будет устойчивой при норме минимального содержания белка от 43,694 г до 44,169 г. Можно сказать, что если потребуется, чтобы в корме содержалось, например, не менее 43,7 г белка, то оптимальная стоимость рациона уменьшится на 0,122*(44 – 43,7) = 0,122*0,3 = 0,0366 (руб.). Если граничное содержание белка составит 44,1 г, то эта стоимость возрастет на 0,122*(44,1 – 44) = 0,0122 (руб.); и т.п. Но если граничное содержание белка составит, например, 43 г, или 10 г, или 44,2 г, или 100 г и т.п., то неизвестно, на сколько изменится оптимум задачи, поскольку при этом исходные данные выйдут за интервал устойчивости двойственных оценок.

 

При изменении нормы содержания клетчатки Db1 = Db2 = Db3 = Db5 = 0, и система примет вид:

0,011 ³ 0

0,489 ³ 0

0,711 ³ 0

1,289 ³ 0

15,222 + Db4 ³ 0

Отсюда Db4 ³ -15,222. Следовательно, двойственная оценка y4 = 0 будет устойчивой при норме содержания клетчатки от 9,778 г до любого значения.

 

И, наконец, при изменении заданной массы корма Db1 = Db2 = Db3 =
= Db4 = 0, и система примет вид:

0,011 + Db5 ³ 0

0,489 ³ 0

0,711 + 380Db5 ³ 0

1,289 - 380Db5 ³ 0

15,222 ³ 0

Отсюда -0,003 £ Db5 £ 0,002. Следовательно, двойственная оценка y5 =
= 4 будет устойчивой при массе рациона от 0,497 г до 0,502 г.

 

Правильность полученных результатов легко проверить, изменяя столбец свободных членов в таблице 22 (ячейки D3:D7). Пока они остаются в пределах найденных интервалов, последняя симплексная таблица (см. таблицу 23) останется допустимой. При этом значения двойственных оценок (коэффициенты в критериальной строке) не изменятся, а оптимум (свободный член в критериальной строке) изменится в соответствии с теоремой об оценке.

Например, подставим в ограничение по содержанию белка вместо свободного члена 44 найденное граничное значение 43,694. Для этого просто введем в ячейку D5 это число. При этом результаты вычислений в других ячейках столбца D изменятся. Новые значения в оптимальной симплексной таблице приведены в таблице 24.

 


Таблица 24 – Проверка результатов анализа устойчивости

  A B C D
  N xб cб B
    x1   0,015
    x5    
    х2   0,485
    x7   15,290
    х4    
  m+1     7,34

 

Итак, решение задачи изменилось. Новое значение оптимума 7,34. Свободный член изменился на (43,694 - 44 = -0,306), поэтому оптимум должен был измениться на у4*(-0,306) = 0,122*(-0,306)» -0,037. При первоначальном варианте исходных данных оптимум составлял 7,378. Поскольку 7,378 - 0,037 = 7,34, проверка подтвердила правильность всех вычислений.

 

Кроме того, результаты вычислений в ячейках D37:I37 таблицы 23 как раз показывают, что оптимум задачи равен 7,378 + 0,122Db3 + 4Db5, т.е. непосредственно иллюстрируют теорему об оценке.

 

Следует обратить внимание, что при подстановке граничных величин свободных членов оптимальный план становится вырожденным (т.е. по крайней мере одна из базисных переменных равна нулю). Например, в таблице 24 x5 = 0. При выходе за пределы найденных интервалов допустимость заключительной таблицы нарушится, т.е. по крайней мере одно из значений в столбце В (ячейках D32:D36) станет отрицательным. Кроме того, следует помнить, что эти интервалы найдены для условий, когда одновременно изменяется только один элемент столбца свободных членов, а остальные остаются неизменными.

 

Отметим, что если сравнить столбец F3:F37 (коэффициенты при Db2) в таблицах 22-23 со столбцом коэффициентов при х5 I3:I37 (после вставки пяти новых столбцов ячейки I3:I37 стали занимать другой диапазон - N3:N37) в таблицах 17-19, можно убедиться, что они полностью совпадают. Это не удивительно, так как во втором ограничении дополнительная переменная х5 была базисной в исходном опорном плане, и содержимое ячеек F3:F7 в таблице 22 полностью совпадает с содержимым ячеек N3:N7. Следовательно, совпадает и результат их линейного преобразования. То же самое можно сказать про Db4 и х7.

 

Вопросы и упражнения

1 Определите понятие двойственности.

2 По каким правилам строится двойственная задача?

3 Сформулируйте основную теорему двойственности.

4 Сформулируйте следствие из основной теоремы двойственности.

5 Если ОДП одной из сопряженных задач пуста, какие выводы можно сделать о другой задаче?

6 Как получить решения двойственной задачи из решения исходной задачи симплекс-методом?

7 Сформулируйте теорему о равновесии.

8 Сформулируйте следствие из теоремы о равновесии.

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

10 В чем заключается экономическая интерпретация теоремы о равновесии?

11 Сформулируйте теорему об оценке.

12 Что такое теневая цена?

13 Что такое сокращенные затраты?

14 В чем заключается понятие устойчивости двойственных оценок?

15 В чем заключается понятие постоптимизационного анализа?

16 Нельзя ли сделать некоторые выводы (и какие) о решении двойственной задачи, если известно, что оптимальный план исходной задачи - вырожденный?

17 Построить и решить задачу, двойственную задаче 11 в разделе 3.4 (пример 3).

18 Построить задачу, двойственную задаче 11 в разделе 3.4 (пример 1), решить ее симплекс-методом и извлечь из заключительной симплексной таблицы решение прямой задачи.

19 Извлечь из решения задачи 11 (раздел 3.4, пример 1) решение двойственной задачи; проверить полученное решение с точки зрения теоремы о равновесии; провести анализ устойчивости двойственных оценок и проверить его результаты.

20 Построить и решить задачу, двойственную задаче 12 в разделе 3.4, решить ее симплекс-методом и извлечь из заключительной симплексной таблицы решение прямой задачи.

21 Извлечь из решения задачи 12 (раздел 3.4) решение двойственной задачи; проверить полученное решение с точки зрения теоремы о равновесии; провести анализ устойчивости двойственных оценок и проверить его результаты.

 


* Не следует путать переменные Y двойственной задачи ни с дополнительными, ни с искусственными переменными, которые в предыдущих разделах иногда обозначались той же буквой.

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

* Это следует из физического смысла производной (скорость изменения функции в точке).

* В разделе 4.2 обозначения y1, y2 и y3 использовались для искусственных переменных при построении расширенной задачи. Поскольку при изучении теории двойственности значения этих переменных использоваться почти не будут, нет необходимости использовать другие буквы латинского алфавита (хотя можно это сделать). Тем не менее, сделаем замечание о том, что не следует путать двойственные переменные с искусственными переменными, которые использовались при решении прямой задачи методом искусственного базиса.

* Если при выполнении студентами упражнений для решения двойственной задачи потребуется применить метод искусственного базиса, рекомендуется обозначить искусственные переменные z1, z2 и т.д.

* В таблице 19 нет оптимального значения двойственной переменной у5. Это объясняется тем, что она соответствует уравнению прямой задачи, в котором, естественно, нет дополнительной переменной.

Однако при решении задачи симплекс-методом в разделе 4.2 в это ограничение вводилась искусственная переменная (ей соответствовал последний столбец N в таблицах 17-19). Отметим, что таблице 19 соответствуют строки 31-37 электронной таблицы. После того, как расширенная задача была решена, все искусственные переменные были исключены из рассмотрения и соответствующие столбцы коэффициентов более не пересчитывались. Если пересчитать их по тем же формулам, что и остальные столбцы (для этого в таблицах 18 и 19 достаточно скопировать формулы в диапазоне К25:К37 на диапазон L25:N37), то в последнем столбце критериальной строки (в N37) коэффициент будет равен 4 (у5 = 4).

* Чтобы вычислить по модулю правые части неравенств во второй системе, нет необходимости проводить расчеты вручную. Можно ввести в любую свободную ячейку, например, в ячейку Т32 формулу =D32/G32, которую затем скопировать на диапазон Т33:T36.

 


Дата добавления: 2015-01-18 | Просмотры: 478 | Нарушение авторских прав



1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |



При использовании материала ссылка на сайт medlec.org обязательна! (0.018 сек.)