
- •Ю.Ю. Герасимов, в.К. Хлюстов
- •Математические методы и модели в расчетах на эвм: применение в лесоуправлении и экологии
- •Часть 1. Вариационная статистика
- •Глава 1.
- •1.1. Общие положения
- •1.2. Основные понятия статистики
- •1.3. Основы теории вероятностей
- •1.3.1. Понятие случайной величины
- •1.3.2. Классическое и статистическое определения вероятности события
- •1.3.3. Основные теоремы теории вероятностей
- •1.4. Контрольные вопросы и задания
- •Глава 2.
- •2.1. Постановка задачи
- •2.2. Классификация и группировка вариант
- •2.3. Графическое представление вариационных рядов
- •2.4.1. Показатели центральной тенденции
- •2.4.2. Показатели вариации
- •2.4.3. Достоверность статистических показателей
- •2.4.4. Показатели скошенности и крутизны
- •2.5. Доверительный интервал
- •2.6. Контрольные вопросы и задания
- •Глава 3.
- •3.1. Постановка задачи
- •3.2. Нормальное распределение
- •3.3. Логнормальное распределение
- •3.4.2. Бета-распределение
- •3.5. Распределение Пуассона
- •3.6. Семейство кривых распределения Джонсона
- •3.7. Семейство кривых Пирсона
- •Контрольные вопросы и задания
- •Глава 4.
- •4.1. Постановка задачи
- •4.3. Сравнение эмпирического распределения с теоретическим (критерий "хи-квадрат")
- •4.5. Сравнение дисперсий двух эмпирических совокупностей
- •4.6. Сравнение частот взвешенных рядов по критерию
- •4.7. Использование пакетов прикладных программ
- •4.8. Контрольные вопросы и задания
- •Глава 5.
- •5.1. Постановка задачи
- •5.2. Однофакторный комплекс
- •5.3. Двухфакторный комплекс
- •5.4. Использование ms Excel для проведения дисперсионного анализа
- •5.4.1. Однофакторный дисперсионный анализ
- •5.4.2. Двухфакторный дисперсионный анализ без повторения
- •5.5. Контрольные вопросы и задания
- •Глава 6.
- •6.1. Постановка задачи
- •6.2. Коэффициент корреляции
- •6.3. Корреляционное отношение
- •6.4. Схема полного корреляционного анализа
- •6.5. Использование пакетов прикладных программ Вычисление коэффициента корреляции с использованием ms Excel
- •Контрольные вопросы и задания
- •Глава 7.
- •7.1. Постановка задачи
- •7.2. Статистический анализ одномерных моделей
- •Уравнение прямой линии
- •Уравнение гиперболы
- •Уравнение показательной кривой
- •Окончательный выбор типа уравнения регрессии
- •7.4. Множественная регрессия
- •7.5. Применение ms Excel для расчета регрессии
- •Часть 2. Исследование операций
- •Глава 8.
- •8.1. Общие положения
- •8.2. Основные понятия системного анализа
- •8.3. Основные понятия исследования операций
- •8.4. Постановка задач принятия оптимальных решений
- •8.5. Контрольные вопросы и задания
- •Глава 9.
- •9.1. Постановка задачи
- •9.2. Графическое решение задачи линейного программирования
- •9.3. Задача линейного программирования в стандартной форме
- •Преобразования неравенств
- •Преобразование неограниченных по знаку переменных
- •2.4. Основы симплекс - метода линейного программирования
- •9.5. Метод искусственных переменных
- •9.6. Анализ чувствительности в линейном программировании
- •9.7. Решение задач линейного программирования на эвм
- •9.8. Контрольные вопросы и задания
- •Глава 10.
- •10.1. Постановка задачи
- •10.2. Метод ветвей и границ
- •10.3. Рекомендации по формулировке и решению задач цп
- •10.4. Задачи оптимизации раскроя
- •XA 0, xB 0, k 0 - целые.
- •XA 0, xB 0, k 0 - целые.
- •10.5. Постановка задачи дискретного программирования
- •Решение задач целочисленного и дискретного программирования на эвм
- •10.7. Контрольные вопросы и задания
- •Глава 11.
- •11.1. Общие понятия
- •11.2. Практические рекомендации при постановке задач динамического программирования
- •11.3. Оптимальное распределение ресурсов
- •11.4. Оптимальное управление запасами
- •11.5. Оптимальная политика замены оборудования
- •11.6. Контрольные вопросы и задания
- •Глава 12.
- •12.1. Постановка задачи
- •12.2. Применение стохастического программирования
- •12.3. Метод статистического моделирования
- •12.4. Контрольные вопросы и задания
- •Глава 13.
- •13.1. Постановка задач нелинейного программирования
- •13.2. Безусловная однопараметрическая оптимизация
- •13.2.1. Методы исключения интервалов
- •13.2.2. Методы полиномиальной аппроксимации
- •13.2.3. Методы с использованием производных
- •13.2.4. Сравнение методов безусловной однопараметрической оптимизации
- •13.3. Безусловная многопараметрическая оптимизация
- •13.3.1. Постановка задачи
- •13.3.2. Методы прямого поиска
- •13.3.3. Градиентные методы
- •13.4. Нелинейная условная оптимизация
- •13.4.1. Постановка задач условной нелинейной оптимизации
- •13.4.2. Методы штрафных функций
- •13.4.3. Методы прямого поиска
- •13.4.4. Методы линеаризации
- •13.5. Решение задач нелинейной оптимизации на эвм
- •13.6. Контрольные вопросы и задания
- •Приложение 1 Значения t - распределения Стьюдента при доверительной вероятности р и числе степеней свободы k
- •Плотность вероятности нормального распределения
- •Приложение 3 Значения χ2 при доверительной вероятности р и числе степеней свободы k
- •Продолжение приложения 3
- •Значения -функции
- •Приложение 5 Значения - в распределении Джонсона
- •Продолжение приложения 5
- •Продолжение приложения 5
- •Продолжение приложения 5
- •Приложение 6
- •Продолжение приложения 6
- •Продолжение приложения 6
- •Продолжение приложения 6
- •Приложение 7
- •Продолжение приложения 7
- •Продолжение приложения 7
- •Продолжение приложения 7
8.5. Контрольные вопросы и задания
Какие научно-технические предпосылки становления науки “Исследование операций” можно выделить?
Изобразите графически связи исследования операций с другими науками.
Приведите основные принципы системного анализа.
Приведите примеры декомпозиции систем в лесном хозяйстве и лесной промышленности.
Приведите примеры иерархии систем в лесном хозяйстве и лесной промышленности.
Дайте определения и поясните суть следующих терминов:
операция;
цель;
оптимальное решение;
показатель эффективности;
множество допустимых решений.
Глава 9.
Линейное программирование
Задачами линейного программирования (ЛП) называются задачи, в которых линейны как целевая функция, так и ограничения в виде равенств и неравенств и для которых методы математического анализа оказываются непригодными. ЛП представляет собой наиболее часто используемый метод оптимизации. В сфере лесного комплекса к их числу относятся задачи по темам:
рациональное использование сырья и материалов, оптимизация раскроя;
оптимизация производственной программы предприятий;
оптимальное размещение и концентрация производства;
составление оптимального плана перевозок, работы транспорта;
управление производственными запасами;
прочие, принадлежащие сфере оптимального планирования.
По оценкам американских экспертов около 75% от общего числа применяемых оптимизационных методов приходится на ЛП. Около четверти машинного времени, затраченного в последние годы на проведение научных исследований, было отведено решению задач ЛП и их многочисленным модификациям.
9.1. Постановка задачи
Постановка практической задачи ЛП включает следующие основные этапы: определение показателя эффективности, переменных задачи, задание линейной целевой функции W(x), подлежащей минимизации или максимизации, функциональных hk(x), gj(x) и областных xli <xi <xui ограничений.
Пример 9.1. Оптимизация размещения побочного производства лесничества
Лесничество имеет 24 га свободной земли под паром и заинтересовано извлечь из нее доход. Оно может выращивать условную растениеводческую продукцию, которая достигает товарного состояния за один год или бычков, отведя часть земли под пастбище. Растениеводческая продукция выращивается и продается партиями по 1000 кг. Требуется 1.5 га для выращивания одной партии продукции и 4 га для вскармливания одного бычка. Лесничество может потратить только 200 часов в год на свое побочное производство. Допустим, что требуется 20 часов для выращивания и заготовки одной партии продукции. Для ухода за одним бычком также требуется 20 часов. Лесничество имеет возможность израсходовать на эти цели 6 тыс. руб. Годовые издержки на одну партию растениеводческой продукции составляют 150 руб. и 1тыс. 200 руб. на одного бычка. Уже заключен контракт на поставку 2 бычков. По сложившимся ценам один килограмм продукции принесет чистый доход 2 руб. 50 коп., один бычок - 5 тыс. руб.
Постановка задачи.
1. В качестве показателя эффективности целесообразно взять доход за операцию (годовой чистый доход с земли в рублях).
2. В качестве управляемых переменных задачи следует взять:
x1 - количество откармливаемых бычков в год;
x2 - количество выращиваемой партии растениеводческой продукции по 1000 кг каждая в год.
3. Целевая функция:
5000 x1 + 2500 x2 max,
где
5000 - чистый доход от одного бычка, тыс. руб.;
2500 - чистый доход от одной партии растениеводческой продукции (1000 кг по 2 руб. 50 коп. за 1 кг).
4. Ограничения:
4.1. По использованию земли, га:
4 x1 + 1,5 x2 24.
4.2. По бюджету, руб.:
1200 x1 + 150 x2 6000.
4.3. По трудовым ресурсам, час:
20 x1 + 20 x2 200.
4.4. Обязательства по контракту, шт.:
x1 2.
4.5. Областные ограничения:
x1 0, x2 0.
Пример 9.2. Оптимизация программы рубок с учетом фактора биоразнообразия
Леса включают в себя приблизительно 12000 га и находятся на побережье. Приблизительно 2400 га в устье реки заняты парками или используются для других рекреационных целей. Остальные 9600 га, которые состоят из 10 отдельных участков, размеры которых колеблются от 110 до 3870 га., доступны для многоцелевого использования и могут подвергаться техногенному воздействию. Участки леса находятся среди земель, на которых производится как сельскохозяйственная, так и лесная продукция. Доминирующая древесная порода - сосна: 85% лесной территории покрыты сосняками или смешанными сосновыми древостоями. Остающиеся 15% покрыты смешанными лиственными древостоями. Леса, предназначенные для многоцелевого использования, разделены на 66 выделов на основе данных последнего лесоустройства. Имеются таблицы хода роста древостоев. Планируемый 60-летний период проведения лесозаготовок разделен на 3 промежутка по 20 лет каждый.
Цель управления многоцелевыми лесами заключается в максимизации выхода ликвидной древесины с учетом необходимости сохранения первозданной лесной среды. Кроме того, объем древесины, приблизительно равный расчетной лесосеке, полагается заданным. Лесозаготовки выполняются по сплошнолесосечной технологии.
Решим данную задачу на примере одного условного лесного массива площадью 360 га., состоящего из трех групп древостоев, характеристики которых представлены в табл. 9.1. Заметим, что рассматриваемая далее технология принятия решений ввиду аналогичности целей и ограничений может быть использована и для планирования лесозаготовок по другим участкам.
Таблица 9.1.
Группа древо-стоев |
Площадь, га |
Запас вырубаемой древесины в будущем через указанный период, кбм/га |
||
|
|
1-20 лет (1) |
21-40 лет (2) |
41-60 лет (3) |
1 |
100 |
3 |
10 |
30 |
2 |
200 |
12 |
17 |
20 |
3 |
60 |
25 |
20 |
18 |
1. В качестве показателя эффективности целесообразно взять годовой выход деловой древесины.
2. В качестве управляемых переменных задачи следует взять площади каждой из трех групп древостоев, каждая из которых может быть вырублена в любой из трех периодов, т.е. имеем 9 переменных xij, где i – группа древостоев под рубку, j - период вырубки:
x11 - площадь древостоев группы 1, вырубаемых в период (1), га;
x12 - площадь древостоев группы 1, вырубаемых в период (2), га;
x13 - площадь древостоев группы 1, вырубаемых в период (3), га;
x21 - площадь древостоев группы 2, вырубаемых в период (1), га;
x22 - площадь древостоев группы 2, вырубаемых в период (2), га;
x23 - площадь древостоев группы 2, вырубаемых в период (3) , га;
x31 - площадь древостоев группы 3, вырубаемых в период (1), га;
x32 - площадь древостоев группы 3, вырубаемых в период (2), га;
x33 - площадь древостоев группы 3, вырубаемых в период (3), га.
3. Целевая функция:
3 x11 + 12 x21 + 25 x31 + 10 x12 + 17 x22 + 20 x32 + 30 x13 + 20 x23 + 18 x33 max.
4. Ограничения:
4.1. По общей вырубаемой площади древостоев, га. Общая вырубаемая площадь древостоев в каждой группе за 3 периода не может превышать площадь группы древостоев. Это очевидное ограничение может быть формализовано следующим образом:
x11 + x12 + x13 100,
x21 + x22 + x23 200,
x31 + x32 + x33 60.
4.2. По равномерности выхода вырубаемой древесины, м3. То же самое требование налагается на объем древесины, который вырубается в каждом периоде. С повторяющимися циклами проблема решается получением потенциально-возможного уровня рубки:
3 x11 + 12 x21 + 25 x31 = 2000,
10 x12 + 17 x22 + 20 x32 = 2000,
30 x13 + 20 x23 + 18 x33 = 2000.
4.3. По расчетной лесосеке, га. Равная периодически вырубаемая лесосека принята во внимание с целью сохранения лесной среды и лесозаготовок. Предполагается, что120 га (360/3) - площадь, пройденная рубкой за период желательна. Это также указывает на то , что весь лесной массив должен быть пройден рубками за 60 лет. Для трех периодов это условие формулируется следующим образом:
x11 + x21 + x31 = 120,
x12 + x22 + x32 = 120,
x13 + x23 + x33 = 120.
4.4. По ограничению в сохранении окружающей среды, га. Условие принимает во внимание ограничения по лесоводственно-биологическим соображениям. Вырубаемая лесосека для каждой группы древостоев в каждый период обеспечивает сохранение подроста и жизнедеятельность обитателей лесных фитоценозов. Это требование обеспечивает каждый вырубаемый период меньше, чем вся площадь в группе древостоев, и помогает обеспечить сбалансированность распределения возрастной и пространственной структуры в лесном массиве. Максимальная вырубаемая площадь по каждой группе древостоев за каждый период: для древостоев 1 группы - 40 га, 2 группы - 90 га, 3 группы - 25 га:
x11 40, x21 90, x31 25,
x12 40, x22 90, x32 25,
x13 40, x23 90, x33 25.
4.5. Областные ограничения:
x11, x21, ... x33 0.
Пример 9.3. Оптимизация плана перевозок лесоматериалов.
Менеджер лесной компании должен решить, как снабжать их три лесозавода древесиной, заготовленной на трех лесосеках. Расстояния между лесозаводами и лесосеками приведены в табл. 9.2. Транспортные затраты на вывозку древесины лесовозами (одной модели) - 10 руб. за км. Каждый завод требует непрерывного снабжения древесиной, причем минимальное ежедневное снабжение каждого из них - 25 лесовозов. Ежедневный максимальный объем вырубаемой древесины по лесосекам (в лесовозах) следующий: первая - 25; вторая - 30; третья - 25.
Таблица 9.2.
Лесосека |
Расстояние между лесозаводами и лесосеками, км |
||
|
Лесозавод 1 |
Лесозавод 2 |
Лесозавод 3 |
1 |
8 |
15 |
50 |
2 |
10 |
17 |
20 |
3 |
30 |
25 |
15 |
Цель принятия решения - минимизация транспортных затрат. Требуется принять решение по количеству ежедневно отгружаемой древесины с лесосек к каждому лесозаводу.
Постановка задачи.
1. В качестве показателя эффективности целесообразно взять транспортные издержки (ежедневные в рублях).
2. В качестве управляемых переменных задачи следует взять число отгружаемых лесовозов в день из трех лесосек на три лесозавода, т.е. 9 переменных xij, где i - номер лесосеки, j - номер лесозавода.
3. Целевая функция:
80 x11 + 150 x21 + 500 x31 + 100 x12 + 170 x22 + 200 x32 + 300 x13 + 250 x23 + 150 x33 min.
4. Ограничения:
4.1. По вырубаемой древесине на лесосеках, лесовозы:
x11 + x12 + x13 25,
x21 + x22 + x23 30,
x31 + x32 + x33 25.
4.2. По потребности лесозаводов, лесовозы:
x11 + x21 + x31 25,
x12 + x22 + x32 25,
x13 + x23 + x33 25.
4.3. Областные ограничения:
x11, x21, ... x33 0.