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

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

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

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

Добавлен: 31.03.2024

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

Скачиваний: 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. Электротепловые аналогии

Литература

Рис. 4.18 Двигатель постоянного тока

Рассмотрим отдельно электрическую и механическую части электрического двигателя. Входной величиной здесь является напряжение, а выходной - частота вращения вала двигателя ω. Под влиянием напряжения, через обмотку якоря протекает ток i, который, взаимодействуя с магнитным полем возбуждения Ф, создает на валу электродвигателя движущий момент:

M = c × i , (4.1)

где с - коэффициент, зависящий от конструкции двигателя.

При вращении якоря в магнитном поле в нем возникает ЭДС:

e = c × i. (4.2)

Она направлена против питающего напряжения U и поэтому вызывает уме-ньшение тока i. Коэффициент пропорциональности С в формулах (4.1) и (4.2) зависит от конструкции двигателя и силы магнитного поля Ф. Рассмотрим схему замещения якорной цепи двигателя при индуктивности якоря Lя = 0 (рис.4.19). Внешнее напряжение U уравновешивается суммой падения напряжения UR на омическом сопротивлении R якоря U и ЭДС e.

Следовательно, можно записать: U - e =UR.

Рис. 4.19. Схема замещения

Зная величину UR, можно определить ток i, а ток i определяет момент M на валу двигателя, последний, в свою очередь, связан с частотой вращения урав-нением движения:

(4.3)

Определив ω - найдем ЭДС е. Рассмотренную причинно-следственную связь физических величин представим сетью связей (рис.4.20).

Рис. 4.20

В соответствии с сетью связей, запишем физические зависимости и следующие из них передаточные функции:

Сопоставив полученные выражения для передаточных функций с сетью связей, легко построить структурную схему двигателя (рис.21).

Рис.4.21. Структурная схема двигателя

Пример 4.3

Рассмотрим регулируемый электропривод на основе двигателя постоянного тока. Требуется построить имитационную модель электропривода с учетом массы якоря электропривода. Индуктивностью якорной цепи электродвигателя пренебречь (рис.4.22).


Рис. 4.22. Регулируемый электропривод

В данной системе входное воздействие – напряжение U1- сравнивается с выходным напряжением тахогенератора (ТГ) U3. Разница этих напряжений ΔU =U1 -U3 усиливается усилителем Y, и подается на якорь двигателя постоянного тока Д. Частота вращения вала электродвигателя с помощью тахогенератора ТГ, преобразуется в пропорциональный сигнал обратной связи U3. Сеть связей рассматриваемой системы имеет следующий вид (рис. 4.23).

Рис.4.23. Сеть связей

где U2- напряжение на якоре; e – ЭДС якоря; UR – падение напряжения на сопротивлении якоря; R - сопротивление якоря; i - ток якоря; M - момент на валу; J - момент инерции якоря; ω - частота вращения; c - постоянный коэффициент; U3- напряжение тахогенератора; kТГ - коэффициент передачи тахогенератора. Структурная схема представлена на рис.4.24.

Рис.4.24. Структурная схема


11.4.5. Преобразование передаточных функций звеньев в дифференциальные уравнения в форме Коши

1) Интегрирующее звено.

Из передаточной функции интегрирующего звена

следует искомое дифференциальное уравнение:

Эквивалентная структурная схема (рис.4.25) состоит из последовательно соединенных элементарных звеньев с передаточными функциями ki и 1/p

Рис.4.25. Структурная схема

2) Апериодическое звено.

Передаточная функция звена:

откуда

или

Полученному дифференциальному уравнению соответствует эквивалентная структурная схема на рис.4.26.

Рис.4.26. Структурная схема

3) Колебательное звено.

Колебательное звено описывается дифференциальным уравнением второго порядка, которому равносильна система из двух уравнений первого порядка. Преобразовывая выражение для передаточной функции звена

получаем

Введем обозначение pYi +1 = Yi в последнее уравнение. Получим:

Система дифференциальных уравнений в форме Коши для колебательного звена имеет вид:

Эквивалентная структурная схема этого звена изображена на рис.4.27.

Рис. 4.27. Структурная схема

4) Дифференцирующее звено с замедлением.

Запишем передаточную функцию звена в виде


Эквивалентная структурная схема (рис.4.28) включает элементарное усилительное звено с передаточной функцией ki /Ti и рассматривавшееся выше апериодическое звено. Таким образом, для описания дифференцирующего звена с замедлением можно использовать уравнение для апериодического звена путем вычитания сигнала Yi из сигнала .

Рис. 4.28. Структурная схема

5) Сложное звено.

В структурной схеме исследуемой динамической системы может встретиться сложное звено с передаточной функцией

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

Введя обозначение , перепишем последнее уравнение в виде:

Раскроем скобки.

Упростив, получаем дифференциальные уравнения первого порядка в форме Коши.

С целью сокращения записи представим эти уравнения в виде:

где .

На рис. 4.29 дана эквивалентная структурная схема рассмотренного звена.

Рис. 4.29. Структурная схема

Если степень полинома в числителе передаточной функции элемента равна степени полинома в её знаменателе, то необходимо на структурной схеме системы представить этот элемент как параллельное соединение двух эле-ментов в соответствии со следующими формулами:


где ;

где .

Если степень полинома в числителе передаточной функции элемента равна степени полинома в её знаменателе, причем знаменатель имеет действительные корни, то знаменатель можно разложить на сомножители. Тогда на структурной схеме системы можно представить этот элемент как параллельное соединение двух элементов в соответствии со следующей схемой:

Требуется найти неизвестные коэффициенты, обозначенные как x и y . Преобразуем последнее уравнение.

(4.4)

При p = 0 y = b2. (4.5)

Решая уравнения (4.4) и (4.5) совместно, определяем другой неизвестный коэффициент x .

Рассмотрим второй пример.

откуда

(4.6)

При p = 0 x + y = b2. (4.7)

Решаем уравнения (4.6) и (4.7) совместно. Получаем:

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

Рис. 4.30 Функциональная схема