- •Волков Сергей Николаевич землеустройство экономико-математические методы и модели
- •Раздел 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
Глава 21
КРИТЕРИИ ОПТИМАЛЬНОСТИ ПРИ РЕШЕНИИ ЗЕМЛЕУСТРОИТЕЛЬНЫХ ЗАДАЧ
21.1. МОДЕЛИРОВАНИЕ ЦЕЛЕВОЙ ФУНКЦИИ
Выбор критерия оптимальности — один из наиболее важных и ответственных этапов моделирования. Даже при самой тщательной постановке и математической формулировке землеустроительной проектной задачи, обосновании системы переменных и условий, адекватно отражающих действительность, неудачно выбранный критерий оптимальности может привести к неудов-
474
летворительным решениям и исказить целевую установку проекта землеустройства.
Возникновение понятия «критерий оптимальности» было обусловлено разработкой оптимизационных моделей, в которых задавалось достижение экстремального (максимального или минимального) значения какого-то экономического результата. В этой связи требовалось не только качественно определить показатели экономической функции, но и аналитически выразить их в виде конкретной математической функции. Поэтому в теории моделирования возникло такое понятие, как целевая функция.
Целевая функция (функционал, целевая установка, функция цели) — это аналитическая форма выражения критерия оптимальности задачи.
В связи с тем что развитие экономико-математического моделирования в землеустройстве началось с постановки и решения частных задач по оптимизации трансформации угодий, проектированию севооборотов, планированию структуры посевов и т.д., в качестве критериев оптимальности использовались простейшие показатели оценки эффективности (эффекта) землеустройства. При этом эффективность землеустройства сводилась к экономической эффективности развития сельскохозяйственного производства на конкретных сельскохозяйственных предприятиях, а в качестве целевых установок задач применялись критерии:
максимизирующие валовую и товарную продукцию в стоимостном и натуральном выражении, валовой, чистый доход, прибыль, рентабельность производства, производительность труда и др.;
минимизирующие приведенные затраты, затраты труда, материально-денежных средств, некоторые виды ресурсов (пашни, кормов и т. д.), себестоимость продукции;
обусловленные порайонными особенностями землеустройства (минимум коэффициента эрозионной опасности культур, максимум оленеемкости пастбищ, максимум проектного покрытия почв растениями, минимум смыва почвы, максимум накопления в почве органического вещества и т. д.), отображаемые, как правило, в натуральном или безразмерном выражении.
При математическом моделировании экономических процессов в целом по народному хозяйству выделяют так называемые глобальный, а также отраслевой и локальные критерии оптимальности.
Глобальный критерий оптимальности — это критерий функционирования народного хозяйства как целостной экономической системы общества. В настоящее время проблема построения количественного значения глобального критерия оптимальности как математической функции цели не решена.
Математики-экономисты выделяют две основные концепции при решении этой задачи:
475
глобальный критерий — это максимум совокупного общественного продукта или важных его составных частей (национального дохода, фонда накопления, фонда потребления и т. д.);
глобальный критерий должен максимизировать благосостояние общества (его материальные и духовные потребности).
Однако сторонники и противники концепций понимают, что выбрать и обосновать глобальный критерий оптимального функционирования народнохозяйственной экономической системы формально, математическими методами невозможно. Это сложная социально-экономическая проблема.
Отраслевой критерий оптимальности характеризует эффективность отрасли или определенной сферы деятельности, каковой и является землеустройство.
Любая отрасль или сфера деятельности как подсистема народнохозяйственного комплекса имеет иерархическую структуру, элементы которой обладают определенной самостоятельностью, специфичностью, локальными целями, поэтому отраслевые критерии оптимальности могут реализовываться через локальные критерии оптимальности.
Экономические процессы, направленные на решение частных технико-экономических задач в землеустройстве, преследуют конкретные цели и оптимизируются с помощью частных критериев оптимальности, подчиняющихся требованиям локальных критериев.
Проблема эффективности землеустройства более подробно рассматривается в курсе «Экономика землеустройства».
В общем виде целевая функция линейной оптимизационной задачи записывается так:
я
2~Р(х)= ХсуХу-»тах(пип)
У = 1
или в расширенной постановке
Г(х) = С|Х) + с2х2 + ... + с„х„ —> тах (тш),
где с, — коэффициент целевой функции при переменных величинах; значение с,-в данной постановке должно быть известно.
Если задача предполагает определение оптимального плана, в котором за счет имеющихся производственных ресурсов должно быть произведено максимальное количество валовой продукции, то коэффициентами в целевой функции с^ будет стоимость валовой продукции, полученной в расчете на единицу принятой размерности для переменных величин. Оптимальное решение системы обеспечит достижение максимально возможного значения избранного показателя решения задачи, то есть максимума производства валовой продукции в стоимостном выражении.
476
Таким же образом в качестве с,- могут использоваться известные, рассчитанные заранее значения стоимости товарной продукции, валового или чистого дохода и т. д., что соответствует таким критериям оптимальности, как максимум стоимости товарной продукции, максимум валового дохода, максимум чистого дохода. То есть коэффициенты с* могут иметь как прямой характер, например стоимость валовой продукции, так и быть расчетными (производными) величинами. Такими будут коэффициенты при решении экономико-математической задачи с целевыми функциями по максимуму прибыли (где с, — прибыль, полученная в расчете на единицу размерности, принятой по объекту, обозначенному х,), минимуму приведенных затрат (где с,- рассчитывается как $ + ЕК/, $ — себестоимость продукции по у'-й переменной; 2?—нормативный коэффициент эффективности капиталовложений; А} —удельные капиталовложения на единицу размерности, принятой по переменной ху).
В последнее время для решения землеустроительных задач, имеющих природоохранный характер, используют и такой критерий оптимальности, как максимум чистого дохода в расчете на единицу приведенных затрат по каждой конкретной переменной Ху Тогда значение су- рассчитывают по следующей формуле:
с-=—^ 1 5^ЕК/
где Р/ — чистый доход по/-му объекту.
Многие землеустроительные задачи решаются по комбинированным (смешанным) критериям оптимальности, то есть по таким, которые предполагают вычисление целевой функции в процессе решения задачи.
Например, значение чистого дохода может вычисляться и так:
п Р(х)= X С'Х;-х(;->тах,
где с7 — стоимость валовой продукции на единицу вводимой переменной; х-, — сум марные производственные затраты (вычисляются в процессе решения задачи, ис ходя из ограничения ^ 5 ос,--:*,■ = 0, где 5) — себестоимость продукции на единицу переменной). ■'=1
В целевую функцию значение х-, вводится с коэффициентом (—1). Таким образом, первая часть функционала представляет собой стоимость валовой продукции, вторая — производственные затраты, а все вместе — чистый доход, который и максимизируется.
477
В ряде случаев используют критерии оптимальности, в которых целевая функция имеет дробно-линейный вид. В числе таких критериев производительность труда и т. д. Так, в качестве оценки переменных величин при названных критериях будет дробь
С -X ■
./ ./; а общая величина показателя качества решения задачи бу-
й:Х; '
дет определена как —п , которая должна достичь макси му*/
У = 1
мального значения.
Решения, которые учитывают одновременно действие нескольких критериев оптимальности, называются субоптимальными.
При решении землеустроительных задач могут применяться и другие критерии оптимальности. Например, при оптимизации структуры посевных площадей в районах водной эрозии почв предлагалось использовать такие критерии оптимальности, как минимальный суммарный коэффициент эрозионной опасности культур (в условиях эрозии, вызываемой весенним снеготаянием), а также максимум проективного покрытия почв растениями (в условиях ливневой эрозии).
Рассмотрим, как строится целевая функция с использованием первого критерия оптимальности.
Обозначим коэффициент эрозионной опасности культур через Кк. Известно, что для зяби (пара) Кк~\, для пропашных культур 0,7—0,85, для яровых зерновых 0,4—0,6, для озимых 0,3, для многолетних трав 0,01—0,06 в зависимости от года использования.
Данные коэффициенты соответствуют участкам с крутизной склона от 3° до 8° (в среднем 6°). На ровной местности опасность смыва при любом составе культур близка к 0. Поэтому в Кк вводится поправка 8Ь учитывающая крутизну склона:
где ^ — коэффициент эрозионной опасности с учетом рельефа местности; Гт — средняя крутизна склона по севообороту.
Опасность эрозионного разрушения зависит также от проти-воэрозионной устойчивости почв, поэтому в Кк1 вводится поправка 82, учитывающая устойчивость почв к смыву:
где Кп = 52 — коэффициент, учитывающий противоэрозионную устойчивость почв.
478
При известных площадях пашни или севооборотов значение Су в целевой функции выглядит следующим образом:
_ Кк5ф2
■V
где Р/— площадь пашни или севооборота.
Целевая функция примет вид
/=1у=1 /=1у=1 ^
Для удобства вычисления значений А* умножают на 10 000. Для сахарной свеклы (А* =0,85), выращиваемой в полевом севообороте площадью Р= 1000 га, со средней крутизной склона
3 3°5]=-=0,5 при противоэрозионной устойчивости почв Кп =
о = 82 = 2 значение сд определится следующим образом:
_0,85-0,5-2-10000 св 1000
Для многолетних трав в том же севообороте
0,03-0,5-210000
смн.трав- ^0д -и>->-
При решении задач на максимум с использованием коэффициента эрозионной опасности культур значение целевой функции будет выглядеть следующим образом:
/=1у=1
V Ъ
Ху—>тах.
Почвозащитное значение сельскохозяйственных культур в период ливневых дождей сказывается иначе, чем в период снеготаяния. В этом случае значительную роль играют пропашные культуры (особенно в период сентябрьских ливней) и многолетние травы, образующие растительный полог, хорошо защищающий почву от эрозии. Поэтому в условиях ливневого стока в качестве критерия оптимальности целесообразно использовать максимум проективного покрытия.
479
Обозначим проективное покрытие культур в 1-й месяц вегетации через а. Средневзвешенное проективное покрытие /-и культуры будет определяться по формуле
** = —'
где /—число месяцев вегетационного периода.
Функция цели будет выглядеть следующим образом: 8 п 8 п а,-
где Р— площадь пашни.
Значение с/у- для пропашных культур при площади пашни 1000 га и значениях проективного покрытия в мае 5%, июне— 20, июле — 50, августе — 80, сентябре — 100 % определяется так:
5+20+50+80+100 255 .,„ а,/= =——=51%.
21.2. СУБОПТИМАЛЬНЫЕ ЗЕМЛЕУСТРОИТЕЛЬНЫЕ РЕШЕНИЯ
Известно, что при землеустроительном проектировании решают многие вопросы. Во-первых, оптимизируется размещение объектов производственной и социальной инфраструктуры хозяйств; во-вторых, устанавливаются объемы производства сельскохозяйственной продукции, сочетание отраслей на предприятиях и размещение производства с учетом особенностей землеуст-раиваемой территории; в-третьих, определяются основные направления использования земель и организуется земельная площадь сельскохозяйственных предприятий. Поэтому составление землеустроительного проекта — многоцелевая задача.
Правильно составленный и экономически обоснованный проект землеустройства должен обеспечивать получение максимального количества валовой и товарной продукции, прибыли, способствовать снижению (минимизации) издержек производства, обеспечивать высокую производительность труда, низкую себестоимость продукции, а также создавать условия для постоянного повышения плодородия почв.
Качество проекта землеустройства оценивается по многим показателям. Поэтому проектное предложение, получаемое в ре-
480
зультате оптимизации только по одному критерию, может оказаться нелучшим. Использование других критериев оптимальности в отдельности создает аналогичную ситуацию, так как в каждом оптимальном плане значение выбранного в качестве целевой функции показателя экстремальное, а значения других —хуже, чем могли бы быть.
Известны примеры, когда план, максимизирующий объемы товарной продукции в стоимостном выражении, дает наибольшие издержки производства и, как следствие, наименьший чистый доход. В то же время план, максимизирующий чистый доход, может значительно снизить объем товарной или валовой продукции. Поэтому оба этих плана могут быть непригодны при составлении проекта землеустройства в реальной экономической ситуации.
В связи с этим возникает задача поиска такого решения, которое было бы наилучшим (компромиссным) по выполнению всех критериев оптимальности. В теории экономико-математических методов и моделирования такое решение называется субоптимальным.
Таким образом, субоптимальиое решение — это план, который учитывает одновременно действие всех критериев оптимальности данной задачи и отражает все реально поставленные условия, то есть субоптимальный план является или может быть неоптимальным по каждому отдельно взятому критерию, но должен быть наилучшим с точки зрения выполнения всех критериев одновременно (рис. 30).
В общей модели проекта землеустройства система ограничений представлена линейными неравенствами и уравнениями, выделяющими в евклидовом «-мерном пространстве некоторый
А Область тимыхр |
Максимум дохода |
|
|
допус- ^•**У*~^ |
|
||
|
*\ • ^\ |
|
|
|
|
>ч Субоптимальный тк план \ |
}§ Максимум валовой '\ продукции |
|
. |
Ф Эмпирический ^ЯЛ план ^>^"^ 1111 0-'^ |
|
Минш затрс |
11П |
|
► |
Рис. 30. Геометрическая интерпретация многокритериальных землеустроительных
задач
481
выпуклый многогранник. Величина этого многогранника сильно влияет на точность планов, полученных разными способами.
Если многогранник мал, любое допустимое решение будет близко по своему экономическому эффекту к оптимальному. В этом случае поиск субоптимального решения необязателен.
Если многогранник вырождается в точку, то план, найденный любым способом, будет оптимальным, поскольку он в задаче единственный.
В случае, если многогранник допустимых решений имеет значительные размеры, произвольный допустимый план, например полученный традиционными методами (эмпирический), может отличаться от оптимального по любому критерию весьма существенно. Будет он отличаться и от субоптимального плана.
При отыскании минимума целевой функции с положительными коэффициентами при неизвестных оптимум получается в вершине, достаточно близкой к началу координат, при нахождении максимума — в наиболее удаленной вершине. Точка, соответствующая субоптимальному плану, располагается в области допустимых решений где-то между минимальной и максимальной вершинами в зависимости от значимости критериев.
Учитывая это, для поиска субоптимальных решений были предложены методы, учитывающие различную предпочтительность критериев оптимальности: последовательных уступок, штрафных функций, равных и наименьших относительных отклонений, линейного мультипрограммирования, выпуклой комбинации. Обзор этих методов дан в научной работе А. М. Они-щенко (Критерии оптимизации сельскохозяйственного производства и методы нахождения наиболее эффективных планов по нескольким критериям. — Киев, 1970), а также в учебнике профессора А. М. Гатаулина (Математическое моделирование экономических процессов в сельском хозяйстве/ Под ред. А. М. Гатаулина. — М.: Агропромиздат, 1990. — С. 134—139).
Получение субоптимальных планов в экономике называют еще решением многокритериальных задач, многоцелевой оптимизацией или решением задач с векторным критерием качества.
Исследования показали, что при решении таких задач необходимы:
обоснование набора (перечня) критериев, подлежащих рассмотрению в данной модели;
оценка относительной предпочтительности критериев или построение некоторой шкалы предпочтительности;
определение условий возможного компромисса (выбор схемы компромисса) и обоснование метода нахождения компромиссного варианта решения (выбор схемы расчета обобщенного критерия).
Набор (перечень) возможных критериев определяется характером исследуемого экономического процесса и устанавливается
482
на основе логического анализа. На практике редко встречаются задачи, когда необходимо одновременно рассматривать более трех-четырех критериев.
При оценке предпочтительности различных критериев оптимальности землеустроительных задач с использованием экспертных оценок можно построить специальную шкалу. В шкале условия предпочтительности могут быть выражены в баллах оценки каждого к-то критерия из некоторого множества ^ или в виде не-
5 которых весовых коэффициентов рк; при рк > О ^рк=1.
к=\
При невозможности установить шкалу предпочтительности исходят из предположения экономической равнозначности критериев, и их ранжирование не производится:
рк=±(к=1,2,...,5).
Условия возможного компромисса определяют путем:
минимизации относительных отклонений от оптимальных значений по всем рассматриваемым критериям;
фиксирования одного из критериев на некотором заданном уровне и оптимизации по следующему критерию и т. д.
В соответствии с различными условиями компромисса разработаны методы нахождения многокритериальных компромиссных или субоптимальных решений.
Метод последовательных уступок состоит в отыскивании оптимума наиболее предпочтительного критерия, затем экстремальная величина уменьшается (или увеличивается) посредством введения в задачу нового ограничения. В расширенной задаче находится экстремум второго критерия, после чего вводится дополнительное ограничение на его величину (делается вторая уступка). В новой задаче оптимизируется третий критерий и т. д., пока все критерии не будут использованы. Метод обладает тем недостатком, что степень приближения окончательного решения к каждому отдельному оптимуму, кроме первого, остается неопределенной, и решение может оказаться ближе к экстремуму по менее важному критерию.
Метод штрафных функций заключается в том, что к основному критерию прибавляются некоторые специальные функции, сформулированные по остальным критериям. Эти так называемые штрафные функции ухудшают значение функционала тем больше, чем больше отклоняются от своих экстремальных значений учитываемые ими показатели. Решив задачу по новому критерию, получим субоптимальный план.
Если связать два критерия в одной целевой функции посред-
483
ством некоторого параметра, то компромиссный план может быть найден путем решения задачи параметрического программирования. Иногда оказывается возможным по тем же исходным данным составить дробно-линейный критерий, который также оптимизируется. Из полученной совокупности оптимальных планов выбирают наилучший по выполнению всех трех критериев. (Применение математических методов при внутрихозяйственном землеустройстве. Заключительный отчет ГИЗР. —Б. Вяземы, 1975.-С. 159-160).
Метод равных и наименьших относительных отклонений, предложенный И. Ныковским (Польша), состоит в следующем. Исходная задача решается по каждому критерию отдельно, для всех критериев вычисляют экстремальные значения. После этого ставится требование, чтобы компромиссному плану соответствовали равные и минимальные относительные отклонения всех критериев от своих экстремальных значений. Равенство отклонений обеспечивается дополнительными ограничениями, вводимыми в задачу, минимизация — новой целевой функцией. Предпочтительность критериев можно учесть путем введения поправочных коэффициентов. Решение такой «замещающей» задачи дает компромиссный план.
Способ линейного мультипрограммирования, разработанный чешским ученым И. Саской, по характеру математических операций лучше назвать способом минимакса. После отыскания оптимума для каждого рассматриваемого критерия в многограннике решений определяется точка, максимальное удаление которой от всех гиперплоскостей, соответствующих экстремальным значениям функционалов, было бы минимально. Эта точка и считается субоптимальной (Математическое моделирование экономических процессов в сельском хозяйстве/ Под ред. А. М. Гатаули-на. — М.: Агропромиздат, 1990. — С. 136—137).
По методу выпуклой комбинации, предложенному немецким исследователем X. Юттлером, сначала отыскивают оптимальные планы по каждому критерию, а затем составляют их выпуклую линейную комбинацию, коэффициенты которой определяют из решения дополнительной задачи (Юттлер X. Линейная модель с несколькими целевыми функциями// Экономика и математические методы, 1967. — № 3).
Рассмотренные методы разработаны для решения только линейных задач. Однако общая модель землеустроительного проекта, в которую целесообразно включение нескольких функционалов, содержит задачу размещения производственных центров, вызывающую разрыв почти всех целевых функций. Кроме того, некоторые важнейшие экономические показатели (например, рентабельность) имеют дробно-линейную структуру. Поэтому в Государственном научно-исследовательском институте земельных ресурсов (ГИЗР) на основе метода выпуклой комбинации в
484
1970—1975гг. И.Ф.Полуниным был разработан и апробирован для землеустроительных задач специальный способ получения субоптимальных решений без привлечений теории игр, как это имело место у X. Юттлера (Полунин И. Ф. Субоптимальные решения при обосновании землеустроительных проектов. — В сб.: Вопросы землепользования и землеустройства. Научные труды ГИЗР. - М., 1974. - № 10. - С. 56-63).
Рассмотрим порядок применения этого метода для решения землеустроительных задач.
Пусть требуется найти максимум (минимум) г показателей:
Ъс=Рк(Х1, х2, ■■; х„), к=\,2, ..., Г
при условиях
Мх[,х2, ..., х„)>0, /= 1, 2, ..., т.
Решим задачу отдельно по каждому критерию и вычислим соответствующие оптимальные планы:
Х[(хр,хр,...,х1Р), /=1,2, ...,г.
Будем считать компромиссным планом вектор Х0(Х[,х2,...,х®),
являющийся выпуклой линейной комбинацией найденных оптимальных решений:
Х§-Х\Х\ +Х2Х2 + ---+ХГХГ',
ЪМ=1, \к>0,к = 1,2,...,г.
к = \
Если область определения исходной задачи выпуклая, то такой план всегда окажется допустимым. Следовательно, данная задача заключается в определении значений X.
Используя значения X, можно вычислить все переменные субоптимального плана следующим образом:
х?=Х1х11)+Х2х[2) + ...+Хгх1г) _х^=Х1х^)+Х2х^) + ...+Хгх{2г)
х®=Х{х^+Х2х^ + ...+Хгх^
Коэффициенты Хк выпуклой комбинации И. Ф. Полунин
485
предложил определить, исходя из условия минимума наибольшего относительного отклонения величины каждого критерия в компромиссном плане от его экстремального значения. Для этого предлагалось использовать следующий путь. Подставим каждый оптимальный план Х[ в каждую целевую функцию гк. Получим численные значения функционалов:
4/=^(*Л к= Ь 2> -■->г' 1=1>2> ■■■>г>
которые при к=1 будут экстремальными. Вычислим коэффициенты аш по формуле
акГ
*№)-*№)
*№>
1к1-1кк
1кк
к=\,2,...,г, /=1,2,..., г.
Они представляют собой модули относительных отклонений величин критериев в разных оптимальных планах от их экстремальных значений. При к = 1 имеем акк=ац=0. Все коэффициенты сведены в таблицу 147.
