Учебная работа. Обоснование принятия оптимальных решений для хлебозавода

Обоснование принятия оптимальных решений для хлебозавода

Обоснование принятия оптимальных решений для хлебозавода

Введение

Сегодня, проблемы рационального выбора, проблемы поиска оптимальных решений, безусловно, увеличивают свою актуальность и значимость, в особенности. существенно возрастает необходимость в принятии решений различной важности, повышается ответственность за принятые решения, и более ощутимыми и значимыми становятся их последствия.

таким образом, возникает потребность в методах, инструментах и, прежде всего, позволяющих человеку находить наиболее приемлемые для него варианты решений среди множества альтернатив.

В настоящее время методы оптимальных решений широко применяются для решения различных экономических задач. Это происходит в связи с необходимостью организаций принимать наиболее выгодные для своей деятельности решении, будь то минимизация издержек производства, или максимизация прибыли.

Целью выполнения курсовой работы является овладение математическими методами решения экономических задач.

задачи выполнения типового расчета:

научиться строить экономико-математические модели;

освоить симплекс-метод табличного решения задачи линейного программирования;

освоить двойственный симплекс-метод решения задачи линейного программирования;

освоить метод потенциалов решения транспортной задачи;

освоить методику решения антагонистических игр;

освоить методику решения задачи динамического программирования.

1. Особенности обоснования принятия оптимальных решений для хлебозавода

Хлебозавод — одно из старейших предприятий хлебопекарной отрасли России. Он ведет свою историю с октября 1934 года, когда с конвейера сошла первая партия продукции.

Уже тогда, почти век назад, в эпоху индустриализации страны хлебозавод был задуман как передовое предприятие пищевой промышленности по производству хлеба и сухарных изделий. На нем впервые были установлены конвейерные печи, в 6 раз, превосходящие по мощности жаровые печи на других хлебозаводах. Вплоть до Великой отечественной войны Хлебозавод непрерывно повышал свои производственные мощности и объемы выпускаемой продукции, обеспечивая растущее население города хлебом отличного качества.

Затем завод стал круглосуточно работать на нужды фронта. В 1942 г. каждые 24 часа с его конвейеров сходила одна тонна сухарей для воинских частей. 24 августа 1942 года во время массированной бомбардировки немецкой авиации завод был частично разрушен. Но уже 1 апреля 1943 г. хлебопекарная база была восстановлена и заработала с еще большей мощностью, выпуская 6 тонн хлеба в сутки.

Хлеб всегда был важным продовольственным запасом страны, а в военные годы его бесперебойные поставки имели еще и стратегическое дальнейшее развитие и совершенствование производственной базы завода шло рука об руку с техническим прогрессом. На хлебозавод внедрялись новые технологии производства, и оставалось неизменным качество хлеба — эталонное, ГОСТовское.

Хлебозавод — предприятие, где все технологические процессы изготовления хлеба полностью автоматизированы. Производство хлеба на этих заводах ведётся на автоматических линиях, рассчитанных на выпуск определённых сортов хлебных изделий.

Мука доставляется в цистернах муковозов, а дрожжи, соль, сахар, жир и т.д. — специализированными автомашинами. На заводе мука хранится в бункерах; её транспортирование осуществляется пневматическим способом. Остальное сырьё хранится в жидком виде в ёмкостях и транспортируется по трубопроводам. Дозирование муки производится дозаторами непрерывного действия, а жидких компонентов — автоматическими дозировочными станциями. Тесто приготовляется с помощью агрегатов, включающих месильные машины непрерывного действия и аппараты для брожения. затем тесто разделяется на куски и формуется на поточных линиях, состоящих из тестоделительных, тестокруглительных, тестозакаточных машин, а также конвейерных установок для расстойки (выдержки) теста в кусках. Хлеб выпекается в автоматизированных конвейерных печах. Режим работы тестоприготовительных агрегатов, тесторазделочных линий и печей автоматически контролируется и регулируется в соответствии с качественными показателями поступающего сырья, изменением условий работы предприятия и т.д. Выпеченный хлеб подаётся на автоматизированный склад (включает хлебохранилище и экспедицию), где хранится в контейнерах до отправки в торговую сеть.

В состав хлебозавода входят 2 производства: хлебобулочное и кондитерское. ОАО » Пышка» — хлебозавод. Он находится в Туле и распространяет свою продукцию по территории, данной области.

Основная задача хлебозавода: обеспечение населения хлебобулочной и кондитерской продукцией.

Принципы и цели:

постоянное обновление и совершенствование ассортимента продукции;

использование только натурального и высококачественного сырья.

Главная задача предприятия — предлагать покупателю вкусный и полезный хлеб, сохраняющий здоровье, помогающий в работе и создающий хорошее настроение.

Хлебозавод выпускает следующие изделия:

П1 — Хлеб

Это пищевой продукт, получаемый путем выпечки, паровой обработки или жарки теста <#"justify">2. задача оптимального распределения ресурсов

Организация имеет возможность выпускать три вида изделий П1, П2, П3, При их изготовлении используется три вида ресурсов Р1, Р2, Р3. Размеры допустимых затрат ресурсов ограничены соответственно величинами b1, b2, b3. Расход ресурса i-го вида (i=1,2,…,m) на единицу изделия j-го вида (j=1,2,…,n) составляет aijден. ед. цена единицы продукции j-го вида равна сj. Требуется найти оптимальный план выпуска изделий, который обеспечивал бы организации максимальный доход.

Построим экономико-математическую модель задачи. С целью ее построения следует ввести переменные и представить исходные данные в табличном виде (Таблица 1).

Таблица 1- Переменные для решения задачи оптимального распределения ресурсов

Норма затрат РесурсыВиды изделийЗапас ресурсовСкрытые цены ресурсовyiyi*a11a12a13b1a21a22a23b2a31a32a33b3цена единицы изделияc1c2c3fmax(х)gmin(у)План выпускаxjxj*

Введем переменные: х1 — объем производства продукции первого вида; х2 — объем производства продукции второго вида; х3 — объем производства продукции третьего вида.

Представим исходные данные варианта 44 в виде таблицы 2.

Таблица 2 — исходные данные для решения задачи оптимального распределения ресурсов

Норма затрат РесурсыВиды изделийЗапас ресурсовСкрытые цены ресурсовyiyi*3323241103124цена единицы изделия341fmax(х)gmin(у)План выпускаxjxj*

целевая функция, отражающая доход от реализации произведенной продукции, представляет собой сумму произведений объема производства каждого вида продукции на

,

где n — количество видов продукции.

поскольку требуется максимизировать доход, то целевая функция стремиться к максимуму. При ресурсных ограничениях, представленных системой неравенств, левые части которых отражают затраты ресурсов каждого вида на Производство продукции соответствующего вида, а правые отражают запасы ресурсов каждого вида. Знак неравенств «меньше или равно», поскольку расход ресурсов не должен превысить имеющихся запасов:

,

где m — количество ресурсов.

Также должно выполняться условие неотрицательности переменных:

.

Таким образом, экономико-математическая модель прямой задачи линейного программирования (ПЗЛП) варианта 44 имеет вид:

при ограничениях:

Данная ПЗЛП имеет стандартную форму записи, поскольку в задаче на максимум все функциональные (ресурсные) ограничения имеют знаки «меньше или равно».

Для построения двойственной задачи линейного программирования (ДЗЛП) следует ввести двойственные переменные: у1 — скрытая цена первого ресурса; у2 — скрытая цена второго ресурса; у3 — скрытая цена третьего ресурса.

Целевая функция ДЗЛП представляет собой совокупные затраты второго предприятия на приобретение всех ресурсов первого предприятия, при этом второе предприятие стремиться, чтобы его затраты на приобретение ресурсов у первого предприятия были минимальными:

.

Ограничениями ДЗЛП является система неравенств, отражающая условия, при которых первому предприятию будет выгодно продать свои ресурсы вместо производства из них продукции, то есть при равенстве или превышении суммы, полученной от второго предприятия, над суммой дохода, полученной от реализации продукции:

.

должно выполняться условие неотрицательности переменных:

.

Также для построения ДЗЛП можно руководствоваться следующими правилами.

. В первой задаче определяется максимум линейной целевой функции, во второй — минимум.

. Коэффициенты при переменных в целевой функции первой задачи являются правыми частями в системе ограничений во второй задаче.

. Каждая из задач задана в стандартной форме, при этом в задаче максимизации все неравенства вида «меньше или равно», а в задаче минимизации все неравенства вида «больше или равно».

. Матрицы коэффициентов при переменных в системах ограничений обеих задач являются транспонированными друг к другу. Для варианта 0 матрица коэффициентов при переменных в системе ограничений ПЗЛП имеет вид:

Тогда транспонированная матрица, отражающая коэффициенты при переменных в системе ограничений ДЗЛП, имеет вид:

.

. Число неравенств в системе ограничений одной задачи совпадает с числом переменных в другой задаче.

. Условие неотрицательности переменных имеются в обеих задачах.

Таким образом, экономико-математическая модель ДЗЛП варианта 0 имеет вид:

при ограничениях

Данная ДЗЛП имеет стандартную форму записи, поскольку в задаче на минимум все функциональные (затратные) ограничения имеют знаки «больше или равно». Для решения задачи линейного программирования необходимо перейти к ее канонической форме записи, то есть перейти в системе ограничений от функциональных неравенств к равенствам посредством включения дополнительных переменных. Для ПЗЛП каноническая форма записи имеет вид:

Для перехода к канонической форме записи ПЗЛП варианта 0 следует добавить дополнительные переменные х4, х5, х6, в соответствующие неравенства со знаком «+», поскольку они отражают возможный остаток неиспользованных ресурсов:

Для ДЗЛП каноническая форма записи имеет вид:

Для перехода к канонической форме записи ДЗЛП варианта 0 следует добавить дополнительные переменные y4, y5, y6, в соответствующие неравенства со знаком «-», поскольку они отражают возможное превышение затрат на приобретение ресурсов над ценой реализации продукции (возможный убыток от производства продукции):

Во взаимодвойственных задачах линейного программирования первоначальным переменным ПЗЛП соответствуют дополнительные переменные ДЗЛП и аналогично первоначальным переменным ДЗЛП соответствуют дополнительные переменные ПЗЛП.

Установим соответствие переменных прямой и двойственной задачи для варианта 0:

Таблица

Решение задачи оптимального распределения ресурсов возможно двумя способами:

А. Одновременное решение ПЗЛП и ДЗЛП с помощью симплекс-таблиц.

Б. двойственный симплекс-метод.

Рассмотрим оба метода.

Решаем ПЗЛП.

Находим исходное опорное решение и проверяем его на оптимальность. Для этого заполняем симплекс-таблицу. Все строки таблицы 1-го шага, за исключением строки Dj (индексная строка), заполняем по данным системы ограничений и целевой функции канонической формы записи ПЗЛП.

В первую строку сjвносятся значения коэффициентов при переменных из целевой функции. В столбец базисных переменных вносим дополнительные переменные хj.

В столбец сi вносят значения коэффициентов при переменных из целевой функции, которые вошли в базис.

Рабочее поле таблицы (столбцы под переменными хjзаполняют коэффициентами при соответствующих переменных из системы равенств канонической формы записи ПЗЛП. Отсутствие той или иной переменной в равенстве означает, что ей соответствует коэффициент, равный нулю.

В столбец bi вносят свободные члены системы равенств.

Индексная строка (Dj) для переменных находится по формуле:

,

где hij — соответствующий элемент, стоящий на пересечении i-й строки и j-го столбца рабочего поля симплекс таблицы.

Индексная строка (Dj) для свободного члена находится по формуле

,

где li — соответствующий элемент, из столбца bi.

В рассматриваемом варианте в первую строку сjвносим значения коэффициентов при переменных из целевой функции 3, 4 и 1, остальные значения равны нулю. В столбец базисных переменных вносим дополнительные переменные х4, х5, х6. В рассматриваемом случае базисные переменные не входят в целевую функцию, поэтому им соответствуют нулевые значения в столбце сi.

Рабочее поле таблицы (столбцы под переменными х1, х2, х3, х4, х5, х6) заполняем коэффициентами при соответствующих переменных из системы равенств канонической формы записи ПЗЛП.

поскольку, в первом равенстве отсутствуют переменные х5 и х6, поэтому в первой рабочей строке симплекс-таблицы этим переменным соответствует коэффициент «0».

Таблица 3 — Шаг симплекс-метода

сiсj341000Базисные пер. (БП)х1х2х3х4х5х6bi0x433210030x5241010100х63120014Dj-3-4-10000

Заполняем индексную строку:

D1 = 0×3 + 0×2 + 0×3 — 3 = -3;

D2 = 0×1 + 0×4 + 0×3 — 4 = -4;

D3 = 0×2 + 0×1 + 0×2 — 1 = -1;

D4 = 0; D5 = 0; D6 = 0;

D7= 0×4 + 0×10 + 0×3 = 0.

Первое опорное решение имеет вид: = (0, 0, 0, 3, 10, 4), = 0. (компоненты опорного решения выписывают для базисных переменных из столбца свободных членов bi; переменные, не входящие в базис, имеют поэтому их компоненты в опорном решении равны нулю, базисные переменные имеют значения х4= 3, х5 = 10, х6 = 4. строки столбца. На данном шаге симплекс-метода = 0.

Следует проверить первое опорное решение на оптимальность.

возможны следующие случаи решения задачи на максимум:

если все оценки Dj³ 0, то найденное решение оптимальное;

если хотя бы одна оценка Dj£ 0, но при соответствующей переменной нет ни одного положительного коэффициента, решение задачи прекращают, так как f(X) ® ¥, т.е. целевая функция не ограничена в области допустимых решений;

если хотя бы одна оценка отрицательная, а при соответствующей переменной есть хотя бы один положительный коэффициент, то нужно перейти к другому опорному решению;

если отрицательных оценок в индексной строке несколько, то в столбец базисной переменной (БП) вводят ту переменную, которой соответствует наибольшая по абсолютной величине отрицательная оценка.

В рассматриваемом случае первое опорное решение не оптимальное, поскольку имеются три отрицательные оценки D1 = — 3, D2 = — 4, D3 = — 1.

следовательно, нужно сделать следующий шаг симплекс-метода: ввести в базис переменную, которую нужно улучшить.

Пусть одна оценка Dk£ 0 или наибольшая по абсолютной величине Dk£ 0, тогда k-й столбец принимают за ключевой (разрешающий).

Такой переменной является переменная из столбца с наибольшей по абсолютной величине оценкой |D2| = 4, т.е. переменная х2.

Таблица 4 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij0x433210030x5241010100х63120014Dj-3-4-10000

теперь следует определить, какую переменную нужно вывести из базиса.

За ключевую (разрешающую) строку принимают ту, которой соответствует минимальное отношение свободных членов (bi) к положительным коэффициентам ключевогоk-го столбца. Элемент, находящийся на пересечении ключевых строки и столбца, называют ключевым (разрешающим) элементом.

В столбце оценочных отношений min{bi/hij} отражены отношения свободных членов к элементам из разрешающего столбца, минимальным является отношение 1, соответствующее строке переменной х4, следовательно, эта строка — разрешающая, а переменная х4 выводится из базисных переменных:

Таблица 5 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij0x433210033/30×52410101010/40х631200144/1Dj-3-4-10000-

На пересечении ключевой строки и ключевого столбца находится ключевой элемент «3».

Переход ко второму шагу симплекс-метода.

В столбец базисных переменных на место выведенной из базисных вносится новая базисная переменная. При этом в столбец сi вносится соответствующее значение коэффициента при переменной из целевой функции.

Ключевая строка переписывается, при этом все ее элементы делят на ключевой элемент. Заполняют столбцы для базисных переменных: на пересечении строки и столбца, в которых находится одна и та же переменная, ставят «1», в остальных клетках столбца нули.

Остальные коэффициенты таблицы находят по правилу «прямоугольника»: для искомого элемента новой таблицы выбирают стоящий на пересечении той же строки и того же столбца элемент из предыдущей таблицы , умножающего на разрешающий, из этого произведения вычитают произведение элементов, расположенных на противоположной диагонали прямоугольника, образуемого искомым и разрешающим элементами, а полученную разность делят на разрешающий элемент

.

Оценки Dj можно считать по приведенным ранее формулам или по правилу «прямоугольника».

Заполним симплекс-таблицу второго шага для варианта 44. Вместо переменной х4 вводим в базис переменную х2, в целевой функции ей соответствует коэффициент 4, который вносится в соответствующую клетку столбца сi. В строке, соответствующей переменной х2 переписываем элементы из ключевой строки предыдущей таблицы, предварительно разделив их на ключевой элемент «3»:

Таблица 6 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij4x 2112/31/30010x 500х60Dj0-

Таблица 7 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij4x 2112/31/30010x 50100х6001Dj000-

остальные клетки заполняем по правилу прямоугольника, отмечая, что ключевой элемент «3» находится в третьей строке и втором столбце предыдущей таблицы:

= ;

= ;

= ;

= ;

= ;

= ;

= ;

= .

Таблица 8 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij4x 2112/31/30010x 5-20-5/3-4/31060х6204/3-1/3013Dj000-

Находим оценки переменных:

D1 = 0×(2) + 0×(-2) + 4×(1) — 3 = 1;

D3 = 0×(4/3) + 0×(-5/3) + 4×(2/3) — 1 = 5/3;

D4 = 0×(-1/3) + 0×(-4/3) + 4×(1/3) — 0 =4/3;

D7= 0×(3) + 0×(6) + 4×(1) = 4.

Таблица 9 — Шаг симплекс-метода

сiсj341000minБазисные пер. (БП)х1х2х3х4х5х6bibi/hij4x 2112/31/30010x 5-20-5/3-4/31060х6204/3-1/3013Dj105/34/3004-

поскольку все оценки положительные Dj, то полученное опорное решение является оптимальным: = (0, 1, 0, 0, 6, 3), а максимальное значение целевой функции равно = 4.

учитывая соответствие переменных взаимодвойственных задач, можно выписать оптимальное решение двойственной задачи. Значения соответствующих переменных берут из последней строки последней симплекс-таблицы:

Таблица 10 — Шаг симплекс-метода

сiсj341000Базисные пер. (БП)х1х2х3х4х5х6bi4x 2112/31/30010x 5-20-5/3-4/31060x 6204/3-1/3013Dj105/34/3004yiy4y5y6y1y2y3

таким образом, = (4/3, 0, 0, 1, 0, 5/3), поскольку y1«x4, y2«x5, y3«x6, y4«x1, y5«x2, y6«x3. Экономический смысл переменных, участвующих в решении отражен в таблице 11.

Таблица 11 — Экономический смысл переменных

компоненты оптимального решения ПЗЛППлан производстваОстатки ресурсов, единиц

Превышение затрат на ресурсы над ценой реал. (возможный убыток от производства продукции)Объективно обусловленные оценки ресурсов (теневые, условные, скрытые цены ресурсов)Компоненты оптимального решения ДЗЛП

Экономический смысл переменных:

x1*, x2*, x3*-основные переменные — оптимальный план производства;

x4*,x5, x6*- дополнительные переменные — остатки ресурсов;

y1,y2,y3-основные переменные — скрытые цены;

y4,y5,y6-дополнительные переменны — превышение затрат на ресурсы над ценой реализации (возможный убыток от производства продукции).

Анализ решения ПЗЛП

Подставим оптимальные значения переменных x*в исходную систему ограничений ПЗЛП:

) 3×х1 + 3×х2 + 2×х3£ 3

×0+ 3×1 + 2×0 = 3

=3, следовательно, х4 = 0, ресурс Р1 используется полностью;

) 2×х1 + 4×х2 + х3£ 10

×0+ 4×1 + 1×0 = 4

< 10, следовательно, х5 = 6, ресурс Р2 используется не полностью;

) 3×х1 + х2 + 2×х3£ 4

×0 + 1×1 + 2×0 = 1

< 4, следовательно, х4 = 3, ресурс Р3 используется не полностью.

анализ решения ДЗЛП

Подставим оптимальные значения переменных у*в исходную систему ограничений ПЗЛП:

1) 3×у1 + 2×у2 + 3×у3³ 3

3×(4/3) + 2×0 + 3×0 = 4

³ 3, на 1, следовательно, у4 = 1, убытки от производства первого вида продукции П1, которая не вошла в оптимальный план производства, в случае ее производства будут составлять 1 ден. ед. с каждого изделия первого вида;

) 3×у1 + 4×у2 + у3³ 4

×(4/3) + 4×0 + 1×0 = 4

= 4, следовательно, у5 = 0, убытки от производства второго вида продукции П2, которая вошла в оптимальный план производства отсутствуют;

) 2×у1 + у2 + 2×у3³ 1

×(4/3) + 1×0 + 2×0 = 8/3

/3 ³ 1 на 5/3, следовательно, у6 = 5/3, убытки от производства третьего вида продукции П3, которая не вошла в оптимальный план производства, в случае ее производства будут составлять 5/3 ден. ед. с каждого изделия третьего вида.

Рассчитаем границы изменения дефицитных ресурсов, в пределах которых не изменится структура оптимального плана. Ресурсы, которые используются полностью, называются дефицитными. Признаком дефицитности ресурсов является отличие от нуля соответствующей данному ресурсу двойственной переменной и равенство нулю соответствующей дополнительной переменной.

В рассматриваемом случае дефицитным является ресурс Р1.

Для исследования границ изменения первого вида ресурса Р1 из последней симплекс-таблицы составляют систему неравенств для базисных переменных ПЗЛП, используя элементы из столбца свободных членов bi и столбца, соответствующего переменной у1. Коэффициенты из столбца «у1» умножают на искомое изменение Db1запаса ресурса Р1:

Þ.

учитывая, что b1 = 3, допустимый интервал изменения границ первого вида ресурса составит или .

Проверим границы с помощью вычислительного эксперимента:

Зададим b1 = 7,48:

рисунок 1 — Прямая задача линейного программирования

структура ПЗЛП не изменилась

Рисунок 2 — двойственная задача линейного программирования

Структура ДЗЛП не изменилась.

Зададим b1 = 0,02:

Рисунок 3 — Прямая задача линейного программирования

структура решения ПЗЛП не изменилась

Рисунок 4 — двойственная задача линейного программирования

Структура ДЗЛП не изменилась.

Зададим b1 = 7,51:

Рисунок 5 — Прямая задача линейного программирования

Получили новую структуру оптимального плана прямой задачи: (х1=0,01; х2=2,5; x3=0)

рисунок 6 — Двойственная задача линейного программирования

Получили новую структуру оптимального плана ДЗЛП.

Можно сделать вывод о правильности установленных границ изменения ресурса Р1: 0≤ b2 ≤ 15/2.

соответствующей дополнительной переменной.

В рассматриваемом случае не дефицитными являются ресурсы Р2 и Р3. оптимальный план производства не изменится при .

Проверим границы с помощью вычислительного эксперимента:

Зададим b2 = 3,98:

Рисунок 7 — Двойственная задача линейного программирования

Получили новый оптимальный план решения ДЗЛП.

Зададим b2 = 4,02:

рисунок 8 — Прямая задача линейного программирования

оптимальный план решения ПЗЛП не изменился

рисунок 9 — Двойственная задача линейного программирования

оптимальный план решения ДЗЛП не изменился.

Можно сделать вывод о правильности установления границы изменения ресурса Р2: b2 ≥ 4.

оптимальный план производства для не дефицитного ресурса Р3 не изменится при .

Проверим границы с помощью вычислительного эксперимента:

Зададим b3 = 0,98:

Рисунок 10 — Двойственная задача линейного программирования

Получили новый оптимальный план решения ДЗЛП.

Зададим b2 = 1,02:

рисунок 11 — Прямая задача линейного программирования

оптимальный план решения ПЗЛП не изменился.

рисунок 12 — Двойственная задача линейного программирования

оптимальный план решения ДЗЛП не изменился.

Можно сделать вывод о правильности установления границы изменения ресурса Р3: b3 ≥ 1.

Рассчитаем границы изменения цены изделия, попавших в оптимальный план производства, в пределах которых оптимальный план не изменится.

В план производства вошел второй вид продукции П2.

Для второго вида продукции П2, которая попала в план производства, из последней симплекс-таблицы составляют систему неравенств для базисных переменных ДЗЛП (они соответствуют свободным переменным ПЗЛП), используя элементы из строки Dj и строки, соответствующей переменной х2. Коэффициенты из строки «х2» умножают на искомое изменение Dс1цены продукции П2

Þ.

учитывая цену второго вида продукции с2 = 4, интервал устойчивости изменения цен составит или .

Проверим границы с помощью вычислительного эксперимента:

Зададим c2 = 2,98:

Рисунок 13 — Прямая задача линейного программирования

Получили новый оптимальный план решения ПЗЛП.

Зададим c2 = 3,02:

рисунок 14 — Прямая задача линейного программирования

структура оптимального плана решения ПЗЛП не изменилась.

Рисунок 15 — Двойственная задача линейного программирования

структура решения ДЗЛП не изменилась.

Можно сделать вывод о правильности определения границы изменения цены с2: с2 3.

Для видов продукции, не попавших в оптимальный план производства, исследование допустимых границ изменения цен не проводится.

Определим величины ∆bs ресурса Рs, введением которого в Производство можно компенсировать убыток и сохранить максимальный Доход на прежнем уровне (ресурсы предполагаются взаимно заменяемыми), получаемый при исключении из производства ∆br единиц ресурса Рr.

В рассматриваемом случае: r = 1; ∆br = 0,2; s = 2.

Для взаимозаменяемых ресурсов (коэффициент взаимозаменяемости >0, но отличен от бесконечности) количество ресурса ∆biвида i, необходимое для замены выбывающего количества ∆bkресурса k, определяется по формуле:

.

таким образом, . Следовательно, замена первого ресурса невозможна.

Оценим целесообразность приобретения ∆bk единиц ресурса Рk по цене сk за единицу. Для оценки целесообразности приобретения дополнительного количества ресурса ∆biвида i по цене сk необходимо сравнить предлагаемую цену с рассчитанной ранее теневой ценой этого ресурса . Приобретение дополнительного количества ресурса целесообразно, если выполняется условие непревышения новой цены над теневой ценой

сk£.

В противном случае приобретение дополнительного количества ресурса нецелесообразно.

В рассматриваемом случае: ∆bk = 0,2; k = 3, ck = 15.

Поскольку < 15, то приобретение дополнительного количества ресурса не целесообразно.

Оценим целесообразность выпуска нового изделия П4, на единицу которого ресурсы Р1, Р2, Р3 расходуются в количествах a14, a24, a34 единиц, а цена единицы изделия составляет с4 денежных единиц. Включение дополнительного вида продукции n+1 в план производства целесообразно, если соотношение дополнительных затрат и цены реализации дополнительного вида продукции удовлетворяет следующему условию

.

Расчет затрат осуществим по формуле

ден. ед.

учитывая, что затраты на ресурсы для производства продукции третьего вида меньше цены реализации с4 = 35 ден. ед., то включение ее в план производства целесообразно.

Решим прямую и двойственную задачу линейного программирования в среде MicrosoftExсel. Решение задач линейного программирования с помощью программы MicrosoftExcel осуществляется через меню Сервис и вкладку Поиск решения. Если данная вкладка не установлена, то ее установка осуществляется следующими действиями:

. Войти в меню Сервис.

. выбрать команду Надстройки.

. В появившемся диалоговом окне Надстройки установить флажок напротив строки Поиск решения и нажать кнопку ОК.

После произведенных действий в меню Сервис появится команда Поиск решения.

Прежде чем начать выполнение команды Поиск решения следует провести подготовительные действия по введению данных задания в таблицу.

выбираем ячейку для введения целевой функции (например, ячейку А1). При записи целевой функции в ячейку А1 вместо значений переменной хi подставляют названия пустых ячеек, в которых хотим получить искомые значения х1, х2, х3, например, С1, С2, С3. Тогда запись целевой функции в ячейке А1 будет иметь вид = 3*С1+4*С2+С3. после введения целевой функции в ячейку А1 нажимаем клавишу «Enter», в ячейке А1 отобразится 0.

Условия ограничений вписываем в ячейки столбца В. В выбранные ячейки записываются только левые части неравенств в следующем виде:

Ячейка В1: =3*С1+3*С2+2*С3 «Enter»

Ячейка В2: =2*С1+4*С2+С3 «Enter»

Ячейка В3: =3*С1+С2+2*С3 «Enter»

Ячейка В4: =С1 «Enter»

Ячейка В5: =С2 «Enter»

Ячейка В6: =С3 «Enter»

Через меню Сервис/Поиск решения открыть окно поиска решения:

рисунок 16 — Решение задачи в среде MicrosoftExсel

В поле ввода Установить целевую ячейку вводим ссылку на ячейку А1.

В поле ввода Изменяя ячейки укажем ссылки на ячейки С1:С3.

Данная операция осуществляется следующими действиями:

.Щелкнуть левой кнопкой мыши в поле ввода Изменяя ячейки.

. Выделить при помощи левой кнопки мыши ячейки, начиная с С1 и до ячейки С3(поле ввода изменения ячейки должно автоматически заполниться).

В поле ввода Ограничения введем ограничения, соответствующие ячейкам В1, В2, В3. Ввод значений осуществляется в следующем порядке:

.нажать кнопку Добавить. Появится диалоговое окно Добавление ограничения.

. Для каждого логического выражения, находящихся в ячейках В1, В2, В3 вводим свое условие и свое ограничение, последовательно нажимая кнопку Добавить. По окончании нажать кнопку ОК.

Рисунок 17 — Решение задачи в среде MicrosoftExсel

. Должен появиться следующий результат:

Рисунок 18 — Результат решения задачи в среде MicrosoftExсel

Для поиска максимального значения устанавливаем флажок.

Для осуществления вычислений нажмем кнопку Выполнить. Откроется окно диалога Результаты поиска решения. В окне Тип отчета выберем строку Результаты и нажмем кнопку ОК. Получим следующие результаты:

В ячейке А1отражено максимально возможное В столбце В отражены значения логических выражений, удовлетворяющих исходным ограничениям.

х1 + 3х2 +2х3 = 3

х1 + 4х2 + х3 = 4

х1 + х2 + 2х3 = 1

В столбце С отражены значения переменных, при которых значение функции принимает максимальную величину:

х1 = 0; х2 = 1; х3 = 0.

На листе «отчет по результатам» появится следующая информация, отражающая результаты расчета:

Рисунок 19 — результат решения задачи в среде MicrosoftExсel

все действия для двойственной задачи, получим следующий отчет по результатам:

Рисунок 20 — Результат решения задачи в среде MicrosoftExсel

симплекс-метод

. Выражение базисных переменных ПЗЛП и ДЗЛП через свободные.

.Определение исходного решения прямой и двойственной задач и проверка его на оптимальность.

Для этого заполняют симплекс-таблицу двойственного симплекс-метода. Таблица имеет два входа: по переменным прямой задачи xjи по переменным двойственной задачи yi. базисные переменные ПЗЛП записывают в столбце xбаз, а свободные переменные в строке xсв со знаком «-». базисные переменные ДЗЛП записывают в строке убаз, а свободные переменные в столбце усв.

Строки таблицы 1-го шага (верхние части клеток) заполняют по данным системы ограничений и целевой функции прямой задачи. Столбцы таблицы 1-го шага (верхние части клеток) заполняют по данным системы ограничений и целевой функции двойственной задачи. Значение целевой функции равно нулю, поскольку базисные переменные не входят в выражение целевой функции.

Решение является оптимальным, если в строке и столбце bi нет отрицательных элементов. В противном случае следует провести изменения в базисных переменных.

Симплексная таблица двойственного симплекс-метода для варианта 40 имеет следующий вид:

Таблица 12 — Шаг двойственного симплекс-метода

Решение ПЗЛП выписывается по строкам, значения базисных переменных берутся из столбца bi, если переменная с соответствующим индексом не входит в базис, то ее Решение ДЗЛП выписывается по столбцам, значения базисных переменных берутся из строки cj, если переменная с соответствующим индексом не входит в базис, то ее Данное решение не является оптимальным, поскольку решение не допустимое (не выполнено условие неотрицательности переменных), ему соответствуют отрицательные элементы в строке . поэтому следует провести замену переменных в базисе.

Для выбора разрешающей строки находят отрицательный элемент в строке .В столбце над этим найденным элементом выбирают любой положительный элемент, эта строка — разрешающая.

Если в столбце над найденным элементом нет положительных элементов, то ПЗЛП не имеет смысла (целевая функция не ограничена в области допустимых решений), а ДЗЛП не имеет решения.

В рассматриваемом примере в строке три отрицательных элемента.

Таблица 13 — Шаг двойственного симплекс-метода

Можно выбрать любой из них, поэтому выбираем «-1», а над ним выбираем элемент из второй строки «2». следовательно, первая строка — разрешающая (выделена жирным шрифтом).

Выберем разрешающий столбец. Элементы строкиделят на соответствующие элементы разрешающей строки под переменными. Из полученных отношений выбирают максимальное отрицательное отношение, этот столбец — разрешающий (максимальным из отрицательных отношений может быть отношение — отрицательный ноль).

Если среди полученных отношений нет отрицательных, то ПЗЛП не имеет решения, ДЗЛП не имеет смысла или решения.

На пересечении разрешающей строки и разрешающего столбца находится разрешающий элемент.

В рассматриваемом примере максимальным среди отношений элементов из строкии элементов разрешающей строки является отношение в третьем столбце . В качестве разрешающего выберем третий столбец (выделен жирным шрифтом).

Таблица 14 — Шаг двойственного симплекс-метода

Разрешающим элементом является элемент «2», который стоит на пересечении первой строки и третьего столбца.

Под разрешающим элементом всегда ставят «1». Остальные элементы разрешающей строки переписываются без изменений, а остальные элементы разрешающего столбца переписываются с противоположным знаком.

Заполнение нижних частей клеток для разрешающей строки и разрешающего столбца варианта 44:

Таблица 15 — Шаг двойственного симплекс-метода

Остальные элементы таблицы в нижних частях клеток находят по правилу прямоугольника: для искомого элемента из нижней части клетки элемент из верхней части этой же клетки умножают на разрешающий элемент и из этого произведения вычитают произведение элементов, расположенных на противоположной диагонали прямоугольника, образуемого искомым и разрешающим элементами (все элементы из верхних клеток).

заполнение нижних частей для остальных клеток варианта 44:

на пересечении второй строки и первого столбца:

Рис. 2×2 — 3×1 = 1;

клетка на пересечении третьей строки и первого столбца:

Рис. 2×3 — 2×3 = 0;

клетка на пересечении четвертой строки и первого столбца:

Рис. — 3×2 — 3×(-1) = -3;

клетка на пересечении второй строки и второго столбца:

Рис. 4×2 — 3×1 = 5;

клетка на пересечении третьей строки и второго столбца:

Рис. 2×1 — 3×2 = -4;

клетка на пересечении четвертой строки и второго столбца:

-4×2 — (-1)×3 = -5;

клетка на пересечении второй строки и четвертого столбца:

Рис.×10 — 1×3 = 17;

клетка на пересечении третьей строки и четвертого столбца:

Рис.×4 — 2×3 = 2;

клетка на пересечении четвертой строки и четвертого столбца:

Рис. 0×2 — (-1)×3= 3.

Таким образом, получаем следующую таблицу:

Таблица 16 — Шаг двойственного симплекс-метода

Изменяют базисные переменные: меняют местами переменные из разрешающей строки и разрешающего столбца. Элементы из нижних клеток предыдущей симплекс-таблицы делят на верхний разрешающий элемент и записывают на соответствующие места в верхние клетки новой симплекс-таблицы.

Если в новой таблице в строке есть отрицательные элементы, то следует сделать следующий шаг симплекс-метода. (Нецелесообразно выбирать за разрешающую строку — те же строки, что и на предыдущих шагах).

Если в новой таблице в строке нет отрицательных элементов, а в столбце свободных членов остались отрицательные элементы, то строка с отрицательным значением выбирается за разрешающую, и выполняется следующий шаг симплекс-метода.

Если в строке есть нулевой элемент, то это признак альтернативного оптимума для ПЗЛП. Для нахождения альтернативного решения выполняется еще один шаг симплекс-метода: Столбец с нулевым элементом в строке выбирается за разрешающий. Находят неотрицательные отношения столбца свободных членов к соответствующим элементам разрешающего столбца. Из полученных отношений выбирают минимальное неотрицательное отношение — это разрешающая строка, разрешающий элемент найден. Затем выполняют еще один шаг симплекс-метода.

Если в столбце есть нулевой элемент, то это признак альтернативного оптимума для ДЗЛП. Для нахождения альтернативного решения выполняется еще один шаг симплекс-метода. При этом строка с нулевым элементом в столбце выбирается за разрешающую.

переходим к следующей симплекс-таблице. При этом в первую строку включается пара переменных х3, у6, соответствующих разрешающему столбцу, а в третий столбец вводится пара переменных х4, у1, соответствующая разрешающей строке. Верхние части клеток заполняются элементами, полученными в результате деления соответствующих (стоящих на аналогичном месте) элементов из предыдущей таблицы в нижних частях клеток на разрешающий элемент «2»:

Таблица 17 — Шаг двойственного симплекс-метода

Решение ПЗЛП на втором шаге двойственного симплекс-метода также выписывается по строкам: = (0, 0, 3/2, 0, 17/2, 1), = 3/2, решение ДЗЛП выписывается по столбцам: = (1/2, 0, 0, -3/2, -5/2, 0), = 3/2. Данное решение также не оптимальное, поскольку в строке еще остались отрицательные элементы. необходимо продолжить поиск оптимального решения.

В строке выбираем отрицательный элемент «-3/2». Над ним выбираем положительный, предпочтительнее «3/2», следовательно, первая строка — разрешающая (выделена жирным шрифтом):

Таблица 18 — Шаг двойственного симплекс-метода

Находим максимальное отношение среди отношений элементов в строке целевой функции к соответствующим элементам разрешающей строки

. Разрешающий столбец — первый, поскольку максимальное отношение соответствует ему.

Разрешающий элемент «3/2» находится на пересечении разрешающей строки и разрешающего столбца.

Заполнение нижних частей клеток осуществляется аналогично рассмотренному выше. Получаем следующую таблицу:

Таблица 19 — Шаг двойственного симплекс-метода

Переходим к следующей симплекс-таблице. При этом в первую строку включается пара переменных х1, у4, соответствующих разрешающему столбцу, а в первый столбец вводится пара переменных х3, у6, соответствующая разрешающей строке. Верхние части клеток заполняются элементами, полученными в результате деления соответствующих (стоящих на аналогичном месте) элементов из предыдущей таблицы в нижних частях клеток на разрешающий элемент «3/2»:

Таблица 20 — Шаг двойственного симплекс-метода

Решение ПЗЛП на третьем шаге двойственного симплекс-метода также выписывается по строкам: = (1, 0, 0, 0, 8, 1), = 3, решение ДЗЛП выписывается по столбцам: = (1, 0, 0, 0, -1, 1), = 3. Данное решение также не оптимальное, поскольку в строке еще остался отрицательный элемент. необходимо продолжить поиск оптимального решения.

В строке выбираем отрицательный элемент «-1». Над ним выбираем положительный, предпочтительнее «1», следовательно, первая строка — разрешающая (выделена жирным шрифтом):

Таблица 21 — Шаг двойственного симплекс-метода

Находим максимальное отношение среди отношений элементов в строке целевой функции к соответствующим элементам разрешающей строки. Разрешающий столбец — второй, поскольку максимальное отношение соответствует ему.

Разрешающий элемент «1» находится на пересечении разрешающей строки и разрешающего столбца.

Заполнение нижних частей клеток осуществляется аналогично рассмотренному выше. Получаем следующую таблицу:

Таблица 22 — Шаг двойственного симплекс-метода

Перейдем к следующей симплекс-таблице. При этом в первую строку включается пара переменных х2, у5, соответствующих разрешающему столбцу, а во второй столбец вводится пара переменных х1, у4, соответствующая разрешающей строке. клетки следующей таблицы заполняются элементами, полученными в результате деления соответствующих элементов из предыдущей таблицы в нижних частях клеток на разрешающий элемент «1». поскольку в нижних частях клеток таблицы все элементы положительные и разрешающий элемент также положительный, то в следующей таблице будет получено оптимальное решение и нет необходимости делить на две части клетки в последней таблице.

Таблица 23 — Шаг двойственного симплекс-метода

Оптимальное решение ПЗЛП: = (0, 1, 0, 0, 6, 3), = 4, оптимальное решение ДЗЛП: = (4/3, 0, 0, 1, 0,5/3), = 4.

Данное решение аналогично решению, полученному в соответствии с методом одновременного решения пары взаимодвойственных задач для ПЗЛП и совпадает с решением ДЗЛП в отношении значения целевой функции.

3. Метод потенциалов решения транспортной задачи

Начальное решение транспортной задачи

Имеется 4 оптовых склада и 4 магазина.

A1, A2, A3, A4 — склады.1 = 19, a2 = 19, a3 = 19, a4 = 19 — соответственно, запасы на складах.

B1, B2, B3, B4 — магазины.

b1 = 15, b2 = 15, b3 = 23, b4 = 23 — соответственно, потребности магазинов.

Найдем исходное опорное решение по методу минимальной стоимости. Для этого заполним распределительную таблицу:

Таблица 24 — Транспортная задача

Таблица 25 — Транспортная задача

Таким образом, все поставки распределены, получено начальное решение транспортной задачи:

Рис.

fo(xo) = 12·19 + 16·15 + 15·4 + 16·19 + 10·15 + 10·4 = 1022

Количество ненулевых элементов 6 < 7 , решение вырожденное. Введем значимый ноль (Õ).

Решение транспортной задачи методом потенциалов

Таблица 26 — Транспортная задача

Вычисляем оценки свободных клеток:

D11 = c11 — u1 -v1 = 21- 0 — 13= 8 > 0,

D12 = c12- u1 — v2 = 17 — 0 -11 = 6 >0,

D14 = c14 — u1 — v4 = 24 — 0 — 11= 13 >0,

D22 = c22 — u2 — v2 = 19 — 3 — 11 = 5 > 0,

D24 = c24 — u2 — v4 = 19 — 3 -11 = 5 > 0,

D31 = c31 — u3 — v1 = 17 — 5 — 13= — 1< 0,

D32 = c32 — u3 — v2 = 15 — 5 — 11= -1< 0,

D33 = c33 — u3 — v3 = 24 — 5- 12 = 7 >0,

D41 = c41 — u4 — v1 = 13+1-13 = 1 >0.

В распределительной таблице оценки клеток проставлены в нижнем левом углу в скобках:

Таблица 27 — Транспортная задача

Переход к следующему опорному решению.

Выбираем клетку, от которой начнем построение цикла перераспределения поставок, по правилу . У вершин цикла с соответствующими значениями поставок по правилу чередования знаков ставим знаки (+) и (-).У вершин со знаком (-) выбираем минимальный груз q = min[15,19] = 15.

Рис.

Проверка решения Х1 на оптимальность выполняется аналогично предыдущему шагу.

Таблица 28 — Транспортная задача

У вершин со знаком (-) выбираем минимальный груз q = min[15,4,0] = 0.

Рис.

Проверка решения Х2 на оптимальность.

Результаты расчета представлены в распределительной таблице:

Таблица 29 — Транспортная задача

Все оценки свободных клеток положительные, следовательно, решение оптимально. Значение целевой функции:

F*(x2) = 12·19 + 16·15 + 15·4 + 15·15 + 16·4 + 10·19 = 1007

таким образом, решение транспортной задачи:

Рис.

Решение транспортной задачи в среде Microsoft Exсel

Ввод исходных данных (в области C3:F6 — тарифы на перевозку продукции; в столбце G3:G6 — запасы; в ячейках С7, D7, E7, F7 — потребности).

В области решения в ячейке G10 введите формулу стоимости перевозок:

=СУММПРОИЗВ(C12:F15;C3:F6). Для этого необходимо нажать на значок f(x) на панели инструментов, выбрать математическую функцию СУММПРОИЗВ и ввести два массива C12:F15 и C3:F6. далее в области C12:F15 проставьте любое первоначальное решение (например, единицы)/

В ячейке С16 записывается формула: =СУММ(C12:C14), т.е. сумма значений по столбцу (можно выделить значения столбца и нажать на знак автосуммы Σ на панели инструментов). Аналогично в D16, E16, F16. автоматически суммируются значения по столбцам.

В ячейке G12 записывается формула: =СУММ(C12:F12), т.е. сумма значений по строке (можно выделить значения строки и нажать на знак автосуммы Σ на панели инструментов). Аналогично в G13, G14, G15. автоматически суммируются значения по строкам:

Рисунок 21 — Поиск оптимального решения

Далее выполняют команду Поиск решения (вкладка Сервис или Данные).

Установить целевую ячейку G10, равной минимальному значению.

В поле ввода Изменяя ячейки установить C12:F15

В поле ввода Ограничения установить C12:F15 >= 0

C16:E16 = C7:E7

G12:G15 = G3:G6

рисунок 22 — Поиск оптимального решения

Вычисления производятся при нажатии кнопки Найти решение. Получим решение задачи:

Рисунок 23 — Поиск оптимального решения

4. Решение игры

рассмотрим игру, представленную платежной матрицей:

a ¹ b.

Так как игра не имеет седловой точки, то она не решается в чистых стратегиях.

Все элементы четвертого столбца не превышают элементов первого столбца. Следовательно, стратегия, соответствующая первому столбцу платежной матрицы заведомо невыгодна для второго игрока, и первый столбец можно исключить. Все элементы второго столбца не превышают элементов пятого столбца. следовательно, стратегия, соответствующая пятому столбцу платежной матрицы заведомо невыгодна для второго игрока, и пятый столбец можно исключить. Получаем матрицу из трех столбцов:

Решим упрощенную игру с помощью графического доминирования.

Рисунок 24 — Геометрическая интерпретация исходной игры

хлебопекарный игра отрасль

С помощью геометрического доминирования исходная игра сведена к игре размерности [2х2]

Решим полученную игру аналитическим методом:

Проверка:

V(H)=x*HY*T=

ЗЛП для игрока 1:

min [ f0 (х)=х1+х2]

ЗЛП для игрока 2:

max [ g0(y)=y1+y2+y3+y4+y5]

Найдем решение задачи линейного программирования для первого игрока с помощью программы Microsoft Excel.

Введем целевую функцию в ячейку А1, ограничения в ячейки В1 — В5, а искомые значения x в ячейки С1 — С2. Далее с помощью поиска решений найдем оптимальное решения. Получим следующее:

Рисунок 25 — Поиск оптимального решения

Найдем решение задачи линейного программирования для второго игрока с помощью программы Microsoft Excel.

Введем целевую функцию в ячейку А1, ограничения в ячейки В1 — В2, а искомые значения y в ячейки С1 — С5. Далее с помощью поиска решений найдем оптимальное решения. Получим следующее:

Рисунок 26 — Поиск оптимального решения

X=(1/7; 1/7); Y=(0; 1/7; 0; 1/7; 0); ,

g0(y*)=f0(x*)=2/7.

Тогда

V(H)=1/(2/7)=7/2;

X*1=(1/7)/(2/7)=1/2;

X*2=(1/7)/(2/7)=1/2;

Y*1=0/(2/7)=0;

Y*2=(1/7)/(2/7)=1/2;

Y*3=0/(2/7)=0;

Y*4=(1/7)/(2/7)=1/2;

Y*5=0/(2/7)=0;

Следовательно:

Х= (1/2; 1/2),

Y= (0; 1/2; 0; 1/2; 0).

V(H) = 7/2

5. Задача динамического программирования

На развитие трех предприятий выделено В млн. руб. Известна эффективность капитальных вложений xi в каждое j-е предприятие, заданная таблично значением нелинейной функции fj(xi), где ,, n — количество предприятий, m — количество возможных сумм капитальных вложений.

необходимо распределить выделенные средства между предприятиями таким образом, чтобы получить максимальный суммарный Доход.

Таблица 30- Данные для решения задачи динамического программирования

Объем капиталовложений Xi (тыс. руб.)Прирост выпуска продукции ѓі (Xi) в зависимости от объема капиталовложений (тыс. руб.)предприятие 1предприятие 2предприятие 30 100 200 300 400 500 600 7000 30 50 90 100 160 170 2100 50 80 90 150 190 210 2200 40 50 110 120 180 220 240

Математическая модель задачи.

Определить х* = (,, …, , …, ), обеспечивающий максимум целевой функции

и удовлетворяющий условиям

,

Математическая модель задачи :

при ограничениях:

,

.

Условная оптимизация.

Максимально возможный доход, который может быть получен с предприятий (с k-го по n-е), определяется с помощью функции Беллмана:

,

где Сk — количество средств, инвестируемых вk-е предприятие, 0≤ Сk≤ В.

На первом шаге условной оптимизации при k = n функция Беллмана представляет собой Прибыль только с n-го предприятия. При этом на его инвестирование может остаться количество средств Сn, 0≤ Сn≤ В. Чтобы получить максимум прибыли с этого предприятия, можно вложить в него все эти средства, т.е.

Fn(Сn) = fn(Сn) и хn = Сn.

Для упрощения расчетов предполагаем, что распределение средств осуществляется в целых числах xi = {0,100,200,300,400,500, 600, 700} тыс. руб.

На первом шаге все капиталовложения идут на 1-е предприятие (Таблица 31).

Таблица 31 — первый шаг задачи динамического программирования

В шапке таблицы отражены варианты значений капиталовложений х1, которые могут быть предоставлены первому предприятию. В столбце C1 отражены варианты значений капиталовложений, которые могут быть выделены всем трем предприятиям в совокупности.

Предположим, что все средства в количестве x1 = 700тыс. руб. отданы первому предприятию. В этом случае максимальный доход составит f1(x1) = 700 тыс. руб., следовательно: F1(C1) = f1(x1) и x1= C1.

На втором шаге определяем оптимальную стратегию при распределении денежных средств между вторым и первым предприятиями. При этом рекуррентное соотношение Беллмана имеет вид:

.

Таблица 32 — второй шаг задачи динамического программирования

В шапке таблицы отражены варианты значений капиталовложений х2, которые могут быть предоставлены второму предприятию при условии, что часть средств выделяется первому предприятию. В клетках таблицы первое слагаемое — это возможный прирост выпуска продукции второго предприятия f2(х2) в результате освоения капиталовложений х2; второе слагаемое — возможный прирост выпуска продукции первого предприятия, если ему будет выделена оставшаяся часть капиталовложений, определяемая как C2 — х2.

На третьем шаге определяем оптимальную стратегию при распределении денежных средств между третьим и двумя другими предприятиями, используя следующую формулу для расчета суммарного дохода:

,

В шапке таблицы отражены варианты значений капиталовложений х1, которые могут быть предоставлены третьему предприятию при условии, что часть средств выделяется второму и первому предприятию. В клетках таблицы первое слагаемое — это возможный прирост выпуска продукции первого предприятия f1(х1) в результате освоения капиталовложений х1; второе слагаемое — возможный прирост выпуска продукции второго и первого предприятий, если им будет выделена оставшаяся часть капиталовложений, определяемая как C1 — х1.

Таблица 33 — третий шаг задачи динамического программирования

Значение функции Беллмана F1(С1) представляет собой максимально возможный доход со всех предприятий, а оптимальным количеством средств, вложенных в третье предприятие.

следовательно, значение целевой функции равно Fmax(x*) = 270 тыс. руб. безусловная оптимизация.

Далее на этапе безусловной оптимизации для всех последующих шагов вычисляется величина Сk = (Сk-1 — хk-1) оптимальным управлением на k-м шаге является то соответствующем состоянии системы Sk. Определяем компоненты оптимальной стратегии. Для этого значения функций Беллмана и соответствующие им оптимальные значения х вносим в итоговую таблицу 34.

Таблица 34 — обратный ход

По данным из таблицы 34 максимальный доход при распределении 700 тыс. руб. между тремя предприятиями составляет: C1 = 700, F3(700) = 270 тыс. руб.

таким образом, возможный вариант оптимального плана инвестирования предприятий:

х* = (0, 100, 600), который обеспечит максимальный доход, равный

F(700) = f1 (0) + f2(100) + f3(600) = 0 + 50 + 220 = 270 тыс. руб.;

Заключение

В первом разделе курсовой работы требовалось найти оптимальный план выпуска изделий, который обеспечивал бы организации максимальный доход. Оптимальный план выпуска = 4 единицы. Были найдены решения В ДЗЛП И ПЗЛП: X* = (0, 1, 0, 0, 6, 3) и Y* = (4/3;0;0;1;0;5/3). Следовало произвести расчет границ изменения дефицитных ресурсов, в пределах которых не изменится структура оптимального плана. В рассмотренном случае лишь один ресурс используется полностью, следовательно, является дефицитным. Это ресурс Р1: 0 £ b1 +Db1 £ 15/ 2. Также был произведен расчет границ изменения цены изделия, попавших в оптимальный план производства, в пределах которых оптимальный план не изменится.

В план производства вошел второй вид продукции П2 и цена на него устанавливается в размере 4 ден. ед. В курсовой работе произведена оценка целесообразности приобретения дополнительного количества ресурса и была доказано, что данная процедура не целесообразна. Произведена оценка целесообразности выпуска нового изделия. Учитывая, что затраты на ресурсы для производства продукции третьего вида меньше цены реализации с4 = 15 ден. ед., то включение ее в план производства целесообразно. В курсовой работе был найден оптимальный план для хлебозавода, оптимальный план перевозки, который равен 1007 денежным единицам, при котором общая стоимость перевозок минимальна, была обоснована ценовая стратегия, цена игры равна 7/2 денежных единиц, было произведено распределение ресурсов между проектами для того, чтобы получить максимальный суммарный доход, который равен 270 тыс. руб.. Ко всем расчетам мною были приложены решения в среде Microsoft Excel. таким образом, все поставленные задачи курсовой работы были выполнены, все методы и методики решения были проанализированы на конкретном виде предприятия.

Список литературы

. «Математические методы в программировании: Учебник» / Агальцов В.П. — М.:ИД «форум», 2013

. «Динамическое программирование» / Окулов С.М., Пестов О.А. — М.: бином. Лаборатория знаний, 2012

. «Компьютерное моделирование математических задач: учебное пособие» / Сулейманов Р.Р. — М.: бином. Лаборатория знаний, 2012

. «Введение в математическое моделирование: учебное пособие» / Под ред. Трусова П.В. — М.: университетская книга, Логос, 2007

. «Математические методы и модели в управлении» / Шикин Е. В., Чхартишвили А. Г. — М: Дело, 2002

. «Математические методы и модели исследования операций»: Учебное пособие / Кутузов А. Л. — издательство СПб ГПУ, 2005

. «Математические методы: Учебник» / Партика Т. Л., Попов И. И. — М: форум: ИНФРА, 2005

. «Математическое программирование» / Костевич Л., издательство «Новое знание», 2003

. «Экономико-математическое моделирование: практическое пособие по решению задач» / Мадера А. Г. — М: ИЭУП, 2004

Учебная работа. Обоснование принятия оптимальных решений для хлебозавода