- •Методы оптимизации
- •Методы оптимизации
- •1. Основные понятия и определения.
- •1.1 Место и роль методов оптимизации при моделировании и решении прикладных задач.
- •1.2 Общая постановка прикладной задачи.
- •1.3 Классификация оптимизационных методов и моделей.
- •1.4 Основные этапы построения оптимизационных моделей.
- •2. Линейное программирование.
- •2.1 Пример постановки задачи линейного программирования.
- •2.2 Общая постановка задачи линейного программирования.
- •2.3 Методика составления экономико-математических моделей.
- •2.4 Геометрическая интерпретация решения злп.
- •Выполнить самостоятельно. Задание №1.
- •2 .6 Симплекс-метод (см) решения злп
- •2.7 Отыскание опорного и оптимального решения злп с использованием табличного алгоритма с заменой базисных переменных.
- •Выполнить самостоятельно.
- •2.8 Двойственные задачи линейного программирования (дзлп).
- •2.9 Алгоритм составления двойственных злп:
- •3. Теорема двойственности
- •3. Целочисленное программирование.
- •3.1 Постановка задачи целочисленного программирования (зцп)
- •3.2 Метод отсечения (метод Гомори).
- •3.3 Алгоритм решения злцп
- •4. Нелинейное программирование.
- •4.1Классические методы оптимизации.
- •4.2 Метод множителей Лагранжа.
- •5. Методы направленного поиска экстремумов.
- •5.1 Методы определения экстремума унимодальной функции.
- •5.2 Методы определения локального экстремума функции нескольких переменных
Выполнить самостоятельно.
В соответствии с индивидуальным заданием №1 решить задачу максимизации с использованием симплексных таблиц. Вариант задания выбирается по номеру зачетной книжки:
-предпоследняя цифра - № столбца
-последняя цифра – № строки
Пример: Симплексным методом решить задачу максимизации.
F(x)=5x1-x2→max
2x1-x2+x3≤3
3x1+2x2≤6
x1≥0; x2≥0
1. Переведем систему ограничений в канонический вид, введя выравнивающие (базисные) неизвестные – x3 и x4.
Задача принимает следующий вид:
2x1-x2+x3=3
3x1+2x2=6
x1=x2=x3=x4=0
F-5x1+x2=0
2. Заполняем первую симплексную таблицу. В ней x3, x4 – основные переменные (базис).
Таблица 1.
Базис |
Свободный член |
Переменные |
Оценочные отношения |
|||
х1 |
х2
|
х3
|
х4
|
|||
x3 |
3 |
2 |
-1 |
1 |
0 |
3/2 |
x5 |
6 |
3 |
2 |
0 |
1 |
2 |
F |
0 |
-5 |
1 |
0 |
0 |
|
← разрешающая строка
↑
разрешающий столбец
3. Проверяем выполнение критерия на max – первый опорный план не оптимальный, т.к. в F коэффициент при х1<0
4. Выбираем наибольший по модулю отрицательный коэффициент F, который определяет разрешающий столбец.
5. Делим свободные члены на коэффициенты разрешающего столбца (определяем оценочные отношения). Выбираем разрешающую строку, где это отношение минимально
Min (3, 2)= 3
2
Разрешающий элемент будет а11=2
Для построения таблицы 2 в качестве основной переменной выбираем х1, т.к. она образует разрешающий столбец таблицы 1.
Таблицы 2.
Базис |
Свободный член |
Переменные |
Оценочные отношения |
|||
х1 |
х2
|
х3
|
х4
|
|||
х1 |
3/2 |
1 |
-1/2 |
1/2 |
0 |
-3 |
х4 |
3/2 |
0 |
7/2 |
-3/2 |
1 |
3/7 |
F |
15/2 |
0 |
-3/2 |
5/2 |
0 |
|
← разрешающая строка
↑
разрешающий столбец
Построение таблицы №2.
1. Заменим переменные в базисе с х3 на х1.
2. Делим элементы разрешающей строки (табл.1) на разрешающий элемент, результаты занесем в строку х1 в табл. №2.
3. В остальных клетках разрешающего столбца (табл. 1) записываем 0.
4. Остальные клетки заполняем по правилу прямоугольника
НЭ=СТЭ-(АхВ)/РЭ
СТЭ А
В РЭ
5. Критерий эффективности опять не выполнен, т.к. F имеет два коэффициента -3 < 0
2
Построение таблицы №3
Таблицы 3.
Базис |
Свободный член |
Переменные |
Оценочные отношения |
|||
х1 |
х2
|
х3
|
х4
|
|||
х1 |
12/7 |
1 |
0 |
2/7 |
1/7 |
|
х2 |
3/7 |
0 |
1 |
-1/7 |
2/7 |
|
F |
57/7 |
0 |
0 |
13/7 |
3/7 |
|
Из таблицы 3 видно – критерий оптимальности выполнен.
Оптимальные базисные решения Х*=(12, 3,0,0)
7
F=5x1-x2=5x12- 3=57
7 7 7
ИНДИВИДУАЛЬНЫЕ ЗАДАНИЯ №1
|
0,1,2 |
3,4,5,6 |
7,8,9 |
0 1 |
f(x)=3x1-2x2; 2x1+x2≤2; x1+2x2≤2; x1,x2≥0; |
f(x)=4x1+8x2; x1+x2≤3; x1+2x2≤2; x1,x2≥0; |
f(x)=-4x1+8x2; 3x1+5x2≤15; x1-x3≤1; x1,x2≥0; |
2 3 |
f(x)=3x1-2x2; -x1+2x2≤2; 2x1-x2≤2; x1,x2≥0; |
f(x)=-x1+6x2; 4x1+3x2≤12; -x1+x2≤1; x1,x2≥0; |
f(x)=-x1+6x2; x1+x2≤3; -2x1+x2≤2; x1,x2≥0; |
4 5 |
f(x)=x1+6x2; 3x1+4x2≤12; -x1+x2≤2; x1,x2≥0; |
f(x)=2x1+6x2; 3x1+4x2≤12; x1+2x2≤2; x1,x2≥0; |
f(x)=8x1+12x2; 2x1+x2≤4; 2x1+5x2≤10; x1,x2≥0; |
6 7 |
f(x)=8x1+12x2; -3x1+2x2≤0; 4x1+3x2≤12; x1,x2≥0; |
f(x)=3x1-2x2; 2x1+x2≤2; 2x1+3x2≤6; x1,x2≥0; |
f(x)=6x1+4x2; x1+2x2≤2; -2x1+x2≤0; x1,x2≥0; |
8 9 |
f(x)=6x1+4x2; 3x1+2x2≤6; 3x1+x2≤3; x1,x2≥0; |
f(x)=8x1+6x2; -x1+x2≤1; 3x1+2x2≤6; x1,x2≥0; |
f(x)=8x1+6x2; -x1+x2≤2; 3x1+4x2≤12; x1,x2≥0; |
Пример. Решить следующую задачу максимизации
F(x)=5x1-x2→max
2x1-x2+x3≤3
3x1+2x2≤6
x1≥0; x2≥0
Запишем систему в каноническом виде
2x1-x2+x3=3
3x1+2x2=6
x1=0; x2=0
Решение. В системе координат Х1ОХ2 построим ОДР.
И з первого уравнения:
При Х1=0 X2=-3 → (X1,X2)=(0,-3) прямая (1)
При Х2=0 Х1=3 → (X1,X2)=(3, 0)
2 2
И з второго уравнения:
При Х1=0 X2=3 → (X1,X2)=(0,3) прямая (2)
При Х2=0 Х1=2 → (X1,X2)=(2, 0)
Из условия (3) (Х1, Х2)=(0,0) Х1=0 – прямая (3)
Х2=0 – прямая (4)
Для определения с какой стороны прямых ОДР все точки удовлетворяют неравенствам возьмем точку (1,1). Таким образом, ОДР лежит внутри АВСD.
Для определения точки выхода, строим прямую:
F(x)=5x1-x2=0
5x1=x2→ x1=x2=0 при x1=1, x2=5
Строим нулевую линию уровня целевой функции.
Строим нормальный вектор q (5.-1)
Передвигая линию уровня целевой функции параллельно себе, по вектору q, фиксируем ее крайнее положение (т.В), которая принадлежит прямым (1) и (2).
2x1-x2+x3=3
3x1+2x2=6
x1=12 x2=3, т.е. Х*=(12, 3)
7 7 7 7
F(Х*)=5х12 -3=57
7 7 7