- •Методические указания
- •Часть 2
- •Одесский припортовый завод
- •Содержание
- •Лабораторная работа №10. Исследование устойчивости Реактора
- •Материальный баланс
- •Уравнение кинетики
- •Лабораторная работа №11. Решение задач оптимизации с использованием аналитических методов
- •Лабораторная работа №12. Методы одномерного поиска для решения задач оптимизации
- •Лабораторная работа №13. Бе3градиентные методы решения задач оптимизации
- •1.2. Метод Хука-Дживса
- •Лабораторная работа №14. Симплексные методы решения задач оптимизации
- •Лабораторная работа №15. Градиентные методы решения задач оптимизации
- •Лабораторная работа №16. Использование методов линейного программирования для решения оптимизационных задач
- •Литература, рекомендуемая при подготовке к лабораторным работам
Лабораторная работа №16. Использование методов линейного программирования для решения оптимизационных задач
Цель работы:
– определить экстремум целевой функции с использованием симплекс-метода данцига.
1. ОПИСАНИЕ МЕТОДИКИ РАСЧЕТОВ
Линейное программирование представляет собой математический аппарат, разработанный для решения оптимизационных задач с линейным выражением для критерия оптимальности и линейными ограничениями на область изменения переменных.
Математически задача линейного программирования формулируется следующим образам: имеется ряд переменных х1, х2, …, хn. Требуется найти такие неотрицательные значения этих переменных, которые бы удовлетворяли системе линейных уравнений (ограничениям):
а11x1 + а12x2 + … + а1nxn ≥ b1
а21x1 + а22x2 + … + а2nxn ≥ b2 (16.1)
………………………………
аm1x1 + аm2x2 + … + аmnxn ≥ bm
и, кроме того, линейная функция достигала экстремума.
R = C1x1 + C2x2 + … + Cnxn . (16.2)
В системе (16.1) могут быть как равенства, так и неравенства. Для решения таких задач используется симплекс-метод данцига
Пусть система ограничений (16.1) записана в виде равенств. Полагая, что ранг матрицы коэффициентов меньше числа неизвестных, выразим r неизвестных системы через остальные:
x1
=
+
+
…+
+
x2
=
+
+
…+
+
(16.3)
……………………………………………
xr
=
+
+
…+
+
Неизвестные x1, x2, …xr называются базисными, остальные переменные – свободные.
Т.к. свободным неизвестным можно придавать
любые неотрицательные значения (xr+1≥0),
то свободные члены
,
,…,
должны быть положительными, поскольку
в противном случае (т.е. если
<0
при xr+1=0, xr+2=
0, … xn=0)
получим xj <0 (j
= 1, 2,…, r).
Положив все свободные неизвестные равными 0 (xr+1=0, xr+2=0,…, xn=0) получим:
x1= ; x2= , …, xr = . (16.4)
При этом вектор ( , ,…, , 0, 0,…, 0) называется базисным вектором.
Полученное решение (16.4) системы (16.1) является допустимым (x1>0; x2>0; …xr>0) и называется базисным (отвечает первому базису).
Подставив в линейную форму критерия оптимальности (16.2) вместо базисных переменных их выражения через свободные переменные, получим
R
=
+
+
+
…+
.
(16.5)
Причем для первого базисного решения (xr+i=0)
R = (16.6)
Далее переходим к другому базисному решению с таким расчетом, чтобы при этом значение линейной формы R не увеличивалось (в случае поиска минимума). Этот переход осуществляется следующим образом: заменяется одна из базисных переменных другой, которая раньше была небазисной (свободной) и т. д. Этот процесс повторяется до тех пор, пока R не достигнет экстремального значения.
2. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ
Работа выполняется в такой последовательности:
определяется ранг матрицы коэффициентов;
отыскивается первое базисное решение;
анализируется целевая функция и находится второе базисное решение;
заменяется одна из базисных переменных на небазисную и определяется новый вид целевой функции;
отыскивается очередное базисное решение и критерий оптимальности и т.д.
3. ИНДИВИДУАЛЬНЫЕ ЗАДАНИЯ
Определить оптимальные значения переменных, при которых целевая функция достигает экстремума и соблюдаются следующие ограничения:
Таблица 16.1
№ вар. |
Целевая функция |
Ограничения |
1 |
R = x4 – x5 Min |
x1 =1 – x4 + 2 x5 ; x2 = 2 + 2 x4 – x5 ; x3 = 6 – 3 x4 –x5 |
2 |
R = 6 – 2 x4 + 2 x5 Min |
2 x2 + 4 x4 + 6 x5 – 12 = 0; 2 x3 – 2 x4 – 6 x5 – 4 = 0; 2 x1 + 2 x4 + 2 x5 – 4 = 0 |
3 |
R = x1 – x2 Min |
x3 = 2 + x1 – x2 ; x4 = 5 + x1 – 2 x2 |
4 |
R = 3 x5 + 2 x4 + 3 Max |
x1 =6 – x4 +2 x5; x2 = 4 – 2 x4 – x5; x3 = 8 – x4 + 3 x5 |
5 |
R = x5 – x4 + 3 Min |
3 x5 + x1 + 2 x2 = 6; x3 – x4 – 3 x5 = 4; x4 + x5 + x1 = 4 |
6 |
R = – x5 + x4 Min |
3 – 3 x4 – x5 – x3 = 0; 1 – x4 + 2 x5 – x1 = 0; 2 + 2 x4 – x5 – x2 = 0 |
7 |
R = x1 +0,5 x2 + x3 Max |
2 x1 + x2 + x3 + x4 = 3; 3 x1 – x2 + x3 + x5 =4; 2 x1 – 2 x2 – x3 + x6 = 2 |
4. ОФОРМЛЕНИЕ ПРОТОКОЛА
В протоколе по лабораторной работе формулируется цель работы, описываются результаты расчетов по определению экстремума целевой функции.
5. КОНТРОЛЬНЫЕ ВОПРОСЫ
Классификация методов решения задач оптимизации.
Общая идея методов линейного программирования.
Симплекс–метод Данцига, алгоритм расчета.
