Линейное программирование: решение задач графическим способом
рать
ограничение [pic]. Оставшаяся часть будет неограниченным выпуклым
многоугольником.
[pic]
Наконец, возможен случай, когда неравенства (1.31) противоречат друг
другу, и допустимая область вообще пуста.
Рассмотрим теорию на конкретном примере:
Найти допустимую область задачи линейного программирования,
определяемую ограничениями
|[pic] |(1.3|
| |2) |
[pic]
Решение:
1. Рассмотрим прямую [pic]. При [pic], а при [pic]. Таким образом, эта
прямая проходит через точки (0,1) и (-1,0). Беря [pic]получим, что
-0+0<1 и поэтому интересующая нас полуплоскость лежит ниже прямой,
изображенной на рис. 4.а.
2. Рассмотрим прямую [pic]. При [pic], а при . Таким образом, эта прямая
проходит через точки (0, -1/2) и (1,0). так как [pic](4.б).
3. Наконец, рассмотри[pic]м прямую [pic]. Она проходит через точки (0,3)
и (3,0) и так как 0+0<3, то интересующая нас полуплоскость лежит ниже
прямой, изображенной на рис. 4.в.
[pic]
Сводя все вместе и добавляя условия [pic]получим рисунок 5, где
выделена область, в которой выполняются одновременно все ограничения
(1.32). Обратите внимание на то, что получившаяся область имеет вид
выпуклого многоугольника.
Этап 2.
Вернёмся теперь к исходной задаче линейного программирования. В ней,
кроме системы неравенств, есть еще целевая функция [pic].
[pic]
Рассмотрим прямую[pic]. Будем увеличивать L. Что будет происходить с
нашей прямой?
Легко догадаться, что прямая будет двигаться параллельно самой себе в
том направлении, которое дается вектором [pic], так как это - вектор
нормали к нашей прямой и одновременно вектор градиента функции [pic].
А теперь сведем всё вместе. Итак, надо решить задачу
[pic]
[pic]
Oграничения задачи вырезают на плоскости некоторый многоугольник. Пусть
при некотором L прямая [pic]пересекает допустимую область. Это пересечение
дает какие-то значения переменных [pic], которые являются планами.
Этап 3
Увеличивая L мы начнем двигать нашу прямую и её пересечение с
допустимой областью будет изменяться (см. рис. 7). В конце концов эта
прямая выйдет награницу допустимой области - как правило, это будет одна из
вершин многоугольника. Дальнейшее увеличение L приведёт к тому, что
пересечение
[pic]
прямой [pic]с допустимой областью будет пустым. Поэтому то положение прямой
[pic], при котором она вышла на граничную точку допустимой области, и даст
решение задачи, а соответствующее значение L и будет оптимальным значением
целевой функции.
1.4 Примеры задач, решаемых графическим методом.
Пример:
Решить задачу
|[pic] |(1.4|
| |1) |
Решение
Допустимую область мы уже строили - она изображена на рис. 5.
Повторим еще раз этот рисунок, оставив только допустимую область и
нарисовав дополнительно прямые [pic] (см. рис. 8).
[pic]
Пусть, например, L=2. Тогда прямая [pic]проходит через точки (2,0) и
(0,1) и изображена на рис. 8. Будем теперь увеличивать L. Тогда прямая
начнёт двигаться параллельно самой себе в направлении, указанном стрелкой.
Легко догадаться, что максимальное значение L получится тогда, когда прямая
пройдет через вершину многоугольника, указанную на рисунке, и дальнейшее
увеличение L приведет к тому, что прямая выйдет за пределы многоугольника и
её пересечение с допустимой областью будет пустым.
Выделенная вершина лежит на пересечении прямых
[pic]
и поэтому имеет координаты [pic]. Это и есть решение нашей задачи, т.е.
[pic]есть оптимальный план задачи (1.41). При этом значение целевой функции
[pic], что и дает её максимальное значение.
Обратите внимание на то, что оптимальный план, как правило,
соответствует какой-то вершине многоугольника, изображающего допустимую
область. И лишь в том случае, когда прямая [pic]случиться так, что решение
не будет единственным. Но и в этом случае вершины, соответствующие границам
этой стороны, дают оптимальные планы нашей задачи линейного
программирования. Таким образом, вершины допустимой области играют в
решении задач линейного программирования особую роль.
[pic]
Ну, а если допустимая область неограничена, то и значение целевой
функции может быть неограниченным.
Подводя итог этим примерам, можно сформулировать следующие положения:
1. допустимая область - это выпуклый многоугольник;
2. оптимум достигается в вершине допустимой области (если допустимая
область ограничена и не пуста);
3. ограниченность целевой функции в допустимой области является
необходимым и достаточным условием разрешимости задачи.
Гл 2 Решение задач линейного программирования графическим способом на ЭВМ
2.1 Описание работы программы
Программа написана с использованием собственных функций и процедур и
трех стандартных модулей System, Crt и Graph.
При запуске программы она проверяет возможно ли использование
графического интерфейса. Если это возможно то программа переходит к
следующему этапу.
Далее процедурой ShowXOY Рисуется на экран координатные оси. На этом
работа этой процедуры заканчивается и пользователь в следующей процедуре
(EnterNerav и в частности в подпроцедуре GetNerav) предлагается ввести
коэффициенты неравенства a1x+a2y=b в следующем порядке: a1 пробел a2 пробел
b. Сразу после ввода всех коэффициентов процедурой ShowLine рисуется нужная
линия. После нажатия [Esc] процедура EnterNerav заканчивается и передает
управление процедуре EnterMainF в которой пользователю предлагается ввести
коэффициенты целевой функции. Далее работа переходит к процедуре GetResult
где идет подсчет оконцательного товета с помощбю процедуры SolveOprtel где
считаетя определитель т. е. точки пересечения целевой функции с каждой
линией ограничения. Далее выводится ответ, если это возможно.
Далее следует описание используемых стандартных процедур и функций.
Процедуры и функции модуля System:
Function Frac(X : Real) : Real;
Возвращает дробную часть аргумента.
Параметр X - выражение вещественного типа. Результат - дробная часть X, то
есть Frac(X) = X-Int(X).
Procedure Str(X [: Width [: Decimals ]]; Var S : String);
Преобразовывает число в строку. Преобразовывает числовое значение X в
строковое представление этого числа, которое можно выводить операторами
типа Write и OutText.
Function Round(X : Real) : Longint;
Округляет значение вещественного типа до значения целочисленного типа.
X - выражение с реальным типом. Round возвращает значение типа Longint,
которое является значением X, округленного к самому близкому целому числу.
Если X – ровно посередине между двумя целыми числами, то результатом будет
число с самой большой абсолютной величиной.
Если округленное значение X ненаходится внутри допустимого диапазона
Longint, то происходит ошибка во время выполнения программы.
Модуль Crt:
В модуле Crt находятся мощные подпрограммы, которые дают вам
возможность полного управления возможностями вашего PC.
Подпрограммы модуля Crt обеспечивают контроль над текстовыми режимами
экрана, расширенными кодами клавиатуры, цветами, окнами и звуком.
Crt может использоваться только в программах, предназначенных для IBM
PC, AT, PS/2 и полностью совместимых.
Procedure Read( [ var F : Text; ] V1 [ , V2, …, VN ]); (текстовые
файлы)
Читает одну или более величин из текстового файла в одну или более
переменных. Параметры: F - необязательная переменная текстового файла, если
не указана, то используется стандартная переменная Input; V1,...,VN -
переменные типа Char, Integer, Real или String.
В случае переменной типа Char процедура Read считывает из файла один
символ и присваивает его переменной. В случае переменной целого типа
процедура Read ожидает поступления последовательности символов, образующих
число со знаком, согласно принятому в Паскале синтаксису. Любые пробелы,
знаки табуляции или метки конца строки, предшествующие числовой строке,
пропускаются. Считывание прекращается при обнаружении первого пробела,
символа табуляции или метки конца строки, которые следуют за числовой
строкой, или в том случае, если функция Eof(F) принимает значение True.
Если числовая строка не соответствует ожидаемому формату, то происходит
ошибка ввода-вывода, в противном случае переменной присваивается считанное
значение. Если Eof(F) принимала значение True перед выполнением процедуры
Read, или Eof(F) приняла значение True при пропуске начальных пробелов,
знаков табуляции или меток конца строки, то переменной присваивается
нулевое значение. Следующая операция Read начнется с пробела, символа
табуляции или метки конца строки, которыми завершилась числовая строка.
В случае переменной вещественного типа процедура Read ожидает
поступления последовательности символов, которые образуют число со знаком в
соответствии с принятым в Паскале синтаксисом за исключением того, что
шестнадцатиричное представление не до
| |  скачать работу |
Линейное программирование: решение задач графическим способом |
