Периодические колебания встречаются в популяционной динамике, химических реакциях, климате и эпидемиологии. Детерминированные модели часто дают устойчивый предельный цикл, но реальные системы подвержены шуму. Шум может сдвигать частоту, расширять спектр и разрушать временную когерентность. Существующие подходы либо слишком частные, либо не дают явных формул. Данная работа заполняет пробел, предлагая общий формализм для анализа зашумлённых предельных циклов. Она стоит рядом с исследованиями квазициклов и стохастических осцилляторов, но отличается универсальностью и применимостью к произвольным системам.
Рассматривается автономная система обыкновенных дифференциальных уравнений с устойчивым предельным циклом, возмущённая аддитивным гауссовым белым шумом. Траектория раскладывается на продольную компоненту вдоль цикла и поперечную компоненту в гиперплоскости, перпендикулярной касательной. Вводится сопутствующий ортонормированный базис, вращающийся вместе с циклом. Линеаризуя динамику вблизи цикла, получают стохастические дифференциальные уравнения для фазы и отклонения. Для случая нормальной формы суперкритической бифуркации Хопфа выводятся явные формулы автоковариации и спектральной плотности. Проверка проводится численным моделированием методом Рунге-Кутты высокого порядка.
Фаза вдоль цикла продвигается с детерминированной скоростью плюс случайный дрейф, похожий на броуновское движение. Отклонение от цикла описывается (n-1)-мерным процессом Орнштейна-Уленбека, который затухает со скоростью, определяемой локальной устойчивостью цикла. Для нормальной формы Хопфа получены явные формулы: автоковариация экспоненциально затухает с декрементом, зависящим только от радиуса цикла и интенсивности шума, а спектр мощности имеет пик, сдвинутый в сторону более высоких частот. Численные симуляции подтверждают точность приближений при слабом шуме (отношение шум/сигнал меньше 0.1) и близких частотах фокуса и цикла. Применение к данным Эль-Ниньо даёт оценку периода около 4.2 года и ляпуновский показатель 0.15-0.17 в год, что согласуется с предыдущими оценками.
Работа даёт универсальный инструмент для анализа временных рядов, порождённых зашумлёнными предельными циклами. Полученные формулы позволяют оценивать параметры системы непосредственно по данным, различая вклад детерминированной динамики и шума. Это отличает подход от моделей с шумоподдерживаемыми осцилляциями, где когерентность возникает только благодаря шуму. Результаты применимы к широкому классу моделей в экологии, эпидемиологии и климатологии. В отличие от соседних работ по квазициклам, здесь рассматриваются истинные предельные циклы, а не затухающие колебания.
Можно обобщить метод на случай неаддитивного шума, а также на системы с несколькими предельными циклами. Перспективно исследовать влияние шума на синхронизацию связанных осцилляторов. Возможно расширение на стохастические дифференциальные уравнения с запаздыванием, часто используемые в моделях климата. Развитие методов подгонки формул к данным позволит применять подход к более широкому кругу временных рядов.
Метод может быть полезен в климатологии для анализа Эль-Ниньо, в экологии для популяционных циклов, в эпидемиологии для периодических вспышек заболеваний.
Следующим шагом может быть применение формул к другим климатическим индексам, таким как Североатлантическое колебание. Также стоит проверить метод на лабораторных данных химических осцилляторов.
В динамических системах остаётся открытым вопрос о том, как шум влияет на бифуркации и устойчивость предельных циклов. Работа вносит вклад в понимание стохастической устойчивости и когерентности колебаний, что связано с проблемой предсказуемости хаотических и квазипериодических систем.
Эль-Ниньо, несмотря на свою нерегулярность, имеет характерный период около 4 лет, что близко к оценке, полученной в работе.