Прикладные задачи
^^^^^^^^^^»нелинейной теории колебаний и вслн
УДК 537.86
ОПРЕДЕЛЕНИЕ ВРЕМЕНИ ЗАДЕРЖКИ ПО ВРЕМЕННЫМ РЯДАМ НА ОСНОВЕ МЕТОДА БЛИЖАЙШИХ СОСЕДЕЙ
М.Д. Прохоров, В.И. Пономаренко} В.С. Хорев2
1 Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
2 Саратовский государственный университет им. Н.Г. Чернышевского
Предложен метод определения времени запаздывания систем с задержкой по их временным рядам, основанный на применении метода ближайших соседей. Метод может быть применен к широкому классу систем с запаздыванием и остается эффективным при высоких уровнях динамического и измерительного шума.
Ключевые слова: Системы с запаздыванием, анализ временных рядов, оценка параметров.
Введение
Автоколебательные системы с запаздывающей обратной связью чрезвычайно широко распространены в реальном мире. Например, в пространственно-развитых системах запаздывание определяется тем, что сигналы распространяются с конечной скоростью и им требуется время на преодоление расстояний [1]. Запаздывание реакции на сигнал и обратная связь с запаздыванием присущи многим физическим [2, 3], химическим [4], климатическим [5] и биологическим [6-8] объектам и процессам. При исследовании систем с запаздыванием важно знать значения временных задержек, величина которых во многом определяет динамику и свойства системы. Знание времен запаздывания имеет большое значение для построения модели системы и предсказания поведения системы во времени и при изменении параметров. Поэтому задача восстановления времени запаздывания по временным рядам наблюдаемых динамических переменных привлекает внимание многих исследователей.
Для ее решения были предложены различные методы, позволяющие восстановить времена запаздывания систем по хаотическим временным рядам их колебаний. Многие из этих методов основаны на проецировании бесконечномерного фазового пространства системы с запаздыванием в подпространства малой размерности [9-14]. При этом используются такие критерии качества реконструкции системы с запаздыванием, как минимальная ошибка прогноза построенной модели [9-11], минимальная величина информационной энтропии [12] или различные меры сложности спроецированного временного ряда [13, 14]. Известны также методы восста-
новления времени запаздывания, основанные на применении регрессионного анализа [15, 16], статистического анализа экстремумов временного ряда [17], информационных подходов [18, 19], метода множественной стрельбы [20], метода оптимизации [21] и адаптивной синхронизации [22, 23]. В отдельную группу можно выделить методы определения времени запаздывания, основанные на возмущении системы внешним воздействием и анализе отклика [24-27], которые могут быть применены к системам не только в хаотических, но и в периодических режимах.
В данной работе мы предлагаем новый метод восстановления времени задержки по временным рядам, основанный на использовании метода ближайших соседей. Метод ближайших соседей широко используется в различных научных дисциплинах при анализе временных рядов [28, 29]. При этом основными областями его применения являются классификация объектов и прогноз динамики исследуемой системы. В задаче классификации объектов идея метода ближайших соседей заключается в том, что объект полагается относящимся к тому классу, который является наиболее распространенным среди его ближайших соседей. Применительно к задаче прогноза временного ряда, суть метода ближайших соседей состоит в том, чтобы использовать для предсказания динамики системы ее состояния в прошлом, наиболее близкие по своим свойствам состоянию в настоящий момент времени. Мы впервые предлагаем использовать метод ближайших соседей для определения по временному ряду системы с запаздыванием ее времени задержки.
1. Восстановление времени задержки в системах первого порядка с запаздыванием
Поясним идею предлагаемого метода на примере одного из наиболее широко используемых дифференциальных уравнений первого порядка с одним временем запаздывания
ех(Ь) = -х{г) + / {х(Ь - т)) , (1)
где е - параметр инерционности, г - время запаздывания, / - нелинейная функция. К уравнению вида (1) сводятся, например, уравнение Маккея-Гласса [6] и уравнение Икеды [1], ставшие эталонными при исследовании систем с запаздыванием.
Поскольку при анализе временных рядов мы всегда имеем дело с переменными, измеренными в дискретные моменты времени, удобно перейти от дифференциального уравнения (1) к разностному уравнению
х(г + Дг) — х(г) .. ,, . .. е -Д-— = —х(Ь) + / (х(г — т)), (2)
где Дг - малое время, равное интервалу выборки. Уравнение (2) можно переписать в виде
х(г + Дг) = ахх(г) + 02/ (х(г — т)), (3)
где 01 = 1—Дг/е, а2 = Дг/е. Запишем уравнение (3) в виде дискретного отображения
Хп+1 = 01Хп + 02/(хп-Л), (4)
где п = г/Дг - дискретное время, а й = т/Дг - дискретное время задержки.
Пусть у нас имеется временной ряд {хп}^=1 системы (1) длиной N точек. Введем вектор Х^ = (х^х^) и найдем для него ближайшего соседа - вектор
X] = [Х],х]-а), где ] = г. Существует несколько критериев выбора ближайшего соседа для заданного вектора, отличающихся выбором метрики [29]. Наиболее часто используемой является Евклидова метрика
Ь (X,, X]) = \](хг — Х])2 + (х—а — Х]-а)2. (5)
Ближайшим соседом для вектора Хг будет такой вектор X], для которого расстояние Ь (Хг, X]) минимально. В общем случае, по временному ряду находят не одного, а к ближайших соседей для заданного вектора.
Идея предлагаемого метода состоит в том, что в системе (4) ближайшим соседним векторам, содержащим динамическую переменную в моменты времени п и п — й, где п € [й + 1, N — 1], будут соответствовать близкие состояния системы в моменты времени п + 1, так как эволюция системы (4) определяется ее текущим состоянием и состоянием в задержанный момент времени. Поскольку время запаздывания нам неизвестно, будем перебирать пробные задержки т из некоторого интервала и для к ближайших соседних векторов каждого вектора Хп = (хп, хп-т) временного ряда оценивать дисперсию оп соответствующих им состояний системы в моменты времени п + 1.
При ошибочном выборе величины т (т = й) дисперсия этих состояний может оказаться большой, так как состояния системы в моменты времени п + 1 не будут зависеть от ее состояний в моменты времени п — т. Истинному времени запаздывания й будет соответствовать положение минимума величины
1 М-1
О(т) =-1- V о2п. (6)
N — т — 2 ^ п v '
п=т+1
Мы применили метод к временным рядам уравнения Маккея-Гласса
*<« = —х('> + т+Х——Т) ■ <7>
Делением на Ь уравнение (7) приводится к виду (1) с е = 1/Ь. Параметры системы (7) а = 0.2, Ь = 0.1, с = 10, т = 300 соответствовали движению на хаотическом аттракторе, шаг выборки точек АЬ = 1, длина ряда N = 10000. Фрагмент временного ряда приведен на рис. 1, а.
На рис. 1, б приведены зависимости величины О от пробного времени запаздывания т при различном количестве к ближайших соседей для вектора Хп = (х
n, хп-т). Величина т перебирается из интервала от 1 до 500 с шагом 1.
Все зависимости О(т) демонстрируют ярко выраженный абсолютный минимум при т = 300, обеспечивая точное восстановление дискретного времени запаздывания й = т/АЬ = 300.
Отметим, что при выбранных значениях параметров процессы в системе Маккея-Гласса происходят с достаточно большой скоростью, и для численного решения системы часто используется шаг интегрирования Н = 0.1. Если при этом для построения зависимости О(т) использовать все точки временного ряда, то есть положить шаг выборки точек АЬ = 0.1, то дискретное время запаздывания будет й = 3000. В этом случае величину т потребуется перебирать в значительно более широком интервале, а длину ряда - увеличить на порядок, чтобы иметь достаточное количество колебаний системы. В результате время численного счета, необходимое
Рис. 1. а - временная реализация уравнения Маккея-Гласса в хаотическом режиме; б - зависимости величины Б от пробного времени запаздывания т при различном количестве к ближайших соседей; в - зависимости Б(т) при различном количестве I исключаемых из рассмотрения близких по времени векторов
для построения зависимости D(m), очень сильно увеличится. Для уменьшения времени счета мы пробовали при h = 0.1 брать лишь каждую десятую точку временного ряда (At = 1 и d = 300) для построения D(m). Полученные при этом для разных к графики D(m) были очень близки соответствующим графикам, представленным на рис. 1, б, и имели такой же ярко выраженный абсолютный минимум при m = 300. Поэтому в последующих примерах мы использовали h = 1 и At = 1.
При слишком частой выборке точек возможна ситуация, когда ближайшими соседями вектора Xj = (xi,xi-d) окажутся близкие ему по времени векторы Xj = (Xj,Xj-d) с j = i±p, где p = 1, 2,... P, которые мешают оценке времени запаздывания. В таких случаях при поиске ближайших соседей вектора Xj = (xi,xj-d) следует исключить из рассмотрения близкие ему по времени векторы Xj = (xj ,Xj-d) в количестве l = 2P.
На рис. 1, в построены зависимости D(m) при различном количестве l близких по времени векторов, не учитываемых при поиске ближайших соседей, для случая к = 10. Как и на рис. 1, б, все графики имеют абсолютный минимум при m = d = 300.
Следует отметить, что вместо поиска фиксированного числа к ближайших соседей вектора Xj = (xj, x—d) можно считать его ближайшими соседями все векторы Xj = (xj, xj-d), для которых L (Xj, Xj) < 8, где 8 - некоторая малая величина. Графики D(m), построенные при таком способе нахождения ближайших соседних векторов, очень похожи на графики, представленные на рис. 1, б. Подобрав параметры к и 8 соответствующим образом, можно добиться практически полного совпа-
дения результатов поиска ближайших соседей обоими способами. Кроме того, нами установлено, что выбор метрики для поиска ближайших векторов почти не влияет на вид зависимостей D(m).
Для оценки возможностей метода при наличии возмущений мы применили его к рядам, полученным при добавлении к временному ряду уравнения (7) гауссова белого шума с нулевым средним значением. Полученные результаты представлены на рис. 2, а для различных уровней измерительного шума при k = 10 и l = 10. Положение минимума D(m) позволяет точно восстановить время запаздывания даже при уровнях шума порядка 30% (отношение сигнал/шум около 10 дБ), что существенно превышает уровень шума, допустимый для применения большинства других методов восстановления времени запаздывания.
По отношению к динамическому шуму предложенный нами метод оказывается еще более устойчивым. На рис. 2, б зависимости D(m) при k = 10 и l = 10 построены для случая, когда гауссов белый шум с нулевым средним значением добавлен в правую часть уравнения (7). Минимум всех графиков D(m), построенных на рис. 2, б при различных уровнях динамического шума, наблюдается при m = 300.
Рассмотрим случай, когда система (7) совершает периодические колебания (a = 0.2, b = 0.1, c = 10, т = 10). Фрагмент временного ряда таких колебаний приведен на рис. 3, а. При построении зависимостей D(m) будем считать ближайшими соседями векторы, для которых L (Xj, Xj) < 0.02. В отсутствие шума на графике D(m) нет выраженного минимума (рис. 3, б). Однако добавление динамического шума оказывается полезным для восстановления времени запаздывания. При 10% шуме (отношение сигнал/шум равно 20 дБ) на графике D(m) появляется отчетливый минимум при m = d = 10 (см. рис. 3, б). Присутствие измерительного шума, разумеется, не оказывает положительного влияния на определение времени запаздывания.
В отличие от большинства других методов восстановления времени запаздывания, предложенный нами подход можно применять для определения времени задержки не только систем вида (1) с линейной зависимостью от текущего состояния и нелинейной зависимостью от состояния в задержанный момент времени, но и систем более общего вида, с нелинейной функцией F, зависящей от обеих переменных x(t) и x(t — т)
X(t) = F (x(t),x(t — т)). (8)
Рассуждения, приведенные выше для системы (1), остаются справедливыми
D 0.12 0.08 0.04 0
а
о
D
30% 0.016
20% 0.012
\г 10% 0.008
0.004
у V ' 0%
^ . I
. I
100 200 300 400
300 400
Рис. 2. Зависимости В(гп) для системы Маккея-Гласса в хаотическом режиме при различных уровнях измерительного шума (а) и динамического шума (б). Уровни шума указаны в процентах на соответствующих кривых
а о
Рис. 3. а — временная реализация уравнения Маккея-Гласса в периодическом режиме; б — зависимости Б(т) при различных уровнях динамического шума, указанного в процентах на соответствующих кривых
для систем вида (8), поскольку в основе метода лежит учет зависимости эволюции исследуемой системы от ее текущего состояния и состояния в задержанный момент времени, а сам вид зависимости не имеет значения.
2. Восстановление времени задержки в системах второго порядка с запаздыванием
Предложенный метод может быть легко развит на системы высокого порядка с задержкой, в частности на системы, описываемые дифференциальными уравнениями с запаздыванием второго порядка
£2Х(Ь) + £1Х(Ь) = ^ (х(Ь),х(Ь - т)), (9)
где £\ и £2 - параметры инерционности. Проиллюстрируем это на примере системы следующего вида:
£2Х{г) + £1 х(Ь) = -х(г) + / {х(Ь - т)) . (10)
Используя описанный выше формализм, можно перейти от дифференциального уравнения (10) к дискретному отображению
Хп+2 = Ь1 Хп+1 + Ь2 Хп + Ьз/(хп-а), (11)
где Ь1 =2 - (^Д*) /£2, Ь2 = -1 + (^Д* - (Д*)2) I£2, Ьз = (М)2 !£2.
Для каждого вектора Хп = (хп+1, хп, хп-т) временного ряда системы (10) найдем к ближайших соседних векторов и оценим дисперсию соответствующих им состояний системы в моменты времени п + 2. Перебирая различные пробные задержки т, будем искать минимум величины
1 N-2
°(т) = м-т-3 £ "«• (12)
п=т+1
положение которого даст оценку дискретного времени запаздывания с! = т/Д*.
Описанные методы могут быть использованы для определения по временному ряду априорно не известного порядка системы с запаздыванием. Идея подхода состоит в том, чтобы сначала восстановить время запаздывания исследуемой системы
в предположении, что она описывается уравнением с запаздыванием первого порядка (1), а затем в предположении, что модельным уравнением системы является уравнение второго порядка (10), и построить зависимости (6) и (12) на одном графике. При правильном выборе порядка модельного дифференциального уравнения зависимость D(m) будет лежать ниже зависимости D(m), построенной при ошибочном выборе порядка модельного уравнения.
Например, пусть у нас имеется временной ряд системы с запаздыванием второго порядка (10) с квадратичной нелинейной функцией f (x) = X — x2, где X - параметр нелинейности, находящейся под действием 10% динамического шума. Параметры системы т = 1000, X = 1.9, £\ = 7, £2 = 10 соответствуют движению на хаотическом аттракторе. Шаг выборки точек At = 1, длина ряда N = 10000. Фрагмент временного ряда приведен на рис. 4, а. Будем считать, что порядок модельного уравнения этой системы нам не известен, и восстановим сначала время запаздывания, полагая, что система описывается уравнением с запаздыванием первого порядка (1). На рис. 4, б построена зависимость (6) при k = 10 и l = 10. Она имеет минимум при m = 1001, что немного больше времени запаздывания d = т/At = 1000.
Восстановим теперь время запаздывания, предположив, что исследуемая система описывается уравнением с запаздыванием второго порядка (10). На рис. 4, б построена зависимость (12) при k = 10 и l = 10. Она располагается ниже зависимости (6), что указывает на то, что система лучше описывается уравнением второго порядка. Минимум зависимости (12) наблюдается при m = d = 1000, то есть при правильном выборе порядка модельного уравнения время запаздывания удается восстановить точно.
Рассмотрим теперь ситуацию, когда исследуемый временной ряд сгенерирован системой с запаздыванием первого порядка (1) с квадратичной нелинейной функцией и параметрами т = 1000, X = 1.9, £ = 10, соответствующими движению на хаотическом аттракторе. Как и в рассмотренном выше примере, на систему действует 10% динамический шум, At = 1, N = 10000.
При реконструкции времени задержки в предположении, что модельное уравнение имеет вид (1), график
D
-2'
1.0
-1.0
-2.0
а
D
0.008
0.004
2\
1 , 1 |
300 600 900 1200 t
б
400
800
1200
m
Рис. 4. а — временная реализация уравнения (10) с квадратичной нелинейностью в хаотическом режиме; б — зависимости В(гп), построенные в предположении, что модельное уравнение имеет первый порядок (1) и второй порядок (2)
0
О(т), построенный при к = 10 и I = 10, демонстрирует минимум при т = 6 = 1000. (Кривая 1 на рис. 5.)
Зависимость О(т), построенная в предположении, что модельное уравнение имеет вид (10), показана на рис. 5 (кривая 2). Она имеет минимум при т = 999 и расположена преимущественно выше кривой 1, что указывает на то, что модельное уравнение системы имеет первый порядок.
Отметим, что при выборе шага интегрирования для рассмотренных в данном разделе систем мы руководствовались такими же соображениями, как и в разделе 1 для системы Маккея-Гласса. При Н = 1 и Д* = 1 количество точек на характерном временном масштабе колебаний оказывается достаточным для успешной реализации предложенного метода и построения зависимостей О(т) .
Предложенный метод может быть распространен и на системы с несколькими временами запаздывания. Рассмотрим систему с запаздыванием с двумя различными временами задержки т1 и т2
Используя описанный выше подход, можно перейти от дифференциального уравнения (13) к дискретному отображению
где а1 = 1 - Д*/£, а2 = Д*/£, 61 = т1/Д*, 62 = т2/Д*.
Из уравнения (14) следует, что ближайшим соседним векторам, содержащим динамическую переменную в моменты времени п, п - 6,1 и п - 62, где п € [62 + 1, N - 1] (62 > 61), будут соответствовать близкие состояния системы в моменты времени п + 1. Поскольку времена запаздывания 61 и 62 нам не известны, будем перебирать пробные задержки т1 и т2 из некоторого интервала и для к ближайших соседних векторов каждого вектора Хп = (хп, хп-т1, хп-т2) временного ряда оценивать дисперсию оп соответствующих им состояний системы в моменты времени п + 1 .
При ошибочном выборе т1 и/или т2 (т1 = 61, т2 = 62) дисперсия этих состояний может оказаться большой. В качестве оценки времен запаздывания 61 и 62 будем использовать положение минимума величины
3. Восстановление времен запаздывания в системах с двумя задержками
£x(t) = -X(t) + fi (x(t - Ti)) + /2 (x(t - T2)) .
(13)
Xn+i = aiXn + a2/i (xn-di) + a2/2(xn-d2),
(14)
1
N-1
D(mi,m2)
N - m2 - 2
(15)
n=m,2+i
Продемонстрируем работоспособность метода, применив его к временному ряду обобщенного уравнения Маккея-Гласса, получаемого введением второго времени запаздывания
, . . 1 ax(t - t1)
x(t) = -bx(t) +--Ц————
w w 2 1+ xc(t - T—)
1 ax(t - t2) + 2 1 + xc(t - T2) '
(16)
Поделив обе части уравнения (16) на Ъ, можно привести его к виду (13) с £ = 1/Ъ. При а = 0.2, Ъ = 0.1, с = 10, т1 = 70, т2 = 300 система (16) демонстрирует хаотические колебания. По временному ряду длиной N = 10000 точек при шаге выборки Дí = 1 построим зависимость 0(т1, т2), представляющую собой двумерную поверхность (рис. 6, а). При ее построении мы считали ближайшими соседями вектора Х^ и Xj, расстояние
L (Xii Xj) — \J\xi xj) + (xi-mi xj-mi) + (xi-m2 xj-m2)
(17)
между которыми Ь (XXj) < 0.02. Зависимость 0(т1,т2) имеет минимум при т1 = 61 = 70, т2 = 62 = 300, обеспечивая точное восстановление обоих времен запаздывания. Так же как при численном интегрировании системы (7), мы выбрали Н = 1. Такой шаг интегрирования позволяет при Дt = 1 точно восстановить времена запаздывания, а уменьшение шага приводит к существенному увеличению вычислительных затрат для построения зависимости 0(т1, т2).
Для сравнения на рис. 6, б приведена зависимость О(т) , полученная при применении к временному ряду системы (16) метода, предложенного в разделе 1 для систем (1) с одной задержкой. График этой зависимости (6) имеет глубокие минимумы при т = 68 и т = 298. Таким образом, без учета вида модельного уравнения (13) оценка времен запаздывания оказывается менее точной. Еще один характерный минимум Б(т) наблюдается на рис. 6, б вблизи т = 61 + 62.
Отметим, что, если применить метод восстановления двух задержек к временному ряду системы (1) с одним временем запаздывания, минимум величины 0(т1, т2) будет наблюдаться при т1 = т2 = 6.
Описанный в данном разделе метод можно применять для определения времен запаздывания не только систем вида (13), но и систем более общего вида с двумя задержками
x(t) = F (x(t),x(t - T—),x(t - T2)).
(18)
0.0002-
0.0001-
а
/ У \1Ь Л
У
290 60
100 200 300 400
Рис. 6. Зависимости П(гп1,гп2) и В(гп) для обобщенной системы Маккея-Гласса в хаотическом режиме
4. Восстановление времени задержки
в системах с запаздыванием с двумя динамическими переменными
Предложенный метод восстановления времени запаздывания по временным рядам, основанный на применении метода ближайших соседей, может быть распространен на системы, описываемые дифференциальными уравнениями с запаздыванием с несколькими динамическими переменными:
х(Ь) = ^ (х(Ь),х(Ь - т),у(Ь - т)),
(19)
у(Ь) = ^ (у(Ь),у(Ь - т),х(Ь - т)).
В этом случае, используя временные ряды обеих переменных х(Ь) и у (г), будем искать ближайших соседей для векторов Хп = (хп,хп-т,уп-т) или V« = (у п, уп- т, хп- т ) при различных пробных задержках т и определять время запаздывания по положению минимума величины (6).
Применим метод к временным рядам системы двух связанных нелинейных дифференциальных уравнений с запаздыванием
х(Ь) = тх(Ь) - ц [х2(Ь - т) + су2(Ь - т) х(Ь),
(20)
у (г) = ту (г) - ц [у2 (г - т) + сх2(г - т) у (г),
введенной в [30]. Выберем следующие значения параметров: т = 4, ц = 4, с = 0.5, т = 0.35. Как показано в [30], при этих параметрах система (20) демонстрирует периодические колебания. Фрагмент временного ряда переменной х(Ь) приведен на рис. 7, а.
По временным рядам переменных х(Ь) и у(Ь) длиной N = 10000 точек при шаге выборки АЬ = 0.01 построим зависимость О(ш) . В отсутствие шума на графике О(ш) нет выраженного минимума. Однако при добавлении в правую часть обоих уравнений системы (20) достаточно больших независимых динамических шумов на графике О(ш) появляется минимум. На рис. 7, б зависимость О(ш) построена при к = 10,1 = 10 и 40% динамических шумах (отношение сигнал/шум около 8 дБ). Она имеет минимум при т = 35, что совпадает с временем запаздывания с! = т/АЬ = 35.
0 5.0 10.0 t х 15 30 45 60 т
Рис. 7. а — временная реализация уравнения (20); б — зависимость Б(т) при 40% динамических шумах
Заключение
Нами предложен метод восстановления времени запаздывания систем с задержкой по их временным рядам, основанный на применении метода ближайших соседей. Метод позволяет определить время запаздывания систем высокого порядка с задержкой, систем с несколькими временами запаздывания и систем, описываемых дифференциальными уравнениями с запаздыванием с несколькими динамическими переменными. Метод может быть применен к системам с запаздыванием с произвольным видом нелинейной функции, в том числе зависящей одновременно от задержанной и незадержанной переменной. Кроме того, метод можно использовать для определения по временному ряду априорно не известного порядка системы с запаздыванием.
Параметры метода могут быть выбраны в широких пределах, а сам метод остается эффективным при очень высоких уровнях как динамического, так и измерительного шума. Показано, что метод позволяет восстановить время запаздывания не только систем с задержкой, совершающих хаотические колебания, но и систем с задержкой в периодических режимах, находящихся под действием динамического шума.
Работа выполнена при государственной поддержке ведущих научных школ, грант НШ-1726.2014.2 и при поддержке РФФИ, грант № 13-02-00227.
Библиографический список
1. Ikeda K. Multiple-valued stationary state and its instability of the transmitted light by a ring cavity system // Opt. Commun. 1979. Vol. 30. P. 257.
2. Lang R., Kobayashi K. External optical feedback effects on semiconductor injection lasers // IEEE J. Quantum Electron. 1980. Vol. 16. P. 347.
3. Erneux T. Applied Delay Differential Equations. New York: Springer, 2009.
4. Epstein I. R. Delay effects and differential delay equations in chemical-kinetics // Int. Rev. in Phys. Chem. 1992. Vol. 11. P. 135.
5. Mokhov I.I., Smirnov D.A. El Nino Southern Oscillation drives North Atlantic Oscillation as revealed with nonlinear techniques from climatic indices // Geophys. Research Lett. 2006. Vol. 33. L03708.
6. Mackey M.C., Glass L. Oscillations and chaos in physiological control systems // Science. 1977. Vol. 197. P. 287.
7. Kuang Y Delay Differential Equations with Applications in Population Dynamics. Boston: Academic Press, 1993.
8. Bocharov G.A., Rihan F.A. Numerical modelling in biosciences using delay differential equations // J. Comp. Appl. Math. 2000. Vol. 125. P. 183.
9. Fowler A.C., Kember G. Delay recognition in chaotic time series // Phys. Lett. A. 1993. Vol. 175. P. 402.
10. Hegger R., Bunner M.J., Kantz H., Giaquinta A. Identifying and modeling delay feedback systems // Phys. Rev. Lett. 1998. Vol. 81. P. 558.
11. Bunner M.J., Ciofini M., Giaquinta A., Hegger R., Kantz H., Meucci R., Politi A. Reconstruction of systems with delayed feedback: (I) Theory // Eur. Phys. J. D. 2000. Vol. 10. P. 165.
12. Tian Y.-C, Gao F. Extraction of delay information from chaotic time series based on information entropy // Physica D. 1997. Vol. 108. P. 113.
13. Kaplan D.T., Glass L. Coarse-grained embeddings of time series: Random walks, gaussian random process, and deterministic chaos // Physica D. 1993. Vol. 64. P. 431.
14. BUnner M.J., Meyer Th., Kittel A., Parisi J. Recovery of the time-evolution equation of time-delay systems from time series // Phys. Rev. E. 1997. Vol. 56. P. 5083.
15. Voss H., Kurths /.Reconstruction of non-linear time delay models from data by the use of optimal transformations // Phys. Lett. A. 1997. Vol. 234. P. 336.
16. Ellner S.P., Kendall B.E., Wood S.N., McCauley E., Briggs C.J.Inferring mechanism from time-series data: Delay differential equations // Physica D. 1997. Vol. 110. P. 182.
17. Пономаренко В.И., Прохоров М.Д., Караваев А.С., Безручко Б.П. Определение параметров систем с запаздывающей обратной связью по хаотическим временным реализациям // ЖЭТФ. 2005. Т. 127. Вып. 3. С. 515.
18. Udaltsov VS., Larger L., Goedgebuer J.P., Locquet A., Citrin D.S. Time delay identification in chaotic cryptosystems ruled by delay- differential equations // J. Opt. Technology. 2005. Vol. 72. P. 373.
19. Zunino L., Soriano M.C., Fischer I., Rosso O.A., Mirasso C.R. Permutation-information-theory approach to unveil delay dynamics from time- series analysis // Phys. Rev. E. 2010. Vol. 82. 046212.
20. Horbelt W., Timmer J., Voss H.U. Parameter estimation in nonlinear delayed feedback systems from noisy data // Phys. Lett. A. 2002. Vol. 299. P. 513.
21. Dai C., Chen W., Li L., Zhu Y., Yang Y Seeker optimization algorithm for parameter estimation of time-delay chaotic systems // Phys. Rev. E. 2011. Vol. 83. 036203.
22. Sorrentino F. Identification of delays and discontinuity points of unknown systems by using synchronization of chaos // Phys. Rev. E. 2010. Vol. 81. 066218.
23. Ma H., Xu B., Lin W., Feng /.Adaptive identification of time delays in nonlinear dynamical models // Phys. Rev. E. 2010. Vol. 82. 066210.
24. SiefertM. Practical criterion for delay estimation using random perturbations // Phys. Rev. E. 2007. Vol. 76. 026215.
25. Ponomarenko V.I., Prokhorov M.D. Recovery of systems with a linear filter and nonlinear delay feedback in periodic regimes // Phys. Rev. E. 2008. Vol. 78. 066207.
26. Yu D., Frasca M., Liu F. Control-based method to identify underlying delays of a nonlinear dynamical system // Phys. Rev. E. 2008. Vol. 78. 046209.
27. Prokhorov M.D., Ponomarenko V.I. Reconstruction of time-delay systems using small impulsive disturbances // Phys. Rev. E. 2009. Vol. 80. 066206.
28. Farmer J.D., Sidorowich J.J. Predicting chaotic time series // Phys. Rev. Lett. 1987. Vol. 59. P. 845.
29. Garcm P., /шйтг J., Marcano A., Moleiro F. Local optimal metrics and nonlinear modeling of chaotic time series // Phys. Rev. Lett. 1996. Vol. 76. P. 1449.
30. Villermaux E. Memory-induced low frequency oscillations in closed convection boxes // Phys. Rev. Lett. 1995. Vol. 75. P. 4618.
Поступила в редакцию 25.06.2013 После доработки 6.12.2013
DELAY TIME ESTIMATION FROM TIME SERIES BASED ON NEAREST NEIGHBOR METHOD
M.D. Prokhorov1 V.I. Ponomarenko} VS. Khorev2
1 Institute of Radio-Engineering and Electronics RAS, Saratov Branch 2 Saratov State University
The method is proposed for delay time estimation in time-delay systems from their time series. The method is based on the nearest neighbor method. It can be applied to a wide class of time-delay systems and it is still efficient under very high levels of dynamical and measurement noise.
Keywords: Time-delay systems, time series analysis, parameter estimation
Прохоров Михаил Дмитриевич - родился в Саратове (1968). Окончил Саратовский государственный университет (1992). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1997) и доктора физико-математических наук (2008). Ведущий научный сотрудник Саратовского филиала Института радиотехники и электроники РАН. Область научных интересов - нелинейная динамика и ее приложения, математическое моделирование, анализ временных рядов. Имеет более 150 научных публикаций.
410019 Саратов, ул. Зеленая, 38
Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН E-mail: mdprokhorov@yandex.ru
■V» |
Пономаренко Владимир Иванович - родился в Саратове (1960). Окончил Саратовский государственный университет (1982). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (1992) и доктора физико-математических наук (2008). Ведущий научный сотрудник Саратовского филиала Института радиотехники и электроники РАН. Область научных интересов - статистическая радиофизика, анализ временных рядов, нелинейная динамика и ее приложения. Автор более 150 научных публикаций.
410019 Саратов, ул. Зеленая, 38
Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН E-mail: ponomarenkovi@gmail.com
Хорев Владимир Сергеевич - родился в Гомеле (1987). Окончил магистратуру в Саратовском государственном университете (2011). Аспирант кафедры динамического моделирования и биомедицинской инженерии СГУ. Область научных интересов - моделирование по временным рядам, нелинейная динамика и ее приложения к задачам биофизики и медицины. Принимал участие в 10 международных и всероссийских конференциях. Имеет 14 научных публикаций.
410012 Саратов, ул. Астраханская, 83
Саратовский государственный университет им. Н.Г. Чернышевского E-mail: zergs@inbox.ru