Файл: Matematicheskoe_modelirovanie-2011-2.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.09.2021

Просмотров: 415

Скачиваний: 4

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

НЭ=СЭ – (А*В)/РЭ (6)

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

Пример №2.

На приобретение оборудования для производственного участка выделено 20тыс. руб. Оборудование может быть размещено на площади, не превышающей 72 кв.м. Можно заказать оборудование двух типов: типа А, требующие производственную площадь 6кв.м и дающие 6 тыс.ед. продукции в смену ( цена 5000 руб.) и типа В, требующие площадь 12 кв.м и дающие 3тыс.ед., (цена 2000 руб.). Каков оптимальный план приобретения оборудования, обеспечивающий максимальную производительность участка?

Обозначим количество приобретаемого оборудования типа А и В через Х1 и Х2 соответственно.

Производительность участка (целевая функция) : Z(X) =6Х1+3Х2.

Основные ограничения связаны

с денежными средствами : 5Х1+2Х2 ≤ 20,

с площадью производственного участка : 6Х1+12Х2 ≤ 72.

Вводим новые базисные переменные Х3 (остаток денежных средств после закупки оборудования) и Х4 (остаток площадей после размещения оборудования) и перепишем ограничения в виде системы уравнений:

5X1+2Х2+X3=20 (X3=20 – 5X1 - 2X2)

1+12Х2+X4 = 72 (X4=72 – 6X1 – 12X2)

При этом функция цели: Z(X) =6Х1+3Х2+0Х3+0Х4.

Составляем опорный (0-ой) план: Х= (0, 0, 20, 72), т.е. пока ничего не приобретено (деньги не потрачены, площади пустуют). Составляем симплексную таблицу

План

Базис

Ci/Cj

Знач. Xi

X1

X2

X3

X4

Qmin

0

X3

0

20

5

2

1

0

20/5=4

X4

0

72

6

12

0

1

72/6=12

Z(X) = 0

- 6

- 3

0

0

Индексная строка

1

X1

6

4

1

0,4

0,2

0

4/0,4=10

X4

0

48

0

9,6

-1,2

1

48/9,6=5

Z(X) = 6*4=24

0

-0,6

1,2

0

Индексная строка

2

X1

6

2

1

0

0,25

-1/24

-

X2

3

5

0

1

-1/8

5/48

-

Z(X) =6*2+3*5=27

0

0

9/8

1/16

Индексная строка















Очевидно, что ведущий столбец соответствует Х1, так как имеет самый большой индекс 6. Находим минимальное значение Qmin=4 ( самое жесткое ограничение ресурса), определяя ведущую строку, показывающую, что из базисных переменных выводится Х3, а вместо нее вводится Х1. Пересчитываем элементы ведущей строки, разделив их на 5, а по формуле (6) определяем элементы второй и индексной строк. Целевая функция для 1-ого плана равна Z(X) =6*4+3*0=24.Однако, один из коэффициентов индексной строки для столбца Х2 остается отрицательным -0,6, следовательно данный план не оптимален, и требуется еще одна итерация (шаг) для его улучшения. Выбираем ведущим 2-ой столбец и по минимальному значению Qmin=5 определяем ведущую строку с базисной переменной Х4. Выполнив те же преобразования, получаем 2-ой план, который будет оптимальным, так как все индексные коэффициенты положительны. Проанализируем полученные результаты. При оптимальном решении целевая функция имеет максимальное значение 27тыс.руб., при этом оба ресурса выведены из базиса, следовательно израсходованы полностью! Убедимся в этом : 5*2+2*5=20 тыс.руб., 6*2+12*5=72 кв.м. Искомое решение Х= (2; 5;0;0;).Так бывает далеко не всегда.


Пример №3.

Коммерческое предприятие реализует три вида продукции А,В, и С. Нормативы затрат материально-денежных ресурсов на реализацию единицы продукции (или 1 тыс. товарооборота) представлены в таблице.

Виды ресурсов

Нормы затрат ресурсов на 1тыс. руб. товарооборота

Объем ресурсов

Группа А

Группа В

Группа С

Рабочее время продавцов, чел.-ч

0,1

0,2

0,4

1100

Площадь торговых залов, м2

0,05

0,02

0,02

120

Площадь складских помещений, м2

3

1

2

8000

Доход, тыс. руб.

3

5

4




В последней строке задан доход от продажи товара каждой группы. Необходимо определить объем и структуру товарооборота, при котором доход предприятия будет максимальным. Математическая модель задачи: определение вектора Х=(Х1; Х2; Х3), где Х – количество проданных изделий А, В,С, который удовлетворяет неравенствам



0,1X1 + 0,2X2 +0,4 X3 ≤ 1100,

0,05X1 + 0,02X2 +0,02 X3 ≤ 120,

3X1+1X2+2X3≤8000. X1≥0 X2≥0 X3≥0

и обеспечивает максимальное значение целевой функции Z(X)=3Х1+5Х2+4Х3.

Для построения опорного 0-ого плана вводим дополнительно базисные переменные Х4, Х56 (определяем объемы неиспользованных ресурсов), преобразовывая систему неравенств в систему уравнений,

0,1X1 + 0,2X2 +0,4 X34= 1100,

0,05X1 + 0,02X2 +0,02 X35= 120,

3X1+1X2+2X36=8000.

решаем полученную систему относительно базисных переменных Х4, Х5, Х6:

Х4=1100-0,1X1 - 0,2X2 -0,4 X3 ,

Х5= 120-0,05X1 - 0,02X2 -0,02 X3,

Х6=8000-3X1 -X2-2X3.

Функцию цели запишем в виде Z(X)=0 – (- 3Х1 – 5Х2 – 4Х3) и полагая, что основные переменные Х123=0, получим опорный план Х=(0; 0; 0; 1100; 120; 8000), при котором целевая функция Z(X)=0 (товары не продаются, ресурсы не используются, доход нулевой). Опорный план запишем в симплексную таблицу:

План

Базис

Ci/Cj

Знач. Xi

X1

X2

X3

X4

X5

X6

Qmin

0

X4

0

1100

0,1

0,2

0,4

1

0

0

5500

X5

0

120

0,05

0,02

0,02

0

1

0

6000

X6

0

8000

3

1

2

0

0

1

8000

Z(X) = 0

- 3

- 5

- 4

0

0

0

Индекс строка

1

X2

5

5500

0,5

1

2

5

0

0

11000

X5

0

10

0,04

0

-0,02

-0,1

1

0

250

X6

0

2500

2,5

0

0

-5

0

1

1000

Z(X) = 27500

- 0,5

0

6

25

0

0

Индекс строка

2

X2

5

5375

0

1

2,25

6,25

-12,5

0

-

X1

3

250

1

0

-0,5

-2,5

25

0

-

X6

0

1875

0

0

1,25

1,25

-62,5

1

-

Z(X) = 27625

0

0

5,75

23,75

12,5

0

Индекс строка




Естественно, 0-ой план не оптимален, так как коэффициенты в индексной строке отрицательные -3, -5, -4.Ведущий столбец содержит Х2, а ведущая стока содержит Х4, РЭ=0,2 находится на их пересечении. Формируем 1-ый план аналогично предыдущему примеру, в индексной строке только коэффициент при Х1 отрицателен -0,5, поэтому ведущий столбец содержит Х1. Строка, содержащая Х5, становится ведущей , и после итерации получаем 2-ой план, все коэффициенты которого в индексной строке положительны. Это означает, что план Х= (250; 5375; 0; 0; 0; 1875) оптимален, и целевая функция при данных значениях Х имеет максимальное значение Z(X) =27625 тыс.руб.

Проанализируем полученное решение: для получения максимальной прибыли необходимо продать 250 ед. товаров группы А, 5375 ед. товаров группы В, товары группы С вообще не реализуются! В оптимальном плане среди базисных переменных находится Х6. Это указывает на то, что 1875м2 складских помещений не используются, то есть данный ресурс в исходном объеме избыточен, и данная его часть (1875м2) может быть применена в других целях.

Метод искусственного базиса.

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

ai1 X1 + ai2 X2 +…ain X n ≥ bi (7)

или уравнений ai1 X1 + ai2 X2 +…ain X n = bi (7*),

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

ai1 X1 + ai2 X2 +…ain X n +Yi = bi (8).

Целевая функция при решении задачи на максимум запишется в виде:

Z(X) =∑CjXj+(-M)∑Yi (9),

при решении аналогичной задач на минимум :

Z(X)=∑CjXj+(M)∑Yi (9*),

где М – очень большое положительное число, своего рода штраф за использование искусственных переменных.

В случае неравенств (7) первоначально вводим дополнительные переменные Хn+i со знаком минус. Их матрица не будет единичной, поэтому в каждое неравенство системы (7) вводим искусственные переменные Уi:

ai1X1+ai2X2+…ainXnXn+i+Yi=bi (10) Целевая функция при этом Z(X)=∑CjXj+0∑Xn+i+(-M)∑Yi (для нахождения максимума). Применение искусственного базиса придает симплексному методу большую гибкость и позволяет использовать его для широкого круга задач.

Пример №4.

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


Виды ресурсов

Нормы затрат на производство 1 т

Запасы ресурсов

Группа А

Группа В

Сырье ,т

1,0

1,0

6

Рабочее время чел.-час

2,0

1,0

8

Доход с 1 т, тыс. руб.

3

2


Математически ограничения выпуска продукции запишутся в виде смешанной системы:

1 + 1Х2≤ 6,

1 + 1Х2 =8.

Введем для первого неравенства базисную переменную Х3, а для второго уравнения искусственную переменную Y1:

1 + 1Х2+ Х3 = 6,

1 + 1Х2 +Y1 =8.

Выразим из полученной системы уравнений Х3 и Y1 и для определения максимума целевую функцию представим:

Z(X)= 3X1+ 2X2+0X3 –MY1= 3X1+ 2X2 –M(8 -2X1 –X2)=

= 3X1+ 2X2 –8M +2MX1 + MX2 = (2M + 3)X1 + (M + 2)X2 -8M

Для опорного плана - Х=(0,0,6,8). Построим симплексную таблицу:

План

Базис

Ci/Cj

Знач. Xi

X1

X2

X3

Y1

Qmin

0


X3

0

6

1

1

1

0

6/1=6

Y1

-M

8

2

1

0

1

8/2=4

Z(X) = -8M

-2M-3

-M-2

0

0

Индексная строка

1

X3

0

2

0

0,5

1

-0,5

2/0,5=4

X1

3

4

1

0,5

0

0,5

4/0,5=8

Z(X) = 3*4=12

0

- 0,5

0

М+1,5

Индексная

строка

2

X2

2

4

0

1

2

-1

-

X1

3

2

1

0

-1

1

-

Z(X) =3*2+2*4=14

0

0

1

М+1

Индексная строка



Как правило, улучшение опорного плана начинается с выведения из базиса искусственной переменной Y1 .Оптимальный план Х=(2,4,0,0) получен на второй итерации, при этом доход максимален 14тыс. руб. , а коэффициенты индексной строки неотрицательны. Легко убедиться, что в данной задаче при оптимальном плане ресурсы использованы полностью (2*1+4*1=6; 2*2+1*4=8).

При нахождении минимальной доходности иначе формулируем целевую функцию ( в качестве слагаемого вводится +MY1 :

Z(X)= 3X1+ 2X2+0X3 +MY1= 3X1+ 2X2 +M(8 -2X1 –X2)=

= 3X1+ 2X2 +8M - 2MX1 - MX2 = (3 - 2M)X1 + (2 - M )X2 +8M

Опорный план тот же , но коэффициенты индексной строки в симплексной таблице иные. Ведущий столбец, по-прежнему, выбираем по наибольшему по абсолютному значению положительному коэффициенту при X1 , ведущая строка определяется по минимальному значению Qmin=4.При первой итерации из базиса выводится искусственная переменная Y1.





План

Базис

Ci/Cj

Знач. Xi

X1

X2

X3

Y1

Qmin

0

X3

0

6

1

1

1

0

6/1=6

Y1

M

8

2

1

0

1

8/2=4

Z(X) = 8М

2M-3

M-2

0

0

Индексная строка

1

X3

0

2

0

0,5

1

-0,5

2/0,5=4

X1

3

4

1

0,5

0

0,5

4/0,5=8

Z(X) = 3*4=12

0

- 0,5

0

-М+1,5

Индексная строка










Полученные отрицательные значения коэффициентов в индексной строке Xi свидетельствуют об оптимальности 1-ого плана, при этом минимальный доход 12 тыс. рублей.

Он обеспечивается только выпуском продукции А (продукция В не выпускается), сырье не используется полностью (остаток Х3=2т), при этом выполнено основное условие - рабочие полностью заняты на производстве.

Пример №5.

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

Виды ресурсов

Нормы затрат на производство 1 т

Запасы ресурсов

Группа А

Группа В

Сырье ,т

3,0

2,0

12

Доход с 1 т, тыс. руб.

2

1


Математически ограничения выпуска продукции запишутся в виде системы неравенств:

3X1 + 2X2 ≤ 12,

X1 + X2 ≥ 5.

Введем для первого неравенства базисную переменную Х3, а для второго уравнения отрицательную Х4 и искусственную переменную Y1:

3X1 + 2X2 + X3 = 12,

X1 + X2X4 + Y1 = 5.

Выразим из полученной системы уравнений Х3 и Y1 и для определения максимума целевую функцию представим:

Z(X) = 2X1+ 1X2+0X3 +0Х4MY1=2X1+ 1X2M*(5 - X1X2 + Х4) =

= (2 + М)X1+ (1 + М)X2MХ4 - 5M

Для нулевого опорного плана Х=( 0,0,12,5) построим симплексную таблицу:

План

Базис

Ci/Cj

Знач. Xi

X1

X2

X3

Х4

Y1

Qmin

0

X3

0

12

3

2

1

0

0

12/3=4

Y1

-M

5

1

1

0

-1

1

5/1=5

Z(X) = - 5М

-M-2

-M-1

0

М

0

Индекс строка

1

X1

2

4

1

2/3

1/3

0

0

4/(2/3)=6

Y1

-M

1

0

1/3

-1/3

-1

1

1/(1/3)=3

Z(X) = 2*4 - M=8-M

0

-M/3+1/3

(M+2)/3

M

0

Индекс строка

2

X1

2

2

1

0

1

2

-2

-

X2

1

3

0

1

-1

-3

3

-

Z(X) = 2*2+1*3=7

0

0

2/3M+1

1

М-1

Индекс строка



Полученный план Х=(2,3,0,0), при этом доход максимален 7тыс. руб., ресурсы использованы полностью (2*3+3*2=12) и выполнено условие маркетинга (2*1+3*1=5).

Получение целочисленных решений

при решении оптимизационных задач.

Метод Гомори.

Значительная часть задач хозяйственной и коммерческой деятельности требует целочисленного решения, например, выпуск и распределение товаров, использование агрегатов при загрузке оборудования и т.д. Как быть, если симплекс-метод дает при оптимальном решении дробные значения основных переменных? В этом случае идет поиск наиболее близкого к оптимальному целочисленного решения. Применяются различные методики: 1) отсечения, 2) комбинированные, 3) приближенные. К числу первых относится известный метод Гомори.