- •Методы оптимальных решений
- •По условию: (1)
- •1. Линейное программирование
- •Решение
- •2. На плоскости Оx1x2 построить область допустимых решений и найти оптимальное решение (оптимальный план производства) геометрическим методом.
- •Задача №3.
- •7. Для определения нового (третьего) опорного плана формируем симплекс-таблицу 2.4.
Задача №3.
В транспортной задаче составить первоначальный план поставок и найти оптимальный план
Решение
С учетом начального условия :
Таблица 3.1.
Поставщики |
Мощность поставщиков |
Потребители и их спрос |
||
В1 |
В2 |
В3 |
||
190 |
120 |
40 |
||
А1 |
100 |
4 |
2 |
4 |
А2 |
200 |
2 |
5 |
3 |
А3 |
60 |
1 |
5 |
6 |
Тогда задача водится к решению транспортной задачи вида
.
Для решения транспортной задачи необходимо
1) составить ее экономико-математическую модель и
2) найти оптимальное распределение поставок и минимальные затраты на перевозку.
Первоначальное распределение выполним методом наименьших затрат.
1. Экономико-математическая модель формулируется следующим образом: необходимо найти объемы перевозок для каждой пары «поставщик-потребитель» так, чтобы:
мощности всех поставщиков были реализованы;
спросы всех потребителей были удовлетворены;
суммарная стоимость затрат была минимальна.
1.
Искомый объем перевозки от i-го
поставщика к j-му
потребителю обозначим
и назовем поставкой клетки (ij).
Заданные мощности поставщиков и спросы потребителей накладывают ограничения на значения . Для реализации мощностей каждого поставщика и удовлетворения спроса потребителей составляем уравнения баланса для каждой строки и для каждого столбца таблицы поставок 3.2:
и
с
ограничениями
т.к. объем перевозок не может быть
отрицательным.
Таким образом, экономико-математическая модель задачи составлена.
2. Первоначальное распределение определяется методом потенциалов.
По
условию задачи (см. таблицу 3.2) транспортная
задача является открытой, т.к. суммарные
мощности поставщиков
больше суммарного спроса потребителей
и часть поставщиков останется
незагруженными.
Для
баланса спроса и предложения и для
приведения системы к закрытому типу
вводим фиктивного потребителя
с нулевыми условными затратами на
перевозку и объемом спроса, равным
разности объемов поставок и потребления
.
Таблица 3.2 – таблица поставок.
Находим в таблице поставок клетки с наименьшей ненулевой стоимостью затрат: (3.1)) и анализируем возможности поставок для этой клетки:
для
клетки (3.1):
;
для
клетки (1.4):
;
для
клетки (1.1):
;
Минимальную
из максимально возможных поставок 60
для потребителя
даем в клетку (3.1), поставку 190-60=130 даем
в клетку (2.1) и исключаем из рассмотрения
третью строку и первый столбец.
В
оставшейся части таблицы поставок
наименьшими затратами, равными 2, обладает
клетка (1.2) с максимальной возможностью
поставок 100. Отдаем поставку 100 для
потребителя
в клетку (1.2), оставшуюся часть поставки
120-100=20 даем в клетку (2.2) и исключаем из
рассмотрения второй столбец
В
оставшейся части таблицы наименьшие
затраты 3 с максимальной (остаточной)
возможностью поставок 200-130-20=50 имеет
клетка (3.2). Для удовлетворения потребителя
отдаем поставку 40 в клетку (3.2) и исключаем
из рассмотрения 3 столбец.
При
этом поставщик
останется недогруженным в объеме
50-40=10 единиц.
Таблица поставок
Минимальные затраты на перевозку согласно таблицы поставок составят
Ответ: Минимальные затраты на перевозку согласно таблицы поставок составят 740 у.е.
4. A={aij}‑ матрица прямых материальных затрат, y – вектор конечного выпуска. Требуется:
1) Построить таблицу межотраслевого баланса в стоимостном выражении.
2) Найти изменение валовых выпусков при увеличении конечного выпуска первой отрасли на 20%, третьей – на 25% и неизменном выпуске второй отрасли.
№ |
а11 |
а12 |
а13 |
а21 |
а22 |
а23 |
а31 |
а32 |
а33 |
y1 |
y2 |
y3 |
22 |
0,1 |
0,2 |
0,1 |
0,2 |
0,1 |
0,0 |
0,0 |
0,2 |
0,1 |
100 |
200 |
300 |
Задача об инвестировании предприятий. Капитал в 5 условных д.е. требуется распределить между четырьмя предприятиями. Номер варианта совпадает с номером студента в списке группы. каждый вариант представляет собой таблицу размера 54. Строки таблицы соответствуют размеру инвестиций х, которые может получить предприятие, х=1, 2, 3, 4 или 5 условных д.е. соответственно. Столбцы таблицы соответствуют прибыли, котрую принесут
1-е, 2-е, 3-е и 4-е предприятие соответственно, при таком объеме инвестиций. (строку с нулевой прибылью, соответствующей х=0, добавить самостоятельно).
Найти распределение инвестиций, обеспечивающее максимальную суммарную прибыль от четырех предприятий вместе
Р
Систему уравнений, представляющих собой соотношения баланса, можно записать в матричном виде:
X = AX + Y,
где A – матрицей прямых затрат, X – вектор валового выпуска, Y – вектор конечного продукта.
Основная цель межотраслевого баланса состоит в том, чтобы отыскать такой вектор X, который при известной матрице прямых затрат A обеспечивает заданный вектор конечного продукта Y.
Система уравнений межотраслевого баланса легко решается методом обратной матрицы. Действительно,
X = AX + Y EX = AX + Y EX - AX = Y (E – A)X = Y
X
= (E – A)-1
Y
(E-A)-1 (E – A)X = (E – A)-1 Y EX = (E – A)-1 Y .
где матрица S = (E – A)-1 называется матрицей полных затрат.
В соответствии с введенными обозначениями имеем:
,
В задаче требуется найти такой вектор X, который при матрице А дал бы вектор конечного продукта.
Найдем матрицу E – A:
Определитель этой матрицы:
.
Найдем алгебраические дополнения Aij к элементам матрицы E – A:
,
,
,
,
,
,
,
,
Присоединенная матрица, т.е. транспонированная матрица алгебраических дополнений:
Тогда обратная матрица (E – A)-1 имеет вид:
Находим искомый вектор X:
.
Проверим правильность нахождения обратной матрицы и вектора Х, применяя функции Microsoft Office System Professional 2003 «МОПРЕД», «МОБР», «МУМНОЖ»:
Таким образом, объем валового продукта для 1 отрасли составит 214,8 условных денежных единиц, 2 отрасли – 269,96 условных денежных единиц, а 3 отрасли – 393,32 условных денежных единиц.
Для построения таблицы межотраслевого баланса в стоимостном выражении необходимо определить величины xij – часть объема продукции i-той отрасли, потребляемая j-той отраслью
,
,
,
,
,
,
,
,
.
Таблица межотраслевого баланса в стоимостном выражении
Отрасль |
Потребление |
Конечный продукт |
Валовой продукт |
|||
1 |
2 |
3 |
||||
Производство |
1 |
21,49 |
53,99 |
39,33 |
100 |
214,80 |
2 |
42,98 |
26,99 |
0,0 |
200 |
269,96 |
|
3 |
0,0 |
53,99 |
39,33 |
300 |
393,32 |
|
Данная таблица означает, что 1 отрасль производит 214,8 условных денежных единиц продукции, из них тратятся на нужды 1 отрасли 21,49 у.е., на нужды 2 отрасли 53,99 у.е., на нужды 3 отрасли 39,33 у.е. и 100 у.е. идут на потреблении; 2 отрасль производит 269,96 у.е. продукции, из них тратится на нужды 1 отрасли 42,98 у.е., на нужды 2 отрасли 26,99 у.е. и 200 у.е. идут на потребление; 32 отрасль производит 393,32 у.е. продукции, из них тратится на нужды 2 отрасли 53,99 у.е., на нужды 3 отрасли 39,33 у.е. и 300 у.е. идут на потребление.
2)
Для того, чтобы вычислить
изменение валовых выпусков при увеличении
конечного выпуска первой отрасли на
20%, третьей – на 25% и неизменном выпуске
второй отрасли,
необходимо найти такой вектор
,
который при матрице А дал бы вектор
конечного продукта
.
.
Отрасль |
Потребление |
Конечный продукт |
Валовой продукт |
Абсолютное изменение ВП |
Относительное изменение ВП |
|||
1 |
2 |
3 |
||||||
Производство |
1 |
21,49 |
53,99 |
39,33 |
120 |
248,11 |
33,31 |
15,51% |
2 |
42,98 |
26,99 |
0,0 |
200 |
277,36 |
7,40 |
2,74% |
|
3 |
0,0 |
53,99 |
39,33 |
375 |
478,30 |
84,98 |
21,61% |
|
Таблица изменений межотраслевого баланса в стоимостном выражении
Ответ:
валовой выпуск увеличится в 1 отрасли
на 15,51%, во 2 отрасли – на 2,74%, а в 3 отрасли
– на 21,61% и составит
Задача №2. Решить задачу симплекс-методом.
Для производства двух видов изделий - изделий А и изделий В - предприятие использует три вида сырья.
Нормы расхода сырья каждого вида на изготовление единицы продукции приведены в таблице. В ней же указана прибыль от реализации одного изделия каждого вида и общее количество сырья данного вида, которое может быть использовано предприятием.
Учитывая, что изделия А и В могут производиться в любых соотношениях (спрос обеспечен), требуется составить такой план их выпуска, при котором прибыль предприятия от реализации всех изделий максимальна.
Таблица 2.1.
Вид сырья |
Нормы расхода сырья (кг) на одно изделие |
Общее количество сырья (кг) |
|
А |
В |
||
1 |
12 |
6 |
300 |
2 |
4 |
6 |
180 |
3 |
4 |
12 |
240 |
Прибыль от реализации одного изделия (руб.) |
30 |
40 |
|
Решение.
Данные, представленные в таблице 1 удобно представить в виде
Задачу можно решить симплекс-методом ЗЛП. Для этого сформулируем задачу и запишем условие в каноническом виде.
1. Из анализа условия задачи следует, что целевая функция ЗЛП имеет вид
Q = 30x1 + 40x2 max (1)
при ограничениях:
(2)
2.
Приводим ЗЛП к каноническому виду, вводя
по числу неравенств неотрицательные
переменные x5
, x6
и x7
со знаком «+», так как все неравенства
системы (1)
:
Q = 30x1 + 40x2 + 0x3 + 0x4 + 0x5 max (3)
при ограничениях:
(4)
Для нахождения начального опорного плана нужно найти общее решение системы ограничений, но такое, чтобы все свободные члены (4) i были неотрицательными.
Второе уравнение имеет предпочтительный вид, в то время как остальные уравнения его не имеют. Для нахождения допустимого общего решения системы ограничений применим метод Жордана. Удобно это делать, записав систему ограничений в виде симплекс-таблицы.
Для
нахождения общего решения методом
Жордана необходимо на каждом шаге
выбирать разрешающий элемент. Чтобы
соблюсти при этом условие неотрицательности
свободных членов, для каждого положительного
коэффициента aij
рассчитываем
.
При этом, если коэффициент aij £ 0, то qij не вычисляем, а в соответствующей клетке таблицы ставим прочерк. Затем произвольно выбираем разрешающую строку (кроме строки, соответствующей уравнению в предпочтительном виде), а разрешающим столбцом выбираем тот, в котором в разрешающей строке находится наименьшее в этом столбце значение qij.
3.
За базисные (основные) переменные
принимаем x5
; x6
; x7
, исходные переменные – за свободные
переменные, причем x1
= 0; x2
= 0; x3
= 0; x4
= 0. Выражая
через свободные переменные
(5)
и подставляя значения x1= 0; x2 = 0; x3 = 0 в систему (5) получим x3 =300; x4 =180; x5 =240. Таким образом, первый опорный план имеет вид:
,
а
функция
.
4. Для проверки оптимальности начального опорного плана составляем симплекс-таблицу 2.2.
Симплекс-таблица 2.2
Базис |
cj |
30 |
40 |
0 |
0 |
0 |
|
|
||
xi
|
i |
ci |
bi |
ai1 |
ai2 |
ai3 |
ai4 |
ai5 |
||
x3 |
3 |
0 |
300 |
12 |
6 |
1 |
0 |
0 |
322 |
50 |
x4 |
4 |
0 |
180 |
4 |
6 |
0 |
1 |
0 |
195 |
30 |
|
5 |
0 |
240 |
4 |
12 |
0 |
0 |
1 |
262 |
20 |
j |
0=0 |
– 30 |
– 40
|
0 |
0 |
0 |
779 |
|
||
Для
базисных неизвестных величины j
= 0. Поэтому подсчитываем
0
= Q
и 1,2,3,4,5.
:
,
Q (Х1)= 0 = 300 0 + 180 0 + 240 0 = 0;
1 = 012 + 0 4 + 0 4 – 30 = – 30 ;
2 = 06 + 0 6 + 0 12 – 40 = – 40 ;
3 = 01 + 0 0 + 0 0 – 0 = 0;
4 = 00 + 0 1 + 0 0 – 0 = 0.
5 = 00 + 0 0 + 0 1 – 0 = 0.
Поскольку в индексной строке j имеются отрицательные числа 1=–30, 2 = – 40, то начальный опорный план не является оптимальным.
5. Находим направляющий столбец и направляющую строку.
Направляющему
столбцу соответствует максимальный по
модулю отрицательный элемент индексной
строки
.
Следовательно, переменную
следует перевести из свободных в
базисные.
Направляющей
строке соответствует минимальное
значение
:
.
Следовательно, переменную
следует перевести из базисных в свободные.
Разрешающим
является элемент
.
Направляющие строка и столбец обозначены
в симплекс-таблице 2.1 стрелками.
6. Для определения нового опорного плана формируем симплекс-таблицу 2.3.
На
место переменной
в столбце «Базис,
» ставим переменную
,
в разрешающую строку в столбце «Базис,
»
записываем значение
коэффициента целевой функции при
переменной
,
остальные элементы строки
получаются делением элементов строки
на разрешающий элемент
.
В
симплекс-таблице 2.3 элемент
,
а остальные элементы разрешающего
столбца полагаем равными 0. Оставшиеся
элементы строк
вычисляем по правилу прямоугольника
.
Симплекс-таблица 2.3
Базис |
cj |
30 |
40 |
0 |
0 |
0 |
|
j*=1 |
||
xi
|
i |
ci |
bi |
ai1 |
ai2 |
ai3 |
ai4 |
ai5 |
||
x3 |
3 |
|
180 |
10 |
0 |
1 |
0 |
-0,5 |
190,5 |
18 |
x4 |
4 |
|
60 |
2 |
0 |
0 |
1 |
-0,5 |
62,5 |
30 |
x2 |
2 |
40 |
20 |
0,3 |
1 |
0 |
0 |
0,08 |
63,38 |
60 |
j |
0=800 |
– 26,(7)
|
0 |
0
|
0 |
3,(3) |
316,38 |
|
||
Новый
опорный план имеет вид:
.
Для
новых базисных неизвестных
величины 3,4,2
= 0. Подсчитываем 0
= Q;
j
, где
,
:
Q (Х2)= 0 =0 180 + 060+ 4020= 800;
1 = 010 + 0 2 + 401/3 – 30 = – 16,(7) ;
2 = 00 + 0 0 + 40 1 – 40 = 0 ;
3 = 01 + 0 0 +4 0 0 – 0 = 0;
4 = 00 + 0 1 + 40 0 – 0 = 0.
5 = 0(-1/12) + 0 (-1/12) + 40 1/12 – 0 = 3,(3).
В индексной строке имеется отрицательное число 1 = – 16,(7), поэтому полученный опорный план не является оптимальным.
