Формирование и решение линейных задач

Вид работы:

Курсовая работа (т)
Предмет:

Менеджмент
Язык:

Русский
,
Формат файла:
MS Word

153,14 Кб
Опубликовано:

2012-10-13

Скачать курсовую работу Читать текст online Заказать курсовую
*Помощь в написании! Посмотреть все курсовые работы

Вы можете узнать стоимость помощи в написании студенческой работы.

Формирование и решение линейных задач

1. Моделирование экономических систем. Основные понятия и определения

1.1 Возникновение и развитие системных представлений

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

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

В различных сферах человеческой деятельности возникли различные подходы и соответствующие методы решения специфических проблем, которые получили различные названия: в военных и экономических вопросах - «исследование операций», в политическом и административном управлении - «системный подход», в философии «диалектический материализм», в прикладных научных исследованиях - «кибернетика». Позже стало ясно, что все эти теоретические и прикладные дисциплины образуют как бы единый поток, «системное движение», которое постепенно оформилось в науку, получившую название «системный анализ». В настоящее время системный анализ является самостоятельной дисциплиной, имеющей свой объект деятельности, свой достаточно мощный арсенал средств и свою прикладную область. Являясь по существу прикладной диалектикой, системный анализ использует все средства современных научных исследований - математику, моделирование, вычислительную технику и натурные эксперименты.

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

Многие довольно часто недооценивают работу, связанную с формулировкой задачи. Однако многие специалисты полагают, что «хорошо поставить задачу - значит на половину ее решить». Хотя в большинстве случаев заказчику кажется, что он уже сформулировал свою проблему, системный аналитик знает, что предлагаемая клиентом постановка задачи является моделью его реальной проблемной ситуации и неизбежно имеет целевой характер, оставаясь приблизительной и упрощенной. Поэтому необходимо проверить эту модель на адекватность, что приводит к развитию и уточнению первоначальной модели. Очень часто первоначальная формулировка изложена в терминах не тех языков, которые необходимы для построения модели.

1.2 Модели и моделирование. Классификация моделей

моделирование линейный программирование задача

Первоначально моделью называли некое вспомогательное средство, объект, который в определенных ситуациях заменял другой объект. Например, манекен в определенном смысле заменяет человека, являясь моделью человеческой фигуры. Древние философы считали, что отобразить природу можно только с помощью логики и правильных рассуждений, т.е. по современной терминологии с помощью языковых моделей. Через несколько столетий девизом английского Научного общества стал лозунг: «Ничего словами!», признавались только выводы, подкрепленные экспериментально или математическими выкладками.

В настоящее время для постижения истины существует 3 пути:

теоретическое исследование;

эксперимент;

моделирование.

Моделью называется объект-заместитель, который в определенных условиях может заменять объект-оригинал, воспроизводя интересующие нас свойства и характеристики оригинала, причем имеет существенные преимущества:

- дешевизну;

- наглядность;

- легкость оперирования и т.п.

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

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

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

Другим принципом классификации целей моделирования может служить деление моделей на статические и динамические.

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

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

В распоряжении человека имеется два типа материалов для построения моделей - средства самого сознания и средства окружающею материального мира. Соответственно этому модели делятся на абстрактные (идеальные) и материальные.

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

1.3 Виды подобия моделей

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

Прежде всего, это подобие, устанавливаемое в процессе создания модели. Назовем такое подобие прямым. Примером такого подобия являются фотографии, масштабированные модели самолетов, кораблей, макеты зданий, выкройки, куклы и т.д.

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

Второй тип подобия между моделью и оригиналом называется косвенным. Косвенное подобие между оригиналом и моделью объективно существует в природе и обнаруживается в виде достаточной близости или совпадения их абстрактных математических моделей и вследствие этого широко используется в практике реального моделирования. Наиболее характерным примером может служить электромеханическая аналогия между маятником и электрическим контуром.

Оказалось, что многие закономерности электрических и механических процессов описываются одинаковыми уравнениями, различие состоит в разной физической интерпретации переменных, входящих в это уравнение. Роль моделей, обладающих косвенным подобием, очень велика и роль аналогий (моделей косвенного подобия) в науке и практике трудно переоценить. Аналоговые вычислительные машины позволяют найти решение почти всякого дифференциального уравнения, представляя собой, таким образом, модель, аналог процесса, описываемого этим уравнением. Использование электронных аналогов в практике определяется тем, что электрические сигналы легко измерить и зафиксировать, что дает известные преимущества модели.

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

Например, сигналом наступления кочевников у древних славян служили костры на курганах. Бумажные денежные знаки могут играть роль модели стоимости только до тех пор, пока в среде их обращения существуют правовые нормы, поддерживающие их функционирование. Керенки в настоящее время имеют только историческую ценность, но это не деньги, в отличие от царских золотых монет, которые представляют материальную ценность из-за наличия благородного металла. Особенно наглядна условность знаковых моделей: цветок в окне явочной квартиры Штирлица означал провал явки, ни сорт, ни цвет не имели никакого отношения к знаковой функции цветка.

1.4 Адекватность моделей

Модель, с помощью которой успешно достигается поставленная цель, будем называть адекватной этой цепи. Адекватность означает, что требования полноты, точности и правильности (истинности) модели выполнены не вообще, а лишь в той мере, которая достаточна достижения поставленной цели.

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

Приближенность модели не следует путать с адекватностью. Приближенность модели может быть очень высокой, но во всех случаях модель - это другой объект и различия неизбежны (единственной совершенной моделью любого объекта является сам объект). Величину, меру, степень приемлемости различия можно ввести, только соотнося его с целью моделирования. Так некоторые подделки произведений искусства даже эксперты не могут отличить от оригинала, но все-таки это всего лишь подделка, и с точки зрения вложения капитала не представляет никакой ценности, хотя для любителей искусства ничем не отличается от оригинала. У английского фельдмаршала Монтгомери во время войны был двойник, появление которого на разных участках фронта намеренно дезинформировало разведку немцев.

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

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

2. Математические модели и методы их расчета

2.1 Понятие операционного исследования

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

Решения человек принимал всегда и во всех сферах своей деятельности. Раньше хотели, чтобы принимаемые решения всегда были правильными. Теперь принято говорить, что решения должны быть оптимальными. Чем сложнее объект управления, тем труднее принять решение, и, следовательно, тем легче допустить ошибку. Вопросам принятия решений на основе применения ЭВМ и математических моделей посвящена новая наука «Исследование операций», приобретающая в последние годы все более обширное поле приложений. Эта наука сравнительно молодая, ее границы и содержание нельзя считать четко определенными.

Предмет под названием «Исследование операций» входит в программу элитарных вузов, но не всегда в этот термин вкладывается одно и то же содержание. Некоторые ученые под «исследованием операций» понимают, главным образом, математические методы оптимизации, такие как линейные, нелинейные, динамическое программирование. Другие к исследованию операций подходят с позиции теории игр и статистических решений. Наконец, некоторые ученые вкладывают в понятие «исследование операций» чрезмерно широкий смысл, считая ее основой системного анализа и «наукой наук».

Под термином «исследование операций» мы будем понимать применение математических, количественных методов для обоснования решений во всех областях целенаправленной человеческой деятельности.

Окончательно термин «исследование операций» закрепился в конце Второй мировой войны, когда в вооруженных силах США были сформированы специальные группы математиков и программистов, в задачу которых входила подготовка решений для командующих боевыми действиями. В дальнейшем исследование операций расширило область своих применений на самые разные области практики: экономика, транспорт, связь и даже охрана природы.

Выбор задачи - важнейший вопрос. Какие основные требования должна удовлетворять задача? Таких требований два:

должно существовать, как минимум, два варианта ее решения (ведь если вариант один, значит и выбирать не из чего);

надо четко знать в каком смысле искомое решение должно быть наилучшим (кто не знает, куда ему плыть - тому нет и попутного ветра).

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

Хорошую модель составить не просто. Известный математик Р. Беллман сказал так: «Если мы попытаемся включить в нашу модель слишком много черт действительности, то захлебнемся в сложных уравнениях; если слишком упростим ее, то она перестанет удовлетворять нашим требованиям». Таким образом, исследователь должен пройти между западнями Переупрощения и болотом Переусложнения. Для выполнения успеха моделирования надо выполнить три правила, которые, по мнению древних, являются признаками мудрости. Эти правила применительно к задачам математического моделирования и формулируются так: учесть главные свойства моделируемого объекта; пренебрегать его второстепенными свойствами; уметь отделить главные свойства от второстепенных.

Составление модели - это искусство, творчество. Древние говорили: «Если двое смотрят на одно и то же, это не означает, что оба видят одно и то же». И слова древних греков: «Если двое делают одно и то же, это не значит, что получится одно и то же». Эти слова в полной мере относятся к составлению математических моделей. Если математическая модель - это диагноз заболевания, то алгоритм - это метод лечения.

Можно выделить следующие основные этапы операционного исследования:

наблюдение явления и сбор исходных данных;

постановка задачи;

 построение математической модели;

 расчет модели;

 тестирование модели и анализ выходных данных. Если полученные результаты не удовлетворяют исследователя, то следует либо вернуться на этап 3, т.e. предложить для решения задачи другую математическую модель; либо вернуться на этап 2, т.e. поставить задачу более корректно;

 применение результатов исследований.

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

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

Экономико-математическая модель - это математическая модель, предназначенная для исследования экономической проблемы.

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

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

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

2.2 Классификация и принципы построения математических моделей

Можно выделить следующие основные этапы построения математической модели:

Определение цели, т.e. чего хотят добиться, решая поставленную задачу.

Определение пapaметров модели, т.е. заранее известных фиксированных факторов, на значения которых исследователь не влияет.

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

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

 Выявление неизвестных факторов, т.е. величин, которые могут изменяться случайным или неопределенным образом.

 Выражение цели через управляющие переменные, параметры и неизвестные факторы, т.e. формирование целевой функции, называемой также критерием эффективности или критерием оптимальности задачи.

Введем следующие условные обозначения:

a - параметры модели;

x - управляющие переменные или решения;

X - область допустимых решений;

x - случайные или неопределенные факторы;

W - целевая функция или критерий эффективности (критерий оптимальности).

W=W (x, a,x)

В соответствии с введенными терминами, математическая модель задачи имеет следующий вид:

W=W (x, a,x)® max (min) (2.1)

x Î X

Решить задачу - это значит найти такое оптимальное решение x^X, чтобы при данных фиксированных параметрах  и с учетом неизвестных  факторов значения критерия эффективности W было по возможности максимальным (минимальным).

W^=W (x^, a,x) = max (min) W (x, a,x)

x Î X

Таким образом, оптимальное решение - это решение, предпочтительное перед другими по определенному критерию эффективности (одному или нескольким).

Перечислим некоторые основные принципы построения математической модели:

Необходимо соизмерять точность и подробность модели, во-первых, с точностью тex исходных данных, которыми располагает исследователь, и, во-вторых, с теми результатами, которые требуется получить.

Математическая модель должна отражать существенные черты исследуемого явления и при этом не должна его сильно упрощать.

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

 Любая сложная система всегда подвергается малым внешним и внутренним воздействиям, следовательно, математическая модель должна быть устойчивой (сохранять свойства и структуру при этих воздействиях).

По числу критериев эффективности математические модели делятся на однокритериальные и многокритериальные. Многокритериальные математические модели содержат два и более критерия.

По учету неизвестных факторов математические модели делятся на детерминированные, стохастические и модели с элементами неопределенности.

В стохастических моделях неизвестные факторы - это случайные величины, для которых известны функции распределения и различные статистические характеристики (математическое ожидание, дисперсия, среднеквадратическое отклонение и т.п.). Среди стохастических характеристик можно выделить:

- модели стохастического программирования, в которых либо в целевую функцию (2.1), либо в ограничения (2.2) входят случайные величины;

- модели теории случайных процессов, предназначенные для изучения процессов, состояние которых в каждый момент времени является случайной величиной;

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

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

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

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

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

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

Hелинейные модели - это модели, в которых либо целевая функция, либо какое-нибудь из ограничений (либо все ограничения) нелинейны по управляющим переменным. Для нелинейных моделей нет единого метода расчета. В зависимости от вида нелинейности, свойств функции и ограничений можно предложить различные способы решения. Однако может случится и так, что для поставленной нелинейной задачи вообще не существует метода расчета. В этом случае задачу следует упростить, либо сведя ее к известным линейным моделям, либо просто линеаризовав модель.

В динамических моделях, в отличие от статических линейных и нелинейных моделей, учитывается фактор времени. Критерий оптимальности в динамических моделях может быть самого общего вида (и даже вообще не быть функцией), однако для него должны выполняться определенные свойства. Расчет динамических моделей сложен, и для каждой конкретной задачи необходимо разрабатывать специальный алгоритм решения.

Графические модели - используются тогда, когда задачу удобно представить в виде графической структуры.

3. Примеры задач линейного программирования

.1 Транспортная задача

уголь, добываемый в нескольких месторождениях, отправляется ряду потребителей. нам известно, сколько угля добывается в каждом из месторождений, скажем за месяц и сколько его требуется на тот же срок каждому из потребителей. Известны расстояния между месторождениями и потребителями, а также условия сообщения между ними. Учитывая эти данные. Можно подсчитать, во что обходится перевозка каждой тонны угля из любого месторождения в любой пункт потребления. Требуется при этих условиях спланировать перевозки угля таким образом, чтобы затраты на них были минимальными.

Пусть для простоты заданы всего 4 месторождения М₁, М², М³, М₄, причем их ежемесячная добыча составляет a₁, а₂, а₃, а₄ тонн угля. Предположим далее, что этот уголь надо доставить в пункты потребления b₁, b₂, b₃, b₄, b₅, соответственно с ежемесячными потребностями этих пунктов. Будем считать, что общее производство угля равно суммарной потребности в нем (сбалансированность планов): a₁, а₂, а₃, а₄ = b₁, b₂, b₃, b₄, b₅. Задача состоит в определении такого плана перевозок, при котором общая стоимость перевозок была бы наименьшей. Обозначим через x₁₁ количество угля (в тоннах), предназначенное к отправлению из M₁ в П₁; вообще через x_ij обозначим количество угля, отправляемого из месторождения M_i в пункт потребления П_j. Схема перевозок примет вид, изображенный в таблице 3.1

Таблица 3.1 Схема перевозок

	ПН	в П₁	в П₂	в П₃	в П₄	в П₅	Всего
ПО							отправлено
из М₁		х₁₁	х₁₂	х₁₃	х₁₄	х₁₅	a₁
из М²		х₂₁	х₂₂	х₂₃	х₂₄	х₂₅	а₂
из М³		х₃₁	х₃₂	х₃₃	х₃₄	х₃₅	а₃
из М₄		х₄₁	х₄₂	х₄₃	х₄₄	х₄₅	а₄
Всего привезено		b₁	b₂	b₃	b₄	b₅

3.1

х₁₁₊х₁₂₊х₁₃₊х₁₄₊х_{15
=}b₁ х₂₁₊х₂₂₊х₂₃₊х₂₄₊х_{25
=}b₂ х₃₁₊х₃₂₊х₃₃₊х₃₄₊х_{35
=}b₃ х₄₁₊х₄₂₊х₄₃₊х₄₄₊х_{45
=}b₄

3.2

х₁₁₊х₁₂₊х₁₃₊х₁₄₊х_{15
=}a₁ х₂₁₊х₂₂₊х₂₃₊х₂₄₊х_{25
=}a₂ х₃₁₊х₃₂₊х₃₃₊х₃₄₊х_{35
=}a₃ х₄₁₊х₄₂₊х₄₃₊х₄₄₊х_{45
=}a₄

Общее количество угля, привозимое в пункт П₁ из всех месторождений, будет х₁₁₊х₁₂₊х₁₃₊х₁₄₊х_15
=b₁; в другие пункты - П₂, П₃ и т.д. и примет вид уравнений 4.1. общее количество угля, вывозимое из М₁, будет: х₁₁₊х₁₂₊х₁₃₊х₁₄₊х_15
=a₁, примет вид 4.2. предполагаем, что стоимость перевозки прямо пропорциональна количеству перевозимого угля, т.е. стоимость перевозки x_ij тонн угля равна:

x _i_j= C _i_j.X _i_j

Общая стоимость S всех перевозок будет равна: (3.3)

S=c₁₁х₁₁+c₁₂х₁₂+c₁₃х₁₃+c₁₄х₁₄+c₁₅х₁₅+… +c₄₁х₄₁+c₄₂х₄₂+c₄₃х₄₃+c₄₄х₄₄+c₄₅х₄₅.

3.2 Общая формулировка задачи линейного программирования

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

эти значения удовлетворяли некоторой системе линейных уравнений или неравенств;

при этих значениях некоторая линейная функция (линейная форма) от этих неизвестных обращалась в минимум (максимум);

 эти значения были неотрицательными.

Задачами такого рода и занимается линейное программирование.

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

Дана система линейных уравнений:

ìа₁₁x₁+а₁₂x₂+… +а₁_nx_n₌b₁

ïа₂₁x₁+а₂₂x₂+… +а_2nx_n
=b₂

(I)í……………………

îа_m1x₁+а_m_²x₂+… +а_mnx_n₌b_m

и линейная функция ¦=c₁x₁+c₂x₂+… +c_nx_n(II).

Требуется найти такие неотрицательные решения х₁  0, х₂  0… х_n  0 (III) системы (I) при которых функция  принимает наименьшее (наибольшее) значение.

Уравнения (I) называются ограничениями данной задачи, уравнение (II) называется линейной формой, а уравнение (III), строго говоря, тоже являются ограничениями, однако их не принято так называть, поскольку они являются общими для всех задач линейного программирования, а не только конкретной задачи. Любое неотрицательное решение системы уравнений называется допустимым. Допустимое решение, дающее минимум функции , оптимальное решение (если оно существует) не обязательно единственно; возможны случаи, когда имеется бесчисленное множество оптимальных решений. Наконец, саму функцию  часто называют линейной формой или функцией цели.

Казалось бы, т.к. задача линейного программирования ставится как задача минимизации некоторой функции , то можно применить классический прием дифференциального исчисления. Однако частные производные равны коэффициентам при неизвестных, которые в «нуль» одновременно не обращаются.

3.3 Графическая интерпретация решения задач линейного программирования

Задачу линейного программирования (ЛП) можно решать аналитическими и графическими методами. Аналитические методы являются основой для решения задачи на ЭВМ. Их единственный недостаток состоит в том, что в отличие от графических методов, они недостаточно наглядны. Графические методы очень наглядны, но они пригодны лишь для решения задач на плоскости, т.е. когда размерность пространства К=2. Однако, учитывая большую наглядность графических методов, с их помощью рассмотрим идею решения задачи ЛП на примере задачи распределения ресурсов.

Однако прежде чем заняться решением, сделаем некоторые замечания. Пусть мы имеем систему m уравнения с n неизвестными (I).

Возможны следующие варианты:

Число неизвестных меньше, чем число уравнений n  m.

например: ì2x₁=4, в этом случае n=1;

îx₁=5, тогда m=2 (число линейно независимых уравнений). (3.4)

Очевидно, что система (3.4) решения не имеет, и она несовместна;

Число неизвестных равно числу уравнений n=m.

В этом случае система имеет единственное решение или не имеет ни одного. Заметим, что m равно числу линейно независимых уравнений.

Для системы: ì2x=10, n=1, m=1;

î6x=30.

Если число неизвестных больше числа уравнений, то система имеет бесчисленное множество решений. Пусть n  m. Например:

2x₁+x₂=2 (3.5)

Очевидно, что это уравнение прямой, и все значения x₁ и x₂, лежащие на этой прямой, являются решением уравнения (4.2). Значит уравнение (4.5) имеет бесчисленное множество решений.

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

a₁x₁+a₂x₂=b (3.6)

Преобразуем его к виду:

(3.7)

Запись (3.7) называют уравнением прямой в отрезках, что изображено на Рис. 3.1. Рассмотрим еще одну форму представления уравнения (3.6). Запишем это уравнение в виде:

a₂x₂=b-a₁x₁

или

x₂=F-kx₁(3.8)

Уравнение (3.8) изображено на рис. 3.2.

Вспомним неравенства. Если линейное уравнение с двумя переменными может быть представлено в виде прямой на плоскости, то неравенство вида:

a₁x₁+a₂x₂ £ b (3.9)

изображается как полуплоскость, показанная на рис. 4.1. На этом рисунке часть плоскости, удовлетворяющая неравенству, заштрихована. Координаты всех точек, принадлежащих заштрихованному участку, имеют такие значения x₁ и x₂, которые удовлетворяют заданному неравенству. Значит, эти значения составляют область допустимых решений (ОДР). Саму прямую считаем принадлежащей каждой из двух указанных полуплоскостей. Предположим теперь, что задано не одно неравенство, а система:

ìа₁₁x₁+а₁₂x₂£ b₁

ïа₂₁x₁+а₂₂x₂£ b₂

(3.10) îа_m₁x₁+а_m_²x₂£ b_m,

где первое неравенство определяет некоторую полуплоскость П₁, второе - полуплоскость П₂ и т.д.

если какая-либо пара чисел (x₁, x₂) удовлетворяет всем неравенствам (4.10), то, соответствующая точка Р(x₁, x₂), принадлежит всем полуплоскостям П₁, П₂,… П_m одновременно. Другими словами, точка Р принадлежит пересечению (общей части) полуплоскостей П₁, П₂,… П_m, т.е. некоторой многоугольной области М (Рис. 4.3), которая является ОДР. Вдоль контура области изображены штрихи, идущие внутрь области. Они одновременно указывают, с какой стороны от данной прямой лежит соответствующая полуплоскость, то же самое указано и с помощью стрелок на каждой линии. Сразу же отметим, что ОДР не всегда бывает, ограничена: в результате пересечения нескольких полуплоскостей может возникнуть и неограниченная область (Рис. 4.4). Возможен и случай, когда область допустимых решений (ОДР) пуста. Это означает, что система (5.7) противоречива (Рис. 4.5). Многоугольник ОДР обладает весьма важным свойством: он является выпуклым.

фигура называется выпуклой, если вместе с любыми двумя своими точками А и В, она содержит и весь отрезок АВ.

В случае трех неизвестных, каждое уравнение представляет собой плоскость в пространстве. Каждая плоскость разбивает все пространство на два полупространства. Система неравенств определяет в пространстве выпуклый объемный многогранник, который представляет ОДР.

4. Методы решения задач линейного программирования

4.1 Общая и основная задачи линейного программирования

К математическим задачам линейного программирования приводят исследования конкретных производственно-хозяйственных ситуаций, которые в том или ином виде интерпретируются как задачи об оптимальном использовании ограниченных ресурсов (задача о раскрое, смесях и т.д.).

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

Oбщей задачей линейного программирования называется задача, которая coстоит в определении максимального (минимального) значения функции:

(4.1)

при условии:

(4.2)

(4.3)

X_j ³ 0 (j=1, 1; 1 £ n) (4.4)

где a_i_j, b_i, с_j - заданные постоянные величины и k  m.

Функция (4.1) называется целевой функцией (или линейной формой) задачи (4.1) - (5.4), а условия (4.2) - (4.4) - ограничениями данной задачи.

Стандартной (или симметричной) задачей линейного программирования называется задача, которая состоит в определении максимального значения функции (4.1) при выполнении условий (4.2) и (4.4), где k=m и 1=n.

Канонической (или основной) задачей линейного программирования называется задача, которая состоит в определении максимального значения функции (4.1) при выполнении условий (4.3) и (4.4), где k=0 и 1=n.

Совокупность чисел Х^ = (x₁, x₂,…, x_n), удовлетворяющих ограничениям задачи (4.2) - (5.4), называется допустимым решением (или планом).

План Х^ = (x₁, x₂,…, x_n), при котором целевая функция задачи (4.1) принимает свое максимальное (минимальное) значение, называется оптимальным.

Значение целевой функции (4.1) при плане X будем обозначать через F(X). Следовательно, Х^ - оптимальный план задачи, если для любого X выполняется неравенство F(X)  F(Х^) (соответственно F(X)  F(Х^)).

4.2 Геометрический метод решения задач линейного программирования

Перепишем основную задачу линейного программирования в векторной форме: найти максимум функции

F=CX (4.5)

при условиях:

x₁P₁+x₂P₂+… +x_nP_n = P₀ (4.6)

Х ³0 (4.7)

где C = (с₁, с₂,…, с_n), Х = (х₁, х₂,…, х_n); СХ - скалярное произведение; Р₁, Р₂,…, Р_n и Р₀ - m-мерные вектор-столбцы, составленные из коэффициентов при неизвестных и свободных членах системы уравнений задачи.

План X = (х₁, х₂,…, х_n) называется опорным планом основной задачи линейного программирования, если система векторов P_j, входящих в разложение (4.6) с положительными коэффициентами x_j, линейно независима.

Непустое множество планов основной задачи линейного программирования образует выпуклый многогранник. Каждая вершина этого многогранника определяет опорный план. В одной из вершин многогранника решений (т.е. для одного из опорных планов) значение целевой функции является максимальным (при условии, что функция ограничена сверху на множестве планов). Если максимальное значение функция принимает более чем в одной вершине, то это же значение она принимает в любой точке, являющейся выпуклой линейной комбинацией данных вершин.

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

F=c₁x₁ + с₂х₂ (4.8)

при условиях:

a_i₁x₁+a_i₂x₂ £ b_i (i=1, k) (4.9)_j ³ 0 (i=1,2) (4.10)

Каждое из неравенств (4.9), (4.10) системы ограничений задачи геометрически определяет полуплоскость соответственно с граничными прямыми a_i₁+a_i₂=b_i (i=1, k), x₁=0 и x₂=0. В этом случае, если система неравенств (4.9), (4.10) совместна, область ее решений есть множество точек, принадлежащих всем указанным полуплоскостям.

Так как множество точек пересечения данных полуплоскостей выпуклое, то областью допустимых решений (ОДР) задачи (4.8) - (4.10) является выпуклое множество, которое называется многоугольником решений. Стороны этого многоугольника лежат на прямых, уравнения которых получаются из исходной системы ограничений заменой знаков неравенств на знаки точных равенств.

Таким образом, исходная задача линейного программирования состоит в нахождении такой точки области допустимых решений, в которой целевая функция F принимает максимальное значение. Эта точка существует тогда, когда многоугольник решений не пуст и на нем целевая функция ограничена сверху. При указанных условиях в одной из вершин ОДР целевая функция принимает максимальное значение. Для определения данной вершины построим линию уровня c₁x₁+c₂x₂=h (где h - некоторая постоянная), проходящую через ОДР, и будем ее передвигать в направлении вектора C=(с₁; с₂) до тех пор, пока она не пройдет через последнюю ее общую точку с многоугольником ОДР. Координаты указанной точки и определяю оптимальный план данной задачи.

При нахождении решения задачи могут встретиться случаи, изображенные на рис. 5.1-5.4. Рис. 5.1 характеризует такой случай, когда целевая функция принимает максимальное значение в единственной точке А. Из рис. 5.2 видно, что максимальное значение целевая функция принимает в любой точке отрезка АВ.тметим, что нахождение минимального значения линейной функции при данной системе ограничений отличается от нахождения ее максимального значения при тex же ограничениях лишь тем, что линия уровня c₁x₁+c₂x₂=h передвигается не в направлении вектора С, а в противоположном направлении.

Тот факт, что оптимальное решение находится в одной из вершин многоугольника ОДР, позволяет сделать еще два важных вывода:

если оптимальным решением являются координаты вершины многоугольника ОДР, значит, сколько вершин имеет ОДР, столько существует целевых функций и столько оптимальных решений по этим функциям может иметь задача.

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

Из рисунков можно сделать вывод, что вершина, координаты которой являются оптимальным решением, определяются углом наклона прямой, описывающей целевую функцию. Значит, каждая вершина будет соответствовать оптимальному решению для некоторой целевой функции. Итак, нахождение решения задачи линейного программирования (5.8) - (5.10) на основе ее геометрической интерпретации включает следующие этапы.

Этапы нахождения решения задачи линейного программирования:

Строят прямые, уравнения которых получаются в результате замены в ограничениях и знаков неравенств на знаки точных равенств.

Находят полуплоскости, определяемые каждым из ограничений задачи.

Находят многоугольник решений (ОДР).

Строят вектор C=(с₁; с₂).

Строят прямую c₁x₁+c₂x₂=h, проходящую через многоугольник решений.

Передвигают прямую c₁x₁+c₂x₂=h в направлении вектора С, в результате чего-либо находят точку (точки), в которой целевая функция принимает максимальное значение, либо устанавливают неограниченность сверху функции на множестве планов.

Определяют координаты точки максимума функции и вычисляют значение целевой функции в этой точке.

4.3 Графическое решение задачи распределения ресурсов

Пусть для двух видов продукции П₁ и П₂ требуются трудовые, материальные и финансовые ресурсы. Наличие ресурсов каждого вида и их нормы расхода, необходимые для выпуска единицы продукции, приведены в табл. 4.1.

Таблица 4.1

Характеристика	Вид продукции		располаг.
	П₁	П₂	ресурс
Резервы:
трудовые	1	4	14
материальные	3	4	18
финансовые	6	2	27
выпуск	1	1	-
прибыль	4	8,5	-
план	х₁	х₂	-

Составим математическую модель задачи

x₁+4x₂£ 14

3x₁+4x₂£ 18

6x₁+2x₂£ 27

x₁ ³ 0, x₂³

Математическая модель представляет собой систему линейных неравенств. Значит ОДР нашей задачи выпуклый многоугольник. Для удобства построения неравенства можно записать в форме, аналогичной уравнениям в отрезках:

x₁/14+x₂/7/2 £ 1

x₁/6+x₂/9/2 £ 1

x₁/9/2+x₂/27/2 £1

x₁ ³0, x₂³0.

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

=x₁+x₂ ® max (5.13)

Эту зависимость представим в виде x₂= F-x₁. Из графика данного уравнения (Рис. 5.1) следует, что tg= 1, при этом  =135^о, а величина F равна отрезку, отсекаемому прямой функции цели на оси координат. Если прямую перемещать параллельно самой себе в направлении, указанном стрелками, то эта величина будет возрастать. Очевидно, что оптимальным решением будут координаты точки С (x^₁; x^₂). При этом F=F^.

На основании рассмотренного, можно сделать исключительно важный вывод: оптимальным решением являются координаты вершин ОДР. А из этого вывода следует метод решения задачи линейного программирования.

Метод решения задачи линейного программирования:

Найти вершины ОДР, как точки пересечения ограничений.

Определить последовательно значения целевой функции в вершинах.

Вершина, в которой целевая функция приобретает оптимальное значение, является оптимальной вершиной.

Координаты оптимальной вершины являются оптимальными значениями искомых переменных.

Если направление целевой функции совпадает с направлением одной из сторон, то у задачи будет, по крайней мере, два решения. В таком случае говорят, что задача имеет альтернативные решения. А это значит, что одно и то же оптимальное значение целевой функции может быть получено при различных значениях переменных.

Тот факт, что оптимальное решение находится на вершине ОДР, дает еще два очень важных вывода:

если оптимальным решением являются координаты вершин ОДР, значит, сколько вершин имеет ОДР, столько оптимальных решений может иметь задача.

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

Как видно на Рис. 5.1, вершина, координаты которой являются оптимальным решением, определяется углом наклона прямой, описывающей целевую функцию. Значит, каждая вершина будет соответствовать оптимальному решению для некоторой целевой функции. Поясним это на рассмотренном ранее примере. Раньше мы находили оптимальное решение по максимизации суммарного выпуска F₁=x₁+x₂  max. Найдем оптимальные решения еще для четырех целевых функций:

F₂=x₂ ® max (максимизация выпуска продукции П₂)

F₃=x₁ ® max (максимизация выпуска продукции П₁)

F₄=4x₁+8,5x₂ ® max (максимизация прибыли)

F₅=(1+3+6) x₁+ (4+4+2) x₂ = 10х₁+10х₂® max (минимизация используемых ресурсов).

Для каждой целевой функции, так же как и для F₁, можно построить линию целевой функции и определить вершину, в которой целевая функция будет иметь оптимальное значение. Результаты решения задачи по пяти целевым функциям сведем в таблицу 4.2, из анализа которой можно сделать вывод: координаты каждой вершины могут быть оптимальным решением в некотором смысле.

Таблица 4.2

Целевая функция	Вершина	x₁	x₂	x₁+x₂	Прибыль	Используемый ресурс
F₁=x₁+x₂  max	C	4	1,5	5,5	28,75	55
F₂=x₂  max	A	0	3,5	3,5	29,75	35	D	4,5	0	4,5	18	45
F₄=4x₁+8,5x₂  max	B	2	3	5	33,5	50
F₅= 10х₁+10х₂	0	0	0	0	0	0

4.4 Симплексный метод

Симплексный метод или метод последовательного улучшения плана является одним из основных методов решения задач ЛП. название симплексный метод берет от слова «симплекс», которым создатель метода Р. Данциг обозначил наложенное на переменные x₁, x₂… x_n ограничение x₁+x₂+… +x_n=1.

В математике симплексом в k-мерном пространстве называется совокупность k+1 вершин.

Так для плоскости при k=2 симплексом будет треугольник; в пространстве при k=3 симплексом будет тетраэдр, имеющий 4 вершины.

С учетом этого понятия аналитический метод решения задачи ЛП называют симплекс-методом. Он основан на алгоритме направленного перебора вершин. Этот алгоритм обеспечивает переход от одной вершины к другой в таком направлении, при котором значение целевой функции от вершины к вершине улучшается.

Определение значения целевой функции и переменных в одной вершине считается итерацией.

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

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

4x₁+9x₂ £ 56

5x₁+3x₂ £ 37

x₁+2x₂ £2

обращающее в максимум линейную форму:

¦=3x₁+4x₂ (5.15)

Вначале перейдем от системы неравенств (5.14) к системе уравнений, добавив к левым частям неравенств неотрицательные переменные x₃, x₄, x₅. Мы получим:

4x₁+9x₂+x₃+0.x₄+0.x₅=56

5x₁+3x₂+0.x₃+x₄+0.x₅=37

îx₁+2x₂+0.x₃+0.x₄+x₅=2

¦=3x₁+4x₂+0.x₃+0.x₄+0.x₅ (5.17)

перепишем теперь систему (5.16) в виде системы 0-уравнений:

0=56 - (4x₁+9x₂+1.x₃+0.x₄+0.x₅)

=37 - (5x₁+3x₂+0.x₃+1.x₄+0.x₅)

=2 - (-x₁+2x₂+0.x₃+0.x₄+1.x₅

¦=0 - (-3x₁-4x₂-0.x₃-0.x₄-0.x₅) (5.19)

заметим, что система (5.18) может быть записана в виде одного векторного равенства:

0=B - (A₁x₁+A₂x₂+A₃x₃+A₄x₄+A₅x₅),

где вектор-столбец В имеет своими компонентами свободные члены, а векторы A₁, A₂,…, A₅ - коэффициенты при соответствующих переменных x₁, x₂, x₃, x₄, x₅.

Линейная форма имеет вид: ¦=3x₁+4x₂+0.x₃+0.x₄+0.x₅.

Векторы A₃, A₄, A₅ образуют базис. Это означает, что, присвоив х₁=0, х₂=0, получаем из (5.16) первое базисное решение: x₁=0; x₂=0; x₃=56; x₄=37; x₅=2.

При этом значение линейной формы =0. На основании (5.18) строим первую симплексную таблицу.

Симплексная таблица строится следующим образом:

В заглавной строке пишем последовательно векторы B, A₁, A₂, A₃, A_4, A₅. Слева добавляем колонку «Базисные векторы», рядом с ней колонку «С», в которой поставлены коэффициенты при базисных переменных в линейной форме, в данном случае величины С₃, С₄, С₅. В последней строке, называемой индексной, и обозначаемой через  j-С_j, проставляются числа, равные значению линейной формы, в соответствием с уравнением (j=1, 2, 3, 4, 5). В итоге мы имеем таблицу 4.1.

Таблица 4.1.

Базисные	Коэффициенты	Вектор свободных	3	4	0	0	0
векторы	линейной формы С	членов В	A₁	A₂	A₃	A₄	A₅
A₃	0	56	4	9	1	0	0
A₄	0	37	5	3	0	1	0
A₅	0	2	-1	2	0	0	1
Индексная строка  j-С_j		0	-3	-4	0	0	0

Это первая симплексная таблица, соответствующая первому базисному решению: x₁=0; x₂=0; x₃=56; x₄=37; x₅=2. Значение линейной формы, равное нулю, мы записываем в первой клетке индексной строки.

Т.к. мы решаем задачу на максимум, то из выражения линейной формы видно, что имеет смысл увеличить x₁ или x₂. Действительно, коэффициенты при этих переменных в скобках отрицательны (а по существу положительны), и если мы положим x₁0 или x₂0, то значение  увеличится. Но эти же коэффициенты с их знаками стоят в индексной строке.

Итак, мы приходим к следующему выводу: наличие в индексной строке отрицательных чисел при решении задачи на максимум свидетельствует о том, что нами оптимальное решение не получено, и то, что от табл. 4.1 надо перейти к следующей.

Переход к новой таблице, т.е. к новой улучшенной программе осуществляем следующим способом: в индексной строке находим наибольшее по абсолютному значению отрицательное (а при задаче на минимум - наибольшее положительное) число. В нашем примере этим числом будет - 4. Найденное число определяет ведущий или ключевой столбец. Затем мы делим свободные члены на положительные элементы ведущего столбца и выбираем из полученных отношений наименьшее. Наименьшее отношение определяет ведущую строку. В данном случае имеем:

Таким образом, ведущей строкой будет строка A₅. На пересечении ведущего столбца и ведущей строки стоит разрешающий элемент. В нашем случае - это число 2.

Теперь мы приступаем к составлению второй таблицы или второго плана. Вместо единичного вектора A₅ мы в базис вводим вектор A₂. Переход к новому базису, как это известно, эквивалентен элементарному преобразованию матрицы, элементами которой служат числа табл. 4.1. А именно: в новой таблице элемент строки, соответствующий элементу ведущей строки прежней таблицы, равен этому элементу ведущей строки, разделенному на разрешающий элемент. чтобы получить любой другой элемент новой симплексной таблицы, нужно от соответствующего элемента прежней таблицы отнять произведение элемента ведущей строки на элемент ведущего столбца, разделенное на разрешающий элемент. Например, элементу 4 (табл. 4.1) будет соответствовать элемент табл. 4.2:

Таким образом, мы переходим ко второй таблице (таблица 4.2). Указанные выше преобразования относятся к столбцам B, A₁, A₂, A₃, A_4, A₅.

Таблица 4.2.

Базисные	Коэффициенты	Вектор свободных	3	4	0	0	0
векторы	линейной формы С	членов В	A₁	A₂	A₃	A₄	A₅
A₃	0	47	17/2	0	1	0	-9/2
A₄	0	34	13/2	0	0	1	-3/2
A₅	4	1	-1/2	1	0	0	1/2
Индексная строка  j-С_j		4	-5	0	0	0	2

Из табл. 4.2 видно, что значение линейной формы возросло и теперь равно 4. Однако наличие в индексной строке отрицательных чисел свидетельствует о том, что это значение еще можно увеличить. Переходим к следующей симплексной таблице. число «5» определяет ведущий столбец. Находим ведущую строку. Для этого определяем:

Итак, разрешающим элементом будет 13/2. Вектор A₄ выводим из базиса и вводим вместо него вектор A₁. Пересчет коэффициентов осуществляем по указанным выше правилам и получаем таблицу 4.3.

Таблица 5.3

Базисные	Коэффициенты	Вектор свободных	3	4	0	0	0
векторы	линейной формы С	членов В	A₁	A₂	A₃	A₄	A₅
A₃	0	33/13	0	0	1	-17/13	-33/13
A₁	3	68/13	1	0	0	2/13	-3/13
A₂	4	47/13	0	1	0	1/13	5/13
Индексная строка  j-С_j		392/13	0	0	0	10/13	11/13

В индексной строке нет отрицательных элементов. Следовательно, мы получим оптимальную программу. Оптимальное решение:

x₁=68/13; x₂=47/13; x₃=33/13; x₄= x₅= 0.

.5 Анализ симплекс-таблиц

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

При оперативном управлении решается достаточно широкий и важный круг вопросов, которые возникают при ежедневном обеспечении производственного процесса. Мы рассмотрим лишь те вопросы оперативного управления, которые могут быть решены с помощью моделей, уже составленных при планировании. «Что будет, если пять человек из числа трудовых ресурсов отвлекут на другие работы? Что будет, если сырья поставят на 20% меньше? Какую продукцию следует выпускать, если изменились цены?» Рассмотрим, как находить ответы на эти вопросы на конкретном примере.

Допустим, предприятие должно выпускать продукцию четырех видов: П₁, П₂, П₃, П₄, используя для этого три вида ресурсов.

Таблица 4.1

Элемент	Вид продукции				Располагаемый
модели	П₁	П₂	П₃	П₄	ресурс
Ресурсы:
трудовые	1	1	1	1	16
сырье	6	5	4	3	110
оборудование	4	6	10	13	100
Прибыль с единицы продукта	60	70	120	130	-
План	х₁	х₂	х₃	х₄	-

x₁+ x₂+ x₃+ x₄£ 16

6x₁+ 5x₂+ 4x₃+ 3x₄£ 110

4x₁+ 6x₂+ 10x₃+ 13x₄ £100

x_j ³ 0, j³4

F = 60x₁+ 70x₂ + 120x₃+ 130x₄ ® max

От системы неравенств (5.20) перейдем к системе уравнений. Для этого, в каждое неравенство добавим по одной дополнительной переменной: y_i 0, i m. Тогда получим систему уравнений:

x₁+ x₂+ x₃+ x₄+ y₁=16

6x₁+ 5x₂+ 4x₃+ 3x₄+ y₂ =110

4x₁+ 6x₂+ 10x₃+ 13x₄ + y₃ =100_j ³ 0, j =4; y_i ³ 0, i= 1,3.

F = 60x₁+ 70x₂ + 120x₃+ 130x₄  max

Затем перепишем систему (5.21) в следующем виде:

₁=16 - x₁+ x₂+ x₃+ x₄)₂ =110 - 6x₁+ 5x₂+ 4x₃+ 3x₄)₃ =100 - 4x₁+ 6x₂+ 10x₃+ 13x₄)= 0 - (-60x₁- 70x₂ - 120x₃- 130x₄)

Систему (5.22) можно представить в виде Таблицы 5В, которую составляют следующим образом: свободные переменные, заключенные в скобки, выносят в верхнюю строку таблицы. В остальные столбцы записывают свободные члены и коэффициенты перед свободными переменными. Эта, так называемая симплекс таблица, служит основой для решения задач линейного программирования. В этой таблице переменные, являющиеся свободными, в данном случае x₁, x₂, x₃, x₄ по условию равны 0. Поскольку свободные переменные равны 0, то из системы (5.22) видно, что базисные переменные y₁, y₂, y₃, а также целевая функция F, которую записывают снизу, равны свободным членам. Значит y₁=16, y₂=110, y₃=100, F=0.

Таблица 4.2

Величина	Свободный	Свободные переменные
	член	х₁	х₂	х₃	х₄
Базисные переменные:
y₁	16	1	1	1	1
y₂	110	6	5	4	3
y₃	100	4	6	10	13
Индексная строка (F)	0	- 60	- 70	- 120	- 130

Напомним, что в системе общее число переменных N=n+m, где n - число основных переменных, m - число дополнительных переменных. Все переменные можно подразделить с одной стороны на основные и дополнительные, а с другой - на базисные и свободные.

Свободными переменными будем называть такие, которые равны 0.

Из теории известно, что n переменных в допустимом решении должны быть равны 0, т.е. столько же переменных, сколько и основных. Однако, из этого ни в коей мере не следует, что нулю равны все основные переменные. Если из общего числа переменных N=n+m будут свободными n переменных, то очевидно, что m переменных будут базисными, т.е. не равны нулю. С учетом введенных терминов можно сказать, что целью решения задачи ЛП является нахождение базисных и свободных переменных.

Для задачи ЛП, записанной в виде симплекс таблицы, можно сформулировать признаки допустимого и оптимального решений. Решение является допустимым, если в симплекс таблице в столбце свободных членов все значения, относящиеся к базисным переменным будут неотрицательными. Оптимальное значение, как мы знаем, может либо минимизировать, либо максимизировать значение целевой функции. В связи с этим, для оптимальных значений есть два признака: один для случая минимизации целевой функции, другой - для случая максимизации.

Целевая функция имеет минимальное значение, если, во-первых, решение является допустимым, т.е. свободные члены будут неотрицательными, а во=вторых, все элементы в строке целевой функции (свободный член не рассматривается) будут неположительными. При этом целевая функция равна свободному члену. Таким образом можно сделать вывод, что в Табл. 5В получено оптимальное решение нашей задачи для случая минимизации целевой функции.

Действительно, если x₁= x₂= x₃= x₄= 0 мы никакой продукции не выпускаем и при этом прибыль F=0. Дополнительные переменные y₁, y₂, y₃, показывающие объем неиспользованных ресурсов, равны, соответственно: 16, 110,100, т.е. тому ресурсу, который имеется в наличии. В самом деле, мы ничего не выпускаем, но не тратим ресурсы. Следовательно, данные в Табл. 5В соответствуют такой вершине ОДР, в которой целевая функция принимает минимальное значение.

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

Поскольку Табл. 4.2 не удовлетворяет данному признаку, то необходимо перейти к другой вершине ОДР. Переход от одной вершины к другой, производится по определенному алгоритму симплекс-метода, который заключается в обмене переменных.

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

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

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

Таким образом, ведущим столбцом будет столбец х₃, а ведущей строкой будет строка y₃. На пересечении ведущего столбца и ведущей строки будет разрешающий элемент. В нашем случае это число 10. Таким образом, ведущей строкой будет строка y₃, обозначим ее через A_r. Ведущим столбцом будет столбец х₄, обозначим его через A_k, и, следовательно, разрешающий элемент A_rk = 10.

Теперь приступаем к составлению второй таблицы. В этой таблице элементы направляющей строки равны этому элементу ведущей строки, деленному на разрешающий элемент, и находятся в соотношении:

Элементы любой другой строки j находят из соотношения:

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

После этих преобразований новая симплексная таблица после первой итерации примет вид Таблицы 4.3.

Таблица 4.3

Величина	Свободный	Свободные переменные
	член	х₁	х₂	х₃	х₄
Базисные переменные:
y₁	6	0,6	0,4	0	- 0,3
y₂	70	4,4	2,6	0	- 2,2
y₃	10	0,4	0,6	1	1,3
Индексная строка (F)	1200	2	0	26

В Табл. 4.3 в индексной строке только один отрицательный элемент, следовательно, ведущим столбцом будет х₁, а ведущую строку определяем по наименьшему отношению:

Т. о. Ведущей строкой будет y₁, а разрешающим элементом число 0,6.

Применительно к нашей задаче последняя симплекс-таблица, полученная после второй итерации, будет иметь вид:

Таблица 4.4

Величина	Свободный	Свободные переменные
	член	х₁	х₂	х₃	х₄
Базисные переменные:
х₁	10	5/3	2/3	- 1/6	- 1/2
y₂	26	-22/3	- 1/3	1/3	0
х₃	6	- 2/3	1/3	1/6	3/2
Индексная строка (F)	1320	20	10	10	20

Из этой таблицы видно, что в столбце свободных членов все элементы положительные. Значит решение является допустимым. В строке целевой функции все элементы тоже положительные. Следовательно, это решение оптимальное и максимизирует целевую функцию. При этом оптимальным планом будут следующие величины: х₁^*=10, х₃^*=6 (значит, они - базисные) и х₂^*=х₄^*=0 (т.к. они свободные). При этом целевая функция F=1320.

Вот результат решения задачи. Однако, с помощью симплекс-таблицы можно узнать еще много полезных сведений. Так их этой же таблицы видим, что свободные переменные y₁=y₃=0, а базисная переменная y₂=26. А это значит, что в оптимальном плане резервы трудовых ресурсов и оборудования равны нулю. Иными словами, эти ресурсы используются полностью. Вместе с тем резерв ресурсов сырья y₂=26, что свидетельствует о том, что имеются излишки сырья. Вот какие полезные сведения можно получить из окончательной симплекс-таблицы.

4.6 Решение транспортных задач

В качестве примера приведем решение транспортной задачи ЛП с помощью таблицы. Транспортная таблица состоит из m строк и n столбцов. В правом верхнем углу каждой клетки будем ставить стоимость С_ij перевозки единицы груза из A_i в B_j, а в центр клетки поместим перевозку X_ij.

Таблица 4.1

	ПН	В₁	В₂	В₃	В₄	В₅	Запасы а_i
ПО
A₁		¹³	⁷	¹⁴	⁷	⁵	30
A₂		¹¹	⁸	¹²	⁶	⁸	48
A₃		⁶	¹⁰	¹⁰	⁸	¹¹	20
A₄		¹⁴	⁸	¹⁰	¹⁰	¹⁵	30
Заявки b_j		18	27	42	26	15	128

Таблица 4.2

	ПН	В₁	В₂	В₃	В₄	В₅	Запасы а_i
ПО
A₁		18 ¹³	12 ⁷	¹⁴	⁷	⁵	30
A₂		¹¹	15 ⁸	33¹²	⁶	⁸	48
A₃		⁶	¹⁰	9¹⁰	11⁸	¹¹	20
A₄		¹⁴	⁸	¹⁰	15¹⁰	15¹⁵	30
Заявки b_j		18	27	42	26	15	128

Составим опорный план. Можно применить метод «северо-западного угля». Пусть пункт В₁ подал заявки на 18 единиц груза. Удовлетворим ее из запасов А₁. После этого в нем остается еще 30-18=12 единиц груза. Отдадим их пункту В₂. Но заявка этого пункта еще не удовлетворена. Выделим остаток 27-12 из запасов А₂ и т.д. рассуждая аналогичным образом, составим таблицу 4.3. Полученный план перевозок является опорным, но вряд ли он является оптимальным в смысле стоимости перевозок.

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

Таким образом, задача ЛП на геометрическом языке может быть сформулирована так: среди прямых уровня функции цели  найти опорную по отношению к ОДР и притом так, чтобы вся область лежала со стороны больших значений . Наш план - не оптимальный. Сразу видно, что его можно улучшить, если произвести в нем «циклическую перестановку», уменьшив перевозки в «дорогой» клетке (2.3) со стоимостью 12. но зато, увеличив перевозки в «дешевой» клетке (2.4) со стоимостью 6. чтобы план оставался опорным, мы должны при этом сделать одну из свободных клеток базисной, а одну из базисных - свободной.

Сколько единиц груза можем мы перенести по циклу следующему циклу: (2.4)(3.4)(3.3)(2.3), увеличивая перевозки в нечетных вершинах цикла и уменьшая в четных? Очевидно, не больше 11 единиц (иначе перевозки в клетке (3.4) стали бы отрицательными). Также очевидно, что в результате циклического переноса допустимый план остается допустимым - баланс заявок и запасов не нарушается. Произведем перенос и запишем улучшенный план в таблицу 4.3.

Таблица 4.3

	ПН	В₁	В₂	В₃	В₄	В₅	Запасы а_i
ПО
A₁		18 ¹³	12 ⁷	¹⁴	⁷	⁵	30
A₂		¹¹	15 ⁸	33¹²	11⁶	⁸	48
A₃		⁶	¹⁰	20¹⁰	⁸	¹¹	20
A₄		¹⁴	⁸	¹⁰	15¹⁰	15¹⁵	30
Заявки b_j		18	27	42	26	15	128

Таблица 4.4

	ПН	В₁	В₂	В₃	В₄	В₅	Запасы а_i
ПО
A₁		^-3 ¹³	12 ⁷	¹⁴	⁷	⁺15⁵	30
A₂		¹¹	15 ⁸	22¹²	11⁶	⁸	48
A₃		⁶	¹⁰	20¹⁰	⁸	¹¹	20
A₄		⁺15¹⁴	⁸	¹⁰	15¹⁰	^{- 15}	30
Заявки b_j		18	27	42	26	15	128

Посмотрим, что мы сэкономили. Общая стоимость плана в табл. 4.2 равна:

‡₁=18´13+12´7+15´8+33´12+9´10+11´8+15´10+15´15=1287.

Общая стоимость плана табл. 4.1 равна:

‡₂=18´13+12´7+15´8+22´12+11´6+20´10+15´10+15´15=1243.

Таким образом, нам удалось уменьшить стоимость перевозок на 44 единицы.

Действительно алгебраическая сумма стоимостей, стоящих в вершинах цикла со знаком «+», если перевозки в этой вершине увеличиваются, и со знаком «-», если уменьшаются (так называемая «цена цикла»). в данном случае равна 6-8+10-12=-4. Значит, при переносе одной величины груза по этому циклу стоимость уменьшается на 4. А мы перенесем 11 единиц. Следовательно, цена цикла 4. 11=44.

Попробуем еще раз улучшить план табл. 4.3 с помощью цикла (табл. 4.4) с ценой: 5-15+14-13=9. Перебрасывая 15 единиц груза, сокращаем стоимость на: 9. 15=135.

5. Практическая робота

.1 Графическое решение задачи распределения ресурсов

Ресурсы	Продукция		Наличие
	П1	П2
трудовые	1	2	12
материальные	2	3	17
финансовые	3	4	20
выпуск	1	1	-
прибыль	2	3	-
план	x₁	x₂	-
Целевая функция	а) x₁+ x₂® max б) прибыль ® max

Найдем наибольшее значение линейной функции графическим методом.

L = 2 x₁+ 3 x₂

при следующих ограничениях

x₁+ 2 x₂12

x₁+ 3 x₂17

x₁+ 4 x₂20

Решение:

В первую очередь, найдем область допустимых значений, т.е. точки x₁ и x₂, которые удовлетворяют системе ограничений. По условию задачи x₁ 0, x₂ 0, т.е. мы рассматриваем только те точки, которые принадлежат первой четверти.

Шаг 1

Рассмотрим неравенство 1 системы ограничений.

₁ + 2 x₂ ≤ 12

) Построим прямую.

Заменим знак неравенства на знак равенства.

₁ + 2 x₂ = 12

Преобразуем уравнение следующим образом.

x₁+x₂= 12

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

На оси X₁ рисуем точку с координатой 12.

На оси X₂ рисуем точку с координатой 6.

Соединяем полученные точки и получаем необходимую прямую.

) Какие точки нас интересуют?

x₁+ 2 x₂ 12

x₂- x₁+ 12₂-1/2 x₁+ 6

Знак неравенства меньше или равно нуля, следовательно, нас интересуют точки лежащие ниже построенной нами прямой.

) Объединим полученную полуплоскость с ранее найденными ограничениями, получим рисунок, приведенный справа.

Область допустимых значений выделена штриховкой.

Точки принадлежащие области допустимых значений:

А (0, 0)

В (12, 0)

С (0, 6)

Шаг 2

Рассмотрим неравенство 2 системы ограничений.

x₁ + 3 x₂ ≤ 17

) Построим прямую.

Заменим знак неравенства на знак равенства.

x₁ + 3 x₂ = 17

Преобразуем уравнение следующим образом.

x₁+x₂= 17

/2 1/3

Каждый член уравнения разделим на 17.

x₁+x₂= 1

17/2 17/3

На оси X₁ рисуем точку с координатой 17/2.

На оси X₂ рисуем точку с координатой 17/3.

Соединяем полученные точки и получаем необходимую прямую.

) Какие точки нас интересуют?

2x₁+ 3 x₂17

x₂-2 x₁+ 17₂-2/3 x₁+ 17/3

) Объединим полученную полуплоскость с ранее найденными ограничениями, получим рисунок, приведенный справа.

Область допустимых значений выделена штриховкой.

Точки принадлежащие области допустимых значений:

А (0, 0)

В (17/2, 0)

С (0, 17/3)

Шаг 3

Рассмотрим неравенство 3 системы ограничений.

x₁ + 4 x₂ ≤ 20

) Построим прямую.

Заменим знак неравенства на знак равенства.

x₁ + 4 x₂ = 20

Преобразуем уравнение следующим образом.

x₁+x₂= 20

/3 1/4

Каждый член уравнения разделим на 20.

x₁+x₂= 1

/35

На оси X₁ рисуем точку с координатой 20/3.

На оси X₂ рисуем точку с координатой 5.

Соединяем полученные точки и получаем необходимую прямую.

) Какие точки нас интересуют?

3 x₁+ 4 x₂ 20

x₂-3 x₁+ 20₂-3/4 x₁+ 5

) Объединим полученную полуплоскость с ранее найденными ограничениями, получим рисунок, приведенный справа.

Область допустимых значений выделена штриховкой.

Точки принадлежащие области допустимых значений:

А (0, 0)

В (20/3, 0)

С (0, 5)

Шаг 4

Вернемся к нашей исходной функции L.

L = 2 x₁+ 3 x₂

Допустим значение функции L равно 1 (абсолютно произвольно выбранное число), тогда

= 2 x₁+ 3 x₂

Построим вектор ON. Причем очевидно, что значение функции будет возрастать при перемещении прямой в направлении вектора ON.

Диапазон перемещения прямой НЕ от точки O до точки N, а именно, в направлении от точки O к точке N.

Будем перемещать прямую, перпендикулярную вектору ON, до тех пор, пока она полностью не пройдет область допустимых решений.

В нашем случае, касание прямой, перед выходом из области допустимых решений, произойдет в точке G (0, 5). В данной точке значение функции будет наибольшим.

Ответ:

Наибольшее значение функции достигает при:

Х1 = 0

Х2 = 5

Значение функции: L = 15

.2 Транспортная задача

Поставщик	Потребитель					Запасы А
	B ₁	B ₂	B ₃	B ₄	B ₅
A ₁	- 10	- 7	- 14	- 8	- 10	40
A ₂	- 12	- 8	- 10	- 8	- 10	38
A ₃	- 8	- 10	- 10	- 12	- 14	42
A ₄	- 16	- 10	- 8	- 12	- 16	30
Заявки В	28	22	25	35	40	150

) Метод северо-западного угла.

Рассмотрим маршрут доставки от поставщика A₄ к потребителю B₅ (ячейка A₄B₅).

Запасы поставщика A₄ составляют 30 единиц продукции. Потребность потребителя B₅ составляет 30 единиц продукции. (см. таблицу).

От поставщика A₄ к потребителю B₅ будем доставлять 30 единиц продукции.

Разместим в ячейку A₄B₅ значение равное 30.

Мы полностью израсходoвали запасы поставщика A₄. Вычеркиваем строку 4 таблицы, т.е. исключаем ее из дальнейшего рассмотрения.

Поставщик	Потребитель					Запасы А
	B ₁	B ₂	B ₃	B ₄	B ₅
A ₁	28 10	12 7	- 14	- 8	- 10	40
A ₂	- 12	10 8	25 10	3 8	- 10	38
A ₃	- 8	- 10	- 10	32 12	10 14	42
A ₄	- 16	- 10	- 8	- 12	30 16	30
Запасы В	28	22	25	35	40	150

Заполненные нами ячейки будем называть базисными, остальные - свободными.

Для решения задачи методом потенциалов, количество базисных ячеек (задействованных маршрутов) должно равняться m + n - 1, где m - количество строк в таблице, n - количество столбцов в таблице.

Количество базисных ячеек (задействованных маршрутов) равно 8, что и требовалось.

Мы нашли начальное решение, т.е. израсходовали все запасы поставщиков и удовлетворили все потребности потребителей.₀ = 10 * 28 + 7 * 12 + 8 * 10 + 10 * 25 + 8 * 3 + 12 * 32 + 14 * 10 + 16 * 30 = 1722 ден. ед.

Общие затраты на доставку всей продукции, для начального решения, составляют 1722 ден. ед.

) Метод наименьших стоимостей.

Минимальный элемент матрицы тарифов находится в ячейке A₄B₅ и равен 16, т.е. из незадействованных маршрутов, маршрут доставки продукции от поставщика A₄ к потребителю B₅ наиболее рентабельный.

Запасы поставщика A₄ составляют 5 единиц продукции. Потребность потребителя B₅ составляет 5 единиц продукции. (см. таблицу).

От поставщика A₄ к потребителю B₅ будем доставлять 5 единиц продукции.

Разместим в ячейку A₄B₅ значение равное 5

Поставщик	Потребитель					Запас
	B ₁	B ₂	B ₃	B ₄	B ₅
A ₁	- 10	22 7	- 14	18 8	- 10	40
A ₂	- 12	- 8	- 10	17 8	21 10	38
A ₃	28 8	- 10	- 10	- 12	14 14	42
A ₄	- 16	- 10	25 8	- 12	5 16	30
Потребность	28	22	25	35	40	150

Заполненные нами ячейки будем называть базисными, остальные - свободными.

Количество базисных ячеек (задействованных маршрутов) равно 8, что и требовалось.

Мы нашли начальное решение, т.е. израсходовали все запасы поставщиков и удовлетворили все потребности потребителей.₀ = 7 * 22 + 8 * 18 + 8 * 17 + 10 * 21 + 8 * 28 + 14 * 14 + 8 * 25 + 16 * 5 = 1344 ден. ед.

Общие затраты на доставку всей продукции, для начального решения, составляют 1344 ден. ед.

Список литературы

1. Мину М. Математическое программирование. Теория и алгоритмы. М. 2004

Орехов Н.А., Левин А.Г., Горбунов Е.А. «Математические методы и модели в экономике». Учебное пособие для вузов / Под ред. проф. Н.А. Орехова - М.: ЮНИТИ-ДАНА, 200

. Самаров К.Л., Шапкин А.С. «Задачи с решениями по высшей математике и математическим методам в экономике»: Учебное пособие - М.: Издательско-торговая корпорация «Дашков и Ко», 2007

. Мацнев А.П., Якушин А.А. «Экономико-математические методы и модели». Учебное пособие для студентов. - М: Новый Центр, 2006.

Формирование и решение линейных задач

Формирование и решение линейных задач

Похожие работы на - Формирование и решение линейных задач