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

ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ  УЧРЕЖДЕНИЕ

«ПРИДНЕСТРОВСКИЙ ГОСУДАРСТВЕННЫЙ  УНИВЕРСИТЕТ им. Т.Г. ШЕВЧЕНКО»

РЫБНИЦКИЙ ФИЛИАЛ

КАФЕДРА «ФИЗИКИ, МАТЕМАТИКИ И ИНФОРМАТИКИ»

 

 

 

 

Курсовая работа

по дисциплине

 

"Исследование операций"

на тему:

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

 

Выполнила:

студентка III курса

специальности “Информатика с доп. спец. английский язык” 

Нистор А.Г.

 

Проверила:

преподаватель Панченко Т.А.

 

 

 

г. Рыбница

2008 г.

 

ОГЛАВЛЕНИЕ

 

 

ВВЕДЕНИЕ

 

Линейное программирование - это наука о методах исследования и

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

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

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

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

2) Разобрать алгоритм решения ЗЛП геометрическим методом.

3) Решить поставленные задачи, используя рассмотренный метод решения задач линейного программирования.

 

I. ТЕОРЕТИЧЕСКИЙ РАЗДЕЛ

1.1 Линейное программирование

 

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

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

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

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

Математическая формулировка задачи линейного программирования

 

 

Нужно максимизировать

 

 

при условиях  при i = 0, 1, 2, . . . , m .

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

Такую задачу называют "основной" или "стандартной" в линейном программировании.

1.2 Формулировка задачи

 

Даны линейная функция Z=С₁х₁+С₂х₂+...+С_Nx_N (1.1)

и система линейных ограничений

 

a₁₁x₁ + a₂₂x₂ + ... + a_1NХ_N = b₁

a₂₁x₁ + a₂₂x₂ + ... + a_2NХ_N = b₂

. . . . . . . . . . . . . . .

a_i1x₁ + a_i2x₂ + ... + a_iNХ_N = b_i (1.2) . . . . . . . . . . . . . . .

a_M1x₁ + a_M2x₂ + ... + a_MNХ_N = b_M

x_j 0 (j = 1, 2, ... ,n) (1.3)

 

где а_ij, b_j и С_j - заданные постоянные величины.

Найти такие неотрицательные  значения х₁, х₂, ..., х_n, которые удовлетворяют системе ограничений (1.2) и доставляют линейной функции (1.1) минимальное значение.

Общая задача имеет несколько  форм записи.

Векторная форма  записи. Минимизировать линейную функцию Z = СХ при ограничениях А₁х₁ + А₂x₂ + ... + А_Nx_N = А_о, X₀ (1.4)

где С = (с₁, с₂, ..., с_N); Х = (х₁, х₂, ..., х_N); СХ - скалярное произведение; векторы A₁ = A₂ = ,..., A_N состоят соответственно из коэффициентов при неизвестных и свободных членах.

Матричная форма записи. Минимизировать линейную функцию, Z = СХ при ограничениях АХ = А₀Х₀, где С = (с₁, с₂, ..., с_N) - матрица-cтрока; А = (а_ij) - матрица системы; Х =(x_ij)- матрица-столбец, А₀ = (а_i) матрица-столбец

Запись с помощью знаков суммирования. Минимизировать линейную функцию Z = С_jх_j при ограничениях

0пределение 1. Планом или допустимым  решением задачи линейного программирования  называется Х = (х₁, х₂, ..., х_N), удовлетворяющий условиям (1.2) и (1.3).

0пределение 2. План Х = (х₁, х₂, ..., х_N) называется опорным, если векторы А (i = 1, 2, ..., N), входящие в разложение (1.4) с положительными коэффициентами х , являются линейно независимыми.

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

0пределение 3. Опорный план называется  невырожденным, если он содержит М положительных компонент, в противном случае опорный план называется вырожденным.

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

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

 

1.3 Основные понятия линейной алгебры и выпуклого анализа, применяемые в теории математического программирования

 

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

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

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

Вектор λ₁a¹ + λ₂a² + …+ λ_ma^m называется линейной комбинацией векторов а¹ а²,..., а^m с коэффициентами λ₁, λ_2, λ_m,

Система векторов линейного пространства а¹ а²,..., а^m называется линейно зависимой, если существуют такие числа λ₁, λ_2, λ_m не равные одновременно нулю, что их линейная комбинация λ₁a¹ + λ₂a² + …+ λ_ma^m равняется нулевому вектору (вектору, все компоненты которого равны нулю). В противном случае систему а¹, а²,..., а^m называют линейно независимой, т. е. линейная комбинация данных векторов может быть равна нулевому вектору только при нулевых коэффициентах λ₁, λ_2, …, λ_m

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

Линейное пространство обычно обозначают как Rⁿ, где n — его размерность.

Любое подмножество данного линейного  пространства, которое само обладает свойствами линейного пространства, называется линейным подпространством. Множество Н, получаемое сдвигом некоторого линейного подпространства L € Rⁿ на вектор a € Rⁿ: H=L+a, называется аффинным множеством (пространством). Если фундаментальным свойством любого линейного пространства или подпространства является принадлежность ему нулевого вектора, то для аффинного множества это не всегда так. На плоскости примером подпространства является прямая, проходящая через начало координат, а аффинного множества — любая прямая на плоскости. Характеристическим свойством аффинного множества является принадлежность ему любой прямой, соединяющей две любые его точки. Размерность аффинного множества совпадает с размерностью того линейного подпространства, сдвигом которого оно получено.

Если рассматривается некоторое  линейное пространство Rⁿ, то принадлежащие ему аффинные множества размерности 1 называются прямыми, а размерности (n-1)—гиперплоскостями. Так, обычная плоскость является гиперплоскостью для трехмерного геометрического пространства R³, а прямая — гиперплоскостью для плоскости R². Всякая гиперплоскость делит линейное пространство на два полупространства.

Множество V векторов (точек) линейного пространства Rⁿ называется выпуклым, если оно содержит отрезок прямой, соединяющей две его любые точки, или, другими словами, из того, что a €V и b€V , следует, что х = (1- λ) х а+ λ х b € V , где 0 ≤ λ ≤ 1.

Линейная комбинация векторов а¹, а²... а^m называется выпуклой, если λ_i ≥0, i €1:m и

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

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

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

Множество К называется конусом  с вершиной в точке x₀, если x₀ € К , и из того, что некоторая точка х принадлежит К ( х € К ), следует, что в К содержится и луч, начинающийся в х₀ и проходящий через х, т. е.

 

 

или

 

 

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

1.4 Математические основы решения задачи линейного программирования графическим способом

1.4.1 Математический аппарат

Для понимания всего  дальнейшего полезно знать и  представлять себе геометрическую интерпретацию задач линейного программирования, которую можно дать для случаев n = 2 и n = 3.

Наиболее наглядна эта  интерпретация для случая n = 2, т.е. для случая двух переменных x₁ и x₂. Пусть нам задана задача линейного программирования в стандартной форме

1.5)

 

Возьмём на плоскости  декартову систему координат  и каждой паре чисел (x₁,x₂)поставим в соответствие точку на этой плоскости.

 

 

Обратим прежде всего  внимание на ограничения x₁ ≥0 и x₂ ≥ 0. Они из всей плоскости вырезают лишь её первую четверть (см. рис. 1). Рассмотрим теперь, какие области соответствуют неравенствам вида a₁ x₁ + a₂ x₂ ≤ b. Сначала рассмотрим область, соответствующую равенству a₁ x₁ + a₂ x₂ = b. Как Вы, конечно, знаете, это прямая линия. Строить её проще всего по двум точкам.

Пусть b ≠ 0. Если взять x₁ = 0, то получится x₂ = b/a₂. Если взять x₂ = 0, то получится x₁ = b/a₁. Таким образом, на прямой лежат две точки (0, b/a₂) и (b/a₁, 0). Дальше через эти две точки можно по линейке провести прямую линию (рисунок 2).

 

 

Если же b=0, то на прямой лежит точка (0,0). Чтобы найти другую точку, можно взять любое отличное от нуля значение x₁ и вычислить соответствующее ему значение x₂.

Эта построенная прямая разбивает всю плоскость на две  полуплоскости. В одной её части a₁x₁ + a₂x₂ < b, а в другой наоборот a₁x₁ + a₂x₂ > b. Узнать, в какой полуплоскости, какой знак имеет место проще всего посмотрев, какому неравенству удовлетворяет какая-то точка плоскости, например, начало координат, т.е. точка (0,0).

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

Рассмотрим задачу ЛП в стандартной форме записи:

 

max f(X) = с₁х₁ + с₂х₂ + ... + с_пх_п (*)

 

при ограничениях

 

а₁₁х₁ + а₁₂х₂ + … + а_1nх_n ≤ b₁