МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
УЧРЕЖДЕНИЕ ОБРАЗОВАНИЯ
«БРЕСТСКИЙ ГОСУДАОСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Кафедра информатики и прикладной математики
КУРСОВАЯ РАБОТА
на тему
«РЕШЕНИЕ ЗАДАЧ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ»
по дисциплине «Системный анализ и исследование операций»
Выполнил студент
Допущен к защите
Результаты защиты
БРЕСТ 2009
СОДЕРЖАНИЕ
Введение………………………………………………………………………………
Постановка задачи оптимизации……………………………………….
ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ………………………………...
Обоснование выбора метода решения задачи……………………
Решение задачи оптимизации…………………………………………….
АНАЛИЗ МОДЕЛИ НА ЧУВСТВИТЕЛЬНОСТЬ И УСТОЙЧИВОСТЬ………
Предварительный анализ оптимального решения………………...
Исследование чувствительности целевой функции………………..
Исследование устойчивости оптимального базисного плана……..
ПОИСК ОПТИМАЛЬНОГО ЦЕЛОЧИСЛЕННОГО РЕШЕНИЯ………………..
ЗАКЛЮЧЕНИЕ………………………………………………………………………….
СПИСОК ИСПОЛЬЗУЕМЫХ ИСТОЧНИКОВ…………………………………….
Введение
Исследование операций изучает хорошо структурированные задачи. Хорошо структурированная задача – задача, для которой существует математическая модель с одним критерием качества. Математическая модель – отражение интересующих нас свойств объектов или явлений с помощью математических символов.
Цель, которую преследуют в процессе исследования операций (ИО), заключается в том, чтобы выявить наилучший (оптимальный) способ действия при решении той или иной задачи организационного управления в условиях, когда имеют место ограничения технико-экономического или какого-либо другого характера. Когда используют термин исследование операций, то почти всегда имеют в виду применение математических методов для моделирования систем и анализа их характеристик. Действительно, математические модели и методы занимают в исследовании операций центральное место.
Исследование операций можно рассматривать и как науку, и как искусство. Правомерность утверждения о научности вытекает из того обстоятельства, что при решении возникающих проблем эффективно используются математические модели и методы. Исследование операций можно рассматривать и как искусство, поскольку успешное выполнение всех этапов исследования во многом определяется творческими способностями и интуицией исследователей. Поэтому при сборе информации, необходимой для построения математической модели, верификации модели и использовании получаемых с помощью модели результатов успех, несомненно, зависит от способности исследователей устанавливать рабочие контакты как с источниками необходимой информации, так и с лицами, ответственными за реализацию принимаемых решений.
Математическая модель – отображение оригинала в виде математических объектов, переменных, функций, уравнений и неравенств.
В простейшем случае математическая модель содержит три объекта:
1. совокупность неизвестных величин x
i
,
i
=1..
n
, значения которых необходимо найти. В совокупности эти значения образуют множество решений.
2. множество допустимых решений которые могут принимать x
i
,
i
=1..
n
. В хорошо формализованных задачах это множество может быть описано с помощью системы неравенств g
j
(
x
1
,
x
2
, .. ,
x
m
)<=0,
j
=1..
m
.
3. оценка качества допустимого решения. В простейшем случае оценка качества осуществляется с помощью специальной функции f(x1,x2,…,xn), которую называют целевой функцией.
В зависимости от ситуации необходимо найти такое допустимое решение, при котором целевая функция принимает максимальное или минимальное значение:
f
(
x
1,
x
2,…,
xn
)
-
extr
(
max
,
min
)
gi
(
x
1,
x
2,…,
xn
)<= 0,
i
=1..
n
,
где g – ограничение целевой функции.
Такая задача в которой множество решений может быть описано с помощью системы функциональных равенств и неравенств, называется задачей математического программирования
. .
Для решения задач этого типа разработан специальный раздел математики, который называется математическим программированием
.
Но не всегда удается сформулировать ограничение задачи в виде совокупности равенств или неравенств. Не всегда удается сформулировать цель в виде одной целевой функции, поэтому в общем случае математическая модель задачи может выглядеть более абстрактно:
Имеется множество допустимых решений X
.
Имеется совокупность критериев K
{
K
1
,
K
2
,...,
Kn
}
с помощью которых производится оценка качества решений X
. Требуется найти такое решение x
Î
X
, которое является наилучшим с точки зрения критериев K
1
,
K
2
,...,
Kn
. Такая задача и будет называться задачей исследования операций.
Классификация задач исследования операций.
Детерминированные задачи
— задачи, которые не содержат в себе элементов случайности.
Задачи дискретного программирования
отличаются от задач математического программирования тем, что в задачах дискретного программирования на значения переменных xi
накладывается требование дискретности, т.е. xi
может принимать не любые значения из диапазона (
a
,
b
),
а из множества конкретных фиксированных значений а1
, а2
, ..., а
k
.
Простейшим примером является требование целочисленности.
К задачам динамического программированию
относятся задачи, которые являются многошаговыми или могут быть сведены к многошаговым задачам. Примером многошаговых задач могут являться задачи планирования на несколько месяцев или лет.
Стохастические задачи
— задачи, которые содержат в своей постановке вероятностные элементы.
До недавнего времени теория очередей
сводилась всего к трем законам:
1).
Соседняя очередь всегда движется быстрее 2).
Как только вы переходите в соседнюю очередь, ваша прежняя начинает двигаться быстрее.
3).
Ваше метание из одной очереди в другую взвинчивает обе очереди. Теперь же теория очередей развилась в самостоятельный раздел математики — теорию массового обслуживания
. В математическом программировании все функции (и целевая, и функция ограничений) являются линейными.
Если хотя бы одна из функций нелинейная, то такая задача относится к нелинейному программированию
.
Выпуклые задачи
— это задачи, в которых целевая функция и функции ограничений обладают специальными свойствами выпуклости или вогнутости. Достоинством задач выпуклого программирования является то, что если решение существует, то оно единственно.
В выпуклом программировании выделяют специальный класс более простых задач, которые называются квадратичным программированием
. Эти задачи близки к линейным задачам. Как правило у них линейная система ограничений, но в целевую функцию могут входить вторые степени переменных, а так же их произведения.
1.
Постановка задачи оптимизации
2.
ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
Для перевозки грузов на трёх линиях могут быть использованы суда трёх типов. Необходимо определить какие суда и в течение какого следует использовать, чтобы обеспечить максимальную загрузку судов с учётом возможного времени их эксплуатации.
Пусть
-количество дней работы судов 1 типа на 1 линии
-количество дней работы судов 2 типа на 1 линии
-количество дней работы судов 3 типа на 1 линии
-количество дней работы судов 1 типа на 2 линии
-количество дней работы судов 2 типа на 2 линии
-количество дней работы судов 3 типа на 2 линии
-количество дней работы судов 1 типа на 3 линии
-количество дней работы судов 2 типа на 3 линии
-количество дней работы судов 3 типа на 3 линии
Общее время эксплуатации судов берётся из таблицы. Общее время эксплуатации должно быть не больше суммарного времени работы судов на всех трёх линиях. Отсюда следуют следующие ограничения:
(2.1)
Так же у каждого из судов есть производительность. Она отличается в зависимости от типов судов и от линий. Так же есть заданный объём перевозок(таблица 1). Запишем следующие ограничения:
(2.2)
Т.к. нам надо определить такое распределение судов при котором загрузка будет максимальной, то целевая функция будет на максимум.
Исходя из этих данных запишем целевую функцию:
L(x)= (2.3)
Математическая модель задачи будет выглядеть следующим образом
(2.4)
3. Обоснование выбора метода решения задачи
Любая задача линейного программирования может быть представлена в канонической и симплексной форме. В симплексной форме задача разрешена относительно базисных переменных.
Анализ задачи в симплексной форме позволяет сделать вывод: если коэффициенты отрицательные, то базисный план соответствует форме задачи линейного программирования, является наиболее оптимальным, поскольку любое изменение небазисных переменных (их увеличение) будет приводить к ухудшению целевой функции.
Достаточное условие оптимальности:
если в задаче на максимум коэффициенты целевой функции отрицательны, то соответствующий базисный план оптимален.
В задаче на минимум достаточное условие — неотрицательные коэффициенты целевой функции.
Если хотя бы один из коэффициентов целевой функции в задаче на минимум положителен, на максимум — отрицателен, то целевую функцию можно улучшить, увеличив значения соответствующих базисных переменных.
Для решения задач симплекс-методом необходимо:
a) привести задачу к канонической форме.
b) преобразовать задачу линейного программирования к симплексной форме, т.е. разделить переменные задачи на две группы: базисные
и небазисные.
Разрешить систему относительно базисных переменных и исключить базисные переменные из целевой функции.
c) Вспомнить одну или более итераций симплекс-метода, которые включают:
(1) проверку достаточных условий оптимальности и выбор ведущей переменной
(2) вычисление максимально допустимого шага вдоль ведущей переменной и выбор разрешающего уравнения
(3) преобразование симплексной формы, связанное с заменой в базисе — переменная, соответствующая разрешающему уравнению, покидает базис, на ее место вводится ведущая переменная.
Рассмотрим задачу линейного программирования в канонической форме:
(3.1)
Симплексной формой
задачи (3.1)
называется такое ее представление, при котором система основных ограничений разрешены относительно m
переменных, которые называются базисными
.
i
1
,
i
2
, …,
im
— индексы базисных переменных.
Тогда система основных ограничений разрешенная относительно базисных переменных будет иметь вид:
(3.2)
J
б
= {
i
1
,
i
2
, …,
im
}
— базисные переменные
J
н
= {
j
1
,
j
2
, …,
jm
}
— небазисные переменные
J
б
È
J
н
=
J
= {1, 2, …,
n
}
Используя уравнение системы (3.2)
мы можем исключить базисные переменные из целевой функции:
(3.3)
xj
³
0
, j
=
J
(3.4)
— известные ограничения.
Задачу линейного программирования с целевой функцией (3.3)
и ограничениями (3.2)
и (3.4)
будем называть задачей линейного программирования в симплексной форме
, если все свободные члены bi
³
0
, .
(3.5)
Такой задаче ставим в соответствие базисный план, который состоит из:
Введем дополнительное обозначение:
Задачу (5)
можно переписать в следующем виде:
(3.6)
1. Симплекс-метод в табличной форме.
Предположим, что задача линейного программирования в канонической форме приведена к симплексной форме (3.6)
. Перенесем коэффициенты целевой функции вектора свободных членов и матрицы R
,
L
,
в специальную симплексную таблицу:
| бн
|
B
|
|
…
|
|
…
|
|
…
|
|
|
|
|
|
r
|
…
|
r
|
…
|
r
|
…
|
r
|
|
|
|
|
|
…
|
|
…
|
|
…
|
|
|
| …
|
…
|
…
|
…
|
…
|
…
|
…
|
…
|
||
|
|
|
|
…
|
|
…
|
|
…
|
|
|
| …
|
…
|
…
|
…
|
…
|
…
|
…
|
…
|
||
|
|
|
|
…
|
|
…
|
|
…
|
|
|
| …
|
…
|
…
|
…
|
…
|
…
|
…
|
…
|
||
|
|
|
|
|
|
|
1 этап:
проверяем достаточное условие оптимальности (таблица) базисного плана.
Если в задаче на максимум коэффициенты целевой функции удовлетворяет условию:
в задаче на минимум:
, то базисный план, соответствующий симплексной таблице является оптимальным и решение ЗПП окончено.
Допустим, что достаточное условие оптимальности не выполняется. В этом случае найдем максимальный по модулю, не удовлетворяющий условию оптимальности, коэффициент целевой функции. Соответствующую переменную назовем ведущей
. Этот индекс обозначен через j
*
.
Соответствующий столбец симплексной таблицы назовем ведущим
.
2 этап
: вычисление оптимально допустимого шага.
Из системы основных ограничений симплексной формы следует:
Из
вытекает, что
(3.7)
Очевидно, что соотношение (7)
может нарушаться только в том случае, если , поэтому для каждой базисной переменной :
, то соответствующая базисная переменная ни при каком значении не обращается в ноль.
Максимально допустимый шаг определим как минимум из чисел .
возможны 2 случая:
1.
Это означает, что шаг может быть бесконечно большим. Соответственно значение целевой функции может неограниченно возрастать, если задача на максимум, или неограниченно убывать, если задача на минимум. С экономической точки зрения такая ситуация свидетельствует о несоответствии математической модели реальной экономической модели (неучтены или неправильно учтены ограничения задачи). Решение ЗПП при этом считается законченным.
2. (конечное число)
называется ведущей
или разрешающей
переменной.
3 этап:
преобразование симплексной таблицы.
Как уже отмечалось в пункте 1, должна выводиться из базиса, на ее место вводится переменная — небазисная.
Нарисуем новую симплексную таблицу.
| бн
|
b
|
|
…
|
|
…
|
|
|
|
|
|
…
|
R
|
…
|
r
|
|
|
|
|
…
|
|
…
|
|
| …
|
…
|
…
|
…
|
…
|
…
|
|
|
|
|
|
…
|
|
…
|
|
| …
|
…
|
…
|
…
|
…
|
…
|
|
|
|
|
|
|
|
(
3.8
)
(
3.9
)
(
3.10
)
(
3.11
)
Правила:
1)
На место ведущего элемента в новой таблице записывается элемент, обратный ведущему.
2)
,
Правило: новые элементы ведущей строки вычисляются по старым элементам делением последних на ведущий элемент, взятый с противоположным знаком.
i
— номер произвольной строки.
3)
Элементы ведущего столбца получаются из старых элементов делением последних на ведущий элемент.
4)
,
,
| j
|
j
|
||
| i
|
rij
|
…
|
rij*
|
| …
|
…
|
…
|
|
| i
|
ri
|
…
|
ri
|
Новый элемент симплексной таблицы, не находящийся в ведущей строке или столбце, получается из старого, если вычесть из последнего произведение его проекций деленное на ведущий элемент.
Двухфазный симплекс-метод.
Если после приведения математической модели к канонической форме не все ограничения находятся в предпочтительном виде, то для решения задачи требуется применить двухфазный симплекс-метод.
Целевая функция вспомогательной задачи представляет собой сумму фиктивных переменных, которую необходимо минимизировать.
Начальный базис вспомогательной задачи составляют переменные:
1. в ограничениях, которые находятся в предпочтительном виде, -те переменные, которые обеспечили предпочтительность этих ограничений.
2. в ограничениях, которые не находятся в предпочтительном виде,- фиктивные переменные
Целевую функцию вспомогательной задачи выражаем через небазисные переменные и решаем систему обычным симплекс-методом(описан выше).
Если все искусственные переменные покинули базис и исключены из таблицы, а строка вспомогательной функции содержит одни нули, то это означает, что решение вспомогательной задачи окончено.
Т.к. как все фиктивные переменные покинули задачу, то система основных ограничений вспомогательной задачи совпадает с исходной системой основных ограничений. Искуственнная целевая функция с исчезновением фиктивных переменных трансформировалась в нулевую. Поэтому для перехода ко второй фазе симплекс-метода необходимо вернуться к исходной целевой функции. Поскольку исходная целевая функция содержит базисные переменные, их необходимо исключить, выразив через небазисные переменные. Далее система решается обычным симплекс-методом, описанным раннее.
4. Решение задачи оптимизации
Для решения задачи необходимо привести её математическую модель(2.4) к канонической форме(введём свободные переменные).
Т.к. не все ограничения находятся в предпочтительном виде, значит будем применять при решении задачи двухфазный симплекс-метод.
Введём фиктивные переменные х16, х17, х18 и решим задачу относительно их:
Будем решать задачу табличным симплекс-методом, как описывалось ранее.
| min |
b |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
x8 |
x9 |
x10 |
x11 |
x12 |
| L |
11700 |
-8 |
-6 |
-12 |
-14 |
-15 |
-12 |
-11 |
-13 |
-4 |
1 |
1 |
1 |
| x16 |
3000 |
-8 |
-6 |
-12 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
| x17 |
5400 |
0 |
0 |
0 |
-14 |
-15 |
-12 |
0 |
0 |
0 |
0 |
1 |
0 |
| x18 |
3300 |
-1 |
0 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
0 |
0 |
| x13 |
300 |
-1 |
0 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
0 |
0 |
| x14 |
300 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
0 |
| x15 |
300 |
0 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
Этой симплексной таблице соответствует начальный базисный план х(0)
= (0,0,0,0,0,0,0,0,0,0,0,0,300,300,300,3000,5400,3300). L (x(
0
)
) = 11700. Этот план не оптимален, так как целевая функция вспомогательной задачи на минимум содержит отрицательные коэффициенты.
Выбираем базисный столбец и разрешающую строку:
| min |
b |
x1 |
x2 |
x3 |
x4 |
x5
|
x6 |
x7 |
x8 |
x9 |
x10 |
x11 |
x12 |
θ |
| L |
11700 |
-8 |
-6 |
-12 |
-14 |
-15
|
-12 |
-11 |
-13 |
-4 |
1 |
1 |
1 |
|
| x16 |
3000 |
-8 |
-6 |
-12 |
0 |
0
|
0 |
0 |
0 |
0 |
1 |
0 |
0 |
∞ |
| x17 |
5400 |
0 |
0 |
0 |
-14 |
-15
|
-12 |
0 |
0 |
0 |
0 |
1 |
0 |
360 |
| x18 |
3300 |
-1 |
0 |
0 |
-1 |
0
|
0 |
-1 |
0 |
0 |
0 |
0 |
0 |
∞ |
| x13 |
300 |
-1 |
0 |
0 |
-1 |
0
|
0 |
-1 |
0 |
0 |
0 |
0 |
0 |
∞ |
| x14
|
300
|
0
|
-1
|
0
|
0
|
-1
|
0
|
0
|
-1
|
0
|
0
|
0
|
0
|
300
|
| x15 |
300 |
0 |
0 |
-1 |
0 |
0
|
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
∞ |
Пересчёт симплексной таблицы происходит по правилам, описанным в разделе 3.
Произведём пересчёт симплексной таблицы и если план получится не оптимальным, то снова выберем ведущий столбец и разрешающую строку.
| min |
b |
x1 |
x2 |
x3 |
x
|
x14 |
x6 |
x7 |
x8 |
x9 |
x10 |
x11 |
x12 |
θ |
| L |
7200 |
-8 |
-6 |
-12 |
-14
|
15 |
-12 |
- 11 |
2 |
-4 |
1 |
1 |
1 |
|
| x16 |
3000 |
-8 |
-6 |
-12 |
0
|
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
∞ |
| x17
|
900
|
0
|
15
|
0
|
-14
|
15
|
-12
|
0
|
15
|
0
|
0
|
1
|
0
|
64.286
|
| x18 |
3300 |
0 |
0 |
0 |
0
|
0 |
0 |
-11 |
-13 |
-4 |
0 |
0 |
1 |
∞ |
| x13 |
300 |
-1 |
0 |
0 |
-1
|
0 |
0 |
-1 |
0 |
0 |
0 |
0 |
0 |
∞ |
| x5 |
300 |
0 |
-1 |
0 |
0
|
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
0 |
300 |
| x15 |
300 |
0 |
0 |
-1 |
0
|
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
∞ |
Произведём пересчёт симплексной таблицы и если план получится не оптимальным, то снова выберем ведущий столбец и разрешающую строку.
| min |
b |
x1 |
x2 |
x3 |
x14 |
x6 |
x7 |
x
|
x9 |
x10 |
x11 |
x12 |
θ |
| L |
6300 |
-8 |
-6 |
-12 |
0 |
0 |
- 11 |
-13
|
-4 |
1 |
0 |
1 |
|
| x16 |
3000 |
-8 |
-6 |
-12 |
0 |
0 |
0 |
0
|
0 |
1 |
0 |
0 |
∞ |
| x4 |
64.286 |
0 |
1.071 |
0 |
1.071 |
-0.857 |
0 |
1.071
|
0 |
0 |
0.071 |
0 |
∞ |
| x18 |
3300 |
0 |
0 |
0 |
0 |
0 |
-11 |
-13
|
-4 |
0 |
0 |
1 |
253.846 |
| x13
|
235.174
|
-1
|
0
|
0
|
-1
|
0
|
-1
|
-1.071
|
0
|
0
|
-0.071
|
0
|
220
|
| x5 |
300 |
0 |
-1 |
0 |
-1 |
0 |
0 |
-1
|
0 |
0 |
0 |
0 |
300 |
| x15 |
300 |
0 |
0 |
-1 |
0 |
-1 |
0 |
0
|
-1 |
0 |
0 |
0 |
∞ |
Произведём пересчёт симплексной таблицы и если план получится не оптимальным, то снова выберем ведущий столбец и разрешающую строку.
| min |
b |
x1 |
x2 |
x3 |
x14 |
x6 |
x7 |
x13 |
x9 |
x10 |
x11 |
x12 |
θ |
| L |
3440 |
4.133 |
7 |
-12 |
13 |
-10.4 |
1.133 |
12.133 |
-4 |
1 |
0.867 |
1 |
|
| x16 |
3000 |
-8 |
-6 |
-12 |
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
250 |
| x4 |
300 |
-1 |
0 |
0 |
0 |
0 |
-1 |
-1 |
0 |
0 |
0 |
0 |
∞ |
| x18 |
440 |
12.133 |
13 |
0 |
13 |
-10.4 |
1.133 |
12.133 |
-4 |
0 |
0.867 |
1 |
∞ |
| x8 |
220 |
-0.933 |
-1 |
0 |
-1 |
0.8 |
-0.933 |
-.933 |
0 |
0 |
-0.067 |
0 |
∞ |
| x5 |
80 |
0.933 |
0 |
0 |
0 |
-0.8 |
0.933 |
0.933 |
0 |
0 |
0.067 |
0 |
∞ |
| x15 |
300 |
0 |
0 |
-1 |
0 |
-1 |
0 |
0 |
-1 |
0 |
0 |
0 |
300 |
Произведём пересчёт симплексной таблицы и если план получится не оптимальным, то снова выберем ведущий столбец и разрешающую строку.
| min |
b |
x1 |
x2 |
x14 |
x6
|
x7 |
x13 |
x9 |
x10 |
x11 |
x12 |
θ |
| L |
440 |
12.133 |
13 |
13 |
-10.4
|
1.133 |
12.133 |
-4 |
0 |
0.867 |
1 |
|
| x3 |
250 |
-0.667 |
-0.5 |
0 |
0
|
0 |
0 |
0 |
0.083 |
0 |
0 |
∞ |
| x4 |
300 |
-1 |
0 |
0 |
0
|
-1 |
-1 |
0 |
0 |
0 |
0 |
∞ |
| x18
|
440
|
12.133
|
13
|
13
|
-10.4
|
1.133
|
12.133
|
-4
|
0
|
0.867
|
1
|
42.308
|
| x8 |
220 |
-0.933 |
-1 |
-1 |
0.8
|
-0.933 |
0.933 |
0 |
0 |
-0.067 |
0 |
∞ |
| x5 |
80 |
0.933 |
0 |
0 |
-0.8
|
0.933 |
0.933 |
0 |
0 |
0.067 |
0 |
100 |
| x15 |
50 |
0.667 |
0.5 |
0 |
-1
|
0 |
0 |
-1 |
-0.083 |
0 |
0 |
50 |
Произведём пересчёт симплексной таблицы
| min |
b |
x1 |
x2 |
x14 |
x7 |
x13 |
x9 |
x10 |
x11 |
x12 |
| L |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
| x3 |
250 |
-0.667 |
-0.5 |
0 |
0 |
0 |
0 |
0.083 |
0 |
0 |
| x4 |
300 |
-1 |
0 |
0 |
-1 |
-1 |
0 |
0 |
0 |
0 |
| x6 |
42.308 |
1.167 |
1.25 |
1.25 |
0.109 |
1.167 |
-0.385 |
0 |
0.083 |
0.096 |
| x8 |
253.846 |
0 |
0 |
0 |
-0.846 |
0 |
-0.308 |
0 |
0 |
0.077 |
| x5 |
46.154 |
0 |
-1 |
-1 |
0.846 |
0 |
0.308 |
0 |
0 |
-0.077 |
| x15 |
7.692 |
-0.5 |
-0.75 |
-1.25 |
-0.109 |
-1.167 |
-0.615 |
-0.083 |
-0.083 |
-0.096 |
Как видно, все искусственные переменные покинули базис и исключены из таблицы. Первая фаза решения завершена, так как все коэффициенты целевой функции равны нулю. Перейдём ко второй фазе(для этого вернёмся к исходной целевой функции):
| max |
b |
x1 |
x2 |
x14 |
x7 |
x13 |
x9 |
x10 |
x11 |
x12 |
| L |
11700 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
1 |
1 |
| x3 |
250 |
-0.667 |
-0.5 |
0 |
0 |
0 |
0 |
0.083 |
0 |
0 |
| x4 |
300 |
-1 |
0 |
0 |
-1 |
-1 |
0 |
0 |
0 |
0 |
| x6 |
42.308 |
1.167 |
1.25 |
1.25 |
0.109 |
1.167 |
-0.385 |
0 |
0.083 |
0.096 |
| x8 |
253.846 |
0 |
0 |
0 |
-0.846 |
0 |
-0.308 |
0 |
0 |
0.077 |
| x5 |
46.154 |
0 |
-1 |
-1 |
0.846 |
0 |
0.308 |
0 |
0 |
-0.077 |
| x15 |
7.692 |
-0.5 |
-0.75 |
-1.25 |
-0.109 |
-1.167 |
-0.615 |
-0.083 |
-0.083 |
-0.096 |
Базисный план не оптимален, так как целевая функция задачи на максимум содержит положительные коэффициенты. Выбираем базисный столбец и разрешающую строку:
| max |
b |
x1 |
x2 |
x14 |
x7 |
x13 |
x9 |
x10
|
x11 |
x12 |
θ |
| L |
11700 |
0 |
0 |
0 |
0 |
0 |
0 |
1
|
1 |
1 |
|
| x3 |
250 |
-0.667 |
-0.5 |
0 |
0 |
0 |
0 |
0.083
|
0 |
0 |
∞ |
| x4 |
300 |
-1 |
0 |
0 |
-1 |
-1 |
0 |
0
|
0 |
0 |
∞ |
| x6 |
42.308 |
1.167 |
1.25 |
1.25 |
0.109 |
1.167 |
-0.385 |
0
|
0.083 |
0.096 |
∞ |
| x8 |
253.846 |
0 |
0 |
0 |
-0.846 |
0 |
-0.308 |
0
|
0 |
0.077 |
∞ |
| x5 |
46.154 |
0 |
-1 |
-1 |
0.846 |
0 |
0.308 |
0
|
0 |
-0.077 |
∞ |
| x15
|
7.692
|
-0.5
|
-0.75
|
-1.25
|
-0.109
|
-1.167
|
-0.615
|
-0.083
|
-0.083
|
-0.096
|
92.308
|
Произведём пересчёт симплексной таблицы:
| max |
b |
x1 |
x2 |
x14 |
x7 |
x13 |
x9 |
x15 |
x11 |
x12 |
| L |
11792.308 |
-6 |
-9 |
-15 |
-1.308 |
-14 |
-7.385 |
-12 |
0 |
-0.154 |
| x3 |
257.692 |
-1.167 |
-1.25 |
-1.25 |
-1.109 |
-1.167 |
-0.615 |
-1 |
-0.083 |
-0.096 |
| x4 |
300 |
-1 |
0 |
0 |
-1 |
-1 |
0 |
0 |
0 |
0 |
| x6 |
42.308 |
1.167 |
1.25 |
1.25 |
0.109 |
1.167 |
-0.385 |
0 |
0.083 |
0.096 |
| x8 |
253.846 |
0 |
0 |
0 |
-0.846 |
0 |
-0.308 |
0 |
0 |
0.077 |
| x5 |
46.154 |
0 |
-1 |
-1 |
0.846 |
0 |
0.308 |
0 |
0 |
-0.077 |
| x10 |
92.308 |
-6 |
-9 |
-15 |
-1.308 |
-14 |
-7.385 |
-12 |
-1 |
-1.154 |
Поскольку строка L таблицы не содержит положительных элементов, то базисный план х(1)
= (0,0,257.693,300,46.154,42.308,0,253.846,0,92.308,0,0,0,0,0). L (x(2)
) = 11792.308 является оптимальным. Решение задачи окончено.
Значения основных переменных задачи х1
, х2
,х3
, х4
, х5,
х6,
х7,
х8,
х9,
обозначают, что:
-суда 1 типа должны работать только на 2 линии
-суда 2 типа должны работать 46.15385 дней на 2 линии и 253.84615 дней на 3 линии
-суда 3 типа должны работать 257.69231 дней на 1 линии и 42.30769 дней на второй линии
При этом максимальная загрузка суден составит 11792.308 млн. тонно-миль.
5.
АНАЛИЗ МОДЕЛИ НА ЧУВСТВИТЕЛЬНОСТЬ И УСТОЙЧИВОСТЬ
5.1. Предварительный анализ оптимального решения
Первые три ограничения являются количественными, т.к. исходят из обязательного заданного объёма перевозок.
Остальные три ограничения являются ресурсными, т.к. исходят из общего времени эксплуатации судов. Очевидно, что при увеличении времени эксплуатации максимальная загрузка судов будет увеличиваться
5.2. Исследование чувствительности целевой функции
Для исследования чувствительности целевой функции к изменениям правых частей основных ограничений необходимо найти оптимальный двойственный план:
(5.1)
и исследовать его компоненты. Согласно полученному решению базисными переменными являются х3
, х4
, х5
,х6
,х8
,х10
.
Вычислим обратную матрицу:
Найдём оптимальный двойственный план:
Дадим экономическую оценку полученным результатам:
1) Увеличение заданного объёма перевозок на 1 линии к увеличению суммарной максимальной загрузки судов не приведёт
2) Увеличение заданного объёма перевозок на 2 линии к увеличению суммарной максимальной загрузки судов не приведёт
3) Увеличение заданного объёма перевозок на 3 линии на 1 единицу приведёт к увеличению суммарной максимальной загрузки судов на 1
4) Увеличение общего времени эксплуатации судов на 1 линии на 1 сутки приведёт к увеличению максимальной загрузки судов на 14 млн.тнно-миль.
5) Увеличение общего времени эксплуатации судов на 2 линии на 1 сутки приведёт к увеличению максимальной загрузки судов на 15 млн.тнно-миль.
6) Увеличение общего времени эксплуатации судов на 3 линии на 1 сутки приведёт к увеличению максимальной загрузки судов на 12 млн.тнно-миль.
5.3. Исследование устойчивости оптимального базисного плана
Определим интервал устойчивости
, (5.2)
где - максимально возможное уменьшение, а - увеличение правой части i-го ограничения.
Найдём интервал устойчивости для первого ограничения:
Таким образом, интервал устойчивости для b1 составляет:
Найдём интервал устойчивости для второго ограничения:
Таким образом, интервал устойчивости для b2 составляет:
Найдём интервал устойчивости для третьего ограничения:
Таким образом, интервал устойчивости для b3 составляет:
Найдём интервал устойчивости для четвёртого ограничения:
Таким образом, интервал устойчивости для b4 составляет:
Найдём интервал устойчивости для пятого ограничения:
Таким образом, интервал устойчивости для b5 составляет:
Найдём интервал устойчивости для шестого ограничения:
Таким образом, интервал устойчивости для b6 составляет:
6.
ПОИСК ОПТИМАЛЬНОГО ЦЕЛОЧИСЛЕННОГО РЕШЕНИЯ
Для поиска оптимального целочисленного решения воспользуемся методом ветвей и границ. Для этого нецелочисленную компоненту оптимального решения x3=257.692. Исходя из этого решаем 2 задачи:
1) x3≤[257.692]=257
2) x3≥[257.692]+1=258
Решив симплекс-методом задачи с дополнительными ограничениями (1) и (2), получим одно решения:
L(x2*
)=11792.308
При ограничении x3≥258 решения не существует.
x2*
не является целочисленным поэтому процесс решения продолжается дальше. Это решение содержит одну нецелочисленную компоненту х5.
Исходя из этого решаем 2 задачи:
1) x4≤[ 46.154]=46
2) x4≥[46.154]+1=47
Решив симплекс-методом задачи с дополнительными ограничениями (1) и (2), получим 2 решения:
L(x3*
)=11792
L(x4*
)= 11791
Найдено оптимальное целочисленное решение – это решение задачи 3.
Поиск решения можно изобразить графически:
рисунок
Значения основных переменных задачи х1
, х2
,х3
, х4
, х5,
х6,
х7,
х8,
х9,
обозначают, что:
-суда 1 типа должны работать 299 дней на 2 линии и 1 день на 3 линии
-суда 2 типа должны работать 47 дней на 2 линии и 253 дня на 3 линии
-суда 3 типа должны работать 250 дней на 1 линии и 50 дней на второй линии
При этом максимальная загрузка суден составит 11792 млн. тонно-миль
ЗАКЛЮЧЕНИЕ
В процессе работы над курсовой работой была успешно решена задача по оптимальному распределению судов 3 типов по различным линиям с максимальной загрузкой судов. С помощью симплекс-метода был найден оптимальный базисный план задачи, проведёно исследование его целевой функции на чувствительность и определены интервалы устойчивости, а так же было найдено оптимальное целочисленное решение.
СПИСОК ИСПОЛЬЗУЕМЫХ ИСТОЧНИКОВ
1. Ракецкий В.М. Методические указания и задания к курсовой работе “Решение задач линейного программирования” по дисциплине “Системный анализ и исследование операций”
Брест: БГТУ 2007
2. Акулич И.Л. Математическое программирование в примерах и задачах. - М.: Высшая школа, 1986.-319с.
3. Альсевич В.В., Габасов Р., Глушенков В.С. Оптимизация линейных экономических моделей: Статистические задачи. – Мн.: БГУ, 2000. -210с.
4. Балашевич В.А. Основы математического программирования. – Мн.: Вышэйшая школа, 1985.-173с.
5. Банди Б. Основы линейного программирования. – М.: Радио и связь, 1989. – 176 с.