
- •Конспект лекций
- •Владикавказ
- •Математическое моделирование элементов сложных экологических систем
- •Лекция 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. Электротепловые аналогии
- •Литература
2.3. Первая теорема подобия
У явлений, подобных в том или ином смысле (физически или математически), можно найти определенные сочетания параметров, называемые критериями подобия, которые имеют одинаковые значения.
Рассмотрим два процесса, описываемых однородными уравнениями.
Для первого процесса
Q1 + Q2 +……+ Qn = 0. (2.3)
Для второго процесса
q1 + q2 + …… + qn = 0 .4)
где
Qi = fi(P1,P2,….,Pm), i= 1…. n (2.5)
qi = fi(R1,R2,…..,Rm), i = 1…. n
В уравнениях (2.3) и (2.4) Qi и qi - однородные функции. В свою очередь P1 и R1, Р2 и R2 . . . Pm и Rm – сходственные переменные и параметры двух процессов. Поскольку Qn и qn не равны нулю, то уравнения (2.3) и (2.4) можно переписать в виде:
(2.3а)
(2.4а)
Для подобных процессов уравнения (2.3а) и (2.4а) совпадают, т.е. они тождественны. Поскольку процессы подобны, то:
;
;
……;
или
P1 = m1R1; P2 = m2R2; ……; Pm = mmRm (2.6)
После подстановки в уравнение (2.5) соотношений (2.6) вследствие однородности функции i Q можно вынести масштабы m1,m2, ….., mm в соответствующих степенях за знак функции в виде общего множителя Ni , т. е.
Qi = fi(P1,…..,Pm) = fi(m1R1,….., mmRm)= Ni fi(R1,…..,Rm) = Niqi
Поскольку
Qi = Ni qi
то имеют место равенства
-
Q1 = N1 q1
Q2 = N2 q2
Q3 = N3 q3
………..
Qn = Nn qn
(2.7)
Теперь подставим равенства (2.7) в уравнение (2.3а)
В соответствии с уравнением Фурье
N1 = N2 =…..= Nn
т.е.
В результате от уравнения (2.3а) приходим к уравнению (2.4а).
Следовательно, уравнения (2.3а) и (2.4а) оказываются тождественными. Это означает, что между соответствующими членами уравнений (2.3а) и (2.4а) существуют соотношения:
;
;
…….
Обобщая полученные результаты на S подобных процессов, получаем
где 1,2, …, S – номер процесса.
Индексы, характеризующие номер процесса, можно опустить и
записать (2.9) в общем виде
,
где idem - означает "соответственно одинаково для всех рассматриваемых процессов". Таким образом, у подобных процессов некоторые соотношения
параметров, называемые критериями подобия, оказались численно
одинаковыми. В сложных явлениях может одновременно протекать несколько различных процессов. Подобие каждого из этих процессов в отдельности обеспечивает подобие всего явления.
Обозначая критерий подобия буквой П (пи), можно дать краткую формулировку первой теоремы: для всех подобных явлений
П = idem.
Следует заметить, что справедливо и обратное положение: если критерии подобия численно одинаковы, то явления подобны. Нужно обратить внимание на практически важное свойство критериев подобия: критерии подобия любого явления могут преобразовываться в критерии другой формы, получаемые за счет операции умножения или деления ранее найденных критериев друг на друга или на какой-либо из них, т.е. если какие-либо критерии
Пk = idem и Пk+j = idem
то очевидно
Пk
Пk+j
= idem;
;
;
kПkj
= idem.
где k - любая постоянная величина.
Заметим, что уравнение (2.3а) могло иметь совершенно иную физическую природу, чем уравнение (2.4а) при одинаковой их математической записи.
Пример 2.1
Пусть имеются два тела массой М1 и М2 , движения которых подобны и описываются следующими однородными уравнениями:
Разделив первое уравнение движения на f1, а второе - соответственно на f2, получим
(2.8)
Поскольку явления подобны, то их параметры связаны соотношениями:
M1=mMM2
f1=mff2
l1=mll2
t1=mtt2
Подставив эти соотношения в уравнение (2.8), получим
Поскольку уравнение движения первого тела однородное, то
(2.9)
Подставив в (2.9) вместо масштабов отношения сходственных параметров, получим:
или
.
Обобщив полученный результат на S подобных систем, получим
Соотношение
одинаковое для всех подобных систем,
является основным критерием механического
подобия:
1
Пример 2.2
Пусть в цепи (рис. 2.3), обладающей сопротивлением R, и индуктивностью L1 при включении ее на постоянное напряжение U1 протекает процесс, описываемый уравнением
(2.10)
Рис 2.3. Электрическая схема
Во второй цепи с параметрами R2, L2 протекает подобный первому процесс, уравнение которого
(2.11)
Определим критерии подобия и покажем, что они численно одинаковы для обоих процессов.
Разделив (2.10) и
(2.11) соответственно на
и
,
получим
(2.12)
Так как процессы подобны, то
(2.13)
Подставив выражения (2.13) в уравнение (2.12), получим
Поскольку уравнение (2.12) однородное, то
Заменяя в последнем уравнении масштабы отношениями сходственных параметров, получаем
Или в критериальной записи