Показатель Ляпунова

редактировать
Пояснения к показателю Ляпунова

В математике, то показатель Ляпунова или Ляпунова характеристический показатель из динамической системы является величиной, которая характеризует скорость разделения, бесконечно близких траекторий. Количественно две траектории в фазовом пространстве с начальным вектором разделения расходятся (при условии, что расходимость можно рассматривать в рамках линеаризованного приближения) со скоростью, определяемой выражением δ Z 0 {\ displaystyle \ delta \ mathbf {Z} _ {0}}

| δ Z ( т ) | е λ т | δ Z 0 | {\ displaystyle | \ delta \ mathbf {Z} (t) | \ приблизительно e ^ {\ lambda t} | \ delta \ mathbf {Z} _ {0} |}

где - показатель Ляпунова. λ {\ displaystyle \ lambda}

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

Экспонента названа в честь Александра Ляпунова.

СОДЕРЖАНИЕ

  • 1 Определение максимального показателя Ляпунова
  • 2 Определение спектра Ляпунова
  • 3 показатель Ляпунова для нестационарной линеаризации
    • 3.1 Эффекты Перрона инверсии знака наибольшего показателя Ляпунова
  • 4 Основные свойства
  • 5 Значение спектра Ляпунова
  • 6 Численный расчет
  • 7 Локальный показатель Ляпунова
  • 8 Условный показатель Ляпунова
  • 9 См. Также
  • 10 Ссылки
  • 11 Дальнейшее чтение
  • 12 Программное обеспечение
  • 13 Внешние ссылки

Определение максимального показателя Ляпунова

Максимальный показатель Ляпунова можно определить следующим образом:

λ знак равно Lim т Lim | δ Z 0 | 0 1 т пер | δ Z ( т ) | | δ Z 0 | {\ displaystyle \ lambda = \ lim _ {t \ to \ infty} \ lim _ {| \ delta \ mathbf {Z} _ {0} | \ to 0} {\ frac {1} {t}} \ ln { \ frac {| \ delta \ mathbf {Z} (t) |} {| \ delta \ mathbf {Z} _ {0} |}}}

Предел гарантирует справедливость линейного приближения в любое время. | δ Z 0 | 0 {\ displaystyle | \ delta \ mathbf {Z} _ {0} | \ to 0}

Для системы с дискретным временем (карты или итерации с фиксированной точкой) для орбиты, начинающейся с этого, это означает: Икс п + 1 знак равно ж ( Икс п ) {\ Displaystyle х_ {п + 1} = е (х_ {п})} Икс 0 {\ displaystyle x_ {0}}

λ ( Икс 0 ) знак равно Lim п 1 п я знак равно 0 п - 1 пер | ж ( Икс я ) | {\ displaystyle \ lambda (x_ {0}) = \ lim _ {n \ to \ infty} {\ frac {1} {n}} \ sum _ {i = 0} ^ {n-1} \ ln | f '(x_ {i}) |}

Определение спектра Ляпунова.

Показатель Ляпунова Ведущий вектор Ляпунова.

Для динамической системы с уравнением эволюции в n -мерном фазовом пространстве спектр показателей Ляпунова Икс ˙ я знак равно ж я ( Икс ) {\ Displaystyle {\ точка {х}} _ {я} = е_ {я} (х)}

{ λ 1 , λ 2 , , λ п } , {\ displaystyle \ {\ lambda _ {1}, \ lambda _ {2}, \ ldots, \ lambda _ {n} \} \,,}

в общем, зависит от отправной точки. Однако нас обычно интересует аттрактор (или аттракторы) динамической системы, и обычно с каждым аттрактором связан один набор показателей. Выбор начальной точки может определить, на каком аттракторе окажется система, если их больше одного. (Для гамильтоновых систем, не имеющих аттракторов, это не проблема.) Показатели Ляпунова описывают поведение векторов в касательном пространстве фазового пространства и определяются из матрицы Якоби Икс 0 {\ displaystyle x_ {0}}

J я j ( т ) знак равно d ж я ( Икс ) d Икс j | Икс ( т ) {\ displaystyle J_ {ij} (t) = \ left. {\ frac {df_ {i} (x)} {dx_ {j}}} \ right | _ {x (t)}}

этот якобиан определяет эволюцию касательных векторов, заданных матрицей, через уравнение Y {\ displaystyle Y}

Y ˙ знак равно J Y {\ displaystyle {\ dot {Y}} = JY}

с начальным условием. Матрица описывает, как небольшое изменение в точке распространяется до конечной точки. Лимит Y я j ( 0 ) знак равно δ я j {\ displaystyle Y_ {ij} (0) = \ delta _ {ij}} Y {\ displaystyle Y} Икс ( 0 ) {\ Displaystyle х (0)} Икс ( т ) {\ Displaystyle х (т)}

Λ знак равно Lim т 1 2 т бревно ( Y ( т ) Y Т ( т ) ) {\ displaystyle \ Lambda = \ lim _ {t \ rightarrow \ infty} {\ frac {1} {2t}} \ log (Y (t) Y ^ {T} (t))}

определяет матрицу (условия существования предела даются теоремой Оселедца ). Показатели Ляпунова определяются собственными значениями. Λ {\ displaystyle \ Lambda} λ я {\ displaystyle \ lambda _ {я}} Λ {\ displaystyle \ Lambda}

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

Показатель Ляпунова для нестационарной линеаризации

Чтобы ввести показатель Ляпунова, рассмотрим фундаментальную матрицу (например, для линеаризации вдоль стационарного решения в непрерывной системе фундаментальная матрица состоит из линейно независимых решений первого приближения системы. Сингулярные значения матрицы равны квадратные корни из собственных значений матрицы.Наибольший показатель Ляпунова имеет вид Икс ( т ) {\ Displaystyle X (т)} Икс 0 {\ displaystyle x_ {0}} exp ( d ж т ( Икс ) d Икс | Икс 0 т ) {\ displaystyle \ exp \ left (\ left. {\ frac {df ^ {t} (x)} {dx}} \ right | _ {x_ {0}} t \ right)} { α j ( Икс ( т ) ) } 1 п {\ displaystyle \ {\ alpha _ {j} {\ big (} X (t) {\ big)} \} _ {1} ^ {n}} Икс ( т ) {\ Displaystyle X (т)} Икс ( т ) * Икс ( т ) {\ Displaystyle Х (т) ^ {*} Х (т)} λ м а Икс {\ displaystyle \ lambda _ {\ mathrm {max}}}

λ м а Икс знак равно Максимум j лим суп т 1 т пер α j ( Икс ( т ) ) . {\ displaystyle \ lambda _ {\ mathrm {max}} = \ max \ limits _ {j} \ limsup _ {t \ rightarrow \ infty} {\ frac {1} {t}} \ ln \ alpha _ {j} {\ big (} X (t) {\ big)}.}

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

Эффекты Перрона инверсии знака наибольшего показателя Ляпунова

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

Контрпример Перрона показывает, что наибольший отрицательный показатель Ляпунова, как правило, не указывает на стабильность, и что наибольший положительный показатель Ляпунова, как правило, не указывает на хаос.

Следовательно, нестационарная линеаризация требует дополнительного обоснования.

Основные свойства

Если система консервативна (т. Е. Отсутствует диссипация ), элемент объема фазового пространства останется неизменным вдоль траектории. Таким образом, сумма всех показателей Ляпунова должна быть равна нулю. Если система диссипативна, сумма показателей Ляпунова отрицательна.

Если система является потоком и траектория не сходится к одной точке, один показатель всегда равен нулю - показатель Ляпунова, соответствующий собственному значению с собственным вектором в направлении потока. L {\ displaystyle L}

Значение спектра Ляпунова

Спектр Ляпунова можно использовать для оценки скорости производства энтропии, фрактальной размерности и размерности Хаусдорфа рассматриваемой динамической системы. В частности, из знания спектра Ляпунова можно получить так называемую размерность Ляпунова (или размерность Каплана – Йорка ), которая определяется следующим образом: D K Y {\ displaystyle D_ {KY}}

D K Y знак равно k + я знак равно 1 k λ я | λ k + 1 | {\ displaystyle D_ {KY} = k + \ sum _ {i = 1} ^ {k} {\ frac {\ lambda _ {i}} {| \ lambda _ {k + 1} |}}}

где - максимальное целое число, при котором сумма наибольших показателей все еще неотрицательна. представляет собой верхнюю границу информационного измерения системы. Более того, сумма всех положительных показателей Ляпунова дает оценку энтропии Колмогорова – Синая согласно теореме Песина. Наряду с широко используемыми численными методами оценки и вычисления размерности Ляпунова существует эффективный аналитический подход, основанный на прямом методе Ляпунова со специальными функциями типа Ляпунова. Показатели Ляпунова ограниченной траектории и ляпуновская размерность аттрактора инвариантны относительно диффеоморфизма фазового пространства. k {\ displaystyle k} k {\ displaystyle k} D K Y {\ displaystyle D_ {KY}}

Мультипликативный обратная наибольшего показателя Ляпунова иногда называют в литературе как время ляпуновской, и определяет характерное электронное -folding время. Для хаотических орбит время Ляпунова будет конечным, тогда как для обычных орбит оно будет бесконечным.

Численный расчет

Показатель Ляпунова Точки внутри и снаружи множества Мандельброта раскрашены показателем Ляпунова.

Как правило, вычисление показателей Ляпунова, как определено выше, не может быть выполнено аналитически, и в большинстве случаев приходится прибегать к численным методам. Ранний пример, который также представляет собой первую демонстрацию экспоненциального расхождения хаотических траекторий, был проведен Р. Х. Миллером в 1964 году. В настоящее время наиболее часто используемая численная процедура оценивает матрицу, основанную на усреднении нескольких приближений определения предела за конечное время. L {\ displaystyle L} L {\ displaystyle L}

Один из наиболее часто используемых и эффективных численных методов вычисления спектра Ляпунова для гладкой динамической системы основан на периодической ортонормировке Грама – Шмидта векторов Ляпунова, чтобы избежать смещения всех векторов в направлении максимального расширения.

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

Местный показатель Ляпунова

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

Условный показатель Ляпунова

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

Смотрите также

использованная литература

дальнейшее чтение

  • М.-Ф. Данка и Н.В. Кузнецовы (2018). "Код Matlab для показателей Ляпунова систем дробного порядка". Международный журнал бифуркаций и хаоса. 25 (5): ст. число 1850067. arXiv : 1804.01143. DOI : 10.1142 / S0218127418500670.
  • X. Zeng; Р. Эйкхольт и Р. А. Пильке (1991). «Оценка спектра показателей Ляпунова по коротким временным рядам низкой точности». Письма с физическим обзором. 66 (25): 3229–3232. Bibcode : 1991PhRvL..66.3229Z. DOI : 10.1103 / PhysRevLett.66.3229. PMID   10043734.

Программное обеспечение

  • [1] Р. Хеггер, Х. Канц и Т. Шрайбер, Нелинейный анализ временных рядов, TISEAN 3.0.1 (март 2007 г.).
  • [2] Продукт Scientio ChaosKit вычисляет показатели Ляпунова среди других мер Хаоса. Доступ предоставляется онлайн через веб-службу и демонстрацию Silverlight.
  • [3] Лаборатория программного обеспечения для математических воссозданий доктора Рональда Джо Рекорда включает графический клиент X11, lyap, для графического исследования показателей Ляпунова принудительной логистической карты и других карт единичного интервала. На содержание и страницы руководства по mathrec лаборатории программного обеспечения также доступны.
  • [4] Программное обеспечение на этой странице было разработано специально для эффективного и точного расчета полного спектра показателей. Это включает LyapOde для случаев, когда уравнения движения известны, а также Lyap для случаев, связанных с данными экспериментальных временных рядов. LyapOde, который включает исходный код, написанный на "C", также может вычислять условные показатели Ляпунова для связанных идентичных систем. Он предназначен для того, чтобы позволить пользователю предоставить свой собственный набор уравнений модели или использовать одно из включенных. Нет никаких внутренних ограничений на количество переменных, параметров и т. Д. Lyap, который включает исходный код, написанный на Fortran, также может вычислять векторы направления Ляпунова и может характеризовать сингулярность аттрактора, что является основной причиной трудностей при вычислении более отрицательные показатели из данных временных рядов. В обоих случаях есть обширная документация и образцы входных файлов. Программное обеспечение может быть скомпилировано для работы в системах Windows, Mac или Linux / Unix. Программное обеспечение работает в текстовом окне и не имеет графических возможностей, но может генерировать выходные файлы, которые можно легко построить с помощью такой программы, как Excel.

внешние ссылки

Последняя правка сделана 2024-01-05 08:14:29
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте