Введение

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

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

Одно из центральных мест в области моделирования процессов на границе водной и воздушной сред занимает двумерное обобщение уравнения Курамото–Цузуки (Гинзбурга–Ландау) [1,2]. При определённых значениях параметров в системе наблюдается так называемая жёсткая турбулентность – хаотический режим с редкими сильными выбросами (пиками), хорошо подходящий для формализованного описания некоторых катастрофических событий.

Для исследования подобных выбросов используется модификация двумерного уравнения Гинзбурга–Ландау – одномерное зависящее от времени уравнение Гинзбурга–Ландау с нелинейностью 5 порядка. Уравнение содержит малый параметр ε, характеризующий интенсивность диссипации. В пределе при ε → 0 данное уравнение переходит в нелинейное уравнение Шрёдингера. Известные на сегодняшний день уравнения, описывающие рост и спад выброса, получены путём исследования последнего [3].

Остаётся открытым вопрос описания фоновой (межпиковой) динамики и перехода от фоновой динамики к пиковой. Для этого было предложено использовать более простые системы, например, системы, демонстрирующие так называемую переключающуюся перемежаемость – хаотический режим с наличием устойчивого инвариантного многообразия, со временем теряющего устойчивость. Для получения эффекта жёсткой турбулентности подобную систему достаточно дополнить механизмом возвращения устойчивости инвариантному многообразию. Примером такой системы является отображение Ершова [4].

Далее обсуждается вопрос моделирования жёсткой турбулентности в условиях неполноты информации о наблюдаемой системе. В качестве последней выбрана система Ершова, причём, по условиям задачи, для наблюдения доступна только одна из трёх переменных.

Для воссоздания наблюдаемого поведения предлагается использовать метод русел и джокеров, основанный на конструировании модели из «разнородных» компонент – детерминированных и стохастических, непрерывных и дискретных, и т.д. [3,4,5,6,7]. Использование метода оправданно, когда фазовое пространство динамической системы неоднородно, и выделяются области сравнительно простого поведения (русла). В пределах русла система может быть описана простой динамической моделью. Там, где это не представляется возможным, для описания сложного поведения используется джокер – некоторый простейший алгоритм, как правило, вероятностный.

В системе Ершова рост и спад выброса представляют собой пример простого поведения, и им сопоставляются русла. В области устойчивого многообразия система демонстрирует сложное (хаотическое) поведение, для моделирования которого используется джокер.

Главное требование, предъявляемое к разрабатываемой модели – соответствие основных характеристик временных рядов, порождаемых ею и исходной моделью. В качестве таких характеристик были выбраны распределения величины выбросов и длительности межпиковых интервалов.

Уравнение Курамото–Цузуки (Гинзбурга–Ландау) и жёсткая турбулентность

Уравнение Курамото–Цузуки (иначе называемое уравнением Гинзбурга–Ландау) представляет собой одну из основных математических моделей, занимающих особое место в нелинейной динамике [1,2,3]. Данное уравнение демонстрирует большое разнообразие решений, включая стационарные, периодические и хаотические, и используется, в частности, для описания явлений, имеющих место в нелинейных средах – перехода к турбулентным режимам в гидродинамике, ионно-звуковых волн в плазме, и т.д. Нас интересует его двумерное обобщение, используемое при моделировании процессов, происходящих на границе сред – ветровых волн на воде, морфогенеза, протекающих на поверхности колебательных химических реакций [8].

Двумерное обобщение уравнения Курамото–Цузуки выглядит следующим образом:

. (1)

Здесь W – некоторая комплекснозначная функция действительных переменных, а c0, c1, c2 – параметры (действительные постоянные). При определённом соотношении значений параметров в системе (1) можно наблюдать жёсткую турбулентность (в англоязычных источниках используется термин hard/strongturbulence) – хаотический режим с редкими и исключительно сильными выбросами (пиками жёсткой турбулентности). Для исследования данного режима используется специальная модификация уравнения Курамото–Цузуки – одномерная, но с более сильной нелинейностью:

. (2)

Данное уравнение известно под названием зависящего от времени уравнения Гинзбурга–Ландау с нелинейностью 5 порядка (quintic time-dependent Ginzburg–Landau equation, сокращённо – QTDGL). При ε = 0 оно переходит в нелинейное уравнение Шрёдингера

. (3)

Важное свойство последнего – существование ряда законов сохранения, существенных для описания жёсткой турбулентности. Это закон сохранения «массы»

,  

«импульса»

 

и «энергии»

.  

На сегодняшний день также известны следующие свойства решений уравнения (2). Во-первых, при малых ε решение уравнения (3) представляет собой хорошую аппроксимацию профиля W(x,t) в области роста больших пиков. Во-вторых, распределение максимумов пиков h описывается (в промежуточном интервале масштабов) степенной асимптотикой p(h) ~ h–α, где α – константа, α ≈ 7÷8. В области малых h имеет место иное распределение. Также максимумы пиков ограничены сверху, что, предположительно, объясняется воздействием множителя ε в уравнении (2).

На практике, как правило, из всех характеристик распределения наибольший интерес представляют моменты первого и второго порядка. Здесь существенна интерпретация W. Если имеют смысл величины |W| и |W|2, то аномально большие пики не оказывают принципиального влияния на среднее и дисперсию, то есть, распределение p(h) можно считать эффективно ограниченным. Если же рассматривается величина |W|4, то имеет место картина, характерная для некоторых разновидностей катастрофических событий, когда редкие аномальные выбросы оказывают существенное влияние на среднее и обращают в бесконечность дисперсию.

В численных экспериментах жёсткую турбулентность можно наблюдать, решая уравнение (2) в области длины L с периодическими граничными условиями, с параметром ε в интервале 0,01÷0,0001. Поведение системы зависит от размера области: при L < 5 пики жёсткой турбулентности отсутствуют, в областях промежуточной длины (приблизительно от 10 до 80) наблюдаются возникновения одиночного пика, в более длинных областях возможно возникновение нескольких пиков одновременно.

Внешне, жизненный цикл одиночного пика жёсткой турбулентности выглядит следующим образом (здесь и далее предполагается рассмотрение в области промежуточной длины). Вначале, в системе наблюдаются медленные нерегулярные пространственно-временные колебания. В каждый момент времени график W(x) выглядит гладким, причём |W| < 1.

Затем, в какой-то момент времени начинается рост пика. Скорость роста значительно превышает характерную величину Wt вне области пика. Ширина пика при этом сокращается.

Далее, после того, как пик достигает максимальной величины (обратно зависящей от ε), начинается его распад. Вначале фазы распада величина пика уменьшается, причём уменьшение происходит быстрее предшествующего роста. Ширина пика при этом увеличивается, а его форма сглаживается.

Распад пика продолжается его превращением в быстро осциллирующий волновой пакет с колоколообразной огибающей. Высота пика при этом имеет величину порядка единиц.

После этого, волновой пакет распространяется в ширину, до заполнения всей области рассмотрения. Амплитуда волн продолжает уменьшаться, а высшие гармоники затухают.

В конце концов, в области рассмотрения вновь наблюдаются медленные колебания. Характерное время сглаживания профиля имеет порядок ε–1. Новый пик может образоваться не прежде, чем происходит сглаживание. Характерный вид профилей W в различных фазах пика жёсткой турбулентности показан, в частности, в [9].

В ходе эволюции пика масса, импульс и энергия ведут себя следующим образом. В начале выброса происходит всплеск энергии – практически мгновенный рост на несколько порядков. Затем, в ходе распада пика и далее, энергия медленно и монотонно убывает. Новые пики не возникают, пока энергия не достигнет величины порядка 1.

Масса после возникновения пика начинает убывать, пока не перестанет выполняться условие E > M. После этого начинается рост. В целом, масса меняется значительно слабее, чем энергия: минимальное и максимальное значения различаются менее чем на порядок.

Импульс практически не отражает развитие пика: в расчётах, где его начальное значение близко к нулю, его дальнейшее изменение незначительно.

В настоящее время не известны простые способы точного предсказания момента и координаты начала следующего пика, а также его высоты. В этой связи представляется естественной попытка дополнить динамическое описание системы статистическим, возможно, с привлечением более простых моделей жёсткой турбулентности, чем уравнение (2).

Переключающаяся перемежаемость и отображение Ершова

Один из вариантов упрощённого (качественного) описания жёсткой турбулентности основан на использовании систем, демонстрирующих так называемую переключающуюся перемежаемость (on-off intermittency). Приведённый термин используется для обозначения разновидности динамического режима, в котором хаотическое движение чередуется с резкими сильными выбросами. В данном случае такое поведение системы обусловлено наличием устойчивого инвариантного многообразия, теряющего устойчивость в результате бифуркации. Кроме того, система дополнена механизмом возвращения устойчивости инвариантному многообразию. Отметим следующее: предполагается, что хаотическое движение длится достаточно долго, чтобы точное определение основных характеристик выброса (момент начала, длительность, высота, направление) было невозможно.

Примером системы, в которой можно наблюдать переключающуюся перемежаемость, является отображение Ершова [3,4]. Данная система строилась с учётом следующих требований:

  1. Наличие одной медленной переменной, играющей роль параметра в отображении для быстрых переменных. Данная переменная, по аналогии с уравнением (2), называется «энергией» и обозначается E, так как предполагается, что именно энергия, среди всех интегралов движения уравнения (2), играет наиболее существенную роль.
  2. Наличие быстрых переменных, описывающих соответствующее межпиковой фазе движение на хаотическом аттракторе.
  3. Изменение свойств аттрактора в зависимости от значения медленной переменной. При переходе медленной переменной некоторого критического значения должен происходить кризис аттрактора – появление неустойчивой области, так называемой «дырки», через которую изображающая точка покидает аттрактор, что соответствует росту пика. Здесь и далее имеется в виду изображающая точка в фазовом пространстве быстрых переменных.
  4. Прекращение роста пика и возвращение изображающей точки на аттрактор (аттрактор вновь становится глобально притягивающим). Существует два очевидных способа определения момента возвращения аттрактору свойства притяжения. Первый, параметрический, предполагает зависимость от значения энергии (рост прекращается при превышении энергией некоторого порогового значения). Второй, динамический – зависимость от значений быстрых переменных (рост прекращается при попадании изображающей точки в заданную область фазового пространства). Динамический способ считается в большей степени соответствующим сути поведения системы, заданной уравнением (2). Здесь следует отметить, что для простоты было бы естественно ограничиться единственной быстрой переменной. Однако, в сочетании с динамическим способом ограничения роста пика, это привело бы к колебаниям вокруг порогового значения (при превышении порогового значения рост прекращается, начинается спад, после обратного перехода через порог вновь начинается рост, и так далее). Поэтому, необходимо наличие как минимум двух быстрых переменных, и несколько более сложный механизм остановки роста пика. Например, следующий: изображающая точка, приближается к аттрактору по первой координате и продолжает удаляться по второй, пока не произойдёт возврат по первой координате в область аттрактора. После этого, начинается возврат по второй координате.
  5. Максимально возможная простота функции отображения быстрых переменных. Целесообразно использование кусочно-линейной функции.

Система строится следующим образом. Есть две быстрые переменные, x и y, и соответствующие им два связанных одномерных отображения с параметром E. Когда E выше порогового значения Ecr, имеется хаотический аттрактор с некоторой конечной областью притяжения. При уменьшении E до порогового значения и далее, происходит кризис аттрактора – границы аттрактора и его области притяжения соприкасаются, после чего появляется «дырка». При попадании в «дырку» изображающей точки начинается уход последней на бесконечность. Затем, при превышении переменной y некоторого порогового значения ycr, начинается возвращение на аттрактор по координате x. Когда изображающая точка достигает аттрактора (и более не покидает его) по координате x, начинается возвращение по координате y.

Далее, энергия E рассматривается не как параметр, но как медленная переменная. Пока изображающая точка движется в пределах аттрактора, E медленно убывает, сходясь к некоторой неподвижной точке E* < Ecr. При этом, когда условие E > Ecr нарушается, происходит кризис аттрактора, изображающая точка начинает уходить на бесконечность, и вместе с этим начинает расти E. Как только условие E > Ecr вновь начинает выполняться, «дырка» закрывается, и аттрактор становится локально притягивающим (росту значений быстрых переменных это не препятствует). После остановки роста пика по координате x, рост энергии также прекращается, и она вновь начинает медленно убывать.

Соотношение быстрого и медленного масштаба времени определяется параметром ε, задающим скорость изменения энергии. Стоит отметить, что при ε → 0 рост пика начинается при значении энергии, очень мало отличающемся от критического (E → Ecr–0). Поскольку «дырка» возникает при E = Ecr и увеличивается в размерах по мере уменьшения E, постольку в данном случае размеры «дырки» оказываются малы. То есть, разброс начальных значений при росте пика мал и для энергии, и для быстрых переменных. Это не соответствует реально наблюдаемому поведению динамической системы, заданной уравнением (2).

Данная проблема решается путём некоторого усложнения механизма образования пика. При переходе энергии в область E < Ecr открывается «дырка» и начинается рост только по одной координат – x. «Дырка» по y открывается только при достижении x некоторого порогового значения xcr. Параметры отображений подбираются таким образом, чтобы переменная x достигала порогового значения за количество итераций, достаточное для образования разброса по y порядка размера аттрактора. Изображающая точка попадает в область открывшейся «дырки» по y спустя некоторое время T1. Разброс по T1 тем больше, чем больше разброс по y к моменту открытия «дырки», и чем меньше размер «дырки».

Время, которое проходит до момента достижения порогового значения по y, обозначается T2. Таким образом, после достижения xcr рост по x продолжается в течение времени T1 + T2. Если размер «дырки» по y мал, то разброс по T2 также мал, и разброс максимальных значений по x определяется разбросом по T1.

Окончательно, с учётом изложенных соображений, записывается трёхмерное отображение следующего вида:

. (4)

Здесь величины γ и γ0 – параметры связи одномерных отображений быстрых переменных, а f – кусочно-линейная функция, определяющая характер этих отображений:

. (5)

График этой функции показан на рис. 1(а). Можно видеть, что одномерное отображение, заданное этой функцией, всегда имеет неустойчивую неподвижную точку. При k > 1 появляется ещё одна неустойчивая неподвижная точка, обозначаемая x*. Эта точка представляет собой границу области притяжения аттрактора данного отображения.

Функция a(E) определяет влияние энергии на параметр a функции f, и, как следствие, на геометрию аттрактора и его области притяжения, ограничивая это влияние в области больших значений E. Величины acr, a – параметры этой функции; её график показан на рис. 1(б).

 (а)  (б)
Рис. 1. Графики функций, входящих в отображение Ершова

(а) – функция f, входящая в одномерные отображения быстрых переменных; x* – неподвижная точка.

(б) – зависимость параметра a функции f от энергии.

Динамика системы (4), по построению, подразделяется на следующие фазы. В межпиковой фазе происходит хаотическое движение на аттракторе; значения обеих быстрых переменных имеют порядок 1. Энергия медленно убывает, соответственно, значение параметра a в отображении для x увеличивается, и, следовательно, наклон правой ветви f возрастает.

Фаза выброса по координате x начинается, когда энергия, а значит, и наклон правой ветви f, достигает порогового значения. В результате этого, значение xn может оказаться в области x < x*, что означает выход изображающей точки из области притяжения аттрактора и экспоненциальный рост: xn ~ (kx+)n. Энергия быстро растёт вместе с ростом |x|, En ~ ε(kx+)n. В это же время динамика по координате y приближённо описывается отображением

,  

пока выполняется условие |xn| << 1/ε. Аттрактор этого отображения – локально притягивающий, уход с него на бесконечность невозможен, несмотря на то, что левая ветвь f стала растягивающей, когда |xn| превысило |xcr|.

Фаза выброса по y наступает, когда координата x изображающей точки достигает величин порядка 1/ε. По мере роста x, появляется необходимость учитывать добавочный член в отображении, описывающем динамику по координате y:

.  

Параметры подобраны таким образом, что отображение

,  

соответствующее xn = 1/ε, уже не имеет аттрактора. «Дырка» в аттракторе возникает при достижении добавочным членом Γn некоторого критического значения

.  

Её размеры невелики, так что изображающая точка попадает в область «дырки» не сразу, но спустя некоторое количество итераций. Всё это время добавочный член и, как следствие, геометрия «дырки», продолжает меняться, так что момент начала экспоненциального роста yn ~ (ky+)n и его начальные условия имеют некоторый разброс.

Фаза спада по x наступает, как только |y| превышает критическое значение |ycr|. Левая ветвь f в отображении по x становится сжимающей, вследствие чего |xn| начинает убывать: xn ~ (kx)n. Параметры подбираются таким образом, что |xn| убывает быстрее, чем рос до того. Рост энергии сменяется медленным уменьшением: En ~ (1–ε)n. Рост по y при этом продолжается.

Фаза распада пика начинается при переходе уменьшающегося |xn| через пороговое значение |xcr|. Левая ветвь f в отображении по y становится сжимающей, рост |yn| сменяется ещё более быстрым уменьшением: yn ~ (ky)n. Пока значение |yn| превышает пороговое, |ycr|, уменьшение |x| продолжается. Параметры должны быть подобраны так, чтобы изображающая точка успела войти по координате x в область аттрактора к моменту перехода yn через порог, когда правая ветвь f в отображении по x станет растягивающей. Уменьшаясь по модулю, yn входит в область аттрактора. Система вновь вступает в межпиковую фазу.

Общий вид происходящего в системе показан на рис. 2(а) (на рис. 2(б) для сравнения приведён пик жёсткой турбулентности, зарегистрированный в описанном в [10] эксперименте по измерению напряжённости электрического поля при локальном солнечном радиовсплеске III типа). Показанное соответствует следующим, приведённым в [3], значениям параметров:

. (6)

Данный набор значений параметров используется и во всех дальнейших расчётах.

 (а)  (б)
Рис. 2. Жёсткая турбулентность в отображении Ершова и в эксперименте

(а) – графики E(n), x(n), y(n) в полулогарифмическом масштабе, для абсолютных значений. Пунктирная линия соответствует lg |E|, сплошная линия – lg |x|, штриховая линия – lg |y|. Цифрами обозначены фазы: 1 – межпиковая, 2 – выброс по x, 3 – выброс по y, 4 – спад по x, 5 – распад пика.

(б) – экспериментальный график зависимости напряжённости электрического поля от времени для ленгмюровского солитона огибающей, связанного с локальным радиовсплеском III типа [10].

Для системы (4) аналитически (в предположении ε → 0) показано, что распределение максимумов [i] пиков по координатам x и y – степенное, то есть, имеет вид

, (7)

а распределение длительностей межпиковых интервалов – экспоненциальное:

. (8)

(здесь ax, bx, Bx, an, bn, Bn – некоторые константы). В дальнейшем, будем рассматривать эти распределения как основные характеристики режима жёсткой турбулентности.

Реконструкция отображения Ершова: случай одной переменной

Итак, пусть имеется система (4), работающая в режиме жёсткой турбулентности. Пусть система рассматривается как «чёрный ящик», скрывающий её устройство, и из трёх порождаемых ею временных рядов, En, xn и yn, делающий доступным для наблюдения только один – xn. Наша задача заключается в построении простой (модельной) системы, порождающей временной ряд с аналогичными характеристиками.

Заметим, что теоретически, в условиях поставленной задачи, можно воспользоваться теоремой Такенса для восстановления фазового портрета по одномерной наблюдаемой [3]. Однако, из-за наличия двух различных масштабов времени (наблюдаются резкие и короткие пики и длинные межпиковые интервалы, в пределах которых энергия E меняется медленно по сравнению с переменными x и y), при таком подходе потребовалась бы выборка неприемлемого объёма. Таким образом, на практике необходимо использование какого-либо упрощения.

Поскольку информация наблюдателя об исходной системе неполна, будем учитывать как динамические, так и статистические характеристики временного ряда. Следовательно, при построении модельной системы можно совместно использовать (хотя бы, для простоты) динамические и вероятностные компоненты, то есть, задействовать метод русел и джокеров. Иными словами, предполагается, что данный метод может упростить построение модельной системы, и в то же время дать возможность использовать максимум доступной информации о наблюдаемой системе. В этом состоит отличие от распространённой практики моделирования, предписывающей игнорировать как можно большую часть доступной информации, с целью упрощения модели.

Первичный анализ информации о наблюдаемой системе

Обсудим, какие характеристики временного ряда xn могли бы быть рассмотрены при построении модели. Главная особенность временного ряда – наличие редких сильных выбросов, чередующихся с длительными межпиковыми интервалами. Следовательно, среди интересующих нас характеристик временного ряда можно выделить характеристики выбросов (распределение максимумов пиков, их продолжительность, динамика роста и спада), и характеристики межпиковых интервалов (продолжительность, межпиковая динамика). Также нас может заинтересовать процесс перехода от межпиковой (фоновой) динамики к росту пика, и обратно, от спада к фоновой динамике.

Прежде всего, нам понадобится правило, позволяющее отличить значения xn, относящиеся к фазе выброса, от значений xn, относящихся к межпиковым интервалам. На рис. 2 можно видеть, что первые превышают вторые на много порядков величины. В то же время, длительность выброса крайне мала по сравнению с длительностью межпикового интервала. Отсюда следует, что медиана xм временного ряда {|xn|} даст оценку характерной величины переменной состояния в межпиковой фазе. Для выделения пиковых значений xn может использоваться пороговое значение – некоторая величина xп > 0, превышающая медиану на один или несколько порядков. Надпороговые значенияxn, то есть, такие, что (|xn| > xп, считаются относящимися к фазе выброса.

Согласно (7), эмпирическая функция плотности распределения максимумов пиков в логарифмических координатах (ln (–x), ln ρ) должна (в идеальном случае) представлять собой прямую. Оценив свободный член c0 и угловой коэффициент c1 соответствующей линейной функции, например, методом наименьших квадратов, мы получим оценки параметров ax и bx:

(9)

(здесь и далее надстрочный знак «» используется для обозначения оценки).

Чтобы оценить Bx, воспользуемся условием нормировки

 

(предполагается выполнение условия сходимости интеграла bx > 1). Решив уравнение относительно Bx и подставив оценки и , имеем

(10)

Точно так же, согласно (8), построив эмпирическую функцию плотности распределения длительностей межпиковых интервалов в полулогарифмических координатах (n, ln ρ) и воспользовавшись формулами, полностью аналогичными (9), получим оценки и . Решив уравнение нормировки

 

относительно Bn и подставив и имеем оценку

. (11)

Рост пика, согласно (4) и (5), в области больших |x| описывается отображением

. (12)

Здесь мы пренебрегли членом, содержащим yn. В силу его сравнительной малости при заданном наборе значений параметров (6), а также ненаблюдаемости переменной y по условиям задачи, воздействие данного члена можно считать слабой флуктуацией. Отображение (12) можно переписать в виде

, (13)

где x* – его неподвижная точка. При отсутствии флуктуаций эта точка ограничивает область роста справа. Таким образом, точки на графике в полулогарифмических координатах (n, ln (–xn)), относящиеся к росту пика, в области больших |x| (порядка максимума пика) должны аппроксимироваться некоторой прямой. Имеем оценку коэффициента kx+:

, (14)

где c1 – угловой коэффициент данной прямой.

Спад пика в области больших |x| описывается уравнением

(15)

(начальное значение для данного отображения – точка максимума пика).

Оценку неподвижной точки x* несложно вычислить, зная оценки свободного члена c0 и углового коэффициента c1 линейной функции (12):

. (16)

Оценки c0 и c1 можно получить, аппроксимировав прямой точки, относящиеся к росту пика, на графике (xnxn+1). Данные оценки будут весьма неточными по ряду причин. Во-первых, в силу экспоненциального роста xn, для построения аппроксимирующей прямой можно использовать очень небольшое количество соседних точек (3 или 4 точки для набора значений параметров (6); некоторого повышения точности можно добиться, накапливая точки, относящиеся к разным пикам, а также заменив в формуле (16) c1 более точной оценкой ). Во-вторых, оценка по точкам |x| ~ |x*| неточна, так как в этой области возрастает влияние слабых флуктуаций в (12). В-третьих, оценка по точкам |x| >> |x*| неточна, поскольку погрешность оценки свободного члена линейной функции нарастает по мере удаления от начала координат.

С другой стороны, нас интересует не столько получение точной оценки x*, сколько воспроизведение динамики роста пика. Поэтому, оставим за собой право определять , исходя из соображений простоты построения модели.

Процедура получения оценки коэффициента kx полностью аналогична описанной выше. Роль неподвижной точки сжимающего отображения не представляется существенной.

Фоновая динамика внешне выглядит как случайное отображение некоторой ограниченной области на себя, при выходе из этой области начинается рост пика. Будем считать эту область отрезком [Dx,minDx,max]. Поскольку пик растёт в отрицательном направлении, величину Dx,max можно оценить как максимальное значение xn за время наблюдения:

. (17)

Величину Dx,min можно оценить как минимальное значение xn в промежутке между какими-либо двумя соседними пиками:

, (18)

где N1 – момент окончания спада первого пика, а N2 – момент начала роста второго.

Рис. 3. Фоновая динамика, оценки границ области фоновой динамики
График xn(n) в межпиковой фазе. Цифрами обозначены точки: 1 – точка с ординатой согласно (17), 2 – точка с ординатой согласно (18), 3 – точка с абсциссой n = N0, 4 – точка с абсциссой n = N1 согласно (19), 5 – точка с абсциссой n = N3, 6 – точка с абсциссой n = N2 согласно (20). Пунктирной линией обозначен порог |x| = xп.

Для определения N1 и N2 можно поступить следующим образом. Рассматривая поведение системы в области больших |x| (|x| >> xп), нетрудно выделить некоторый момент N0, относящийся к спаду первого пика. Считая, что в процессе спада значение xn непрерывно возрастает, N1 можно оценить как первый, начиная с N0, момент нарушения этой закономерности:

. (19)

Аналогично, выделяется некоторый момент N3, относящийся к росту второго пика, и определяется N2:

. (20)

Фоновая динамика и процедура получения оценок и показана на рис. 3. Для повышения точности оценки Dx,min можно повторить описанную процедуру для нескольких межпиковых интервалов.

Предварительные соображения по схеме русел и джокеров

Итак, анализируя временной ряд xn, можно выделить:

  • фазу фоновой динамики,
  • событие перехода от фоновой динамики к росту пика,
  • фазу роста,
  • событие остановки роста (достижения максимума пика),
  • фазу спада,
  • событие перехода от спада к фоновой динамике.

Динамика в фазах роста и спада описывается простыми уравнениями (13) и (15). Исходя из этого, имеет смысл назначить для моделирования роста и спада два русла, C1 и C2, соответственно.

Фоновая динамика хаотична, а момент перехода от фоновой динамики к росту пика – случайная величина (переменная E, инициирующая переход, не наблюдается). Следовательно, имеет смысл использовать джокер для их описания. Обозначим данный джокер J1. Как будет видно далее, этот же джокер можно использовать и для описания события перехода от спада к фоновой динамике.

Достижение максимума пика – также случайное событие (смена динамики не следует из уравнения роста; инициирующая её переменная y не наблюдается). Для его описания используем джокер, обозначенный J2.

Попутно заметим, что алгоритмы моделирования случайных событий должны быть рассчитаны на использование стандартного датчика случайных чисел, показания которого считаются случайной величиной, равномерно распределённой на отрезке [0; 1]. В дальнейших выкладках мы будем обозначать её Z, а соответствующую функцию плотности распределения –

. (21)

Русла C1 и C2

Динамика на русле C1 задаётся уравнением роста пика (13). Из этого уравнения следует, что рост пика не может начинаться с точки x = x*. Следовательно, имеет смысл ограничить область русла неравенством

.  

В некоторой окрестности x* влияние слабых флуктуаций является определяющим, и динамика хаотична. Область роста пика примыкает к области фоновой динамики слева, и с точки зрения наблюдателя, данная окрестность x* относится к области фоновой динамики. Фактически, точка x = Ax ограничивает окрестность слева, поэтому можно считать, что

. (22)

Будем полагать, что начальное значение x0, согласно намеченной схеме русел и джокеров, представляет собой результат работы джокера J1 – случайную величину. Функцию плотности распределения x0 обозначим .

Динамика на русле C2 задаётся уравнением спада (15). Поскольку рост и спад пика происходят в одной и той же области фазового пространства, область русла C2 будет совпадать с областью русла C1. Таким образом, требуется уточнять, какое из русел в данный момент считается действующим.

Джокер J2

Разработка алгоритма джокера подразумевает решение двух задач. Во-первых, необходимо задать правило выбора момента срабатывания джокера. Во-вторых, требуется определить саму процедуру срабатывания.

Вторая задача в данном конкретном случае представляется более простой. Результат срабатывания джокера J2 должен заключаться в том, что экспоненциальный рост значения наблюдаемой переменной сменяется экспоненциальным спадом. Таким образом, процедура срабатывания представляет собой изменение уравнения отображения (в терминах намеченной схемы русел и джокеров – переключение между руслами C1 и C2). Здесь просматривается аналогия с предложенным в [7] джокером типа «шов», также меняющим уравнение движения без перемещения изображающей точки.

Разница между джокером типа «шов» и джокером J2 должна состоять в определении момента срабатывания. Первый срабатывает при пересечении фазовой траекторией линии, представляющей собой область джокера. В случае второго, это простое правило неприменимо. Координата изображающей точки в момент срабатывания джокера J2 представляет собой максимум пика – случайную величину. Продолжительность роста пика – также случайная величина, что не даёт возможности «настроить» джокер на срабатывание в фиксированные моменты времени.

Таким образом, момент срабатывания джокера J2 должен определятся некоторым вероятностным предикатом. Его аргументами могут быть как время n, так и координата x. Остановимся на определении момента срабатывания исходя из координаты, так как этот способ в большей степени соответствует «традиционному» понятию области джокера.

В результате, мы пришли к следующему алгоритму. В процессе роста пика на каждой итерации принимается решение: либо срабатывает джокер и рост пика останавливается (пик достигает максимального значения), либо рост продолжается в соответствии с уравнением русла C1. Вероятность срабатывания джокера зависит от значения x на данной итерации, p = p2(x). Наша задача сводится к поиску функции p2(x), обеспечивающей распределение максимумов, соответствующее распределению (7). По этой причине мы будем ссылаться на (7) как на целевое распределение.

Рис. 4. Рост пика, срабатывание джокера и спад пика в логарифмической системе координат

Штриховая линия обозначает отображение, описывающее рост пика (движение в русле C1).

Штрихпунктирная линия обозначает отображение, описывающее спад пика (движение в русле C2).

Верхним треугольником отмечены значения переменной состояния в процессе роста.

Квадратом отмечена точка срабатывания джокера.

Нижним треугольником отмечены значения переменной состояния в процессе спада.

Для удобства дальнейших выкладок, перейдём к логарифмической системе координат:

. (23)

Искомую зависимость вероятности срабатывания джокера от координаты далее будем обозначать p̃2(x̃). Удобство (23) состоит в том, что отображение, описывающее рост пика, примет очень простой вид (рис. 4)

. (24)

Найдём вид целевого распределения в новых координатах. Для этого нам понадобится формула распределения функции случайной величины. Пусть имеется некоторая случайная величина V с функцией плотности распределения ρV(v). Пусть задана некоторая детерминированная функция g, непрерывная, строго монотонная и имеющая обратную функцию g–1 на всей области рассмотрения. Функция от случайной величины W = g(V) будет случайной величиной с функцией плотности распределения

. (25)

Сопоставим xmax величину V, x̃maxW, а функцию (23) – g. Отметим, что в области x ≤ Bx ≤ x*, во-первых, функция g отвечает всем предъявляемым требованиям, а во-вторых, величиной x* в выражении (23) можно пренебречь в силу неравенства |x*| << |Bx|. Отсюда:

.  

Далее, можно видеть, что мы перешли от степенного распределения к экспоненциальному:

. (26)

В дальнейшем нам понадобится соотношение коэффициентов, которое нетрудно получить из уравнения нормировки:

. (27)

Отметим следующее важное свойство отображения (24). Пусть дан некоторый интервал [Q̃xQ̃x+s) , Q̃x = const, s ≤ 1. Очевидным образом, он не пересекается со своим образом [Q̃x+1; Q̃x+s+1) и прообразом [Q̃x–1; Q̃x+s–1). Обозначим условную функцию плотности распределения значений x̃ (условием является попадание в данный интервал в процессе роста пика). Из (25) следует, что условная функция плотности распределения на i-образах(прообразах) [Q̃x+iQ̃x+s+i) данного интервала

.  

Интервал [Q̃x+uQ̃x+s+u), u ∈ R, в общем случае, не является образом (прообразом) [Q̃xQ̃x+s). Однако, при условии s = 1 первый интервал можно разбить на два подынтервала – образы(прообразы) соответствующих подынтервалов второго интервала с числом итераций, различающимся на единицу. Следовательно, обозначив , можно записать соотношение

. (28)

В дальнейших выкладках будем предполагать, что в процессе роста пика изображающая точка гарантированно попадает в интервал [B̃xB̃x+1), и, как следствие, функция плотности распределения значений x̃ на данном интервале совпадает с . Для этого необходимо, чтобы начальные значения находились полностью в области x̃ < B̃x, то есть,

 

(хотя бы, с некоторой приемлемой точностью). Тогда функция плотности распределения максимумов пиков в интервале x̃ ∈ [B̃xB̃x+1) будет равна

. (29)

Согласно (28), для образов указанного интервала можно записать последовательность уравнений следующего вида:

. (30)

Здесь мы учли то обстоятельство, что изображающая точка попадает в следующий интервал, только если джокер не срабатывает на данном интервале. Общий вид функции показан на рис. 5.

Рис. 5. Типичный вид распределения максимумов пиков для джокера J2
На графике показан общий вид распределения максимумов пиков, обеспечиваемый джокером J2 в соответствии с формулами (29) и (30). Для наглядности, в качестве и p̃2(x̃) взяты произвольные функции.

Отсюда, имеем уравнение

.

Учитывая, что из (26) следует

,  

получаем для p̃2(x̃) рекуррентное соотношение:

,  

или, что эквивалентно,

. (31)

Будем искать решение p̃2() этого уравнения в простейшем виде:

.  

Подставив в (31), получим выражение для искомой константы:

. (32)

Проверим, что найденное нами решение не нарушает условие нормировки

. (33)

Из (26) и (29) следует выражение для функции плотности распределения

. (34)

Подставив его в (33), получим уравнение:

.  

Преобразовав с учётом (27), получим необходимое условие для (33)

.  

Можно видеть, что (32) не противоречит данному условию.

Теперь, рассмотрим процесс роста пика в целом. В силу определённой свободы выбора границы Ãx русла C1, а также распределения начальных значений x̃0, положим, что все возможные начальные значения находятся в интервале [Ãxx+1). Отсюда, с учётом (28), (27) и (34), следует выражение для распределения начальных значений:

. (35)

Записав для интервала [ÃxÃx+1) последовательность уравнений, аналогичную (29)–(30), с ограничением x̃ + n < B̃x, и сравнив с (26), можно доопределить p̃2() на интервале [ÃxB̃x) нулём. Окончательно, формально распространив область джокера J2 на всё русло C1 и выполнив обратное преобразование координат, имеем

. (36)

Джокер J1

Прежде всего, на джокер J1 возлагается ответственность за моделирование фоновой динамики и за переход от фоновой динамики к росту пика. Сначала, рассмотрим данные задачи по отдельности.

Пусть моделирование фоновой динамики заключается в случайном отображении отрезка [Dx,minDx,max] на себя. Для простоты положим, что xn на каждой итерации присваивается значение некоторой случайной величины, функция плотности распределения которой полностью находится в данном отрезке (то есть, значения xn независимы в совокупности). С той же целью используем случайную величину Z с распределением (21). Итак, в межпиковой фазе

. (37)

Следующая задача джокера J1 – обеспечение длительности межпиковой фазы в соответствии с экспоненциальным распределением (8). Пусть на каждой n-й итерации возможно прекращение межпиковой фазы джокером, с вероятностью, определяемой некоторой функцией p1(n). Чтобы найти эту функцию, рассмотрим пуассоновское приближение биномиального распределения

,  

где Pm(n) – вероятность ровно m успехов в n испытаниях Бернулли, а p – вероятность успеха в отдельном испытании. Известно, что промежуток (количество испытаний) между успешными испытаниями распределён экспоненциально:

.  

Следовательно, взяв вероятность прекращения джокером межпиковой фазы на данной итерации p = bn, мы получим искомый вид распределения межпиковых интервалов [ii]. Окончательно, доопределив функцию p1(n) нулём в начале межпиковой фазы, получим

. (38)

Теперь, рассмотрим переход от фоновой динамики к росту пика. В этот момент, джокер J1, вместо выполнения очередной итерации согласно (37), должен перевести изображающую точку в исток русла C1. Таким образом, имеет смысл возложить на J1 обязанность формирования в истоке русла распределения начальных значений вида (35), необходимого для получения заданного распределения максимумов пиков.

Для этого требуется решить частную задачу поиска зависимости x̃0 = gZ(Z), позволяющей получить распределение (35) на основе (21). Согласно (25), имеем уравнение

.  

Предположив, что выражение под знаком модуля положительно на всём интервале [ÃxÃx + 1), перейдём к следующему уравнению:

.  

Проинтегрировав, получим

.  

Можно видеть, что на интервале [ÃxÃx+1) полученная функция строго возрастает, что соответствует сделанному выше предположению относительно положительности выражения под знаком модуля (рис. 6).

Рис. 6. Получение распределения начальных значений на основе равномерного распределения

1 – функция плотности распределения (21) случайной величины Z.

2 – искомая зависимость (39) начального значения x̃0 от случайной величины Z.

3 – заданная функция плотности распределения (35) начального значения x̃0.

Чтобы найти константу интегрирования C, запишем граничные условия:

.  

Отметим, что одно из граничных условий избыточно, и, следовательно, может использоваться для проверки непротиворечивости решения. Имеем:

.  

Окончательно, выполнив обращение найденной функции, получим

,  

или, после обратного преобразования координат и подстановки определения α̃x из (35),

. (39)

Теперь можно определить значение таким образом, чтобы воспроизвести динамику в начале фазы роста пика. Как было отмечено выше, получение точной оценки x* затруднено, и задача заключается скорее в воспроизведении внешнего вида фронта пика. Поэтому, имеет смысл при определении оценки пользоваться наиболее грубыми и простыми методами.

Рассмотрим начало роста некоторого пика. Момент выхода xn из области [Dx,minDx,max] обозначим Nстарт:

. (40)

Момент превышения абсолютной величиной xn порогового значения xп обозначим Nп:

. (41)

Число итераций между этими двумя моментами обозначим Nфронт:

. (42)

Среднее по всем наблюдаемым пикам значение Nфронт обозначим ⟨Nфронт⟩. Воспользуемся обозначением ⟨Nфронт,м⟩ для модельной системы и ⟨Nфронт,н⟩ для наблюдаемой системы. Будем считать, что внешний вид фронта пика воспроизводится приемлемо, если

.  

Значение ⟨Nфронт⟩ модельной системы можно грубо оценить, воспользовавшись аналогичным (13) равенством

, (43)

где ⟨x0⟩ обозначает математическое ожидание начального значения. Чтобы найти последнее, выполним обратное (23) преобразования координат для распределения начальных значений (35). Согласно (25), имеем

.  

Отсюда,

, (44)

где

. (45)

Приравняв ⟨Nфронт,н⟩ и ⟨Nфронт,м⟩, подставив (22) и (44) в (43) и решив уравнение относительно x*, получим оценку

. (46)

Наконец, ещё одна потенциальная обязанность джокера J1 – прекращение спада пика согласно уравнению русла C2. Фаза спада считается законченной, когда изображающая точка достигает области фоновой динамики. По этой причине, а также с учётом написанного выше, имеет смысл объявить отрезок [Dx,minDx,max] областью джокера J1. Таким образом, двигаясь в соответствии с уравнением русла C2, изображающая точка покинет русло и перейдёт в область джокера, после чего произойдёт активизация последнего.

Построение системы русел и джокеров

Окончательно, цикл работы системы русел и джокеров выглядит следующим образом. В межпиковом интервале в области джокера J1 выполняются итерации случайного отображения (37) в течение времени, определяемого вероятностью начала выброса (38). Выброс начинается с задания в истоке русла C1 джокером J1 начального значения согласно (39). После этого выполняются итерации соответствующего руслу C1 отображения (13). При этом, вероятность прекращения роста определяется джокером J2 согласно (36). В момент прекращения роста русло C1 перестаёт считаться действующим, активизируется русло C2, и начинают выполняться итерации соответствующего ему отображения (15). Это продолжается до тех пор, пока изображающая точка не достигнет области джокера J1, после чего начинается новый межпиковый интервал (рис. 7).

Рис. 7. Система русел и джокеров

Прямой штриховкой обозначены совпадающие области русел C1 и C2 и джокера J2.

Косой штриховкой обозначена область джокера J1.

Верхним треугольником отмечены значения переменной состояния в процессе роста пика.

Квадратиком отмечена точка срабатывания джокера J2.

Нижним треугольником отмечены значения переменной состояния в процессе спада пика.

Цифрами обозначены: 1 – срабатывание джокера J1: случайное перемещение изображающей точки в соответствии с уравнением (37), 2 – срабатывание джокера J1: перевод изображающей точки в исток русла C1, 3 – движение в русле C1 в соответствии с его уравнением (13), 4 – срабатывание джокера J2: перевод изображающей точки в русло C2, 5 – движение в русле C2 в соответствии с его уравнением (15), 6 – изображающая точка вновь попадает в область джокера J1, движение в соответствии с уравнением русла C2 прекращается.

В соответствии с постановкой задачи, исходные данные для построения системы представляют собой временной ряд xn (предполагается, что данных достаточно для получения всех необходимых статистических оценок).

Ниже приводится предлагаемый порядок построения системы русел и джокеров.

  1. Вычисляется медиана xм временного ряда {|xn|} и пороговое значение xп.
  2. Определяются границы областей русел и джокеров:
    1. по формуле (17) вычисляется оценка (используется весь имеющийся временной ряд, либо данные по одному или нескольким межпиковым интервалам);
    2. по надпороговым значениям xn выделяются два следующих непосредственно друг за другом пика;
    3. выделяются соответствующие надпороговым значениям моменты N0 и N3 в фазах спада первого пика и роста второго, соответственно;
    4. по формуле (19), с использованием значений xn, начиная с n = N0, определяется момент N1 начала межпикового интервала;
    5. по формуле (20), с использованием значений xn, начиная с n = N3, в сторону уменьшения n, определяется момент N2 конца межпикового интервала;
    6. по формуле (18) вычисляется оценка . Это же значение, согласно (22), является правой границей области русел Ax;
    7. возможен повтор данной процедуры для остальных межпиковых интервалов, с целью уточнения оценок;
  3. Вычисляются оценки , и :
    1. Строится гистограмма длительностей межпиковых интервалов в полулогарифмических координатах (n, ln r);
    2. Для построенной гистограммы прокладывается аппроксимирующая прямая;
    3. Вычисляются оценки и путём подстановки коэффициентов аппроксимирующей прямой в формулу, аналогичную (9);
    4. Вычисляется оценка путём подстановки полученных оценок и в (11).
  4. Вычисляются оценки , и :
    1. строится гистограмма максимумов пиков xmax в логарифмических координатах (ln (–x), ln ρ);
    2. для построенной гистограммы прокладывается аппроксимирующая прямая;
    3. вычисляются оценки и путём подстановки коэффициентов аппроксимирующей прямой в (9);
    4. вычисляется оценка путём подстановки полученных оценок и в (10).
  5. Вычисляются оценки и :
    1. для одного из пиков строится график в полулогарифмических координатах (n, ln (–xn));
    2. прокладывается аппроксимирующая прямая по точкам в области больших |x| (порядка ), относящимся к фазе роста пика;
    3. вычисляется оценка путём подстановки углового коэффициента аппроксимирующей прямой в (14);
    4. прокладывается аппроксимирующая прямая по точкам в области больших |x|, относящимся к фазе спада пика;
    5. вычисляется оценка путём подстановки углового коэффициента аппроксимирующей прямой в формулу (14);
    6. Возможен повтор данной процедуры для других пиков, с целью уточнения оценок.
  6. Вычисляется оценка :
    1. выполняется подстановка оценки в правило (40). Согласно данному правилу, на основе значений xn начиная с момента n = N2 конца одного из межпиковых интервалов, определяется значение Nстарт;
    2. определяется значение Nп по правилу (41) и вычисляется значение Nфронт по формуле (42);
    3. аналогичным образом, определяются значения Nфронт для остальных пиков. Вычисляется среднее значение ⟨Nфронт,н⟩;
    4. вычисляется значение c0 путём подстановки оценок и в формулу (45);
    5. вычисляется оценка путём подстановки значений c0, ⟨Nфронт,н⟩ и xп и оценок и в формулу (46).
  7. Вычисляются коэффициенты Ãx, , , и вероятность :
    1. Ãx и – путём подстановки x = Ax и в (23);
    2. и – путём подстановки , и в (26);
    3. – путём подстановки оценки в (32).
  8. Окончательно определяются уравнения схемы русел и джокеров:
    1. уравнение (13) русла C1 – путём подстановки и ;
    2. уравнение (15) русла C2 – путём подстановки ;
    3. функция отображения (37) джокера J1 – путём подстановки и ;
    4. функция вероятности прекращения межпиковой фазы (38) джокера J1 – путём подстановки и ;
    5. функция формирования джокером J1 начального значения для роста пика (39) – путём подстановки , Ãx, и ;
    6. функция вероятности срабатывания джокера J2 (36) – путём подстановки Ax и .

Результаты моделирования, сравнение

Для проверки предложенной модели использовался временной ряд xn, порождаемый системой (4) с набором параметров (6) на протяжении 5·107 итераций. За время наблюдения произошёл 8221 выброс.

На основе собранной статистики были получены следующие оценки:

.

Вероятностные параметры джокеров J1 и J2 составили pn = 2,5·10–4 и px = 5,35·10–2, соответственно.

Общий вид динамики модельной системы J1C1J2C2 показан на рис. 8 (ср. рис. 2).

 (а)  (б)
Рис. 8. Динамика системы русел и джокеров

(а) – фоновая динамика и одиночный пик в системе J1C1J2C2, полулогарифмический масштаб.

(б) – фронт пика, подробно. Сплошная линия соответствует системе J1C1J2C2, пунктирная линия – наблюдаемой системе. Штрихпунктирной линией показан уровень , штриховой – уровень x = xп.

Гистограммы абсолютных величин максимумов пиков и длительностей межпиковых интервалов в системах (4) и J1C1J2C2 показаны на рис. 9. Можно видеть, что модельная система соответствует наблюдаемой (исходной) системе по выбранным для сравнения статистическим характеристикам.

 (а)  (б)
Рис. 9. Сравнение статистических характеристик исходной и модельной систем

(а) – гистограммы максимумов пиков, логарифмический масштаб.

(б) – гистограммы длительностей межпиковых интервалов, полулогарифмический масштаб.

Сплошная линия соответствует системе J1C1J2C2, пунктирная линия – наблюдаемой системе.

Заключение

Итак, рассмотрен пример построения модели системы с жёсткой турбулентностью в условиях неполноты доступной информации. Исходные данные представляют собой один из трёх временных рядов, порождаемых трёхмерным отображением Ершова (4). Характерная особенность рассмотренной системы заключается в наличии двух различных масштабов времени, что делает использование теоремы Такенса невозможным на практике. С учётом этого, при решении задачи использовался метод русел и джокеров, заключающийся в совмещении динамических и вероятностных методов моделирования.

Наблюдаемый временной ряд отличается наличием длительных промежутков хаотической динамики, чередующихся с сильными выбросами (экспоненциальный рост, сменяющийся экспоненциальным спадом). В качестве главного показателя соответствия системы русел и джокеров исходной системе выбрано сходство распределений максимумов пиков и длительностей межпиковых интервалов порождаемых временных рядов.

Предложена схема, включающая два русла, C1 и C2, и два джокера, J1 и J2. Уравнения русел описывают экспоненциальный рост и спад. Области русел совпадают, поэтому специально оговаривается, какое из русел в данный момент считается действующим.

Джокер J1 отвечает за моделирование фоновой динамики и за переключение между различными типами динамики. Фоновая динамика моделируется путём случайного отображения области джокера на себя (используется равномерно распределённая случайная величина). Выброс начинается в результате перевода джокером изображающей точки в исток русла C1 (после этого русло C1 объявляется действующим). Спад заканчивается после достижения изображающей точкой, следующей в русле C2, области джокера J1. Область джокера смежна с областью русел.

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

Следует подчеркнуть две особенности предложенной схемы. Во-первых, она предусматривает совпадение областей русел и джокеров, и в этом смысле представляет собой дальнейшее развитие схемы, предложенной в [7]. Во-вторых, в отличие от предлагавшихся ранее схем, наибольшую часть времени изображающая точка проводит не в области русел, а в области джокера.

В дальнейшем возможна доработка джокера J1 для более точного моделирования фоновой динамики (например, воспроизведение распределения значений наблюдаемой величины, её функции автокорреляции,…). Также возможно рассмотрение более сложной задачи, когда для наблюдения доступны две переменные из трёх, x и y (здесь потребуется не только обеспечение заданного вида распределения для каждой переменной, но и соблюдение соотношения между параметрами распределений, учёт распределения запаздывания выбросов по переменной y относительно выбросов по переменной x, учёт кросс-корреляции переменных в межпиковой фазе, и т.д.).

Литература

  1. Kuramoto Y. Chemical Oscillations, Waves, and Turbulence. – Berlin: Springer-Verlag, 1984. – 156 p.
  2. Ахромеева Т.С., Курдюмов С.П., Малинецкий Г.Г., Самарский А.А. Нестационарные структуры и диффузионный хаос. – М.: Наука, 1992. – 544 с.
  3. Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. – М.: УРСС, 2002. – 360 с.
  4. Малинецкий Г.Г., Потапов А.Б., Подлазов А.В. Нелинейная динамика: подходы, результаты, надежды. – М.: УРСС, 2006. – 280 с.
  5. Белайчук Л.В., Малинецкий Г.Г. Проделки джокеров на одномерных отображениях. – М.: ИПМ им. М.В. Келдыша, 1997. Препринт №24. – 29 с.
  6. Малинецкий Г.Г., Потапов А.Б. Русла и джокеры: о новых методах прогноза поведения сложных систем. – М.: ИПМ им. М.В. Келдыша, 1998. Препринт №32. – 32 с.
  7. Зульпукаров М.-Г.М., Малинецкий Г.Г., Подлазов А.В. Метод русел и джокеров на примере исследования системы Розенцвейга–Макартура. – М.: ИПМ им. М.В. Келдыша, 2006. Препринт №21. – 32 с.
  8. Гарел Д., Гарел О. Колебательные химические реакции. – М.: Мир, 1986. – 148 с.
  9. Владимиров В.А., Воробьев Ю.Л., Малинецкий Г.Г. и др. Управление риском. Риск, устойчивое развитие, синергетика/ Кибернетика: неограниченные возможности и возможные ограничения. – М.: Наука, 2000. – 432 с.
  10. Thejappa G., MacDowall R.J. Ulysses Observations of Nonlinear Wave-wave Interactions in the Source Regions of Type III Solar Radio Bursts // Journal of Astrophysics and Astronomy. – 2000, September, December. – Vol.21. – №3,4. – P.447–450.

 



[i] Здесь и далее слова «максимум пика» используются для обозначения максимума по абсолютной величине.

[ii] Замечание. В на первый взгляд аналогичном случае джокера J2 целевое распределение зависит не только от p2(x), но и от распределения начальных значений x0. Поэтому, данная схема, предположительно, не даст существенного уменьшения объёма выкладок.