- •Глава 1. Краткие сведения из элементарной математики. 6
- •Глава 2. Элементы высшей математики. 8
- •Глава 3. Применение дифференциальных уравнений первого порядка для решения прикладных задач физики, биологии, медицины 25
- •Глава 1. Краткие сведения из элементарной математики.
- •1. Функция
- •1.1. Способы задания функции
- •1.2. Основные элементарные функции.
- •2. Логарифмы и их свойства
- •Глава 2. Элементы высшей математики.
- •2.1. Механический смысл производной.
- •4. Функция нескольких переменных.
- •4.1. Частные производные.
- •4.2. Полное приращение и полный дифференциал.
- •4.3. Примеры для самостоятельной работы
- •5. Неопределенный интеграл.
- •5.1. Первообразная и неопределенный интеграл
- •5.2. Свойства неопределенного интеграла.
- •5.3.Таблица интегралов
- •5.4. Интегрирование методом замены переменного или способом подстановки
- •5.5.Интегрирование по частям
- •6. Определенный интеграл
- •6.1. Основные свойства определенного интеграла.
- •6.2. Вычисление определенного интеграла. Формула Ньютона – Лейбница.
- •6.3. Применение интегралов для решения количественных медицинских задач
- •7. Дифференциальные уравнения. Введение. Постановка задачи
- •7.1. Дифференциальные уравнения первого порядка (общие понятия)
- •7.2. Дифференциальные уравнения с разделенными и разделяющимися переменными
- •Глава 3. Применение дифференциальных уравнений первого порядка для решения прикладных задач физики, биологии, медицины
- •Глава 4. Краткое введение в теорию комплексных чисел
- •Глава 5. Дифференциальные уравнения второго порядка
- •1. Линейное однородное дифференциальное уравнение
- •1.1. Алгоритм решения дифференциального уравнения
- •1.2.Примеры решения дифференциальных уравнений
- •Глава 6. Применение дифференциальных уравнений для исследования колебательных процессов
- •1. Состояние динамических систем вблизи положения равновесия
- •2. Дифференциальное уравнение механических колебаний
- •Глава 7. Математическое моделирование в биологии и медицине
- •1. Модель Вольтерра
- •2. Фармакокинетическая модель
- •3. Простейшая математическая модель эпидемии
- •4. Простейшая модель инфекционного заболевания
3. Простейшая математическая модель эпидемии
За многие годы существования человечества огромное число людей погибло от различных эпидемий. Чума, холера, грипп и др. нередко поражали значительные массы людей. Для того, чтобы иметь возможность бороться с эпидемиями, т.е. своевременно применять те или иные медицинские мероприятия (карантины, вакцинации и т.п.),необходимо уметь сравнивать эффективность этих мероприятий. Сравнить же их можно лишь в том случае, если есть возможность предсказать, как при том или ином мероприятии будет меняться ход эпидемии, т.е. как будет меняться число заболевших.
Отсюда возникает необходимость в построении модели, которая могла бы служить целям прогноза.
Для простоты рассмотрим естественный ход эпидемии без какого-либо вмешательства и попробуем спрогнозировать последствия.
Так как нашей целью является лишь создание иллюстративной модели, то здесь мы абстрагируемся от очень многих факторов ( условия размножения бактериальных клеток, степень восприимчивости к инфекции отдельных людей, вероятность встречи носителя инфекции со здоровым человеком и т.д.)
Итак, пусть имеется Nздоровых людей, и в момент времениt=0в эту группу попадает один заболевший человек – источник инфекции. Будем предполагать, что никакого удаления заболевших из группы не происходит. Считаем также, что человек становится источником инфекции сразу же после того, как он сам заразится.
Допустим в некоторый момент времени t
X(t)– число источников инфекции,
Y(t) – число людей, могущих заболеть,
Тогда имеем X(t)+Y(t)=N+1в любой момент времени.
При t=0 выполняется условие X(0)=1
Количество новых больных dX, появившихся за промежуток времени dt, будет пропорционально числу встреч здоровых и заболевших людей, т.е. произведению величин XY. Следовательно, можно записать
или ,
где a – коэффициент пропорциональности.
Полученное дифференциальное уравнение вместе с условием X(0)=1 определяет функцию X(t), т.е. численность заболевших в момент времени t.
Найдем общее решение, предварительно разделив переменные.
Чтобы взять первый интеграл, разделим числитель и знаменатель дроби на и введем новую переменную, тогда, и отсюда.
Заменим под интегралом переменную, получим
После интегрирования имеем
, и, следовательно
С учетом того, что ,получим искомую функциюX(t)
Так как при t=0значениеX(t)=1, то для определения величиныСимеем уравнение
, откудаC=N
Окончательно получаем
Итак, мы знаем число заболевших как функцию времени. Проанализируем полученную формулу. При возрастании tзнаменатель дроби убывает, т. е.X(t)увеличивается. Это соответствует нашим представлениям, так как, согласно им, число заболевших может только увеличиваться.
Интересно выяснить, как меняется скорость увеличения числа больных. Для решения этого вопроса нужно изучить величину .
Дифференцируя, получаем
Числитель дроби обращается в нуль при .
Таким образом, когда , величина, а когда, величина.
Следовательно, функция - скорость возрастания числа заболевших – растет вплоть до момента, а затем убывает.
Этот результат хорошо согласуется с экспериментальными данными, так как известно, что в начале эпидемии число заболевших резко растет, а впоследствии скорость распространения инфекции снижается.