Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

08-МУ к лаб_занятиям

.pdf
Скачиваний:
74
Добавлен:
04.06.2015
Размер:
1.74 Mб
Скачать

11

Для этого надо нажать кнопку Параметры и заполнить некоторые поля окна

Параметры поиска решения (рисунок 8).

Рисунок 8 - Параметры поиска решения, подходящие для большинства задач линейного программирования

Параметр Максимальное время служит для назначения времени (в секундах), выделяемого на решение задачи.

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

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

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

Параметр Сходимость применяется при решении нелинейных задач. Установка флажка Линейная модель обеспечивает ускорение поиска

решения линейной задачи за счет применения симплекс-метода.

Запуск задачи на решение

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

После запуска на решение на экране появляется окно Результаты поиска решения, в котором решения представлены названия трех типов отчетов: "Результаты", "Устойчивость", "Пределы". Они необходимы при анализе полученного решения. Если при вводе условий задачи в Excel были допущены ошибки, то после запуска задачи на решение на экран будет выдано соответствующее сообщение с указанием причины, по которой решение не найдено.

12

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

Рисунок 9 - Экранная форма задачи (8) после получения решения

Порядок выполнения работы

1.Для задачи линейного программирования, соответствующей заданному варианту, найти оптимальное решение в табличном процессоре

Microsoft Excel

2.Оформить отчет по лабораторной работе, который должен содержать:

название работы;

исходные данные варианта;

результаты решения задачи.

Задание

 

 

 

Вариант 1

 

 

 

 

 

 

 

 

Вариант 2

 

 

F ( X ) = 5x1 + 7x2 − 6x3 + 9x4 + 8x5 → max;

F ( X ) = -45x1 + 65x2 + 2x4 - 3x5 ® max;

0,7x

+ 0,9x

+1,5x

+ 2,3x

+1,8x £ 50000,

15x

1

+18x

2

+ 34x

4

- 22x

5

= 56,

 

1

2

3

4

5

 

 

 

 

 

 

 

0,4x

1+1,1x2 - 0,5x3 +1,3x4 - 2,8x5 ³ 32000,

2x 1

+ 7x3 - 4x4 + 3x5 ³ 91,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,5x 1+1,8x3 + 0,7x4 + 2x5 £ 40000,

0,2x 1 + 0,8x2 +1,5x3 + 0,9x4 + 4x5 £ 26,

2,2x

1-1,4x2 - 0,8x3 + 0,9x4 = 15000,

1,8x

1

- 42x

2

+ 6,4x

3

+ 3x

5

= 15,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xj ³ 0 (j = 1,5).

 

 

x j ³ 0 (j = 1,5).

 

 

 

 

 

13

Вариант 3

F(X ) = 46x1 + 2,3x2 + 9,4x3 − 4x5 → max;

3x

1

+ 7,8x +12x

+ 9x ³ 49,

 

 

 

3

4

 

5

2,3x2 + 5x3

+ 5,6x4 - x5 £ 86,

 

 

 

- 40x4 + 29x5

= 50,

16x1

190x

 

-98x - 4x

+150x ³ 300,

 

 

 

1

 

2

 

4

5

 

 

 

 

 

 

 

 

 

³ 0 (j =1,5).

 

 

xj

 

 

Вариант 5

F(X ) = 10x1 + 40x3 + 13x4 + 56x5 → min;

7x

 

+16x

 

 

+ 5x

+ 25x

 

£ 600,

 

1

 

3

 

 

4

 

 

5

 

 

 

 

 

 

 

 

 

 

8x

1+1,7x2 - 0,5x4 + 4,7x5 = 890,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6x 1+ 4x3 - 7x4 + 6,3x5 £ 270,

84x

 

+ 62x + 80x

+14x ³ 2300,

 

 

1

 

 

 

 

2

 

 

 

 

 

3

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j =1,5).

 

 

 

 

 

 

xj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 7

 

F ( X ) = x1 + 4x3 + 8x4

-12x5 ® min;

x

1

+

9x

2

+ 2x

3

- 4x

4

= 250,

 

 

 

 

 

 

 

 

 

 

 

 

 

0,4x1 + x2 - 5x3 + 3x4

+ 8x5 £ 460,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,5x1 + 10x 2 - 8x3 + 6x4 + 2x5 £ 190,

 

 

 

 

-8,5x3 + 3x4 + 2x 5

= 210,

11x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j = 1,5).

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 9

 

F(X ) = 12x2 + 89x3 − 5x5 → max;

2x

1

+9,6x +15,7x + 22x -8x £ 73,

 

 

 

 

2

 

 

 

 

 

3

 

 

4

5

0,9x 1+11,1x2 - 4,3x3 +1,5x4 + 6,4x5 =19,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14x1 + 45x2 -38x4 + 26x5 £ 49,

220x

 

-150x + 3x +95x

=133,

 

 

 

1

 

 

 

2

 

 

3

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j =1,5).

 

 

 

 

 

 

xj

 

 

 

 

 

 

Вариант 4

F ( X ) = 14x1 - 9x2 - x4 + 6,4x5 ® min;

0,9x

1 +10x2

- 28x4 + 5x5 £ 245,

 

1 +1,7x2

 

- 0,2x3 - 0,5x4 = 9,

0,8x

 

 

+ 4x3 - 7x4 + 6,3x5 £ 54,

6x 1

 

 

 

 

 

8x 1 + 6,2x2 - 4,8x4 + 2,9x5 ³ 17,

 

 

 

 

 

 

 

 

 

x j ³ 0 (j = 1,5).

Вариант 6

F(X ) = 0,5x1 +1,8x3 − 9,2x4 +14x5 → min;

9,6x

 

+15,7x + 24x

 

-8x

£ 74,

 

 

 

2

 

 

 

 

 

 

 

 

3

 

 

 

 

 

4

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,8x 1+11,1x 2 - 4,5x3 +1,5x4 - 6,3x5 = 22,

 

 

 

 

 

 

 

 

 

 

-38x4 + 26x5 £ 46,

14x 1+ 45x 2

220x -148x

 

 

 

- 7x

 

 

+ 95x ³150,

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

3

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j =1,5).

 

 

 

 

 

 

 

 

 

 

 

xj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 8

 

 

 

 

F ( X ) = 4x1 + 6x2

− 14x3 + 49x5 → min;

21x

 

+ 9x

 

 

 

 

- 2x

 

-12x

³ 58,

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

4

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

110x2 - 60x3 + 80x4 - 45x5 = 290,

 

 

 

 

 

 

 

 

 

 

-14x4 + x5 £ 72,

5x2 + 27x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 130x4 = 140,

87x 1- 6,4x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j = 1,5).

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 10

 

 

 

F ( X ) = 84x1 + 5,7x2

 

+ 10x4

 

− 3x5 → max;

4x

 

+ 8,5x

 

 

 

+ 16x

 

 

+ 10x

³ 50,

 

1

 

 

 

 

 

2

 

 

 

 

 

 

3

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10,4x

1 + 6x3 + 2x4

 

+ 4x5

£ 120,

 

 

 

 

 

 

 

 

 

 

- 20x4

 

+ 30x5

 

= 600,

19x 1+18x2

 

 

200x

 

+ 45x

 

 

 

 

- 8x

 

 

+ 3,4x

 

 

³ 210,

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

³ 0 (j = 1,5).

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

Вариант 11

F(X ) = -38x 1 + 60x2 + 4x4 +8x5 ® max;

18x

1+ 4x2

+ 2x3

-12x5

£ 86,

 

+19x3

-7x4

+10x5

=130,

2x2

 

 

 

 

+ 2x4 -5x5 £ 34,

0,4x 1+ 3x2 - 4,2x3

 

 

 

 

+ 6x4 =18,

2,1x 1+13x2 - 20x3

x j ³ 0 (j =1,5).

14

Вариант 12

F(X ) = 0,84x2 − 4x3 + 3,8x4 +12x5 → min;

15x + 9,6x

+ 34x

- 8x £180,

 

1

2

 

 

4

5

0,6x 1+11,1x2 - 2,6x3 +1,5x4 - 6,3x5 = 68,

 

 

 

 

- 38x4 +12x5 £ 81,

14x 1+ 64x3

190x -148x - 7x

+ 84x ³ 230,

 

 

1

 

2

3

5

 

 

 

 

 

 

 

 

³ 0 (j =1,5).

 

 

xj

 

 

Контрольные вопросы

1.Как формулируется общая задача линейного программирования?

2.Как формулируется задача линейного программирования в стандартной форме?

3.Какова постановка канонической задачи линейного программирования?

4.Что называется допустимым решением задачи?

5.Какое решение называется оптимальным?

6.

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

8.Каким образом в MS Excel задается направление оптимизации целевой функции?

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

10.Поясните общий порядок работы с окном Поиск решения.

11.Каким образом можно изменять, добавлять, удалять ограничения в окне Поиск решения?

12.Какие сообщения выдаются в MS Excel в случаях: успешного решения задачи; несовместности системы ограничений задачи; неограниченности целевой функции?

13.Объясните смысл параметров, задаваемых в окне Параметры поиска решения.

Лабораторная работа №2 Графический анализ чувствительности оптимального решения задачи оптимизации плана производства

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

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

15

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

Выделяют следующие три задачи анализа на чувствительность.

1.Анализ сокращения или увеличения ресурсов:

1)на сколько можно увеличить запас дефицитного ресурса для улучшения оптимального значения целевой функции?

2)на сколько можно уменьшить запас недефицитного ресурса при сохранении полученного оптимального значения целевой функции?

2.Увеличение (уменьшение) запаса какого из ресурсов наиболее выгодно?

3.Анализ изменения коэффициентов целевой функции: каков диапазон изменения коэффициентов целевой функции, при котором не меняется оптимальное решение?

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

Этап 1. Сначала на координатной плоскости x1 0x2 строится область допустимых решений (ОДР) – совокупность точек, координаты каждой из которых составляют решение системы ограничений задачи. Далее строится

вектор-градиент линейной функции

 

F(X) в

какой-нибудь точке X 0 ,

принадлежащей ОДР:

 

 

 

 

 

 

 

Ñ =

 

F

= c1

 

F

= c2

 

 

,

 

 

 

 

x

x

.

 

 

 

 

 

 

 

1

 

 

2

 

Этап 2. Прямая (линия уровня)

c1 x1 + c2 x2 = F (X 0 ) (F (X ) = const ),

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

Этап 3. Для нахождения координат точки максимума надо решить совместно два уравнения прямых, получаемых из соответствующих ограничений и дающих в пересечении точку максимума. Значение F(X), найденное в получаемой точке, является максимальным.

В случае минимизации целевой функции прямую c1 x1 + c2 x2 = F (X ) надо двигать в направлении, противоположном вектору-градиенту.

Анализ оптимального решения на чувствительность

Рассмотрим графическое решение ЗЛП, представленное на рисунке 10. Ограничения задачи классифицируются на связывающие и не связывающие. Связывающие ограничения проходят через оптимальную точку, например, (1) и (2) на рисунке 10. Не связывающие ограничения не проходят через оптимальную точку, например, (3), (4) и (5) на рисунке 10. Аналогично, ресурс, представляемый связывающим ограничением,

16

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

(3)

(5)

(2)

(4)

В

АС

F

D (1)

О

Е

F (X ) = const

Рисунок 10 - Графическое решение задачи линейного программирования

Область допустимых решений задачи на рисунке 10 – многоугольник ОABCDE. Если прямую (2) передвигать до точки F, то это приведет к расширению области допустимых решений до многоугольника ОABCFE и к получению нового оптимального решения в точке F (рисунок 11).

(5)

(3)

(2)

(4)

В

АС

F

(1)

О

E

Fmax (X )

Рисунок 11 - Анализ максимального изменения запаса дефицитного ресурса

(2) с целью улучшения оптимального решения

17

При этом ограничение (2) станет избыточным. Новое решение (F) лучше прежнего (C), поскольку для пересечения с точкой F линия уровня F (x) = const должна пройти по направлению вектора, выходящего из начала координат и показывающего направление максимизации целевой функции, дальше точки С (рисунок 11).

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

Графический анализ максимально возможного изменения запаса недефицитного ресурса показан на рисунке 12. Передвинем не связывающее ограничение (3) до пересечения с оптимальным решением в точке С. Это соответствует уменьшению запаса недефицитного ресурса (3), который в оптимальной точке С исходной задачи расходовался не полностью. Областью допустимых решений станет многоугольник OGCDE. Оптимальное решение останется прежним (точка С). Таким образом, чтобы графически определить максимальное изменение запаса недефицитного ресурса, не меняющее оптимального решения, необходимо передвигать соответствующую прямую до пересечения с оптимальной точкой.

(5)

 

(4)

(2)

(3)

 

С

 

F

G

(1)

 

D

О

E

F (X ) = const

Рисунок 12 - Анализ максимального изменения запаса недефицитного ресурса (3), не изменяющего оптимальное решение

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

Fmax (X )
Fmax (X ) ,

 

18

 

yi

=

max приращение оптимального значения F (X )

.

 

 

 

max допустимый прирост объема i - го ресурса

Сначала наращивается запас ресурса, имеющего максимальное значение yi , затем второе по величине и т.д.

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

увеличении коэффициента целевой функции c1 или уменьшении коэффициента c2 прямая, соответствующая максимальному значению целевой функции на графике, вращается вокруг оптимальной точки по часовой стрелке. Если c1 уменьшается или увеличивается c2 , то эта прямая вращается вокруг оптимальной точки против часовой стрелки (рисунок 13).

Зафиксируем значение c2 . Оптимальное решение в точке С не будет меняться при увеличении c1 до тех пор, пока прямая соответствующая максимальному значению целевой функции, не совпадет с прямой (2). Аналогично, оптимальное решение в точке С не будет меняться при уменьшении c1 до тех пор, пока прямая не совпадет с прямой

(1). При таких поворотах точка С будет оставаться оптимальной до тех пор,

пока прямая Fmax (X ) не выйдет за границы, определяемые прямыми (1) и

(2). Если прямая Fmax (X )

выйдет за границы (1) или (2), то оптимальной

станет соответственно точка H или J.

Уменьшение цены 1 ( c1 )

(3)

или увеличение

 

цены 2 ( c2 )

 

(2)

(4)

Н

С

(5)

 

c

(1)

 

R

 

J

Fmax (X )

Увеличение

цены 1 ( c1 )

 

или уменьшение

 

цены 2 ( c2 )

Рисунок 13 - Анализ изменения коэффициентов целевой функции

19

Таким образом, нижний и верхний пределы изменения цены 1

определяются значениями коэффициента c1 , при которых прямая Fmax (X ) совпадает соответственно с прямыми (1) и (2).

Аналогичные выводы можно сделать при фиксированном значении коэффициента c1 .

Порядок выполнения работы

1.Рассчитать исходные данные задачи (таблица 2). Они определяются цифрами номера зачетки. Обозначим их буквами А, В, С, D, Е, F, где А- первая цифра зачетки, В вторая, С третья и т. д. Если в номере зачетки есть нули, то соответствующие значения в исходных данных следует принять равными 10.

2.Составить математическую модель задачи и решить ее графически.

3.Проанализировать полученное решение.

4.Оформить отчет по лабораторной работе, который должен содержать:

название работы;

исходные данные варианта;

результаты и анализ решения задачи.

Задание

В результате анализа хозяйственной деятельности предприятия выявлены неиспользованные ресурсы, позволяющие освоить выпуск двух новых видов продукции: А и В. В число этих ресурсов входит сырье двух видов в количестве С1 и С2, соответственно, имеются необходимые для этого фрезерные и строгальные станки, причем, ресурс их свободного машинного времени, составляет Мф и Мс, соответственно. Кроме того, имеются свободные трудовые ресурсы в количестве Т.

На изготовление единицы продукции вида А расходуется С1а сырья первого вида и С2а сырья второго вида, на единицу продукции вида В - С1в и С2в, соответственно. Нормы расхода машинного времени фрезерных станков на производство единицы продукции вида А - Мфа, на производство единицы продукции вида В - Мфв. Для строгальных станков эти нормы составляют соответственно Мса и Мсв. Нормы расхода трудовых ресурсов: Та и Тв .

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

Таблица 2 - Исходные данные задачи

Нормы расхода ресурсов

 

Прибыль от реализации

на единицу продукции

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

единицы продукции

вида А

вида В

 

вида А

вида В

С1а = F

C1в = А

С1=(А+В)·(В+С)·(С+Д)

 

 

С2а = Е

С2в = В

С2=(В+С)·(С+Д)·(Д+Е)

Па=В+С+Д

Пв=Д+Е+F

Мфа = Д

Мфв = С

Мф=(С+Д)·(Д+Е)·(Е+F)

Мса = С

Мсв = Д

Мс=(Д+Е)·(Е+F)·(F+А)

 

 

Та = В

Тв = Е

Та=(Е+F)·(F+А)·(А+В)

 

 

20

Контрольные вопросы

1.Как построить область допустимых решений ЗЛП?

2.Как найти оптимальное решение задачи?

3.Как можно графически определить максимальное изменение запаса дефицитного ресурса, улучшающее оптимальное решение?

4.Как графически определить максимальное изменение запаса недефицитного ресурса, не меняющее оптимального решения?

5.Как определить, запас какого из дефицитных ресурсов выгоднее увеличивать в первую очередь?

6.Объяснить построение математической модели производственной задачи.

7.Какие ограничения задачи называются связывающими, какие - не связывающими?

Лабораторная работа №3 Анализ чувствительности оптимального решения задачи оптимизации плана производства в табличном процессоре Microsoft Excel

Цель работы: приобретение навыков анализа чувствительности оптимального решения задачи линейного программирования, полученного в табличном процессоре MS Excel.

Проведем анализ чувствительности на примере решения производственной задачи.

Постановка задачи. Предприятие выпускает продукцию двух видов: А и В. На производство единицы продукции вида А расходуется 5 кг материала первого вида, 3 кг второго вида и 2 кг третьего вида. На производство единицы продукции вида В расходуется 2 кг материала первого вида и по 3 кг второго и третьего вида. Запасы этих материалов составляют соответственно 505, 393 и 348 кг. От реализации единицы готовой продукции вида А предприятие получит прибыль в размере 7 тыс. руб., а от реализации единицы продукции вида В – 4 тыс. руб.

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

В качестве элементов решения задачи примем x1 - количество единиц

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

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

F ( X ) = 7x1 + 4x2 → max ,