47285 (571981), страница 3
Текст из файла (страница 3)
Тогда сумма всех перевозок:
L=18400
Ответ:
| B1 | B2 | B3 | B4 | B5 | ai | |
| A1 | 25 | 21 200 | 20 | 50 | 18 | 200 |
| A2 | 15 100 | 30 100 | 32 | 25 | 40 | 200 |
| A3 | 23 100 | 40 | 10 200 | 12 100 | 21 200 | 600 |
| bj | 200 | 300 | 200 | 100 | 200 | 1000 |
Задача 4 (№53)
Условие:
Определить экстремум целевой функции вида
= 1112+2222+1212+11+22
при условиях:
111+1221
211+2222.
-
Найти стационарную точку целевой функции и исследовать ее (функцию) на выпуклость (вогнутость) в окрестностях стационарной точки.
-
Составить функцию Лагранжа.
-
Получить систему неравенств в соответствии с теоремой Куна-Таккера.
-
Используя метод искусственных переменных составить симплекс-таблицу и найти решение полученной задачи линейного программирования.
-
Дать ответ с учетом условий дополняющей нежесткости.
| № | b1 | b2 | c11 | c12 | c22 | extr | a11 | a12 | a21 | a22 | p1 | p2 | Знаки огр. 1 2 | |
| 53 | 6 | 1,5 | -2 | -4 | –1 | max | 2,5 | -1 | 3 | 2,5 | 7 | 13 | ||
Решение:
Целевая функция:
F= -212-22-412+61+1,52→max
Ограничения g1(x) и g2(x): 2,51-27 2,51-2–70
31+2,5213 31+2,52-130
1) определим относительный максимум функции, для этого определим стационарную точку (х10, х20):
→
2) Исследуем стационарную точку на максимум, для чего определяем выпуклость или вогнутость функции
F11 (х10, х20) = -4 < 0
F12 (х10, х20)=-4
F21 (х10, х20)=-4
F22 (х10, х20)=-2
F11 F12 -4 -4
F21 F22 -4 -2
Т.к. условие выполняется, то целевая функция является строго выпуклой в окрестности стационарной точки
3) Составляем функцию Лагранжа:
L(x,u)=F(x)+u1g1(x)+u2g2(x)=-212-22-412+61+1,52+u1 (2,51-2–7)+ u2 (31+2,52-13).
Получим уравнения седловой точки, применяя теорему Куна-Таккера:
i=1;2
Объединим неравенства в систему А, а равенства в систему В:
Система А:
Система В:
Перепишем систему А:
6-41-42+2,5u1+3u2 <0
1,5-41-22-u1+2,5u2 <0
2,51-2–70
31+2,52–130
4)Введем новые переменные
V={v1,v2}≥0; W={w1,w2}≥0
в систему А для того, чтобы неравенства превратить в равенства:
6-41-42+2,5u1+3u2 + v1=0
1,5-41-22-u1+2,5u2 + v2=0
2,51-2–7- w1=0
31+2,52–13- w2=0
Тогда
- v1=6-41-42+2,5u1+3u2
- v2=1,5-41-22-u1+2,5u2
w1=2,51-2–7
w2=31+2,52–13
Следовательно, система В примет вид:
- это условия дополняющей нежесткости.
5) Решим систему А с помощью метода искусственных переменных.
Введем переменные Y={y1; y2} в 1 и 2 уравнения системы
6-41-42+2,5u1+3u2 + v1 -y1=0
1,5-41-22-u1+2,5u2 + v2 -y2=0
2,51-2–7- w1=0
31+2,52–13- w2=0
и создадим псевдоцелевую функцию Y=My1+My2→min
Y’=-Y= -My1-My2→max.
В качестве свободных выберем х1, х2, v1, v2, u1, u2;
а в качестве базисных y1, y2, w1, w2.
Приведем систему и целевую функцию к стандартному виду, для построения симплекс-таблицы:
y1=6-(41+42-2,5u1-3u2 - v1)
y2=1,5-(41+22+u1-2,5u2 -v2)
w1=-7-(-2,51+2)
w2=-13-(-31-2,52)
Y’=-Y=-My1-My2=-7,5M-(-81-62+1,5u1+5,5u2+ v1+v2) M
Решим с помощью симплекс-таблицы. Найдем опорное решение:
|
|
|
|
|
|
|
| |
|
| -7,5M 4,5M | -8M 12M | -6M 3M | 1,5M 3M | 5,5M -7,5M | M 0 | M -3M |
|
| 6 -3 | 4 -8 | 4 -2 | -2,5 -2 | -3 5 | -1 0 | 0 2 |
|
| 1,5 3/4 | 4 2 | 2 0,5 | 1 0,5 | -2,5 -5/4 | 0 0 | -1 -0,5 |
|
| -7 -3/4 | -2,5 -2 | 1 -0,5 | 0 -0,5 | 0 5/4 | 0 0 | 0 0,5 |
|
| -13 15/8 | -3 5 | -2,5 5/4 | 0 5/4 | 0 -25/16 | 0 0 | 0 -5/4 |
Меняем
и
|
|
|
|
|
|
|
| |
|
| -3M 3M | 4M -4M | 3M -2M | 4,5M -4,5M | -2M M | M -M | -2M 2M |
|
| 3 3/2 | -4 -2 | -2 -1 | -4,5 -9/4 | 2 0,5 | -1 -0,5 | 2 1 |
|
| 3/4 15/8 | 2 -2,5 | 0,5 -5/4 | 0,5 -45/16 | -5/4 5/8 | 0 -5/8 | -0,5 5/4 |
|
| -31/4 -15/8 | -4,5 2,5 | -0,5 5/4 | -0,5 45/16 | 5/4 -5/8 | 0 5/8 | 0,5 -5/4 |
|
| -89/8 75/32 | 2 -25/8 | 5/4 -25/16 | 5/4 -225/64 | -25/16 25/32 | 0 -25/32 | -5/4 25/16 |
Меняем
и
|
|
|
|
|
|
|
| |
|
| 0 0 | 0 0 | M 0 | 0 0 | M 0 | 0 0 | 0 0 |
|
| 3/2 77/8 | -2 -1 | -1 -3/4 | -9/4 -37/16 | 0,5 5/8 | -0,5 -5/8 | 1 3/4 |
|
| 21/8 77/32 | -0,5 -1/4 | -3/4 -3/16 | -37/16 -37/64 | 5/8 5/32 | -5/8 -5/32 | 3/4 -3/16 |
|
| -77/8 77/16 | -2 -0,5 | 3/4 -3/8 | 37/16 -37/32 | -5/8 5/16 | 5/8 -5/16 | -3/4 3/8 |
|
| -281/32 693/128 | -9/8 -9/16 | -5/16 -27/64 | -145/64 -333/256 | 25/32 45/128 | -25/32 -45/128 | 5/16 27/64 |
Меняем
и
|
|
|
|
|
|
|
| |
|
| 0 0 | 0 0 | M 0 | 0 0 | M 0 | 0 0 | 0 0 |
|
| 89/8 431/18 | -1 -16/9 | -7/4 | -73/16 | 9/8 | -9/8 | 7/4 |
|
| 161/32 431/72 | -1/4 -4/9 | -15/16 | -185/64 | 25/32 | -25/32 | 9/16 |
|
| 77/16 431/36 | -0,5 -8/9 | -3/8 | -37/32 | 5/16 | -5/16 | 3/8 |
|
| -431/32 431/18 | -9/16 -16/9 | -47/64 | -913/256 | 145/128 | -145/128 | 47/64 |
Меняем
и
|
|
|
|
|
|
|
| |
|
| 0 | 0 | M | 0 | M | 0 | 0 |
|
| 2525/72 | ||||||
|
| 3173/288 | ||||||
|
| 2417/144 | ||||||
|
| 431/18 |
Итак,
=
=
=
=
=
,
=16,785,
=11,017,
=23,944,
=35,07
6) Условия дополняющей нежесткости выполняются
,значит, решения исходной задачи квадратичного программирования существует.
Ответ: существует.
Литература.
1) Курс лекций Плотникова Н.В.
2) Пантелеев А.В., Летова Т.А. «Методы оптимизации в примерах и задачах».
0>0>













