Знакомство с нелинейной динамикой: Лекции соросовского профессора - Анищенко В.С.
Скачать (прямая ссылка):
Лекция 10. Реконструкция динамических систем
Введение
Важнейший метод исследования эволюционных процессов в естествознании состоит в построении математических моделей изучаемых систем и их анализе. Как сказал один из великих мыслителей прошлого века, любое утверждение истинно настолько, насколько оно базируется на математике. Наличие математической модели исследуемой системы существенно расширяет возможности ее изучения, позволяя решать задачи предсказания поведения системы во времени и зависимости режимов ее функционирования от параметров. Таким образом, одной из центральных является задача математического моделирования, решение которой дает возможность осуществления научного прогноза функционирования системы во времени, являющегося одной из главных проблем в естествознании.
Решение задачи моделирования теоретически не содержит проблем, если реальная система задана. Хорошо известный пример — колебательный LC-контур. На основе знания схемы контура и электрических законов не представляет труда записать основополагающие соотношения и получить уравнения классического осциллятора:
x + cj$x = Q,uj$ = 1/LC. (10.1)
Решением уравнения (9.1) является гармонический колебательный во времени процесс, частота которого определяется параметрами контура L и С. При заданных начальных условиях х(to) и x(to) состояние системы (9.1) будет однозначно известно для любого времени t ^ to-
Однако очень часто приходится сталкиваться с более сложной ситуацией, когда детальные сведения о реальной системе либо отсутствуют вовсе, либо явно недостаточны. Единственная информация о свойствах системы содержится лишь в экспериментальной зависимости одной из координат состояния системы во времени. Такая зависимость a(t), измеренная в течение конечного времени to, называется наблюдаемой (или реализацией) системы, а при дискретизации с шагом At: а (г At) = сц, і = 1,..., TV; N = [to/At], она носит название одномерного временного ряда. Делается предположение о том, что наблюдаемая a(t) является детерминировано определенной, то есть представляет собой одномерную проекцию фазовой траектории, порождаемой некоторой динамической системой (ДС). Задачей реконструкции динамической системы является восстановление модельной ДС, решение которой с известной степенью точности воспроизводит одномерную наблюдаемую a(t) на заданном интервале времени to и для t > to - Проблема реконструкции ДС, Введение
131
таким образом, относится к классу обратных задач, решение которых не может быть однозначным.
В рамках настоящей лекции мы ограничимся рассмотрением проблемы восстановления (реконструкции) ДС в виде системы обыкновенных дифференциальных уравнений (ОДУ) или дискретных отображений по одномерному временному ряду. Важно при этом не забывать, что временной ряд а (і At) предполагается детерминировано определенным, то есть отражает эволюционный процесс реальной ДС, управляемой детерминированными динамическими законами. Естественно, что такое предположение сужает класс рассматриваемых сигналов. В частности, если временной ряд есть следствие абсолютно случайного (шумового) процесса, то говорить о реконструкции не имеет смысла [4,25].
Попытаемся понять основные проблемы, с которыми связано решение задачи реконструкции. Первая обусловлена необходимостью введения каким-либо образом координат состояния системы. Ведь нам известна зависимость во времени (на конечном интервале!) лишь одной из координат реальной системы a(t). Как ввести новые координаты и сколько их должно быть? Какова размерность ДС, которую мы хотим восстановить?
Предположим, что нам удалось как-то решить эту проблему и мы знаем размерность п модельной ДС. Сразу возникает второй, не менее важный вопрос: как записать сами уравнения? Каков вид модельного оператора эволюции, который в случае ОДУ определяется правыми частями системы п дифференциальных уравнений первого порядка? Дать обоснованные ответы на эти два вопроса по сути дела и есть содержание раздела теории динамических систем, рассматривающего проблему реконструкции ДС по одномерным временным рядам.
В 1980 г. была опубликована работа Н. Пакарда, в которой показано, что фазовый портрет динамической системы может быть восстановлен по скалярному временному ряду сц, если в качестве недостающих координат вектора состояния используется тот же самый ряд сц, взятый с некоторым запаздыванием. В 1981 г. была доказана теорема, утверждающая, что по одномерной реализации a(t) ДС, обладающей аттрактором A7 принадлежащим гладкому d-мерному многообразию, методом задержки можно получить n-мерную реконструкцию Ar исходного аттрактора как множество векторов x(t) в Rn при п ^ 2d + 1 (теорема Такенса):
x(t) = An(a(t)) = (a(t), a(t + г),..., a(t + (тг - 1 )т)) =
= (жі,ж2, • • -,?). (10.2) 132
Лекция 10. Реконструкция динамических систем
Согласно теореме, отображение An : А —Ar является гладким и обратимым на Ar почти при любой задержке т (если N —ос).
Попытаемся разобраться в содержании теоремы Такенса. Она обосновывает введение в качестве новых координат состояния системы значений одномерного временного ряда а (г At), взятых через некоторый интервал времени т:
