Об одном индикаторе, применимом для наглядной оценки быстро возрастающих функций

от автора

image
Для многих моделей эпидемий – SIR, SEIR и подобных (детали математического описания см., например, в www.idmod.org/docs/hiv/model-compartments.html) справедливо следующее утверждение: на начальном этапе эпидемии, когда количество заразившихся (I) много меньше размера популяции, скорость роста количества заболевших пропорциональна количеству заболевших:
$∂I/∂t=βI $, где β – коэффициент, характеризующий скорость заражения.

Решением данного уравнения является показательная функция. Для показательной функции $f(t)=a^t$ справедливо следующее функциональное уравнение:

$f(t+log_a⁡2 )=2f(t)$

Т.о. число $log_a⁡2 $ является периодом удвоения для функции $f(t)=a^t$. По определению, если период удвоения для гладкой неубывающей функции является константой, то функция – показательная.

Как и многие другие в это интересное время, я слежу за графиками роста заболеваемости, публикуемыми, к примеру, на сайте www.worldometers.info/coronavirus.
Уже довольно давно графики стали напоминать нечто, похожее на бумеранг или хоккейную клюшку:
image
Рисунок 1
Те же графики в логарифмическом масштабе дают чуть больше информации:

image
Рисунок 2
Видно, что темпы роста имеют тенденцию к замедлению, поскольку наклон логарифма от соответствующей функции уменьшается, но все же есть неудовлетворенность от непонимания, насколько всё-таки действенны применяемые меры по сдерживанию эпидемии.

Реальная динамика количества заражённых даже в условиях применимости приближения $∂I/∂t=βI $отличается от экспоненциальной, что обусловлено в первую очередь мерами по сдерживанию эпидемии, которые приводят к тому, что $β$ перестаёт быть константой и становится убывающей (если меры эффективные, разумеется ) функцией от времени.

В связи с вышеизложенным в качестве индикатора, применимого для наглядной оценки функций, похожих на показательную, предлагается использовать период удвоения. В общем случае для монотонно возрастающей функции $f(t)$ период удвоения $D(t)$ можно определить из следующего функционального уравнения:

$f(t+D(t))=2f(t)$

Отличие $D(t) $ от константы свидетельствует об отличии $f(t)$ от экспоненты. Применительно к динамике показателей заболеваемости, рост $D(t) $ (в идеале – к бесконечности) свидетельствует об эффективности принимаемых мер по сдерживанию эпидемии.

В случае функций, заданных таблично на дискретном множестве, например, в виде таблицы зависимости количества заболевших от даты, существует произвол в определении $D(t) $. В качестве простейшего способа определения $D(t) $ можно предложить следующий:
Пусть t∈{0;1;…;N}-дискретное время,I(t) – количество заболевших в зависимости от времени t. Тогда
image
Также можно определить «пессимистичный» период удвоения
image
«Пессимизм» в данном случае связан с тем, что сравнение I(t) производится всегда с I(o), т.е. с «низкой» по определению базой. Но мы ведь предполагаем, что с течением времени ситуация должна улучшаться? Для оптимистов есть своё определение:
image
Согласно вышеприведённым определениям image.
Ниже приведены примеры применения вышеописанного индикатора:
image
Рисунок 3
image
Рисунок 4
Данные по Испании, описанной в прессе как пример головотяпства
image
Рисунок 5
Несмотря на очевидное головотяпство на начальном этапе, Испания всё же выглядит небезнадежно.

И в заключение – родные пенаты
image
Рисунок 6
Утешает, что в Роспотребнадзоре темпы прироста заболеваемости COVID-19 в РФ сочли медленными

Файл с исходными данными, формулами и графиками можно взять здесь

Домашние задания:
1. Решить относительно $D(t)$ уравнение $f(t+D(t))=2f(t)$ для следующих функций
$f(t)=t^t$
$f(t)=Γ(t)$, где $Γ(t) $– гамма функция
$f(t)=t^n$
$f(t)=ln⁡(t)$
Также для $f(t)=ln⁡(t)$ решить уравнение $f(t+D(t))=mf(t)$
2. Ответить на вопрос: как связаны определённый выше период удвоения функции и логарифмическая производная функции?

Прошу читателей в течение недели не публиковать решения в комментариях

ссылка на оригинал статьи https://habr.com/ru/post/498020/


Комментарии

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *