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

Курсовая работа

по дисциплине

Исследование операций

Руководитель:

Плотникова Н. В.             

«____» ___________ 2005 г.

Автор:

Студент группы ПС-346

Попов А. Е..                      

«____» ___________ 2005 г.

Работа защищена

с оценкой                          

«____» ___________ 2005 г.


Оглавление

1 Условия задач. 3

2 Решение задач исследования операций. 4

2.1 Решение задачи 1. 4

2.2 Решение задачи 2. 8

2.3 Решение задачи 3. 12

2.4 Решение задачи 4. 17


1 Условия задач
2 Решение задач исследования операций

 

2.1 Решение задачи 1

Для составления математической модели задачи введём переменные:

 – количество горючего, доставляемое со склада A на бензоколонку 1

– количество горючего, доставляемое со склада A на бензоколонку 2

x3a – количество горючего, доставляемое со склада A на бензоколонку 3

x1b – количество горючего, доставляемое со склада B на бензоколонку 1

x2b – количество горючего, доставляемое со склада B на бензоколонку 2

x3b – количество горючего, доставляемое со склада B на бензоколонку 3

x1c – количество горючего, доставляемое со склада C на бензоколонку 1

x2c – количество горючего, доставляемое со склада C на бензоколонку 2

x3c – количество горючего, доставляемое со склада C на бензоколонку 3

На складах A, B, C  находится 90, 60, 90 тонн горючего соответственно, следовательно, можно записать:

На каждую заправку нужно оправить одинаковое количество горючего, равное (90+60+90)/3:

 

В соответствии со стоимостями перевозок запишем целевую функцию, которую необходимо минимизировать:

Имеем классическую транспортную задачу с числом базисных переменных, равным n+m–1 , где m–число пунктов отправления, а n – пунктов назначения. В решаемой задаче число базисных переменных равно 3+3-1=5.

Число свободных переменных соответственно 9-4=4.

Примем переменные x1a, x1b, x2a, x2с, x3с в качестве базисных, а переменные x1c, x2b, x3а,  x3b  в качестве свободных (данный выбор позволяет легко выразить базисные переменные через свободные).

Далее в соответствии с алгоритмом Симплекс метода необходимо выразить базисные переменные через свободные:

Следующий шаг решения – представление целевой функции через свободные переменные:

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

Составим Симплекс таблицу:


bi x3a x2b x3b x1c
L

630

-10

-3

   1       

-1

                0

-4         

                4

1

                -1

x1a

20

       -10

0

                1

-1

                0

-1

                1

1

          -1

x1b

60

                0

0

         0

1

                0

1

                0

0

                0

x2a

70

                10

1

                -1

1

                0

1

         -1

-1

                1

x2c

10

                10

-1

          -1

0

                0

-1

                -1

1

                1

x3c

80

                0

1

                0

0

                0

1

                0

0

                0

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

bi x3a x2b x3b x2c
L 620 -2 -1 0 -1
x1a 10 1 -1 0 -1
x1b 60 0 1 1 0
x2a 80 0 1 0 1
x1c 10 -1 0 -1 1
x3c 80 1 0 1 0

Все коэффициенты при свободных переменных неположительные, следовательно, найденное решение является оптимальным. Запишем его:

x1a=10;     x1b=60;    x1c=10;

x2a=80;     x2b=0;      x2c=0;

x3a=0;     x3b=0;      x3c=80;

L=620;

Для проверки правильности вычислений можно составить транспортную таблицу:

A B C
1 10 60 10 80
2 80 0 0 80
3 0 0 80 80
90 60 90

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

Ответ:

x1a=10     x1b=60    x1c=10

x2a=80     x2b=0      x2c=0

x3a=0       x3b=0      x3c=80

L=620

2.2 Решение задачи 2

Составим систему ограничений исходя из условия задачи

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

Пусть переменные x1 и x2  - свободные, а переменные x3, x4 и x5 – базисные.

Далее необходимо представить систему ограничений в стандартном виде. Для этого проведем ряд преобразований:

 

Подставим выражения для x3 и x4 в третье уравнение системы ограничений:

Упростим полученное выражение и выразим x5:

Теперь можно представить систему ограничений в стандартном виде:

Необходимо также выразить целевую функцию через свободные переменные:

Теперь можно заполнить Симплекс-таблицу

bi x1 x2
L 1 -1 -3
x3 2 -1 2
x4 2 1 1
x5 1 1 -1

Исходя из того, что все свободные члены положительны, можно сделать вывод о том принятое решение является опорным.

Далее нужно выбрать разрешающий элемент. В качестве разрешающего столбца целесообразно принять столбец x1, так как коэффициент при x1 в целевой функции меньше коэффициента при x2. Разрешающей строкой будет строка x5­, так как отношение свободного члена этой строки к коэффициенту при x1 минимально. Отметим найденный разрешающий элемент в таблице, а также заполним необходимые клетки:

bi x1 x2
L

1              

                   1

-1

               1    

-3

                 -1

x3

2

                    1

-1

              1

2

                 -1

x4

2

                -1

1                         

            -1

1

           1

x5

1

                1

1

                 1                                                                                                           

-1

              -1

Перерисуем таблицу с учётом замены x2 на x3:

bi x5 x2
L 2 1 -4
x3 3 1 1
x4 1 -1 2
x1 1 1 -1

Коэффициент при х2 в целевой функции отрицателен, значит необходимо произвести ещё одну замену. В качестве разрешающей строки примем x3.  Таким образом, разрешающим будет элемент, стоящий на пересечении строки x3 и столбца x2.

bi x5 x2
L

2

               12

1

           4

-4

          4

x3

3

                  3

1

           1

1

         1

x4

1                       

               -6

-1

                 -2

2

         -2

x1

1

               3

1

            1

-1

         1

В итоге получим:

bi x5 x3
L 14 5 4
x2 3 1 1
x4 -5 -1 0
x1 4 2 1

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

Ответ:

x1=4

x2=3

x3=0

x4=-5

x5=0

L=14

 

2.3 Решение задачи 3

Условие задачи задано в виде транспортной таблицы:

   ПН

ПО

B1 B2 B3 запасы
A1 50 15 10 300
A2 21 30 20 100
A3 18 40 25 200
A4 23 22 12 800
A5 25 32 45 200
заявки 500 300 800

Применим к задаче метод «Северо-Западного угла». Для  этого заполним таблицу начиная с левого верхнего угла без учёта стоимости перевозок:

   ПН

ПО

B1 B2 B3 запасы
A1 300 300
A2 100 100
A3 100 100 200
A4 200 600 800
A5 200 200
заявки 500 300 800

В таблице заполнено n+m-1=7 клеток, значит найденное решение является опорным. Далее необходимо улучшить план перевозок в соответствии со стоимостями доставки грузов. Для этого используем циклические перестановки в тех циклах, где цена отрицательна.

   ПН

ПО

B1 B2 B3 запасы
  A1

    50

300

    15     10 300
A2

     21

100

    30      20 100
A3

     18

100

    40

100

     25

200
A4      23

    22

200

     12

600

800
A5     25     32

     45

200

200
заявки 500 300 800

В данной таблице в верхней части ячейки указана стоимость перевозки, а в нижней количество перевозимого груза. Прямоугольником выделен отрицательный цикл  γ1=25+22-40-12=-5. Минимальное значение перевозок, стоящих в отрицательных вершинах равно k1=100. В итоге получим уменьшение стоимости перевозки: 

ΔL1=-5*100=-500

Транспортная таблица примет следующий вид:

   ПН

ПО

B1 B2 B3 запасы
  A1

    50

300

    15     10 300
A2

     21

100

    30      20 100
A3

     18

100

    40

     25

100

200
A4      23

    22

300

     12

500

800
A5     25

    32

     45

200

200
заявки 500 300 800

γ2=12+32-45-22=-23            k2=200            ΔL2=-23*200=-4600

   ПН

ПО

B1 B2 B3 запасы
  A1

    50

300

    15

    10

300
A2

     21

100

    30      20 100
A3

     18

100

    40

     25

100

200
A4      23

    22

100

     12

700

800
A5     25

    32

200

     45 200
заявки 500 300 800

γ3=10+18-50-25=-47            k3=100            ΔL3=-47*100=-4700

   ПН

ПО

B1 B2 B3 запасы
  A1

    50

200

    15

    10

100

300
A2

     21

100

    30      20 100
A3

     18

200

    40      25 200
A4

     23

    22

100

     12

700

800
A5     25

    32

200

     45 200
заявки 500 300 800

γ4=10+23-12-50=-29            k4=200            ΔL4=-29*200=-6800

   ПН

ПО

B1 B2 B3 запасы
  A1     50     15

    10

300

300
A2

     21

100

    30      20 100
A3

     18

200

    40      25 200
A4

     23

200

    22

100

     12

500

800
A5     25

    32

200

     45 200
заявки 500 300 800

Отрицательных циклов в транспортной таблице больше нет. Следовательно, можно предположить, что найденное решение является оптимальным. Для проверки применим метод потенциалов.

Составим систему:

Положим β2=0, тогда α4=-22

β1=1,        α2=-20

β3=-10,     α2=-22

α1=-20,     α5=-32

Все коэффициенты α отрицательны, значит, найденный план перевозок является оптимальным.

Ответ:

x21=100;

x31=200;

x41=200;

x42=100;

x52=200;

x13=300;

x43=500.

 

2.4 Решение задачи 4

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

Найти минимум функции f(x1,x2)

 

При ограничениях

Заменив знак функции f(x1,x2) на противоположный, перейдем к поиску максимума функции:

Теперь задача приведена к стандартному виду задачи квадратичного программирования. Приступим к решению.

1) Определим стационарную точку

Решив систему, получим:

x1=10

x2=7

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

2) Составим функцию Лагранжа:

Применив к функции Лагранжа теорему Куна-Таккера, будем иметь систему:

3) Преобразуем полученную систему:

 

Из уравнения 3 системы следует, что x2=6-x1:

Для обращения неравенств системы в равенства введём V1, V2, W и преобразуем систему:

4) Запишем условия дополняющей нежесткости:

5) Введем в систему уравнений искусственные переменные z1,z2:

Поставим задачу максимизации функции  .

Для решения этой задачи воспользуемся Симплекс-методом. Примем переменные z1 и  z2 в качестве базисных:

Составим Симплекс таблицу:

bi x1 U1 U2 V1 V2
φ

-17M

        0

-5M

                0

0

                0

   M             

        0

M

        0 

-M

                0

z1

9

      8

2

      3

-1

       1

 2                

       -3

-1

       0

0

      1

z2

8

        8  

3

       3               

1

      1

-3          

      -3

0

        0

1

        1

W

0

       0

-1

       0

0

        0

0

       0

0

        0

0

        0

bi x1 z2 U2 V1 V2
φ

-17M

    17M

-5M

        M

0

       M

   M             

      -M

M

      -M

-M

        M

z1

17

    17/5

5

      1/5

1

      1/5

 -1                

     -1/5  

-1

     -1/5

1

      1/5     

U1

8

   -51/5   

3

     -3/5

1

     -3/5

-3          

     3/5    

0

      3/5     

1

     -3/5     

W

0

    17/5

-1

      1/5

0

      1/5

0

     -1/5   

0

    -1/5

0

     1/5

bi z1 z2 U2 V1 V2
φ 0 M M 0 0 0
x1 17/5 1/5 1/5 -1/5 -1/5 1/5
U1 -11/5 -3/5 -2/5 1/2 3/5 -2/5
W 17/5 1/5 1/5 -1/5 -1/5 1/5

В итоге получим

x1=17/5

x2=6-x1=13/5

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

Условия дополняющей нежесткости

 выполняются.

Следовательно, найденное решение является оптимальным.

Найдем значения целевой функции:

=- 51/5 - 52/5 + 289/50 – 221/25 + 169/25 =

= -16.9

Ответ:

x1 = 17/5

x2 = 13/5

f(x1,x2) = -16.9


Еще из раздела Информатика, программирование:


 Это интересно
 Реклама
 Поиск рефератов
 
 Афоризм
Странные наблюдения: ошибочно набранный номер телефона никогда не бывает занят…
 Гороскоп
Гороскопы
 Счётчики
bigmir)net TOP 100