- •Волков Сергей Николаевич землеустройство экономико-математические методы и модели
- •Раздел 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
136. Вероятностный и детерминированный планы
Показатель |
Вероятностный план |
Детерминированный план |
Переменные величины
х{ — картофель |
10,31 |
22,14 |
х2 — зерно |
26,75 |
0 |
хъ — мясо |
2,68 |
11,62 |
х4 — осенняя капуста |
32,35 |
57,54 |
Вектор цен |
|
|
у\ — земля (период 1) |
0 |
0 |
у2 — (период 2) |
32,93 |
34,72 |
у3 — капитал (период 1) |
65,96 |
94,01 |
у4 — (период 2) |
0 |
0 |
у5 — (период 3) |
0 |
0 |
Уб — трудовые ресурсы (период 1) |
0 |
0 |
у1 — (период 2) |
0 |
0 |
у8 — (период 3) |
0 |
6,10 |
Ожидаемое значение дохода, долл. |
7209 |
9131 |
Среднее квадратическое отклонение |
2195 |
4225 |
Из таблицы видно, что при учете вероятностного характера вектора выпуска х ожидаемый чистый доход уменьшается по сравнению с детерминированным аналогом примерно на 20 %, но возрастает гарантия получения этого дохода. Коэффициент вариации для вероятностной модели Кв в = 30,4%, для детерминированной УЯВ=46,2 %. Очевидно, что для собственника земли лучше иметь меньший, но устойчивый доход, что обеспечит стабильность производства и его гарантированную эффективность.
Классификация задач и методов стохастического программирования. В настоящее время разработано большое число математических моделей задач стохастического программирования, что требует их классификации. В агроэкономических исследованиях наиболее распространенной является классификация стохастических оптимизационных моделей, разработанная Ю. И. Копенкиным, которая не претерпела существенных изменений с середины 70-х годов до настоящего времени (Крав-
412
ченко Р. Г. Математическое моделирование экономических процессов в сельском хозяйстве. — М.: Колос, 1978. — С. 407—411; Математическое моделирование экономических процессов в сельском хозяйстве/ Под ред. А. М. Гатаулина. — М.: Агропром-издат, 1990.-С. 384-389).
Согласно этой классификации стохастические модели задач делятся на два основных класса: одноэтапные и двухэтапные.
Одноэтапные задачи характеризуются тем, что принимается только одно решение (априорное), которое затем не корректируется. Это решение имеет вид детерминированного вектора:
х = (хь х2, ..., х„).
Оно не поддается корректировке при уточнении производственной ситуации. Например, априорным (предварительным) может быть решение, которое получают при усредненных значениях случайных параметров.
Двухэтапные задачи основаны на применении двух последовательно получаемых решений. На первом этапе определяется априорное (предварительное) решение, которое уточняется на втором этапе в зависимости от изменения параметров случайных условий производства (исходов).
Решение, получаемое на втором этапе, называется апостериорным. В моделях эти решения имеют вид стохастического вектора уг = (у\г, У2г, Укг)> применяемого с вероятностью рп где г— номер исходов (реализации случайных условий, г= 1, 2, ..., Л).
Термин «двухэтапная задача» имеет условный характер, так как разделение на этапы относится лишь к ее постановке. С точки зрения математического программирования это единственная задача, включающая одновременно переменные двух этапов, объединенные единой целевой функцией и системой ограничений.
В настоящее время выделяют следующие постановки одно-этапных задач стохастического программирования.
1. Примитивная постановка задачи заключается в определении оптимального плана на основе решения обычной детерминированной задачи линейного программирования, коэффициенты которой рассчитаны путем усреднения случайных параметров и считаются условно детерминированными.
Если случайным является только вектор (О-коэффициентов целевой функции, то для решения задачи достаточно заменить
коэффициенты с,- их средними значениями су. Тогда целевая функция задачи примет вид
Р(х)=М{Сх)=М(Сх)-> тах,
где М— символ математического ожидания эффекта.
413
2. Жесткая постановка задачи требует, чтобы при любых конкретных реализациях (исходах) случайных параметров выполнялись все ограничения, то есть
Агх <Д.для всех ге Я; х>0;
Р{х)=Сх-^тзх,
где г — индекс возможной реализации (исхода) случайных параметров задачи; К — множество значений г, С — вектор средних значений целевой функции.
Таким образом, решение одноэтапной стохастической задачи сводится к решению громоздкой детерминированной задачи линейного программирования при условиях:
А{х<Ву А2х<В2
Арс< Вг
А^х < 5дг х>0
и целевой функции Р(х)=Сх->тах,
где г — номер возможной комбинации значений А и В, появляющихся с некоторой вероятностью (г= 1, 2, ..., Л^.
Исходя из этого следует, что решение данной задачи в стохастической постановке может быть сведено к задаче линейного программирования с «пессимистическими» значениями, то есть наихудшими величинами параметров технико- экономических коэффициентов и объемов ограничений.
Жесткая постановка должна применяться тогда, когда нужно исключить риск. Этого требуют ситуации, когда выполнение хотя бы одного из ограничений задачи приводит к катастрофическим последствиям (например, к разорению и ликвидации сельскохозяйственного предприятия).
3. Вероятностная постановка задачи (с вероятностными ограничениями) — вероятность выполнения /-го
ограничения должна быть не менее заданной величины й}'', то есть:
г0
^ауХ]<Ь1
и=1
>А(/\0<А(/)<1, /=1,2 /я.
Задача с вероятностными ограничениями обычно решается для ситуаций, когда случайным является вектор В, а матрица А и вектор С — детерминированы.
414
Эта задача может быть также сведена к детерминированной задаче линейного программирования следующего вида:.
Р(х)=Сх-^тах, Ах<В, х>0.
Элемент случайного вектора В=Ь^1=1,2,...,т) вычисляют по уравнению вида
Ь; где ф/(6,) — плотность распределения случайной величины Ь1.
Если Р^-\, данная постановка переходит в жесткую.
Методы постановки двухэтапных задач стохастического программирования делятся на непрямые и прямые.
Непрямыми методами стараются свести зависимость Р(х, со), где со —вероятностный элемент, к зависимости Р(х), а затем применить один из известных методов математического программирования. Иными словами, стохастическую задачу сводят к ее детерминированному аналогу, а последний решают известными методами линейного или нелинейного программирования.
К прямым методам относят методы, основанные на информации только о функции Р(х, со): методы стохастической аппроксимации, проектирования стохастических квазиградиентов и др. (Ермольев Ю. М. Стохастические модели и методы оптимизаций. — Киев: Кибернетика, 1974. — № 4).
Рассмотрим рекомендуемые непрямые приемы сведения двухэтапных стохастических задач к задачам линейного программирования (Копенкин Ю. И. Стохастические задачи математического программирования. — В кн.: Математическое моделирование экономических процессов в сельском хозяйстве/ Под ред. А. М. Гатаулина. — М.: Агропромиздат, 1990. — С. 388—389).
1. Если случайным является вектор В (свободные члены ограничений Ьь обозначающие наличие ресурсов), а матрица А детерминированная, то двухэтапная стохастическая постановка допускает возникновение невязок в ограничениях (нехватку ресурсов) при отдельных исходах, однако им соответствует определенный штраф, который вычитается из целевой функции. Задача линейного программирования при этом имеет вид
415
Ах -В^
Ах +Ду1 =51
Ах +ДУ2 =^2
Ах +Оух=Вх
х>0, У1>0, ..., ум>0 Р{х,у) = Сх + р1уу1 + р2уу2 + ... + рмуум,
где х = {х\, х2, ..., х„) — вектор основных переменных (управляющих переменных первого этапа); уг = {уХп у2г, ..., утг] — вектор вспомогательных переменных (управляющих переменных второго этапа), обозначающих невязки в /-м ограничении при г-м исходе (/- = 1, 2, ..., ЛО; Вг = {Ь1п Ь2п ..., Ьтг) — вектор свободных членов ограничений, обозначающих наличие /-го ресурса при г-м исходе (г = 1, 2, ..., ЛО; А — матрица размерностью тп неизменяющихся технико-экономических коэффициентов ау, обозначающих затраты /-го ресурса за единицу у-го вида деятельности; О — вспомогательная матрица, элементы которой служат для ввода в ограничения невязок (например, диагональная матрица, элементы главной диагонали которой равны —1); У = (Уь Уъ —, Ут} — вектор, элементами которого являются размеры штрафа за единицу невязки по /-му ресурсу; рг — вероятность /--го исхода (г = 1, 2, ..., Л); С= {си с2,.,., с„] — вектор коэффициентов целевой функции; блок Ах = В0 предназначен для отображения всех условий, не зависящих от случайностей.
2. Если случайной является матрица А технико-экономических коэффициентов ау, обозначающих размеры затрат ресурсов и выход продукции на единицу ./-го вида деятельности, а вектор В детерминированный, то соответствующая задача линейного программирования также имеет блочную структуру:
Аох = В0
А\Х+ А^1 = В\ }
А2х + Б2у2 = В2
Амх+ 0мун= Вм
х>0, у,>0, ...,ук>0
Дх, У) = СХ + ряу{ +р2т + ...+Р№м,
где х = [хь хъ ..., хп) — вектор управляющих переменных первого этапа; уг= {у1п у2п •••> Утг) — вектор управляющих переменных второго этапа, соответствующих г-му исходу (г=1, 2, .,., Л7); А\, А2, ..., Ац— матрицы {тп) технико-экономических коэффициентов, соответствующие возможным исходам; Оь Л, ..., ^—вспомогательные матрицы технико-экономических коэффициентов, необходимые для отображения связи управляющих переменных г-го исхода с управляющими переменными первого этапа; В = {йь Ь2, ..., Ьт] — вектор неизменяющихся свободных членов ограничений; С= {с,, с2,..., с„} — вектор коэффициентов целевой функции при переменных первого этапа; у = {уь у2, ..., у} — вектор коэффициентов целевой функции при переменных второго этапа; д. — вероятность г-го исхода (г=\, 2, ..., IV).
Блок Л(>х:=50 предназначен для отображения всех условий, не зависящих от случайностей.
416
Второй случай имеет большое значение для тех стохастических землеустроительных моделей, технико-экономические коэффициенты которых связаны с урожайностью сельскохозяйственных культур, обусловливающей случайный характер матрицы А. К этим моделям относятся задачи оптимизации трансформации земель, определения оптимального сочетания отраслей и структуры земельных угодий, соотношения орошаемых и богарных земель и т. д. Элементы вектора В в таких задачах могут рассматриваться в большинстве случаев как детерминированные величины (площади сельскохозяйственных угодий, трудовые ресурсы, планы производства и реализации продукции). Эти элементы не зависят от конкретных результатов года по урожайности.
Если случайный характер имеют одновременно параметры матрицы А и вектора В, то задачу линейного программирования можно построить путем комбинирования описанных выше приемов.
В агроэкономических исследованиях рекомендуется использовать определенные критерии оптимальности для решения задач стохастического программирования. Эти же критерии могут применяться и для поиска оптимальных планов в землеустроительных задачах. В качестве основных целевых функций (критериев) могут быть использованы следующие.
1. Максимум или минимум математического ожидания эф фекта:
р(х) = М{Сх) -> тах (тт),
где М— символ математического ожидания.
При этом дисперсия эффекта не учитывается. Данный критерий соответствует задачам, в которых критерием оптимальности является максимум стоимости валовой и товарной продукции, максимум прибыли и т. д. В стохастических задачах показателями качества решения будут соответственно математические ожидания названных величин. Этот критерий является наиболее распространенным.
2. Минимум дисперсии эффекта:
ДСх)=1>
г „ л ЕС/*/
->тт.
Задача сводится к выбору такого плана, при котором дисперсия эффекта будет минимальной. Для многих землеустро-
417
ительных задач уменьшение колеблемости результативного показателя весьма важно. Например, желательно выбирать такой план кормопроизводства, при котором колеблемость (дисперсия) объема производимых кормов является минимальной. При использовании этого критерия оптимальности устанавливают некоторое пороговое значение математического ожидания эффекта т0, условие достижения которого вводится в модель в качестве ограничения. В этом случае требуется определить
1шп ДСх) при М(Сх) > т0.
Пороговое значение устанавливается на основе соответствующих нормативов, плановых заданий и т. д.
3. Линейная комбинация математического ожидания и дис персии линейной формы:
Сх-АхДх-> тах,
где Л — штраф за единицу дисперсии; й — квадратичная матрица, элементами которой являются дисперсия и ковариация, например выходов продукции с 1 га по культурам; Зс —транспонированный вектор х переменных задач.
При использовании данного критерия трудно объективно оценить влияние дисперсии на выбор плана, то есть определить размер штрафа за дисперсию.
4. Максимальная вероятность превышения некоторого фик сированного значения эффекта:
Р(Сх> К) -» тах,
где К— заданное пороговое значение эффекта, превышение которого желательно.
Примерами могут служить сведения к максимуму вероятности того, что прибыль будет не менее заданного значения К или объем произведенной продукции заданного вида будет не менее необходимой потребности.
Таким образом, из всего многообразия постановок задач стохастического программирования необходимо выбрать ту, которая в наибольшей степени соответствует характеру и цели исследования.
Стохастическое моделирование оптимального соотношения площадей озимых и яровых зерновых культур (одноэтапная модель).
Постановка задачи. Задача определения оптимального соотношения площадей озимых и яровых зерновых культур в агрономической науке является классической. Это объясняется тем, что озимые зерновые культуры в большинстве районов являются бо-
418
лее урожайными, чем яровые. Однако вероятность того, что озимые погибнут или дадут небольшой урожай вследствие заморозков или вымокания посевов весной, очень велика. Заменять их менее урожайными яровыми культурами часто невыгодно, так как для ряда культур (прежде всего пропашных) яровые являются менее предпочтительным предшественником и реализуются на рынке по более низким ценам, что дает хозяйству меньший доход.
В связи с этим определение оптимального соотношения площадей озимых и яровых зерновых культур — важная задача, в ходе решения которой с учетом стохастических характеристик урожайности культур необходимо получить максимальную выгоду.
Математическая формулировка задачи. Примем, что в данной задаче вероятностный характер имеют коэффициенты целевой функции — вектора с линейной формы и вектора Ь — свободных членов ограничений.
Будем считать матрицу технико-экономических коэффициентов а детерминированной.
Математическую формулировку данной одноэтапной задачи рассмотрим на примере ее постановки, данной в Государственном научно-исследовательском институте земельных ресурсов Т. Я. Перингером и И.Ф.Полуниным (Стохастическое моделирование размещения посевов озимых и яровых культур. Научный отчет «Математические методы в организации использования земель». - М.: ГИЗР, 1977. - С. 53-63).
При этом примем третий случай постановки: для задач с вероятностными ограничениями.
Определим систему переменных. Для этого примем, что ху — площадь озимой ржи, га; х2 — площадь овса, га.
В задаче необходимо определить оптимальное соотношение площадей посева двух культур: озимой ржи и овса.
На неизвестные наложим следующие ограничения.
1. Баланс площадей посевов:
Р(х, + х2=6(1))>Р0(1).
2. Необходимые трудовые затраты:
Р(а\2)х1+а^)х2<Ь{2))>Р^).
3. Объем полевых механизированных работ:
Р(а[\+а^х2<Ь^)>Р^.
419
4. Потребность в денежно-материальных средствах:
Р(а[%1+а^х2<Ь^)>Р{4\
5. Баланс минеральных удобрений:
Р(а15)х1+а(25)х2<Ь(-5))>Р^\
6. Неотрицательность переменных:
х1>0, х2>0, 0<Р0(/)<1.
Данные условия читаются следующим образом: вероятность того, что будет выполнено К-е ограничение, должна быть не
меньше Р0(/)(/=1,2,3,4,5).
Чтобы решить эту задачу при заданных вероятностных ограничениях, необходимо определить составляющие вектора Ъ.
Известно, что по закону распределения компонент вектора Ь
и вероятности Р^ можно найти составляющие вектора ограничений и от вероятностной системы ограничений перейти к ее детерминированному аналогу.
Поскольку составляющие случайного вектора ограничений Ъ независимы и каждая из них представляет собой сумму элементарных независимых слагаемых, можно считать, что компоненты вектора Ь распределены по нормальному закону, характеризуемому плотностью вероятности вида
где а — среднеквадратическое отклонение величины Ь\т — математическое ожидание величины Ь.
Графически это может быть показано так (рис. 26). Допустим, стохастическая задача решена. Тогда после подста-
п новки в левую часть неравенств вида ^а;уху<^ оптимальных
7 = 1
значений х получим некоторую величину Ь. С вероятностью Р0 можно утверждать, что Ь <Ь, так как только в этом случае требуе-420
мое условие выполняется. Причем точка_ Ъ должна быть
правее точки Ь.
Рис. 26. Плотность распределения составляющих вектора ограничений (Ь)
Чтобы по плотности распределения подсчитать вероятность попадания точки Ь в интервал \Ь,с° , надо взять интеграл от плотности в соответствующих границах:
Р0=]ДЬ)С1Ь. ь
В полученном интегральным уравнении есть только одно неизвестное Ь — такое значение свободного члена, которое позволяет перейти к детерминированному ограничению. Для определения Ь подставим в вышеприведенное уравнение вместо ЛЬ) ее выражение:
(Ь-т)1
оо | —
Не останавливаясь на выводе, приведем формулу для вычисления значения Ь, полученную Т. Я. Перингером и И. Ф. Полуниным:
Ь=т+ка,
где т — математическое ожидание (среднее арифметическое значение) величины Ь; а — квадратическое отклонение величины Ь; к— коэффициент, учитывающий изменение параметра Ь в зависимости от различной вероятности событий (табл. 137).
Зная значения ь, мы переводим стохастическую задачу к
обыкновенному детерминированному виду.
Обоснуем значение целевой функции. В связи с тем что в данной задаче важно получить гарантированный суммарный урожай зерновых (на продажу, подстилку и корм скоту), в качестве критерия оптимальности примем минимум набора урожая.
421
Исходя из рекомендуемых выше критериев оптимальности, применяемых в задачах стохастического программирования, проанализируем и выберем наиболее подходящий из них.
Пусть в качестве целевой функции выбрано математическое ожидание величины линейной формы. Тогда функцию цели математически запишем в виде
2{=М
г „ ^
7Л7
и=1
► ГП1П,
где М— знак математического ожидания; су — случайный коэффициент, учитывающий недобор урожая с 1 га посевов (/= 1, 2, ..., л).
Пользуясь правилом вынесения неслучайной величины за знак математического ожидания, а также тем, что математическое ожидание линейной функции равно той же линейной функции от математических ожиданий аргументов, функцию перепишем как
2Х=М
п \ п п п
.7=1 ) 7 = 1 7 = 1 7 = 1
Здесь С] — математическое ожидание (среднее арифметическое значение) коэффициентов линейной формы, то есть матема-
422
Г" Я3Я3Я-5.° „°Я5 Я-5.° Я-1 .Р.0 „° Р .°^° -0«°«° ЯЯ
|—ООООООООООООООООООО —
I I I I М I I I I и>р— оооооооо о о о о о о о о о- — (_^>
I—
(__)
(_^
^_; ^— (ч^ <_*-* +^ 1-П 1
"<->
^— (.
о^^и»^^
— ^'^о^о^омэ^*—-"-лол^с^-^о
Чч-^
р~~"
*4—' ^^Л *ч* I
-•—*
Ч*-* I
Г-4* 1*44 ^■*■—% / ъ I У\ .- * ^ I - « ^ 1 Г*~**\ Г*~\ Г\Л ^ \ - » - * - ■11*^1 Ъ -—^- Ль ,~*.
13 I
а? ю
3"0
1 1 I 1 I I 1 I ! I
^ ^ Г"* .Г"Я5 -°«° Я5 Я-5»° -° Я „°„° Я3 Я5 Я5-"- .г" .Г" ^ ю <?\К)Ъсо о\^л"и> го "*— о ^-1о"о^ <-/»ел "со оюомо
4^4^00(^4^---4ГО00С^ЫОК)и»1>0(О--д-р^и>004й.4^
!1
тическое ожидание недобора урожая в планируемом году, которое определяется как разность между многолетней средней урожайностью и прогнозируемой урожайностью на планируемый год.
Поясним физический и экономический смысл данного функционала (рис. 27).
Линии 1 и 2 показывают значения соответственно прогнозируемого и среднего многолетнего валового сбора озимых и яровых культур. Из рисунка видно, что в любой планируемый год может возникнуть одна из трех ситуаций: недобор урожая озимых и яровых положителен (случай 1), отрицателен (случай 2) и равен 0 (случай 3). В первых двух случаях условие задачи минимизировать величину недобора урожая озимых и яровых культур приводит к тому, что прогнозируемая кривая 1 переместится в положение 3, то есть произойдет увеличение фактического валового сбора озимых и яровых культур.
Целесообразно не только «поднять» кривую валового сбора, но и уменьшить ее колебания по годам, то есть добиться стабильности валовых сборов зерновых.
Математически величину колебаний характеризует дисперсия. Покажем, как изменится вид кривой прогнозируемого валового сбора урожая озимых и яровых культур, если показателем качества решения будет выбрана дисперсия линейной формы:
2-) = Б
и=1
► тт.
В этом выражении Б — знак дисперсии, а с,- и х] имеют тот же физический смысл, что и в целевой функции 2\.
'°ти
Й : » : •3 \
уР0.
3 <
г ;
о^""--&**•
**>
А«!
мгь
А^'
окдайнсх'
С
V Л О 1)
с
Годы
Рис. 27. Возможные пути получения урожая озимых и яровых зерновых культур
423
Учитывая правило нахождения дисперсии линейной функции, приведенное уравнение 22 можно переписать в виде
Г „ Л
ч>1 ;
г2=о
" г п2 = X [х]\ Л/+2 ^Х]Х1Ку ->1ШП,'
У=1 У'<1
где Кц — корреляционный момент величин с,- и с,-.
Последнее выражение можно записать в более удобной форме, заменив /)су на о] и Щ на /ут/Оу (г/,- — коэффициент корреляции величин с,- и с,):
%2 = X ару + 2 X ^^^■а^с^x^x^ -> тт.
Минимизируя выражение 2"2, мы уменьшаем размах колебаний, и следовательно, кривая прогнозируемого валового сбора озимых и яровых культур будет сглажена и займет положение 4 (см. рис. 27).
С экономической точки зрения наиболее выгодным представляется вариант выбора такой целевой функции, которая одновременно уменьшала бы недобор урожая и сглаживала его колебания. Этого можно достичь, если показателем качества решения задачи выступит линейная комбинация математического ожидания и дисперсии линейной формы.
Функция цели в этом случае будет состоять из двух слагаемых:
2з = (а2{ + $22) -> гшп, где а и р — некоторые коэффициенты, не равные нулю.
Преобразуем это уравнение, для чего разделим обе его части на а. Получим
|
а а |
Введем обозначения |
|
|
а а |
Тогда имеем |
|
2о = (2,\ + Х22) -» тт. 424
Выбор неизвестного коэффициента X зависит от конкретных условий задачи. В ряде исследований он называется штрафом за единицу дисперсии.
Подставляя в последнее уравнение значения 2Х и 2Ъ получим целевую функцию в окончательном виде:
,2,,2
^
Хо)х;+2Х^ст(оуХ/ху-
= 1
У<1
Для нее надо найти минимум.
Таким образом, задача стохастического программирования сведена к детерминированной задаче нелинейного программирования.
Рассмотрим пример по одному из районов Ивановской области за пять лет (исходные данные и решение Т. Я. Перингера и И. Ф. Полунина, изложение и обозначения наши).
Построение системы ограничений. Учитывая, что матрица технико-экономических коэффициентов данной задачи детерминирована, их определяют обычным способом. Например, затраты труда в расчете на 1 га посевов зерновых определяют путем деления объема затраченных трудовых ресурсов при возделывании культуры на ее общую площадь (табл. 138).
138. Исходные значения технико-экономических коэффициентов
^Ограничения Культура
Площадь посева, га
Потребность в трудовых ресурсах, чел.-ч
Объем механизированных работ, усл. эт. га
Потребность в денежно-материальных средствах, руб. на I га
Потребность в
минеральных
удобрениях,
кг на 1 га
Озимая рожь, х{ |
1 |
8,02 |
9,64 |
107,0 |
110,0 |
Овес, х2 |
1 |
8,63 |
7,25 |
82,0 |
78,0 |
По условию задачи размещения посевов озимой ржи и овса правая часть системы ограничений считается случайной, то есть
она записывается в виде Р
>о(0
и не может быть за-
,7 = 1
ранее точно определена. Рассмотрим подробно процедуру вычислений технико-экономических коэффициентов в правой части и переход от вероятностной формы ограничений к детерминированной. Из годовых отчетов за пять рассматриваемых лет находим общую площадь, занятую посевами озимой ржи и овса. За этот же период по годовым отчетам устанавливаем наличие трудовых ресурсов (с учетом привлеченной рабочей силы), техники, выделенных денежных средств и наличие минеральных удобрений. Вычисляем среднее арифметическое значение ресур-
425
са Ь;-, а также среднеквадратическое отклонение о,- по формуле
1
Исходные данные и расчеты заносим в таблицу 139. Для того чтобы свести систему вероятностных ограничений к детерминированным, зададим вероятность выполнения каждого
неравенства: например, Р^>0,$5; />0(2)>0,40; Р0(3)>0,60; ^о(4)-°>70' Р0(5)>0,35.
Тогда систему ограничений для района в вероятностной форме можно записать следующим образом:
Дх, + х2 = 15 032) > 0,85;
Р(8,02х! + 8,63х2 < 132 354) > 0,40;
Р(9,64х! + 7,25х2< 131 665) > 0,60;
Д107Х) + 82х2< 1 543 430) > 0,70;
Р(110х,+ 78х2< 1710 000) > 0,35;
Х!>0;
х2>0.
Используя вышеприведенную формулу Ь=т+ка=Ь+ка, а также данные таблицы 139, перейдем от системы вероятностных ограничений к ее детерминированному аналогу.
Например, для ограничения по трудовым ресурсам:
6=132 354; к = 0,2531 (при вероятности 0,40, табл. 139), а = 5169,
Ь =132 354+0,2531-5169=133 662. Тогда общий вид этого ограничения в детерминируемом выражении будет следующим:
8,02x1+ 8,63х2< 133 662.
Учитывая, что первое ограничение по площади посева имеет вид равенства, в детерминированный вид, исходя из закона плотности распределения (рис. 28), оно трансформируется в два условия с параметрами
^=15032-1,0366-947=14 050;
^'=15 032+1,0366-947=16 014. Тогда баланс площадей посева будет выглядеть так:
1. Х[+х2> 14 050;
2. Х!+х2<16 014.
426
