Файл: Математическое моделирование в экологии.doc

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 31.03.2024

Просмотров: 1124

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

СОДЕРЖАНИЕ

Конспект лекций

Владикавказ

Математическое моделирование элементов сложных экологических систем

Лекция 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. Электротепловые аналогии

Литература

Зависимость числа интервалов от объема выборки

n

До 50

50 — 100

100 — 500

500 — 1000

Свыше 1000

k

6 — 8

8-10

10-13

13-15

15-20


При этом желательно, чтобы в каждом интервале было не менее восьми наблюдений (в крайних интервалах число наблюдений может быть меньше восьми).

Для нашего примера принимаем k = 8. Определяем ширину интервала (h) по формуле

h === 15,8/8= 1,975

Принимаем h = 2,0, условно считаем Хmax = 16,0. Записываем статистический ряд по интервалам (табл. 3.3). Эмпирической функцией распределения Fn(х) называют функцию накопленных (кумулятивных) частостей, для нашего примера она имеет вид Рис. 3.3.

Функция Fn(х) может быть изображена графически (см. рис. 3.3) в координатах (х, Fn(х)). Функция Fn(х) служит оценкой неизвестной функции распре- деления F(x) для случайной величины х.

Таблица 3.3

Статистический ряд по интервалам

Номер интер- вала i

Интервал значений для хi.

Среднее значение х в ин- тервале, хi*.

Частота наблюде- ний в интервале, mi.

Частость,

wi

Накоп- ленная частость, wn,i

1

0 ≤ х1 ≤ 2

1

6

0,09375

0,09375

2

2 ≤ х2 ≤ 4

3

8

0,12500

0,21875

3

4 ≤ x3 ≤ 6

5

10

0,15625

0,37500

4

6 ≤ х4 ≤ 8

7

12

0,18750

0,56250

5

8 ≤ x5 ≤ 10

9

10

0,15625

0,71875

6

10 ≤ х6 ≤12

11

8

0,12500

0,84375

7

12 ≤ х7 ≤ 14

13

6

0,09375

0,93750

8

14 ≤ х8 ≤ 16

15

4

0,06250

1,00000

Сумма

64

1,00000

1,00000


Можно получить выборочную функцию распределения плотности вероятностей fn(х) в сере- дине интервалов хi* (см. ниже рис. 2.3)

fn(хi*)= mi /nh,

Рис. 3.3. Эмпирическая функция распределения нитратов в огурцах.

Например, для первого интервала

fn(xi*) = = 0,046875

Тогда f'n(х) для нашего примера будет иметь следующий вид:

⌠ 0,00000 при х < 0

│ 0,04687 при 0 < x1 < 2

│ 0,06250 при 2 < х2 < 4

│ 0,07823 при 4 < х3 < 6 130

fn(x) = { 0,09375 при 6 < х4 < 8

│ 0,07812 при 8 < х5 <10

│ 0,06250 при 10 < х6 < 12

│ 0,04687 при 12 < х7 <14

⌡ 0,03125 при 14 < х8 <16

Если генеральная совокупность N обладает двумерным признаком Х и Y, где Х и Y представляют собой случайные зависимые величины, то статистический ряд может иметь вид, приведенный:

в табл. 2.4 (когда обе величины Х и Y непрерывны); табл. 2.5 (Х непрерывна, а Y дискретна); табл. 2.6 (обе Х и Y дискретны).

В табл. (2.4 - 2.6) частота наблюдений mij показывает, сколько встречается в опыте пар (хi,уj).

Таблица 3.4

Статистический ряд при непрерывных величинах Х и Y

Х

Y

Границы интервала

0,1<у1<2

0,2<у2<0,3

0,3<у3<0,4

0,4<у4<0,5

Грани- цы

интер- вала

0 ≤ х1< 1,1

m11= 5

m12 = 4

m13= 3

m14= 1

1,1≤х2< 2,2

m21= 1

m22= 7

m23= 6

m24= 4

2,2 ≤x3<3,3

m31 = 0

m32 = 8

m33= 5

m34 = 2

2,2 ≤x3<3,3

m41= 2

m42= 4

m43= 1

m44= 3


Таблица 3.5

Статистический ряд, когда Х непрерывна, а Y дискретна

Y X

Значение у,.

y1= 2

y2 = 4

y3 = 6

y4 = 8

Грани-

цы интер-

вала

0,5≤ х2<1,5

m11= 4

m12=6

m13=4

m14= 2

1,5≤ х2<2,5

m21= 3

m22=5

m23=6

m24=5

2,5≤ х2<3,5

m31= 1

m32= 5

m33=4

m34=7


Таблица 3.6

Статистический ряд, когда Х и Y дискретны

У Х

Значение у,.

1

2

3

4

5

.10

m»= 5

т,,=4

m»= 8

т — 7

т =2 1з

Значе-

15

т,,=2

изб

и,,=4

т4= 3

и,,= 2

ние х,.

20

т =1 зз

т,,=6

m»= 5

и =1 з~

тзз 2

25

m4) 3

т =3 42

т,,=4

т =5 44

т,,=6




Если генеральная совокупность N обладает р-мерным признаком

(L = l, 2,...,1,..., р), где x1, x2, x3, …, xn, …, хp — случайные величины, то статистический ряд будет состоять при выборке n из n векторов

(x11, x21, x31, …, xn1)

(x12, x22, x32, …, xn2)

…………………….

(x1n, x2n, x3n, …, xnm)

где хij — случайная величина хj в i-ом испытании (наблюдении).

Результаты наблюдений могут быть также представлены в виде матрицы

X =

назад


Лекция 4. Статистические оценки параметров распределения случайных величин по выборкам

Основная задача выборочного метода заключается в том, что- бы на основе изучения выборочной совокупности получить та- кие выборочные характеристики, которые как можно более точно отражали бы характеристики генеральной совокупности. Пусть изучается некоторый признак Θ. По результатам выборки определяется оценка этого признака Θb. Как бы тщательно ни была организована выборка, будем иметь некоторую ошибку = | Θb — Θ|, которая будет отличаться от нуля. С этой точки зрения основная задача выборочного метода будет заключаться в том, чтобы величина была как можно меньше.

Чтобы выборочную оценку можно было считать, доброкачественной и пригодной для решения поставленных задач, она должна обладать определенными свойствами. Наилучшие оценки обладают такими свойствами, как несмещенность, состоятельность, эффективность и достаточность.

Выборочная оценка называется несмещенной, если ее математическое ожидание при любом объеме выборки равно значению пара- метра в генеральной совокупности, т.е.

М(Θb) — Θ = 0.

Если же имеется смещение, т.е.

М(Θb) — Θ = β,

то величина β отражает это смещение.

Пусть исследуется признак Х в генеральной совокупности, из которой сделана выборка х1, х2, ..., хn. Средняя арифметическая выборочная

является оценкой генеральной средней M(x), так как

M(x)=M

Математическое ожидание выборочной средней равно средней генеральной совокупности: M(x) = (x).

Выборочная оценка Θb параметра Θ, полученная на основе n независимых наблюдений, называется состоятельной, если предел вероятности

Pb - Θ < ) =1.