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

Пример решения задачи методом искусственного базиса

Прочитайте:
  1. HLA- DR, DQ, DP. В этой же зоне находятся и другие гены: например, DN, DO, продукты которых пока не известны.
  2. I. Иммунология. Определение, задачи, методы. История развитии иммунологии.
  3. II -А. Задачи СИТУАЦИОННЫЕ по диагностике в
  4. II. Основные задачи
  5. II. Целевые задачи
  6. II. Целевые задачи
  7. II. Целевые задачи
  8. II. Целевые задачи
  9. II. Целевые задачи
  10. II. Целевые задачи

Решим двухэтапным методом искусственного базиса задачу из раздела 1.3.1. Для этого ее надо вначале привести к канонической форме:

min 4х1 + 15х2 + 40х3

380х1 + х2 + 2х3 – х4 = 4

380х1 + х2 + 2х3 + х5 = 6

90х2 + 50х3 – х6 = 110

20х2 + 80х3 + х7 = 25

х1 + х2 + х3 = 0,5

х1-7 ³ 0

Здесь дополнительная переменная х4 представляет собой превышение содержания кальция в рационе над минимально допустимым, а х5 показывает, на сколько оно ниже максимально допустимого. Переменная х6 показывает, на сколько содержание белка в рационе больше минимально допустимого, а х7 – на сколько меньше содержание клетчатки по сравнению с максимально допустимым. Все эти величины измеряются в граммах.

 

В этой задаче свободные члены неотрицательны. Имеется часть базиса, а именно единичные столбцы стоят при двух переменных – х5 и х7 (во втором и четвертом ограничениях). В остальных трех уравнениях базисных переменных нет, поэтому введем три искусственные переменные и построим расширенную задачу для реализации двухэтапного симплекс-метода следующим образом:

min у1 + у2 + у3

380х1 + х2 + 2х3 – х4 + у1 = 4

380х1 + х2 + 2х3 + х5 = 6

90х2 + 50х3 – х6 + у2 = 110

20х2 + 80х3 + х7 = 25

х1 + х2 + х3 + у3= 0,5

х1-7 ³ 0

у1-3 ³ 0

Решим расширенную задачу обычным симплекс-методом. Тогда исходная таблица примет вид таблицы 14.

 


Таблица 14 – Исходная симплексная таблица для расширенной задачи

                           
N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
  у1           -1            
  x5                        
  у2               -1        
  x7                        
  у3   0,5                    
m+1     114,5       -1   -1        

 

Так как задача на минимум, здесь критерий оптимальности нарушен в столбцах при трех первых переменных (в критериальной строке в них стоят положительные коэффициенты). В качестве разрешающего можно выбрать любой из них. Выберем, например, столбец x2. Тогда из базиса выйдет переменная у3. Следующая таблица будет иметь вид таблицы 15.

 

Таблица 15 – Вторая симплексная таблица для расширенной задачи

N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
  у1   3,5       -1           -1
  x5   5,5                   -1
  у2     -90   -40     -1       -90
  x7     -20                 -20
  x2   0,5                    
m+1     68,5     -39 -1   -1       -92

 

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

 

Таблица 16 – Оптимальная симплексная таблица для расширенной задачи

N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
  x1   0,009     0,003 -0,003       0,003   -0,003
  x5                   -1    
  у2   65,831     -39,763 -0,237   -1   0,237   -90,237
  x7   15,185     60,053 -0,053       0,053   -20,053
  x2   0,491     0,997 0,003       -0,003   1,003
m+1     65,831     -39,763 -0,237   -1   -0,763   -91,237

 

Эта таблица оптимальна, так как положительных коэффициентов в критериальном ограничении нет. Оптимальный план расширенной задачи
(x1, x2, x3, x4, x5, x6, x7, у1, у2, у3) = (0,009; 0,491; 0; 0; 2; 0; 15,185; 0; 65,831; 0). Оптимум расширенной задачи равен у1 + у2 + у3 = 65,831 > 0.

Так как оптимум положителен, можно сделать вывод, что исходная задача неразрешима, так как ее ОДП пуста.

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

 

Предположим теперь, что требования к рациону изменились, и теперь предельное содержание белка в нем составляет не 110 г, а всего 44 г. Тогда модель в канонической форме примет вид:

min 4х1 + 15х2 + 40х3

380х1 + х2 + 2х3 – х4 = 4

380х1 + х2 + 2х3 + х5 = 6

90х2 + 50х3 – х6 = 44

20х2 + 80х3 + х7 = 25

х1 + х2 + х3 = 0,5

х1-7 ³ 0

Решим новую задачу симплекс-методом. Для этого расширенную задачу построим следующим образом:

min у1 + у2 + у3

380х1 + х2 + 2х3 – х4 + у1 = 4

380х1 + х2 + 2х3 + х5 = 6

90х2 + 50х3 – х6 + у2 = 44

20х2 + 80х3 + х7 = 25

х1 + х2 + х3 + у3= 0,5

х1-7 ³ 0

у1-3 ³ 0

Решение расширенной задачи представлено в таблице 17 (слева и сверху указаны номера строк и столбцов электронной таблицы Microsoft Excel).

 

В результате решения расширенной задачи был получен ее оптимум, равный нулю, и оптимальный план (x1, x2, x3, x4, x5, x6, x7, у1, у2, у3) = (0,009; 0,487; 0,004; 0; 2; 0; 14,93; 0; 0; 0). Следовательно, можно перейти ко второму этапу двухэтапного симплекс-метода, т.е. начать решать исходную задачу, начиная с опорного плана (x1, x2, x3, x4, x5, x6, x7) = (0,009; 0,487; 0,004; 0; 2; 0; 14,93).

 

 


Таблица 17 – Решение расширенной задачи

  A B C D E F G H I J K L M N
                             
  N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
    у1           -1            
    x5                        
    у2               -1        
    x7                        
    у3   0,5                    
  m+1     48,5       -1   -1        
  N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
    у1   3,511     1,444 -1   0,011     -0,011  
    x5   5,511     1,444     0,011     -0,011  
    x2   0,489     0,556     -0,011     0,011  
    x7   15,222     68,889     0,222     -0,222  
    у3   0,011     0,444     0,011     -0,011  
  m+1     3,522     1,889 -1   0,022     -1,022  
  N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
    x1   0,009     0,004 -0,003       0,003    
    x5                   -1    
    x2   0,489     0,556     -0,011     0,011  
    x7   15,222     68,889     0,222     -0,222  
    у3   0,002     0,441 0,003   0,011   -0,003 -0,011  
  m+1     0,002     0,441 0,003   0,011   -1,003 -1,011  
  N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
    x1   0,009       -0,003       0,003   -0,009
    x5                   -1    
    x2   0,487       -0,003   -0,025   0,003 0,025 -1,261
    x7   14,93       -0,411   -1,51   0,411 1,51 -156,337
    х3   0,004       0,006   0,025   -0,006 -0,025 2,269
  m+1                     -1 -1 -1

 

Иными словами, переходят к решению следующей задачи:

min 4х1 + 15х2 + 40х3

х1 – 0,003х4 = 0,009

х4 + х5 = 2

х2 – 0,003х4 – 0,025х6 = 0,487

-0,411х4 - 1,51х6 + х7 = 14,93

х3 + 0,006х4 + 0,025х6 = 0,004

х1-7 ³ 0

Именно такая система уравнений записана в последней симплексной таблице (она представляет собой систему уравнений исходной задачи, подвергнутую линейным преобразованиям). Здесь присутствует полный набор базисных переменных: x1, x5, x2, x7 и х3.

 

Однако последняя строка таблицы 17, полученная методом Гаусса, соответствует критериальному ограничению расширенной задачи, т.е. другой целевой функции. Чтобы перейти к исходной задаче, это ограничение надо построить заново, т.е. пересчитать его свободный член (значение целевой функции) и коэффициенты по формулам, выведенным в разделе 3.2.2 (в самом деле, ведь базисные переменные опорного плана здесь входят в целевую функцию).

Для этого в столбец cб записывают коэффициенты целевой функции исходной задачи при переменных x1, x5, x2, x7 и х3: 4, 0, 15, 0 и 40. Значение целевой функции составит 4*0,009 + 15*0,487 + 40*0,004» 7,5. Первые три коэффициента в критериальном ограничении будут нулевыми, так как переменные x1-3 – базисные, D4 = 4*(-0,003) + 15*(-0,003) + 40*0,006 – 0» 0,2* и так далее.

 

При использовании для расчетов электронной таблицы рекомендуется вставить перед последней симплексной таблицей (перед строкой 23) еще одну строку, и ввести в нее коэффициенты целевой функции исходной задачи. При этом все последняя таблица сместится на строку вниз и займет диапазон ячеек А24:N30 (см. таблица 18). Новая строка станет строкой № 23, и коэффициентами целевой функции будет заполнен диапазон Е23:К23 (диапазон L23:N23 можно очистить).

В столбец (диапазон С25:С29) также надо ввести новые значения коэффициентов целевой функции.

Затем изменим критериальное ограничение, которое записано в 30-й строке. Для этого в D30 введем формулу =СУММПРОИЗВ($C25:$C29;D25:D29)-D23, которую затем скопируем по строке на диапазон ячеек Е30:К30. Диапазон L30:N30 также можно очистить (как, впрочем, и диапазон L25:N29, поскольку искусственные переменные в дальнейшем не понадобятся).

 

Таблица 18 – Переход к решению исходной задачи

  A B C D E F G H I J K L M N
                             
  N xб cб B x1 x2 x3 x4 x5 x6 x7 у1 у2 у3
    x1   0,009       -0,003       0,003   -0,009
    x5                   -1    
    x2   0,487       -0,003   -0,025   0,003 0,025 -1,261
    x7   14,93       -0,411   -1,51   0,411 1,51 -156,337
    х3   0,004       0,006   0,025   -0,006 -0,025 2,269
  m+1     7,505       0,179   0,629        

 

Исходная задача здесь также поставлена на минимум. После пересчета критериальной строки становится ясно, что оптимальный план на минимум еще не найден, так как критерий оптимальности нарушается. Положительные коэффициенты 0,179 и 0,629 стоят в обоих небазисных столбцах. Для преобразования таблицы в базис можно ввести либо x4, либо x6. Если провести предварительные расчеты и определить, на сколько уменьшится стоимость рациона в том и другом случае (см. раздел 3.3), станет понятно, что выгоднее ввести в базис x4.

Новая таблица примет вид таблицы 19.

 

Таблица 19 – Оптимальная симплексная таблица

N xб cб B x1 x2 x3 x4 x5 x6 x7
  x1   0,011     0,444     0,011  
  x5   1,289     -167,444     -4,211  
  х2   0,489     0,556     -0,011  
  x7   15,222     68,889     0,222  
  х4   0,711     167,444     4,211  
m+1     7,378     -29,889     -0,122  

 

Теперь положительных коэффициентов в критериальной строке нет, т.е. задача на минимум решена. Оптимальный план Х* = (0,011; 0,489; 0; 0,711; 1,289; 0; 15,222), оптимум 7,378.

Ответ: для составления рациона следует взять 489 г зерна и 11 г известняка. При этом содержание кальция в рационе будет примерно на 0,7 г выше и на 1,3 г ниже, чем предельно допустимые нормы (от 4 до 6 г). Содержание белка будет на границе нормы, а содержание клетчатки – на 15,2 г меньше максимально допустимого. Стоимость рациона составит 7 руб. 38 коп. в неделю на одного цыпленка.

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

1 В чем заключается общая идея методов искусственного базиса?

2 Что можно сказать о разрешимости расширенной задачи?

3 Какова связь между решениями исходной и расширенной задач?

4 В чем сущность двухэтапного симплекс-метода?

5 Решите с помощью Microsoft Excel следующие задачи:

Пример 1. max х1 - 2х2 1 - 3х2 - х3 + х4 + х5 = 6 х1 + 4х2 + х3 + х5 = 15 2х1 - 4х2 - х3 + х4 = -3 х1-5 ³ 0 Пример 2. min 2х1 + 3х2 - 4х3 + х5 1 - 3х2 - х3 + х4 + х5 = -10 х1 + 4х2 + х3 + х5 = 15 2х1 - 4х2 - х3 + х4 = -3 х1-5 ³ 0

 


* Обратите внимание, что целевая функция задачи (19) всегда минимизируется, независимо от того, на минимум или на максимум поставлена задача (18).

* Не следует путать искусственные переменные с дополнительными, рассмотренными в разделе 1.4.1. Дополнительные переменные не нужно сводить к нулю: в допустимом плане они вполне могут быть положительными, так как вводились в неравенства (левая и правая части неравенства могут отличаться друг от друга на величину дополнительной переменной).

* Результат точных расчетов по этим формулам не совпадает, например, со значением D4, приведенным в таблице 17: D4 = 4*(-0,003) + 15*(-0,003) + 40*0,006 – 0 = 0,183 ¹ 0,179. Такое расхождение объясняется погрешностью округлений. Однако правильным для задачи в целом является именно результат 0,179, полученный с помощью электронной таблицы.

Не смотря на то, что Microsoft Excel показывает пользователю лишь первые значащие цифры числа, при осуществлении расчетов эта программа оперирует числами со значительно большей точностью. Поэтому проверка расчетов «по последней цифре» в данном случае неприменима. В самом деле, если ввести в две ячейки число 2,4, а в третью – формулу для суммы этих двух ячеек, то сумма составит 4,8. Если при этом потребовать, чтобы числа были показаны с точностью до нуля знаков после запятой, то в первых двух ячейках будет стоять число 2, а в третьей число 5. Казалось бы, результат «2 + 2 = 5» нелеп. Однако на самом деле электронная таблица «помнит» о том, что число 2,4 лишь приближенно равно двум. Поэтому сумма на самом деле приближенно равна 5, а не 4. И в ячейке с суммой на самом деле содержится не число 5, а число 4,8.


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



1 | 2 |



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