Учебное пособие: Методические указания по изучению дисциплины и выполнению контрольной работы
Название: Методические указания по изучению дисциплины и выполнению контрольной работы Раздел: Остальные рефераты Тип: учебное пособие | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ЭКОНОМИКО – МАТЕМАТИЧЕСКИЕ МЕТОДЫ Методические указания по изучению дисциплины и выполнению контрольной работы Методические указания по изучению дисциплины подготовил к.т.н., проф. кафедры математических методов обработки информации, зав. каф. ММОИ Клетин В.А. Москва, 2007 Экономико-математические методы Введение Экономические проблемы, возникающие перед специалистами, в большинстве своем сложные. Они зависят от множества различных, иногда противоречащих друг другу факторов, изменяются с течением времени и влияют на другие проблемы и процессы. Вследствие этого исследование экономической проблемы целесообразно проводить на адекватной математической модели. Математическая модель отражает проблему в абстрактной форме и позволяет учесть большое число разнообразных характеристик, от которых эта проблема зависит. Анализ и расчет математической модели позволяют выбрать оптимальные решения поставленной задачи и обосновать этот выбор. Впервые математические модели были использованы для решения практической задачи в 30-х годах в Великобритании при создании системы противовоздушной обороны. Для разработки данной системы были привлечены ученые различных специальностей. Система создавалась в условиях неопределенности относительно возможных действий противника, поэтому исследования проводились на адекватных математических моделях. В это время впервые был применен термин: «операционное исследование», подразумевающий исследования военной операции. В последующие годы операционные исследования или исследования операций развиваются как наука, результаты которой применяются для выбора оптимальных решений при управлении реальными процессами и системами. Математические модели, созданные для целей экономики, изучаются специальной научной дисциплиной, получившей название «экономико-математические методы». ЭММ представляют собой своеобразный инструментальный набор, с помощью которого экономисты, бизнесмены, менеджеры, стремясь добиться наилучшего эффекта «обрабатывают» свой материал. Научной основой ЭММ стали методы исследования операций. 1. Основные понятия и определения исследования операций 1.1. Цель и задачи исследования операций Исследование операций - научная дисциплина, занимающаяся разработкой и практическим применением методов наиболее эффективного управления различными организационными системами (система - это совокупность взаимосвязанных, взаимодействующих элементов (людей, машин ...), выполняющих определенную задачу). Управление любой системой реализуется как процесс, подчиняющийся определенным закономерностям. Их знание помогает определить условия, необходимые и достаточные для осуществления данного процесса. Для этого все параметры, характеризующие процесс и внешние условия, должны быть количественно определены, измерены. Т. о., цель исследования - количественное обоснование принимаемых решений по организации управления. При решении конкретной задачи управления применение методов ИО предполагает: - построение экономических и математических моделей для задач принятия решения в сложных ситуациях или в условиях неопределенности; - изучение взаимосвязей, определяющих впоследствии принятие решений, и установление критериев эффективности, позволяющих оценивать преимущество того или иного варианта действия. Операция - любое управляемое мероприятие, направленное на достижение цели. Результат операции зависит от способа ее проведения, организации, иначе - от выбора некоторых параметров. Всякий определенный выбор параметров называется решением. Оптимальными считают те решения, которые по тем или иным соображениям предпочтительнее других. Поэтому основной задачей исследования операций является предварительное количественное обоснование оптимальных решений. Для применения количественных методов исследования требуется построить математическую модель операции. Составление модели операции требует понимания сущности описываемого явления и знания математического аппарата. 1.2. Модели и моделирование. Одним из основных методов научного познания является эксперимент, а самой распространенной его разновидностью - метод моделирования систем. В процессе создания систем приходится проводить многочисленные исследования, эксперименты и расчеты, связанные с оценкой качества функционирования систем, с выбором лучшего варианта для ее создания. Выполнять их непосредственно на реальной системе очень сложно, иногда занимает много времени и экономически невыгодно. Существуют системы (экономика страны), на которых просто невозможно ставить эксперименты с познавательной целью. Значительно проще и дешевле создать модель системы и проводить на ней эксперименты. Под моделью принято понимать систему, способную замещать оригинал так, что ее изучение дает новую информацию об оригинале. Модель должна частично или полностью воспроизводить структуру моделируемой системы, ее функции. Под моделированием понимается процесс построения и исследования модели, способной заменить реальную систему и дать о ней новую информацию. Модели, используемые на практике, условно можно разделить на два типа: физические и символические. Символические модели описывают структуру и функции оригинала с помощью символов и соотношений между ними, выражающих определенные зависимости, присущие оригиналу. Большое место среди символических моделей занимают математические модели (уравнения, неравенства, функции, алгоритмы и т.д.), отражающие математические или логические зависимости. Математическая модель представляет собой систему математических и логических соотношений, описывающих структуру и функции реальной системы. Математическая модель отличается по своей физической природе от оригинала. Исследование свойств оригинала с помощью математической модели значительно удобнее, дешевле и занимает меньше времени по сравнению с физическим моделированием. Многие математические модели являются универсальными, т.е. могут использоваться для исследования различных систем. Целый ряд систем, в том числе экономических, либо трудно, либо вообще невозможно представить с помощью физических моделей. Существенную роль в развитии математического моделирования сыграли ЭВМ, способные выполнять различные по сложности вычисления и логические операции с большой скоростью. Среди математических моделей важное место занимают ЭММ, представляющие собой математическое описание экономических процессов и явлений. Большинство ЭММ включает в себя систему уравнений и неравенств, состоящих из набора переменных и параметров. Переменные величины характеризуют, например, объем производимой продукции, капитальных вложений, перевозок и т.п., а параметры - нормы расхода сырья, материалов, времени на производство определенной продукции. Практически в каждой модели можно выделить две группы переменных: 1) внешние переменные - их значения определяются вне данной модели и считаются в данной модели заданными; 2) внутренние переменные, значения которых определяются в результате исследования данной модели. ЭММ используются преимущественно для планирования или прогнозирования состояния системы на будущее. Наряду с использованием в предсказательных целях ЭММ применяются для описания реально существовавших или существующих экономических процессов. Выделяют описательные и оптимизационные ЭММ, которые используются на любых уровнях народнохозяйственной иерархии. Описательные модели экономических систем представляют собой формализованную с помощью математического аппарата экономическую задачу и используются для более глубокого изучения состояния системы и взаимосвязи ее элементов. К ним относятся матричные модели межотраслевых балансов народного хозяйства и экономического района, производственные функции и др. При определении исходных данных задачи модели данного типа позволяют получить единственное решение. Основной недостаток этих моделей - отсутствие условия нахождения оптимального решения. Оптимизационные модели отражают в математической форме смысл экономической задачи. Отличительной особенностью этих моделей является наличие условия нахождения оптимального решения (критерия оптимальности), которое записывается в виде функционала. Эти модели при определенных исходных данных задачи позволяют получать множество решений, удовлетворяющих условиям задачи, и обеспечивают выбор оптимального решения, отвечающего критерию оптимальности: модели определения оптимальной производственной программы, модели оптимального смешивания компонентов, оптимального раскроя, оптимального размещения предприятий некоторой отрасли на определенной территории, модели транспортной задачи. Большинство существующих оптимизационных моделей являются моделями планирования и имеют один критерий оптимальности. 1.3. Процесс экономико-математического моделирования. Этот процесс состоит из нескольких взаимосвязанных этапов. Разбиение на этапы и выделение на каждом этапе присущих ему процессов условно: на одном из выделенных этапов возможно совмещение процессов, относящихся к разным этапам. Первый этап - постановка задачи. Данный этап начинается с выработки цели исследования. Для конкретной экономической системы цели исследования могут быть различными, например, для предприятия можно задаться целью составить оптимальный план выпуска продукции или перевозок грузов, либо найти оптимальный вариант раскроя исходных материалов и т.д. Исходя из цели исследования необходимо провести подробный анализ системы, выявить ее структуру и функции, изучить особенности. В процессе постановки задачи необходимо помнить, что модель должна, во-первых, правильно воспроизводить действительность, во-вторых, быть доступной для исследования. Эти два обстоятельства оказывают существенное влияние на выбор исходных предпосылок. При моделировании экономических систем, исходя из цели исследования, с одной стороны, необходимо выбрать самые важные в условиях данной задачи факторы и ввести в модель только те, которые самым существенным образом влияют на результат решения, на достижение поставленной цели. Учет в модели несущественных факторов приводит к тому, что модель становится сложной для понимания моделируемой системы и для решения. С другой стороны, игнорирование многих факторов может привести к чрезмерному упрощению модели, нарушению соответствия ее действительности. Компромисс между этими двумя требованиями достигается методом проб и ошибок. Эйнштейн утверждал, что правильная постановка задачи более важна, чем ее решение. Второй этап - построение математической модели На этом этапе проводится формализация задачи - построение математических зависимостей в виде уравнений, неравенств, функций и т.п. Формализованную с помощью математического аппарата запись экономической задачи называют моделью задачи. Приступая к формализации экономического процесса, необходимо проанализировать, подходит ли для его описания одна из ранее созданных ЭММ. К настоящему моменту создано несколько десятков так называемых универсальных, или типовых, моделей (модель транспортной задачи, модели задачи о ранце, диете, раскрое и т.п.), которые используются на практике для описания различных экономических процессов. Самой универсальной моделью считается модель транспортной задачи, с помощью которой формализуется не только процесс перевозки грузов, но и процесс размещения предприятий отрасли на определенной территории, процесс назначения работников на работы и др. Третий этап - получение решения с помощью построенной модели. Основные задачи данного этапа. Первая задача - сбор и обработка необходимой для модели достоверной исходной информации, определение числовых значений параметров и внешних переменных. На практике не всегда удается собрать требуемую информацию, что приводит к невозможности использования модели в полученном виде. Тогда приходится возвращаться к постановке задачи и приспосабливать ее к имеющимся исходным данным. Вторая задача - выбор метода получения решения: используются аналитические (формульные) и численные экономико-математические методы: симплекс-метод, метод потенциалов и др. Экономико-математические методы в определенной степени универсальны и используются для решения различных экономических задач. Однако не любая задача укладывается в рамки модели, для которой уже разработаны эффективные аналитические или численные методы решения. В этом случае пользуются другими методами получения решения, в частности эвристическими и имитационными методами исследования систем. Эвристика (в переводе с греческого - нахожу, придумываю, открываю) - это совокупность неформальных методов решения задач (эвристических методов), основанных на прошлом опыте, интуиции решающего. Эвристические методы в общем случае не гарантируют получение наилучшего решения, поскольку они опираются не на доказательства, а на так называемые правдоподобные рассуждения. Имитационное моделирование следует рассматривать как новую методологию, новое направление в моделировании, позволяющее расширить его возможности. Под имитационным моделированием понимается экспериментирование с моделью реальной системы, в частности, вычислительный эксперимент, проводимый с помощью математической модели путем изменения различных исходных предпосылок. Поскольку вручную такие эксперименты просто невозможны, ИМ получило развитие только с появлением ЭВМ. Имитация (в переводе с латинского - подражание) - это воспроизведение чего-либо искусственными средствами, что позволяет постичь суть явления, не прибегая к экспериментам на реальном объекте. Имитационные модели служат для анализа поведения системы в условиях, определяемых экспериментатором. Четвертый этап - применение полученных с помощью модели результатов на практике. Сложность экономических процессов и явлений, другие особенности экономических систем затрудняют не только построение моделей, но и проверку их адекватности - соответствия ЭММ рассматриваемой экономической системе, цели ее исследования. Любая модель любой системы предполагает абстрагирование от некоторых реальных свойств объекта и отражает лишь основные его свойства. На данном этапе проверяется, насколько принятые допущения правомерны и, следовательно, применима ли построенная модель для исследования моделируемой системы. В случае необходимости модель корректируется. С целью обоснования пригодности модели для конкретных исследований проводится так называемый анализ модели на чувствительность. Полученное с помощью ЭММ решение анализируется на чувствительность путем изменения исходной информации в определенных пределах. Важность данной задачи состоит в том, что исходная информация со временем может меняться и необходимо знать, как будут влиять эти изменения на получаемое решение. 1.4. Общая постановка задачи исследования операций. Все факторы, входящие в описание операции, можно разделить на две группы: - постоянные факторы
(условия проведения операции), на которые мы влиять не можем. Обозначим их через - зависимые факторы (элементы решения) x1 ,x2 ,..., которые в известных пределах мы можем выбирать по своему усмотрению. Критерий эффективности, выражаемый некоторой функцией, называемой целевой , зависит от факторов обеих групп, поэтому целевую функцию Z можно записать в виде Z = f(x1
,x2,...,
Все модели исследования операций могут быть классифицированы в зависимости от природы и свойств операции, характера решаемых задач, особенностей применяемых математических методов. Следует отметить прежде всего большой класс оптимизационных моделей. Такие задачи возникают при попытке оптимизировать планирование и управление сложными системами, в первую очередь экономическими системами. Оптимизационную задачу можно сформулировать в общем виде: найти переменные x1 , x2 ,..., xn , удовлетворяющие системе неравенств (уравнений) и обращающие в максимум (или минимум) целевую функцию, т.е. Z =
f(
(Условия неотрицательности переменных, если они есть, входят в ограничения (1.1)). В тех случаях, когда функции f и Математическое программирование - область математики, разрабатывающая теорию и численные методы решения многомерных экстремальных задач с ограничениями, т.е. задач на экстремум функции многих переменных с ограничениями на область изменения этих переменных. Функцию, экстремальное значение которой нужно найти в условиях экономических возможностей, называют целевой, показателем эффективности или критерием оптимальности. Экономические возможности формализуются в виде системы ограничений. Все это составляет математическую модель. Математическая модель задачи - это отражение оригинала в виде функций, уравнений, неравенств, цифр и т.д. Если критерий эффективности Z = f(x1
,x2
,...,xn
) представляет линейную функцию, а функции Если в задаче математического программирования имеется переменная времени и критерий эффективности (1.2) выражается не в явном виде как функция переменных, а косвенно - через уравнения, описывающие протекание операций во времени, то такая задача является задачей динамического программирования. Из перечисленных методов математического программирования наиболее распространенном и разработанным является линейное программирование. В его рамки укладывается широкий круг задач исследования операций. В создание современного математического аппарата и развитие многих направлений исследования операций большой вклад внесли российские ученые Л.В. Канторович, Н.П. Бусленко, Е.С. Вентцель, Н.Н. Воробьев, Н.Н. Моисеев, Д.Б. Юдин и многие другие. Значительный вклад в формирование и развитие исследования операций внесли зарубежные ученые Р. Акоф, Р. Беллман, Г. Данциг, Г. Кун, Дж. Нейман, Т. Саати, Р. Черчмен, А. Кофман и др. Методы исследования операций, как и любые мат. методы, всегда в той или иной мере упрощают, огрубляют задачу, отражая порой нелинейные процессы линейными моделями, динамические процессы - статическими моделями и т.д. Жизнь богаче любой схемы. Поэтому не следует ни преувеличивать значение количественных методов исследования операций, ни преуменьшать его, ссылаясь на примеры неудачных решений. Уместно привести в связи с этим шутливо-парадоксальное определение исследования операций, сделанное одним из его создателей Т.Саати, как «искусства давать плохие ответы на те практические вопросы, на которые даются еще худшие ответы другими методами». 2. Линейное программирование Термин «линейное программирование» впервые появился в 1951 г. в работах американских ученых (Дж. Данциг, Т. Купманс), а первые исследования по линейному программированию (основные задачи и приложения, критерий оптимальности, экономическая интерпретация, методы решения, геометрическая интерпретация результатов решения) были проведены в конце 30-х годов в СССР в Ленинградском университете Л. В. Канторовичем. Под линейным программированием понимается линейное планирование, т.е. получение оптимального плана - решения в задачах с линейной структурой. Линейное программирование широко применяется в сфере военной деятельности, сельском хозяйстве, промышленности, управлении производственными процессами и запасами, в экономике и на транспорте. 2.1. Общая задача линейного программирования. Общей задачей линейного программирования (ОЗЛП) называют задачу: Максимизировать (минимизировать) функцию при ограничениях: где cj
,
aij
, bi
-заданные действительные числа, (2.1) - целевая функция, (2.2) - ограничения, Экономическая интерпретация модели ЛП состоит в следующем. Моделируемая система характеризуется наличием нескольких видов «производственной деятельности» Цель построения модели состоит в определении уровней (объемов производства) каждого вида производственной деятельности xj , при которых оптимизируется (максимизируется или минимизируется) общий результат производственной деятельности системы в целом без нарушения ограничений, накладываемых на использование ресурсов. Оптимальным решением
(или оптимальным планом
) задачи ЛП называется решение Термины «решение» и «план» - синонимы, однако первый используется чаще, когда речь идет о формальной стороне задачи (ее математическом решении), а второй - о содержательной стороне (экономической интерпретации). Симметричной формой записи ЗЛП называют задачу 2.2. Линейное программирование в экономике 1. Задача о наилучшем использовании ресурсов.
Пусть некоторая производственная единица (цех, завод, фирма и т.д.), исходя из конъюнктуры рынка, технических или технологических возможностей и имеющихся ресурсов, может выпускать n различных видов продукции (товаров) Пj
, Математическая модель задачи имеет следующий вид: Так как переменные xj
входят в целевую функцию Z( 2. Задача о смесях. В различных отраслях народного хозяйства возникает проблема составления таких рабочих смесей на основе исходных материалов, которые обеспечивали бы получение конечного продукта, обладающего определенными свойствами. К этой группе задач относятся задачи формировании минимальной потребительской продовольственной корзины, составлении кормового рациона в животноводстве, шихт в металлургии, горючих и смазочных смесей в нефтеперерабатывающей промышленности, смесей для получения бетона в строительстве и т.д. Высокий уровень затрат на исходные сырьевые материалы и необходимость повышения эффективности производства выдвигает на первый план следующую задачу: получить продукцию с заданными свойствами при наименьших затратах на исходные сырьевые материалы. Модель задачи о наилучшем составе смеси рассмотрим на примере задачи формирования минимальной потребительской продовольственной корзины. Задан ассортимент продуктов Решение задачи - это количества xj продуктов каждого вида, обеспечивающие необходимое количество питательных веществ при минимальных затратах на исходные продукты. Математическая модель задачи имеет следующий вид: 3. Задача о раскрое материалов. Сущность задачи об оптимальном раскрое состоит в разработке таких технологически допустимых планов раскроя, при которых получается необходимый комплект заготовок, а отходы ( по длине, площади, объему, массе или стоимости) сводятся к минимуму. На раскрой (распил, обработку) поступает материал нескольких видов в определенном количестве. Из этого материала необходимо изготовить различные изделия. Материал может быть раскроен разными способами. Каждый способ имеет свою себестоимость и позволяет получить разное количество изделий каждого вида. Определить способ раскроя, при котором суммарная себестоимость минимальна. Пусть n - число различных видов материала, поступающего на раскрой; dj
- количество материала j-го вида, Обозначим через xjk
- количество единиц материала j-го вида, раскраиваемых k-м способом, Математическая модель задачи имеет следующий вид:
Вместо критерия минимизации себестоимости в задаче может быть взят, например, критерий минимизации отходов. В этом случае в условии должно быть задано количество отходов, получаемых при каждом способе раскроя для единицы материала каждого вида. 2.3. Геометрическая интерпретация и графическое решение ЗЛП . Геометрическая интерпретация экономических задач дает возможность наглядно представить их структуру, выявить особенности и открывает пути исследования более сложных свойств. ЗЛП с двумя переменными всегда можно решить графически. Случай двух переменных не имеет особого практического значения, но его рассмотрение проясняет свойства ОЗЛП, приводит к идее ее решения, делает геометрически наглядными способы решения и пути их практической реализации. Пусть дана задача: Z=c1
x1
+c2
x2
Дадим геометрическую интерпретацию элементов этой задачи. Введем на плоскости декартову прямоугольную систему координат x1 Ox2 и сопоставим каждой паре чисел (x1 ,x2 ) точку плоскости с координатами x1 и х2. Выясним сначала, что представляет собой множество точек, соответствующих допустимым решениям данной задачи. Рассмотрим одно линейное неравенство Оно определяет на плоскости одну из двух полуплоскостей, на которые прямая Убедиться в том, с какой стороны от прямой лежит полуплоскость, точки которой удовлетворяют заданному неравенству, можно путем подстановки координат точек одной или другой полуплоскости в неравенство. Если координаты точки удовлетворяют неравенству, то эта точка лежит в полуплоскости, соответствующей данному неравенству. В противном случае неравенству соответствует другая плоскость. Каждое из ограничений (2.7), (2.8) задает на плоскости х1 Oх2 некоторую полуплоскость. Нас интересуют те точки плоскости, координаты которых принадлежат всем полуплоскостям. Следовательно, допустимое множество ЗЛП геометрически изображается пересечением (общей частью) полуплоскостей, определяемых отдельными ограничениями. Полуплоскость - выпуклое множество. Множество называется выпуклым, если ему вместе с двумя произвольными точками принадлежит и прямолинейный отрезок, их соединяющий. Пересечение любого числа выпуклых множеств является выпуклым множеством. Таким образом, область допустимых решений задачи (2.6) - (2.8) есть выпуклое множество. На рис.2.1 представлены возможные ситуации, когда область допустимых решений ЗЛП - выпуклый многоугольник (а), неограниченная выпуклая многоугольная область (б), единственная точка (в), прямая линия (г), луч (д), отрезок (е), пустое множество (ж).
0 а) X1 0 X1
![]()
г) 0 Х1 0 Х1 0 Х 1 Х2 Х2
Х1
Рис.2.1 Перейдем к геометрической интерпретации целевой функции. Пусть область допустимых решений ЗЛП - непустое множество, например, многоугольник ABCDE рис.2.2.
E
0 X1 Рис.2.2 N Выберем произвольное значение целевой функции Z=Z0 . Получим c1 x1 +c2 x2 =Z0 Это уравнение прямой линии (рис.2.2 - прямая MN). В точках прямой MN целевая функция сохраняет одно и то же постоянное значение Z0 . Считая Z0 параметром, получим семейство параллельных прямых, называемых линиями уровня целевой функции (линиями постоянного значения). Нас интересуют те точки области допустимых решений, которые принадлежат линии уровня с наибольшим значением Z0 по сравнению с его значениями для всех других линий уровня, пересекающихся с областью допустимых решений. Возникает вопрос: как установить направление возрастания (убывания) целевой функции по x1 и х2 : Частная производная (2.9) и (2.10) функции показывает скорость ее возрастания вдоль данной оси. Следовательно, с1
и с2
- скорости возрастания вдоль осей Oх1
и Oх2.
Вектор Из геометрической интерпретации элементов ЗЛП вытекает следующий порядок ее графического решения. 1, С учетом системы ограничений строим область допустимых решений. 2. Строим вектор 3. Проводим произвольную линию уровня Z=Z0,
перпендикулярную к вектору 4. При решении задач на max перемещаем линю уровня Z=Z0
в направлении вектора Определяем оптимальный план Пример 2. 1. Предприятию необходимо изготовить два вида продукции А и В, с использованием трех видов ресурсов R1, R2 , R3 количество которых ограничено. Исходные данные задачи представлены в таблице:
Требуется составить такой план выпуска продукции, чтобы при ее реализации получить максимальный доход. Решение. Обозначим через х1 и х2 количество единиц продукции видов А и В, планируемых к выпуску. Известно, что доход от реализации единицы продукции А составляет 12 усл. ед. и количество этой продукции - х1 . Следовательно, доход от реализации всей продукции А составит 12х1 усл. ед. Аналогично, доход от реализации всей продукции В составит 15х2 усл. ед. Учитывая, что доход от реализации продукции должен быть максимальным, целевая функция задачи будет иметь вид: Известно также, что имеющиеся на предприятии ресурсы ограничены. Это обстоятельство в свою очередь необходимо отразить в модели. Предприятие производит продукцию, используя три вида ресурсов. Естественно, что фактический расход никакого вида ресурса не должен превышать запасов соответствующего вида ресурсов на предприятии. Поскольку расход каждого вида ресурсов на единицу каждого вида продукции и запасы ресурсов известны, это обстоятельство отражается в следующих ограничениях: Смысл первого ограничения состоит в том, что фактический расход ресурса R1 на производство продукции А и В, равный 6х1 +6х2 (здесь 6х1 - количество единиц ресурса R1 , идущего на изготовление х1 единиц продукции A; 6х2 - количество единиц ресурса R1, идущее на изготовление х2 единиц продукции В) не должен превышать запаса этого ресурса на предприятии, равного 36 ед. Аналогичный смысл имеют 2-е и 3-е ограничения только для ресурсов R2 и R3 соответственно. Количество продукции, выпускаемое предприятием, должно быть величиной положительной или равной нулю (если предприятие определенный вид продукции не производит). Следовательно, в модели должны присутствовать ограничения неотрицательности переменных: Таким образом, построена математическая модель нашей задачи как задачи линейного программирования: Начнем решение задачи с построения области допустимых решений (рис.2.3) В первую очередь отобразим в прямоугольной системе координат условия неотрицательности переменных. В двумерном пространстве уравнению соответствует прямая линия, а неравенству - полуплоскость, лежащая по одну сторону от прямой. Прямые х1
=0 и х2
=0 совпадают с осями координат. Полуплоскости x1
>0,x2
>0 лежат соответственно справа от оси Oх2
и выше оси Oх1
. Множество точек, удовлетворяющих одновременно неравенствам Теперь рассмотрим ограничения задачи. Построим по порядку прямые: ` 6x1 +6x2 =36 или х1 +х2 =6 (а) 4х1 +2х2 =20 или 2х2 +х2 =10 (б) 4х1 +8х2 =40 или х1 +2х2 =10 (в) и определяем, с какой стороны от этих прямых лежат полуплоскости, точки которых удовлетворяют строгим неравенствам: 6x1 +6x2 <36 4x1 +2x2 <20 4x1 +8x2 <40 Для определения полуплоскости решений первого неравенства возьмем произвольную точку плоскости, не лежащую на прямой 6х1
+6х2
=36, например (0;0), и подставим ее координаты в неравенство В результате подстановки получили верное числовое неравенство 0< 36, и это означает, что начало координат лежит в полуплоскости решений первого неравенства. Сторона, в которой располагается полуплоскость от прямой, указывается штриховой. Аналогично строим полуплоскость решений остальных неравенств.
X2
Заштрихованная часть плоскости и представляет собой искомый многоугольник допустимых решений задачи (рис.2.3) Теперь нужно среди точек построенного многоугольника найти такую, в которой целевая функция Z=12x1 +15x2 достигает максимального значения. Для этого построим прямую, заданную уравнением 12х1 +15х2 =0, которая является линией нулевого уровня функции Z. Направление возрастания линейной функции Z=12x1
+15x2
указывает вектор Для нахождения оптимального плана нужно «передвигать» линию нулевого уровня Z параллельно самой себе в направлении вектора Найдем соответствующее значение целевой функции: Ответ . Для обеспечения максимального дохода от реализации готовой продукции предприятию необходимо выпустить 2 ед. продукции вида А и 4 ед. продукции вида В. При таком плане доход от реализации составит 84 усл. ед. Геометрический метод решения ЗЛП обладает рядом достоинств: он прост, нагляден, позволяет быстро и легко получить ответ. Однако только геометрический метод решения никак не может удовлетворить ни математиков, ни экономистов. Возможны «технические» погрешности, которые неизбежно возникают при приближенном построении графиков. Второй недостаток геометрического метода заключается в том, что многие величины, имеющие четкий экономический смысл (такие, как остатки ресурсов производства, избыток питательных веществ и т.п.), не выявляются при геометрическом решении задач. Но самое главное - геометрический метод неприемлем для решения практических задач. Поэтому необходимы аналитические методы, позволяющие решать ЗЛП с любым числом переменных и выявить экономический смысл входящих в них величин. 3. Транспортная задача Методы линейного программирования, являются хорошим инструментом для решения ряда проблем распределения ресурсов. Применение пакетов прикладных программ позволяет значительно упростить решение задачи. Поэтому лицо, принимающее решение, получает возможность уделить большее внимание интерпретации и оценке решения задачи. Однако применение ППП предполагает предварительную формализацию модели линейного программирования. В процессе решения большинства проблем эта задача является основной. При построении модели необходимо идентифицировать ее переменные и сформулировать систему ограничений. При решении некоторых видов проблем распределения ресурсов использование специально созданных для этих целей алгоритмов упрощает процесс построения исходной модели. В дальнейшем рассмотрим два примера таких алгоритмов, созданных для решения транспортной задачи и задачи о назначениях. Эти задачи используются для моделирования и оптимизации экономических проблем, связанных с формированием оптимального плана перевозок, оптимального распределения индивидуальных контрактов на транспортировки, составление оптимального штатного расписания, определения оптимальной специализации предприятий, рабочих участков, оптимального назначения кандидатов на работы, оптимального использования торговых агентов. 3.1. Постановка задачи Пусть имеется m поставщиков А1
, А2
, ...,Аm
однородного груза в количествах соответственно а1,
а2
,...,аm
единиц и n потребителей В1,
В2
,...,Вn
этого груза, потребность которых составляет соответственно b1
, b2,
...,bn
единиц. Известна стоимость (тариф) cij
перевозки единицы груза от i-го (i= Обычно исходные данные транспортной задачи представляют в виде табл. 1: табл.1
При постановке конкретных задач перевозки грузов может возникнуть одна из трех ситуаций: Каждой ситуации соответствует определенная модель ТЗ: ситуации (3.1) (суммарный объем груза, имеющегося у поставщиков, равен суммарному спросу потребителей) отвечает закрытая модель ТЗ (сбалансированная транспортная модель), а ситуациям (3.2) и (3.3) - отвечает открытая модель ТЗ (несбалансированная транспортная модель). 3.2. Экономико-математическая модель ТЗ Рассмотрим ситуацию (3.1). Обозначим через Экономико-математическая модель ТЗ должна отражать все условия и цель задачи в математической форме. Переменные
Цель ТЗ - минимизировать общие затраты на реализацию плана перевозок, которые можно представить функцией:
Итак, математически ТЗ ставится так. Даны система ограничений (3.4) и (3.5) (ограничения (3.5) отражают тот факт, что весь груз от всех поставщиков должен быть вывезен, а ограничения (3.4) отражают тот факт, что каждый потребитель должен получить ровно столько груза, сколько ему необходимо) при условии (3.6) и линейная функция (3.7). Найти такое неотрицательное решение, при котором линейная функция (4.7) принимает минимальное значение. Для того, чтобы ТЗ (3.4)-(3.7) имела допустимые планы, необходимо и достаточно выполнение равенства (3.1) Решение транспортной задачи состоит из двух этапов: 1 этап. Нахождение начального плана перевозок (xij
), 2 этап. Улучшение начального плана перевозок и получение оптимального плана перевозок (xij
), 3.3. Построение исходного опорного плана. Построение опорных планов, а также их преобразование будем производить непосредственно в распределительной таблице (табл.1). Если в плане перевозок переменная Метод «северо-западного угла». Определение значений xij начинается с левой верхней, условно называемой северо-западной, клетки (1,1) табл.1. Находим x11 =min(a1 , b1 ). 1) если а1 <b1 , то x11 =a1 , строка 1 исключается из дальнейшего рассмотрения, а потребность первого потребителя b1 уменьшается на а1 ; 2) если а1 >b1 , то х11 =b1 , столбец 1 исключается из дальнейшего рассмотрения (первый потребитель В1 будет полностью удовлетворен), а наличие груза у первого поставщика a1 уменьшается на b1 ; 3) если a1 =b1 , то x11 = a1 =b1 , первая строка и первый столбец исключаются из дальнейшего рассмотрения. Эта ситуация приводит к вырождению исходного решения. Затем аналогичные операции проделывают с оставшейся частью таблицы, начиная с ее северо-западного угла. На последнем шаге процесса остается одна строка (последняя) и один столбец (последний). После заполнения клетки, стоящей на их пересечении, т.е. клетки (m, n), процесс завершается. После завершения описанного процесса необходимо провести проверку полученного плана (решения) на вырожденность. Если количество заполненных (занятых) клеток равно m+n-1, то план является невырожденным, в противном случае - вырожденным. Если план вырожденный, то незаполненные клетки с минимальными стоимостями перевозок заполняются нулями, чтобы общее число заполненных клеток стало равным m+n-1. Однако, при расстановке нулей необходимо помнить, что в таблице не должно быть ни одного прямоугольника, все вершины которого являются заполненными клетками. Например, переменные x11 ,x12 ,x21 ,x22 не могут быть одновременно базисными. Метод минимального элемента. В отличии от метода северо-западного угла данный метод учитывает при построении исходного плана стоимости перевозок. В ряде случаев он позволяет получить лучшее с точки зрения критерия оптимальности решение, сокращая количество итераций для получения оптимального плана. Определение значений xij начиная с клетки, имеющей минимальную стоимость перевозки. Если в таблице имеется несколько клеток с одинаковыми минимальными стоимостями, то заполняется прежде та клетка, в которую можно вписать большую поставку. Переменной, отвечающей выбранной клетке, присваивается минимальное из двух возможных значений: xij =min(ai ,bj ). Соответствующая строка или столбец исключаются из дальнейшего рассмотрения, а потребность потребителя или наличие груза у поставщика уменьшается на выбранную величину. Если для выбранной клетки с минимальной стоимостью перевозки наличие груза у поставщика равно потребности потребителя, то из дальнейшего рассмотрения исключаются и строка и столбец (это приводит к вырождению исходного плана). Затем в оставшейся части таблицы проделывают аналогичные операции, опять начиная с клетки, имеющей минимальную стоимость перевозки. Проверка плана на вырожденность и расстановка ( в случае вырожденности плана) нулей осуществляется так же, как это описано для метода северо-западного угла. Для нахождения оптимального плана перевозок необходимо уметь оценивать полученный план на оптимальность. Как это сделать, не имея в распоряжении всех возможных планов перевозок, которые можно было бы сравнить между собой? Для оценки плана на оптимальность вводится понятие косвенных затрат. Косвенные затраты - это затраты, получаемые для маршрутов, по которым не осуществляются перевозки при данном плане. Рассчитанные косвенные затраты сравниваются с реальными затратами, которые имели бы место, если бы перевозки по данным маршрутам осуществлялись. Если для всех невыбранных маршрутов косвенные затраты не меньше реальных, то данный план перевозок является оптимальным. Если хотя бы для одного маршрута косвенные затраты меньше реальных, то план перевозок может быть улучшен путем введения в него данного маршрута. Ввод нового маршрута в план перевозок соответствует вводу в список базисных переменных переменной транспортной задачи, соответствующей этому маршруту. Эти рассуждения лежат в основе ряда методов, применяемых для нахождения оптимального плана перевозок. Рассмотрим один из них - метод потенциалов. 3.4 . Алгоритм решения ТЗ методом потенциалов 1. Построить опорный план по одному из правил. 2. Проверить план на невырожденность. Если полученный план вырожденный, формально заполняют нулями некоторые из свободных клеток так, чтобы общее число занятых клеток было равно m+n-1. Нули надо расставлять так, чтобы не образовался замкнутый цикл из занятых клеток. 3. Проверка плана на оптимальность. 3.1. Определение потенциалов поставщиков и потребителей. Для всех занятых клеток рассчитываются потенциалы по формуле где ui - потенциалы поставщиков (строк), vj - потенциалы потребителей (столбцов). Так как всех занятых клеток должно быть m+n-1, т.е. на единицу меньше числа потенциалов (их m+n), то для определения чисел ui , vj необходимо решить систему из m+n-1 уравнений ui +vj =cij c m+n неизвестными. Система неопределенная, и, чтобы найти частные решения, одному из потенциалов придают произвольное числовое значение (обычно полагают u1 =0), тогда остальные потенциалы определяются однозначно. 3.2. Для всех свободных клеток рассчитываются оценки по формуле Здесь Если все s 4. Среди отрицательных оценок выбирается минимальная. Например, для клеток (i, j) и (k,l) имеем оценки: sij
=-2, skl
=-4. В этом случае наиболее потенциальной (перспективной) является клетка (k,l). Экономически оценка показывает, на сколько денежных единиц уменьшатся транспортные расходы от загрузки данной клетки единицей груза. Так, от загрузки клетки (i,j) единицей груза транспортные расходы уменьшаются на Для наиболее перспективной свободной клетки (клетки с min отрицательной оценкой) строится замкнутый цикл. Для этого из выбранной свободной клетки проводится по строке или столбцу прямая линия до любой занятой клетки, затем под углом 90 5. В вершинах цепи, чередуя проставляются знаки «+» и «-», начиная с выбранной свободной клетки (в нее помещают знак «+»). В клетках со знаком «-» выбирается минимальная поставка (
Рис.3.1 Может случится, что найденное наименьшее число (
30 10
![]()
20 40 Рис.3.2 Составляется новый план, рассчитывается значение целевой функции. Переходим к шагу 2. 3.5. Пример решения транспортной задачи На четырех строительных площадках В1 , В2 , В3 , В4 монтируется в день соответственно 20,120,20 60 м3 сборных плит перекрытий. Производство этих плит налажено на трех заводах А1, А2 , А3 в размере соответственно 100,70 и 50 м3 . Известны стоимости перевозки (табл.2) 1м3 сборных плит из каждого пункта производства в каждый пункт потребления (ден. ед./ м3 ). Требуется так закрепить строительные площадки за заводами, чтобы при полном обеспечении сборными плитами перекрытий затраты на перевозку были минимальными. РЕШЕНИЕ. Исходное опорное решение получим по методу «северо-западного угла» (табл.3) и по методу min элемента (табл.4). Здесь табл.2
табл.3
Транспортные расходы f= табл.4
f= Транспортные расходы для опорного плана, построенного по методу min элемента, меньше. Поэтому за исходное решение возьмем то, которое получено по методу min элемента (табл.4). Количество заполненных клеток в табл.4 равно 6: m+n-1=3+4-1=6. Следовательно, полученный план невырожденный. Для определения потенциалов (см. формула 3.8) составляем уравнения: u1 +v2 =4, u2 +v1 =1, u2 +v4 =2, u3 +v2 =4, u3 +v3 =1, u3 +v4 =4. Положим u1 =0, тогда v2 =4, u3 =0, v3 =1, v4 =4, u2 =-2, v1 =3. Потенциалы проставлены в табл.5 (последняя строка и последний столбец). Их можно вычислять и непосредственно в таблице не выписывая систему уравнений. Т.к. если известны потенциал и тариф (стоимость перевозки) занятой клетки, то из соотношения ui +vj =cij легко определить неизвестный потенциал (из суммы вычесть известное слагаемое, получим неизвестное слагаемое. Роль суммы в данном равенстве играет тариф cij ) Определим по формуле (3.9) оценки свободных клеток: s11 =6-(0+3)=3>0, s13 =2-(0+1)=1>0, s14 =4-(0+4)=0 s22 =2-(-2+4)=0, s23 =7-(-2+1)=8>0, s31 =2-(0+3)= -1<0
План неоптимальный, т.к. s31
<0. Строим для клетки (3;1) цикл непосредственно в табл.5. В цикл войдут клетки (3;1), (3;4), (2;4), (2;1). Наименьшее количество груза, стоящее в вершинах цикла с отрицательным знаком, Будет ли полученный план оптимальным? План невырожденный. Определим для него новые потенциалы: u1 +v2 =4, u2 +v1 =1, u2 +v4 = 2, u3 +v1 =2, u3 +v2 =4, u3 +v3 =1. Положим u1 =0, тогда v2 =4, u3 =0, v1 =2, v3 =1, u2 = -1, v4 =3. Проставим значения найденных потенциалов в табл.6. Определим оценки свободных клеток: s11 =6-(0+2)=4>0, s13 =2-(0+1)=1>0, s14 =4-(0+3)=1>0 s22 =2-(-1+4)=-1<0, s23 =7-(-1+1)=7>0, s44 =4-(0+3)=1>0. План (табл.6) неоптимальный, т.к. s22
<0. Строим для клетки (2;2) цикл непосредственно в табл.6. В цикл войдут клетки (2;2), (2;1), (3;1), (3;2). Наименьшее количество груза, стоящее в вершинах цикла с отрицательным знаком, f =
Будет ли полученный план оптимальным? Определим для него новые потенциалы: u1 +v2 =4, u2 +v2 =2, u2 +v4 =2, u3 +v1 =2, u3 +v2 =4, u3 +v3 =1. Положим u1 =0, тогда v2 =4, u2 =-2, v4 =4, u3 =0, v1 = 2, v3 =1. Проставим значения найденных потенциалов в табл.7. Определим оценки свободных клеток: s11 =6-(0+2)=4>0, s13 =2-(0+1)=1>0, s14 =4-(0+4)=0 s21 =1-(-2+2)=1>0, s23 =7-(-2+1)=8>0, s34 =4-(0+4)=0 Оценки свободных клеток неотрицательны, следовательно, полученный план является оптимальным: x12 =100, x22 =10, x24 =60, x31 =20, x32 =10, x33 =20 Минимальные транспортные расходы для этого плана f=640. Все оценки свободных клеток равны нулю. Это свидетельствует о неединственности оптимального плана. ОТВЕТ. Согласно оптимальному плану, с первого завода A1 нужно поставить 100м3 перекрытый на вторую строительную площадку B2 , с завода А2 - 10м3 на строительную площадку В2 и 60м3 на строительную площадку В4 , с завода А3 - на 20 м3 на строительную площадку В1 , 10 м3 на строительную площадку В2 и 20 м3 на строительную площадку В3. 3.7. Транспортные задачи, имеющие некоторые усложнения в постановке 1. Транспортная задача с избытком запасов : Для отыскания оптимального плана вводят фиктивный (n+1)-й пункт назначения Bn
+1
с потребностью bn
+1
и полагают стоимости перевозок грузов в Bn
+1
равными нулю. Полученная новая транспортная задача является закрытой. Найдя оптимальный план xij
2. Транспортная задача с избытком заявок : Эта задача несколько сложнее предыдущей: дело в том, что как ни развози грузы, все потребности удовлетворить все равно не удается. Поэтому постановка задачи нуждается в уточнении. 2.1. Все пункты назначения требуется удовлетворить пропорционально поданным заявкам. В этом случае, подсчитав коэффициент пропорциональности и «подправив» заявки: 2.2. Если не заботиться о «справедливости» удовлетворения заявок, а по-прежнему интересоваться лишь минимизацией транспортных расходов, то для отыскания оптимального плана вводят фиктивный пункт (m+1)-й отправления Am
+1
с запасом груза am
+1
и полагают стоимости перевозок грузов из Am
+1
равными нулю. Полученная транспортная задача является закрытой. Найдя оптимальный план xij
И в случае 1 и в случае 2.2 при преобразовании открытой задачи в закрытую, целевая функция не меняется, так как все слагаемые, соответствующие дополнительным перевозкам, равны нулю. 3. Транспортная задача с поиском максимума целевой функции . Перейдем к целевой функции Z=-F и найдем план ( 4. Транспортная задача с запретными маршрутами. Речь идет о задачах, в которых нельзя перевозить груз из некоторых пунктов отправления Ai в некоторые пункты назначения Bj . В этом случае стоимости соответствующих перевозок полагаем равными достаточно большому числу. Тогда при отыскании оптимального плана соответствующие перевозки будут блокированы. 5. Транспортная задача с обязательными поставками. Иногда приходится решать транспортную задачу, в которой дополнительным условием в ограничениях является обязательное обеспечение конкретных перевозок по определенным маршрутам. В этом случае каждую обязательную перевозку xij =dij реализуем условно, уменьшая на dij запасы в Ai и потребности в Bj . Если это не удается сделать, то исходная задача не имеет решения. В противном случае стоимости обязательных поставок полагаем равными достаточно большому числу, решаем полученную задачу и от ее оптимального плана переходим к оптимальному плану исходной задачи. 6. Транспортная задача с ограничениями снизу. Пусть требуется решить транспортную задачу, в которой некоторые из перевозок ограничены снизу xij 3.8. Экономические задачи, сводящиеся к транспортным моделям Алгоритмы и методы решения транспортной задачи могут быть использованы при решении некоторых экономических задач, не имеющих ничего общего с транспортировкой груза. В этом случае величины тарифов cij имеют различный смысл в зависимости от конкретной экономической задачи. К таким задачам относятся следующие: - оптимальное закрепление за станками операций по обработке деталей. В них cij является таким экономическим показателем, как производительность. Задача позволяет определить, сколько времени и на какой операции нужно использовать каждый из станков, чтобы обработать максимальное количество деталей. Так как транспортная задача требует нахождения минимума, то значения cij берутся с отрицательным знаком; - оптимальные назначения, или проблема выбора. Имеется m механизмов, которые могут выполнять m различных работ с производительностью cij . Задача позволяет определить, какой механизм и на какую работу надо назначить, чтобы добиться максимальной производительности; - задача о сокращении производства с учетом суммарных расходов на изготовление и транспортировку продукции; - увеличение производительности автомобильного транспорта за счет минимизации порожнего пробега. Уменьшение порожнего пробега сократит количество автомобилей для перевозок, увеличив их производительность. ПРИМЕР. Выбор оптимального варианта использования производственного оборудования. На предприятии имеются три группы станков, каждая из которых может выполнять пять операций по обработке деталей (операции могут выполняться в любом порядке). Максимальное время работы каждой группы станков соответственно равно 100, 250, 180 часов. Каждая операция должна выполняться соответственно 100, 120, 70, 130 часов. Определить, сколько времени и на какую операцию нужно использовать каждую группу станков, чтобы обработать максимальное количество деталей. Производительность каждой группы станков на каждую операцию задана матрицей 3 5 11 10 5 5 10 15 3 2 4 8 6 12 10 1. Оптимальное распределение оборудования Оборудование m различных видов нужно распределить между n рабочими участками. Производительность одной единицы оборудования i-го вида на j-м рабочем участке равна pij
; Данная задача относится к классу транспортных задач при условии, что производительность линейно зависит от количества используемого оборудования. Поставщиками в задаче являются различные виды оборудования, потребителями - рабочие участки. Предложение определяется запасом оборудования каждого вида, спрос - потребностью в нем на рабочем участке. Обозначим через xij
число единиц оборудования i-го вида, выделенное на j-й рабочий участок, Построенная модель является сбалансированной. Если запас оборудования и потребность в нем не равны, то переход к сбалансированной модели осуществляется с помощью преобразований, изложенных в пункте 4.7. В данной задаче требуется максимизировать целевую функцию P, представляющую суммарную производительность. Для перехода к стандартной транспортной модели надо заменить функцию P на противоположную функцию Z=-P, которую нужно будет минимизировать. При решении в транспортной таблице вместо тарифов на перевозки запишутся производительности pij , взятые с противоположным знаком. Далее задача решается как обычно. 2. Формирование оптимального штата фирмы. Фирма набирает штат сотрудников. Она располагает n группами различных должностей по bj
вакантных единиц в каждой группе, Предположим, что общее число кандидатов соответствует числу вакантных должностей. (если это не так, то следует просто проделать преобразования пункта 4.7.) Тогда данная задача соответствует транспортной модели. В роли поставщиков выступают группы кандидатов, а в роли потребителей - группы должностей. Предложением является число кандидатов в каждой группе, спросом - число вакансий в каждой группе должностей. В качестве тарифов на перевозки рассматриваются затраты на переобучение. Обозначим через xij
число кандидатов из i-й группы, назначенных на j-ю должность, Методы решения этой задачи такие же, как и транспортной задачи. ПРИМЕР 2. Имеется m видов сельскохозяйственных культур и n хозяйств, где их можно выращивать. Из-за различных условий доход, получаемый с 1 га посева одной и той же культуры, в различных хозяйствах неодинаков. Обозначим через cij доход для i-й культуры и j-го хозяйства. Общие площади посева культур ai (i=1,m) и площади пашни в хозяйствах bj (j=1,n) известны. Требуется составить такой план размещения сельскохозяйственных культур по хозяйствам, чтобы общий доход был максимальным. Обозначим площадь посева i-й культуры в j-м хозяйстве через xij . Тогда получаем задачу: Найти план X=(xij
) такой, что При условиях: а) план посева по каждой культуре должен быть выполнен б) пашня в хозяйствах должна быть использована полностью в) 4. Нелинейное программирование Во многих экономических моделях исследования операций зависимости между постоянными и переменными факторами лишь в первом приближении можно считать линейными, более детальное рассмотрение позволяет обнаружить их нелинейность. Как правило, такие показатели, как прибыль, себестоимость, капитальные затраты на производство и др., в действительности зависят от объема производства, расхода ресурсов и т.п. нелинейно. В этом случае возникает задача нелинейного программирования (ЗНП). 4.1 Общая задача нелинейного программирования В общем виде ЗНП формулируется следующим образом: где Решить задачу нелинейного программирования - это значит найти такие значения управляющих переменных Для задачи нелинейного программирования, в отличие от линейных задач, нет единого метода решения. В зависимости от вида целевой функции (4.1) и ограничений (4.2) разработано несколько специальных методов решения, к которым относятся методы множителей Лагранжа, квадратичное и выпуклое программирование, градиентные методы, ряд приближенных методов решения, графический метод. Процесс составления математической модели ЗНП принципиально не отличается от составления модели ЗЛП. Рассмотрим несколько примеров. Задача
4.1.
На m предприятиях выпускается некоторый продукт. Себестоимость единицы этого продукта на каждом из указанных предприятий есть Предприятия должны обеспечить n потребителей с потребностями Требуется определить такой план распределения выпуска продукта предприятиями и план перевозок его потребителям, чтобы суммарная себестоимость выпуска и стоимость перевозки была минимальной. Математическую модель задачи
. Пусть Для удобства запишем данные и искомые величины задач в виде таблицы: Табл.4.1
Система ограничений задачи: Целевая функция f запишется так: Вместо Надо найти минимальное значение функции (6.5) на множестве допустимых решений (4.4). Задача
4.2.
На производство некоторого продукта расходуется два вида ресурсов. Определите оптимальное распределение величин затрачиваемых ресурсов, если цена ресурса первого вида 30 рублей, второго - 40 рублей, а всего выделено на производство 120 рублей. Известно, что из количества x1
первого ресурса и x2
второго ресурса можно получить Вообще, функция где Функция y выведена в предположении, что существует только два ресурса: x1
- труд, x2
- капитал, где Математическая модель задачи . Пусть x1 - количество ресурсов вида 1, x2 - количество ресурсов вида 2. Система ограничений: Целевая функция: Требуется найти наибольшее значение функции (4.7) на множестве решений системы (4.6). 4.2. Некоторые особенности решения задач нелинейного программирования Для решения ЗНП существенно знать: 1) выпукло или не выпукло множество допустимых решений задачи; 2) является ли целевая функция выпуклой или вогнутой или она не относится ни к тому, ни к другому классу. Говорят, что множество выпукло , если оно вместе с любыми своими точками А и В содержит и все точки отрезка АВ. Примерами выпуклых множеств в пространстве могут служить сфера, пирамида, призма и т.д. В невыпуклом множестве можно указать хотя бы две точки, такие, что не все точки отрезка АВ принадлежат рассматриваемому множеству. Примером невыпуклого множества в пространстве является тор. Функцию y = f(x) одной переменной будем называть выпуклой , если отрезок, соединяющий две любые точки ее графика, принадлежит графику или расположен выше его. Функция вогнута , если отрезок, соединяющий две любые точки графика, принадлежит ему или расположен ниже его. Аналогично можно сформулировать определения понятий вогнутой и выпуклой функций нескольких переменных. Говорят, что гиперповерхность z=f(x1 ,x2 ,...,xn ) выпуклая, если отрезок, соединяющий две ее любые точки, лежит на поверхности или выше ее. Гиперповерхность z=f(x1 ,x2 ,...,xn ) вогнута, если отрезок, соединяющий две ее любые точки, лежит на поверхности или ниже ее. Пусть дана функция Говорят, что функция Точка Пусть функция Точка Необходимое условие экстремума.
Если в точке Как и в случае одной переменной, необходимое условие не является достаточным для того, чтобы стационарная точка была точкой экстремума. Для получения достаточных условий следует определить в стационарной точке знак дифференциала второго порядка. Дифференциал второго порядка обозначается Достаточные условия экстремума: а) в стационарной точке б) если в) если Для функции двух переменных Найдем значения частных производных второго порядка в стационарной точке Обозначим через Тогда достаточные условия экстремума функции двух переменных имеют вид: а) если б) если в) если Если область Ф замкнута и ограничена, то дифференцируемая функция Таким образом, чтобы найти наибольшее (наименьшее) значение функции 1) найти все стационарные точки внутри области Ф и вычислить значения функции в них: 2) исследовать функцию на экстремум на границе области Ф; 3) сравнить значения функции, полученные в п.1 и п.2: наибольшее (наименьшее) из этих чисел и будет наибольшим (наименьшим) значением функции во всей области. Граница области Ф аналитически может быть задана системой уравнений (условий) относительно переменных Условный экстремум.
Пусть необходимо найти экстремум функции Предполагается, что функции Говорят, что в точке Для задач линейного программирования характерно следующее: - Множество допустимых решений выпукло. Это выпуклое множество имеет конечное число вершин, называемых обычно крайними (угловыми( точками. - Множество всех точек - Локальный max или min является также глобальным max или min целевой функции на множестве допустимых решений, т.е. не существует локального оптимума целевой функции, отличного от глобального оптимума. - Если оптимальное значение целевой функции ограничено, то по крайней мере одна крайняя точка множества допустимых решений является оптимальным решением. Кроме того, начав с произвольной вершины множества допустимых решений, можно достичь оптимума за конечное число шагов, причем на каждом шаге совершается переход только в соседнюю вершину. Окончательно данная вершина является оптимальной тогда и только тогда, когда значение целевой функции в ней по крайней мере не меньше, чем значения целевой функции во всех примыкающих вершинах. У произвольной задачи нелинейного программирования некоторые или все приведенные выше свойства ЗЛП отсутствуют. Поэтому ЗНП гораздо сложнее задач линейного программирования и для них не существует общего универсального метода решения (аналогично симплексному методу). 4.4. Метод множителей Лагранжа Среди задач (4.1)-(4.3) особое место занимают задачи типа для решения которых можно воспользоваться классическим методом оптимизации Лагранжа, или методом разрешающих множителей. При этом предполагают, что функции Можно указать следующий порядок решения задачи (6.10)-(6.11) методом множителей Лагранжа: 1) составить функцию Лагранжа: L здесь 2)Найти частные производные функции Лагранжа по всем переменным Эта система состоит из m + n уравнений. Решить эту систему (если это окажется возможным) и найти таким образом все стационарные точки функции Лагранжа; 3) из стационарных точек, взятых без координат В основе метода Лагранжа решения классической задачи оптимизации (6.10)-(6.11) лежит утверждение, что если является решением системы (6.13). Следовательно, решая систему (6.13), получаем множество стационарных точек, в которых функция Z может иметь экстремальные значения. При этом, как правило, неизвестен способ определения точек глобального максимума или минимума. Однако, если решения системы найдены, то для определения глобального экстремума достаточно найти значения функции в соответствующих точках области определения. Множителям Лагранжа можно придать следующий экономический смысл. Если Пример 6.4. По плану производства продукции предприятию необходимо изготовить 200 изделий. Эти изделия могут быть изготовлены двумя технологическими способами. Производственные затраты на изготовление n изделий первым способом равны 4n+n2 , а для второго способа - 8n+n2 . Сколько изделий надо изготовить каждым способом, чтобы общие затраты на производство продукции были бы минимальными. Решение. Обозначим число изделий, изготовленных первым способом через х1 , вторым способом - х2. Тогда суммарные затраты на изготовление продукции по плану составят: f (x1 ,x2 )=4x1 +x1 2 +8x2 +x2 2 При этом общее число изделий должно быть равно 200, т.е. x1 +x2 =200 Получим математическую модель задачи: f(x1
,x2
)=4x1
+x1
2
+8x2
+x2
2
x1 +x2 =200 x1
Для ее расчета применим метод множителей Лагранжа. Составим функцию Лагранжа L(x1,
x2
, Найдем частные производные функции L по x1
,x2
, имеем систему: Исключим из этой системы Таким образом, по необходимому условию экстремума дифференцируемой функции получим стационарную точку М(101,99) возможного условного экстремума функции f(x1 ,x2 ). Дальнейшее исследование этой точки будем проводить как и в случае безусловного экстремума. Найдем частные производные второго порядка: Для рассматриваемого примера производные постоянны и не зависят от значений х1 и х2 . Поэтому a11 =2, a12 =a21 =0, a22 =2 (см. 6.8). Вычислим определитель Так как определитель больше нуля и a11 =2>0, следовательно, в точке М(101,99) функция f(x1 ,x2 ) имеет минимум: Ответ. При изготовлении 101 детали первым способом и 99 деталей вторым способом затраты на производство будут минимальными и равными 21198 денежным единицам. Литература. 1. Исследование операций в экономике. Под редакцией проф. Н.Ш. Кремера. М., «Банки и биржи», ЮНИТИ, 1997. 2. Хэмди А. Таха. Введение в исследование операций. Издательский дом «Вильямс»,М., С-П, Киев, 2001. 3. Акулич И.Л. Математическое программирование в примерах и задачах. М., 1986. 4. Эддоус М., Стэнфилд Р. Методы принятия решений. М., ЮНИТИ, 1997. 5. Банди Б. Основы линейного программирования. М., 1989. 6. Фомин Г.П. математические методы и модели в коммерческой деятельности. М., Финансы и статистика, 2001. 7. Вентцель Е.С. Исследование операций. М., Наука, 1980. КОНТРОЛЬНЫЕ ВОПРОСЫ 1. Цель и задачи исследования операций. 2. Модели и моделирование. Экономико-математические модели. 3. Математическое программирование. Модель задачи математического программирования. 4. Назовите основные разделы математического программирования. Дайте их краткую характеристику. 5. Приведите математическую формулировку основной задачи линейного программирования. 6. Какое решение называется допустимым и какое оптимальным? 7. Приведите примеры экономических задач. 8. Симметричная форма записи. Как от общей ЗЛП перейти к симметричной задаче? 9. Запишите закрытую модель транспортной задачи. 10. Запишите открытую модель транспортной задачи. 11. Дайте определение цикла свободной клетки. Сколько циклов существует у одной свободной клетки? 12. Что называется оценкой свободной клетки? 13. В каком случае план можно улучшать и как это сделать? 14. Как свести открытую модель транспортной задачи к закрытой? 15. В каком случае план транспортной задачи считается вырожденным? Как с этим бороться? 16. Как можно построить начальное опорное решение транспортной задачи? 17. Дайте экономическую интерпретацию метода потенциалов решения транспортной задачи. 18. Приведите примеры экономических задач, сводящихся к транспортным моделям. 19. Как решаются транспортные задачи, имеющие некоторые усложнения в постановке? Контрольная работа по экономико - математическим методам ТАБЛИЦА для определения индивидуального задания контрольной работы Последняя цифра номера зачетной книжки 1 2 3 4 5 6 7 8 9 0 11 12 13 14 15 16 17 18 19 20 1 36 37 38 39 40 21 22 23 24 25 60 41 42 43 44 45 46 47 48 49 П р 01 02 03 04 05 06 07 08 09 10 е 2 26 27 28 29 30 31 32 33 34 35 д 50 51 52 53 54 55 56 57 58 59 п о 01 02 03 04 05 06 07 08 09 10 с 3 27 28 29 30 31 32 33 34 35 36 л 52 53 54 55 56 57 58 59 60 41 е д 11 12 13 14 15 16 17 18 19 20 н 4 37 38 39 40 21 22 23 24 25 26 я 42 43 44 45 46 47 48 49 50 51 я 01 02 03 04 05 06 07 08 09 10 5 28 29 30 31 32 33 34 35 36 37 54 55 56 57 58 59 60 41 42 43 11 12 13 14 15 16 17 18 19 20 ц 6 38 39 40 21 22 23 24 25 26 27 и 44 45 46 47 48 49 50 51 52 53 ф р 01 02 03 04 05 06 07 08 09 10 а 7 23 24 25 26 27 28 29 30 31 32 45 46 47 48 49 50 51 52 53 54 11 12 13 14 15 16 17 18 19 20 8 25 26 27 28 29 30 31 32 33 34 60 41 42 43 44 45 46 47 48 49 01 02 03 04 05 06 07 08 09 10 9 35 36 37 38 39 40 21 22 23 24 50 51 52 53 54 55 56 57 58 59 11 12 13 14 15 16 17 18 19 20 0 26 27 28 29 30 31 32 33 34 35 42 43 44 45 46 47 48 49 50 51 Номера задач контрольной работы определяются по соответствующей таблице с помощью двух последних цифр номера зачетной книжки студента. Например, для студента, имеющего зачетную книжку с номером 87128, на пересечении горизонтальной колонки 2 и столбца 8 таблицы указаны следующие номера задач его индивидуального задания контрольной работы: 08, 33, 57. ЗАДАЧИ КОНТРОЛЬНОЙ РАБОТЫ Задачи 01 – 10 Используя графический метод, найти решение следующей задачи линейного программирования: Значения параметров a, b, c приведены в таблице 1. Задачи 11 – 20 Используя графический метод, найти решение следующей задачи линейного программирования: Значения параметров a, b, c приведены в таблице 1. Табл. 1
Задачи № 21 – 40 Ниже приведена таблица, в которой указаны запасы Табл. 2
Задача 21 Задача 22 Задача 23
Задача 24 Задача 25 Задача 26
Задача 27 Задача 28 Задача 29
Задача 30 Задача 31 Задача 32
Задача 33 Задача 34 Задача 35
Задача 36 Задача 37 Задача 38
Задача 39 Задача 40
Задачи № 41 – 60 По плану производства продукции предприятию необходимо изготовить в изделий. Эти изделия могут быть изготовлены двумя технологическими способами. Производственные затраты на изготовление n изделий первым способом равны Табл. 3
|