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

Введение в математическое программирование

..pdf
Скачиваний:
14
Добавлен:
15.11.2022
Размер:
1.81 Mб
Скачать

Алгоритм метода множителей Лагранжа для задачи (4.2.1)–(4.2.2)

1. Составляем функцию Лагранжа

m

F(x1, x2, ..., λ1, λ2 , ..., λm) = f (x1, x2 , ..., xn ) + λi (bi −ϕi (x1, ..., xn )).

i=1

2. Находим частные производные функции Лагранжа по x j

и λi ( j =

1, n

),

(i =

1, m

) и приравниваем их к нулю:

 

 

 

F

 

f

m

∂ϕ

 

 

 

 

 

 

 

 

 

i = 0, ( j =1, n),

 

 

 

 

 

=

 

 

λi

 

 

 

x j

x j

 

 

 

 

 

i=1

x j

(4.2.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(i =1, m).

 

 

ϕi (x1, x2 , ..., xn ) = bi

 

3.Решаем систему (4.2.3) и определяем точки, в которых функция f (x1, x2 , ..., xn ) может иметь экстремум.

4.Проверяем полученные точки на экстремум (например,

сприменением достаточных условий локального экстремума)

и определяем экстремальное значение функции f (x1, x2 , ..., xn )

в найденной точке.

Пример 13

Фирма реализует автомобили двумя способами: через магазин и через торговых агентов. При реализации x1 автомобилей

через магазин расходы на реализацию составляют 4x1 + x12 у.е., а при продаже x2 автомобилей через торговых агентов расходы

составляют x22 у.е. Найти оптимальный способ реализации ав-

томобилей, минимизирующий суммарные расходы, если общее число предназначенных для продажи автомобилей составляет

200 штук.

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

Пусть x j – число автомобилей, реализуемых первым и вто-

рым способом( j =1,2) .

91

Целевая функция Z = 4x1 + x12 + x22 (min)

(суммарные издержки минимизируются). Ограничения задачи:

x1 + x2 = 200,

x1 0,x2 0.

Итак, математическая модель задачи имеет вид

Z = 4x1 + x12 + x22 (min) ;

x1 + x2 = 200,

x1 0,x2 0.

Для решения задачи применим метод множителей Лагранжа. Функция Лагранжа имеет вид

F(x1, x2 , λ)= 4x1 + x12 + õ22 + λ(200 – x1 x2 ).

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

F

 

= 4 + 2x λ = 0,

 

x

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

F

 

= 2x2

λ = 0,

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

F

= 200

x

x

 

= 0.

 

 

2

λ

 

 

1

 

 

 

 

 

 

 

 

Решая систему, найдем

x1 =99, x2 =101, λ = 202; F(x1, x2 ) = 20 398.

Составим определитель из частных производных 2-го порядка для функции F(x1, x2 ) по x1 и x2 :

92

2F = 2 > 0 , то в полученной точке функция
x12

2F

2F

 

 

 

 

 

 

x2

x1x2

 

 

2

0

 

 

 

 

1

 

 

=

 

 

 

 

= 4 > 0.

2F 2F

 

0

2

 

x2x1

 

x2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

Следовательно, в силу достаточного условия существования условного экстремума, функция F(x1, x2 ) в точке x1 =99 ,

x2 =101 имеет экстремум.

Так как

F(x1, x2 ) имеет условный минимум.

Итак, оптимальное решение: X * = (99,101), экстремальное значение целевой функции:

Z(X *) = Zmin = 20 398 у.е.

Таким образом, нужно продать 99 автомобилей через магазин и 101 автомобиль через торговых агентов, при этом расходы минимизируются и составляют 20 398 у.е.

ВАРИАНТЫ ЗАДАНИЙ ДЛЯ КОНТРОЛЬНЫХ РАБОТ

1. В табл. 4.2.1 приведены данные о предприятии, производящем продукцию двух видов Р1 и Р2. Для изготовления этой продукции использованы три вида ресурсов S1, S2, S3. Запасы ресурсов равны соответственно b1, b2, b3. Расход i-го ресурса Si на изготовление единицы продукции вида Рj равен

aij (i =1,3; j =1,2) .

Доход, получаемый предприятием от реализации единицы j-го вида продукции, равен cj ( j =1,2) .

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

Решитьзадачугеометрическимспособомисимплекс-методом.

93

Таблица 4.2.1

 

 

Виды

Запасы

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

продукции

(b )

 

P1

 

P2

i

S1

à11

 

à12

b1

S2

à21

 

à22

b2

S3

à31

 

à32

b3

Доход (ñj )

ñ1

 

ñ2

 

Варианты заданий

1.

 

P1

P2

bi

 

S1

3

5

55

 

S2

1

5

45

 

S3

5

2

60

 

cj

5

6

 

3.

 

 

 

 

 

P1

P2

bi

 

S1

3

5

65

 

S2

5

2

55

 

S3

5

4

65

 

cj

6

7

 

5.

 

 

 

 

 

P1

P2

bi

 

S1

5

6

68

 

S2

3

2

36

 

S3

1

4

36

 

cj

5

8

 

7.

 

 

 

 

 

P1

P2

bi

 

S1

1

4

48

 

S2

5

3

75

 

S3

5

4

80

 

cj

4

7

 

2.

 

P1

P2

bi

 

S1

4

1

56

 

S2

2

7

77

 

S3

5

6

89

 

cj

3

5

 

4.

 

 

 

 

 

P1

P2

bi

 

S1

1

5

55

 

S2

4

5

70

 

S3

3

2

42

 

cj

6

5

 

6.

 

 

 

 

 

P1

P2

bi

 

S1

5

3

49

 

S2

3

5

55

 

S3

3

1

27

 

cj

7

6

 

8.

 

 

 

 

 

P1

P2

bi

 

S1

4

5

59

 

S2

1

6

48

 

S3

3

2

39

 

cj

3

7

 

94

9.

 

P1

P2

bi

 

S1

2

5

65

 

S2

5

2

75

 

S3

4

3

67

 

cj

6

7

 

11.

 

 

 

 

 

P1

P2

bi

 

S1

1

7

70

 

S2

3

2

54

 

S3

2

3

41

 

cj

3

7

 

13.

 

 

 

 

 

P1

P2

bi

 

S1

2

7

98

 

S2

3

2

62

 

S3

4

1

76

 

cj

2

5

 

15.

 

 

 

 

 

P1

P2

bi

 

S1

1

3

45

 

S2

3

5

83

 

S3

5

2

75

 

cj

7

5

 

17.

 

 

 

 

 

P1

P2

bi

 

S1

5

2

65

 

S2

5

4

75

 

S3

5

7

105

 

cj

7

4

 

19.

 

 

 

 

 

P1

P2

bi

 

S1

1

4

32

 

S2

2

3

29

 

S3

5

2

45

 

cj

8

7

 

10.

 

P1

P2

bi

 

S1

5

2

60

 

S2

5

6

80

 

S3

1

4

44

 

cj

7

5

 

12.

 

 

 

 

 

P1

P2

bi

 

S1

2

3

41

 

S2

2

7

77

 

S3

5

2

75

 

cj

7

6

 

14.

 

 

 

 

 

P1

P2

bi

 

S1

1

5

60

 

S2

3

7

92

 

S3

4

3

72

 

cj

3

5

 

16.

 

 

 

 

 

P1

P2

bi

 

S1

1

2

32

 

S2

4

1

52

 

S3

3

2

44

 

cj

5

6

 

18.

 

 

 

 

 

P1

P2

bi

 

S1

2

5

95

 

S2

5

2

80

 

S3

5

1

75

 

cj

4

5

 

20.

 

 

 

 

 

P1

P2

bi

 

S1

2

5

75

 

S2

3

1

39

 

S3

7

6

113

 

cj

5

8

 

95

21.

 

P1

P2

bi

 

S1

4

3

49

 

S2

3

7

70

 

S3

3

1

33

 

cj

3

5

 

23.

 

 

 

 

 

P1

P2

bi

 

S1

5

2

65

 

S2

6

5

91

 

S3

1

3

39

 

cj

3

4

 

25.

 

 

 

 

 

P1

P2

bi

 

S1

7

2

91

 

S2

3

5

68

 

S3

1

6

66

 

cj

4

3

 

27.

 

 

 

 

 

P1

P2

bi

 

S1

5

4

66

 

S2

1

3

33

 

S3

2

1

24

 

cj

3

4

 

29.

 

 

 

 

 

P1

P2

bi

 

S1

1

5

50

 

S2

2

5

55

 

S3

5

2

85

 

cj

3

4

 

22.

 

P1

P2

bi

 

S1

5

3

75

 

S2

4

7

83

 

S3

1

5

50

 

cj

4

5

 

24.

 

 

 

 

 

P1

P2

bi

 

S1

5

2

65

 

S2

2

5

60

 

S3

5

6

85

 

cj

4

3

 

26.

 

 

 

 

 

P1

P2

bi

 

S1

5

2

60

 

S2

4

5

65

 

S3

2

5

55

 

cj

4

3

 

28.

 

 

 

 

 

P1

P2

bi

 

S1

5

4

61

 

S2

4

1

40

 

S3

2

5

55

 

cj

3

4

 

30.

 

 

 

 

 

P1

P2

bi

 

S1

5

2

55

 

S2

1

6

60

 

S3

4

3

51

 

cj

3

5

 

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

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

96

Варианты заданий

2x1 +3x2 x3 +2x4 = 6, 1. x1 +2x2 3x3 + x4 = 7.

x1 + 2x2 x3 + 4x4 = 6, 3. 3x1 x2 + 2x3 x4 = 3.

2x1 2x2 +3x3 4x4 = 7, 5. 2x1 4x2 x3 +2x4 = −1.

2x1 3x2 + 4x3 5x4 = −2, 7. 3x1 + 4x2 + x3 2x4 = 6.

3x1 + 4x2 + x3 x4 = 7, 9. 2x1 6x2 + 2x3 +3x4 =1.

3x1 4x2 4x3 + 2x4 = −3, 11. 2x1 + x2 +3x3 x4 = 5.

4x1 5x2 3x3 + 2x4 = −2, 13. 3x1 2x2 +5x3 +3x4 = 9.

3x1 4x2 6x3 + 2x4 = −5, 15. 5x1 5x2 + 2x3 +3x4 = 5.

4x1 2x2 + 4x3 + x4 = 7, 17. 6x1 3x2 5x3 2x4 = −4.

3x1 x2 + 2x3 + 4x4 =8, 19. 6x1 4x2 5x3 2x4 = −5.

3x1 4x2 + x3 2x4 = −2, 2. 2x1 + 2x2 2x3 + x4 = 3.

x1 3x2 3x3 + x4 = −4, 4. x1 + 2x2 + 4x3 x4 = 6.

2x1 4x2 +5x3 8x4 = −5, 6. 3x1 x2 2x3 + x4 =3.

3x1 +3x2 4x3 6x4 = −4, 8. 2x1 2x2 +5x3 +3x4 =8.

4x1 3x2 +5x3 2x4 = 4, 10. x1 4x2 + 2x3 5x4 = −6.

5x1 + 2x2 x3 x4 = 5, 12. 4x1 + x2 3x3 4x4 = −2.

x1 3x2 4x3 2x4 = −8, 14. 3x1 5x2 + x3 + 4x4 = 3.

6x1 5x2 3x3 x4 = −3, 16. 2x1 x2 + 4x3 +5x4 =10.

5x1 +3x2 6x3 4x4 = −2, 18. 6x1 + 4x2 2x3 +3x4 =11.

2x1 + x2 + 2x3 + 4x4 = 9, 20. 5x1 3x2 + x3 6x4 = −3.

97

21.

2x

+ 2x

 

 

x

x

 

=1,

22.

3x

+5x

 

 

+ x

+ x

 

= 32,

 

x

1

 

 

 

 

2

3

 

 

4

 

 

 

x

1

 

 

2

 

 

3

 

 

4

 

 

 

 

 

 

x

2

 

3x x

4

= −1.

 

 

3x

2

 

x + x

4

= −8.

 

1

 

 

 

 

 

 

3

 

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

23.

3x1 +8x2 + x3 + x4 = 50,

24.

4x1

 

 

 

 

+ x3 + x4 =16,

 

5x1 4x2 x3 + x4 =14.

 

6x1 4x2 x3 + x4 = 4.

25.

x

2x

2

+ 2x

+ x

4

= −1,

26.

2x + 2x

2

2x + x

4

 

= 3,

 

1

 

 

 

 

 

 

3

x

= 4.

 

 

1

 

 

 

 

 

 

3

 

 

 

 

= 2.

 

 

2x

4x

2

+ x

4

 

 

3x + 4x

2

x + 2x

4

 

 

1

 

 

 

 

3

 

 

 

 

 

 

 

1

 

 

 

 

 

3

 

 

 

 

 

27.

x

+3x

 

+3x

x

 

= 4,

28.

3x + x

 

 

2x + x

 

= 3,

 

x

1

 

 

 

 

 

2

 

3

 

 

4

 

 

 

1

 

2

 

 

 

 

3

 

 

4

 

 

 

= −5.

 

 

+ 2x

2

+ 4x

x

4

= 6.

 

 

2x 4x

+5x 8x

4

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

1

 

 

 

2

 

 

 

3

 

 

 

 

29.

2x1 2x2 +5x3 +3x4 =8,

30.

6x1 + 4x2 2x3 +3x4 =11,

 

3x1 +3x2 4x3 6x4 = −4.

 

5x1 3x2 +6x3 + 4x4 = 2.

 

4. Решить транспортную задачу, данные для которой при-

ведены

 

в

 

табл.

4.2.2.

В ней

a1, a2 ,

a3, a4

 

запасы груза

в пунктах отправления А1, А2, А3, А4; b1, b2 , b3, b4 – потребности в грузе в пунктах назначения B1, B2, B3, B4; cij – стоимость

перевозки единицы груза из i-го пункта отправления в j-й пункт назначения.

Требуется составить план перевозок, при котором суммарные транспортные издержки минимизируются. Решить транспортную задачу методом потенциалов.

 

 

 

 

 

 

 

Таблица 4.2.2

Поставщик

 

Потребитель

 

Запас груза (аi)

В1

 

В2

В3

 

В4

А1

c11

 

с12

с13

 

с14

a1

A2

c21

 

с22

с23

 

с24

a2

A3

c31

 

с32

с33

 

с34

a2

A4

c41

 

с42

с43

 

с44

а4

Потребность в грузе (bj)

b1

 

b2

b3

 

b4

 

 

 

 

 

 

 

 

 

 

 

98

 

 

 

 

Варианты заданий

1.

 

В1

В2

В3

В4

аi

 

А1

9

3

6

5

35

 

А2

4

10

11

8

70

 

А3

2

2

3

9

65

 

А4

3

4

9

12

30

 

bj

40

85

25

50

 

3.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

7

6

9

8

130

 

А2

6

9

11

12

80

 

А3

5

5

7

9

40

 

А4

4

8

10

9

35

 

bj

90

25

105

65

 

5.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

10

9

3

4

110

 

А2

5

6

9

6

75

 

А3

3

9

2

5

60

 

А4

6

8

7

7

100

 

bj

80

125

90

50

 

7.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

4

8

8

10

90

 

А2

5

6

9

5

35

 

А3

3

4

10

11

70

 

А4

5

7

5

8

90

 

bj

80

65

55

85

 

9.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

11

6

4

3

100

 

А2

5

8

4

5

25

 

А3

10

10

11

6

55

 

А4

7

9

5

5

70

 

bj

30

65

80

75

 

2.

 

В1

В2

В3

В4

аi

 

А1

11

2

10

3

20

 

А2

6

6

9

10

90

 

А3

4

8

3

8

95

 

А4

9

11

2

5

95

 

bj

75

80

45

100

 

4.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

9

4

2

2

65

 

А2

3

8

5

7

10

 

А3

2

7

3

4

110

 

А4

4

5

9

6

25

 

bj

40

45

95

30

 

6.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

10

4

5

9

60

 

А2

10

5

3

5

80

 

А3

5

9

10

4

45

 

А4

7

6

4

5

95

 

bj

30

20

110

120

 

8.

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

10

9

3

35

 

А2

5

6

6

2

100

 

А3

4

6

8

7

30

 

А4

7

8

5

3

85

 

bj

105

20

60

65

 

10.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

7

10

6

5

85

 

А2

6

11

7

4

115

 

А3

4

5

6

8

40

 

А4

8

5

4

6

60

 

bj

90

75

120

15

 

99

11.

 

В1

В2

В3

В4

аi

 

А1

6

3

7

10

45

 

А2

10

4

12

10

100

 

А3

5

9

8

11

20

 

А4

4

2

4

8

75

 

bj

30

80

95

35

 

13.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

10

6

5

65

 

А2

4

3

5

9

80

 

А3

11

4

4

8

35

 

А4

5

5

3

6

90

 

bj

60

50

85

75

 

15.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

7

6

5

8

80

 

А2

8

8

9

7

100

 

А3

12

10

9

8

60

 

А4

8

9

6

10

70

 

bj

45

90

65

110

 

17.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

5

5

4

6

40

 

А2

3

7

9

6

45

 

А3

7

8

5

9

80

 

А4

6

4

6

8

60

 

bj

90

75

50

10

 

19.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

6

7

4

8

70

 

А2

6

4

5

3

80

 

А3

9

8

6

10

40

 

А4

7

4

3

5

60

 

bj

20

45

95

90

 

12.

 

В1

В2

В3

В4

аi

 

А1

7

3

6

9

45

 

А2

6

12

7

5

85

 

А3

9

9

5

4

20

 

А4

8

10

12

6

90

 

bj

25

60

35

120

 

14.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

4

5

4

80

 

А2

4

8

7

3

25

 

А3

9

9

6

6

60

 

А4

5

5

3

7

45

 

bj

30

35

85

60

 

16.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

8

6

4

20

 

А2

3

10

11

6

90

 

А3

4

6

9

4

65

 

А4

9

7

8

7

45

 

bj

30

70

50

70

 

18.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

6

4

9

95

 

А2

7

4

6

6

70

 

А3

12

5

5

11

60

 

А4

11

12

9

10

85

 

bj

40

100

40

130

 

20.

 

 

 

 

 

 

 

В1

В2

В3

В4

аi

 

А1

8

3

4

9

40

 

А2

8

7

5

6

80

 

А3

10

12

9

9

120

 

А4

8

11

5

10

30

 

bj

50

100

60

60

 

100