Автор работы: Пользователь скрыл имя, 05 Февраля 2014 в 14:16, контрольная работа
Задание 2 Для производства 4-х видов продукции используется 3 вида сырья. Нормы расхода сырья (кг) запасы (кг) его ценность от реализации единицы продукции заданы таблицей. Составим план выпуска продукции, обеспечивающий получение максимальной прибыли, используя симплексный метод, а также построить двойственную задачу и решить ее симплекс-методом.
Задание 4 Решить задачу целочисленного программирования геометрическим методом.
Задание 1
Решить графически: maxF=2x1+4x2;
3x1+6x2≤12
2x1-x2≥-2
-x1+3x2≥0
x1≥0, x2≥0
Решение:
Строим прямые по уравнениям:
3x1+6x2=12
2x1-x2=-2
-x1+3x2=0
x1=0
x2=0
x2
2x1-x2=-2
4
C
3
2 A 3x1+6x2=12
1
B
x1
1 2 3 4
Целевая функция принимает максимальное значение в любой точке отрезка АВ.
Fmax=8
Задание 2
Для производства 4-х видов продукции используется 3 вида сырья. Нормы расхода сырья (кг) запасы (кг) его ценность от реализации единицы продукции заданы таблицей.
Составим план выпуска продукции, обеспечивающий получение максимальной прибыли, используя симплексный метод, а также построить двойственную задачу и решить ее симплекс-методом.
Нормы расхода ресурсов на единичное изделие  | 
  Запас ресурсов  | ||||
Изделие 1  | 
  Изделие 2  | 
  Изделие 3  | 
  Изделие 4  | ||
Ресурс 1  | 
  4  | 
  6  | 
  8  | 
  10  | 
  72  | 
Ресурс 2  | 
  10  | 
  3  | 
  5  | 
  9  | 
  90  | 
Ресурс 3  | 
  1  | 
  5  | 
  7  | 
  2  | 
  50  | 
Ценность  | 
  5  | 
  7  | 
  3  | 
  8  | 
  |
Решение:
Fmax=5x1+7x2+3x3+8x4
ограничения:
4x1+6x2+8x3+10x4≤72
10x1+3x2+5x3+9x4≤90
x1+5x2+7x3+2x4≤50
Приведем систему неравенств к 
системе уравнений путем 
4x1+6x2+8x3+10x4+x5=72
10x1+3x2+5x3+9x4+x6=90
x1+5x2+7x3+2x4+x7=50
Запишем преобразованную систему уравнений в векторной форме:
x1*P1+x2*P2+x3*P3+x4*P4+x5*P5+
P1= ; P2=; P3= ; P4= ; P5= ; P6= ; P7= ; P0=
Опорный план: X=(0;0;0;72;90;50)
Таблица I итерации:
Базис  | 
  P0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x5  | 
  72  | 
  4  | 
  6  | 
  8  | 
  10  | 
  1  | 
  0  | 
  0  | 
x6  | 
  90  | 
  10  | 
  3  | 
  5  | 
  9  | 
  0  | 
  1  | 
  0  | 
x7  | 
  50  | 
  1  | 
  5  | 
  7  | 
  2  | 
  0  | 
  0  | 
  1  | 
zj-cj  | 
  0  | 
  -5  | 
  -7  | 
  -3  | 
  -8  | 
  0  | 
  0  | 
  0  | 
Опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.
В качестве ведущего выберем столбец X4, так как это больший коэффициент по модулю.
Определим вектор подлежащий исключению из базиса:
Ѳ0=min(72/10; 90/9; 50/2;)=72/10=7.2
Следовательно, первая строка является ведущей. Разрешающий элемент равен 10.
Получаем новую симплекс таблицу:
Итерация II
Базис  | 
  P0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   X4  | 
  7,2  | 
  0,4  | 
  0,6  | 
  0,8  | 
  1  | 
  0,1  | 
  0  | 
  0  | 
x6  | 
  25,5  | 
  6,4  | 
  -2,4  | 
  -2,2  | 
  0  | 
  -0,9  | 
  1  | 
  0  | 
x7  | 
  35,6  | 
  0,2  | 
  3,8  | 
  5,4  | 
  0  | 
  -0,2  | 
  0  | 
  1  | 
zj-cj  | 
  57,6  | 
  -1,8  | 
  -2,2  | 
  3,4  | 
  0  | 
  0,8  | 
  0  | 
  0  | 
Текущий опорный план неоптимален.
Ведущий столбец – x2
Ѳ0=min(7.2/0.6; 25.5/-2.4; 35.6/3.8)=35.6/3.8=9
Третья строка является ведущей. Разрешающий элемент - 3,8.
Итерация III
Базис  | 
  P0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x4  | 
  1  | 
  0  | 
  -  | 
  1  | 
  0  | 
  -  | ||
x6  | 
  47  | 
  6  | 
  0  | 
  1  | 
  0  | 
  -1  | 
  1  | 
  |
x2  | 
  9  | 
  1  | 
  1  | 
  0  | 
  -  | 
  0  | 
  ||
zj-cj  | 
  78  | 
  -1,  | 
  0  | 
  6  | 
  0  | 
  0  | 
  
Текущий опорный план неоптимален.
Ведущий столбец – x1
Ѳ0=min=1 / = 4
Первая строка является ведущей. Разрешающий элемент -
Итерация IV
Базис  | 
  P0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x1  | 
  4  | 
  1  | 
  0  | 
  -  | 
  2  | 
  0  | 
  -  | |
x6  | 
  19  | 
  0  | 
  0  | 
  2  | 
  -17  | 
  -3  | 
  1  | 
  3  | 
x2  | 
  9  | 
  0  | 
  1  | 
  1  | 
  -  | 
  -  | 
  0  | 
  |
zj-cj  | 
  85  | 
  0  | 
  0  | 
  6  | 
  4  | 
  1  | 
  0  | 
  -  | 
Текущий опорный план неоптимален.
Ведущий столбец – x7
Ѳ0=min=19/3=5
Разрешающий элемент равен 3
Базис  | 
  P0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x1  | 
  6  | 
  1  | 
  0  | 
  -  | 
  0  | |||
x7  | 
  5  | 
  0  | 
  0  | 
  -5  | 
  -  | 
  1  | ||
x2  | 
  7  | 
  0  | 
  1  | 
  1  | 
  1  | 
  -  | 
  0  | |
zj-cj  | 
  86  | 
  0  | 
  0  | 
  6  | 
  3  | 
  1  | 
  0  | 
Индексная строка не содержит отрицательных элементов. Найден оптимальный план.
x1=6
x2=7
Fmax=5*6+7*7=86
Двойственная задача:
Составим двойственную задачу к прямой задаче:
4y1+10y2+y3≥5
6y1+3y2+5y3≥7
8y1+5y2+7y3≥3
10y1+9y2+2y3≥8
72y1+90y2+50y3→min
y1≥0; y2≥0; y3≥0
Приведем систему ограничений к системе неравенств смысла ≤ , умножив соответствующие строки на (-1)
-4x1-10x2-x3≤-5
-6x1-3x2-5x3≤-7
-8x1-5x2-7x3≤-3
-10x1-9x2-2x3≤-8
F(min)=72x1+90x2+50x3
Приведем систему неравенств к 
системе уравнений путем 
-4x1-10x2-x3+x4=-5
-6x1-3x2-5x3+x5=-7
-8x1-5x2-7x3+x6=-3
-10x1-9x2-2x3+x7=-8
Матрица коэффициентов этой системы имеет вид:
-4  | 
  -10  | 
  -1  | 
  1  | 
  0  | 
  0  | 
  0  | 
-6  | 
  -3  | 
  -5  | 
  0  | 
  1  | 
  0  | 
  0  | 
-8  | 
  -5  | 
  -7  | 
  0  | 
  0  | 
  1  | 
  0  | 
-10  | 
  -9  | 
  -2  | 
  0  | 
  0  | 
  0  | 
  1  | 
Полагая что свободные переменные равны 0, получим первый опорный план:
X=(0,0,0,-5,-7,-3,-8)
Базис  | 
  Р0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x4  | 
  -5  | 
  -4  | 
  -10  | 
  -1  | 
  1  | 
  0  | 
  0  | 
  0  | 
x5  | 
  -7  | 
  -6  | 
  -3  | 
  -5  | 
  0  | 
  1  | 
  0  | 
  0  | 
x6  | 
  -3  | 
  -8  | 
  -5  | 
  -7  | 
  0  | 
  0  | 
  1  | 
  0  | 
x7  | 
  -8  | 
  -10  | 
  -9  | 
  -2  | 
  0  | 
  0  | 
  0  | 
  1  | 
F(X)  | 
  0  | 
  -72  | 
  -90  | 
  -50  | 
  0  | 
  0  | 
  0  | 
  0  | 
План в симплексной таблице является псевдопланом.
Ведущая строка – x7
Ведущий столбец – x1
Разрешающий элемент: -10
Выполняем преобразования симплексной таблицы методом Жордано-Гаусса.
Итерация I
Базис  | 
  Р0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x4  | 
  -1  | 
  0  | 
  -6  | 
  -  | 
  1  | 
  0  | 
  0  | 
  -  | 
x5  | 
  -2  | 
  0  | 
  2  | 
  -3  | 
  0  | 
  1  | 
  0  | 
  -  | 
x6  | 
  3  | 
  0  | 
  2  | 
  -5  | 
  0  | 
  0  | 
  1  | 
  -  | 
X1  | 
  1  | 
  0  | 
  0  | 
  0  | 
  -  | |||
F(X)  | 
  57  | 
  0  | 
  -25  | 
  -35  | 
  0  | 
  0  | 
  0  | 
  -7  | 
План в симплексной таблице является псевдопланом, так как в базисном столбце имеются отрицательные элементы.
Ведущая строка – x5
Ведущий столбец – x3
Разрешающий элемент: -3
Итерация II
Базис  | 
  Р0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x4  | 
  -1  | 
  0  | 
  -6  | 
  0  | 
  1  | 
  -  | 
  0  | 
  -  | 
x3  | 
  0  | 
  -  | 
  1  | 
  0  | 
  -  | 
  0  | 
  ||
x6  | 
  6  | 
  0  | 
  -1  | 
  0  | 
  0  | 
  -1  | 
  1  | 
  |
x1  | 
  1  | 
  1  | 
  0  | 
  0  | 
  0  | 
  -  | ||
F(X)  | 
  78  | 
  0  | 
  -47  | 
  0  | 
  0  | 
  -9  | 
  0  | 
  -1  | 
План является псевдопланом.
Ведущая строка – х4
Ведущий столбец – х7
Разрешающий элемент: -
Итерация III
Базис  | 
  Р0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x7  | 
  4  | 
  0  | 
  17  | 
  0  | 
  -2  | 
  0  | 
  1  | |
x3  | 
  -  | 
  0  | 
  -  | 
  1  | 
  -  | 
  0  | 
  0  | |
x6  | 
  6  | 
  0  | 
  -2  | 
  0  | 
  -1  | 
  1  | 
  0  | |
x1  | 
  1  | 
  1  | 
  3  | 
  0  | 
  -  | 
  0  | 
  0  | |
F(X)  | 
  85  | 
  0  | 
  -19  | 
  0  | 
  -4  | 
  -9  | 
  0  | 
  0  | 
План является псевдопланом.
Ведущая строка – х3
Ведущий столбец – х2
Разрешающий элемент: -
Итерация IV
Базис  | 
  Р0  | 
  x1  | 
  x2  | 
  x3  | 
  x4  | 
  x5  | 
  x6  | 
  x7  | 
| 
   x7  | 
  3  | 
  0  | 
  0  | 
  5  | 
  -  | 
  -1  | 
  0  | 
  1  | 
x2  | 
  0  | 
  1  | 
  -  | 
  -  | 
  0  | 
  0  | ||
x6  | 
  6  | 
  0  | 
  0  | 
  -  | 
  -  | 
  -1  | 
  1  | 
  0  | 
x1  | 
  1  | 
  1  | 
  0  | 
  -  | 
  0  | 
  0  | ||
F(X)  | 
  86  | 
  0  | 
  0  | 
  -5  | 
  -6  | 
  -7  | 
  0  | 
  0  | 
Информация о работе Контрольная работа по "Методам принятия оптимальных решений"