Построение математической модели и разработка программного обеспечения для решения задачи организационного управления. 2
Государственное образовательное учреждение
профессионального образования
«ТЮМЕНСКИЙ ГОСУДАРСТВЕННЫЙ НЕФТЕГАЗОВЫЙ УНИВЕРСИТЕТ»
КАФЕДРА ИНФОРМАТИКИ И ВЫЧИСЛИТЕЛЬНОЙ ТЕХНИКИ
КУРСОВАЯ РАБОТА
по курсу «Теория принятия решений»
Тема: «Построение математической модели и разработка программного обеспечения для решения задачи организационного управления»
Вариант
№29
Тюмень 2010
Содержание
Постановка задачи
На
производство поступила партия стержней
длиной 250 и 190 см, причем количество стержней
длиной 250 см ограничено и равно 200 шт. Из
этих стержней необходимо получить 470
заготовок длиной 120 см и 450 заготовок длиной
80 см. Как следует разрезать стержни, чтобы
количество отходов было минимальным?
Построение математической модели
Математическая формулировка
j —индекс вида стержня, j = 1,..., n;
k —индекс вида заготовки, k = 1, ..., q;
i — индекс способа раскроя стержня, i = 1,..., p;
аijk — количество (целое число) заготовок вида k, полученных при раскрое j-го типа стержня i-м способом;
bk — число заготовок вида k которое требуется получить;
dj — количество стержней j-го вида;
xij — количество стержней j-го типа, раскраиваемых по i-му способу (интенсивность использования способа раскроя);
cij
— величина отхода, полученного при
раскрое стержня j-го типа по i-му
способу;
Целевая функция. Необходимо минимизировать отходы, полученые при раскрое стержней.
min
F = xji cji, j=1...n,
k=1...q
Ограничения. При решении рассматриваемой задачи должны быть учтены ограничения на расход исходных стержней имеющихся в наличии и количество заготовок необходимых для выполнения заказа.
Ограничение на расход исходных стержней можно записать в следующем виде:
Это приводит к следующему ограничению:
(x11i)
≤ 200;
Ограничения на количество заготовок необходимых для выполнения заказа принимает следующий математический вид
(аij1 xij) = 470;
(аij2
xij) = 450;
Неявное ограничение заключается в том, что количества затраченных стержней и произведенных заготовок не могут принимать отрицательных значений. Чтобы предотвратить получение таких недопустимых решений, потребуем выполнения условия неотрицательности переменных, т.е. введем ограничения:
xij
≥ 0, аijk
≥ 0;
Итак, математическую модель можно записать следующим образом.
Определить интенсивности использования для каждого метода раскроя, при которых достигается минимальное количество отходов:
min
F =
xji cji, j=1...n,
k=1...q (целевая функция)
и
которые удовлетворяют условиям
(xi1) ≤ d1;
(аij1 xij) = b1;
(аij2
xij) = b2;
xij
≥ 0, аijk
≥ 0;
Выбор и обоснование метода решений поставленной задачи.
Т.к. все входящие в модель функции (ограничения и целевая функция) являются линейными, то данная задача относится к классу задач линейного программирования (ЛП), поэтому для ее решения необходимо применить один из методов решения задач ЛП. Универсальный метод решения таких задач – симплекс-метод.
Решение задачи
Заполним
таблицу методов раскроя
a для стержней типа j = 1 (250 см)
| Вариант раскроя i | Количество полученных стержней | Полученный
отход c | |
| тип k = 1 (120 см) | тип k = 2 (80см) | ||
| 1 | 2 | 0 | 10 |
| 2 | 1 | 1 | 50 |
| 3 | 0 | 3 | 10 |
a для стержней типа j = 2 (190 см)
| Вариант раскроя i | Количество полученных стержней | Полученный
отход c | |
| тип k = 1(120 см) | тип k = 2(80см) | ||
| 1 | 1 | 0 | 70 |
| 2 | 0 | 2 | 30 |
развернем целевую
функцию в соответствии с таблицей раскроя:
min F = 10x11
+ 50x21 + 10x31 + 70x12 + 30x22
развернем ограничения
в соответствии с таблицей раскроя:
x11 + x21 + x31 ≤ 200;
2x11 + x21
+ x12 = 470
x21 + 3x31
+ 2x22 = 450
xji
≥ 0, аijk
≥ 0;
Решим сформулированную задачу с помощью симплекс-метода. Для этого запишем целевую функцию в виде: Z - 10x11 - 50x21 - 10x31 - 70x12 - 30x22 = 0 и запишем исходные данные в симплекс-таблицу. Процесс нахождения оптимального решения приведен в таблице 2.
Таблица 2
Начальная симплекс таблица задачи на минимум, приведенной к каноническому виду.
Возможно улучшение плана.
Разрешающий столбец определяется по двум последним строкам таблицы.
В пересечении колонок Х0-Х8(выбирается максимальное положительное число).
Сначала просматривается строка помеченная знаком "M-->"
и если в ней нет положительныхчисел, просматривается последняя строка.
Если разрешающий столбец не нашли, то в таблице представлен оптимальный план.
Разрешающая строка определяется по минимальному не отрицательному отношению
коэффициентов столбца
Х0 и разрешающего столбца(что представлено
в столбце Alfa).
| F(Min) | 10 | 50 | 10 | 70 | 30 | 0 | M | M | |||
| Сi | P0 | X0 | X11 | X21 | X31 | X12 | X22 | X6 | X7 | X8 | Alfa |
| 0 | 6 | 200,00 | 1,00 | 1,00 | 1,00 | 0,00 | 0,00 | 1,00 | 0,00 | 0,00 | 200 |
| M | 7 | 470,00 | 2,00 | 1,00 | 0,00 | 1,00 | 0,00 | 0,00 | 1,00 | 0,00 | -1 |
| M | 8 | 450,00 | 0,00 | 1,00 | 3,00 | 0,00 | 2,00 | 0,00 | 0,00 | 1,00 | 150 |
| M--> | 920 | 2 | 2 | 3 | 1 | 2 | 0 | 0 | 0 | ||
| 0,00 | -10 | -50 | -10 | -70 | -30 | 0 | 0 | 0 |
| После итерации №1(возможно улучшение плана) | |||||||||||
|
Так как в предпоследнейстроке( |
|||||||||||
| (имеется максимальное положительное число). | |||||||||||
| А в столбце "Alfa" имеется не отрицательное число(выбрано минимальное) | |||||||||||
| F(Min) | 10 | 50 | 10 | 70 | 30 | 0 | M | M | |||
| Сi | P1 | X0 | X11 | X21 | X31 | X12 | X22 | X6 | X7 | X8 | Alfa |
| 0 | 6 | 50,00 | 1,00 | 0,67 | 0,00 | 0,00 | -0,67 | 1,00 | 0,00 | -0,33 | 50 |
| M | 7 | 470,00 | 2,00 | 1,00 | 0,00 | 1,00 | 0,00 | 0,00 | 1,00 | 0,00 | 235 |
| 10 | 3 | 150,00 | 0,00 | 0,33 | 1,00 | 0,00 | 0,67 | 0,00 | 0,00 | 0,33 | -1 |
| M--> | 470 | 2 | 1 | 0 | 1 | 0 | 0 | 0 | -1 | ||
| 1500,00 | -10 | -46,7 | 0 | -70 | -23,3 | 0 | 0 | 3,33 | |||
| После итерации №2(возможно улучшение плана) | |||||||||||
|
Так как в предпоследнейстроке( |
|||||||||||
| (имеется
максимальное положительное |
|||||||||||
| А в столбце "Alfa" имеется не отрицательное число(выбрано минимальное) | |||||||||||
| F(Min) | 10 | 50 | 10 | 70 | 30 | 0 | M | M | |||
| Сi | P2 | X0 | X11 | X21 | X31 | X12 | X22 | X6 | X7 | X8 | Alfa |
| 10 | 1 | 50,00 | 1,00 | 0,67 | 0,00 | 0,00 | -0,67 | 1,00 | 0,00 | -0,33 | -75 |
| M | 7 | 370,00 | 0,00 | -0,33 | 0,00 | 1,00 | 1,33 | -2,00 | 1,00 | 0,67 | 278 |
| 10 | 3 | 150,00 | 0,00 | 0,33 | 1,00 | 0,00 | 0,67 | 0,00 | 0,00 | 0,33 | 225 |
| M--> | 370 | 0 | -0,33 | 0 | 1 | 1,33 | -2 | 0 | -0,33 | ||
| 2000,00 | 0 | -40 | 0 | -70 | -30 | 10 | 0 | 0 | |||
| После итерации №3(возможно улучшение плана) | |||||||||||
|
Так как в предпоследнейстроке( |
|||||||||||
| (имеется
максимальное положительное |
|||||||||||
| А в столбце "Alfa" имеется не отрицательное число(выбрано минимальное) | |||||||||||
| F(Min) | 10 | 50 | 10 | 70 | 30 | 0 | M | M | |||
| Сi | P3 | X0 | X11 | X21 | X31 | X12 | X22 | X6 | X7 | X8 | Alfa |
| 10 | 1 | 200,00 | 1,00 | 1,00 | 1,00 | 0,00 | 0,00 | 1,00 | 0,00 | 0,00 | -1 |
| M | 7 | 70,00 | 0,00 | -1,00 | -2,00 | 1,00 | 0,00 | -2,00 | 1,00 | 0,00 | 70 |
| 30 | 5 | 225,00 | 0,00 | 0,50 | 1,50 | 0,00 | 1,00 | 0,00 | 0,00 | 0,50 | -1 |
| M--> | 70 | 0 | -1 | -2 | 1 | 0 | -2 | 0 | -1 | ||
| 8750,00 | 0 | -25 | 45 | -70 | 0 | 10 | 0 | 15 | |||
| После итерации №4 получили оптимальный план. | |||||||||||
|
Для функции 10x11+50x21+10x31+70x12+30x22= |
|||||||||||
| Достигается при X11=200;X12=70;X22=225; | |||||||||||
| Номера переменных Х и их значения находятся, соответственно, во второй и третьей колонках симплекс таблицы | |||||||||||
| Оптимальный план находится в первой ячейке последней строки симплекс таблицы | |||||||||||
| F(Min) | 10 | 50 | 10 | 70 | 30 | 0 | M | M | |||
| Сi | P4 | X0 | X11 | X21 | X31 | X12 | X22 | X6 | X7 | X8 | Alfa |
| 10 | 1 | 200,00 | 1,00 | 1,00 | 1,00 | 0,00 | 0,00 | 1,00 | 0,00 | 0,00 | -1 |
| 70 | 4 | 70,00 | 0,00 | -1,00 | -2,00 | 1,00 | 0,00 | -2,00 | 1,00 | 0,00 | -1 |
| 30 | 5 | 225,00 | 0,00 | 0,50 | 1,50 | 0,00 | 1,00 | 0,00 | 0,00 | 0,50 | -1 |
| M--> | 0 | 0 | 0 | 0 | 0 | 0 | 0 | -1 | -1 | ||
| 13650,00 | 0 | -95 | -95 | 0 | 0 | -130 | 70 | 15 | |||
Выполнив четыре итерации для получения оптимального решения, получили результирующую симплекс-таблицу, из которой следует, что оптимальное решение имеет вид:
x11 = 200, x12 = 70, x22 = 225.
Анализ модели на чувствительность.
Проведем анализ полученного решения. Результирующая симплекс-таблица «насыщена» весьма важными данными, лишь небольшую часть которых составляют оптимальные значения переменных. Из симплекс-таблицы непосредственно, либо при помощи простых дополнительных вычислений можно получить информацию относительно
- оптимального решения,
- статуса ресурсов,
- ценности каждого ресурса,
- чувствительности оптимального решения к изменению запасов ресурсов, вариациям коэффициентов целевой функции.
Оптимальное решение
Используя данные, содержащиеся в симплекс-таблице для оптимального решения, основные результаты можно представить так:
Таблица
3
| Управляемые
переменные |
Оптимальное значение | Решение |
| x11 | 200 | 200 стержней первого типа (250 см) должны быть порезаны на заготовки первого типа (120 см) |
| x12 | 70 | 70 стержней второго типа (190 см) должны быть порезаны на заготовки первого типа (120 см) |
| x22 | 225 | 225 стержней второго типа (190 см) должны быть порезаны на заготовки второго типа (80 см) |
Статус ресурсов
Виды стержней могут быть разделены на дефицитные и недефицитные в зависимости от того, полное или частичное их использование предусматривает оптимальное решение задачи.
Применительно к рассматриваемой задаче можно привести следующую сводку результатов.
Таблица 4
| Ресурс | Остаток | Статус ресурса |
| Стержни 250 см. | 0 | Дефицитный |
| Стержни 190 см. | ∞ | Недефицитный |
Положительное значение остаточной переменной указывает на неполное использование соответствующего ресурса, т.е. данный ресурс не является дефицитным. Если же остаточная переменная равна нулю, то это свидетельствует о полном потреблении соответствующего ресурса.
Ресурс,
увеличение запасов которых позволяет
улучшить решение (увеличить доход) –
это стержни длинной 250 см. т.к. они дефицитные.
Ценность ресурсов
Ценность
ресурса характеризуется
Максимальное изменение запаса ресурса
При
решении вопроса о том, запас
какого из ресурсов следует увеличить
в первую очередь, обычно используются
теневые цены. Чтобы определить интервал
значений изменения запаса ресурса, при
которых теневая цена данного ресурса,
фигурирующая в последней симплекс-таблице,
остается неизменной, необходимо выполнить
ряд дополнительных вычислений.
Максимальное изменение коэффициентов стоимости
Наряду с определением допустимых изменений запасов ресурсов представляет интерес и установление интервала допустимых изменений коэффициентов прибыли или стоимости.
Замечание. Любое изменение коэффициента целевой функции при небазисной в оптимальном решении переменной приводит лишь к тому, что в заключительной симплекс-таблице в Z-уравнении изменяется только коэффициент, соответствующий этой переменной. Причем коэффициент при небазисной переменной в результирующем Z-уравнении нужно уменьшить на ту величину, на которую он увеличивается в исходном Z-уравнении.
В заключении заметим, что разобранный пример является простейшим. Он приведен в данных указаниях, чтобы продемонстрировать выполнение отдельных этапов исследования операций и, безусловно, не охватывает все возможные ситуации, которые могут возникнуть при решении задач.
Список литературы
- Курс лекций по Теории Принятия Решений //Гапанович И.В. 2010
- Основы теории принятия решений //Орлов А.И. 2002 - 51 с