Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

3376

.pdf
Скачиваний:
1
Добавлен:
08.01.2021
Размер:
542.46 Кб
Скачать

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ОБРАЗОВАНИЯ «ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ ЛЕСОТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

ИМЕНИ Г.Ф. МОРОЗОВА»

МЕТОДЫ ОПТИМАЛЬНЫХ РЕШЕНИЙ

Методические указания для самостоятельной работы студентов

по направлению подготовки 38.03.01 – Экономика

Воронеж 2016

3

УДК 512.8

Раецкая, Е. В. Методы оптимальных решений [Текст] : методические указания для самостоятельной работы студентов по направлению подготовки 38.03.01 – Экономика / Е. В. Раецкая, И.В. Сапронов, Н.М. Спирина; М-во образования и науки РФ, ФГБОУ ВО «ВГЛТУ». – Воронеж, 2016. – 29 с.

Печатается по решению учебно-методического совета ФГБОУ ВО «ВГЛТУ» (протокол № 5 от 22 апреля 2016 г.)

 

Рецензент д-р физ.-мат. наук, доцента кафедры математического анализа

ВГУ

Зубова С.П.

4

Введение

Целью изучения дисциплины «Методы оптимальных решений» является воспитание достаточно высокой математической культуры, привитие навыков современных видов математического мышления, ознакомление с математическими свойствами моделей и методами оптимизации; применению методов оптимизации, которые могут использоваться при анализе и решении широкого спектра экономических задач.

Для достижения поставленной цели, при самостоятельной работе решаются следующие задачи:

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

-закрепление теоретического материала и выработка умения самостоятельно

применять

математические методы в технических и в экономических

приложениях;

 

Врезультате самостоятельного освоения дисциплины студент должен:

-знать основные понятия, определения и методы исследования объектов с помощью теорем и формул различных разделов курса математики;

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

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

Студент по результатам освоения дисциплины «Методы оптимальных

решений» должен обладать способностями:

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

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

5

Содержание

Введение……………………………………………………………………………..3

1.Графический метод решения задач линейного программирования……...5

2.Метод Гаусса-Жордана решения системы линейных уравнений…………6

3.Базисные решения системы линейных уравнений…….……………….…..7

4.Допустимые базисные решения системы линейных уравнений………....9

5.Симплексный метод решения задачи линейного программирования....10

6.Метод потенциалов решения транспортной задачи…………………….....14

7.Метод отсечений решения задач целочисленного программирования…23

Вопросы для контроля.………………………………… ……………………….29

Библиографический список …………………………………………………….29

6

1.ГРАФИЧЕСКИЙ МЕТОД РЕШЕНИЯ ЗАДАЧ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ

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

 

F c1x1 c2 x2

(1.1)

для переменных x1 , x2 удовлетворяющих системе ограничений-неравенств

a11x1 a12 x2 b1,

 

 

 

 

 

 

b2

,

a21x1 a22 x2

 

 

 

 

 

 

 

 

(1.2)

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

 

 

 

 

 

 

 

 

a

m1

x a

m2

x

2

b

 

1

 

 

m

и условиям неотрицательности

 

 

 

 

 

 

x1 0,

x2

0 .

(1.3)

З а м е ч а н и е. Можно рассматривать стандартную задачу линейного

программирования с двумя переменными, в которой требуется найти минимум линейной функции (1.1). Тогда знаки всех неравенств в системе ограниченийнеравенств (1.2) заменяются знаками , а условия (1.3) остаются теми же.

Областью допустимых решений (ОДР) системы неравенств (1.2)

назовем множество всех точек координатной плоскости, координаты каждой из которых являются решением системы неравенств (1.2) и удовлетворяют условиям неотрицательности (1.3). Каждое неравенство ai1x1 ai 2 x2 bi

i 1,2,..., m системы неравенств (1.2) определяет полуплоскость, границей которой является прямая с соответствующим уравнением ai1x1 ai 2 x2 bi , поэтому ОДР системы неравенств (1.2) является пересечением всех этих полуплоскостей и первого координатного угла. Будем считать, что ОДР системы (1.2) не пуста.

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

1) строим область допустимых решений системы неравенств (1.2);

7

2)

строим вектор c c1, c2 , где c1 , c2 – коэффициенты при

соответствующих переменных в функции (1.1);

3)перпендикулярно вектору c строим прямую (линию уровня) с уравнением c1x1 c2 x2 m , где m – какое-либо число (удобно брать m 0 );

4)перемещаем линию уровня параллельно еѐ начальному положению (в направлении вектора c для задачи "на максимум", в противоположном направлении – для задачи "на минимум") до такого положения, при котором она будет иметь общую точку с ОДР и ОДР целиком будет находиться в одной из полуплоскостей, на которые она разбивает всю плоскость; на прямой, имеющей такое положение (опорной прямой), и будет находиться точка (точки), в которой (которых) функция (7.1) достигает максимума или минимума.

2. МЕТОД ГАУССА-ЖОРДАНА РЕШЕНИЯ СИСТЕМ ЛИНЕЙНЫХ УРАВНЕНИЙ

Рассмотрим систему линейных уравнений

a11x1 a12 x2 ...

a1n xn b1,

 

 

 

 

 

 

 

a2n xn

b2 ,

 

a21x1 a22 x2

(2.1)

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

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x a

n 2

x

2

...

a

nn

x

n

b ,

 

 

n1 1

 

 

 

 

n

 

где x1 , x2 , …, xn – неизвестные, aij ,

bi ( i 1, 2, ..., n ;

j 1, 2, ..., n ) – заданные

числа.

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

8

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

Этот метод удобно применять в матричной форме, записав расширенную матрицу

a11

a12

...

a1n

 

b1

 

 

a

 

a

 

...

a

 

 

b

 

 

21

 

22

 

 

2n

 

2

 

... ...

... ...

 

...

 

 

 

 

 

 

 

 

 

 

an1

an 2

...

ann

 

bn

 

 

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

3. БАЗИСНЫЕ РЕШЕНИЯ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Рассмотрим систему линейных уравнений

a11x1 a12 x2

... a1n xn

b1,

 

 

 

 

 

 

... a2n xn

b2

,

a21x1 a22 x2

 

 

 

 

 

 

 

 

 

 

 

(3.1)

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

 

 

 

 

 

 

 

 

 

 

 

 

a

m1

x a

m2

x

2

... a

mn

x

n

b ,

 

1

 

 

 

 

m

в которой число неизвестных больше числа уравнений ( n m). Предположим, что эта система уравнений может быть разрешена относительно каких-то m

неизвестных, например, относительно x1, x2 ,..., xm . Это означает,

что от неѐ

можно перейти (например, применив метод Гаусса-Жордана

и выразив

9

 

неизвестные x1, x2 ,..., xn

через

остальные

неизвестные) к эквивалентной

системе уравнений вида

 

 

 

x1

 

 

 

f1,

a1 m 1 xm 1

... a1n xn

 

 

 

 

 

 

 

 

f2 ,

x2

 

a2 m 1 xm 1

... a2n xn

 

 

 

 

(3.2)

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

 

 

 

 

 

 

 

 

 

 

xm

am m 1 xm 1 ... amn xn fm ,

разрешенной относительно неизвестных x1, x2 ,..., xm .

Неизвестные, относительно которых разрешена система уравнений, называются базисными. Совокупность базисных неизвестных образует

базисную систему неизвестных или, короче, базис. Если выбрана базисная система неизвестных, то остальные неизвестные называются свободными. Свободным неизвестным можно давать произвольные значения и находить соответствующие значения базисных неизвестных. Таким образом, система уравнений (3.1) имеет бесчисленное множество решений. Решение системы (3.1), которое получается при нулевых значениях свободных неизвестных, называется базисным. В нашем случае x1, x2 ,..., xm – базис для системы уравнений (3.1). Придав свободным неизвестным (в системе (9.2)) нулевые значения, получаем базисное решение: x1 f1 , x2 f2 ,…, xm fm ,

xm 1 0 ,…, xn 0 .

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

Система уравнений (3.1) имеет лишь конечное число базисных решений,

не превосходящее Cnm .

Для нахождения всех базисных решений системы (3.1) можно воспользоваться методом последовательного замещения. Опишем кратко этот метод.

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

10

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

4. ДОПУСТИМЫЕ БАЗИСНЫЕ РЕШЕНИЯ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ

Рассмотрим систему линейных уравнений (9.1), для которой выполнены предположения, сделанные в теоретической части предыдущей лабораторной работы.

Решение системы уравнений (3.1) называется допустимым, если значения всех неизвестных в нѐм неотрицательны.

Пусть нам известно невырожденное допустимое базисное решение

системы

уравнений (3.1),

например, решение x0

f 0

,

x0

f

2

0, …,

 

 

 

 

 

 

 

1

1

 

2

 

 

x0

f

m

0,

x0

0 , …,

x0 0 , и система (3.1) приведена методом Гаусса-

m

 

 

m 1

 

n

 

 

 

 

 

 

 

Жордана к соответствующему виду

 

 

 

 

 

 

 

 

 

 

x1

 

h1 m 1 xm 1 ...

h1n xn f1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

h2 m 1 xm 1

h2n xn

f2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.1)

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xm

hm m 1 xm 1

hmn xn fm .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

11

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

отношение fi является наименьшим, и будет номером ключевой строки. Тогда

his

преобразование Гаусса-Жордана приведет к новому допустимому базисному решению.

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

5. СИМПЛЕКСНЫЙ МЕТОД РЕШЕНИЯ ЗАДАЧИ ЛИНЕЙНОГО

ПРОГРАММИРОВАНИЯ

Каноническая задача линейного программирования (КЗЛП) состоит в нахождении максимума (или минимума) линейной функции

F X c1x1

c2 x2 ... cn xn

(5.1)

для переменных x1 , x2 , …, xn ,

удовлетворяющих

системе ограничений-

равенств (9.1) и условиям неотрицательности

 

x j 0,

j 1, 2, ..., n .

(5.2)

Здесь X обозначает набор переменных x1 , x2 , …, xn , функцию (5.1) можно

называть целевой функцией.

Любое решение системы уравнений (3.1), удовлетворяющее условиям неотрицательности (5.2), называется планом канонической задачи линейного программирования. Любое допустимое базисное решение системы уравнений (3.1) называется опорным планом данной задачи. План КЗЛП, при котором линейная функция (5.1) принимает максимальное (или минимальное для задачи на минимум) значение, называется оптимальным.

12

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]