- •Волков Сергей Николаевич землеустройство экономико-математические методы и модели
- •Раздел I
- •Глава 1
- •1.2. Математические методы, применяемые в экономических расчетах
- •Глава 2
- •1. Классификация математических моделей, применяемых в землеустройстве
- •Вычислении площадей треугольников и четырехугольников
- •4.3. Метод решения задачи на условный экстремум лагранжа
- •Глава 5
- •4. Годовые затраты на 1 га пашни в зависимости от размера территории (зерно-свекловичный тип хозяйства с развитым молочно-мясным скотоводством)
- •5.3. Определение оптимальных размеров полей севооборотов
- •5. Расчет оптимального размера поля севооборота
- •Глава 6 итерационные методы
- •6.1. Постановка и математическая формулировка
- •8. Расчет координат оптимального размещения ферм методом последовательных приближений
- •9. Расчет координат животноводческих комплексов
- •10. Расчет оптимальных координат молочного комплекса при селешш Большая Вруда (итерационный метод)
- •11. Анализ эффективности размещения животноводческих комплексов (тыс. Руб., в ценах 1990 г.)
- •12. Расчет значения предельной ошибки целевой функции
- •Раздел III
- •Глава 7
- •Глава 8 расчет параметров производственных функций
- •13. Исходные данные к задаче 8.1
- •8.2. Принцип наименьших квадратов
- •8.3. Системы нормальных уравнений
- •16. Исходные данные к задаче 8.2
- •18. Исходные данные к задаче 8.3
- •8.5. Применение линейных моделей регрессии
- •22. Исходные данные к задаче 8.5
- •Глава 9
- •9.2. Оценка погрешностей определения коэффициентов корреляции
- •25. Корреляционные и дисперсионные характеристики демонстрационных задач
- •26. Формулы для расчета экономических характеристик некоторых однофакторных производственных функций
- •27. Формулы для расчета предельных норм заменяемости для некоторых двухфакторных производственных функций
- •10.2. Примеры расчета экономических характеристик
- •28. Зависимость коэффициента эластичности Ег от стоимости животноводческих построек (х2)
- •Раздел IV
- •Глава 11
- •33. Расчет бета-коэффициентов уравнения регрессии
- •11.3. Обоснование укрупнения (разукрупнения)
- •Контрольные вопросы и задания
- •Глава 12
- •35. Расчет урожайности зерновых культур на землях различных категорий с учетом изменения факторов интенсификации
- •36. Динамика урожайности сельскохозяйственных культур в совхозе «40 лет Октября»
- •38. Планируемая урожайность, рассчитанная с использованием производственных функций, ц с 1 га
- •39. Основные показатели эффективности внесения минеральных удобрений на черноземных почвах
- •40. Расчет планируемой урожайности основных сельскохозяйственных культур по факторам интенсификации расчетпо конструктивным методом
- •41. Возможные урожаи полевых культур при 3 % использования фар
- •42. Сводные данные по планированию урожайности в совхозе «40 лет Октября»
- •Глава 13
- •43. Зависимость между размерами молочной фермы и удельными капиталовложениями
- •44. Расчет параметров уравнения гиперболы
- •46. Влияние механического состава почв на величину удельного сопротивления, кг/см2
- •47. Расчет значений коэффициента а2
- •54. Потери на холостые повороты агрегатов при поперечных работах и от недобора продукции пропашных культур в полосе поворотов при различной густоте лесных полос
- •55. Экономическая эффективность продольных (основных) полезащитных лесных полос при различной густоте посадок
- •57. Чистый доход в расчете на 1 га полевого севооборота
- •58. Влияние густоты сети полевых дорог на общую величину транспортных затрат и потерь от недобора продукции
- •1 Га пашни,
- •1 Га пашни,
- •60. Определение коэффициентов значимости факторов производственной функции
- •Раздел V
- •Глава 14 общая модель линейного программирования
- •63. Исходные данные к задаче 14.2
- •64. Исходные данные к задаче 14.3
- •65. Исходные данные к задаче 14.4
- •67. Характеристика вершин области допустимых значений задачи 14.5
- •68. Первая симплекс-таблица задачи 14.5
- •69. Вторая симплекс-таблица задачи 14.5
- •70. Третья симплекс-таблица задачи 14.5
- •71. Четвертая симплекс-таблица задачи 14.5
- •72. Пятая (последняя) симплекс-таблица задачи 14.5
- •74. Исходные данные к задаче 14.6
- •75. Исходные данные к задаче 14.7
- •76. Оптимальное решение прямой задачи 4.6
- •78. Оптимальное решение прямой задачи 4.7
- •79. Оптимальное решение двойственной задачи
- •Глава 15 распределительная (транспортная) модель
- •80. Исходные данные к задаче 15.1
- •81. Исходные данные к задаче 15.2
- •82. Исходные данные к задаче 15.3
- •83. Исходные данные к задаче 15.4
- •84. Табличная форма представления транспортной модели
- •85. Исходные данные к задаче 15.5
- •86. Нахождение опорного решения задачи 15.5 методом минимального
- •87. Нахождение опорного решения задачи 15.5 методом аппроксимации*
- •15.3. Метод потенциалов
- •88. Цикл испытуемой клетки (3,5)
- •89. Цикл испытуемой клетки (2,5)
- •90. Потенциалы и оценки* для опорного решения задачи 15.5, полученного методом аппроксимации
- •91. Потенциалы и оценки* для опорного решения задачи 15.5, полученного методом минимального элемента
- •92. Потенциалы и оценки на втором шаге решения задачи 15.5
- •93. Потенциалы и оценки на третьем шаге решения задачи 15.5
- •94. Оптимальный план закрепления источников кормов за фермами
- •15.4. Особые случаи постановки и решения распределительных задач
- •95. Исходные данные к задаче 15.6
- •96. Сбалансированная исходная транспортная таблица задачи 15.6
- •97. Исходная транспортная таблица задачи 15.6, в которой учтены требование сбалансированности задачи и первые три дополнительных условия*
- •101. Оптимальный план распределения кормовых культур по участкам
- •102. Исходные данные к задаче 15.7
- •103. Опорный план задачи 15.7
- •104. Второй (оптимальный) план задачи 15.7
- •106. Опорный план задачи 15.3
- •107. Оптимальный план задачи 15.3
- •108. Средние значения коэффициентов урожайности культур в зависимости
- •Глава 16
- •109. Последняя симплекс-таблица для задачи 14.4
- •16.2. Коэффициенты замещения
- •16.3. Использование коэффициентов замещения
- •113. Исходные данные к задаче 16.1
- •114. Оптимальное решение прямой задачи 16.1
- •115. Оптимальное решение двойственной задачи 16.1*
- •16.6. Альтернативные решения распределительных задач
- •117. Исходные данные к задаче 16.2
- •121. Оптимальное решение задачи 16.2 с дополнительным условием (пример 1)*
- •Глава 17
- •122. Исходные данные к задаче 17.1
- •123. Исходная симплекс-таблица задачи 17.1
- •124. Первая расчетная симплекс-таблица 17.1
- •17.3. Роль ограничений в формировании облика
- •Глава 18
- •130. Исходная таблица
- •131. Оптимальный план
- •132. Исходная таблица
- •133. Первый оптимальный план
- •134. Промежуточный опорный план
- •135. Последний оптимальный план
- •136. Вероятностный и детерминированный планы
- •139. Исходные и расчетные данные для вычисления значений ресурсов в ограничениях
- •140. Исходные данные для вычисления коэффициентов целевой функции
- •141. Исходные данные для расчета гц
- •142. Схема двухэтапной стохастической задачи
- •143. Схема числовой стохастической модели оптимизации производственной структуры
- •Контрольные вопросы и задания
- •Раздел VI
- •Глава 19 информационное обеспечение моделирования
- •19.3. Построение матрицы экономико-математической
- •Контрольные вопросы и задания
- •Глава 20
- •146. Баланс гумуса в почве под посевами различных сельскохозяйственных культур
- •Глава 21
- •147. Вычисление значений ак1
- •148. Числовые значения ак1*
- •Раздел VII
- •152. Сведения о максимально возможных объемах и эффективности различных мероприятий по освоению и интенсификации использования земель
- •153. Матрица экономико-математической модели задачи оптимизации мероприятий по освоению и интенсификации использования
- •Контрольные вопросы и задания
- •Глава 23
- •154. Основные переменные
- •155. Исходные данные
- •23.2. Оптимизация трансформации
- •157. Качественная характеристика участков
- •158. Расчет значения Сд для полевого севооборота № 1 по 1-му участку
- •159. Оптимизация трансформации угодий
- •161. Расчет капитальных затрат на трансформацию угодий
- •162. Сводная таблица оценки вариантов, тыс. Руб.
- •Глава 24
- •163. Исходные данные для системы ограничений
- •164. Ориентировочные коэффициенты изменения урожайности культур в зависимости от их предшественников по отношению к средней урожайности хозяйства (зона неустойчивого увлажнения)
- •165. Расчет с,- по полевому севообороту
- •166. Доля сельскохозяйственных культур в рекомендуемых к освоению севооборотах
- •167. Исходные данные для построения экономико-математической модели задачи
- •168. Матрица задачи по проектированию системы севооборотов хозяйства
- •24.2. Размещение севооборотов и сельскохозяйственных
- •170. Фрагмент матрицы оптимального размещения культур (севооборотов) по участкам с различным плодородием
- •171. Фактическое размещение посевов сельскохозяйственных культур
- •172. Оценка предшественников сельскохозяйственных культур
- •173. Матрица задачи по оптимизации плана перехода к запроектированным севооборотам
- •174. Корректировка плана перехода к запроектированным севооборотам
- •175. Структура посевов после корректировки, га
- •176. Окончательный план перехода к запроектированным севооборотам
- •Глава 25
- •25.2. Особенности подготовки исходной информации и пример решения
- •178. Состав и площадь сельскохозяйственных угодий на год землеустройства и по проекту
- •179. Число работников и общий объем трудовых ресурсов
- •Глава 26
- •181. Матрица экономико-математической модели задачи проектирования противоэрозионных мероприятий
- •182. Результаты решения задачи проектирования противоэрозионных мероприятий
- •184. Расчет допустимого слоя стока
- •185. Зависимость площадей линейных элементов организации территории
- •26.3. Оптимизация размещения посевов
- •188. Исходная матрица задачи
- •Глава 27
- •27.2. Особенности подготовки
- •191. Результаты решения задачи организации территории плодовых и ягодных многолетних насаждений
- •Глава 28
- •28.2. Особенности подготовки исходной информации и пример решения
- •192. Допустимые площади кормовых культур и пастбищ*
- •193. Расчет потребности в зеленом корме
- •194. Расчет потребности в кормах с пашни
- •Контрольные вопросы и задания
- •Глава 29
- •199. Результаты решения экономико-математической задачи
- •Контрольные вопросы и задания
- •Раздел VIII
- •Глава 30
- •30.2. Особенности подготовки исходной информации и пример решения
- •200. Продолжительность рабочего периода в крестьянском хозяйстве
- •201. Нормы внесения минеральных удобрений под сельскохозяйственные культуры, кг д. В. На 1 т продукции*
- •202. Технико-экономические характеристики животноводческих хозяйств*
- •204. Математическая модель молочно-картофелеводческого крестьянского хозяйства
- •205. Выход кормов с 1 га культурных пастбищ
- •206. Оптимальные размеры землевладений и структура производства в крестьянском хозяйстве молочно-картофелеводческого направления
- •30.3. Автоматизация расчетов модели на эвм
- •207. Значения переменных задачи
- •31.2. Особенности подготовки
- •209. Состав земельных угодий до и после перераспределения земель
- •Глава 32
- •32.1. Экономико-математическая модель
- •32.2. Экономико-математическая модель
- •210. Исходная матрица задачи
- •211. Оптимальный план формирования сырьевых зон перерабатывающих предприятий
- •Контрольные вопросы и задания
- •Литература
- •Раздел I. Общие сведения об экономико-математических методах и моделировании в землеустройстве 9
- •Глава 1. Моделирование и современные методы вычислений 9
- •Глава 2. Основные этапы развития математического моделирования в аграрно- экономической и землеустроительной науке 32
- •Глава 3. Классификация математических моделей, применяемых в земле устройстве 57
- •Раздел II. Аналитическое моделирование в земле устройстве 72
- •Глава 4. Построение и исследование аналитических моделей 72
- •Глава 5. Применение дифференциального и интегрального исчисления при построении оптимизационных аналитических моделей 92
- •Глава 9. Оценка производственных функций с использованием методов корре ляционно-регрессионного анализа 161
- •Глава 10. Экономические характеристики производственных функций и
- •Раздел IV. Применение производственных функций
- •Глава 11. Оптимизация интенсивности использования земли при землеуст ройстве 197
- •Глава 12. Планирование урожайности сельскохозяйственных культур 209
- •Глава 13. Разработка землеустроительных нормативов и решение нестан дартных задач 234
- •Раздел V. Методы математического программирования
- •Глава 14. Общая модель линейного программирования 261
- •Глава 15. Распределительная (транспортная) модель 303
- •Глава 16. Анализ и корректировка оптимальных решений 344
- •Глава 17. Дополнительные аспекты решения задач линейного программиро вания 383
- •Глава 18. Некоторые виды задач математического программирования 398
- •Раздел VI. Основы экономико-математического моделирования 436
- •Глава 19. Информационное обеспечение моделирования 436
- •Глава 20. Выбор переменных и построение ограничений задачи 451
- •Глава 21. Критерии оптимальности при решении землеустроительных задач 474
- •Раздел VII. Экономико-математические модели
- •Глава 22. Экономико-математическая модель оптимизации мероприятий
- •Глава 23. Экономико-математическая модель трансформации угодий 506
- •Глава 24. Экономико-математическая модель организации системы сево оборотов хозяйства 519
- •Глава 25. Экономико-математическая модель оптимизации структуры посевных площадей при агроэкономическом обосновании проектов внутрихо зяйственного землеустройства 553
- •Глава 26. Экономико-математическая модель проектирования комплекса противоэрозионных мероприятий в условиях развитой водной эрозии почв 566
- •Глава 31. Экономико-математическая модель оптимизации перераспреде ления земель сельскохозяйственных предприятий 614
- •Глава 32. Экономико-математические модели в схемах землеустройства 659
115. Оптимальное решение двойственной задачи 16.1*
Лз/л)
(осн.)
Базисные переменные
Номера ограничений
(для дополнительных
переменных)
(значение
базисных
переменных)
(осп.)
Коэффициенты замещения
А\}(у'4) (осн.)
(изб. в огр. 1)
(изб. в огр. 2)
у2(осн.) 5000 -
у7 (изб.) — 3
уц (изб.) - 4
у; (осн.) -24000 -Индексная строка (Щ — В^
80 150 1,25
1,25 770000
0,04
-0,8
-0,0033
-0,0033 -520
0,4
-2
-0,0021
-0,0021 -920
0
0
0,0042
0,0042 -100
-0,04
-0,2
-0,0008
-0,0008 -380
*В соответствии с (16.15) и значениями у'4,у% из таблицы получаем у4 — -1,25.
В п. 14.8 мы уже указывали на соответствие между элементами индексной строки прямой задачи и значениями основных переменных, попавших в базис оптимального решения двойственной задачи. Проведем теперь более полный анализ именно с этой точки зрения.
Примем во внимание следующие факты:
элементы индексной строки последней симплекс-таблицы прямой задачи, соответствующие базисным переменным (в табл. 114 не показаны), равны нулю;
значения небазисных переменных двойственной задачи также равны нулю;
в каждое ограничение прямой задачи в стандартной форме входит точно одна дополнительная переменная (остаточная или избыточная) и с этим же ограничением ассоциирована одна основная переменная двойственной задачи (см. систему 16.12). Следовательно, существует попарное соответствие между дополнительными переменными прямой задачи и основными переменными двойственной задачи: х4—>у4 (4-е ограничение прямой задачи); х5->У] (1-е ограничение); х6 ->у2 (2-е ограничение); х7 -»у3 (3-е ограничение);
значения целевых функций 7, (прямой) и Ж (двойственной) задач в оптимуме совпадают.
С учетом сказанного сопоставление таблиц 114 и 115 показывает, что все элементы индексной строки прямой задачи, соответствующие дополнительным переменным, точно равны ассоциированным с ними двойственным переменным с учетом знака, с которым дополнительная переменная входит в ограничение прямой задачи:
(24 - С4) = -у4 = 1,25; (2Г5 - С5) = у, = 0; (26-С6)=у2 = №, (2"7-С7)=>;з = 0. (16.16)
Дальнейшие выводы строятся на анализе выражения для целевой функции двойственной задачи. Запишем его для задачи 16.1 в общем виде:
У/= Ьху{ + Ь2у2 + Ьъуъ + Ь4у4, (16.17)
где в качестве коэффициентов при основных переменных двойственной задачи используются правые части ограничений прямой задачи — Ь„ /'= 1,...,4, то есть уровни ресурсов, которыми обладает хозяйство, и уровни плановых заданий по производству отдельных видов продукции.
Для того чтобы прямая и двойственная целевые функции совпадали в оптимуме не только по численному значению, но и по единицам измерения, необходимо, чтобы в выражении (16.17) каждая двойственная переменная имела размерность [у,] =
361
= [руб.]/[единица измерения соответствующего ресурса или планового задания], то есть фактически являлась «ценой» данного ресурса или планового задания. Например, у2 является «ценой» ресурса Ь2, то есть органических удобрений, у4 является «ценой» планового задания по производству зеленой массы многолетних трав (Ь4).
В литературе такие «цены» получили название скрытых или теневых. Впервые они были исследованы Л. В. Канторовичем, который называл их объективно обусловленными оценками.
В целом с учетом (16.16) можно констатировать, что:
скрытые цены ресурсов равны элементам индексной строки последней симплекс-таблицы прямой задачи, соответствующим остаточным переменным;
скрытые цены плановых заданий равны взятым с обратным знаком элементам индексной строки последней симплекс-таблицы прямой задачи, соответствующим избыточным переменным.
В соответствии с (16.16) элементы индексной строки прямой задачи часто называют двойственными оценками. Рассмотрим содержание этих оценок более детально.
В соответствии с соотношением (16.2) и рассмотренной выше экономической интерпретацией знака дополнительных переменных различных типов, вводимых в базис, двойственная оценка определяет изменение целевой функции при изменении соответствующего ресурса на единицу, то есть характеризует эффективность использования (ценность) данного ресурса и смысле достижения требуемого эффекта (увеличения целевой функции). Аналогично, если анализируется роль планового задания, то двойственная оценка может рассматриваться как цена увеличения задания в плане потерь в достижении требуемого эффекта.
В рассматриваемой задаче 16.1 двойственная оценка ресурса органических удобрений составляет 80 руб. за 1 т, то есть при увеличении ресурса органических удобрений на 1 т прибылI, увеличивается на 80 руб. Это полностью согласуется с тем фактом, что органические удобрения являются дефицитным ресурсом, так как соответствующая остаточная переменная является небазисной. В то же время недефицитные ресурсы (пашня и трудовые ресурсы) имеют нулевые двойственные оценки, то л есть изменения в наличии этих ресурсов (в определенных пределах) не сказываются на значениях целевой функции. Двойственная оценка планового задания по производству зеленой массы многолетних трав составляет 1,25 руб. за 1 ц, то есть уие личение их производства на 1 ц приведет к снижению прибыли на 1,25 руб.
Двойственные оценки основных небазисных переменных прямой задачи (то есть соответствующие этим переменным элементы индексной строки) можно также назвать скрытыми ценами не
362
эффективных отраслей производства. Такая цена показывает, чего будет стоить развитие неэффективной отрасли с точки зрения достижения поставленной цели. В рассматриваемой задаче 16.1 двойственная оценка производства зеленого горошка составляет 150 руб. на 1га, то есть возделывание каждого гектара пашни под зеленый горошек будет снижать общую прибыль на 150 руб.
В практическом отношении важно, что двойственные оценки прямой задачи позволяют выявить наличие альтернативных оптимальных решений. Действительно, из соотношения (16.2) прямо следует, что если двойственная оценка какой-то небазисной переменной равна нулю, то введение ее в оптимальный план не изменит значения целевой функции, а новое решение, соответствующее любому допустимому значению вводимой переменной, также будет оптимальным. Здесь можно провести аналогию с транспортными задачами, в которых признаком наличия альтернативных оптимальных решений является равенство нулю оценки какой-либо свободной клетки (см. анализ задачи 15.7).
Наличие альтернативных оптимальных решений расширяет возможности землеустроителя при выборе рационального проекта с учетом дополнительных факторов, прямо не учтенных в постановке задачи линейного программирования.
Принципиальная особенность любой двойственной оценки заключается в том, что она отражает влияние на процесс производства не только данного ресурса, но и представляет собой интегрированный показатель его влияния в комплексе со всеми остальными. Фактически двойственные оценки формируются под влиянием всех факторов, определяющих постановку задачи, — типов ограничений, коэффициентов при неизвестных в целевой функции и в ограничениях, правых частей ограничений. Можно сказать, что реальные рыночные цены на данные виды ресурсов или производимую продукцию являются только одной — «видимой» составной частью факторов формирования двойственных оценок. Это объясняет другое наименование двойственных оценок—«скрытые цены».
Из определения двойственных оценок следует, что их экономическое содержание распространяется только на данную задачу. При изменении исходных данных (например, коэффициентов целевой функции) двойственные оценки могут существенно измениться.
Абсолютная величина двойственных оценок, соответствующих остаточным переменным, характеризует степень эффективности ресурсов, находящихся в дефиците. Так, в задачах на максимум чем выше значение двойственной оценки, тем данный ресурс дефицитнее.
363
Как известно, в сельском хозяйстве используются самые разнообразные ресурсы — земля, труд, машины, топливо, семена, удобрения и т. д. При существенном ограничении объемов данных ресурсов они получают в оптимальном плане ненулевые двойственные оценки. При этом экономическое содержание двойственных оценок зависит от степени детальности выбора переменных, а также от структуры отдельных ограничений и их системы в целом. Например, если дифференциация переменных в процессе постановки задачи ведется по ряду земельных участков различного плодородия, то двойственные оценки показывают эффективность применения различных средств производства на этих землях. Они могут отражать, в частности, размер дифференциальной земельной ренты по плодородию, местоположению и по фактору повышения интенсивности производства.
При построении ограничений по затратам ручного и механизированного труда двойственные оценки будут указывать на резервы повышения производительности труда в конкретном хозяйстве. Например, более высокая двойственная оценка по затратам механизированного труда укажет на необходимость более широкого использования, улучшения подготовки механизаторских кадров, повышения их квалификации.
Двойственные оценки дефицитных ресурсов позволяют выявить различные направления улучшения хозяйства. Так, высокая двойственная оценка по трудовым ресурсам свидетельствует о необходимости первоочередного привлечения рабочей силы со стороны.
Двойственные оценки могут применяться также при планировании объема производственных ресурсов и денежных средств на перспективу, а также указывать на наиболее целесообразное распределение ресурсов между производственными подразделениями хозяйства.
В заключение данного раздела дополнительно проанализируем экономическое содержание двойственных оценок, основываясь на теореме равновесия (Ларченко Е. Г. Вычислительная техника и экономико-математические методы в землеустройстве. — М.: Недра, 1978.— С. 295). Рассмотрим случай, когда в задаче имеются только ресурсные ограничения (в соответствии с принятыми обозначениями заданные объемы ресурсов — это А,, /= \,...,т). Согласно теореме равновесия в оптимальных решениях прямой и двойственной задач соответствующие двойствен и 1.1 с оценки, то есть значения двойственных переменных у, по отношению к прямой задаче и значения переменных ху по отноше нию к двойственной задаче, равны нулю при выполнении следу ющих условий:
364
у(°)=0, как только X а»1 ° <А/, и наоборот;
х^=0, как только Ё %у^<с,-, и наоборот.
Символ (о) при переменных указывает на оптимальность решений. Обратите внимание на то, что в неравенствах использованы знаки строгого неравенства.
В соответствии со смыслом переменных прямой и двойственной задач, а также коэффициентов ау, Ьь с,-, формирующих прямую задачу, указанные выше суммы следует интерпретировать так:
т . .
2., ОуХу _ общий объем использования /-го ресурса во всех отраслях хозяйства;
ЬиуУ{ _ затраты на единицу продукции ву-и отрасли хозяи-
1=1
ства.
Экономически интерпретируя данную теорему, можно сделать следующие выводы.
Во-первых, если по оптимальному плану имеется недоисполь-
т , .
чованный ресурс (X аух) <^) и, значит, соответствующая остаточная переменная находится в базисе, то цена на этот ресурс в <кшном хозяйстве падает до нуля. Расходование остатка этого ресурса может иметь смысл только при увеличении ресурсов, находящихся в дефиците. Положительные цены имеют лишь те ресурсы, которые лимитируют дальнейшее развитие производства. Например, если недоиспользуются пашня и трудовые ресурсы, а денежно-материальные средства израсходованы полностью, то положительную цену будет иметь только ресурс, указанный пос-иедним. В рассмотренном выше примере (см. задачу 16.1, табл. I 14) положительную цену имел только один ресурс — органические удобрения.
Во-вторых, если затраты на производство у'-го продукта пре-
пышают получаемый от него доход (Х0//), >с/)> то соответству-
/=1 „
ющаяу-я отрасль исключается из оптимального плана. В указанном примере такой отраслью является производство зеленого горошка (/'= 3, см. табл. 114).
365
Таким образом, использование двойственной задачи позволяет:
осуществить анализ и контроль решения прямой задачи;
оценить оптимальность любого базисного решения (по данным индексной строки);
облегчить решение задачи при числе ограничений, значительно превышающем число переменных (заменив прямую задачу двойственной);
оценить дефицитность и эффективность использования различных ресурсов путем сравнения двойственных оценок;
находить альтернативные оптимальные решения по признаку равенства нулю двойственных оценок при переменных, не вошедших в оптимальный план в последней симплекс-таблице.
Тем самым двойственные оценки становятся важным инструментом анализа оптимальных решений и в этом качестве могут служить для обоснования действий, направленных на повышение эффективности производства.
16.5. ПРЕДЕЛЫ УСТОЙЧИВОСТИ ОПТИМАЛЬНОГО РЕШЕНИЯ ПРИ ИЗМЕНЕНИИ КОЭФФИЦИЕНТОВ ЦЕЛЕВОЙ ФУНКЦИИ
Фактически мы уже рассмотрели вопрос об устойчивости структуры оптимального решения (то есть списка базисных переменных) и элементов индексной строки (двойственных оценок) при изменении дефицитных ресурсов и критических плановых заданий, а также при развитии неэффективных отраслей хозяйства. Было показано, что при изменении этих факторов в некоторых пределах структура решения и двойственные оценки сохраняются. Конкретно указанные пределы определяются в рамках алгоритмов введения небазисных переменных в оптимальный план. При этом, однако, значения базисных переменных и целевой функции меняются.
Практически полезным свойством оптимальных решений задач линейного программирования является определенная устойчивость этих решений к изменению коэффициентов целевой функции (их иногда называют коэффициентами удельной прибыли). А именно при изменении этих коэффициентов в определенных пределах решение сохраняется как по составу базисных переменных, так и по их значениям, однако значение целевой функции и двойственные оценки при этом меняются. Таким образом, использование указанного свойства расширяет возможности по корректировке оптимального решения задачи в случае относительно малых изменений коэффициентов целевой функции При этом повторного применения симплекс-метода, так же кпк и при корректировке решения с помощью коэффициентов замещения, не требуется.
366
Алгоритмы корректировки значения целевой функции и двойственных оценок небазисных переменных прямой задачи при изменении коэффициентов целевой функции 2Г основываются на детальном сопоставлении решений прямой и двойственной задач; для иллюстрации мы воспользуемся уже имеющимися решениями задачи 16.1 в прямой и двойственной постановках (см. табл. 114 и 115). Все расчеты ведутся только для случая максимизации целевой функции; формулы для случая минимизации легко получить по аналогии.
Сначала рассмотрим изменение коэффициентов целевой функции при переменных, входящих в базис оптимального решения прямой задачи. Допустим, например, что коэффициент при переменной х{ в целевой функции I (16.11) изменяется следующим образом:
С[=С,+ДС, (16.18)
где Ас — изменение коэффициента, которое может быть как положительным, так и отрицательным.
Если решить задачу с новым значением коэффициента, то изменение претерпит только индексная строка, причем новое значение у'-го элемента индексной строки можно рассчитать, основываясь на уже имеющемся решении прямой задачи, по формуле
(4- С}у = (3- С,-) + АуДс, у > 0, (16.19)
где Л4у—коэффициенты замещения, расположенные в симплекс-таблице (см. табл. 114) в строке, соответствующей базисной переменной хь то есть в 4-й строке.
Таким образом, формула (16.19) позволяет определить новые значения как целевой функции 2'0^={2^-С^)', так и всех двойственных оценок небазисных переменных: {2} — С])' при у > 1. При этом так как рассматриваемое изменение задачи не затрагивает структуры последней симплекс-таблицы, элементы индексной строки, соответствующие базисным переменным (в последней симплекс-таблице они не показаны), сохраняют нулевые тачения.
Схему рассуждений, приводящих к формуле (16.19), кратко поясним следующим образом.
Изменение коэффициента Сх в целевой функции 2 одновременно означает идентичное изменение правой части первого ограничения системы ограничений (16.14) двойственной задачи. Следовательно, вариацию Асиз формулы (16.18) можно рассматривать как изменение «планового задания» в двойственной зада-
367
че, соответствующего этому ограничению (тип ограничения — «>»). В то же время если решение двойственной задачи рассматривать как обычное решение симплексной задачи (16.13)— (16.16), то результатом указанного изменения правой части первого ограничения будет обычная корректировка оптимального решения двойственной задачи за счет введения в оптимальный план избыточной переменной у5 = Ас (см. табл. 115). Как известно, в процессе такой корректировки будут меняться базисные переменные двойственной задачи и значение целевой функции Ж При этом новые значения базисных переменных и целевой функции будут определяться с использованием коэффициентов замещения и элемента индексной строки таблицы 115, соответствующих избыточной переменной у5. Теперь остается вспомнить, что:
каждая базисная переменная двойственной задачи соответствует определенной двойственной оценке прямой задачи;
в оптимуме значения целевых функций прямой и двойственной задач совпадают.
Если учесть, что между множеством чисел столбца у5 последней симплекс-таблицы двойственной задачи, включая элемент индексной строки, и числами 4-й строки последней симплекс-таблицы прямой задачи, включая элемент из столбца А®, существует взаимно однозначное соответствие с точностью до изменения знака на обратный (ср. табл. 114 и 115), то из формул (16.1) и (16.2) непосредственно вытекает формула (16.19).
Допустимые пределы изменений величины Ас (в рамках которых структура симплекс-таблицы не меняется) можно определить из тех же соображений, воспользовавшись уже рассмотренным выше алгоритмом ввода в оптимальный план избыточной переменной. Однако возможен и более прямой путь, которым мы и воспользуемся.
Из правила определения оптимальности симплекс-таблицы (и задачах на максимум — неотрицательность элементов индексной строки) следует, что при применении соотношения (16.19) новые значения (2} — С,)' приу> 1 не должны быть отрицательными. Таким образом, в случае максимизации целевой функции для любого/^ 1 должно выполняться условие
(2}- С]) + М} Ас^ 0. (16.20)
Следовательно, для определения пределов допустимых изменений Дс необходимо разделить (2} — С-) на соответствующие ненулевые значения Л4у и частные взять с обратным знаком. Среди положительных частных от такого деления нужно найти наименьшее (обозначим его 0^т), а среди отрицательных — наи меньшее по модулю (обозначим его Цп\п)- Затем следует прове-
368
рить, не превышает ли значение Втт по модулю коэффициент
Сь и если так, положить Д~;п=-С1. Такое действие необходимо, чтобы предотвратить появление отрицательного коэффициента удельной прибыли по отрасли хозяйства, соответствующей переменной х{ (см. 16.18). В результате получим следующее ограничение:
Очевидно, что если среди коэффициентов Ац нет отрицательных, то вариация Дс может быть сколь угодно большой положительной, а если нет положительных — должно выполняться условие Ят1П=-С,.
Пусть в задаче 16.1 чистый доход от реализации томатов уменьшается с 2000 до 1500 руб. с 1 га, то есть Ас= —500. Определить, как при этом изменяются оптимальное значение чистого дохода, двойственные оценки неэффективной отрасли (производство зеленого горошка), дефицитного ресурса (органические удобрения) и критического планового задания (по выращиванию многолетних трав).
Учитывая, что меняется коэффициент С2 при переменной х2 (см. целевую функцию 16.11), для определения допустимых значений Ас разделим элементы индексной строки на ненулевые коэффициенты замещения, стоящие в строке, соответствующей переменной х2 (см. табл. 114). Используя полученные результаты, на основании рассмотренного выше алгоритма оценим допустимые пределы изменений С2:
-2000 <ДС< 750 (руб. с 1 га).
Таким образом, предполагаемое изменение чистого дохода от реализации томатов находится в допустимых пределах. Используя формулу (16.19), получим:
чистый доход по хозяйству в целом уменьшится на 190 000 руб.;
скрытая цена (двойственная оценка) производства зеленого горошка уменьшится на 100 руб. с 1 га;
скрытая цена органических удобрений уменьшится на 20 руб. за 1 т;
скрытая цена критического планового задания по выращиванию многолетних трав уменьшится на 0,4 руб. за 1 ц.
Теперь следует рассмотреть случай, когда изменяется коэффициент в целевой функции при небазисной переменной — например, переменной хз в задаче 16.1 (см. табл. 114). Пусть
369
Сз=с3+лс,
(16.21)
где Дс— вариация коэффициента, которая может принимать как положительные, так и отрицательные значения.
В данном случае это приведет к изменению только элемента индексной строки, соответствующего рассматриваемой переменной:
(23-С3)' = (23-С3)-ДС. (16.22)
Потребовав, чтобы новое значение этого элемента не было отрицательным, получим следующее ограничение на допустимые изменения С3 (для случая максимизации целевой функции):
-С3<Дс<(2з-С3). (16.23)
Пусть в задаче 16.1 чистый доход от реализации зеленого горошка увеличивается с 250 до 350 руб. с 1 га, то есть Дс= 100. Определить, как при этом изменится скрытая цена этой неэффективной отрасли.
Учитывая, что меняется коэффициент при небазисной переменной х3 (см. целевую функцию 16.11), получим следующие допустимые пределы его вариации:
-250<Лс<150(руб. с 1га).
Таким образом, предполагаемое изменение чистого дохода от реализации зеленого горошка находится в допустимых пределах. Используя формулу (16.22), получим, что скрытая цена производства зеленого горошка уменьшится на 100 руб. с 1 га. При этом другие двойственные оценки и чистый доход по хозяйству в целом не изменятся.
Заметим, что изменение коэффициента при небазисной переменной в допустимых пределах в принципе не может перевести соответствующую неэффективную отрасль в число эффективных. В то же время в соотношениях вида (16.23) содержится важная информация экономического характера: например, из них следует, что при прочих неизменных исходных данных производство зеленого горошка станет эффективным, если будет выполнено условие
АС>(73-С3),
то есть если в данном конкретном примере коэффициент удельной прибыли по зеленому горошку (С3) будет больше 400 руб. с 1 га. Для иллюстрации этого утверждения в таблице 116 приведс-
370
но оптимальное решение задачи 16.1 при значении коэффициента С3 = 410.
116. Оптимальное решение задачи 16.1 при увеличении удельной прибыли по зеленому
горошку до 410 руб. с 1 га
Л,7
(*т)
(ост.
в огр.
3)
(изб. в огр. 4)
Базисные переменные
С,-
Номера ограничений (для дополнительных переменных)
А® (значения базисных переменных)
Коэффициенты замещения
Л (*й) (ост. в огр. 2)
х5 (ост.) — 1
л:, (осн.) 2000 -
х3(осн.) 410 -
х, (осн.) 100 -Индексная строка (2^— С,)
140 |
0,0025 |
0,12 |
-0,4 |
285 |
0,0006 |
0,08 |
-0,1 |
475 |
0,001 |
-0,2 |
0,5 |
100 |
-0,004 |
0 |
0 |
774750 |
1,26 |
78 |
5 |
Возрастание коэффициента удельной прибыли по зеленому горошку привело к переходу данной отрасли в число эффективных. Одновременно произошли следующие изменения (ср. табл. 116 и 114):
увеличился чистый доход хозяйства в целом на 4750 руб.;
более полно стала использоваться пашня (уменьшилась остаточная переменная х$, соответствующая 1-му ограничению);
теперь полностью исчерпаны трудовые ресурсы (остаточная переменная х-], соответствующая 3-му ограничению, перешла в число небазисных);
уменьшилась площадь пашни под томатами (переменная х2).
