Прикладная математика. Задача 1 (Линейная производственная задача) Задача 2 Двойственная задача
Скачать 439 Kb.
|
Т аблица 7.
Таблица 9.
Таблица 10.
Таблица 11.
Проверим выполнение равенства f1(x*1) + f2(x*2) + f3(x*3) + f4(x*4) = z max 58 + 49 + 22 + 68 = 197 Равенство выполняется. Задача №5. Динамическая задача управления производством и Рассмотрим трехэтапную систему управления запасами с дискретной продукцией и динамическим детерминированным спросом. Пусть спрос (заявки) потребителей на нашу продукцию составляют: на первый этап d1=2 единицы, на второй – d2=3, на третий - d3=3 единицы. К началу первого этапа на складе имеется только 2 единицы продукции, т.е. начальный уровень запаса равен y1=2. Затраты на хранение единицы продукции на разных этапах различны и составляют соответственно h1=3, h2=2, h3=2. Затраты на производство xj единиц продукции на j-м этапе определяются функцией j(xj) =2xj2 + 3xj + 4 (1) т.е. а=2; b=3; с=4. Требуется указать, сколько единиц продукции на отдельных этапах следует производить, чтобы заявки потребителей были удовлетворены, а наши общие затраты на производство и хранение за все три этапа были наименьшими. Решение. Представим исходные данные задачи в виде таблицы: Таблица 12.
Воспользовавшись рекуррентными соотношениями, последовательно вычисляем F1 ( = y2), F2 ( = y3), F3 ( = y4), и соответственно находим 1 (= y2), 2 ( = y3 ), 3 ( = y4) Положим k = 1. Тогда имеем (2) Учтем, что параметр состояния = у2 может принимать целые значения на отрезке 0 у2 d2 + d3 (3) 0 y2 3 + 3 т.е. у2 = 0, 1, 2, 3, 4, 5, 6. При этом каждому значению параметра состояния должна отвечать определенная область изменения переменной x1, а именно 0 х1 2 + у2 (4) На первом этапе спрос d1 = 2, исходный запас у1 = 2. Из балансового уравнения х1 + у1 - d1 = у2 (5) непосредственно следует, что объем производства связан со значением параметра состояния = у2 соотношением x1 = y2 + d1 - y1 = y2 + 2 - 2 = y2 (6) Если задан уровень запаса к началу первого этапа, то каждому значению у2 отвечает единственное значение х1 и потому F1( = y2) = 1 (x1, y2) = ах12 + bx1 + c + h1y2 (7) П y2 = 0, x1 = 0, 1 (0;0) = 2 ∙ 02 + 3 0 + 4 + 3 0 = 4 y2 = 1, x1 = 1, 1 (1;1) = 2 ∙ 12 + 3 1 + 4 + 3 1 = 12 y2 = 2, x1 = 2, 1 (2;2) = 2 ∙ 22 + 3 ∙ 2 + 4 + 3 ∙ 2 = 24 y2 = 3, x1 = 3, 1 (3;3) = 2 ∙ 32 + 3 ∙ 3 + 4 + 3 ∙ 3 = 40 y2 = 4, x1 = 4, 1 (4;4) = 2 ∙ 42 + 3 ∙ 4 + 4 + 3 ∙ 4 = 60 y2 = 5, x1 = 5, 1 (5;5) = 2 ∙ 52 + 3 ∙ 5 + 4 + 3 ∙ 5 = 84 y2 = 6, x1 = 6, 1 (6;6) = 2 ∙ 62 + 3 ∙ 6 + 4 + 3 ∙ 6 = 112 Значения функции состояния F1( ) представлены в табл. 13. Таблица 13.
Переходим ко второму этапу. Полагаем k = 2 и табулируем функцию F2( = y3). (8) Здесь минимум берется по единственной переменной х2, которая может изменяться в пределах 0 x2 d2 + y3 или 0 x2 3 + y3 (9) где верхняя граница зависит от параметра состояния = у3, который, принимает значения на отрезке 0 y3 d3 , т.е. 0 y3 3 (10) а аргумент у2 в последнем слагаемом справа в соотношении (8) связан с х2 и у3 балансовым уравнением x2 + y2 - d2 = y3 откуда следует y2 = y3 + d2 - x2 = y3 + 3 - x2 (11) Придавая параметру состояния различные значения от 0 до 3, будем последовательно вычислять 2 (x2, ), а затем определять F2( ) и 2( ). Положим, например = у3 = 2. Тогда, согласно (9), 0 x2 5, т.е. переменная х2 может принимать значения: 0, 1, 2, 3, 4, 5, а каждому значению х2 отвечает определенное значение у2, вычисляемое по формуле (11): у2 = 5 - х2 Последовательно находим: если x2 = 0, то y2 = 5-0 = 5, 2 (0,2) = 2∙02 + 30 + 4 + 22 + F1(5) = 8 + 84 = 92, x2 = 1, y2 = 5-1 = 4, 2 (1,2) =2∙12 + 31+ 4 + 22 + F1(4) = 13 + 60 = 73, x2 = 2, y2 = 5-2 = 3, 2 (2,2) =2∙22 +32 + 4 + 22 + F1(3) = 22 + 40 = 62, x2 = 3, y2 = 5-3 = 2, 2 (3,2) =2∙32+33 + 4 + 22 + F1(2) = 35 + 24 = 59*, x2 = 4, y2 = 5-4 = 1, 2 (4,2) = 2∙42 + 34 + 4 + 22 + F1(1) = 52+12 = 64, x2 = 5, y2 = 5-5 = 0, 2 (5,2) =2∙ 52 + 35 + 4 + 22 + F1(0) = 73 + 4 =77 . Н F2 ( = y3 = 2) = min 2 (x2,2) = min (92, 73, 62, 59, 64, 77) = 59, x2 причем минимум достигается при значении х2, равном 2 ( = y3 = 2) = 3 Аналогично для всех остальных значений у3 = 0, 1, 3. F2 ( = y3 = 0) = 30; 2 ( = y3 = 0) = 2. F2 ( = y3 = 1) = 44; 2 ( = y3 = 1) = 2. F2 ( = y3 = 3) = 77; 2 ( = y3 = 3) = 3. Процесс табулирования функции F2 ( = y3) приведен в табл. 14, а результаты табулирования сведены в табл. 15. |