
- •Лабораторный практикум
- •По дисциплине «Математические методы в экономике»
- •Часть 1.Линейное и дискретное программирование
- •Для подготовки студентов по направлениям 080500 - менеджмент и 080100 - экономика
- •Лабораторный практикум по дисциплине «Математические методы в экономике» . Часть 1.Линейное и дискретное программирование. Уч. Пособие. М.: фгоу впо ргау - мсха им. К.А. Тимирязева, - 124 с.
- •Содержание
- •Введение
- •Лабораторная работа № 1 «Запись условий задач линейного программирования»
- •Соответствие решений канонической, исходной и однородной форм задач линейного программирования
- •Алгоритм перехода от исходной формы к канонической форме
- •Лабораторная работа № 2 «Переход от задачи в канонической форме к задаче в однородной форме»
- •Эквивалентность различных форм задач линейного программирования
- •Лабораторная работа № 3 «Графический метод решения задач линейного программирования»
- •2 Область допустимых решений и допустимое решение
- •3 Линия уровня
- •4 Вектор – градиент
- •5 Оптимальное решение
- •6 Возможные варианты графического решения для двух переменных
- •Лабораторная работа № 4-5 «Решение задач линейного программирования на основе фундаментальной теоремы»
- •1 Теоремы о выпуклых множествах и области допустимых решений
- •2 Фундаментальная теорема
- •3 Применение фундаментальной теоремы для неограниченной области допустимых решений
- •2 Случай
- •4 Обобщенная фундаментальная теорема линейного программирования
- •5 Решение задач на основе фундаментальной теоремы
- •Лабораторная работа № 6 «Получение исходного опорного решения»
- •1 Общие сведения
- •2 Переход от однородной формы или исходной формы общей задачи линейного программирования к канонической форме
- •3 Переход от исходного базисного решения к исходному опорному решению
- •4 Получение исходного опорного решения м-задачи.
- •Правила перехода к м-задаче от исходной формы (основной) задачи
- •Контрольные вопросы
- •Лабораторная работа №7 - 12 «Решение общей задачи линейного программирования симплексным методом»
- •Алгоритм симплексного метода в полных таблицах при решении задач на максимум целевой функции.
- •Глоссарий основных понятий
- •Рекомендуемая литература
4 Получение исходного опорного решения м-задачи.
В задаче, записанной в канонической форме с неотрицательными свободными членами, но не имеющей в каждом ограничении базисную переменную, для получения исходного опорного решения разрешающую строку можно выбирать по минимальному симплексному отношению, то есть отношению неотрицательных свободных членов к строго положительным элементам разрешающего столбца, причём разрешающий столбец выбирается произвольно среди столбцов свободных переменных. При больших размерностях задачи эти преобразования громоздки, поэтому к таким задачам применяют правила перехода к М-задаче, всегда имеющей исходное опорное решение.
Правила перехода к м-задаче от исходной формы (основной) задачи
В ограничения типа меньше или равно (≤) вводим неотрицательные дополнительные переменные с коэффициентом плюс единица.
В ограничения типа больше или равно (≥) вводим неотрицательные дополнительные переменные с коэффициентом минус единица.
В ограничения, не содержащие базисные переменные, вводим неотрицательные искусственные базисные переменные Yi.
В целевую функцию дополнительные переменные вводятся с нулевым коэффициентом.
В целевую функцию искусственные базисные переменные Yi при решении задачи на максимум целевой функции вводятся с коэффициентом "-М", а при решении задачи на минимум целевой функции вводятся с коэффициентом "+М", где "М" большое положительное число: М ≥ 104.
Правила перехода к М-задаче от канонической формы (основной) задачи без единичного базиса с неотрицательными свободными членами
В ограничения, не содержащие базисные переменные, вводим неотрицательные искусственные базисные переменные Yi.
В целевую функцию искусственные базисные переменные Yi при решении задачи на максимум целевой функции вводятся с коэффициентом "-М", а при решении задачи на минимум целевой функции вводятся с коэффициентом "+М", где "М" большое положительное число: М ≥ 104.
Таким образом, получаем следующую М-задачу:
Задача приведена к единичному базису. Каждое ограничение имеет базисные переменные. Это соответственно: Xn+1, Xn+2,..., Xn+k, Y1, ... , Y m-k . Свободным переменным присвоим значение ноль, тогда базисные переменные равны неотрицательным свободным членам. Задача имеет следующее опорное решение
X(0,0, ... , 0, a10, a20, ... , ak0, 0, ... , 0, a(k+1)0, ... , al0, a(l+1)0, ... , am0).
Значение целевой функции для исходного опорного решения М-задачи имеет вид: Z= C0 -Ma(k+1)0- ... -Mal0-Ma(l+1)0- ...-Mam0.
Пример выполнения работы:
Получить исходное опорное решение в примерах 1-3.
Пример 1.
5X1+3X2-X3 10
6X1-5X2+2X3 15
X10, X20, X30
maxZ=4X1+3X2+2X3
Перейдём к канонической форме, получим
5X1+3X2-X3 +X4= 10
6X1-5X2+2X3 +X5=15
Xj0, j=15
maxZ=4X1+3X2+2X3+0X4+0X5
Переменные X4 и X5 - базисные переменные, X1, X2 и X3 - свободные переменные, равные нулю, поэтому задача приведена к исходному опорному решению.
Исходное опорное решение имеет вид Xисх.опорное(0,0,0,10,15).
Пример2. 3X1+2X2 6
X1-X2 -9 .
-X1-2X2 -4
Xj0, j=12
maxZ=4X1+3X2
Перейдём к канонической форме.
3X1+2X2+X3 = 6
X1-X2 +X4= -9 - выделенная строка.
-X1-2X2 +X5= -4
Xj0, j=15
maxZ=4X1+3X2+0X3+0X4+0X5
Получили каноническую форму, каждое ограничение которой содержит базисную переменную, то есть система приведена к единичному базису, но свободные члены произвольные (и положительные, и отрицательные числа). Поэтому система приведена к исходному базисному решению Xисх.базисное(0,0,6,-9,-4). В выделенной строке есть отрицательный коэффициент при свободной переменной, поэтому исходное базисное решение записываем в таблицу Гаусса.
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
|
X3 |
3 |
2 |
1 |
0 |
0 |
6 |
6/2=3 |
|
X4 |
1 |
-1 |
0 |
1 |
0 |
-9 |
-9/ (-1)=9 |
выдел.строка |
X5 |
-1 |
-2 |
0 |
0 |
1 |
-4 |
-4/ (-2)=2 |
I |
|
|
j |
|
|
|
|
|
|
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
|
X3 |
2 |
0 |
1 |
0 |
1 |
2 |
2/1=2 |
I |
X4 |
3/2 |
0 |
0 |
1 |
-1/2 |
-7 |
-7/ (-1/2)=14 |
выдел.строка |
X2 |
1/2 |
1 |
0 |
0 |
-1/2 |
2 |
- |
|
|
|
|
|
|
j |
|
|
|
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
|
X5 |
2 |
0 |
1 |
0 |
1 |
2 |
- |
|
X4 |
5/2 |
0 |
1/2 |
1 |
0 |
-6 |
- |
выдел.строка |
X2 |
3/2 |
1 |
1/2 |
0 |
0 |
3 |
- |
|
Выделенная строка с минимальным отрицательным свободным членом «-6». В выделенной строке нет ни одного отрицательного коэффициента при переменных, следовательно, система не совместна в области опорных (допустимых) решений. Можно это проверить графическим методом.
Пример 3. X1+2X2 8
-3X1-2X2 -6 .
-X1+0,5X2 -3
Xj0, j=12
maxZ=X1+6X2
Перейдём к канонической форме.
X1+2X2+X3 = 8
-3X1-2X2 +X4= -6 - выделенная строка.
-X1-0,5X2 +X5= -3
Xj0, j=15
maxZ=X1+6X2+0X3+0X4+0X5
Получили каноническую форму, каждое ограничение которой содержит базисную переменную, то есть система приведена к единичному базису, но свободные члены произвольные (и положительные, и отрицательные числа). Поэтому система приведена к исходному базисному решению Xисх.базисное(0,0,8,-6,-3). В выделенной строке есть отрицательные коэффициенты при свободных переменных, поэтому исходное базисное решение записываем в таблицу Гаусса.
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
|
X3 |
1 |
2 |
1 |
0 |
0 |
8 |
8/1=4 |
|
X4 |
-3 |
-2 |
0 |
1 |
0 |
-6 |
-6/ (-3)=2 |
выдел.строка I |
X5 |
-1 |
0,5 |
0 |
0 |
1 |
-3 |
-3/ (-1)=3 |
|
|
j |
|
|
|
|
|
|
|
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
|
X3 |
0 |
4/.3 |
1 |
1/3 |
0 |
6 |
6/(1/3)=18 |
|
X1 |
1 |
2/3 |
0 |
-1/3 |
0 |
2 |
- |
|
X5 |
0 |
7/6 |
0 |
-1/3 |
1 |
-1 |
-1/(-1/3)=3 |
выдел.строка I |
|
|
|
|
j |
|
|
|
|
Баз\ пер |
X1 |
X2 |
X3 |
X4 |
X5 |
a i0 |
РСО |
X3 |
0 |
5/2 |
1 |
0 |
1 |
5 |
- |
X1 |
1 |
-1/2 |
0 |
0 |
-1 |
3 |
- |
X4 |
0 |
-7/2 |
0 |
1 |
-3 |
3 |
- |
Система приведена к единичному базису, все свободные члены неотрицательные, следовательно, система совместна в области опорных (допустимых) решений. Получили исходное опорное решение X(3,0,5,3,0), в котором свободные переменные X2 и X5 равны нулю, а базисные переменные неотрицательные. Можно это проверить графическим методом, найдя на графике вершину с координатами (3, 0), однозначно соответствующую исходному опорному решению.
Пример 3а.
Получить исходное опорное решение М задачи.
X1+2X2 8
-3X1-2X2 -6 .
-X1+0,5X2 -3
Xj0, j=12
maxZ=X1+6X2
Задача дана в однородной форме, но не все свободные члены неотрицательные. Приведём задачу к исходной форме с неотрицательными свободными членами.
X1+2X2 8
-3X1+2X2 6 .
X1-0,5X2 3
Xj0, j=12
maxZ=X1+6X2
Перейдём к записи задачи в канонической форме. Получим
X1+2X2+X3 = 8 , X3 - базисная переменная;
3X1+2X2 -X4 = 6 .
X1-0,5X2 -X5= 3
Xj0, j=15
maxZ=X1+6X2+0X3+0X4+0X5
Получили каноническую форму, но не каждое ограничение её содержит базисную переменную, то есть система не приведена к единичному базису, причём свободные члены неотрицательные числа. Поэтому систему приведем к исходному опорному решению М задачи, введя во второе и третье ограничения искусственные неотрицательные базисные переменные Yi. Получим
X1+2X2+X3 = 8 , X3 - базисная переменная;
3X1+2X2 -X4 +Y1 = 6, Y1 - искусственная базисная переменная;
X1-0,5X2 -X5 +Y2= 3 , Y2 - искусственная базисная переменная;
Xj0, j=15
maxZ=X1+6X2+0X3+0X4+0X5 -M*Y1 –M*Y2
Система приведена к единичному базису с неотрицательными свободными членами. Пусть свободные переменные равны нулю, тогда исходное опорное решение М задачи имеет вид Хисходн. опорное решение М задачи(0,0,8,0,0,6,3),
Z(Хисходное опорное решение М задачи)= -M*6 –M*3= - 9М. Оно не совпадает с исходным опорным решением основной задачи.
Список индивидуальных данных
1). Получить исходное опорное решение с помощью расширенных симплексных отношений или установить несовместность системы ограничений в области допустимых (опорных) решений.
2). Получить исходное опорное решение М задачи и подсчитать Z(Xисх.опор.М задачи).
3) Проверить результат графическим методом.
Z=C1*X1+C2*X2+C0
ai1* Х1+ ai2*Х2( ai0, I=1m
Х1 0, Х2 0
№ вари- анта |
Х1 |
Х2 |
Тип ограничения |
|
|
№ варианта |
Х1 |
Х2 |
Тип ограничения |
|
|
ai1 |
ai2 |
|
ai0 |
|
|
ai1 |
ai2 |
|
ai0 |
1 |
|
|
|
|
|
2 |
|
|
|
|
maxZ |
-1 |
-2 |
|
-3 |
|
minZ |
1 |
2 |
|
5 |
|
1 |
1 |
|
4 |
|
|
-1 |
1 |
|
4 |
|
2 |
1 |
|
2 |
|
|
2 |
-1 |
|
2 |
|
1 |
2 |
|
2 |
|
|
-1 |
2 |
|
2 |
3 |
|
|
|
|
|
4 |
|
|
|
|
maxZ |
-1 |
-2 |
|
4 |
|
minZ |
1 |
1 |
|
6 |
|
1 |
1 |
|
7 |
|
|
1 |
1 |
|
7 |
|
2 |
1 |
|
2 |
|
|
-2 |
1 |
|
2 |
|
1 |
2 |
|
2 |
|
|
1 |
-2 |
|
2 |
5 |
|
|
|
|
|
6 |
|
|
|
|
minZ |
3 |
1 |
|
2 |
|
maxZ |
1 |
2 |
|
4 |
|
1 |
1 |
|
3 |
|
|
1 |
1 |
|
5 |
|
2 |
-3 |
|
6 |
|
|
-1 |
1 |
|
3 |
|
3 |
–2 |
|
6 |
|
|
1 |
-1 |
|
3 |
7 |
|
|
|
|
|
8 |
|
|
|
|
maxZ |
-1 |
2 |
|
-5 |
|
maxZ |
1 |
-3 |
|
-9 |
|
1 |
1 |
|
5 |
|
|
1 |
1 |
|
6 |
|
2 |
-1 |
|
6 |
|
|
2 |
-1 |
|
4 |
|
-1 |
2 |
|
6 |
|
|
-1 |
2 |
|
4 |
9 |
|
|
|
|
|
10 |
|
|
|
|
maxZ |
2 |
-4 |
|
3 |
|
minZ |
-2 |
-2 |
|
3 |
|
5 |
2 |
|
10 |
|
|
2 |
5 |
|
10 |
|
2 |
-1 |
|
4 |
|
|
2 |
-3 |
|
6 |
|
-1 |
2 |
|
4 |
|
|
-3 |
2 |
|
6 |
11 |
|
|
|
|
|
12 |
|
|
|
|
maxZ |
3 |
6 |
|
7 |
|
minZ |
-3 |
-3 |
|
-7 |
|
3 |
1 |
|
6 |
|
|
3 |
1 |
|
6 |
|
2 |
-1 |
|
10 |
|
|
2 |
-1 |
|
5 |
|
-1 |
2 |
|
10 |
|
|
-1 |
2 |
|
5 |
13 |
|
|
|
|
|
14 |
|
|
|
|
maxZ |
-3 |
4 |
|
7 |
|
minZ |
-4 |
2 |
|
-7 |
|
-1 |
-2 |
|
-6 |
|
|
-1 |
-1 |
|
-4 |
|
2 |
-3 |
|
-6 |
|
|
2 |
1 |
|
2 |
|
1 |
-2 |
|
2 |
|
|
1 |
2 |
|
2 |
15 |
|
|
|
|
|
16 |
|
|
|
|
minZ |
-3 |
1 |
|
2 |
|
maxZ |
1 |
-2 |
|
4 |
|
1 |
1 |
|
-3 |
|
|
1 |
1 |
|
-5 |
|
2 |
-3 |
|
6 |
|
|
-1 |
1 |
|
3 |
|
3 |
–2 |
|
12 |
|
|
1 |
-1 |
|
8 |
17 |
|
|
|
|
|
18 |
|
|
|
|
maxZ |
-1 |
2 |
|
-5 |
|
maxZ |
-1 |
-3 |
|
9 |
|
1 |
1 |
|
5 |
|
|
1 |
1 |
|
6 |
|
2 |
-1 |
|
6 |
|
|
2 |
-1 |
|
4 |
|
-1 |
2 |
|
6 |
|
|
-1 |
2 |
|
4 |
19 |
|
|
|
|
|
20 |
|
|
|
|
maxZ |
2 |
5 |
|
3 |
|
minZ |
-5 |
-2 |
|
-3 |
|
5 |
2 |
|
10 |
|
|
2 |
5 |
|
10 |
|
2 |
-1 |
|
-4 |
|
|
2 |
-3 |
|
6 |
|
-1 |
2 |
|
14 |
|
|
-3 |
2 |
|
-6 |
21 |
|
|
|
|
|
22 |
|
|
|
|
maxZ |
3 |
-4 |
|
8 |
|
minZ |
-3 |
4 |
|
-8 |
|
3 |
1 |
|
6 |
|
|
3 |
1 |
|
6 |
|
2 |
-1 |
|
1 |
|
|
2 |
-1 |
|
15 |
|
-1 |
2 |
|
1 |
|
|
-1 |
2 |
|
15 |
23 |
|
|
|
|
|
24 |
|
|
|
|
maxZ |
-1 |
2 |
|
7 |
|
minZ |
-8 |
2 |
|
-7 |
|
-1 |
-2 |
|
6 |
|
|
-1 |
1 |
|
-4 |
|
2 |
-3 |
|
-6 |
|
|
2 |
-1 |
|
-2 |
|
3 |
-2 |
|
12 |
|
|
1 |
-2 |
|
2 |
25 |
|
|
|
|
|
26 |
|
|
|
|
minZ |
3 |
-1 |
|
12 |
|
maxZ |
-1 |
2 |
|
14 |
|
1 |
1 |
|
3 |
|
|
1 |
1 |
|
6 |
|
2 |
-3 |
|
12 |
|
|
-1 |
1 |
|
3 |
|
3 |
–2 |
|
6 |
|
|
1 |
-1 |
|
3 |
27 |
|
|
|
|
|
28 |
|
|
|
|
maxZ |
1 |
-2 |
|
-5 |
|
maxZ |
-1 |
-3 |
|
9 |
|
1 |
1 |
|
5 |
|
|
1 |
1 |
|
6 |
|
2 |
-1 |
|
6 |
|
|
-2 |
-1 |
|
4 |
|
-1 |
-2 |
|
6 |
|
|
-1 |
-2 |
|
4 |
29 |
|
|
|
|
|
30 |
|
|
|
|
maxZ |
2 |
3 |
|
3 |
|
minZ |
-2 |
-3 |
|
-3 |
|
5 |
2 |
|
10 |
|
|
2 |
5 |
|
10 |
|
-2 |
-1 |
|
4 |
|
|
-2 |
-3 |
|
6 |
|
-1 |
-2 |
|
4 |
|
|
-3 |
-2 |
|
6 |
31 |
|
|
|
|
|
32 |
|
|
|
|
maxZ |
3 |
5 |
|
7 |
|
minZ |
-5 |
-2 |
|
-7 |
|
3 |
1 |
|
6 |
|
|
3 |
1 |
|
6 |
|
-2 |
-1 |
|
10 |
|
|
-2 |
-1 |
|
5 |
|
-1 |
-2 |
|
10 |
|
|
-1 |
-2 |
|
5 |
33 |
|
|
|
|
|
34 |
|
|
|
|
maxZ |
-3 |
-2 |
|
2 |
|
minZ |
3 |
2 |
|
-6 |
|
-1 |
-2 |
|
-6 |
|
|
-1 |
-1 |
|
-4 |
|
-2 |
-3 |
|
6 |
|
|
2 |
1 |
|
2 |
|
-1 |
-2 |
|
2 |
|
|
-1 |
-2 |
|
2 |