Детерминированные модели

Общий обзор

Курс для магистрантов и аспирантов
Математическое моделирование, численные методы и комплексы программ
Тема 1.2

Содержание лекции

  • Определение и характеристики детерминированных моделей
  • Классификация детерминированных моделей
  • Обыкновенные дифференциальные уравнения в моделировании
  • Системы ОДУ: фазовые портреты и устойчивость
  • Уравнения в частных производных: классификация и типы
  • Краевые задачи для уравнений в частных производных
  • Интегральные уравнения в математическом моделировании
  • Примеры применения в различных областях
  • Численные методы решения

Определение детерминированных моделей

Детерминированная модель

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

Основные характеристики:

  • Отсутствие случайных факторов
  • Однозначная связь «вход-выход»
  • Воспроизводимость результатов
  • Предсказуемость поведения системы

Важно: В отличие от стохастических моделей, детерминированные не учитывают случайные воздействия и неопределенности.

Классификация детерминированных моделей

По отношению ко времени:

  • Статические — описывают состояние системы в фиксированный момент времени
  • Динамические — описывают эволюцию системы во времени

По типу переменных:

  • Дискретные — переменные принимают конечное множество значений
  • Непрерывные — переменные изменяются непрерывно

По математическому аппарату:

  • Алгебраические модели
  • Дифференциальные модели (ОДУ, УЧП)
  • Интегральные модели
  • Комбинированные модели

Обыкновенные дифференциальные уравнения

ОДУ — уравнения, содержащие неизвестную функцию одной переменной и её производные.

$$F(x, y, y', y'', ..., y^{(n)}) = 0$$

Применения:

  • Моделирование динамических процессов
  • Описание колебательных систем
  • Кинетика химических реакций
  • Модели роста популяций
  • Задачи механики и физики

Радиоактивный распад

$$\frac{dy}{dt} = -\lambda y$$
$$y(t) = y_0 e^{-\lambda t}$$

Гармонический осциллятор

$$\frac{d^2x}{dt^2} + \omega^2 x = 0$$
$$x(t) = A \cos(\omega t + \phi)$$

Системы ОДУ: фазовые портреты и устойчивость

$$\frac{dx}{dt} = f(x, y), \quad \frac{dy}{dt} = g(x, y)$$

Фазовый портрет — геометрическое представление всех возможных траекторий системы в фазовом пространстве.

Типы критических точек:

Узел (устойчивый/неустойчивый)

Все траектории сходятся к точке или расходятся от неё

Седло

Траектории сходятся по одним направлениям и расходятся по другим

Фокус (устойчивый/неустойчивый)

Траектории спирально приближаются к точке или удаляются от неё

Центр

Замкнутые траектории вокруг точки равновесия

Фазовые портреты и типы особых точек

Устойчивость определяется собственными значениями матрицы линеаризации системы в точке равновесия.

Уравнения в частных производных: классификация

УЧП — дифференциальные уравнения, содержащие неизвестные функции нескольких переменных и их частные производные.

$$F\left(x, y, u, \frac{\partial u}{\partial x}, \frac{\partial u}{\partial y}, \frac{\partial^2 u}{\partial x^2}, ...\right) = 0$$

Классификация линейных УЧП второго порядка:

$$\Delta = B^2 - 4AC$$

Эллиптический тип

Условие: $\Delta < 0$

Пример: Уравнение Лапласа: $\nabla^2 u = 0$

Применение: Стационарные процессы, задачи равновесия

Параболический тип

Условие: $\Delta = 0$

Пример: Уравнение теплопроводности: $\frac{\partial u}{\partial t} = a^2 \nabla^2 u$

Применение: Процессы диффузии, теплопроводность

Гиперболический тип

Условие: $\Delta > 0$

Пример: Волновое уравнение: $\frac{\partial^2 u}{\partial t^2} = a^2 \nabla^2 u$

Применение: Колебания, распространение волн

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

Краевые задачи для УЧП

Краевая задача — дифференциальное уравнение с заданными условиями на границе области или в начальный момент времени.

Типы краевых условий:

Краевые условия I рода (Дирихле)

$$u|_{\partial D} = \phi(x, y)$$

Задаются значения функции на границе

Краевые условия II рода (Неймана)

$$\frac{\partial u}{\partial n}\bigg|_{\partial D} = \psi(x, y)$$

Задаются значения производной по нормали на границе

Краевые условия III рода

$$\alpha u + \beta \frac{\partial u}{\partial n}\bigg|_{\partial D} = \chi(x, y)$$

Смешанные условия

Типы задач:

  • Задача Коши — задаются начальные условия
  • Краевая задача — задаются граничные условия
  • Смешанная задача — задаются начальные и граничные условия

Интегральные уравнения в моделировании

Интегральное уравнение — уравнение, в котором неизвестная функция находится под знаком интеграла.

Уравнения Фредгольма I рода

$$\int_a^b K(x,t)y(t)dt = f(x)$$

Неизвестная функция только под интегралом

Уравнения Фредгольма II рода

$$y(x) = f(x) + \lambda\int_a^b K(x,t)y(t)dt$$

Неизвестная функция внутри и вне интеграла

Уравнения Вольтерра

$$y(x) = f(x) + \lambda\int_a^x K(x,t)y(t)dt$$

Верхний предел интегрирования — переменная

Применения:

  • Задачи рассеяния в физике
  • Теория переноса излучения
  • Обращение интегральных преобразований
  • Краевые задачи математической физики

Примеры применения детерминированных моделей

Физика

  • Механика: движение материальных точек и твердых тел
  • Электродинамика: уравнения Максвелла
  • Квантовая механика: уравнение Шрёдингера
  • Термодинамика: уравнения состояния

Химическая технология

  • Кинетика химических реакций: $\frac{dk}{dt} = k_0 e^{-E/RT}$
  • Массоперенос: уравнения диффузии
  • Теплоперенос в реакторах
  • Гидродинамика химических аппаратов

Биология и медицина

  • Модели роста популяций: логистическая модель
  • Модель хищник-жертва: уравнения Лотки-Вольтерра
  • Распространение инфекций: SIR-модели
  • Фармакокинетика: модели распределения лекарств

Инженерия

  • Теория автоматического управления
  • Расчет конструкций: теория упругости
  • Гидравлика: уравнения Навье-Стокса
  • Электротехника: анализ цепей

Численные методы решения

Для ОДУ:

  • Методы Рунге-Кутта различных порядков
  • Многошаговые методы (Адамса, Гира)
  • Методы для жестких систем
  • Симплектические методы для гамильтоновых систем

Для УЧП:

  • Метод конечных разностей (МКР)
  • Метод конечных элементов (МКЭ)
  • Метод конечных объемов (МКО)
  • Спектральные методы

Для интегральных уравнений:

  • Метод последовательных приближений
  • Метод коллокации
  • Метод Галёркина
  • Квадратурные методы

Важные соображения:

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

Заключение

Ключевые выводы:

  • Детерминированные модели — основа математического описания природных и технических процессов
  • ОДУ эффективны для моделирования динамических систем с одной независимой переменной
  • УЧП необходимы для описания пространственно-распределенных процессов
  • Интегральные уравнения возникают при решении краевых задач и задач рассеяния
  • Фазовые портреты позволяют качественно исследовать поведение динамических систем

Преимущества:

  • Точность и воспроизводимость результатов
  • Возможность аналитического исследования
  • Хорошо развитые численные методы
  • Физическая интерпретируемость

Ограничения:

  • Не учитывают случайные факторы
  • Требуют точного знания параметров
  • Могут быть чувствительны к начальным условиям
  • Ограниченная применимость к сложным системам

Следующая тема:

Стохастические модели