Выполнить самостоятельно.

В соответствии с индивидуальным заданием №1 решить задачу максимизации с использованием симплексных таблиц. Вариант задания выбирается по номеру зачетной книжки:

-предпоследняя цифра - № столбца

-последняя цифра – № строки

Пример: Симплексным методом решить задачу максимизации.

F(x)=5x1-x2→max

2x1-x2+x3≤3

3x1+2x2≤6

x1≥0; x2≥0

1. Переведем систему ограничений в канонический вид, введя выравнивающие (базисные) неизвестные – x3 и x4.

Задача принимает следующий вид:

2x1-x2+x3=3

3x1+2x2=6

x1=x2=x3=x4=0

F-5x1+x2=0

2. Заполняем первую симплексную таблицу. В ней x3, x4 – основные переменные (базис).

Таблица 1.

  Базис Свободный член Переменные Оценочные отношения
х1 х2   х3   х4  
x3 2 -1 3/2
x5
F -5  

 

 

← разрешающая строка

 

 

разрешающий столбец

 

3. Проверяем выполнение критерия на max – первый опорный план не оптимальный, т.к. в F коэффициент при х1<0

4. Выбираем наибольший по модулю отрицательный коэффициент F, который определяет разрешающий столбец.

5. Делим свободные члены на коэффициенты разрешающего столбца (определяем оценочные отношения). Выбираем разрешающую строку, где это отношение минимально

 

Min (3, 2)= 3

2 2

Разрешающий элемент будет а11=2

Для построения таблицы 2 в качестве основной переменной выбираем х1, т.к. она образует разрешающий столбец таблицы 1.

 

Таблицы 2.

  Базис Свободный член Переменные Оценочные отношения
х1 х2   х3   х4  
х1 3/2 -1/2 1/2 -3
х4 3/2 7/2 -3/2 3/7
F 15/2 -3/2 5/2  

 

 

← разрешающая строка

 

разрешающий столбец

Построение таблицы №2.

1. Заменим переменные в базисе с х3 на х1.

2. Делим элементы разрешающей строки (табл.1) на разрешающий элемент, результаты занесем в строку х1 в табл. №2.

3. В остальных клетках разрешающего столбца (табл. 1) записываем 0.

4. Остальные клетки заполняем по правилу прямоугольника

НЭ=СТЭ-(АхВ)/РЭ

 

СТЭ А

 


 

В РЭ

 

5. Критерий эффективности опять не выполнен, т.к. F имеет два коэффициента -3 < 0

Построение таблицы №3

Таблицы 3.

  Базис Свободный член Переменные Оценочные отношения
х1 х2   х3   х4  
х1 12/7 2/7 1/7  
х2 3/7 -1/7 2/7  
F 57/7 13/7 3/7  

 

Из таблицы 3 видно – критерий оптимальности выполнен.

Оптимальные базисные решения Х*=(12, 3,0,0)

7 7

F=5x1-x2=5x12- 3=57

7 7 7

 

 

ИНДИВИДУАЛЬНЫЕ ЗАДАНИЯ №1

  0,1,2 3,4,5,6 7,8,9
f(x)=3x1-2x2; 2x1+x2≤2; x1+2x2≤2; x1,x2≥0; f(x)=4x1+8x2; x1+x2≤3; x1+2x2≤2; x1,x2≥0; f(x)=-4x1+8x2; 3x1+5x2≤15; x1-x3≤1; x1,x2≥0;
f(x)=3x1-2x2; -x1+2x2≤2; 2x1-x2≤2; x1,x2≥0; f(x)=-x1+6x2; 4x1+3x2≤12; -x1+x2≤1; x1,x2≥0; f(x)=-x1+6x2; x1+x2≤3; -2x1+x2≤2; x1,x2≥0;
f(x)=x1+6x2; 3x1+4x2≤12; -x1+x2≤2; x1,x2≥0; f(x)=2x1+6x2; 3x1+4x2≤12; x1+2x2≤2; x1,x2≥0; f(x)=8x1+12x2; 2x1+x2≤4; 2x1+5x2≤10; x1,x2≥0;
f(x)=8x1+12x2; -3x1+2x2≤0; 4x1+3x2≤12; x1,x2≥0; f(x)=3x1-2x2; 2x1+x2≤2; 2x1+3x2≤6; x1,x2≥0; f(x)=6x1+4x2; x1+2x2≤2; -2x1+x2≤0; x1,x2≥0;
f(x)=6x1+4x2; 3x1+2x2≤6; 3x1+x2≤3; x1,x2≥0; f(x)=8x1+6x2; -x1+x2≤1; 3x1+2x2≤6; x1,x2≥0; f(x)=8x1+6x2; -x1+x2≤2; 3x1+4x2≤12; x1,x2≥0;

 

Пример. Решить следующую задачу максимизации

F(x)=5x1-x2→max

2x1-x2+x3≤3

3x1+2x2≤6

x1≥0; x2≥0

 

Запишем систему в каноническом виде

1) 2x1-x2+x3=3

2) 3x1+2x2=6

3) x1=0; x2=0

 

Решение. В системе координат Х1ОХ2 построим ОДР.

Из первого уравнения:

При Х1=0 X2=-3 → (X1,X2)=(0,-3) прямая (1)

При Х2=0 Х1=3 → (X1,X2)=(3, 0)

2 2

 

Из второго уравнения:

При Х1=0 X2=3 → (X1,X2)=(0,3) прямая (2)

При Х2=0 Х1=2 → (X1,X2)=(2, 0)

 

Из условия (3) (Х1, Х2)=(0,0) Х1=0 – прямая (3)

Х2=0 – прямая (4)

 

Для определения с какой стороны прямых ОДР все точки удовлетворяют неравенствам возьмем точку (1,1). Таким образом, ОДР лежит внутри АВСD.

Для определения точки выхода, строим прямую:

F(x)=5x1-x2=0

5x1=x2→ x1=x2=0 при x1=1, x2=5

Строим нулевую линию уровня целевой функции.

Строим нормальный вектор q (5.-1)

Передвигая линию уровня целевой функции параллельно себе, по вектору q, фиксируем ее крайнее положение (т.В), которая принадлежит прямым (1) и (2).


2x1-x2+x3=3

3x1+2x2=6

 

x1=12 x2=3, т.е. Х*=(12, 3)

7 7 7 7

F(Х*)=5х12 -3=57

7 7 7

 

 

2.8Двойственные задачи линейного программирования (ДЗЛП).

Рассмотрим пример 2.6 и пример 2.5 и представим их в виде таблице

 

Задача 1 (исходная) Задача 2 (двойственная)
F=2x1+3x2→max при ограничениях: x1+3x2≤18 2x1+x2≤16 x2≤5 3x1≤21 x1, x2≥0 Составить такой план выпуска продукции X=(x1, х2) при котором прибыль при реализации будет max. В общем виде задача запишется: n F=∑Cjxj→max, j=1 при ограничениях: n j=1,2…m ∑aijхi≤ bi i=1,2…n j=1 Z=18 y1+16y2+5y3+21y4→min при ограничениях: y1+2y2+3y4≥2 3y1+y2+y3≥3 y1, y2, y3, y4≥0 Найти такой набор цен ресурсов Y=(y1,y2,y3,y4), при котором общие затраты на ресурсы будут минимальны. В общем виде задачу запишем: m Z=∑ biyi →min, i=1 при ограничениях: m j=1,2…m ∑ aijyi≤ cj i=1,2…n i=1  

 

Из таблице видно, что задачи 1 и 2 обладают следующими свойствами:

1) В задаче 1 ищут max, во 2 ищут min.

2) Коэффициенты целевой ф-ции 1-ой задачи являются свободными членами системы ограничения 2-ой задачи.

3) Каждая из задач задана в стандартной форме, причем в задаче максимизации все неравенства вида ≤, а в задаче минимизации вида ≥.

4) Матрицы коэффициентов в системе ограничения обеих задач являются транспонированными друг другу.

Для 1 3 18

задачи 1 A1= 2 1 16

0 1 5

3 0 21

2 3 F

 

 

Для

задачи 2 A’1= 1 2 0 3 2

3 1 1 0 21

18 16 5 21 Z

 

 

 

 

5) Число неравенств в системе ограничений 1-ой задачи совпадает с числом переменных 2-ой задачи.

6) условие не отрицательности переменных есть в обеих задачах.

Две задачи, обладающие указанными свойствами, называются симметричными взаимно – двойственными задачами.

2.9Алгоритм составления двойственных ЗЛП:

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

1) Привести все неравенства системы ограничений исходной задачи к одному смыслу: если в исходной задаче ищут max целевой функции, то все неравенства системы ограничений привести к виду ≤, а если ищут min – к виду ≥. Для этого неравенства, в которых это требование не выполняется, умножить на (-1).

2) Составить расширенную матрицу системы A1.

3) Найти матрицу A1, транспонированную к матрице A1.

4) Сформулировать двойственную ЗЛП

Этот алгоритм можно проиллюстрировать следующим примером:

Составить задачу двойственную следующей задаче:

F=14x1+10x2+14x3+11x4→max

Система ограничений:

1+2х2+2х3+3х4≤35

х12+2х3+3х4≤30

12+2х34≥40*

х1≥0

х2≥0

х3≥0

х4≥0

 

Решение:

Алгоритм Конкретные ситуации данного алгoритма
1.Все неравенства системы ограничений привести к виду ≤;   2.Составить расширенную матрицу системы А1, в которую включить: - матрицу коэффициентов переменных системы ограничений. - столбец свободных членов. - строку коэффициентов при переменных линейной функции.   3.Найти матрицу A'1 транспонированную к матрице А1.   4.Сформулировать двойственную задачу на основании полученной матрицы A'1 Приведем третье неравенство к виду: -3х12-2х34≤-40    
       
   
 
 


4 2 2 3 35

А1= 1 1 2 3 30

-3 -1 -2 -1 -40

14 10 14 11 F

 

4 1 -3 14

2 1 -1 10

A'1 = 2 2 -2 14

3 3 -1 11

35 30 -40 Z

 

Z=35y1+30y2-40y3→min


4y1+y2-3y3≥14

2y1+y2-y3≥10

2y1+2y2-2y3≥14

3y1+3y2-y3≥11

y1≥0; y2≥0; y3≥0

 

Выполнить самостоятельно:

1. Составить ДЗЛП следующей задачи

F=-x1+2x2→max

При ограничении: 2x1-x2≥1

-x1+4x2≤24

x1 -x2≤3

x1 +x2≥5

2. Составить ДЗЛП, используя исходные данные таблицы №1. Вариант задания выбирается по номеру зачетной книжки

–предпоследняя цифра - № строки

–последняя цифра - № столбца.

3.Теорема двойственности

Первая теорема двойственности

Устанавливает связь между оптимальными решениями ДЗЛП.

 

Если одна из взаимно свойственных задач имеет оптимальное решение, то его имеет и другие, причем оптимальные значения их линейных функций равны:

Fmax=Zmin или F(X*)=Z(Y*)

Если линейная функция одной из задач не ограничена (F=∞), то по условию другой задачи противоречивы.

Рассмотрим пример подтверждающий справедливость I теоремы.

Задача 1. Даны 2 взаимно двойственные задачи

I F=2x1-3x2→max При ограничениях: x1+3x2≤ 18 2x1+x2≤16 x2≤5 3x1≤21 x1≥0 x2≥0 II Z=18y1+16y2+5y3+21y4→min При ограничениях: y1+2y2+3y4 ≥2 3y1+y2+y3≥3   yi≥0 (i=1,2,3,4)

 

Задача I об использовании ресурсов и двойственная ей задача II были решены ранее (§ 2.6, 2.7) получены оптимумы, Fmax=24 и Zmin=24, т.е. заключение I теоремы двойственности верно.

 

Экономический смысл теоремы

План производства X*=(x1*,x2*…xn*) и набор цен Y*=(y1*,y2*…ym*) оказываются оптимальными только тогда, когда прибыль от продукции, найденная при “внешних” ценах c1,c2…cn, равны затратам на ресурсы по “внутренним ” (определяемым из решения задачи) ценам, для всех же других планов, прибыль от продукции всегда меньше (или равна) затрат на ресурсы, т.е. F(X)≤Z(Y).

 

Вторая теорема двойственности

Пусть даны две взаимно двойственные задачи (табл. 1). Преобразуем эти две задачи для решения симплексным методом. Тогда системы ограничений каждой из задач будут

 
 


 

aij xj + xn+i = bi j=1,2…n

 

aij yj - ym+j = cji=1,2…m

установим соответствие между первоначальными переменными одной из двойственных задач и дополнительными переменными другой задачи. См табл.1

Переменные исходной задачи I
первоначальные дополнительные
x1 x2 … xi … xn ↕ ↕ ↕ ↕ ym+1 ym+2 … ym+i …ym+n xn+1 xn+2 … xn+j … xn+m ↕ ↕ ↕ ↕ y1 y2 … yj … ym
дополнительные первоначальные
Переменные двойственной задачи II

Вторая теорема двойственности:

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

 

Справедливость I теоремы двойственности подтверждается решением задачи 1

На основании соответствий в табл. 1 установим следующие соотношения

x1 x2 x3 x4 x5 x6

↕ ↕ ↕ ↕ ↕ ↕

y1 y2 y3 y4 y5 y6

Обе задачи решены симплексным методом. На последнем шаге решения каждой задачи получено:

 

 

Компоненты оптимального решения двойственной задачи yi*=4/5, y2*=3/5, y3*=0, y4*=0, y5* =0, y6*=0 равны (по абсолютной величине) коэффициентом при соответствующих переменных линейной функции

F=24-4/5x3-3/5x4-0x5-0x6-0x1-0x2,

а компоненты оптимального решения исходной задачи xi*=6, x2*=4, x3*=0, x4*=0, x5* =1, x6*=3 равны коэффициентам при соответствующих переменных линейной функции II:

F=24+6y5+4y6+0y1+0y2+1y3+3y4

Компоненты оптимального решения двойственной задачи называются оптимальным оценкам исходной задачи.