- •Л.В. Водолазская, в.С. Пецевич математическое моделирование социально-экономических процессов
- •Оглавление
- •Предисловие
- •Методические рекомендации
- •Введение
- •Лекция 1 Основные понятия и определения
- •1.1 Основные понятия и определения математического программирования
- •1.2. Основные понятия и определения математического моделирования
- •Контрольные вопросы
- •Лекция 2 Симплексный метод линейного программирования
- •2.1. Общая характеристика симплексного метода
- •2.2. Решение задачи линейного программирования в симплексных таблицах. Правила построения симплексных таблиц
- •Определение оптимальности плана. Построение новой симплексной таблицы
- •2.3. Альтернативный оптимум
- •2.4. Вырождение основной задачи линейного программирования
- •Контрольные вопросы
- •Лекция 3 Метод искусственного базиса или м - метод
- •Контрольные вопросы
- •Лекция 4 Транспортная задача
- •4.1. Нахождение опорного плана транспортной задачи
- •4.2. Нахождение оптимального плана методом потенциалов
- •Контрольные вопросы
- •Лекция 5 Оптимизация структуры посевных площадей овощных культур.
- •5.1. Постановка задачи.
- •5.2. Состав переменных и ограничений
- •5.3. Структурная экономико-математическая модель
- •5.4. Исходная информация
- •5.5. Разработка числовой экономико-математической задачи
- •5.6. Анализ оптимального решения.
- •Контрольные вопросы
- •Лекция 6 Оптимизация структуры посевных площадей зерновых культур с учетом предшественников
- •6.1 Постановка задачи
- •6.2. Состав переменных и ограничений
- •6.3. Исходная информация
- •6.4. Разработка числовой экономико-математической модели
- •6.4. Анализ оптимального решения
- •Контрольные вопросы
- •7.Оптимизация рационов кормления животных
- •7.2. Состав переменных и ограничений задачи.
- •7.3. Исходная информация Для составления экономико-математической модели оптимального рациона кормления скота необходимы следующие данные:
- •7.4. Разработка числовой экономико-математической модели
- •Питательная ценность и стоимость кормов (в расчете на 1 кг корма)
- •7.5. Анализ оптимального решения
- •Контрольные вопросы
- •8. Оптимизация производственной структуры сельскохозяйственного предприятия
- •8.1.Постановка задачи
- •8.2. Система переменных и ограничений
- •8.3. Подготовка исходной информации и составление числовой экономико-математической модели
- •8.3.Анализ оптимального решения
- •Контрольные вопросы
- •9. Оптимизация плана производства кормов
- •9.1. Постановка задачи
- •9.2. Состав переменных и ограничений
- •Подготовка исходной информации и составление числовой экономико-математической модели
- •Анализ оптимального решения
- •Контрольные вопросы
- •Заключение
- •Библиографический список литературы
- •Типография издательства ОмГау, Омск-8, Сибаковская, 4
Определение оптимальности плана. Построение новой симплексной таблицы
В симплексных таблицах формальным признаком оптимальности является содержание оценочной строки.
Ключевая теорема симплексного метода (Z на максимум): Если после выполнения очередной итерации:
|
Согласно данной теоремы:
1. Просматриваются оценки Z - строки. Если они все неотрицательны, то получено оптимальное решение при решении на максимум целевой функции, если все оценки 0, то получено оптимальное решении при решении на минимум целевой функции.
2. Если оптимальное решение не получено, то выбирается разрешающий столбец по наибольшей по абсолютной величине отрицательной оценке Z - строки при решении на максимум и по наибольшей положительной оценке при решении на минимум целевой функции.
3. Составляются симплексные отношения - отношения свободных членов к положительным коэффициентам разрешающего столбца. Минимальное из этих отношений ( min ai 0/ aip) определит разрешающую строку. Коэффициент, находящийся на пересечении разрешающей строки и разрешающего столбца, называется разрешающим элементом.
Если в разрешающем столбце нет ни одного положительного коэффициента , то задача не имеет решения по причине неограниченности целевой функции (Zmax+ ,Zmin- ).
4. Осуществляется переход к новой таблице, где базисная переменная заменяется на свободную, соответствующую разрешающему столбцу.
5. Бывший разрешающий элемент заменяется обратным по величине (1/ аqp).
6. Все остальные элементы бывшего разрешающего столбца делятся на разрешающий элемент и меняют знаки на противоположный.
Если в бывшем разрешающем столбце в какой-то строке стоял “0”, то эта строка переносится в следующую таблицу без изменений.
7. Все остальные элементы бывшей разрешающей строки делятся на разрешающий элемент .
Если в бывшей разрешающей строке в каком-то столбце стоял “0”, то этот столбец переносится в следующую таблицу без изменений.
8. Остальные элементы таблицы пересчитываются по правилу “прямоугольника”:
aij= |
Исходный коэффициент |
- |
Коэффициент разрешающей строки |
х |
Коэффициент разрешающего столбца |
Разрешающий элемент
|
aij . . . . aip
. . . . . . . . . aij = aij - aqj *aip/ aqp , где
aqj. . . (aqp) i =1,2,..... j=0,1,..... , т.е. и столбцов свободных
членов.
9. Осуществляется контроль за правильностью расчетов для элементов Z-строки по формуле (2).
10. Если ошибок нет, то алгоритм повторяется с пункта 1.
В первой симплексной таблице нашего примера все коэффициенты оценочной строки отрицательны. Следовательно, согласно теореме, план (на максимум целевой функции) не является оптимальным и требует улучшения.
В последнем столбце рассчитывают симплексные отношения.
Таблица 2.1
Симплексная таблица (исходное опорное решение)
-
Св.П
Cj
-2
1
2
Б.П.
Ci
ai0
X1
X2
ai0/aip
X3
0
4
3
1
4
X4
0
4
1
(3)
4/3
Z
2
-1
-2
Исходное опорное решение: Х1(0,0,4,4); Zmax=2
Во второй симплексной таблице меняем местами базисную переменную X4 и свободную X2.
На месте бывшего разрешающего элемента запишем обратную величину, т.е. 1/3.
Все остальные элементы разрешающей строки в новой таблице вычислим путем деления на разрешающий элемент: 4/3; 1/3.
Все остальные элементы бывшего разрешающего столбца определим путем деления на разрешающий элемент. Результат запишем в новую таблицу с противоположным знаком: -1/3; 2/3.
Все остальные элементы таблицы вычислим по правилу «прямоугольника».
-
Первая строка:
4 – 4 1/3 =8/3
3 – 1 1/3 =8/3
Оценочная строка:
2 – 4 (-2)/3 =14/3
-1 – 1 (-2)/3 = -1/3
Таблица 2.2.
Вторая симплексная таблица
-
Св.П
Cj
-2
1
0
Б.П.
Ci
ai0
X1
X4
ai0/aip
X3
0
8/3
(8/3)
-1/3
1
X2
2
4/3
1/3
1/3
4
Z
14/3
-1/3
2/3
Проверим правильность заполнения Z – строки по формуле (2).
а00 = 0 8/3 + 2 4/3 – (-2) = 14/3
а01 =0 8/3 + 2 1/3 – 1 = -1/3
а02 =0 (-1/3) + 2 1/3 – 0 = 2/3
Так как в оценочной строке есть еще отрицательные оценки, оптимальное решение не получено, можем записать только опорное решение:
X2( 0,4/3,8/3,0) Zmax=14/3
По тому же алгоритму пересчитываем коэффициенты четвертой таблицы.
Таблица 2.3.
Третья симплексная таблица
-
Св.П
Cj
-2
0
0
Б.П.
Ci
ai0
X3
X4
X1
1
1
3/8
-1/8
X2
2
1
-1/8
3/8
Z
5
1/8
5/8
Все оценки неотрицательны, следовательно, получено оптимальное решение Хопт(1,1,0,0) Zmax=5
Пример2
Задача решалась на максимум функции (Z max). На последнем шаге была получена следующая таблица.
Таблица 2.4.
Последняя симплексная таблица
-
Св.П
Cj
0
1
0
Б.П.
Ci
ai0
X3
X2
X1
2
2
1
-1
X4
0
1
1
-1/2
Z
4
1
-2
В разрешающем столбце нет ни одного положительного элемента, следовательно, Zmax ,т.е. задача линейного программирования не имеет решения по причине неограниченности целевой функции Z сверху.