Задача по линейному программированию

Автор работы: Пользователь скрыл имя, 14 Сентября 2011 в 09:57, задача

Краткое описание

решение одной задачи.

Содержимое работы - 1 файл

Задача 1.doc

— 110.00 Кб (Скачать файл)

  Программированные задания

Вариант № 7 

Задание 1 

Решить задачу линейного программирования симплекс-методом. 

f(x) = 4 x1 + 5x2 2 x3 + x4 + 2x5 → min
 
  -3x1 + 5 x2 + 4 x3 + 2 x4 + 2 x5 =   1   (1)
    -4x1 - 6 x2 -   x3 +   x4 + 3 x5 =   -1   (2)

   x1, x2, x3, x4, x≥ 0 

Решение:

 
I этап. Избавимся от отрицательных свободных членов, умножив на -1 ограничения 2.

  -3x1 + 5 x2 + 4 x3 + 2 x4 + 2 x5 =   1               (1)
     4x1 + 6 x2 +   x3 -   x4 - 3 x5 =   1               (2)

  x1, x2, x3, x4, x5 ≥ 0 

          В исходной задаче базисные переменные отсутствуют, это значит, что исходная задача не содержит в себе допустимого базисного решения. Для его нахождения составим расширенную задачу и будем рассматривать ее в симметричной форме (F=- f, тогда поиск min f эквивалентен поиску max F, причем min f = max F).

        Для этого в левые части уравнений системы ограничений введем дополнительные неотрицательные неизвестные x6 и х7 . Получим следующую систему ограничений:

 
 
  -3x1 + 5 x2 + 4 x3 + 2 x4 + 2 x5 + x6         =   1               (1)
     4x1 + 6 x2 +   x3 -   x4 - 3 x5         + x7 =   1               (2)

  x1, x2, x3, x4, x5, x6, x7 ≥ 0 
 
с базисными переменными x6, x7.
 

I I этап. Целью решения данной системы является получение допустимого базисного решения, не содержащего искусственных переменных (x6, x7).

 Приравняем свободные переменные системы уравнений (ограничений) к  нулю: х1= х2= х3= х4= х5=0, получаем начальное опорное решение расширенной задачи х1 = (0,0,0,0,0,1,1). Вычисляем оценки разложений векторов условий по базису опорного решения и записываем в симплексную таблицу (табл. 1). При этом оценки ∆ј и F(x1) для удобства  вычислений записываем в две строки: в первую - слагаемые ∆'ј, независящие от коэффициентов x6 и x7, во вторую - слагаемые ∆"ј, зависящие от коэффициентов x6 и x7.

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

      Составим первую симплекс-таблицу. 

                                                                                                                    Табл.1 
Симплекс-таблица № 1

cі БП -4 -5 -2 -1 -2 0 0 bі Пояснения
х1 х2 х3 х4 х5 х6 х7
0 х6 -3 5 4 2 2 1 0 1 max (׀1-׀, ׀-11׀, ׀-5׀, ׀-1׀, ׀1׀) = 11,
0 х7 4 6 1 -1 -3 0 1 1 х2 в базис. Второй столбец – ключевой.
  'ј 4 5 2 1 2 0 0 0 min (1/5, 1/6) = 1/6,
  ∆"ј -1 -11 -5 -1 1 0 0 -2 х7 из базиса. Вторая строка – ключевая.
 

Индексная строка находится по формулам:

 
 

 
 
 
 

      Проверяем критерий оптимальности: начальное  опорное решение не является оптимальным, так как в индексной строке ∆"ј имеются отрицательные оценки. 

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

IV этап. Построение нового базиса.

      Выбираем  номер вектора хj вводимого в базис опорного решения, и вектора хі , выводимого из базиса. Для этого вычисляем приращения целевой функции ∆Fј при введении в базис каждого из векторов с положительной оценкой и находим минимум этого приращения. При этом слагаемые оценок ∆'ј пренебрегаем до тех пор, пока хотя бы одно слагаемое ∆"ј отлично от нуля. 

min     bi      = min  (1/5; 1/6) = 1/6.

         ai2>0

V этап. Определение нового опорного плана.

 Ключевое число 6. Вектор х7 , выводимый из базиса, исключаем из рассмотрения. Получаем опорное решение х2 = (0, 1/6, 0, 0, 0, 1/6, 0) с базисом х2 , х6 . Составим вторую симплексную таблицу (табл.2) 

                                                                                                                    Табл.2 
Симплекс-таблица № 2

cі БП -4 -5 -2 -1 -2 0 bі Пояснения
х1 х2 х3 х4 х5 х6
5 х2 2/3 1 1/6 -1/6 -1/2 0 1/6 max = 9/2,
0 х6 -19/3 0 19/6 17/6 9/2 1 1/6 х5 в базис. Пятый столбец – ключевой.
  'ј 2/3 0 7/6 11/6 9/2 0 -5/6 min = 1/27,
  ∆"ј 19/3 0 -19/6 -17/6 -9/2 0 -1/6 х6 из базиса. Вторая строка – ключевая.
 

      Данное  решение не является оптимальным, так  как векторы х3, х4 и х5 имеют отрицательные оценки. Проделываем со второй симплексной таблицей то же, что и с первой. Получаем новую третью симплекс-таблицу (табл.3). 

                                                                                                                    Табл.3 
Симплекс-таблица № 3

cі БП -4 -5 -2 -1 -2 bі Пояснения
х1 х2 х3 х4 х5
2 х5 -38/27 0 19/27 17/27 1 1/27 max = 9/2,
5 х2 -1/27 1 14/27 4/27 0 5/27 х5 в базис. Пятый столбец – ключевой.
  'ј 7 0 -2 -1 0 -1 min = 1/19,
  ∆"ј 0 0 0 0 0 0 х6 из базиса. Вторая строка – ключевая.
 

      Получено  оптимальное решение расширенной задачи. Все искусственные переменные вышли из базиса и поэтому можно приступить к решению исходной задачи, приняв полученное базисное решение в качестве опорного. Строка ∆"ј больше не нужна, принятие решения о ключевом столбце, во всех последующих симплекс-таблицах, будем принимать по строке ∆'ј . Составим четвертую симплексную таблицу (табл.4).

Информация о работе Задача по линейному программированию