Анализ бифуркаций для модели SIRS с нелинейной скоростью заражения

Основная информация

• ID статьи: 2510.27292
• Название: Bifurcation analysis for a SIRS model with a nonlinear incidence rate
• Авторы: Xiaoling Wang, Kuilin Wu (Математический и статистический факультет Гуйчжоуского университета)
• Классификация: math.DS (Динамические системы)
• Дата публикации: 31 октября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2510.27292

Аннотация

В данной работе проводится глубокое исследование модели эпидемии SIRS с общей нелинейной скоростью заражения $f(I)S=\beta I(1+\upsilon I^{k-1})S$ ($k>0$). Авторы систематически анализируют существование и устойчивость точек равновесия, применяя локальную теорию бифуркаций для исследования богатого динамического поведения модели. В статье выводятся нормальные формы различных типов бифуркаций, включая бифуркацию Богданова-Такенса, бифуркацию нильпотентного фокуса и бифуркацию Хопфа. Путём вычисления первых четырёх фокусных величин определяется коразмерность бифуркации Хопфа, что раскрывает сложное динамическое явление возможного возникновения нескольких предельных циклов в системе. Численное моделирование подтверждает корректность теоретического анализа.

Исследовательский контекст и мотивация

Исследуемая проблема

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

Значимость проблемы

1. Практическое значение: Глобализация и экологические трансформации приводят к возрастающей сложности моделей передачи инфекционных заболеваний. Анализ бифуркационных явлений критически важен для прогнозирования тенденций эпидемий, разработки эффективных стратегий контроля и оптимального распределения ресурсов общественного здравоохранения.
2. Теоретическая ценность: Классическая билинейная модель скорости заражения $\beta IS$ может описывать только идеальные и простые процессы передачи, неспособна отражать сложное динамическое поведение реального мира (такое как внезапные вспышки, устойчивые эпидемии и периодические колебания).

Ограничения существующих методов

1. Ограничения простых моделей: Классические модели с билинейной скоростью заражения обычно имеют не более одной эндемической точки равновесия, неспособны представлять бистабильность или периодичность и другие сложные явления.
2. Исследования с конкретными параметрами: Предыдущие исследования скорости заражения $f(I)S = \beta I(1+\upsilon I^{k-1})S$ в основном сосредоточены на конкретных значениях $k$ (например, $k=1$ или $k=2$).
3. Неполный анализ бифуркаций: Отсутствует систематическое исследование полной структуры бифуркаций системы при общих параметрах $k>0$.

Исследовательская мотивация

Путём введения насыщения и нелинейных характеристик нелинейная скорость заражения может отражать богатое динамическое поведение, вызванное сложными факторами, такими как психологическое поведение, ограничения ресурсов и суперраспространение. Данная работа направлена на проведение комплексного анализа бифуркаций при общих параметрах $k>0$, раскрывая более богатые динамические явления, чем в предыдущих исследованиях.

Основные вклады

1. Систематический анализ точек равновесия: Для общих параметров $k>0$ систематически анализируется существование, типы и устойчивость точек равновесия, обнаруживается, что система может иметь не более 3 положительных точек равновесия.
2. Анализ бифуркаций высокой коразмерности:
   • Доказано, что система может претерпевать бифуркацию Богданова-Такенса коразмерности 3 при $1<k\leq 2$
   • Доказано, что система может претерпевать бифуркацию Богданова-Такенса коразмерности 4 при $k>2$
   • Проанализированы бифуркации нильпотентного фокуса (коразмерности 3 и 4)
3. Тонкий анализ бифуркации Хопфа:
   • Вычислены первые четыре фокусные величины
   • Доказано, что система может порождать не более 4 предельных циклов
   • Определено количество и устойчивость предельных циклов в различных областях параметров
4. Полная диаграмма бифуркаций: Построена детальная диаграмма бифуркаций, показывающая сосуществование и переходы между различными бифуркационными явлениями, такими как бифуркация седло-узел, бифуркация Хопфа и гомоклинические бифуркации.
5. Численная верификация: Предоставлены обширные результаты численного моделирования, наглядно демонстрирующие различные динамические поведения, предсказанные теорией.

Подробное описание методов

Определение задачи

Исследование модели эпидемии SIRS с нелинейной скоростью заражения:

где:

• $S(t), I(t), R(t)$: количество восприимчивых, инфицированных и выздоровевших особей
• $\beta>0$: коэффициент передачи
• $\upsilon>0, k>0$: параметры нелинейной скорости заражения
• $b>0$: коэффициент пополнения популяции
• $d>0$: коэффициент естественной смертности
• $\mu>0$: коэффициент выздоровления
• $\delta\geq 0$: коэффициент потери иммунитета

Упрощение модели

Путём анализа общей численности популяции $N(t)=S(t)+I(t)+R(t)$ получаем $\lim_{t\to\infty}N(t)=\frac{b}{d}\equiv\Lambda$.

Упрощённая двумерная система:

Путём замены переменных $I=\frac{d+\delta}{\beta}x$, $R=\frac{d+\delta}{\beta}y$, $t=\frac{1}{d+\delta}\tau$ получаем стандартизированную систему:

где $p=(d+\delta)^{k-1}\upsilon/\beta^{k-1}$, $\Lambda_0=\beta\Lambda/(d+\delta)$, $\gamma=(d+\mu)/(d+\delta)$, $\eta=\mu/(d+\delta)$.

Базовое число воспроизведения: $R_0 := \Lambda_0/\gamma$

Методология анализа точек равновесия

Точка равновесия без болезни $E_0(0,0)$:

• При $R_0<1$ является устойчивым узлом
• При $R_0>1$ является седлом
• При $R_0=1$ является седло-узлом

Эндемические точки равновесия: Определяются путём анализа нулей функции $H(x) = (1+px^{k-1})\left(1-\frac{1+\eta}{\Lambda_0}x\right) - \frac{1}{R_0}$

В зависимости от различных диапазонов значений $k$ система может иметь 0-3 положительные точки равновесия.

Методы анализа бифуркаций

1. Бифуркация седло-узел

Когда $H(x)=0$ и $H'(x)=0$, точка равновесия вырождается в седло-узел. Используя теорему о центральном многообразии и теорию нормальных форм, доказывается, что система претерпевает бифуркацию седло-узел при изменении параметров.

2. Бифуркация Богданова-Такенса

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

и проверяется условие невырожденности $\frac{\partial(\eta_1,\eta_2,\eta_3,\eta_4)}{\partial(\lambda_1,\lambda_2,\lambda_3,\lambda_4)}\neq 0$.

3. Анализ бифуркации Хопфа

Когда $\text{Tr}(J(E))=0$ и $\text{Det}(J(E))>0$, точка равновесия может быть слабым фокусом. Используется метод формальных рядов для вычисления первых четырёх фокусных величин:

и так далее.

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

Технические инновации

1. Единая аналитическая схема для общих параметров: В отличие от предыдущих исследований, сосредоточенных на конкретных значениях $k$, в работе установлена единая аналитическая схема, применимая к общим $k>0$.
2. Полная характеризация бифуркаций высокой коразмерности:
   • Впервые доказано, что модель может претерпевать бифуркацию Богданова-Такенса коразмерности 4 при $k>2$
   • Обнаружены бифуркации нильпотентного фокуса коразмерности 4
3. Точное определение количества предельных циклов: Путём вычисления четырёх фокусных величин доказано, что система может иметь не более 4 предельных циклов, что является более точным результатом, чем известные в литературе (2 предельных цикла).
4. Систематическая классификация по параметрам: В зависимости от различных диапазонов $k$ ($0<k<1$, $k=1$, $1<k<2$, $k=2$, $k>2$) предоставляется полная динамическая классификация.

Экспериментальная установка

Численные методы

• Численное моделирование проводится в MATLAB
• Используется четырёхпорядковый метод Рунге-Кутты для решения дифференциальных уравнений
• Применяется метод продолжения для отслеживания кривых бифуркаций

Выбор параметров

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

Бифуркация Богданова-Такенса ($k=2$):

• $\Lambda_0=1.8993$, $p=5.7966$, $\eta=2.3072$
• Путём изменения $\gamma$ наблюдается последовательность бифуркаций

Бифуркации высокой коразмерности ($k=3$):

• $\Lambda_0=3.1832$, $p=0.9331$, $\eta=0.6355$
• Демонстрируется явление сосуществования трёх предельных циклов

Бифуркация Хопфа:

• Множественные наборы параметров демонстрируют случаи с 1-4 предельными циклами

Методы визуализации

1. Фазовые портреты: Демонстрируют поведение траекторий в фазовом пространстве
2. Диаграммы бифуркаций: Отмечают различные кривые бифуркаций в пространстве параметров
3. Предельные циклы: Наглядно показывают сосуществование нескольких предельных циклов

Результаты экспериментов

Основные результаты

1. Существование и устойчивость точек равновесия (раздел 2)

Теорема 2.1: Свойства точки равновесия без болезни $E_0(0,0)$ полностью определяются базовым числом воспроизведения $R_0$.

Леммы 2.2-2.4: В зависимости от различных диапазонов $k$ система может иметь:

• $k\leq 1$: не более 1 положительной точки равновесия
• $1<k\leq 2$: не более 2 положительных точек равновесия
• $k>2$: не более 3 положительных точек равновесия

2. Анализ вырожденных точек равновесия (раздел 3)

Теорема 3.1: Условия существования седло-узла

• При $1<k<2$ точка $E_2^*$ является седло-узлом
• При $k=2$ точка $\bar{E}_0^*$ является седло-узлом
• При $k>2$ множественные вырожденные точки равновесия могут быть седло-узлами

Теорема 3.2: Существование бифуркации каспа При выполнении определённых условий на параметры вырожденная точка равновесия является каспом коразмерности 2.

Теорема 3.3: Касп высокой коразмерности

• $1<k\leq 2$: касп коразмерности 3
• $k>2$ и $z=\tilde{z}$: касп коразмерности 4

Теорема 3.4: Бифуркация нильпотентного фокуса При определённых условиях на параметры система может претерпевать бифуркацию нильпотентного фокуса коразмерности 3 или 4.

3. Анализ бифуркаций (раздел 4)

Теорема 4.1: Бифуркация седло-узел Система претерпевает бифуркацию седло-узел вблизи параметрических поверхностей $SN_1, SN_2, SN_3$.

Теорема 4.2: Бифуркация Богданова-Такенса ($k>2$) Система претерпевает бифуркацию Богданова-Такенса коразмерности 4 при приближении $(\Lambda_0,\gamma,\eta,p)$ к $(\tilde{\Lambda}_0,\tilde{\gamma},\tilde{\eta},\check{p})$, включающую:

• Бифуркацию каспа коразмерности 2
• Бифуркацию каспа коразмерности 3
• Бифуркацию седло-узел, бифуркацию Хопфа, гомоклинические бифуркации и т.д.

Теорема 4.3: Бифуркация Богданова-Такенса ($1<k\leq 2$) Система претерпевает бифуркацию Богданова-Такенса коразмерности 3.

4. Бифуркация Хопфа (раздел 5)

Теорема 5.1: Когда $R_i=0$ ($i=0,1,2,3$), точка равновесия является слабым фокусом не более чем второго порядка.

Теорема 5.2: Когда $l_1=0$, $\frac{3}{2}<k<2$ и $k\neq 1.9839$, точка равновесия может быть слабым фокусом третьего порядка.

Лемма 5.1: Когда $k=2$, система может иметь не более 2 предельных циклов.

Теорема 5.3: Когда $k>2$ и $R_i=0$, система может иметь не более 3 предельных циклов.

Теорема 5.4: Когда $k>2$ и $l_2=0$ или $l_4=0$, система может иметь не более 3 предельных циклов.

Замечание 5.1: Численные результаты показывают, что система может иметь 4 предельных цикла (см. рис. 5.5).

Численная верификация последовательностей бифуркаций

Рисунок 4.2 ($k=2$): Демонстрирует последовательность бифуркаций при уменьшении параметра $\gamma$:

1. Бифуркация седло-узел ($\gamma=2.6731$)
2. Бифуркация Хопфа ($\gamma=2.6717$)
3. Гомоклиническая бифуркация ($\gamma=2.6712$)
4. Два предельных цикла ($\gamma=2.6687$)

Рисунок 4.3 ($k=3$): Демонстрирует более сложную последовательность бифуркаций:

1. Бифуркация седло-узел
2. Два предельных цикла
3. Три предельных цикла
4. Бифуркация Хопфа

Численные свидетельства существования предельных циклов

Рисунок 5.1: Демонстрирует, что при $\frac{3}{2}<k<2$ и $l_1=0$ система может иметь 3 предельных цикла.

Рисунок 5.2: Демонстрирует 1-2 предельных цикла при $k=2$.

Рисунок 5.3: Демонстрирует 3 предельных цикла при $k>2$.

Рисунок 5.4: Демонстрирует 3 предельных цикла при $k>2$ и $l_2=0$.

Рисунок 5.5: Демонстрирует сосуществование 4 предельных циклов (параметры: $k=3$).

Экспериментальные находки

1. Ключевая роль параметра $k$: Диапазон значений $k$ определяет сложность динамики системы:
   • $k\leq 1$: относительно простая динамика
   • $1<k\leq 2$: появляются бифуркации коразмерности 3
   • $k>2$: появляются бифуркации коразмерности 4 и больше предельных циклов
2. Разнообразие предельных циклов: Система может демонстрировать 1-4 предельных цикла, что редко встречается в эпидемиологических моделях.
3. Богатство бифуркационных явлений: Система демонстрирует почти все классические типы бифуркаций: седло-узел, Хопф, гомоклинические, Богданова-Такенса и т.д.
4. Влияние базового числа воспроизведения: Изменение $R_0$ приводит к переходам устойчивости между точкой равновесия без болезни и эндемическими точками равновесия.

Связанные работы

Классические модели SIRS

• Kermack-McKendrick (1927): Предложили классическую схему модели SIR
• Hethcote & van den Driessche (1991): Исследовали модели с нелинейной скоростью заражения

Модели с насыщенной скоростью заражения

• Capasso & Serio (1978): Введена насыщенная скорость заражения $\frac{\beta IS}{1+\alpha I}$
• Xiao & Ruan (2007): Исследовали немонотонную скорость заражения $\frac{\beta IS}{1+\alpha I^2}$
• Zhou et al. (2007): Исследовали $\frac{kIS}{1+\beta I+\alpha I^2}$

Модели с нелинейной скоростью заражения

• Ruan & Wang (2003): Исследовали $\frac{\beta I^2 S}{1+\alpha I^2}$, обнаружили бифуркацию Богданова-Такенса коразмерности 2 и 2 предельных цикла
• Zhang et al. (2022): Исследовали общую насыщенную скорость заражения $\frac{\beta I^p S}{1+\alpha I^q}$
• Cui & Zhao (2024): Исследовали бифуркацию Богданова-Такенса при общих параметрах $p, q$

Скорость заражения, исследуемая в данной работе

• Lu et al. (2023): Исследовали случай $k=1$, обнаружили бифуркацию Богданова-Такенса коразмерности 3
• Jin et al. (2007): Исследовали случай $k=2$, обнаружили бифуркацию Богданова-Такенса коразмерности 2

Преимущества данной работы

1. Общность: Единый анализ для произвольных $k>0$, а не для конкретных значений
2. Полнота: Обнаружены бифуркации более высокой коразмерности (коразмерность 4)
3. Точность: Доказано существование не более 4 предельных циклов, что точнее известных результатов
4. Систематичность: Предоставлена полная диаграмма бифуркаций и динамическая классификация

Выводы и обсуждение

Основные выводы

1. Структура точек равновесия:
   • Система всегда имеет единственную точку равновесия без болезни $E_0$
   • При $R_0>1$ может существовать 1-3 эндемические точки равновесия
   • Количество и тип точек равновесия определяются диапазоном параметра $k$
2. Бифуркационные явления:
   • Система может претерпевать бифуркации седло-узел, Хопфа, гомоклинические
   • Существует бифуркация Богданова-Такенса коразмерности 3 ($1<k\leq 2$)
   • Существует бифуркация Богданова-Такенса коразмерности 4 ($k>2$)
   • Существуют бифуркации нильпотентного фокуса коразмерности 3 и 4
3. Предельные циклы:
   • Система может иметь не более 4 предельных циклов
   • Количество предельных циклов определяется коразмерностью бифуркации Хопфа
   • В различных областях параметров можно наблюдать 1-4 сосуществующих предельных цикла
4. Эпидемиологическое значение:
   • Нелинейная скорость заражения может порождать сложное динамическое поведение
   • Система может демонстрировать бистабильность, периодические колебания и другие явления
   • Параметр $\beta$ (коэффициент передачи) играет ключевую роль в $R_0$

Ограничения

1. Теоретическая полнота:
   • Для некоторых областей параметров (например, когда $R_i\neq 0$ и $l_1\neq 0$) полный анализ фокусных величин не завершён
   • Существование 4 предельных циклов подтверждено только численно, отсутствует строгое теоретическое доказательство
2. Практическое применение:
   • Модель предполагает, что общая численность популяции стремится к константе, что может быть неприменимо при быстрых изменениях численности
   • Не учитываются пространственная неоднородность и возрастная структура
   • Не обсуждаются практические трудности оценки параметров
3. Численная верификация:
   • Численная верификация некоторых бифуркаций высокой коразмерности является сложной
   • Отсутствует систематический анализ чувствительности к параметрам
4. Предположения модели:
   • Предполагается, что выздоровевшие в конечном итоге теряют иммунитет, что может быть неприменимо для всех заболеваний
   • Игнорируется дополнительная смертность, вызванная болезнью

Направления будущих исследований

1. Теоретическое совершенствование:
   • Завершить анализ фокусных величин для всех областей параметров
   • Строго доказать существование 4 предельных циклов
   • Исследовать глобальное динамическое поведение, включая глобальную устойчивость и глобальные бифуркации
2. Расширение модели:
   • Учесть пространственную диффузию (системы реакции-диффузии)
   • Введение эффектов запаздывания
   • Рассмотрение случайных возмущений (стохастические дифференциальные уравнения)
   • Включение возрастной структуры или неоднородности
3. Прикладные исследования:
   • Применение модели к конкретным заболеваниям (COVID-19, грипп и т.д.)
   • Исследование эффективности стратегий контроля (вакцинация, изоляция и т.д.)
   • Оценка параметров и верификация модели
4. Вычислительные методы:
   • Разработка более эффективных алгоритмов анализа бифуркаций
   • Построение полной диаграммы бифуркаций в пространстве параметров
   • Исследование хаотической и сложной динамики

Глубокая оценка

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

1. Теоретическая глубина:
   • Применены продвинутые методы теории динамических систем (теорема о центральном многообразии, теория нормальных форм, теория бифуркаций)
   • Вычислены первые четыре фокусные величины, что является технически очень сложным
   • Обнаружены бифуркации коразмерности 4, что крайне редко встречается в эпидемиологических моделях
2. Систематичность:
   • Проведён единый анализ для общих параметров $k>0$
   • Предоставлена полная классификация точек равновесия и анализ устойчивости
   • Построена детальная диаграмма бифуркаций
3. Новизна:
   • Впервые доказано, что модель может иметь 4 предельных цикла
   • Впервые обнаружена бифуркация Богданова-Такенса коразмерности 4
   • Обнаружено сосуществование множественных бифуркаций высокой коразмерности
4. Численная верификация:
   • Предоставлены высокого качества фазовые портреты и диаграммы бифуркаций
   • Численные результаты высоко согласуются с теоретическими предсказаниями
   • Визуализация ясная и легко понимаемая
5. Качество изложения:
   • Логика ясна, структура рациональна
   • Математические выводы детальны (хотя некоторые длинные выражения опущены)
   • Богатые и информативные рисунки

Недостатки

1. Теоретические пробелы:
   • Анализ фокусных величин для некоторых областей параметров не завершён (например, раздел 5.1.2)
   • Отсутствует строгое доказательство существования 4 предельных циклов
   • Недостаточен анализ глобальной динамики
2. Практическое применение:
   • Отсутствует сравнение с реальными эпидемиологическими данными
   • Не обсуждается биологическое значение параметров и практическая измеримость
   • Отсутствуют исследования применения стратегий контроля
3. Технические детали:
   • Многие промежуточные результаты вычислений опущены ("omitted for brevity")
   • Некоторые обозначения определены недостаточно ясно
   • Недостаточно деталей в описании численных методов
4. Читаемость:
   • Для неспециалистов некоторые части чрезмерно технизированы
   • Отсутствует интуитивное объяснение биологического смысла
   • Условия некоторых теорем чрезмерно сложны
5. Дизайн экспериментов:
   • Отсутствует анализ чувствительности к параметрам
   • Не предоставлены материалы для воспроизведения (код, данные)
   • Выбор параметров для некоторых численных результатов недостаточно обоснован

Влияние

1. Научный вклад:
   • Важное дополнение к теории бифуркаций эпидемиологических моделей
   • Обнаруженные явления высокой коразмерности и множественные предельные циклы имеют теоретическую ценность
   • Предоставлена новая аналитическая схема для связанных исследований
2. Практическая ценность:
   • Раскрывает сложное эпидемиологическое поведение, вызванное нелинейной скоростью заражения
   • Предоставляет теоретическую поддержку для принятия решений в области общественного здравоохранения
   • Указывает на то, что простые стратегии контроля могут быть неэффективны (из-за бистабильности и периодических колебаний)
3. Воспроизводимость:
   • Теоретические выводы детальны и могут быть проверены
   • Параметры численного моделирования ясно указаны
   • Однако отсутствуют код и подробное описание численных методов
4. Цитируемость:
   • Важный справочный материал для исследователей, работающих над нелинейными эпидемиологическими моделями
   • Методологическая ценность для математиков, изучающих бифуркации высокой коразмерности
   • Вероятно, станет важной литературой в этой области

Области применения

1. Теоретические исследования:
   • Применение теории динамических систем
   • Развитие теории бифуркаций
   • Математический анализ эпидемиологических моделей
2. Эпидемиологические приложения:
   • Заболевания со сложными механизмами передачи (например, события суперраспространения)
   • Заболевания, требующие учёта психологических факторов и поведенческих изменений
   • Заболевания, проявляющие периодические вспышки
3. Разработка политики:
   • Понимание сложной динамики передачи заболеваний
   • Оценка потенциального эффекта мер вмешательства
   • Выявление возможных критических точек и внезапных изменений
4. Образовательные цели:
   • Примеры применения теории динамических систем
   • Иллюстрация применения теории бифуркаций
   • Материалы для курсов математического моделирования

Избранные ссылки

В статье цитируются 37 важных работ, основные из которых:

Классические эпидемиологические модели:

• 28 Kermack & McKendrick (1927): Классическая модель SIR
• 18 Levin et al. (1989): Прикладная математическая экология

Нелинейная скорость заражения:

• 14 Capasso & Serio (1978): Насыщенная скорость заражения
• 15 Xiao & Ruan (2007): Немонотонная скорость заражения
• 9 Ruan & Wang (2003): Бифуркации при нелинейной скорости заражения

Теория бифуркаций:

• 32 Zhang et al. (1992): Качественная теория дифференциальных уравнений
• 33 Khibnik et al. (1998): Глобальное исследование кубического уравнения Лиенара
• 34 Dumortier et al. (2001): Развёртывание нильпотентного седла коразмерности 4

Связанные работы:

• 30 Lu et al. (2023): Исследование случая $k=1$
• 31 Jin et al. (2007): Исследование случая $k=2$
• 20-22 Zhang, Cui, Hu и др.: Исследования общей насыщенной скорости заражения

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