Новосибирский Государственный Университет
Опубликован: 06.11.2008 | Доступ: свободный | Студентов: 3338 / 1210 | Оценка: 4.24 / 3.67 | Длительность: 14:37:00
Специальности: Математик
Лекция 6:

Дискретное программирование

< Лекция 5 || Лекция 6: 1234 || Лекция 7 >

Преобразуем матрицу (6.11.) в соответствии с указанным выше требованием получения базисных значений переменных величин. Для этого необходимо выполнить над ней такие преобразования, чтобы базисные переменные остались по одной в каждом из уравнений (строке матрицы), а коэффициенты при них были равны единице. Начинаем с коэффициента при x_1 в первом уравнении. Чтобы сделать его равным единице, делим все коэффициенты первого уравнения на четыре. Для исключения переменной x_1 из остальных уравнений отнимаем от каждого из них первое уравнение, умноженное на такое число, при котором разность коэффициентов при x_1 была бы равна нулю. Например, второе и третье уравнения (строки) нужно умножить на нуль, четвертое – на единицу.

В результате преобразований получим

\begin{pmatrix}
1&0&0&\frac{1}{4}&0&0&4&\\
0&2&0&0&1&0&10&\\
0&0&1&2&6&0&76&\\
0&3&0&-1&0&1&8&\\
\end{pmatrix}

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

\begin{pmatrix}
1&0&0&\frac{1}{4}&0&0&4&\\
0&1&0&0&\frac{1}{2}&0&5&\\
0&0&1&2&6&0&76&\\
0&0&0&-1&-\frac{3}{2}&1&-7&\\
\end{pmatrix}

Для переменной x_3 — в третьей строке и x_4 — в четвертой:

\begin{pmatrix}
1&0&0&0&-\frac{3}{8}&\frac{1}{4}&\frac{9}{4}&\\
0&1&0&0&\frac{1}{2}&0&5&\\
0&0&1&0&3&2&62&\\
0&0&0&1&\frac{3}{2}&-1&7&\\
\end{pmatrix} ( 6.12)

Выполненная процедура носит название метода полного исключения (так называемое Жарданово исключение ) .

Теперь, приравнивая переменные x_5 и x_6 (соответственно пятый и шестой столбцы матрицы) нулю, можем написать значение базисных переменных, которые будут в этом случае равны свободным членам соответствующих уравнений:

x_1=\frac{9}{4};\\ 
x_2=5;\\
x_3=62;\\
x_4=7;\\

Обращаясь к геометрической интерпретации (рис. 6.1), можно убедиться, что полученные координаты (x_1=2\frac{1}{4}, x_2=5, x_5=x_6=0) соответствуют вершине A многоугольника ABCDEF — области допустимых планов. Это и есть первый допустимый план.

Теперь можно перейти ко второму шагу симплекс-метода – установлению того, является ли допустимый план, соответствующий найденной вершине A, оптимальным.

Наиболее естественным путем решения этой задачи был бы сплошной перебор всех вершин области допустимых планов, определение для каждой из них значений переменных x_j(j=1,2,...,6) и вычисление по ним в каждой вершине величины целевой функции.

Та вершина, в которой величина y окажется минимальной, и даст искомый оптимальный план.

Но этот путь весьма неэкономный, ибо требует просчитать большое количество планов, в том числе и явно неоптимальных.

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

В чем сущность направленного перебора?

В первом допустимом плане, соответствующем вершине A, целевая функция в соответствии с формулой (6.7.) равна:

y_A=-2,4\cdot 2\frac{1}{4}+0,8\cdot 5+22,8=21,4

Мы уже знаем из (6.10.), что y_{min}=13,2. Следовательно, целевая функция в точке A значительно больше минимума, и необходимо продолжать перебор вершин-планов до тех пор, пока не придем к оптимальному.

Из вершины A можно перейти к соседним вершинам F и B (рис. 6.1), двигаясь по сторонам многоугольника AF и AB соответственно. Видимо, нужно избрать такое направление перехода к соседней вершине, которое приведет к наибольшему уменьшению целевой функции.

Рассчитаем значения целевой функции для соседних вершин F,B. Пользуясь формулой (6.7.) и подставляя соответствующие значения x_1, x_2, получим:

y_F= -2,4\cdot 0+0,8\cdot 5=26,8;\\ y_B= -2,4\cdot 4+0,8\cdot 2\frac{2}{3}+22,8=15,33.

Сопоставляя два последних выражения, нетрудно убедиться, что минимизация функции цели достигается при движении к точке B по стороне AB. Это означает, что в базис вводится переменная x_5, которая в вершине A была равна нулю.

Поскольку при отсутствии наглядного геометрического представления заранее нельзя располагать значениями переменных в вершинах многоугольника, то для установления необходимости и направления перебора планов пользуются специальным критерием \delta_j:

\delta_j=\sum\limits_{i=1}^{m}c_ia_{ij}-c_j ( 6.13)
,

где индекс j приписывается небазисным (нулевым) переменным, а индекс i — базисным. Имеется доказательство того, что в случае оптимальности полученного плана все \delta_j становятся равными нулю или меньшими нуля. Включению в базис подлежит та переменная, для которой \delta_j принимает наибольшее положительное значение. В нашем примере это:

\delta_5=-\frac{8}{5}\cdot (\frac{3}{8})+(-\frac{16}{10})\cdot \frac{1}{2}+\frac{2}{10}\cdot 3+\frac{8}{10}\cdot\frac{3}{2}-0=\frac{13}{10};

\delta_6=-\frac{8}{10}\cdot \frac{1}{4}-\frac{16}{10}\cdot 0+\frac{2}{10}\cdot 2-\frac{8}{10}\cdot 1-0=-\frac{3}{5}.

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

Далее необходимо установить, какая переменная должна быть выведена из базиса при введении в него переменной x_5. Чтобы ответить на этот вопрос, будем рассуждать так.

Очевидно, следует переместиться по стороне AB как можно дальше от точки A, чтобы как можно больше уменьшить целевую функцию. Стало быть, можно взять в качестве координаты x_1 точки B ее максимальное возможное значение, которое допускается системой уравнений , соответствующей матрице (6.12.), то есть такое, при котором ни одна из переменных не становится отрицательной. Можно показать, что это достигается в том случае, если вывести из базиса переменную, которой соответствует минимальное положительное значение отношения свободного члена уравнения к коэффициенту при x_5 в соответствующем столбце матрицы (6.12.).

Поэтому избирается четвертая строка матрицы и соответственно переменная x_4, подлежащая исключению из базиса.

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

Получим

\begin{pmatrix}
1&0&0&0&0&0&4\\
0&1&0&0&0&\frac{1}{3}&\frac{8}{3}\\
0&0&1&0&0&4&48\\
0&0&0&\frac{2}{3}&1&-\frac{2}{3}&\frac{14}{3}
\end{pmatrix} ( 6.14)

Данной матрице отвечает допустимый план в вершине B. Приравнивая небазисные переменные нулю (x_4=x_6=0), получаем значения остальных переменных, соответствующих второму плану:

x_1=4; x_3=48;\\ x_2=\frac{8}{3};x_5=\frac{14}{3}.

Как уже было показано, y_B=15\frac{1}{3}. Итак, получено существенное сокращение целевой функции, однако критерий \delta_6 продолжает оставаться положительным, что говорит о необходимости дальнейшего улучшения плана:

\delta_4=-\frac{2}{5}; \delta_6=\frac{3}{10}.

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

После преобразования матрицы (6.14.) получаем матрицу (6.15.), отвечающую третьему плану:

\begin{pmatrix}
1&0&0&0&0&0&4\\
0&3&0&0&0&1&8\\
0&-12&1&0&0&0&16\\
0&2&0&\frac{2}{3}&1&0&10
\end{pmatrix} ( 6.15)
.

Данный план соответствует вершине C:

x_1=4; x_4=0;\\
x_2=0; x_5=10;\\
x_3=16; x_6=8.

Критерии для данного плана равны:

\delta_2=-\frac{4}{5}; \delta_4 =-\frac{2}{5}.

Поскольку нет ни одного положительного критерия, то третий план и является оптимальным, что следует из рис. 6.1.

Целевая функция при данном плане равна:

y_C=13,2

Итак, мы пришли аналитическим путем к тому же оптимальному плану, который был ранее получен геометрическим способом.

Решение примера 6.1. можно сформулировать следующим образом.

Чтобы общие потери были минимальны, количество носителей первого типа должно быть равно 4, второго – 0, третьего – 16, четвертого 0, пятого 10, шестого – 8.

При этом потери в носителях будут составлять 13,2 единицы.

Важно отметить, что наихудший план распределения, соответствующий точке F, приводит к потерям 26,8 единиц. Таким образом, без дополнительного расходования носителей и овощей (только за счет их рационального распределения) улучшен результат решения задачи перевозки более чем в два раза.

Полученное решение – целочисленное, что, однако, не является обязательным для задач линейного программирования. В ряде случаев для получения целочисленного решения приходится прибегать к искусственным приемам, что существенно усложняет вычислительную процедуру.

< Лекция 5 || Лекция 6: 1234 || Лекция 7 >
Данил Комардин
Данил Комардин

мне задали дистанционное задание на сертификат,но я не могу его найти