- •Общая характеристика симплекс-метода
- •Требования к задачам, решаемым симплекс-методом
- •В качестве критерия оптимизации может выступать:
- •Математическая формулировка задачи
- •Информационное обеспечение моделирования
- •Требования, предъявляемые к информации
- •Подготовка исходных данных для составления матрицы эмм и решения задачи на эвм
- •Технолого-экономические коэффициенты
- •Классификация технико-экономических коэффициентов
- •Моделирование системных ограничений. Формирование ограничений по земельным ресурсам
- •Моделирование использования сельскохозяйственных угодий с учетом трансформации
- •Моделирование использования пашни и сельскохозяйственных угодий с учетом структуры угодий или посевных площадей
- •По потребности в семенах и их производству
- •К ресурсным ограничениям относятся условия по использованию трудовых, денежно-материальных средств, минеральных удобрений, машин и механизмов, оросительной воде. Общий вид
- •2. Приведение задач линейного программирования к каноническому представлению
- •Алгоритм симплекс-метода
- •Расчет всех элементов новой симплекс-таблицы
- •К основным блокам информации относятся
- •Дополнительные переменные, попавшие в базис
- •Дополнительные переменные, не попавшие в базис
- •Введение в план дополнительной переменной
- •Двойственные задачи линейного программирования
- •Тогда структурный вид двойственной задачи будет иметь вид:
- •Изменение коэффициентов в целевой функции при переменной, вошедшей в базисное решение.
- •Изменение коэффициента в целевой функции при переменной, не вошедшей в базисное решение
- •Математическая модель задачи дз
- •Последняя симплекс-таблица задачи дз-1
- •Пределы устойчивости оптимального решения при изменении коэффициентов целевой функции
Последняя симплекс-таблица задачи дз-1
(Zmах = 240538)
Номер строки |
Базисные переменные |
Ном.огр. для допол. переменных |
Аiо |
|
||||
Х5 основн. |
Х7 изб.огр.8 |
Х8 ост.огр.1 |
Х11 ост.огр4 |
Х12 ост. огр.5 |
||||
1 |
Х13(ОСТ.) |
6 |
3812 |
6.2 |
0.2831 |
6.15 |
-0.877 |
2.338 |
2 |
Х9 (ОСТ.) |
2 |
28620 |
390 |
0.86 |
-10 |
0.2 |
-2.2 |
3 |
Х10 (ОСТ) |
3 |
9121 |
410 |
1.478 |
-89.62 |
4.208 |
-12.05 |
4 |
Х2 (ОСН.) |
- |
60.15 |
0.08 |
0.0035 |
0.0769 |
0.0015 |
0.029 |
5 |
Х3 (ОСН.) |
- |
30 |
0 |
-0.01 |
0 |
0 |
0 |
6 |
Х1 (ОСН.) |
- |
841.8 |
0.92 |
0.0025 |
0.923 |
0.0184 |
-0.049 |
7 |
Х14(ОСТ.) |
7 |
49.85 |
-0.08 |
-0.0035 |
-0.077 |
-0.015 |
-0.029 |
8 |
Х6 (ОСН.) |
- |
80880 |
-410 |
-1.478 |
89.62 |
-4.2 |
12.05 |
9 |
Х4 (ОСН.) |
- |
53.15 |
0.08 |
-0.0025 |
0.0769 |
-0.018 |
0.049 |
Анализ табл. 4 показывает, что производство сахарной свеклы хозяйству невыгодно: соответствующая основная переменная Х5 не попала в базис, то есть в оптимальном сочетании отраслей площадь пашни под свеклу X5 = 0.
Предположим теперь, что по каким-либо внешним причинам хозяйство вынуждено выращивать сахарную свеклу в определенном объеме. Пусть такой заданный объем равен 4800 га. Такое задание эквивалентно ограничению:
240*Х5 = 4800 ,
что соответствует требованию отвести под свеклу
Х5 = 4800/240 = 20 га.
Возникает естественный вопрос можно ли, не проводя полного решения задачи симплекс-методом с учетом указанного ограничения, а только используя уже полученное решение (табл.4), найти, тем не менее, новое оптимальное решение.
Введение в базис основной переменной, не вошедшей в него в последней симплекс-таблице, приведет к уменьшению значения целевой функции (задача на максимизацию).
Количественно уменьшение будет определяться значением элемента индексной строки, соответствующего вводимой переменной. Для рассматриваемой задачи это означает, что отведение под свеклу каждого гектара пашни будет приводить к уменьшению чистого дохода хозяйства (целевой функции) на 149 тыс.руб. Следовательно, при отведении под свеклу 20 га целевая функция примет значение:
Z = Zmax - 149*20 = 240538 - 2980 = 237558 тыс.руб.
Таким образом, информация, сосредоточенная в индексной строке, позволяет рассчитать новое значение целевой функции. Поскольку с экономической точки зрения проведенное изменение oптимального решения означает выделение определенного ресурса для производства свеклы, то это определенным образом может сказаться на других отраслях хозяйства, а также на объеме недоиспользованных ранее ресурсов, т.е. (с математической точки зрения) на значениях базисных переменных. Изменение самого решения, то есть значений базисных переменных, можно оценить с помощью коэффициентов замещения, стоящих в столбце вводимой в план переменной.
Новое значение любой базисной переменной из последней симплекс-таблицы (табл.4) может быть определено по формуле:
Х'jб = Xjб - К*Х5
где К - коэффициент замещения, стоящий на пересечении столбца, соответствующего вводимой переменной, и строки, соответствующей рассматриваемой базисной переменной. Таким образом, если коэффициент замещения К, соответствующий данной базисной переменной, отрицателен, то при введении в план переменной Х5 базисная переменная будет увеличиваться, в противном случае - уменьшаться. Указанное обстоятельство позволяет установить правило определения допустимого предела увеличения переменной X5: увеличение значения переменной Х5 не должно приводить к отрицательным значениям базисных переменных. Иначе говоря, если рассмотреть базисную переменную, которой соответствует положительный коэффициент замещения, и новое значение базисной переменной Х'jб принять равным нулю, то есть минимальному допустимому значению, то на последнего соотношения получим следующее ограничение на величину Х5:
Хjб – К*Х5 = 0
или
Х5max = Хjб/К.
Резюмируя изложенное, придем к следующему алгоритму введения в план основной переменной, не вошедшей в оптимальный базис:
1.разделить значения базисных переменных из последней симплекс-таблицы на соответствующие положительные коэффициенты замещения. Для рассматриваемого примера (табл. 4) получим (числа округленные):
Х13: 3810/6.2 = 614
X9: 28600/390 = 73
X10: 9120/410 = 22
X2: 60.2/0.08 = 753
X1: 842/0.92 = 915
X14: 581/0.08 = 7262
Сравнить все полученные числа и выбрать наименьшее. Это и будет максимально возможное значение вводимой переменной. Для рассматриваемого примера такое значение равно Х5max = 22.
Для понимания влияния полученного ограничения на допустимое значение вводимой переменной необходимо учесть следующее. Именно не превышение ограничения Х5max = 22 обеспечивает возможность построения нового оптимального плана из старого (табл.4) без полного решения задачи симплекс-методом. В принципе допустимо задать требование, нарушающее указанное ограничение, например, положить Х5 = 30, но в этом случае придется решать новую задачу (с ограничением Х5 = 30) самостоятельно.
(Терминология: иногда отношение, породившее ограничение на допустимые значения вводимой переменной, называет "узким местом").
2.Расчет новых значений целевой функции и базисных переменных (далее положим, что необходимо ввести в план Х5 = 20):
Z’ = Zmax – (Z5 – C5)* X5 (12)
Х’jб = Kjб – К*Х5 (13)
где (Z5-C5) - элемент индексной строки, соответствующей переменной Х5;
К - коэффициент замещения, расположенный в столбце, соответствующем небазисной переменной Х5, и в строке, соответствующей базисной переменной Xj6.
Результаты расчета заносят в таблицу следующего вида:
Таблица 5
Введение в оптимальный план задачи ДЗ-1 основной переменной Х5 = 20
Базисные Переменные |
Значение базисной переменной в оптимальном плане |
Коэффициенты замещения при Х5 |
Произведение к-в замещения на вводимую переменную (К*20) |
Новый план при Х5 = 20 |
Х13 (ост) |
3810 |
6.2 |
124 |
3686 |
Х9 (ост) |
2860 |
390 |
7800 |
20800 |
Х10 (ост) |
9120 |
410 |
8200 |
920 |
Х2 (осн) |
60 |
0.08 |
1.6 |
58.4 |
Х3 (осн) |
30 |
0 |
0 |
30 |
Х1 (осн) |
842 |
0.92 |
18.4 |
823.6 |
Х14 (ост) |
49.8 |
-0.08 |
-1.6 |
51.4 |
Х6 (осн) |
80900 |
-410 |
-8200 |
89100 |
Х4 (осн) |
58 |
0.08 |
1.6 |
56.4 |
Х5 |
0 |
1 |
-20 |
20 |
Z |
241000 |
149 |
2980 |
23758 |
Подчеркнем, что новый план также можно рассматривать как оптимальный, но для задачи ДЗ-1 с дополнительным условием Х5 = 20га.
Сравнивая таблицы 4 и 5, придем к выводу, что требование занять 20 гектаров пашни под сахарную свеклу привело к уменьшению:
- площади зерновых культур (X1) - на 18.4 гa;
- площади кормовых культур (Х4) - на 1.6 гa;
- поголовья коров (Х2) - на 2 головы.
При этом поголовье свиней сохранялось, что и должно быть, поскольку в условиях задачи сформулировано жесткое плановое задание по свинине. Общие денежные расходы хозяйства (X5) увеличились на 82тыс. руб., а чистый доход хозяйства (целевая функция) уменьшилась на 2980 тыс.руб.
Введение в план дополнительной переменной.
Введение в план отрицательного значения остаточной переменной эквивалентно увеличению соответствующего ресурса и приводит (в задачах на максимум) к возрастанию целевой функция. Соответствует предельное по модулю отрицательное значение переменной, при не превышении которого состав базисных переменных не меняется.
Введение в план положительного значения остаточной переменной эквивалентно уменьшению соответствующего ресурса со всеми вытекающими последствиями.
Аналитический алгоритм преобразования оптимального плана при введении в план остаточной переменной заключается в следующем (на примере задачи ДЗ - см. Ta6n.3):
1.Разделить значения базисных переменных из последней симплекс-таблицы на коэффициенты замещения, соответствующей вводимой в план остаточной переменной (причем, в отличие от случая введения в план основной переменной, необходимо делить как на положительные, так и на отрицательные коэффициенты).
Для рассматриваемого примера получим (числа округлены):
X13 : 3812/6.15 = 620
X9 : 28620/-10 = -2862
X10 : 9121/-89.62 = -102
X2 : 60.15/0.0769 = 782
X3 : 30/0 = --
X1 : 841.8/0.923 = 912
X14 : 49.85/-0.077 = -647
X6 : 80880/89.62 = 902
X4 : 53.15/0.0769 = 691
Из найденных чисел выбрать наименьшее положительное и наименьшее (по модулю) отрицательное. Тем самым будут установлены пределы, в которых можно изменять значение остаточной переменной без изменения структуры оптимального решения. Для рассматриваемого примера допустимые пределы изменения переменной Х8 составляют (-102, +620).
Рассмотрим далее случай введения в план переменной Х8 = -100.
2. Расчет новых значений целевой функции и базисных переменных по формулам аналогичным (12) и (13):
Z’ = Zmax – (Z8 – C8)*X8;
X’jб = Xjб – K*X8,
где (Z8-C8) – элемент индексной строки, соответствующий переменной Х38.|
К - коэффициент замещения, расположенный в столбце, соответствующем небазисной переменной Х8, и в строке, соответствующей базисной переменной Хjб.
Результаты расчетов заносят в таблицу 6.
Таблица 6
Введение в оптимальный план задачи ДЗ-1 остаточной переменной Х8 = -100
Базисные Переменные |
Значение базисной переменной в оптимальном плане |
Коэффициенты замещения при Х8 |
Произведение к-в замещения на вводимую переменную (К* (-100)) |
Новый план при Х8 = -100 |
Х13 (ост) |
3812 |
6.15 |
-615 |
4427 |
Х9 (ост) |
28620 |
-10 |
1000 |
27620 |
Х10 (ост) |
9121 |
-89.62 |
8962 |
159 |
Х2 (осн) |
60.15 |
0.0769 |
-7.69 |
67.84 |
Х3 (осн) |
30 |
0 |
0 |
30 |
Х1 (осн) |
841.8 |
0.923 |
-92.3 |
934.1 |
Х14 (ост) |
49.85 |
-0.077 |
7.7 |
42.15 |
Х6 (осн) |
80880 |
89.62 |
-8962 |
89842 |
Х4 (осн) |
53.15 |
0.0769 |
-7.69 |
60.84 |
Х8 |
0 |
-1 |
100 |
-100 |
Z |
240538 |
269.2 |
-26920 |
267458 |
Таким образом, в рассмотренном пример введение в план остаточной переменной Х8 =-100, то есть увеличение ресурса пашни на 100 гектаров, приводит к увеличению площади посевов зерна до 934 гектаров (X1) и к увеличению поголовья коров - 50 голов (Х2).
В случае введения в план избыточной переменной алгоритм вычислений идентичен рассмотренному.