
- •Конспект лекций
- •Владикавказ
- •Математическое моделирование элементов сложных экологических систем
- •Лекция 1. Введение в моделирование. Исторический экскурс.
- •1. Основы моделирования в экологии 1.1. Общие принципы построения моделей в экологии
- •Лекция 2.
- •2.1. Элементы моделирования
- •2.2. Этапы построения математической модели
- •1.4. Элементы теории подобия, применяемые в моделировании
- •Лекция 3
- •3.2. Экологические модели
- •3.2.1. Основы экологометрики
- •3.2.2. Выборочный метод в экологометрике.
- •Зависимость числа интервалов от объема выборки
- •Статистический ряд по интервалам
- •Лекция 4. Статистические оценки параметров распределения случайных величин по выборкам
- •4.4. Статистические оценки гипотез об экологических моделях
- •Определение вариантов выборок
- •Выборка из генеральной совокупности
- •Статистическая таблица
- •Лекция 5.
- •Результаты эксперимента
- •Статистическая таблица эксперимента
- •Пример преобразования членов уравнения регрессии
- •Вычисление данных для линеаризации уравнения регрессии
- •Нормальные уравнения мнк для некоторых функций
- •Статистическое оценивание уравнения регрессии и парной корреляции.
- •Обработка результатов наблюдений
- •Лекция 6.
- •Рекомендации по выбору вида функции
- •3.4. Динамические статистические модели
- •Посадка леса
- •Данные по объему сброса качественных сточных вод
- •Данные по объему сброса сточных вод за 5-летие
- •Пример расчета 5-летних средних
- •Условное обозначение времени
- •Расчетные значения для определения уравнения динамики
- •Ряд динамики для определения сезонных колебаний
- •Лекция 7. Многофакторные эколого-математические модели. Анализ влияния отдельных факторов в экологической модели.
- •Эксперименталъный материал исследования
- •Результаты проведенных опытов
- •8.1. Анализ влияния отдельных факторов в экологической модели.
- •Лекция 9. Методы оптимизации. Метод Лагранжа
- •Лекция 10. Метод линейного программирования.
- •Лекция 11. Функциональные модели.
- •Лекция 12. Модели процессов содержащие обыкновенные дифференциальные уравнения.
- •Численные ошибки использованных для вычисления данных
- •Лекция 13. Статистические модели динамики.
- •Лекция 14. Балансовые модели.
- •Лекция 15.
- •Лекция 16. Информационные технологии в экологии. Экологические информационные системы.
- •1 6.1. Экологические информационные системы
- •1. Какова область значения для числовых характеристик?
- •Лекция 17. Использование информационных технологий для решения задач экологии.
- •Специальные приложения.
- •Значение функции
- •Значение критерия
- •Значение критерия
- •Критические значения коэффициента корреляции rk;α
- •2. Основы теории подобия
- •2.1. Подобие физических явлений и его признаки
- •2.2. Анализ размерностей
- •2.3. Первая теорема подобия
- •2.4. Применение методов подобия в математическом
- •11.3. Численные методы решения дифференциальных уравнений
- •11.3.1. Постановка задачи
- •11.3.2. Процесс численного решения
- •11.3.3. Метод Эйлера
- •11.3.4. Модифицированный метод Эйлера
- •11.3.5. Метод Рунге – Кутта
- •11.3.6. Метод Рунге – Кутта для систем дифференциальных уравнений
- •11.3.7. Общая характеристика одношаговых методов
- •3.8. Многошаговые методы
- •11.3.9. Методы прогноза и коррекции
- •11.3.10. Краткая характеристика методов прогноза и коррекции.
- •11.3.11. Выбор шага и погрешность решения.
- •11.3.12. Жесткие задачи
- •11.4. Имитационное моделирование систем
- •11.4.1. Принципы имитационного моделирования
- •11.4.2. Объекты моделирования
- •11.4.3. Динамическая модель исследуемого объекта
- •11.4.4. Построение имитационных моделей динамических систем
- •11.4.5. Преобразование передаточных функций звеньев в дифференциальные уравнения в форме Коши
- •11.4.6. Синтез имитационной модели на основе структурной схемы
- •11.5. Теоретические основы построения математических моделей систем
- •11.5.1. Компонентные и топологические уравнения
- •11.5.2. Компонентные и топологические уравнения механической системы
- •11.5.3. Компонентные и топологические уравнения электрической системы
- •11.5.4. Компонентные и топологические уравнения гидравлической системы
- •11.5.5. Компонентные и топологические уравнения тепловой системы
- •11.6. Метод электроаналогий
- •11.6.1. Сущность метода электроаналогий.
- •11.6.2. Электромеханические аналогии
- •11.6.3. Построение имитационных моделей методом электроаналогий
- •11.6.4. Плоское прямолинейное движение звеньев
- •11.6.5. Электрогидравлические аналогии
- •11.6.6. Электротепловые аналогии
- •Литература
Лекция 15.
Балансовые модели (продолжение)
для
каждого
s
S,
0
<
α
<
1
(α
—
риск),
где
f(∆R-,
∆X)
—
ожидаемые
потери,
определяемые оптимальной
стратегией
на
бесконечном
плановом
периоде
при
заданном текущем
состоянии
S.
математического программирования.
Рис. 15. 13. Схема материальных потоков
В другом случае используется модель межотраслевого баланса с учетом затрат на ликвидацию загрязнений. Здесь применяется модель Леонтьева — Форда.
Модель показывает, каким образом материальные потоки, поступающие в экономико-экологическую систему, распределяются затем по разным видам деятельности. Все, что поступает в систему в виде сырья и материалов, либо преобразуется в готовые изделия, либо идет в отходы производства. Представление об объемах материальных потоков, поступающих в окружающую среду, дает следующая условная схема (рис. 15.13), составленная для города с населением в 1 млн. человек. На схеме дано изображение трех основных входных потоков (вода, пища и топливо) и трех выходных потоков (сточные воды, твердые отходы и загрязнения воздуха), которые являются общими для всех городов. В этой модели появляются величины, измеренные в натуральных единицах, а именно отходы производства по каждому виду загрязнителей. Это обстоятельство существенно меняет привычные свойства модели межотраслевого баланса, в которой все величины выражены в стоимостной форме.
Основные условия модели Леонтьева — Форда в матрично-векторной форме можно записать следующим образом
,
(5.38)
где Х1 - вектор валовых выпусков продукции размерности т,
Х2 - вектор объемов загрязнений, подлежащих ликвидации, размерности n - т,
,
;
Y1 — вектор конечной продукции размерности n;
Y2— вектор объемов загрязнений, которые в настоящее время не могут быть ликвидированы (например, из-за нехватки средств), размерности n - т,
,
;
А11= (аij)n-m — матрица прямых затрат (матрица А в первоначальной модели межотраслевого баланса, которая имеет вид
где аij, — определяет какой объем i-гo ресурса необходим для производства единицы продукции j-й отрасли);
А12=(аiq)m(n-m) — матрица прямых затрат продукта i на уничто- жение единицы загрязнения вида q,
А21=(аki)m(n-m) — матрица коэффициентов, характеризующих количество поступающих в окружающую среду отходов по каждо- му виду загрязнителей k в расчете на единицу валового выпуска продукции j каждой из отраслей;
А22=(аkq)(n-m)(n-m)— матрица коэффициентов выброса загряз- нений k-гo вида при уничтожении единицы загрязнения вида q, т.е. учитывает вторичный эффект загрязнений.
Перепишем модель Леонтьева — Форда (5.38) в виде:
X1=A11X1+A12X2+Y1; (5.39)
X2=A21X1+A22X2+Y2. (5.40)
С помощью этой модели может быть решена задача определе- ния валовых выпусков продукции отраслей с учетом затрат на лик- видацию загрязнений, которая была ранее описана, но лишь с пред- положением, что отходы по каждому виду загрязнителей пропорци- ональны валовым выпускам продукции отраслей, а затраты на лик- видацию загрязнений пропорциональны объемам загрязнений, под- лежащих ликвидации. В этой модели учтены только производственные загрязнения A12 и вторичные загрязнения А22. Вектор Х2+Y2 характеризует общие объемы отходов по каждому виду загрязне- ний, образовавшихся в течение года в результате производственной деятельности. Загрязнения, которые образуются в сфере конечного потребления, могут быть отражены в модели также, как это сделано для производственных загрязнений. В системе уравнений (5.39 - 5.40) неизвестными являются вектор валовых выпусков продукции отраслей Х1 и вектор подлежащих ликвидации объемов загрязнений Х2. Они могут быть найдены, если заданы вектор конечных выпусков продукции отраслей Y1 и вектор объемов загрязнений, которые в на- стоящее время не могут быть ликвидированы Y2. Но сразу по данной системе уравнений этого сделать нельзя. Преобразуем эту систему.
Если пренебречь вторичными загрязнениями, связанными с деятельностью предприятий, ликвидирующих загрязнения, т.е. если считать А22 = 0, то получим:
X1=A11X1+A12X2+Y1;
X2=A21X1+0 - Y2.
Подставив в первое уравнение значение Х2 получим:
X1=A11X1+A12(A21X1-Y2)+Y1;
X1=A11X1+A12A21X1-A12Y2+Y1;
(I-A11-A12A21)X1 = -A12Y2+Y1;
(I-A11-A12A21)-1(I-A11-A12A21)X1 = (I-A11-A12A21)-1(-A12Y2+Y1).
Пусть В11= (I - А11 - А12А21)-1 и В12 = В11А12,
тогда
X1 = B11Y1-B12Y2;
X2 = A21(B11Y1-B12Y2)-Y2;
X2=A21B11Y1-(A21B12+I)Y2.
Обозначим В21 = А21В11 и В22 = А21В 12+ I,
тогда
Х2 = В21Y1 — В22Y2.
Итак, мы получили следующую систему уравнений:
X1 = B11Y1 – B12Y2; (5.41)
X2 = B21Y1 – B22Y2. (5.42)
где для матриц В11, В12, В21 и В22 имеем:
В11= (I — А11 — А12А21)-1; В12 = В11А12;
В21 = А21В11; В22 = А21В12 +I,
где
.
В решении (5.41 — 5.42) вектор (В11Y1) характеризует валовые вы- пуски продукции отраслей при условии ликвидации всех загрязнений, образовавшихся в течение данного года в результате производственной деятельности (здесь и далее полагаем А22 = 0). Иначе (В11Y1) — это вектор валовых выпусков продукции отраслей в случае, если бы все образующиеся в течение года загрязнения были бы ликвидированы; (В12Y2) — вектор потенциальных затрат промежуточной продукции. Промежуточная продукция — это часть всей продукции, представляющая закупки данного вида продукции отраслями-потребителями в качестве исходных материалов для производства их продукции на ликвидацию не ликвидируемых загрязнений. Этот вектор является частичной характеристикой обычно не учитываемых издержек производства. (B21Y1) — вектор потенциальных отходов производства в случае ликвидации всех загрязнений.
Таким образом, матрица В11 характеризует затраты промежуточ- ной продукции в расчете на единицу конечного выпуска продукции каждой из отраслей с учетом ликвидации всех загрязнений (матрица полных затрат), а матрица В21 характеризует количество отходов в расчете на единицу конечного выпуска продукции каждой из отрас- лей с учетом всех загрязнений.
Модель межотраслевого баланса с учетом затрат на ликвидацию загрязнений иллюстрирует, насколько сложными становятся хозяй- ственные взаимосвязи. Сокращение не ликвидируемых загрязнений Y2 приведет к росту объемов загрязнений, подлежащих ликвидации Х2, а это в свою очередь вызовет рост расходов на ликвидацию загрязнений (А12Х2) и, следовательно, приведет к росту валовых выпусков продукции отраслей Х1,что вызовет увеличение объемов загрязнений (А21Х1).
Уравнения межотраслевого баланса с учетом затрат на ликвидацию загрязнений (5.39 — 5.40) могут быть дополнены ограничения- ми по ресурсам: во-первых на трудовые ресурсы и во-вторых на факторы производства, объем вовлечения которых ограничен.
Пусть имеется один вид ограниченного ресурса — трудовые ресурсы. Обозначим вектор прямых затрат труда на единицу валового выпуска продукции (затрат труда в отрасли в расчете на единицу ее валового выпуска) через l1. и вектор прямых затрат труда на ликви- дацию загрязнений через l2. Обозначим далее через L ресурсы труда, тогда
L1Х1 + 12Х2 = L.
Подставив значения Х1 и Х2, получим:
q1Y1 – q2Y2 = L
где q1 = l1В11 +l2В21 и q2=l1В12+l2В22.
Вектор q1 состоит из компонент, характеризующих затраты труда на производ-ство единицы конечной продукции Y1 при ус- ловии ликвидации всех загрязнений, образовавшихся в данном году. Аналогично q2 — вектор затрат труда (потенциальных) на ликвидацию единицы не ликвидируемых загрязнений. Уменьшение, Y2 делает баланс трудовых и прочих ресурсов более напряженным, а его увеличение — менее напряженным.
Следует, однако, помнить, что помимо ресурсов, воспроизводимых внутри системы и являющихся обычными продуктами производства, существуют также ресурсы, которые вводятся в систему извне, и их максимальный объем вовлечения ограничен (так называемые ограниченные факторы производства). Характерно, что уро- вень деятельности ограничен не только трудом, но в зависимости от выбора продолжительности периода производства также и основными фондами, главными составными элементами которых являются производственные здания и станки, а также землей и многими другими важными ресурсами. Если через В обозначить матриц„ коэффициентов затрат ограниченных ресурсов на единицу валового выпуска продукции каждой из отраслей и через b — вектор ограни- ченных ресурсов, то реально достижимый объем выпуска продук- ции должен отвечать следующему условию
ВХ1= b.
Наиболее упрощенной модификацией модели Леонтьева — Фор- да является задача нахождения Х1 и Y2 при заданных Y1 и Х2. В этом случае вектор валовых выпусков Х1 равен:
Х1= (I-A11)-1(А12Х2+Y1)
и, следовательно, вектор не ликвидируемых объемов загрязнений Y равен
Y2=A21X1+A22X2 – X2
Чем больше воздействие производства на окружающую среду и больше средств выделяется для природоохранных целей, тем больше становятся коэффициенты матрицы (А12А21) и, следовательно, коэффициенты полных затрат В11. В крайнем предельном случае ко- эффициенты полных затрат могут стать сколь угодно большими. Рассмотрение этого крайнего, но теоретически допустимого случая полезно, так как в различных методах описания межотраслевых взаимосвязей он может давать в итогах значительные расхождения.
Использование модели Леонтьева — Форда для вариантных рас- четов позволяет получить информацию на макроуровне относитель- но отраслевой структуры затрат на охрану окружающей среды, вли- яния их на величину конечного или общего выпуска, изменения цен в зависимости от предлагаемого уровня загрязнения среды и других показателей.
Пример.
В
регионе
действуют
два
предприятия.
Они
выпускают
промежуточную
и
конечную
продукцию.
Матрица
прямых
затрат
на
производ-ство
этой
продукции
А11
имеет
вид
.
Производство
продукции
двумя
отраслями
сопровождается
выбросом
загрязнений.
Матрица
коэффициентов,
характеризующих
количество
поступающих в
окружающую
среду
отходов,
в
расчете
на
единицу
валового
выпуска продукции
каждой
из
отраслей
А21
имеет
вид
.Предприятия
создают
очистные
сооружения
для
уничтожения
загрязнений. Матрица
прямых
затрат
предприятий
на
уничтожение
единицы
загрязнения
А12
имеет
вид
.
Характерно,
что
в
результате
деятельности
предприятия
по
уничтожению
заг-
рязнений
от
производства продукции
загрязнения
не
выделяются.
Определить
валовый
выпуск продукции
предприятий
Х1
и
объем
загрязнений, подлежащих
ликвидации,
Х2,
если
известно,
что
потребность
региона в
конечной
продукции
предприятий
Y1
имеет
вид
,а
объем
загрязнений
Y2,
который
в
настоящий
момент
не может
быть
уничтожен
—
.
Решение.
X1 = B11Y1 – B12Y2;
где В11= (I — А11 — А12А21)-1, а В12 = В11А12.
А12А21=.
I —
А11
—
А12А21=.
det(I — А11 — А12А21)=0,76x0,67-0,57x0,59=0,1729.
(I — А11 — А12А21)(1,1) = 0,67 (I — А11 — А12А21)(1,2) = 0,57
(I — А11 — А12А21)(2,1) = 0,59 (I — А11 — А12А21)(2,2) = 0,76
(I
—
А11
—
А12А21V)T
=
.
B11
=
(I
—
А11
—
А12А21)-1
=
.
B12
= B11
А12
=.
X1
= B11Y1
– B12Y2
=
.
X2=B21Y1+B22Y2,
где В21 = А21В11; В22 = А21В12 +I.
В21
=
А21В11
=
.
В22
=
А21В12
+I
=
.
X2=B21Y1+B22Y2
=
.
Ответ: X1
=
;
X2
=
.
назад