- •Раздел 1. Модели линейного программирования и его приложения Тема 1. Математическое моделирование экономических систем
- •1.1. Описание процесса математического моделирования.
- •1.2. Примеры построения математических моделей простейших экономических задач
- •1. Задача об использовании ресурсов (планирование производства, задача оптимального использования удобрений и т.Д.).
- •2. Оптимальное смешивание (составление рациона питания, смесей, задача о диете и т.Д.).
- •3. Задача об использовании мощностей (задача о загрузке оборудования)
- •4. Задача о раскрое материалов
- •Задача о закреплении самолетов за воздушными линиями.
- •Контрольные вопросы
- •Рекомендованная литература: [ 3, 8, 11, 12] Тема 2. Линейные векторные пространства
- •Понятие n-мерного пространства
- •Линейная зависимость векторов
- •2.3. Базис n-мерного векторного пространства
- •2.4. Система единичных векторов n-мерного векторного пространства
- •2.5. Решение системы линейных уравнений методом Жордана – Гаусса
- •Контрольные вопросы
- •Для приведения неравенства к равенству необходимо к его левой части прибавить неотрицательную величину . (2)
- •Пусть удовлетворяет уравнению (3) и неравенству (2), т.Е.
- •3.2. Выпуклые множества
- •3.3. Геометрическая интерпретация задачи линейного программирования
- •3.4. Свойства решений задачи линейного программирования
- •Контрольные вопросы
- •Примеры задач, решаемых графическим методом.
- •4.2. Графический метод решения злп с n переменными
- •Контрольные вопросы
- •5.2. Алгоритм симплекс-метода
- •5.3. Пример отыскания максимума линейной функции
- •5.4. Пример отыскания минимума линейной функции
- •Решение. Введем дополнительные неотрицательные переменные y5, y6 со знаком “-”, т.К. Неравенства имеют вид ””:
- •Симплексные таблицы
- •В общем случае достаточно воспользоваться правилом:
- •5.6. Метод искусственного базиса
- •Контрольные вопросы
- •Предположим, что некоторая организация решила закупить ресурсы s1, s2,…, Sm предприятия и необходимо установить оптимальные цены на эти ресурсы y1,y2,…,ym.
- •Аналогично стоимость всех затраченных ресурсов, идущих на изготовление единицы j-ой продукции, не может быть меньше стоимости окончательного продукта, т.Е.
- •6.2. Взаимно двойственные задачи линейного программирования и их свойства
- •6.3. Первая теорема двойственности
- •6.4. Вторая теорема двойственности
- •Рассмотренная теорема является следствием следующей теоремы.
- •6.5. Объективно обусловленные оценки и их смысл
- •Контрольные вопросы
- •7.2. Нахождение опорного плана
- •Пример 1. Исходные данные приведены в таблице 1.
- •7.3. Метод последовательного улучшения плана перевозок, цикл пересчета
- •7.4. Решение транспортной задачи методом потенциалов
- •Приложения транспортной задачи к решению некоторых экономических задач
- •Контрольные вопросы
- •Понятие цикла. Рекомендованная литература: [ 3, 5, 6, 7, 8, 11]
- •Тема 8. Элементы теории матричных игр
- •8.1. Предмет теории игр, основные понятия
- •8.2. Платежная матрица. Нижняя и верхняя цена игры
- •8.3. Решение игры в смешанных стратегиях
- •Решая эту систему, получим оптимальную стратегию:
- •8.4. Геометрическая интерпретация игры 2 х 2
- •8.5. Приведение матричной игры к задаче линейного программирования
- •Составив расширенные матрицы для задач, убеждаемся, что одна матрица получилась из другой транспонированием:
- •Контрольные вопросы
- •Графический метод решения игры mx2. Рекомендованная литература: [ 1, 2, 4, 9, 12] раздел 2. Модели нелинейного программирования
- •Тема 9. Нелинейное программирование
- •9.1. Общая задача нелинейного программирования
- •9.2. Метод множителей Лагранжа
- •9.3. Обобщёние метода множителей Лагранжа
- •9.4. Теорема Куна-Таккера
- •9.5. Модели выпуклого программирования
- •9.6. Приближенное решение задач выпуклого программирования методом кусочно-линейной аппроксимации
- •Контрольные вопросы
- •Задачі про мінімізацію розходу горючого літаком за набором висоти та швидкості. Рекомендованная литература: [ 3, 6, 7, 8, 11] Бібліографічний список
3.4. Свойства решений задачи линейного программирования
Теорема 1. Множество всех допустимых решений системы ограничений ЗЛП выпукло.
Доказательство.
Дана ЗЛП в матричном виде: Z=CXmax(min) при ограничениях AX = B, X 0 .
Пусть Х1=(х11, х21,…, хn1) и Х2=(х12, х22,…, хn2) – два допустимых решения ЗЛП, соответственно выполняются соотношения AX1 = B, X1 0; AX2 = B, X2 0.
Необходимо доказать, что выпуклая комбинация решений Х1 и Х2
Х=1Х1+2Х2, 1 0, 20, 1+2=1 также допустимое решение задачи.
Действительно, АХ=А(1Х1+2Х2)= 1АХ1+2АХ2=1В+2В=(1+2)В=В,
т.е. получаем, что Х удовлетворяет системе AX = B, следовательно является решением.
Но т.к. X1 0; X2 0, 1 0,2 0, то и Х 0, т.е. удовлетворяет и условию X 0, значит является допустимым решением ЗЛП.
Ответ на вопрос, в какой точке многогранника решений возможно оптимальное решение ЗЛП, дается в следующей фундаментальной теореме.
Теорема 2. Если ЗЛП имеет оптимальное решение, то линейная функция достигает своего минимального значения в угловой точке многогранника решений. Если линейная функция принимает минимальное значение более чем в одной угловой точке, то она принимает его в любой точке, являющейся выпуклой линейной комбинацией таких угловых точек.
Доказательство.
Предположим, что многогранник решений ограниченный, имеет конечное число угловых точек. Обозначим его через К. В двумерном пространстве К имеет вид многоугольника, изображенного на рисунке.
О
бозначим
угловые точки через
х2
Х2
Х3,
Х1, Х2,…,Хр, Х1
а
оптимальный план – через Х*.
Xi
Хi
Тогда Z(X*) Z(X) 0 Хр х1
для всех точек Х многогранника решений К. Если Х* - угловая точка, то первая часть теоремы доказана.
Предположим, что Х* не является угловой точкой. Тогда Х* на основании теоремы 2 (о представлении выпуклого многогранника) можно представить как выпуклую линейную комбинацию угловых точек К, т.е.
Х*=1Х1+2Х2+…+рХр,
i
0 ( i=1,2,…,p
),
=1.
Так как Z (X) – линейная функция, получаем
Z (X*)= Z (1Х1+2Х2+…+рХр)= 1 Z (X1)+ 2 Z (X2)+…+ p Z (Xp). (*)
В этом разложении среди значений Z(Xi) ( i=1,2,..,p ) выберем наименьшее (пусть оно соответствует угловой точке Хк ( 1кр ) ) и обозначим его через m, т.е. Z (Xк)=m. Заменим в (*) каждое значение Z (Xi) этим наименьшим значением. Тогда, так как i 0 , =1, то
Z(X*)1m+2m+…+рm= m =m.
По предположению, Х* - оптимальный план, поэтому, с одной стороны, Z(X*)m, но, с другой стороны, доказано, что Z(X*)m, значит, Z(X*)=m= Z(Xк), где Хк – угловая точка.
Итак, существует угловая точка Хк, в которой линейная функция принимает минимальное значение.
Для доказательства второй части теоремы допустим, что Z (X) принимает минимальное значение более чем в одной угловой точке, например в точках Х1, Х2,…,Хq,1<qp; тогда
Z(X1)= Z(X2)=…= Z(Xq)=m.
Пусть Х – выпуклая линейная комбинация этих угловых точек:
Х=1Х1+2Х2+…+qХq,
i
0 ( i=1,2,…,q
),
=1.
Тогда, учитывая, что Z (X) – линейная функция, получаем
Z(X)=Z(1Х1+2Х2+…+qХq)= 1Z(X1)+ 2Z(X2)+…+ qZ(Xq)=1m+2m+…+qm= =m =m,
т.е. линейная функция Z принимает минимальное значение в произвольной точке Х, являющейся выпуклой линейной комбинацией угловых точек Х1,Х2,…,Хq.
Замечание. Требование ограниченности многогранника решений в теореме является существенным, так как в случае неограниченной многогранной области, как отмечалось в теореме 2 пункта 3.2., не каждую точку такой области можно представить выпуклой комбинацией ее угловых точек.
Доказанная теорема является фундаментальной, т.к. она указывает принципиальный путь решения ЗЛП. Действительно, согласно теореме вместо исследования бесконечного множества допустимых решений для нахождения среди них оптимального решения необходимо исследовать лишь конечное число угловых точек многогранника решений.
Следующая теорема посвящена аналитическому методу нахождения угловых точек.
Теорема 3. Каждому допустимому базисному решению (опорному плану) ЗЛП соответствует угловая точка многогранника решений, и наоборот, каждой угловой точке многогранника решений соответствует допустимое базисное решение.
Доказательство.
Пусть Х= (x1, x2, ... , xm,0,…,0) - допустимое базисное решение системы ограничений A1x1+A2x2+…+Anxn= B, в которой первые m компонент – базисные переменные, а остальные n-m компонент – свободные, равные нулю в базисном решении (если это не так, то соответствующие переменные можно переменуровать).
Покажем, что Х – угловая точка многогранника решений.
Предположим противное, т.е. что точка Х не является угловой. Тогда она может быть представлена внутренней точкой отрезка, соединяющего две различные, не совпадающие с Х, точки
Х1= (x1(1), x2(1), ... , xn(1)), Х2= (x1(2), x2(2), ... , xn(2)),
другими словами, - выпуклой линейной комбинацией точек Х1 и Х2 многогранника решений, т.е. (*) Х=1Х1+2Х2, 1> 0, 2>0, 1+2=1.
(Полагаем, что 10, 20, ибо в противном случае точка Х совпадет с точкой Х1 или Х2).
Запишем векторное равенство (*) в координатной форме:
Поскольку xj(1)0, xj(2)0 (j=1,2,…,n), 1>0 и 2>0 , то из последних n-m равенств следует, что xm+1(1)=0, xm+1(2)=0,…, xn(1)=0, xn(2)=0, т.е. в решениях Х1, Х2 и Х системы уравнений A1x1+A2x2+…+Anxn= B значения n-m компонент равны в данном случае нулю:
Х1= (x1(1), x2(1), ... , xm(1),0,…,0), Х2= (x1(2), x2(2), ... , xm(2),0,…,0).
Эти компоненты можно считать значениями свободных переменных.
Поскольку Х1 и Х2 – планы, то
A1x1(1)+A2x2(1)+…+Amxm(1)= B, A1x1(2)+A2x2(2)+…+Amxm(2)= B.
Вычитая из первого соотношения второе, получаем
(x1(1)- x1(2))A1+(x2(1)- x2(2))A2+…+(xm(1)- xm(2))Am= 0.
По предположению, векторы А1, А2, ... ,Аm линейно независимы, поэтому последнее соотношение выполняется, если
x1(1)- x1(2)=0; x2(1)- x2(2)=0;…;xm(1)- xm(2)= 0.
Отсюда x1(1)=x1(2); x2(1)=x2(2);…;xm(1)=xm(2).
Итак, Х невозможно представить как выпуклую линейную комбинацию двух точек многогранника решений. Следовательно, Х – угловая точка.
Докажем обратное утверждение. Пусть Х = (x1, x2, ... , xm,0,…,0) – угловая точка многогранника решений и первые ее m координат положительны. Покажем, что Х – допустимое базисное решение (опорный план).
Предположим противное, т.е. векторы А1, А2, ... ,Аm линейно зависимы. Тогда существуют такие числа l1, l2, ... ,lk, не все равные нулю, при которых выполняется соотношение
l1А1+ l2А2+ ... + lmАm=0. (*)
Зададим некоторое число >0, умножим на него равенство (*):
l1А1+ l2А2+ ... + lmАm=0. (**)
По условию А1x1+ А2x2+ ... +Аmxm=B. (***)
Равенство (**) почленно сложим с равенством (***):
(x1 + l1 )А1+ (x2 + l2 )А2+ ... + (xm + lm )Аm=B,
Теперь равенство (**) почленно вычтем из (***), получаем:
(x1 - l1 )А1+ (x2 - l2 )А2+ ... + (xm- lm )Аm=B.
Сравнивая полученные последних два равенства с равенством (***), заключаем, что при любом исходная система ограничений А1x1+ А2x2+ ... +Аnxn=B имеет два решения:
X1=( x1 + l1; x2 + l2; …; xm + lm; 0;…;0),
X2=( x1 - l1; x2 - l2; …; xm - lm; 0;…;0).
Все хi0 (j=1,2,…,n), поэтому число можно подобрать настолько малым, что все компоненты решений Х1 и Х2 примут неотрицательные значения, тогда Х1 и Х2 будут различными допустимыми решениями ЗЛП.
При этом, как легко видеть, решение ½ (Х1 + Х2) =X (показать самостоятельно), т.е. точка Х лежит на отрезке (в данном случае на его середине), расположенном в многограннике решений, а значит является выпуклой линейной комбинацией точек Х1, Х2. А значит не является угловой точкой, что противоречит условию теоремы о том, что Х – угловая точка.
Итак, наше допущение, что система векторов А1, А2, ... ,Аm линейно зависима, привело к противоречию. Следовательно, сделанное допущение неверно и система векторов линейно независима.
Следствие 1. Так как векторы А1, А2, ... , Аn имеют размерность m, то угловая точка многогранника решений имеет не более чем m положительных компонент xi > 0 (i=1,2, ..., m).
Следствие 2. Каждой угловой точке многогранника решений соответствует km линейно независимых векторов системы А1, А2, ... , Аn.
Можно предположить, что система векторов А1, А2, ... , Аn задачи линейного программирования всегда содержит m линейно независимых векторов. Если при решении частной задачи это не очевидно, то первоначальную систему векторов дополняют m линейно независимыми векторами, затем находят решение расширенной задачи.
Следствие 3. Если ЗЛП имеет оптимальное решение, то оно совпадает, по крайней мере, с одним из ее допустимых базисных решений.
Итак, если линейная функция ЗЛП ограничена на многограннике решений, то:
- существует такая угловая точка многогранника решений, в которой линейная функция ЗЛП достигает своего оптимального значения;
- каждый опорный план (допустимое базисное решение, каждое из которых определяется системой m линейно-независимых векторов, выбираемой из системы n векторов А1, А2, ... ,Аn) соответствует угловой точке многогранника решений.
