
- •Лабораторный практикум
- •По дисциплине «Математические методы в экономике»
- •Часть 2 .Линейное и дискретное программирование
- •Для подготовки студентов по направлениям 080500 - менеджмент и 080100 - экономика
- •Лабораторный практикум по дисциплине «Математические методы в экономике» . Часть 1.Линейное и дискретное программирование. Уч. Пособие. М.: фгоу впо ргау - мсха им. К.А. Тимирязева, - 127 с.
- •Содержание
- •Введение
- •Лабораторная работа № 13 «Решение задач линейного программирования в приложении ms Excel «Поиск решения»»
- •Последовательность выполнения
- •Лабораторная работа № 14 «Метод искусственного базиса или м- метод»
- •Лабораторная работа №15-17 «Основы теории двойственности»
- •1 Экономический смысл двойственных задач
- •2 Правило записи двойственных задач. Симметричные двойственные задачи
- •3 Свойства двойственных задач. Основные теоремы двойственности
- •4 Запись оптимального решение двойственной задачи по оптимальному решению прямой задачи. Краткий анализ оптимального решения прямой задачи
- •Лабораторная работа №18 «Транспортная задача на минимум целевой функции»
- •Лабораторная работа № 19 «Особенности решения транспортной задачи на максимум»
- •Лабораторная работа № 20 Видоизменения транспортной задачи (блокировка перевозок, ограничение пропускной способности).
- •I. Общая постановка задачи
- •II. Общая постановка задачи
- •Лабораторная работа № 21 «Транспортная задача с учетом производственных затрат»
- •Лабораторная работа № 22 «Решение задач о назначениях»
- •1 Постановка задачи и математическая формализация условий
- •2 Венгерский алгоритм решения задачи на минимум целевой функции
- •3 Особенности решения задачи на максимум
- •Лабораторная работа №23 «Целочисленное программирование»
- •Лабораторная работа № 24 «Параметрическое программирование с параметром в целевой функции»
- •Алгоритм решения задач в полных симплексных таблицах на максимум целевой функции
- •Лабораторная работа № 25. «Параметрическая транспортная задача с параметром в целевой функции»
- •Алгоритм решения параметрической транспортной задачи с параметром в целевой функции
- •Глоссарий основных понятий
- •Рекомендуемая литература
На пересечении разрешающей строки и разрешающего столбца определяется разрешающий элемент и выполняются преобразования однократного замещения (аналогично симплексному методу).
Далее алгоритм повторяется с пункта 2 до получения оптимального целочисленного решения. Алгоритм метода Гомори конечен.
Для решения задачи на минимум целевая функция умножается на "-1" и в дальнейшем используется алгоритм решения задачи на максимум. При записи оптимального решения значение целевой функции умножается на минус единицу, то есть переходим от максимума целевой функции к минимуму.
Пример
Найти maxZ=x1+3x2, при условиях
x1+x2 1
x1+2x2 4
-x1+2x2 2
x1 0; x2 0
x1, x2 - целые
Сначала необходимо найти оптимальный план, не принимая во внимание условие целочисленности. Для этого приводим задачу к канонической форме, исходное опорное решение получаем методом искусственного базиса (М-методом), заносим в симплексную таблицу.
maxZ=x1+3x2 +0x3 +0x4 +0x5 –My1
x1+x2 – x3+y1 =1
x1+2x2 +x4 = 4
-x1+2x2 +x5 = 2
xj 0, j=1 5, y1 0
№1 |
Cj |
0 |
1 |
3 |
0 |
0 |
0 |
-M |
Сi |
П Б |
ai0 |
x1 |
x2 |
x3 |
x4 |
x5 |
y1 |
-M |
y1 |
1 |
1 |
1 |
-1 |
0 |
0 |
1 |
0 |
x4 |
4 |
1 |
2 |
0 |
1 |
0 |
0 |
0 |
x5 |
2 |
-1 |
2 |
0 |
0 |
1 |
0 |
Z |
j |
-M |
-M-1 |
-M-3 |
M |
0 |
0 |
0 |
В четвертой симплексной таблице получено оптимальное решение задачи.
№4 |
Cj |
0 |
1 |
3 |
0 |
0 |
0 |
-M |
Сi |
П Б |
ai0 |
x1 |
x2 |
x3 |
x4 |
x5 |
y1 |
3 |
x2 |
3/2 |
0 |
1 |
0 |
1/4 |
1/4 |
0 |
1 |
x1 |
1 |
1 |
0 |
0 |
1/2 |
-1/2 |
0 |
0 |
x3 |
3/2 |
0 |
0 |
1 |
3/4 |
-1/4 |
-1 |
Z |
j |
11/2 |
0 |
0 |
0 |
5/4 |
1/4 |
M |
Так как искусственная переменная y1=0, то можно исключить ее из таблицы и в дальнейшем не рассматривать.
Переменная x2 и x3 имеют дробные значения, т.е. условия целочисленности не выполняются. Так как дробные части значений переменных x2 и x3 равны, то можно выбрать любую переменную. Введем ограничение на целочисленность переменной x2.
Определим коэффициенты дополнительного ограничения:
аi0 = 3/2 – 1 =1/2 ( 1 – максимальное целое число, не превышающее 3/2);
аi1 = аi2 = аi3 = 0 (0-0=0; 1-1=0)
аi4 = 1/4-0=1/4
аi5 =1/4-0=1/4
Получим дополнительное ограничение: 1/4х4+1/4х5≥1/2
Приводим это ограничение к каноническому виду:
1/4х4+1/4х5 – х6 =1/2;
умножаем на "-1": -1/4х4-1/4х5 + х6 = -1/2 и заносим в последнюю симплексную таблицу.
№4а |
Cj |
0 |
1 |
3 |
0 |
0 |
0 |
0 |
|
Сi |
П Б |
ai0 |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
|
3 |
x2 |
3/2 |
0 |
1 |
0 |
1/4 |
1/4 |
0 |
|
1 |
x1 |
1 |
1 |
0 |
0 |
1/2 |
-1/2 |
0 |
|
0 |
x3 |
3/2 |
0 |
0 |
1 |
3/4 |
-1/4 |
0 |
|
0 |
x6 |
-1/2 |
0 |
0 |
0 |
-1/4 |
-1/4 |
1 |
i* |
Z |
j |
11/2 |
0 |
0 |
0 |
5/4 |
1/4 |
0 |
|
ДСО
|
- |
- |
- |
-5 |
-1
|
- |
|
||
|
|
|
|
|
j* |
|
|
Принимаем новую строку за разрешающую и подсчитываем отношение неотрицательных оценок к строго отрицательным коэффициентам разрешающей строки. По наибольшему отрицательному отношению выбираем разрешающий столбец (это столбец x5). Проводим преобразования однократного замещения.
№5 |
Cj |
0 |
1 |
3 |
0 |
0 |
0 |
0 |
Сi |
П Б |
ai0 |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
3 |
x2 |
1 |
0 |
1 |
0 |
0 |
0 |
1 |
1 |
x1 |
2 |
1 |
0 |
0 |
1 |
0 |
-2 |
0 |
x3 |
2 |
0 |
0 |
1 |
1 |
0 |
-1 |
0 |
x5 |
2 |
0 |
0 |
0 |
1 |
1 |
1 |
Z |
j |
5 |
0 |
0 |
0 |
1 |
0 |
1 |
В таблице 5 получено оптимальное решение (все j 0) с целочисленными значениями переменных: Х*(2,1,2,0,2,0). MaxZ(X*)=5.
Общая постановка задачи
Найти оптимальное решение задачи симплексным методом без условий целочисленности. Проверить: является ли найденное решение целочисленным и если нет, то найти целочисленное решение методом Гомори.
Список индивидуальных данных
1. max Z= x1+ 2x2 1,5x1 + 2x2 3 2x1+ x2 3 x1+ x2 4 x1 0; x2 0 x1, x2-целые |
2. max Z= x1+x2 1,5x1 + 2x2 3 2x1+ x2 3 x1+ 1,5x2 4 x1 0; x2 0 x1, x2-целые |
3. max Z= 2x1+3x2 x1 + 2x2 3 -x1+ 2x2 2 x1 - x2 1 x1 0; x2 0 x1, x2-целые |
4. max Z= 2x1+x2 1,5x1 + 2x2 5 2x1+ x2 3 x1 0; x2 0 x1, x2-целые
|
5. max Z= x1+x2 3x1 + x2 5 x1+ x2 4 x2 4 x1 0; x2 0; x1, x2; -целые |
6. max Z= x1+x2 1,5x1 + 2x2 5 2x1+ x2 3 x1 + 1,5x2 3 x1 0; x2 0 x1, x2-целые |
7. max Z= 1,5x1+ 2x2 x1 + 2,5x2 6 x1+ x2 4 x1 0; x2 0 x1, x2-целые
|
8. min Z= x1+x2 4x1 + x2 2 x1+x2 6 x1 3 x1 0; x2 0 x1, x2-целые |
9. min Z= x1+x2 x1 + 4x2 2 x1+x2 6 x1 3 x1 0; x2 0; x1, x2 -целые |
10. min Z= x1-x2 x1 + 2x2 3 2x1- x2 2 x1 0; x2 0 x1, x2-целые |
11. max Z= 2x1-x2 x1 + x2 2 x1- x2 3 2x1 + x2 1 x1 0; x2 0 x1, x2-целые
|
12. max Z= 2x1-x2 x1 + x2 4 2x1- x2 1 x1 0; x2 0 x1, x2-целые
|
13. max Z= 2x1 x1 + 2x2 3 2x1+ x2 3 x1 0; x2 0 x1, x2-целые
|
14. max Z= 11x1+9x2 x1 + 3x2 3 3x1+ x2 9 x1 + 0,5x2 2 x1 0; x2 0 x1, x2-целые |
15. min Z= x1 + x2 x1+ 4x2 2 x1+ x2 6 x1 0; x2 0; x1, x2-целые
|
16. max Z= 2x1+x2 -x1 + 2x2 3 x1+ 2x2 7 x1 0; x2 0 x1, x2-целые |
17. min Z= x1 + x2 x1+ 2x2 3 2x1+ x2 3 x1 0; x2 0; x1, x2-целые |
18. max Z= 2x1 x1 + 2x2 4 2x1+ x2 3 x1 0; x2 0 x1, x2-целые |
Задача 19. Необходимо забороновать не менее 375 га зяби. На эту работу можно выделить не более 7 тракторов двух марок А и В. Производительность трактора марки А - 80 га/смена, В – 40 га/смена. Затраты на работу на тракторе марки А – 50 ден.ед/смена –В – 20 ден.ед./смена. Необходимо направить на работу такое количество тракторов, чтобы затраты были минимальными , а работа выполнена за одну смену.
Контрольные вопросы
Сформулируйте общую постановку задачи целочисленного программирования.
Чем вызвана необходимость введения на переменные условия целочисленности?
На чем основывается метод Гомори?
Как реализуется принцип Гомори в процессе решения задачи симплексным методом?
Как определяется строка симплексной таблицы, по которой вводится дополнительное ограничение?
Как рассчитываются коэффициенты нового ограничения?
Какой тип ( , , =) имеет новое ограничение?
Как занести новое ограничение в симплексную таблицу?
Как выбирается разрешающая строка?
Как выбирается разрешающий столбец?
Лабораторная работа № 24 «Параметрическое программирование с параметром в целевой функции»
Теоретическая часть
Общая задача линейного программирования имеет вид:
В ней коэффициенты и свободный член целевой функции (сj,c0), коэффициенты системы ограничений (аij), свободные члены (аi0), являются постоянными величинами. Однако на практике встречаются случаи, когда эти величины не являются постоянными, а их значения изменяются в некоторых интервалах. Кроме того, возникает необходимость не только получить конкретное оптимальное решение экономической задачи при фиксированных значениях сj, аij, аi0, но и исследовать зависимость оптимального решения от изменяющихся параметров, то есть установить, в каких допустимых пределах можно их изменять, чтобы план оставался оптимальным.
Исследование зависимости оптимального решения от изменяющихся параметров составляет предмет параметрического программирования.
Рассмотрим задачу линейного программирования с параметром в целевой функции.
Пусть
коэффициент целевой функции Сj
может изменяться в некоторых допустимых
пределах, т.е.
,
где
постоянные значения,
а t
– параметр, изменяющийся в некоторых
пределах
.
В этом случае задачу математически можно сформулировать следующим образом:
Дана линейная функция
Для
каждого значения t в интервале
,
где
произвольные действительные числа,
найти неотрицательный вектор
удовлетворяющий
системе ограничений и обеспечивающий
оптимальное значение целевой функции.
Методом параметрического программирования может решаться такая экономическая задача, как определение диапазона оптимального объема производства продукции при изменении условий (цен) реализации, при изменении себестоимости производства.
Алгоритм решения задач в полных симплексных таблицах на максимум целевой функции
1. Записать исходное опорное решение задачи в полную симплексную таблицу.
Форма симплексной таблицы несколько изменяется.
Таблица 1. Симплексная таблица для решения задач с параметром в целевой функции
№ табл. |
Cj' |
-C0 |
-C1' |
… |
-Cm' |
… |
-Cn' |
|
Cj" |
-C0 |
-C1" |
… |
-Cm" |
… |
-Cn" |
||
Ci' |
Ci" |
П БП |
аi0 |
X1 |
… |
Xm |
… |
Xn |
C1' |
C1" |
X1 |
а10 |
1 |
… |
0 |
… |
а1n |
... |
… |
… |
… |
0 |
… |
0 |
… |
… |
Cm' |
Cm" |
Xm |
аm0 |
0 |
… |
1 |
… |
аmn |
|
j' |
∆0' |
0 |
… |
0 |
… |
∆n' |
|
j" |
∆0" |
0 |
… |
0 |
… |
∆n" |
||
t = |
j |
∆0 |
0 |
… |
0 |
… |
n |
В первых двух столбцах таблицы записываются коэффициенты целевой функции при базисных переменных, учитывая, что . В первых двух строках записываются коэффициенты целевой функции при всех переменных и свободный член. Дополнительно вводятся две строки, в которых подсчитываются j', j" по формулам:
2. Берём t = , где - начальное значение параметра t (или t = tнижн., т.е.нижней границе отрезка), при условии, что - ∞.
Находим j по формуле: j = j'+t j" и заполняем последнюю строку таблицы.
3. Проверяем решение на оптимальность (при решении на максимум все j при переменных должны быть неотрицательные).
Если решение неоптимальное, то решаем задачу симплексным методом, выбирая разрешающий столбец по наименишей отрицательной общей оценке j, полученной при t = , (или t = tнижн., т.е.нижней границе отрезка), до получения оптимального плана (решения).
Если решение оптимальное, то переходим к пункту 4.
4. Если получено оптимальное решение, то определяем отрезок значений параметра t, для которых полученное решение сохраняет оптимальность.
Оптимальное решение, записанное в таблице, будет оптимальным до тех пор, пока неотрицательны все оценки j при переменных, т.е. при определении отрезка значений параметра t составляем и решаем систему линейных неравенств:
1'+t
1"
0
2'+t 2" 0
…
j'+t j" 0
Находим пределы изменения параметра t:
tнижн. t tверхн.
5. Выписываем оптимальное решение на найденном отрезке из таблицы.
Если tверхн. , то решение закончено, если tверхн. , то переходим к пункту 6.
6. Для нахождения оптимального решения на следующем отрезке определяем разрешающий столбец.
Берем t = tверхн. , которое является верхней границей предыдущего (найденного) интервала и нижней границей искомого (нового) интервала изменения параметра t.
Находим столбец, по которому j = j'+t j"= 0. Этот столбец является разрешающим. При небольшом увеличении t (t > tверхн.) оценка j* станет отрицательной, найденное оптимальное решение потеряет свою оптимальность.
7. Выполняем симплексное преобразование однократного замещения с выбранным разрешающим столбцом и получим новое оптимальное решение, в котором j вычисляются для t = tнижн., искомого отрезка. Далее переходим к пункту 4.
Процесс решения задачи продолжается до тех пор, пока не закончится разбиение всего заданного интервала [;] на отрезки.
Замечание. Если задано t, то для начала исследования берём любую точку числовой оси ( принадлежащую заданному интервалу) и проводим исследования в обе стороны от первоначально полученного отрезка.
Пример выполнения работы
Найти max Z = (1+t) x1+(1-t) x2
x
1+x2
2
x1 1
x1 0, x2 0
0 t 20
Последовательность решения.
Переходим к канонической форме записи:
x 1+x2+x3 =2
x1+x4 = 1
x1 0, x2 0, x3 0, x4 0
max Z = (1+t) x1+(1-t) x2+(0+t0) x3+(0+t0) x4
1. Записываем исходное опорное решение в симплексную таблицу.
№ 1 |
Cj' |
0 |
1 |
1 |
0 |
0 |
|
|
Cj" |
0 |
1 |
-1 |
0 |
0 |
|
||
Ci' |
Ci" |
П БП |
аi0 |
X1 |
X2 |
X3 |
X4 |
СО |
0 |
0 |
X3 |
2 |
1 |
1 |
1 |
0 |
2:1=2 |
0 |
0 |
X4 |
1 |
1 |
0 |
0 |
1 |
1:1=1 * |
|
|
j' |
0 |
-1 |
-1 |
0 |
0 |
|
|
|
j" |
0 |
-1 |
1 |
0 |
0 |
|
t = 0 |
j |
0 |
-1 * |
-1 |
0 |
0 |
|
СО – симплексные отношения
2.Берем t = 0, т.е. начальному значению интервала изменения параметра. Находим j по формуле: j = j'+t j"
0 = 0+0 0=0; 1= -1+0 (-1)=-1; 2= -1+0 1=-1 и т.д.
3. Проверяем решение на оптимальность (при решении на максимум все j при переменных должны быть неотрицательные).
В нашем примере еще не получено оптимальное решение, т.к. 1, и 2 отрицательные. Решаем задачу симплексным методом до получения оптимального решения, которое получается в третьей симплексной таблице.
№ 2 |
Cj' |
0 |
1 |
1 |
0 |
0 |
|
|
Cj" |
0 |
1 |
-1 |
0 |
0 |
|
||
Ci' |
Ci" |
П БП |
аi0 |
X1 |
X2 |
X3 |
X4 |
СО |
0 |
0 |
X3 |
1 |
0 |
1 |
1 |
-1 |
1:1=1 * |
1 |
1 |
X1 |
1 |
1 |
0 |
0 |
1 |
- |
|
|
j' |
1 |
0 |
-1 |
0 |
1 |
|
|
|
j" |
1 |
0 |
1 |
0 |
1 |
|
t = 0 |
j |
1 |
0 |
-1 * |
0 |
1 |
|
№ 3 |
Cj' |
0 |
1 |
1 |
0 |
0 |
|
|
Cj" |
0 |
1 |
-1 |
0 |
0 |
|
||
Ci' |
Ci" |
П БП |
аi0 |
X1 |
X2 |
X3 |
X4 |
СО |
1 |
-1 |
X2 |
1 |
0 |
1 |
1 |
-1 |
1:1 =1 * |
1 |
1 |
X1 |
1 |
1 |
0 |
0 |
1 |
- |
|
|
j' |
2 |
0 |
0 |
1 |
0 |
|
|
|
j" |
0 |
0 |
0 |
-1 |
2 |
|
t = 0 |
j |
2 |
0 |
0
|
1 * |
0 |
|
4. Определяем отрезок значений параметра t, для которых полученное решение сохраняет оптимальность. Для этого составляем систему неравенств:
3=
1+ t (-1)
0 1- t
0
4= 0+ t 2 0 2t 0
0 t 1
5. Выписываем оптимальное решение на отрезке 0 t 1.
,
Z1*=2
+ t0
= 2
tверхн.=1. Так как tверхн.=1 меньше 20 (конечного значения параметра t в заданном интервале), то решение параметрической задачи еще не закончено.
6. Для нахождения оптимального решения на следующем отрезке определяем разрешающий столбец.
Берем t = 1. Находим столбец, по которому j = j'+t j"= 0. Этот столбец является разрешающим.
3 = 1+ 1 (-1) =0
4 = 0+1 2 = 2
Так как 3 =0, то разрешающим столбцом будет столбец X3 (см. табл. 3).
7. Выполняем симплексное преобразование однократного замещения с выбранным разрешающим столбцом и получим новое оптимальное решение.
№ 4 |
Cj' |
0 |
1 |
1 |
0 |
0 |
|
Cj" |
0 |
1 |
-1 |
0 |
0 |
||
Ci' |
Ci" |
П БП |
аi0 |
X1 |
X2 |
X3 |
X4 |
0 |
0 |
X3 |
1 |
0 |
1 |
1 |
-1 |
1 |
1 |
X1 |
1 |
1 |
0 |
0 |
1 |
|
|
j' |
1 |
0 |
-1 |
0 |
1 |
|
|
j" |
1 |
0 |
1 |
0 |
1 |
t = 1 |
j |
2 |
0 |
0 |
0 |
2 |
8. Определяем отрезок значений параметра t, для которых полученное решение сохраняет оптимальность. Для этого составляем систему неравенств:
2= -1+ t 1 0 -1+ t 0
4= 1+ t 1 0 1+t 0
1 t ∞
9. Выписываем оптимальное решение на отрезке 1 t ∞
,
Z2*=1
+ t1
=1+t
10. Так как последний интервал изменения параметра равен 1 t ∞, то = 20 входит в найденный интервал. Решение закончено.
Ответ: , Z1*=2 на отрезке 0 t 1.
, Z2*=1+t на отрезке 1 t ∞, в т.ч. на отрезке 1 t 20
Индивидуальное задание
Общая постановка задачи
Задача определения оптимального объема производства продукции при изменении условий реализации.
Предприятие выпускает два вида продукции, для изготовления которых используются различные виды сырья. Известны расходы сырья каждого вида на производство единицы данного вида продукции и запасы сырья на предприятии. Установлено, что цена единицы продукции может изменяться в определенных пределах. Сырье может быть недоиспользовано.
Найти оптимальный план производства продукции в условиях изменения цен реализации такой, чтобы обеспечить максимум общей стоимости произведенной продукции.
Задание. По исходным данным записать задачу параметрического программирования в исходной форме и выполнить решение согласно алгоритму.
Список индивидуальных данных:
Таблица 2. Варианты индивидуальных заданий.
№ варианта |
Вид сырья |
Вид продукции |
№ варианта |
Вид сырья |
Вид продукции |
1 |
А,В,С |
1,2 |
14 |
А,В,С |
2,4 |
2 |
А,С,Д |
1,2 |
15 |
В,С,Д |
2,4 |
3 |
А,В,С |
2,3 |
16 |
А,С,Д |
1,3 |
4 |
А,С,Д |
2,3 |
17 |
А,В,Д |
2,5 |
5 |
А,В,С |
3,4 |
18 |
В,С,Д |
2,5 |
6 |
А,С,Д |
3,4 |
19 |
А,В,С |
1,4 |
7 |
А,В,С |
1,3 |
20 |
В,С,Д |
1,4 |
8 |
А,С,Д |
1,4 |
21 |
А,В,Д |
3,4 |
9 |
А,В,С |
2,5 |
22 |
В,С,Д |
3,4 |
10 |
А,С,Д |
2,5 |
23 |
А,В,Д |
2,3 |
11 |
В,С,Д |
1,3 |
24 |
В,С,Д |
2,3 |
12 |
А,В,Д |
1,3 |
25 |
А,В,Д |
1,2 |
13 |
А,С,Д |
2,4 |
26 |
В,С,Д |
1,2 |
Таблица 3. Исходные данные
Вид сырья |
Нормы расхода сырья на единицу продукции, усл.ед. |
Запасы сырья, усл.ед. |
||||
виды продукции |
||||||
1 |
2 |
3 |
4 |
5 |
||
А |
1 |
2 |
3 |
1 |
2 |
4 |
В |
1 |
1 |
2 |
3 |
1,5 |
3 |
С |
2 |
2 |
1 |
3 |
1 |
8 |
Д |
1,5 |
1 |
1 |
2 |
1,5 |
6 |
Цена изделия, ден.ед. |
2+t |
20- t |
1+2 t |
15- t |
3+ t |
0 t 10 |
Контрольные вопросы
Какие изменения вносятся в симплексную таблицу при решении задачи линейного программирования с параметром в целевой функции?
По какой формуле вычисляются оценки j?
Каким образом получают оптимальное решение в симплексной таблице?
Как определяют отрезок значений параметра t, для которого в таблице получено оптимальное решение?
Как выписывается из таблицы оптимальное решение на найденном отрезке?
Как выбирается разрешающий столбец для получения оптимального решения на следующем отрезке?
Какие действия выполняются после определения разрешающего столбца?
Как долго продолжается процесс решения задачи?
Лабораторная работа № 25. «Параметрическая транспортная задача с параметром в целевой функции»
Теоретическая часть
При решении транспортных задач мы рассматривали перевозку однородного груза. Часто на практике решают параметрические транспортные задачи с параметром в целевой функции, то есть коэффициенты транспортных затрат заданы с параметром. Так при перевозке груза разными марками автомашин с неодинаковой грузоподъёмностью затраты на перевозку одной тонны груза могут изменяться. В этом случае возможность изменения коэффициентов транспортных затрат учитывается с помощью параметра. Параметр t задаётся на интервале [,], – нижняя (начальная) граница интервала, – верхняя (конечная) граница данного интервала.
Алгоритм решения параметрической транспортной задачи с параметром в целевой функции
Записывают условие транспортной задачи в распределительную таблицу, причём коэффициенты транспортных затрат - с параметром t.
Проверяют, является ли модель задачи закрытой. Если да, то переходят к пункту 3, если нет, то введением фиктивного отправителя или фиктивного потребителя сводят задачу к закрытой модели.
Если навык решения транспортной задачи есть, то мысленно, если нет, то письменно полагают t=t нижней границе рассматриваемого интервала.
Любым известным методом находят оптимальный план (оптимальное решение) грузоперевозок при этом t.
Выписывают характеристики свободных переменных (клеток) для оптимального решения и составляют неравенства, наложив на них условие неположительности.
Из решения системы неравенств неположительности характеристик свободных клеток определяют интервал, на котором найденный план грузоперевозок останется оптимальным (t нижняя граница t t верхняя граница).
Проводят преобразование однократного замещения по циклу свободной клетки, характеристика которой при подстановке t=tверхняя граница обратится в ноль, а при подстановке ttверхняя граница станет положительной.
Переход к пункту 5. Процесс решения выполняется до тех пор, пока не будет исследован весь заданный интервал изменения параметра t.
Замечание. Преобразование однократного замещения по циклу свободной клетки – это нахождение = min{xij} в четных клетках цикла, перемещение по вершинам цикла (прибавление к переменным в нечётных вершинах цикла и вычитание из переменных в чётных вершинах цикла). Переменная, находящаяся в свободной клетке цикла, становится базисной, свободная клетка цикла становится занятой, а одна из базисных переменных (из клеток) цикла становится свободной. При этом ранг задачи не меняется. Баланс строчек и столбцов таблицы не нарушается.
Пример выполнения работы
Задание. Записать транспортную задачу с параметром в целевой функции в распределительную таблицу и решить её.
Необходимо перевезти груз от трёх поставщиков пяти потребителям с минимальными затратами на грузоперевозку. Ресурсы, потребности и коэффициенты транспортных затрат даны в таблице. Параметр t принадлежит интервалу [0,1].
-
Потребности
Ресурсы
b1
10
b2
20
b3
25
b4
15
b5
30
а1 20
1+t
12-t
2t
1+2t
2+t
а2 30
2
1+t
2t+1
4-t
6-t
а3 50
t+1
4-t
3-t
t+2
2t-1
Найдём исходное опорное решение методом минимального элемента в таблице, приняв t=0.
-
Потребности
Ресурсы
10
20
25
15
30
Ui
20
1+t
12-t
2t
20
1+2t
2-t
0
30
2
1+t
20
2t+1
5
4-t
5
6-t
1
50
t+1
10
4-t
3-t
t+2
10
2t-1
30
2t-1
Vj
-t+2
t
2t
3-t
0
Проверяем решение на оптимальность. Находим потенциалы по формуле: Cij= Ui +Vj , где U1 =0, Cij - коэффициенты транспортных затрат базисных клеток (переменных).
Вычисляем характеристики свободных клеток по формуле:
ij Ui +Vj)- Cij. Тогда 0-t+2)-(1+t)=1-2t; 2(0+t)- (12-t)=2t-12;
(0+3-t)- (1+2t)=2-3t; (0+0)- (2-t)=t-2; 21(1-t+2)- 2=1-t;
25(0+1)- (6-t)=t-5; 32(2t-1+t)-(4-t)=4t-5; (2t-1+2t)- (3-t)=5t-4.
Решение неоптимальное, так как при t=0 есть неотрицательные характеристики (оценки) свободных клеток. Это =1; 2; 211; выбираем среди них максимальную. 2 - максимальная неотрицательная характеристика. Для этой клетки строим замкнутый цикл: [1,4]-[1,3]-[2,3]-[2,4]-[1,4], имеющий соответственно знаки вершин: +, -,+, -. Перемещаем по циклу. Строим новую таблицу
-
Потребности
Ресурсы
10
20
25
15
30
Ui
20
1+t
12-t
2t
15
1+2t
5
2-t
0
30
2
1+t
20
2t+1
10
4-t
6-t
1
50
t+1
10
4-t
3-t
t+2
10
2t-1
30
-t+1
Vj
2t
t
2t
1+2t
3t-2
Проверяем решение на оптимальность. Находим потенциалы по формуле: Cij= Ui +Vj , где U1 =0, Cij - коэффициенты транспортных затрат базисных клеток (переменных).
Вычисляем характеристики свободных клеток по формуле:
ij Ui +Vj)- Cij. Тогда 0+2t)-(1+t)= -1+t; 2(0+t)- (12-t)=2t-12;
(0+3t-2)- (2-t)=4t-4; 21(1+2t)- 2= -1+2t; (1+1+2t)- (4-t)=--2+3t;
25(1+3t-2)- (6-t)=4t-7; 32(-t+1+t)-(4-t)=t-3; (2t-t+1)- (3-t)=2t-2.
Решение оптимальное, так как при t=0 все характеристики (оценки) свободных клеток неположительные.
Находим интервал изменения t, на котором решение остаётся оптимальным. Выписываем характеристики свободных переменных (клеток) для оптимального решения и составляем неравенства, наложив на них условие неположительности. ij Ui +Vj)- Cij 0.
Получаем систему:
-1+t0; 2t-120; 4t-40; -1+2t0; -2+3t0; 4t-70; t-30; 2t-20. Откуда: t1; t; t1; t0,5; t2/3; t7/4; t3; t1.
Решением системы является t(-; ½]. По условию t[0; 1]. Поэтому суживаем интервал изменения t. Получаем при t[0; ½] сохраняется оптимальное решение
X=
min Z= 2t1+2t)1+t)2t+1)t+1)10 +(t+2)10 +(2t-1) 30=
=30t +5+10t+20+20t+20t+10+10 t+10+10 t+20+60t-30=160 t+35.
Переходим к новой распределительной таблице. Проводим преобразование однократного замещения по циклу свободной клетки, характеристика которой при подстановке t=tверхняя граница обратится в ноль.
Это клетка [2,1]. Строим для неё цикл: [2,1]-[2,3]-[1,3]-[1,4]-[3,4]-[3,1]-[2,1]. Из переменных в нечётных клетках (с отрицательными вершинами цикла) выбираем Перемещаем по циклу. Заполняем следующую таблицу.
-
Потребности
Ресурсы
10
20
25
15
30
Ui
20
1+t
12-t
2t
20
1+2t
2-t
0
30
2
5
1+t
20
2t+1
5
4-t
6-t
1
50
t+1
5
4-t
3-t
t+2
15
2t-1
30
t
Vj
1
t
2t
2
t-1
Проверяем балансы строк, балансы столбцов и ранг. Всё совпадает; таблица построена верно. Проверяем решение на оптимальность. Подсчитываем потенциалы строк и столбцов и вносим их в таблицу. Находим характеристики (оценки) свободных клеток.
ij Ui +Vj)- Cij. Тогда 0-1)-(1+t)= -t; 2(0+t)- (12-t)=2t-12;
(0+2)- (1+2t)=1-2t; (0+ t-1)- (2-t)=2t-3; 24(1+2)- (4-t )=-1+ t;
25(1+ t-1)- (6-t)=2t-6; 32(t+t)-(4-t)=3t-4; (t+2t)- (3-t)=4t-3.
При t=1/2 все характеристики неположительные. Следовательно, решение оптимальное.
Находим интервал изменения t, на котором решение остаётся оптимальным. Выписываем характеристики свободных переменных (клеток) для оптимального решения и составляем неравенства, наложив на них условие неположительности. ij Ui +Vj)- Cij 0.
Получаем систему:
ij Ui +Vj)- Cij. Тогда 0-1)-(1+t)= -t; 2(0+t)- (12-t)=2t-12;
(0+2)- (1+2t)=1-2t; (0+ t-1)- (2-t)=2t-3; 24(1+2)- (4-t )=-1+ t;
25(1+ t-1)- (6-t)=2t-6; 32(t+t)-(4-t)=3t-4; (t+2t)- (3-t)=4t-3.
-t0; 2t-120; 1-2t0; 2t-30; -1+ t0; 2t-60; 3t-40; 4t-30. Откуда: t0; t; t1/2; t3/; t1; t3; t3. t3/.
Решением системы является t[1/2; 3/4]. По условию t[0; 1].
Поэтому весь интервал не исследован.
Получаем при t[1/2; 3/4] сохраняется оптимальное решение
X=
min Z= 2t21+t)2t+1)t+1) +(t+2)15 +(2t-1) 30=
=40t +10+20+20t+10t+5+5t+5+15t+30+60t-30=150 t+40.
Переходим к новой распределительной таблице. Проводим преобразование однократного замещения по циклу свободной клетки, характеристика которой при подстановке t=tверхняя граница обратится в ноль.
Это клетка [3,3]. Строим для неё цикл: [3,3]-[2,3]-[2,1]-[3,1]-[3,3]. Из переменных в нечётных клетках (с отрицательными вершинами цикла) выбираем Перемещаем по циклу. Заполняем следующую таблицу.
-
Потребности
Ресурсы
10
20
25
15
30
Ui
20
1+t
12-t
2t
20
1+2t
2-t
0
30
2
10
1+t
20
2t+1
4-t
6-t
4-4t
50
t+1
0
4-t
3-t
5
t+2
15
2t-1
30
3-3t
Vj
4t-2
5t-3
2t
4t-1
5t-4
Проверяем балансы строк, балансы столбцов и ранг. Всё совпадает; таблица построена верно. Проверяем решение на оптимальность.
Находим интервал изменения t, на котором решение остаётся оптимальным. Выписываем характеристики свободных переменных (клеток) для оптимального решения и составляем неравенства, наложив на них условие неположительности. ij = (Ui +Vj)- Cij 0.
Получаем систему:
ij Ui +Vj)- Cij. Тогда 0+4t-2)-(1+t)= 3t-3; 2(0+5t-3)- (12-t)=6t-15; (0+4t-1)- (1+2t)= -2+2t; (0+ 5t-4)- (2-t)=6t-6; (4-4t+2t)- (2t+1)= -4t+3. 24(4-4t+4t-1)- (4-t )= -1+ t; 25(4-4t+5t-4)-(6-t)=2t-6; 32(3-3t+5t-3)-(4-t)=3t-4;
3t-30; 6t -150; -2+2t0; 6t-60; -4t+30; -1+t0; 2t-60; 3t-40; Откуда: t1; t5/2; t1; t1; t3/4; t1; t3; t3.
Решением системы является t[3/4;1]. По условию t[0; 1]. Так как t=tверхняя граница=то поэтому весь интервал исследован.
Получаем, что при t[3/4;1] сохраняется оптимальное решение. Кроме того решение вырожденное, так как одна из базисных переменных (x31 =0) равна нулю.
X=
min Z= 2t21+t)+( t+1)0 +(3-t)5t+2)1 +(2t-1) 30=
=40t +20+20+20t+0+15-5t+15t+30+60t-30=130t+55.
Ответ. Заданный интервал изменения параметра t (t[0; 1]) разбит на три интервала:
на 1-ом интервале t[0; ½] получено оптимальное решение
X=
min Z= 160 t+35;
на 2-ом интервале t[1/2; 3/4] получено оптимальное решение
X=
min Z= 150 t+40;
на 3-ем интервале t[3/4;1] найдено вырожденное оптимальное решение
X=
min Z= 130t+55.
Индивидуальное задание
Общая постановка задачи
Записать транспортную задачу с параметром в целевой функции в распределительную таблицу. Найти оптимальные планы перевозки однородного груза от поставщиков к потребителям, обеспечивающие минимальные транспортные расходы и учитывающие параметры в целевой функции. Исследовать весь интервал изменения параметра.
Список индивидуальных данных
Получить исходное опорное решение методами: или «северо-западного» угла, или наименьшего элемента в строке, в столбце, в таблице. Подсчитать значение целевой функции при полученном исходном опорном решении.
Получить оптимальное решение методом потенциалов, улучшая исходное опорное решение, полученное в п.1.
Таблица 1. Варианты заданий к транспортной задаче с параметром в целевой функции
Поставщики |
Ресурсы |
Потребители |
||||||
а |
б |
в |
г |
д |
е |
ж |
||
Потребности |
||||||||
50 |
100 |
70 |
80 |
50 |
120 |
150 |
||
А |
100 |
4-2t |
3+t |
5-t |
6+3t |
2-4t |
1+2t |
8-6t |
Б |
150 |
2-t |
1+3t |
8-5t |
3+4t |
7-t |
6-3t |
4+2t |
В |
50 |
7-5t |
3+t |
9-6t |
10-7t |
8-4t |
4+2t |
7-5t |
Г |
20 |
3+2t |
4+t |
1+5t |
5-2t |
10-8t |
9-6t |
6-3t |
Д |
50 |
6-5t |
7-4t |
2+t |
3+2t |
4-t |
8-3t |
5+t |
Е |
180 |
5-3t |
9-5t |
4+2t |
6-3t |
7-4t |
2+2t |
3+4t |
Таблица 2. Значение параметров для целевой функции
Номер варианта, [,] |
Поставщики |
Потребители |
Номер варианта, [,] |
Поставщики |
Потребители |
1, [-2,20] |
А, Б, В |
а, б, в, г |
13, [0,22] |
А, Б, В |
б, д, ж |
2, [-5,15] |
Б, В, Г |
а, б, в |
14, [-1,12] |
В, Г, Д, Е |
б, д, ж |
3, [-3,15] |
Б, В, Д |
б, в, г |
15, [-2,18] |
В, Г, Д, Е |
в, г, ж |
4, [-4,20] |
В, Г, Д, Е |
а, б, в, г |
16, [-3,20] |
А, Б, Д |
а, б, ж |
5, [-1,22] |
А, В, Г, Д |
а, б, в |
17, [-3,21] |
А, Г, Е |
б, в, г, д |
6, [0,25] |
А, Б, Г |
а, б, в, д |
18, [-1,14] |
Б, В, Г, Е |
б, в, г, ж |
7, [-3,23] |
Г, Д, Е |
б, в, г |
19, [-2,19] |
А, Б, В |
а, г, д, е |
8, [-1,21] |
Б, В, Д |
а, в, г, д |
20, [-3,20] |
А, Б, В |
а, г, д, е |
9, [-2,23] |
Г, Д, Е |
а, в, г, д |
21, [-1,21] |
В, Г, Е |
а, д, ж |
10, [-3,22] |
Б, В, Г, Д |
а, б, в, д |
22, [-4,18] |
А, Г, Д, Е |
а, б, д, ж |
11, [-4,25] |
В, Г, Д, Е |
б, г, е |
23, [0,25] |
Б, В, Г |
а, д, е |
12, [-2,24] |
А, Б, В |
б, г, е |
24, [-5,15] |
В, Д, Е |
б, г, д, ж |
Контрольные вопросы
Всегда ли параметрическая транспортная задача с параметром в целевой функции имеет оптимальное решение (план)?
Какая параметрическая транспортная задача называется: закрытой, открытой?
Как свести открытую параметрическую транспортную задачу к закрытой?
Чему равно число основных (базисных) переменных закрытой параметрической транспортной задачи?
Как найти исходное опорное решение параметрической транспортной задачи методом «северо-западного угла»?
Как найти исходное опорное решение параметрической транспортной задачи методом наименьшего элемента?
Как проверить найденное опорное решение параметрической транспортной задачи на оптимальность?
Опишите метод потенциалов, используемый для нахождения оптимального плана перевозок параметрической транспортной задачи .
Как осуществляется выбор переменной, выводимой из списка базисных переменных? Правила построения цикла. Как осуществляется выбор ?
Как осуществляется выбор переменной, вводимой в список базисных переменных?
Как находиться интервал изменения параметра t, в котором сохраняется оптимальность решения?
Глоссарий основных понятий
А |
Алгоритм симплекс-метода - алгоритм последовательного улучшения плана, позволяющий осуществлять переход от одного допустимого базисного решения к другому таким образом, что значение целевой функции непрерывно возрастает (при решении на максимум) и за конечное число шагов находится оптимальное решение. Алгоритм улучшения плана транспортной задачи - алгоритм перехода к новому опорному плану транспортной задачи, дающему меньшее значение функции потерь, до обнаружения оптимального плана. Альтернативный оптимум - если задача имеет несколько оптимальных решений с разным значением переменных, но с одинаковым значением целевой функции, то она имеет альтернативные оптимальные решения (альтернативный оптимум).Антагонистические игры |
Б |
Базисная переменная - переменная, коэффициенты которой образуют единичный столбец. Базисное решение - система ограничений приведена к единичному базису, а все свободные переменные равны нулю. |
В |
Вектор коэффициентов - вектор, компонентами которого являются коэффициенты целевой функции задачи линейного программирования. Вектор ограничений - вектор, компонентами которого являются ограничения выражений, определяющих допустимую область задачи линейного программирования. Вершина выпуклого многогранника - это любая точка выпуклого многогранника, которая не является внутренней ни для какого отрезка целиком принадлежащего этому многограннику. Возможное решение - решение, удовлетворяющее системе ограничений задачи. Выпуклая комбинация точек - точка, компоненты которой представлены суммой произведений неотрицательных коэффициентов не больших единицы и соответствующих компонент данных точек, при этом сумма всех коэффициентов равна единице. Выпуклое множество - множество, которое вместе с двумя принадлежащими ему точками обязательно содержит отрезок, соединяющий эти точки. Вырожденный опорный план - опорный план, число ненулевых компонент которого меньше числа ограничений. |
Г |
Геометрическая интерпретация задачи линейного программирования - интерпретация зависимостей, имеющих место в задаче линейного программирования в виде геометрических фигур (точек, прямых, полуплоскостей, многоугольников) в декартовой системе координат. |
Д |
Двойственные задачи линейного программирования - задачи линейного программирования, которые могут быть составлены из исходных задач линейного программирования. Дельта-метод - один из методов проверки опорного плана транспортной задачи на оптимальность. Допустимая область задачи линейного программирования - множество опорных планов задачи линейного программирования. Допустимое решение - допустимым решением называется совокупность значений n переменных, удовлетворяющая системе ограничений и условиям неотрицательности. |
Е |
Единичный базис системы - базисные переменные образуют единичный базис системы. |
З |
Задача линейного программирования - характеризуется тем, что целевая функция является линейной функцией переменных, а область допустимых значений определяется системой линейных равенств или неравенств. Задача математического программирования
В
общей постановке задачи этого раздела
выглядят следующим образом. Имеются
какие-то переменные
Задача математического программирования (МП) состоит из системы ограничений, условий неотрицательности переменных и целевой функции. Задача о составлении плана производства - возникает при необходимости максимизации дохода от реализации продукции, производимой некоторой организацией, при этом производство ограничено имеющимися сырьевыми ресурсами. |
И |
Исследование операций - наука, занимающаяся разработкой и практическим применением методов наиболее оптимального управления организационными системами. |
К |
Каноническая форма задачи линейного программирования - форма задачи линейного программирования, в которой все ограничения имеют тип "="и переменные неотрицательные. Критерий оптимальности - показатель, количественно определяющий цель решения модели. |
Л |
Линейное программирование - часть математического программирования, задачами которой является нахождение экстремума линейной целевой функции на допустимом множестве значений аргументов. |
М |
Математическое программирование - раздел современной математики, задачами которого является нахождение экстремума функции при условии принадлежности переменных определенному множеству. Матрица коэффициентов - матрица, элементами которой являются коэффициенты системы линейных равенств или неравенств определенного типа. Матричная форма задачи линейного программирования - форма задачи линейного программирования, когда все элементы задачи представлены в матричных и векторных обозначениях. Метод искусственного базиса (М-метод) - один из методов, упрощающий определение исходного опорного плана задачи линейного программирования и симплекс-таблицы. Метод минимального элемента - один из группы методов определения первоначального опорного плана транспортной задачи. Метод потенциалов - один из методов проверки опорного плана транспортной задачи на оптимальность. Метод северо-западного угла - один из группы методов определения первоначального опорного плана транспортной задачи. |
Н |
Невырожденный опорный план - план, соответствующий вершине допустимой области, который имеет m отличных от нуля компонент, где m есть количество ограничений задачи линейного программирования. Неопределенная система -совместная система называется неопределенной, если она имеет более одного решения. Несовместная система ограничений - система ограничений называется несовместной, если она не имеет ни одного решения. |
О |
Общее решение - решение называется общим, если базисные переменные выражены через свободные члены и свободные переменные. Опорное решение - опорным решением называется базисное решение, в котором базисные переменные неотрицательные или - решение называется опорным, если система ограничений приведена к единичному базису, все свободные переменные равны нулю, а базисные переменные неотрицательные. Определенная совместная система Совместная система называется определенной, если она имеет только одно решение. Оптимальный план ЗЛП - решение задачи линейного программирования, т. е. такой план, который входит в допустимую область и доставляет экстремум целевой функции. Основная теорема линейного программирования - если целевая функция принимает максимальное значение в некоторой точке допустимой области, то она принимает это же значение в крайней точке допустимой области. Если целевая функция принимает максимальное значение более, чем в одной крайней точке, то она принимает это же значение в любой их выпуклой комбинации. Открытая транспортная задача - несбалансированная транспортная задача. Отрезок - множество точек, которые могут быть представлены в виде выпуклой комбинации данных двух точек. |
П |
План - набор чисел, удовлетворяющий ограничениям задачи линейного программирования. Признак вершины допустимой области Если система из k ненулевых векторов-столбцов, образованных соответствующими столбцами матрицы ограничений является линейно независимой и ненулевые координаты точки X, удовлетворяют ограничениям, то эта точка является вершиной допустимой области. Принцип последовательного «улучшения» опорных решений
|
Р |
Ранг совместной системы -рангом совместной системы линейных ограничений называется ранг её матрицы. Решение задачи линейного программирования - это план, доставляющий экстремальное значение целевой функции. |
С |
Свободные переменные -переменные, не входящие в единичный базис системы. Сбалансированная транспортная задача -транспортная задача, в которой сумма ресурсов равняется сумме потребностей Симплекс-метод - последовательное улучшение плана задачи линейного программирования, позволяющее осуществлять переход от одного допустимого базисного решения к другому, причем так, что значения целевой функции непрерывно возрастают и за конечное число шагов находится оптимальное решение. Совместная система ограничений - система ограничений называется совместной, если она имеет хотя бы одно решение. |
Т |
Теорема двойственности Первая теорема двойственности: Если одна из взаимно двойственных задач имеет оптимальное решение, то его имеет и другая, причем оптимальные значения их линейных функций равны, т.е. max Z = min F. Если линейная функция одной из задач не ограничена, то условия другой задачи противоречивы. Вторая теорема двойственности: Положительным (ненулевым) компонентам оптимального решения одной из взаимно двойственных задач соответствуют нулевые компоненты оптимального решения другой задачи, т.е. для любых i= 1,2,…, m и j=1,2,…,n: если Xj*>0, то ym+j*=0; если Xn+I*>0, то yi*=0. Теорема о выпуклом множестве и выпуклой комбинации этого множества Пусть G - выпуклое множество. Тогда любая выпуклая комбинация точек, принадлежащих этому множеству, также принадлежит этому множеству. Теорема о выпуклости допустимого множества ЗЛП Допустимая область задачи линейного программирования является выпуклым множеством. Теорема о выпуклости оптимальных планов ЗЛП Множество оптимальных планов задачи линейного программирования выпукло (если оно не пусто). Теорема о том, что любая точка выпуклого многогранника является выпуклой комбинацией вершин Любая точка выпуклого многогранника является выпуклой комбинацией его вершин. Транспортная задача Пусть имеется однородный продукт, распределенный в определенных количествах (не обязательно одинаковых) в m складах. Этот продукт необходимо доставить в n пунктов потребления, причем в каждый пункт установленное количество. Запасы и потребности сбалансированы. Стоимость перевозки из конкретного склада в конкретный пункт индивидуальна. Товар должен быть вывезен из всех складов и доставлен в требуемом количестве в каждый пункт. Задача заключается в минимизации транспортных расходов. |
Ц |
Целевая функция - функция в математическом программировании, для которой требуется найти экстремум. Математическое выражение критерия оптимальности. |
Ч |
Частное решение - решение называется частным, если в общем решении свободные переменные принимают любые значения. Число базисных решений Число базисных решений не превышает числа сочетаний из n по m, где n - число переменных, а m - число линейно независимых ограничений Число опорных решений - число опорных решений не превышает числа базисных решений. |
Рекомендуемая литература
а) Основная литература
Абрамов Л.М., Капустин В.Ф. Математическое программирование. Л.: Изд-во ЛГУ, 1981
Калихман И.Л. Линейная алгебра и программирование. - М.: Высшая школа, 1967
Красс М.С.Математические методы и модели. М.: Финансы и статистика, 2007
Красс М.С., Чупрынов Б.П. Основы математики и ее приложения в экономическом образовании. М.: ДЕЛО, 2001
Красс М.С.,Чупрынов Б.П. Основы математики. М.: Финансы и статистика, 2005
Кремер Н.Ш. и др. Исследование операций в экономике. М.: ЮНИТИ, 1997, 2003, 2004
Кузнецов Б.Т. Математические методы и модели исследования операций. Гриф УМЦ "Профессиональный учебник".-М.: ЮНИТИ, 2005
Кузнецов Ю.Н., Кузубов В.И., Волощенко А.Б. Математическое программирование. М.: Высшая школа, 1980,1986
Математическое программирование. /Под ред. Н.Ш. Кремера. – М.: Финстатинформ, 1995
Математические методы в экономике: учеб. пособие/Л.Э. Хазанова. – 3-е изд., стереотип. – М., Волтерс Клувер, 2005. – 144с.
Шапкин А.С., Мазаева Н.П. Математические методы и модели исследования операций. Учебник.-М.: Дашков и К°, 2005
б) Дополнительная литература
Акулич И.Л. Математическое программирование в примерах и задачах. – М.: Высшая школа, 1986
Астафьев Н.Н. Линейные неравенства и выпуклость. М.: Наука, 1982
Банди Б. Основы линейного программирования. М.: Радио и связь, 1989
Волков Е.А. Численные методы. М.: Наука, 1982
Гольштейн Е.Г., Юдин Д.Б. Линейное программирование (теория, методы и приложения). – М.: Наука, 1969
Карпилевич Ф.И., Садовский Л.Е. Элементы линейной алгебры и линейного программирования. – М.: Наука, 1967
Кофман А., Фор Р. Займемся исследованием операций. –М.: Мир, 1966
Максимилиан С.В., Замбицкий М.К. Линейное программирование. Кишинев: Штиинца, 1981
Мустаф, Брюс. Современное линейное программирование: Теория и практика. М.: Мир, 1984
Поляк Б.Т. Введение в оптимизацию. М.: Наука, 1983
Пшеничный Б.Н. Необходимые условия экстремума. М.: Наука, 1982
Солодовников А.С., Бабайцев В.А., Браилов А.В. Математика в экономике. Ч.1.- М.: Финансы и статистика, 2001,2007
1 Рекомендуется перед началом оптимизации вводить в эти ячейки начальные значения, максимально близкие к предполагаемому оптимальному значению — это ускоряет поиск решения и повышает шансы на его успешное завершение.