Линейное программирование. 8

Министерство  образования РФ

ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ УПРАВЛЕНИЯ

Институт  информационных систем управления

Кафедра прикладной математики 
 
 
 
 
 
 
 
 
 

Курсовой  проект

по  дисциплине «Прикладная  математика»

Вариант №17 
 
 
 
 
 
 
 
 
 
 

Выполнила:

студентка II курса 1 группы

Молодова  А.М.

Проверил  преподаватель:

Никифорова  О.П.  
 

Москва 2011

Содержание

   Задание  на курсовую работу……………………………………………………................3

1. Линейная производственная задача……………………………………….......................6

2. Двойственная задача линейного программирования…………......................................12

3. Задача о "расшивке узких мест производства"………………………………................14

4.Транспортная задача линейного программирования…………………….......................17

5. Динамическое программирование. Задача распределения капитальных вложений....20

6. Динамическая  задача управления производственными  запасами.................................23

7. Анализ доходности и риска финансовых операций...............…………….....................27

8. Матричная модель производственной программы………………………......................30

    Список  использованной литературы……………………………………………….........33 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Задание на курсовую работу

 
  1. Составить математическую модель линейной производственно  задачи, имея следующие исходные данные:

         ,

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

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

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

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

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

    Применить наёденные двойственные оценки ресурсов к решению следующей задачи.  

    3.       Сформулировать задачу о «расшивке узких мест производства» и составить математическую модель. Определить область устойчивости двойственных     оценок, где сохраняется структура программы производства. Решить задачу о «расшивке узких мест производства» при условии, что дополнительно можно получить от поставщиков не более одной трети первоначально выделенного объёма ресурса любого вида (если задача окажется с двумя переменными, то только графически); найти план приобретения дополнительных объёмов ресурсов, дополнительную возможную прибыль.

    Составить сводку результатов.

      4.     Составить математическую модель транспортной задачи по следующим      данным:

    Вектор  объёма производства

    Вектор  объёма потребления 

    Матрица транспортных издержек:

    Если  полученная модель окажется открытой, то свести её к замкнутой и найти  оптимальное решение транспортной задачи методом потенциалов.  

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

     

    6.       Рассмотреть динамическую задачу  управления производством и запасами с исходными данными:

     d1 d2 d3  a b c  h1 h2 h3      y1

           2 2 2  1 1 5  1 2 4         2 ,

    где  d1,d2,d3 - заказ на изделия на каждый месяц, a, b, с – коэффициенты в квадратичной функции, h1,h2, h3 – затраты на хранение единицы изделия для каждого месяца, y1  - к началу первого месяца на складе имеется 2 единицы изделия. 

    7.       Провести анализ доходности и  риска финансовых операций по исходным    данным:

      17. (0,1/2) (4,1/4) (8,1/8) (32, 1/2)

      18. (-6,1/2) (-4,1/4) (-2,1/8) (10, 1/8)

      19. (0,1/4) (8,1/4) (12,1/3) (24, 1/6)

      20. (-6,1/4) (-2,1/4) (0,1/3) (-6, 1/6) 
     
     
     
     
     
     
     

    8.       Составить матричную модель производственной программы предприятия   по следующим исходным данным:

              структурная матрица производства  А =

     матрица коэффициентов прямых затрат (затраты  на физический выпуск)

     B = 

      матрица коэффициентов прямых затрат внешних ресурсов  Y =

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

    

 
 
 
 
    1. Линейная  производственная задача

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

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

      Применим следующие обозначения:

  i – номер группы оборудования (i=1,2, … , m);

  j – номер вида изделия (j=1,2, … , n);

  aij – норма времени на обработку единицы i-го изделия на j-ой группе оборудования;

  bi – действительный фонд времени работы i-й группы оборудования;

  xi – планируемое количество единиц  j-го изделия;

  (x1, x2, … , xn) – искомый план производства. 

Дано:

                  1 4 3 4  120

            A= 3 0 2 2 B= 168 C= (31,10,14,20)

                  2 5 0 3  80 

   Необходимо  сформулировать линейную производственную задачу и составить ее математическую модель, (исходные данные взяты из приложения I, где технологическая матрица А затрат различных ресурсов на единицу каждой продукции, вектор объемов ресурсов В и вектор удельной прибыли С при возможном выпуске четырех видов продукции с использованием трех видов ресурсов

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

   В последней симплексной таблице  указать обращенный базис Q-1, соответствующий оптимальному набору базисных неизвестных. Проверить выполнение соотношения

H = Q-1*B

и решить графически.

Решение

 

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

найти производственную программу

          (Х1,Х2,Х3,Х4) 

максимизирующую прибыль

          Z=31Х1+10Х2+14Х3+20Х4    

при ограничениях по ресурсам 

   Х1+4Х2+3Х3+4X4  120

 3Х1              +2Х3+2X4  168

 2Х1+5Х2               +3X4  80 

Где X1 0, Х2 0,Х3 0, Х4 0.   

Получили задачу на условный экстремум. Для ее решения систему неравенств при помощи дополнительных неотрицательных неизвестных X5, Х6,Х7 заменим системой линейных алгебраических уравнений 

Z=31Х1+10Х2+14Х3+20Х4+0Х5+0Х6+0Х7     =>    max 

  Х1+4Х2+3Х3+4X4 +X5= 120

3Х1             +2Х3+2 X4          +X6  = 168

2Х1+5Х2             +3X4                    +X7  = 80 

Где X1 0, Х2 0,Х3 0, Х4 0, Х5 0, Х6 0, Х7 0. 

Дополнительные  переменные имеют смысл остатков соответствующих ресурсов.  

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

              Таблица 1

 
Б Н 31 10 14 20 0 0 0 Пояснения
      X1 X2 X3 X4 X5 X6 X7 minj(∆j<0) = -31 

  minj(bi / ai1) =

  minj(120/1;168/3;80/2)=

  80/2 

  X7 из базиса

  Х1 в базис

0 X5 120 1 4 3 4 1 0 0
0 Х6 168 3 0 2 2 0 1 0
0 Х7 80 2 5 0 3 0 0 1
  Z0-Z 0-Z -31 -10 -14 -20 0 0 0
                     
0 Х5 80 0 15/10 3 25/10 1 0 -5/10 minj(∆j<0) = -14

   minj(bi / ai3) =

  minj(80/3;48/2;40/0)=

  48/2 

  X6из базиса

  Х3 в базис

0 Х6 48 0 -75/10 2 -25/10 0 1 -15/10
31 Х1 40 1 25/10 0 15/10 0 0 5/10
  Z0-Z 1240-Z 0 670/10 -14 265/10 0 0 155/10
                     
0 Х5 8 0 1275/100 0 625/100 1 -15/10 175/100  
 
 
 
 
              ∆j  ≥0
14 Х3 24 0 -375/100 1 -125/100 0 5/10 -75/100
31 Х1 40 1 25/10 0 15/10 0 0 5/10
  Z0-Z 1576-Z 0 15 0 9 0 7 5
 
 

Симплексный  метод   решения   линейной  производственной  задачи  представлен  в  Таблице 1. 

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

Х1=40,  Х2=0,  Х3=24,  Х4=0,

а максимальная прибыль равна 

ZМАХ = 1576 

При этом  второй и третий  ресурс будут исчерпаны  полностью, а первый ресурс имеет остаток Х5= 8 единиц т. е. узкие места производства Х6=0 и Х7=0 

Проверим выполнение соотношения H=Q-1 В, где Q-1 базис в последней симплексной таблице, соответствующий набору базисных неизвестных.  
 
 
 
 

Проверка: Q-1*B = H 
 

              1    -15/10   175/100                        120             8

   Q-1 =  0       5/10   -75/100   B=    168            H=      24

         0        0       5/10               80                       40

                                     

           1   -15/10   175/100         120        1*120+(-15/10)*168+175/100*80          8                     H =     0     5/10     -75/100   *    168  =   0*120+5/10*168+(-75/100)*80        =  24

           0         0          5/10          80         0*120+0*168+5/10*80                          40 
 
 

Вывод: при решении  симплексной таблицы ошибок допущено не было.

Графическое решение 

Проверим полученный результат.

Воспользуемся тем, что в оптимальной производственной программе Х2=0, Х4=0. Предположим, что первую и вторую продукции мы не намеревались выпускать с самого начала. Рассмотрим задачу с оставшимися двумя переменными, сохранив их нумерацию. Математическая модель задачи будет выглядеть следующим образом:

Z = 31X1 + 14X3 max 

X1  + 3X3    120  (1)

3X1  + 2X3   168  (2)

2X1  + 0X3   80    (3) 

X1 0, X3 0 

1) X1  + 3X3 = 120

     X1=0    X3=40

     X3=0    X1=120 

2) 3X1  + 2X3=168

      X1=0   X3=84

      X3=0   X1=56    

3) 2X1=80

      X1 =40

Решая эту задачу графически, получим следующее: 

Точка M лежит на пересечении прямых (2) и (3) 

3X1  + 2X3 =168

2X1   = 80 

2X1 =80  => X1 = 40 

3*40+2X3   = 168 

2X3   = 168-120

2X3   = 48

X3   = 24 

Ответ Х1=40,  Х3=24,   

Вывод:  решение  в симплексной таблице верно 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

2. Двойственная задача линейного программирования 

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

Теперь представим себе что возникла новая ситуация. Некий предприниматель, занимающийся производством каких-то других видов продукции, но с использованием таких же трех видов ресурсов, какие имеются у нас, предлагает нам «уступить» по определенным ценам все имеющиеся у нас ресурсы и обещает платить y1 руб. за каждую единицу первого ресурса, у2 руб. – второго, у3 руб. – за единицу третьего ресурса. Возникает вопрос: при каких ценах y1, у2, у3 мы можем согласиться с предложение предпринимателя?

Величины y1, у2, у3  принято называть расчетными или двойственными оценками ресурсов. Они прямо зависят от условий, в которых действует наше предприятие. 

Дано:

                   1 4 3 4  120

            A=  3 0 2 2 B= 168 C= (31,10,14,20)

                   2 5 0 3   80 

A – технологическая матрица

B – вектор объема ресурсов

C – вектор удельной прибыли 

Для производства единицы продукции первого вида мы должны затратить, как видно из матрицы А,  1 единицу ресурса первого вида, 3 единицы ресурса второго вида и 2 единицы третьего (элементы первого столбца матрицы). В ценах y1, у2, у3 наши затраты составят y1 + 3у2 +2y3,   т.е. столько заплатит предприниматель за все ресурсы, идущие на производство единицы первой продукции. На рынке за единицу первой продукции мы получили бы прибыль 31 руб.  Следовательно,  мы  можем  согласиться  с предложением предпринимателя только в том случае, если он заплатит не меньше  

   y1 + 3 у2 + 2 y3  31

    4 y1 +          5 y3  10

    3 y1 + 2 у2               14

    4 y1 + 2 у2 + 3 y3  20 

Поэтому перед предпринимателем мы ставим такие условие по всем видам продукции. Учтем, что за все имеющиеся у нас ресурсы нам должны заплатить

120у1 + 168у2 +80у3 рублей. При поставленных нами условиях предприниматель будет искать такие значения величин y1, у2, у3, чтобы эта сумма была как можно меньше. Таким образом, проблема определения расчетных оценок ресурсов приводит к задаче линейного программирования: найти вектор двойственных оценок У (У1У2,У3) минимизирующий общую оценку всех ресурсов

  f =120у1 + 168у2 +80у3 -> min

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

у1 >=0    у2>=0   у3>=0 

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

Х (Х1, Х2, Х3, Х4) и у (y1, у2, у3) пары двойственных задач необходимо и достаточно выполнение условий 

     X1( y1+3 y2+2 y3 -31 )=0  y1 ( Х1+4Х2+3Х3+4 Х4- 120)=0

    X2(4 y1 + 5 y3 - 10)=0  y2 (3Х1+ 2Х3+2Х4- 168)=0

    X3(3 y1 +2 y2 - 14)=0  y3 (2 Х1+5 Х2+ 3Х4- 80)=0

    X4(4 y1 +2 y2+3 y3 - 20)=0 

Ранее было найдено, что в решении исходной задачи X1 >0, X3>0  => (X1 =40,

X2=0, Х3=24, Х4=0).

Поэтому

          y1+3 y2+2 y3 - 31=0

          3y1+2 y2 - 14=0 

Если же учесть, что первый   ресурс был избыточным и, согласно той же теореме двойственности, его двойственная оценка равна нулю: y1 =0, то приходим к системе уравнений

          3 y2+2 y3 - 31=0

         2 y2 - 14=0

откуда следует

         y2=7, у3=5 

y*= (0;7;5) 

Таким образом, получили двойственные оценки, причем общая оценка всех ресурсов равна 1576: 

f=120*0+168*7+80*5=1576 

Так как первый ресурс был избыточным, согласно теории свойств, его двойственная оценка равна 0 (y1=0) 

Экономический смысл двойственных оценок. Например, двойственная оценка второго ресурса  у2=7 показывает, что добавление одной единицы этого ресурса обеспечит прирост прибыли в 7 единиц.

3. Задача о  расшивке узких мест производства 

       При выполнении оптимальной производственной программы второй и третий ресурсы используются полностью, т.е. образуют "узкие места производства". Будем их заказывать дополнительно. Пусть T(t1,t2,t3)- вектор дополнительных объемов ресурсов. Так как мы будем использовать найденные двойственные оценки ресурсов, то должно выполняться условие Н + Q-1T>=0.

Задача состоит  в том, чтобы найти вектор Т (0, t2 , t3) максимизирующий суммарный прирост прибыли (т. к. y1=0 то t1=0)

       W = 7 t2 +5 t3   → max

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

1 -1,5 1,75
0 0,5 -0,75
0 0 0,5

       8   0                  0

       24       +                           * t2   0

        40  t3  0 
 

         0                             120

         t2 ≤    1/3     168

         t3                             80 

причем по смыслу задачи t2 0  , t3 0 

Составим систему  уравнений для графического решения

          1,5* t2 -1,75* t3 8  (1)  точки (0; -4 ); (5 ;0)

      -0,5* t2 +0,75* t3 24 (2)  точки (0;32); (-48;0)

                    -0,5* t3 40 (3) 

             t2 (4) 

                    t3 (5)

приходим к  задаче линейного программирования 

Эту задачу решим графически:  

1,5* t2 -1,75* t3 =8 

t3

1,5* t2 -1,75* =8

1,5* t2  =      =>  t2 =16 

Программа "расшивки" имеет вид

Tопт = 0;16;26

Wmax= 7 t2 +5 t3 =7*16+5*26 =388 и прирост прибыли составит 388

При новом плане  выпуска продукции

X5 н =-1,5* 16 +1,75* 26 +8=0

X3 н =0,5* 16 -0,75* 26 +24=22

X1 н =1,5* 16 +40=53

Z н =31*53 +14*22 =1964

Z н =1576+388 =1964

Графическое решение

 
 
 
 
 
 

Сводка  результатов приведена в Таблице1  

                                                                                    Таблица1