Реферат: Содержание введение
СОДЕРЖАНИЕ
Введение………………………………………………………………………....3
1. Основные понятия теории оптимизации……………………………...…5
1.1. Общая постановка задачи оптимизации…………………………………..5
1.2. Ограничения на допустимое множество……………………………...…6
1.3. Классическая задача оптимизации……………………………………….6
1.4. Функция Лагранжа…………………………………………………………6
2. Линейное программирование: формулировка задач и их
графическое решение………………………………………………………….8
2.1. Задача ЛП……………………………………………………………………8
2.2. Графическое решение задачи ЛП………………………………………….9
3. Алгебраический метод решения задач……………………………………11
3.1. Стандартная форма линейных оптимизационных моделей……………...11
3.2. Симплекс-метод………………………………………………………….....12
4. Двойственность………………………………………………………………25
Заключение……………………………………………………………………...29
Библиографический список ……………………………………………...…..30
Введение
Линейное программирование представляет собой наиболее часто используемый метод оптимизации. К числу задач линейного программирования можно отнести задачи: ^ · рационального использования сырья и материалов; задачи оптимального раскроя; · оптимизации производственной программы предприятий; · оптимального размещения и концентрации производства; · составления оптимального плана перевозок, работы транспорта; ^ · управления производственными запасами; · и многие другие, принадлежащие сфере оптимального планирования.
Для большого количества практически интересных задач целевая функция выражается линейно - через характеристики плана, причем допустимые значения параметров подчинены линейным равенствам или неравенствам. Нахождение при данных условиях абсолютного экстремума целевой функции носит название линейного программирования.
Прямая задача линейного программирования является математической формулировкой проблемы составления такого плана использования различных способов производства, который позволяет получить максимальное количество однородного продукта при имеющихся в наличии ресурсах.
Математическое программирование - это прикладная отрасль математики, которая является теоретической основой решения задач оптимального планирования.
Существуют следующие разделы математического программирования: линейное, параметрическое, нелинейное и динамическое программирование. Наиболее разработанным и широко применяемым разделом математического программирования является линейное программирование, целью которого служит отыскивание оптимума (max, min)заданной линейной функции при наличии ограничений в виде линейных уравнений или неравенств.
Наиболее известным и широко применяемым на практике для решения общей задачи линейного программирования (ЛП) является симплекс-метод. Несмотря на то, что симплекс-метод является достаточно эффективным алгоритмом, показавшим хорошие результаты при решении прикладных задач ЛП, он является алгоритмом с экспоненциальной сложностью. Причина этого состоит в комбинаторном характере симплекс-метода, последовательно перебирающего вершины многогранника допустимых решений при поиске оптимального решения.
Симплекс метод - является универсальным методам, которым можно решить любую задачу линейного программирования.
^ 1. Основные понятия теории оптимизации
1.1. Общая постановка задачи оптимизации
В общей задаче требуется найти вектор
из допустимой области , который обращает в минимум целевую функцию q(x), т.е. такой , для которого
(1)
Если существует, то он определяет слабый, глобальный (абсолютный) минимум q*(x) в допустимой . Слабый, т.к. удовлетворяет нестрогому неравенству. Глобальный, т.к. неравенство справедливо для любых x из области X. Минимум при сильный, если для. Если поменять знаки неравенств – получим сильный и слабый максимумы. Минимум в точке называется локальным (относительным), если найдётся такая окрестность O(x*) точки , что для всех имеет место . Если дифференцируема, то задача отыскания локальных минимумов сводится к нахождению стационарных точек, в которых обращаются в ноль частные производные q(x):
(2)
(2) – необходимое, но не достаточное условие. Достаточным условием существования в стационарной точке относительного минимума является положительная определённость квадратичной формы.
^ 1.2. Ограничения на допустимое множество
Теорема Вейерштрасса: непрерывная функция, определённая на непустом замкнутом ограниченном множестве, достигает минимума (максимума) по крайней мере в одной из точек этого множества.
^ 1.3. Классическая задача оптимизации
Состоит в нахождении минимума целевой функции, где – точка в пространстве при начальных ограничениях типа равенств
(3)
Если (3) имеют место, то минимум q(x) называется условным минимумом. Если ограничения (3) отсутствуют, то говорят о безусловном минимуме.
Классический способ решения данной задачи состоит в том, что (3) используют для исключения из рассмотрения переменных. При этом целевая функция приводится к виду
(4)
,где через обозначены неисключаемые переменные. Задача теперь состоит в нахождении значений , которые обращают в минимум q1 и на которые не наложено ограничений (задача на безусловный экстремум).
^ 1.4. Функция Лагранжа
Введём в рассмотрение вектор и исследуем свойства функции
(5)
– функция Лагранжа, - множители Лагранжа.
– функция n+m переменных .
Рассмотрим стационарные точки функции , которые получим, приравняв к нулю частные производные по и по :
(6)
(7)
Если в стационарной точке (x*, y*) функция достигает минимума, то обеспечивает минимум функции q(x) и при выполнении ограничений (3), т.е. даёт решение задачи.
Задача на условный минимум целевой функции q(x) при наличии ограничений типа равенств сводится к задаче на определение стационарных точек функции Лагранжа .
^ 2. Линейное программирование: формулировка задач и их графическое решение
2.1. Задача ЛП
Рассмотрим на примере задачи фирмы. Небольшая фабрика изготовляет два вида красок: для наружных (E) и внутренних (I) работ. Продукция поступает в оптовую продажу. Для производства красок используется два исходных продукта – A и B. Максимально возможные суточные запасы этих продуктов составляют 6т и 8т соответственно. Расходы A и B на производство 1т соответствующих красок приведены в таблице.
Исходный продукт
Расход на тонну краски
Максимальный запас, т.
краска E
краска I
A
1
2
6
B
2
1
8
Суточный спрос на краску I никогда не превышает спроса на краску E более чем на 1т. Спрос на I не превышает 2т. Оптовая цена за 1т краски E – 3000$, I – 2000$. Какое количество краски каждого вида фабрика должна производить, чтобы доход от реализации продуктов был максимальным?
Так как нужно определить объём производства каждого вида краски, переменными в модели являются:
xE – суточный объём производства краски E (в тоннах);
xI – суточный объём производства краски I (в тоннах).
Обозначив доход (в тыс. $) через , можно дать математическую формулировку целевой функции: определить (допустимые) значения xE и xI, максимизирующие величину общего дохода
Ограничения на расход исходных продуктов:
(для A)
(для B)
Ограничения на величину спроса на продукцию:
Потребуем выполнения условия неотрицательности переменных:
Получили математическую модель:
Определить суточные объёмы производства (в т.) краски I и E, при которых достигается
(целевая функция)
при ограничениях
^ 2.2. Графическое решение задачи ЛП
Построим область допустимых решений, в которой одновременно выполняются все ограничения. Искомое пространство решений – многоугольник ABCDEF. Пространство решений содержит бесконечное число точек, являющихся допустимыми решениями, но, несмотря на это, можно найти оптимальное решение, если выяснить, в каком направлении возрастает целевая функция модели z=3xE+2xI. На график наносят ряд параллельных линий, соответствующих уравнению целевой функции при нескольких произвольно выбранных и последовательно возрастающих значениях , что позволяет определить наклон целевой функции и направление её увеличения. На видно, что оптимальному решению соответствует точка C, являющаяся пересечением прямых
Решив систему, получим
Тогда получаемый доход
тыс $.
Оптимальному решению всегда соответствует одна из допустимых угловых точек пространства решений. Какая из этих точек окажется оптимальной, зависит от наклона прямой, представляющей целевую функцию (т.е. от коэффициентов целевой функции).
^ 3. Алгебраический метод решения задач
Использование графического метода удобно при решении задач ЛП с двумя переменными. При большем их числе необходимо применение алгебраичского аппарата.
Процесс решения задачи ЛП симплекс-методом носит итерационный характер: однотипные вычислительные процедуры в определённой последовательности выполняются до тех пор, пока не будет получено оптимальное решение.
^ 3.1. Стандартная форма линейных оптимизационных моделей
Все ограничения записываются в виде равенств с неотрицательной правой частью. Исходное ограничение можно представить в виде равенства, прибавляя остаточную переменную к левой части ограничения (вычитая избыточную переменную из левой части).
Например,
Введём остаточную переменную s1>0, тогда
Правую часть равенства можно сделать неотрицательной, умножив обе части на –1.
Значения всех переменных модели неотрицательны.
Любую переменную yi, не имеющую ограничения в знаке, можно представить как разность двух неотрицательных переменных:
При любом допустимом решении только одна из этих переменных может принимать положительное значение, т.к. если , то , и наоборот. Это позволяет рассматривать как остаточную переменную, а – как избыточную.
Целевая функция подлежит максимизации или минимизации.
Максимизация некоторой функции эквивалентна минимизации той же функции, взятой с противоположным знаком, и наоборот.
Эквивалентность означает, что при одной и той же совокупности ограничений оптимальные значения переменных в обоих случаях будут одинаковы.
3.2. Симплекс-метод
Общую идею симплекс-метода проиллюстрируем на примере. На исходная точка алгоритма – начало координат (т. A) – начальное решение. От исходной точки осуществляется переход к некоторой смежной угловой точке (т. B или т. F). Её выбор зависит от коэффициентов целевой функции. Т.к. коэффициент при xE больше коэффициента при xI, а целевая функция подлежит максимизации, требуемое направление перехода соответствует увеличению xE (т. B). Далее указанный процесс повторяется для выяснения, существует ли другая экстремальная точка, соответствующая лучшему допустимому решению.
Правила выбора экстремальной точки:
Каждая последующая угловая точка должна быть смежной с предыдущей.
Обратный переход к предшествующей экстремальной точке не может производиться.
Чтобы описать рассмотренные процедуры формальными способами, необходимо определить пространство решений и угловые точки алгебраически. Требуемые соотношения устанавливаются по таблице:
Геометрическое определение (графический метод)
Алгебраическое определение (симплекс-метод)
Пространство решений
Ограничения модели стандартной формы
Угловые точки
Базисные решения задачи в стандартном виде
с
Модель:
максимизировать целевую функцию
при ограничениях
На – пространство решений. Каждую точку можно определить с помощью
Для идентификации нужной точки воспользуемся тем, что при ограничения модели эквивалентны равенствам, которые представляются соответствующими рёбрами пространства решений.
Анализируя , заметим, что
можно упорядочить, исходя из того, какое значение (нулевое или ненулевое) имеет данная переменная в экстремальной точке.
Экстр.
переменные
точка
нулевые
ненулевые
A
B
C
D
E
F
Из таблицы выделим закономерности:
Стандартная модель содержит четыре уравнения и шесть неизвестных, поэтому в каждой из экстремальных точек (6–4) переменные должны иметь нулевые значения.
Смежные экстремальные точки отличаются только одной переменной в каждой группе (нулевых и ненулевых переменных).
Если линейная модель стандартной формы содержит уравнений и неизвестных, то все допустимые экстремальные точки определяются как все однозначные неотрицательные решения системы уравнений, в которых m-n переменных равны нулю. Однозначные решения такой системы – базисные решения. Если они удовлетворяют требованию неотрицательности правых частей, то это – допустимое базисное решение. Переменные, равные нулю – небазисные, остальные – базисные. Каждую следующую экстремальную точку можно определить определить путём замены одной из текущих небазисных переменных текущей базисной переменной. В нашем примере при переходе из т. A в т. B необходимо увеличивать небазисную переменную от исходного нулевого значения до значения, соответствующего т. B. В т. B s2 обращается в нуль (становится небазисной). Т.о., происходит взаимообмен xE и s2 между небазисными и базисными переменными.
Включаемая переменная – небазисная в данный момент переменная, которая будет включена в множество базисных переменных на следующей итерации. Исключаемая переменная – базисная в данный момент переменная, которая на следующей итерации подлежит исключению из множества базисных переменных.
^ Вычислительные процедуры симплекс-метода
Симплекс-алгоритм:
Шаг 0: Используя линейную модель стандартной формы, определяют НДБР путём приравнивания к нулю n-m (небазисных) переменных.
Шаг 1: Из числа текущих небазисных переменных выбирается включаемая в новый базис переменная, увеличение которой обеспечивает улучшение значения целевой функции. Если её нет -- текущее базисное решение оптимально, иначе переход к Шагу 2.
Шаг 2: Из числа переменных текущего базиса выбирается исключаемая переменная, которая должна стать небазисной при введении в состав базиса новой переменной.
Шаг 3: Находится новое базисное решение, соответствующее новым составам базисных и небазисных переменных. Переход к Шагу 1.
Если xE=xI=0, то
(соответствует точке A Ha ) – начальное допустимое решение.
Решение
-3
-2
2
6
2
8
-1
2
Если в задаче максимизации все небазисные переменные в -уравнении имеют неотрицательные коэффициенты, полученное пробное решение является оптимальным. Иначе в качестве новой базисной переменной следует выбрать ту, которая имеет наибольший по абсолютной величине отрицательный коэффициент. Применяя это условие к исходной таблице – переменная, включаемая в базис.
Процедура выбора подключаемой переменной предполагает проверку условия допустимости, требующего, чтобы в качестве исключаемой переменной выбиралась та (из текущего базиса), которая первой обращается в нуль при увеличении включаемой переменной вплоть до значения, соответствующего смежной экстремальной точке.
Отношение, идентифицирующее исключаемую переменную, можно определить из симплекс-таблице. Для этого в столбце вводимой переменной вычёркиваются отрицательные и нулевые элементы ограничений. Затем вычисляются отношения постоянных из правых частей ограничений к оставшимся элементам столбца. Исключаемая переменная – та, для которой это отношение минимально.
Решение
Отношение
-3
-2
-
2
6
2
8
-1
-
2
-
Столбец, ассоциированный с вводимой переменной – ведущий столбец; строка, соответствующая исключаемой переменной – ведущая строка; на их пересечении – ведущий элемент.
Поиск нового базисного решения осуществляется методом исключения переменных (метод Жордана-Гаусса). Этот процесс включает в себя вычислительные процедуры двух типов.
Тип 1. Формирование ведущего уравнения
Новая ведущая строка = предыдущая ведущая строка/ведущий элемент
Тип 2. Формирование остальных уравнений
Новое уравнение = предыдущее уравнение – (коэффициент ведущего столбца предыдущего уравнения)*(новая ведущая строка)
Новая симплекс-таблица, полученная после проведения рассмотренных операций:
Решение
Отношение
-
-
-
-
-
-
xI – вводимая переменная (т.к. коэффициент в -уравнении -1/2). Исключаемая переменная s1, (отношение 4/3 – минимальное). Снова проведём вычисления двух типов. Последняя симплекс-таблица соответствует оптимальному решению задачи, т.к. в -уравнении ни одна из небазисных переменных не фигурирует с отрицательными коэффициентами.
В случае минимизации целевой функции в этом алгоритме необходимо изменить только условие оптимальности: в качестве новой базисной переменной следует выбирать переменную, которая в -уравнении имеет наибольший положительный коэффициент. Условия допустимости в обоих случаях одинаковы.
^ Искусственное начальное решение
Для получения начального базисного решения мы использовали остаточные переменные. Однако когда исходное ограничение записано в виде равенства или имеет знак, нельзя сразу же получить НДБР. Поэтому были разработаны два метода, в которых используется «штрафование» искусственных переменных.
Метод больших штрафов (М-метод)
Рассмотрим линейную модель, приведённую к стандартной форме:
минимизировать
при ограничениях
В первом и втором уравнениях нет переменных, исполняющих роль остаточных. Поэтому введём в каждое из этих уравнений по одной из искусственных переменных R1 и R2:
За использование этих переменных в составе целевой функции можно ввести штраф, приписывая им достаточно большой положительный коэфффициент . Получим линейную модель:
минимизировать
при ограничениях
Теперь если
,то начальное допустимое решение:
Метод оптимизации, направленный на нахождение минимального значения , приведёт к тому, что переменные R1 и R2 в оптимальном решении обратятся в нуль.
Необходимо переформулировать условие задачи, чтобы представить процесс решения в удобной табличной форме. Подставив в целевую функцию полученные из соответствующих ограничений выражения для искусственных переменных
получим выражение для :
Решение представлено в сводной таблице:
Итерация
Решение
Отношение
-
-
-
-
-
-
4. Двойственность
Двойственная задача – вспомогательная задача ЛП, формулируемая с помощью определённых правил непосредственно из исходной, или прямой задачи.
Прямая задача ЛП в стандартной форме:
максимизировать (минимизировать)
при ограничениях
В состав включаются избыточные и остаточные переменные.
Для формулировки двойственной задачи расположим коэффициенты прямой задачи согласно схеме:
каждому ограничению прямой задачи соответствует переменная двойственной задачи;
каждому переменной прямой задачи соответствует ограничение двойственной задачи;
коэффициенты при некоторой переменной, фигурирующие в ограничения прямой задачи, становятся коэффициентами левой части соответствующего ограничения двойственной задачи, а коэффициент, фигурирующий при той же переменной в выражении для целевой функции прямой задачи, становится постоянной правой части этого же ограничения двойственной задачи.
Информация о других условиях двойственной задачи (направление оптимизации, ограничения и знаки двойственных переменных) представлена в таблице:
Прямая задача в стандартной форме.
Двойственная задача
Целевая функция
Целевая функция
Ограничения
Переменные
Максимизация
Минимизация
Не ограничены в знаке
Минимизация
Максимизация
Не ограничены в знаке
Рассмотрим пример:
Прямая задача:
максимизировать
при ограничениях
Прямая задача в стандартной форме:
максимизировать
при ограничениях
Двойственная задача:
минимизировать
при ограничениях:
(означает, что y1>0). y1, y2 не ограничены в знаке.
Заключение
В настоящее время линейное программирование является одним из наиболее употребительных аппаратов математической теории оптимального принятия решений, в том числе и в финансовой математике. Для решения задач линейного программирования разработано сложное программное обеспечение, дающее возможность эффективно и надежно решать практические задачи больших объемов. Эти программы и системы снабжены развитыми системами подготовки исходных данных, средствами их анализа и представления полученных результатов. В развитие и совершенствование этих систем вложен труд и талант многих математиков, аккумулирован опыт решения тысяч задач. Владение аппаратом линейного программирования необходимо каждому специалисту в области прикладной математики.
Библиографический список
1.Онегов, В.А. Исследование операций. Задачи, методы, алгоритмы. – Киров: ВГПУ, 2001.
2. Сборник задач и упражнений по высшей математике: мат. программирование: Учеб. Пособие/ А.В. Кузнецов, В.А. Сакович, Н.И. Холод; Под. общ. ред. А.В. Кузнецова – Мн.: Выш. шк., 2002
3. Ф. П. Васильев, А. Ю. Иваницкий Линейное программирование - Издательство: Факториал Пресс, 2003.
4. Н. И. Глебов, Ю. А. Кочетов, А. В. Плясунов, Методы оптимизации учебное пособие , Новосибирск: Новосибирский государственный университет, 2000
5. http://simplex-metod.narod.ru/
еще рефераты
Еще работы по разное