Lyapunov eksponenti - Lyapunov exponent

Yilda matematika The Lyapunov eksponenti yoki Lyapunovning o'ziga xos xususiyati a dinamik tizim cheksiz yaqin bo'linish tezligini tavsiflovchi miqdor traektoriyalar. Miqdoriy ravishda, ikkita traektoriya fazaviy bo'shliq dastlabki ajratish vektori bilan diverge (divergensiyani chiziqli yaqinlashuv ichida davolash mumkin bo'lsa) tomonidan berilgan tezlik bilan

qayerda Lyapunov eksponentidir.

Boshlang'ich ajratish vektorining turli yo'nalishlari uchun ajratish tezligi har xil bo'lishi mumkin. Shunday qilib, a Lyapunov eksponentlarining spektri- son fazasi fazasining o'lchamiga teng. Ulardan eng kattasini "." Deb atash odatiy holdir Maksimal Lyapunov eksponenti (MLE), chunki u tushunchasini belgilaydi bashorat qilish dinamik tizim uchun. Odatda tizimning belgisi sifatida ijobiy MLE olinadi tartibsiz (ba'zi bir boshqa shartlar bajarilgan taqdirda, masalan, fazoviy bo'shliqning ixchamligi). Ixtiyoriy ravishda boshlang'ich ajratish vektori odatda MLE bilan bog'liq yo'nalishdagi ba'zi tarkibiy qismlardan iborat bo'lishiga e'tibor bering va eksponent o'sish sur'ati tufayli boshqa ko'rsatkichlarning ta'siri vaqt o'tishi bilan yo'q bo'lib ketadi.

Eksponent nomi berilgan Aleksandr Lyapunov.

Maksimal Lyapunov ko'rsatkichining ta'rifi

Lyapunovning maksimal ko'rsatkichini quyidagicha aniqlash mumkin:

Chegara har qanday vaqtda chiziqli yaqinlashuvning haqiqiyligini ta'minlaydi.[1]

Diskret vaqt tizimi uchun (xaritalar yoki sobit nuqta takrorlashlari) , bilan boshlanadigan orbitaga bu quyidagicha tarjima qilinadi:

Lyapunov spektrining ta'rifi

Evolyutsiya tenglamasiga ega bo'lgan dinamik tizim uchun ichida n- o'lchovli faza maydoni, Lyapunov eksponentlari spektri

umuman, boshlang'ich nuqtaga bog'liq . Biroq, biz odatda jalb qiluvchi (yoki attraktorlar) dinamik tizim, va odatda har bir attraktor bilan bog'liq bo'lgan bitta ko'rsatkich mavjud. Boshlang'ich nuqtani tanlash, agar tizim bir nechta bo'lsa, tizim qaysi jalb qiluvchi bilan yakunlanishini aniqlashi mumkin. (Attraktsionerlari bo'lmagan Hamilton tizimlari uchun bu tashvish tug'dirmaydi.) Lyapunov eksponentlari vektorlarning faza fazosining teginish fazosidagi xatti-harakatlarini tavsiflaydi va Yakobian matritsasi

bu Jacobian matritsa tomonidan berilgan tangens vektorlarning evolyutsiyasini belgilaydi , tenglama orqali

dastlabki shart bilan . Matritsa nuqtada qanday kichik o'zgarishlarni tasvirlaydi yakuniy nuqtaga qadar tarqaladi . Chegara

matritsani belgilaydi (chegaraning mavjudligi uchun shartlar Oseledets teoremasi ). Lyapunov eksponentlari ning xos qiymatlari bilan aniqlanadi .

Lyapunov eksponentlari to'plami anning deyarli barcha boshlang'ich nuqtalari uchun bir xil bo'ladi ergodik dinamik tizimning tarkibiy qismi.

Lyapunov vaqt bo'yicha o'zgaruvchan chiziqlash uchun ko'rsatkich

Lyapunov ko'rsatkichini tanishtirish uchun asosiy matritsani ko'rib chiqing(masalan, statsionar eritma bo'ylab chiziqlash uchun uzluksiz tizimda asosiy matritsa hisoblanaditizimning birinchi darajali yaqinlashuvining chiziqli mustaqil echimlaridan iborat.Yakkalik qiymatlarimatritsaning matritsaning o'ziga xos qiymatlarining kvadrat ildizlari .Lyapunovning eng katta eksponati quyidagicha[2]

A.M. Lyapunov birinchi yaqinlashuv tizimi muntazam bo'lsa (masalan, doimiy va davriy koeffitsientli barcha tizimlar muntazam bo'lsa) va uning eng katta Lyapunov ko'rsatkichi manfiy bo'lsa, u holda asl tizimning echimi asimptotik ravishda Lyapunov barqaror.Keyingi vaqtlarda O. Perron birinchi yaqinlashuvning muntazamligi talabi katta ahamiyatga ega deb aytgan edi.

Lyapunovning eng katta ko'rsatkichli inversiyasining perron effektlari

1930 yilda O. Perron ikkinchi tartibli tizimning namunasini yaratdi, bu erda birinchi yaqinlashuv dastlabki tizimning nol eritmasi bo'yicha salbiy Lyapunov ko'rsatkichlariga ega, ammo shu bilan birga, asl chiziqli bo'lmagan tizimning bu nol eritmasi Lyapunov beqaror. Bundan tashqari, ushbu nol eritmaning ma'lum bir mahallasida asl tizimning deyarli barcha echimlari Lyapunovning ijobiy ko'rsatkichlariga ega. Shuningdek, teskari misolni yaratish mumkin, bunda birinchi yaqinlashuv dastlabki tizimning nolli eritmasi bo'yicha ijobiy Lyapunov ko'rsatkichlariga ega, ammo shu bilan birga, bu chiziqli bo'lmagan tizimning nol eritmasi Lyapunov barqaror.[3][4]Lyapunovning dastlabki tizim echimlari ko'rsatkichlari va bir xil dastlabki ma'lumotlarga birinchi yaqinlashish tizimining belgilar inversiyasining ta'siri keyinchalik Perron effekti deb nomlandi.[3][4]

Perronning qarama-qarshi namunasi shuni ko'rsatadiki, Lyapunovning eng katta salbiy ko'rsatkichi umuman barqarorlikni anglatmaydi va ijobiy eng katta Lyapunov ko'rsatkichi, umuman, betartiblikni ko'rsatmaydi.

Shuning uchun vaqt bo'yicha o'zgaruvchan chiziqlash qo'shimcha asoslashni talab qiladi.[4]

Asosiy xususiyatlar

Agar tizim konservativ bo'lsa (ya'ni, yo'q) tarqalish ), faza fazosining hajm elementi traektoriya bo'ylab bir xil bo'ladi. Shunday qilib, Lyapunovning barcha eksponentlari yig'indisi nolga teng bo'lishi kerak. Agar tizim dissipativ bo'lsa, Lyapunov ko'rsatkichlari yig'indisi salbiy bo'ladi.

Agar tizim oqim bo'lsa va traektoriya bitta nuqtaga yaqinlashmasa, bitta ko'rsatkich har doim nolga teng bo'ladi - Lyapunov ko'rsatkichi o'z qiymatiga mos keladi oqim yo'nalishi bo'yicha xususiy vektor bilan.

Lyapunov spektrining ahamiyati

Lyapunov spektridan entropiya hosil bo'lish tezligini baholash uchun foydalanish mumkin fraktal o'lchov va of Hausdorff o'lchovi ko'rib chiqilgan dinamik tizim[5]. Xususan, Lyapunov spektridan olingan ma'lumotdan shunday deb atash mumkin Lyapunov o'lchovi (yoki Kaplan-York o'lchovi ) quyidagicha belgilanadi:

qayerda ning yig'indisi bo'ladigan maksimal butun son eng katta ko'rsatkichlar hali ham salbiy emas. uchun yuqori chegarani bildiradi axborot o'lchovi tizimning.[6] Bundan tashqari, Lyapunovning barcha ijobiy ko'rsatkichlari yig'indisi $ ga baho beradi Kolmogorov - Sinay entropiyasi Pesin teoremasiga muvofiq.[7]Hisoblash va hisoblash uchun keng qo'llaniladigan raqamli usullar bilan bir qatorda Lyapunov o'lchovi Lyapunovga o'xshash maxsus funktsiyalarga ega to'g'ridan-to'g'ri Lyapunov uslubiga asoslangan samarali analitik yondashuv mavjud.[8]Lyapunov chegaralangan traektoriyaning eksponentlari va Lyapunov o'lchovi attraktor ostida o'zgarmasdir diffeomorfizm faza makonining[9]

The multiplikativ teskari Lyapunovning eng yirik eksponentlaridan biri ba'zan adabiyotda shunday nomlanadi Lyapunov vaqti va xarakteristikasini belgilaydi e- katlama vaqti. Xaotik orbitalar uchun Lyapunov vaqti cheklangan bo'ladi, doimiy orbitalar uchun esa cheksiz bo'ladi.

Raqamli hisoblash

Odatda Lyapunov eksponentlarini hisoblash, yuqorida ta'riflanganidek, analitik tarzda amalga oshirilmaydi va aksariyat hollarda raqamli usullarga murojaat qilish kerak. Xaotik traektoriyalarning eksponent divergentsiyasining birinchi namoyishini tashkil etgan dastlabki misol, tomonidan amalga oshirildi R. H. Miller 1964 yilda.[10] Hozirgi kunda eng ko'p ishlatiladigan raqamli protsedura matritsa, chegarani belgilashning bir necha cheklangan vaqt yaqinlashuvlarini o'rtacha hisobiga asoslangan .

Silliq dinamik tizim uchun Lyapunov spektrini hisoblashda eng ko'p ishlatiladigan va samarali sonli texnikalardan biri davriylikka bog'liq.Gram-Shmidt orthonormalizatsiya Lyapunov vektorlari maksimal kengayish yo'nalishi bo'yicha barcha vektorlarning noto'g'riligini oldini olish uchun.[11][12][13][14]

Cheklangan eksperimental ma'lumotlardan Lyapunov ko'rsatkichlarini hisoblash uchun turli usullar taklif qilingan. Biroq, ushbu usullarni qo'llashda juda ko'p qiyinchiliklar mavjud va bunday muammolarga ehtiyotkorlik bilan murojaat qilish kerak. Asosiy qiyinchilik shundaki, ma'lumotlar fazaviy makonni to'liq o'rgana olmaydi, aksincha ular ma'lum yo'nalishlar bo'yicha kengaytirilgan (agar mavjud bo'lsa) chekuvchi bilan cheklangan bo'lsa. Ma'lumotlar to'plamidagi ushbu ingichka yoki alohida yo'nalishlar ko'proq salbiy ko'rsatkichlar bilan bog'liq bo'lgan yo'nalishlardir. Attraktordan kichik siljishlar evolyutsiyasini modellashtirish uchun chiziqli bo'lmagan xaritalardan foydalanish Lyapunov spektrini tiklash qobiliyatini keskin yaxshilaganligi,[15][16] ma'lumotlar shovqin darajasi juda past bo'lsa. Ma'lumotlarning o'ziga xos xususiyati va uning salbiy ko'rsatkichlar bilan aloqasi ham o'rganildi.[17]

Mahalliy Lyapunov eksponenti

Lyapunov (global) eksponenti tizimning umumiy taxminiyligi uchun o'lchov berar ekan, ba'zida lokal prognozni nuqta atrofida baholash qiziq bo'ladi. x0 fazaviy fazoda. Bu orqali amalga oshirilishi mumkin o'zgacha qiymatlar ning Jacobian matritsa J 0(x0). Ushbu o'ziga xos qiymatlar Lyapunovning mahalliy eksponatlari deb ham ataladi.[18] (E'tiborli narsa: global eksponentlardan farqli o'laroq, bu mahalliy ko'rsatkichlar koordinatalarning chiziqli bo'lmagan o'zgarishi ostida o'zgarmas emas).

Shartli Lyapunov ko'rsatkichi

Ushbu atama odatda nisbatan qo'llaniladi betartiblikni sinxronlashtirish, unda birlashtirilgan ikkita tizim mavjud bo'lib, odatda bir yo'nalishda, haydovchi (yoki master) tizimi va javob (yoki qul) tizimi mavjud. Shartli ko'rsatkichlar - bu (xaotik) qo'zg'alish signalining manbai sifatida qabul qilingan haydovchi tizimi bilan javob berish tizimining ko'rsatkichlari. Sinxronizatsiya barcha shartli ko'rsatkichlar manfiy bo'lganda sodir bo'ladi.[19]

Shuningdek qarang

Adabiyotlar

  1. ^ Cencini, M .; va boshq. (2010). Jahon ilmiy (tahrir). Xaos Oddiy modellardan murakkab tizimlarga. ISBN  978-981-4277-65-5.
  2. ^ Temam, R. (1988). Mexanika va fizikadagi cheksiz o'lchovli dinamik tizimlar. Kembrij: Springer-Verlag.
  3. ^ a b N.V.Kuznetsov; G.A. Leonov (2005). Diskret tizimlar uchun birinchi taxmin bo'yicha barqarorlik to'g'risida (PDF). 2005 yil Fizika va boshqaruv bo'yicha xalqaro konferentsiya, PhysCon 2005 yil. Ishlar jildi 2005. 596-599 betlar. doi:10.1109 / PHYCON.2005.1514053. ISBN  978-0-7803-9235-9. S2CID  31746738.
  4. ^ a b v G.A. Leonov; N.V.Kuznetsov (2007). "Vaqt bo'yicha o'zgaruvchan chiziqlash va Perron effektlari" (PDF). Xalqaro bifurkatsiya va betartiblik jurnali. 17 (4): 1079–1107. Bibcode:2007IJBC ... 17.1079L. CiteSeerX  10.1.1.660.43. doi:10.1142 / S0218127407017732.
  5. ^ Kuznetsov, Nikolay; Reitmann, Volker (2020). Dinamik tizimlar uchun attraktor o'lchamlari taxminlari: nazariya va hisoblash. Cham: Springer.
  6. ^ Kaplan, J. va York, J. (1979). "Ko'p o'lchovli farq tenglamalarining xaotik harakati". Peitgen, H. O. & Walther, H. O. (tahr.). Funktsional differentsial tenglamalar va sobit nuqtalarni yaqinlashishi. Nyu-York: Springer. ISBN  978-3-540-09518-7.
  7. ^ Pesin, Y. B. (1977). "Lyapunovning o'ziga xos xususiyatlari va yumshoq ergodik nazariya". Rus matematikasi. So'rovnomalar. 32 (4): 55–114. Bibcode:1977RuMaS..32 ... 55P. doi:10.1070 / RM1977v032n04ABEH001639.
  8. ^ Kuznetsov, N.V. (2016). "Lyapunov o'lchovi va uni Leonov usuli bo'yicha baholash". Fizika xatlari. 380 (25–26): 2142–2149. arXiv:1602.05410. Bibcode:2016PhLA..380.2142K. doi:10.1016 / j.physleta.2016.04.036. S2CID  118467839.
  9. ^ Kuznetsov, N.V .; Alekseyeva, T.A .; Leonov, G.A. (2016). "Lyapunov eksponentlarining o'zgaruvchanligi va muntazam va tartibsiz chiziqli chiziqlar uchun Lyapunov o'lchovi". Lineer bo'lmagan dinamikalar. 85 (1): 195–201. arXiv:1410.2016. doi:10.1007 / s11071-016-2678-4. S2CID  119650438.
  10. ^ Miller, R. H. (1964). "Kichik yulduzlar dinamik tizimlarida qaytarilmaslik". Astrofizika jurnali. 140: 250. Bibcode:1964ApJ ... 140..250M. doi:10.1086/147911.
  11. ^ Benettin, G.; Galgani, L .; Giorgilli, A .; Strelcyn, J. M. (1980). "Yalang'och dinamik tizimlar va hamilton tizimlar uchun Lyapunovga xos xususiyatlar; ularning barchasini hisoblash usuli. 1-qism: Nazariya". Makkanika. 15: 9–20. doi:10.1007 / BF02128236. S2CID  123085922.
  12. ^ Benettin, G.; Galgani, L .; Giorgilli, A .; Strelcyn, J. M. (1980). "Lyapunov silliq dinamik tizimlar va gamilton tizimlar uchun xarakterli ko'rsatkichlar; ularning barchasini hisoblash usuli. 2-qism: Raqamli qo'llanma". Makkanika. 15: 21–30. doi:10.1007 / BF02128237. S2CID  117095512.
  13. ^ Shimada, I .; Nagashima, T. (1979). "Dissipativ dinamik tizimlarning ergodik muammosiga raqamli yondashuv". Nazariy fizikaning taraqqiyoti. 61 (6): 1605–1616. Bibcode:1979PhPh..61.1605S. doi:10.1143 / PTP.61.1605.
  14. ^ Ekman, J. -P .; Ruelle, D. (1985). "Xaos va g'alati attraksionlarning ergodik nazariyasi". Zamonaviy fizika sharhlari. 57 (3): 617–656. Bibcode:1985RvMP ... 57..617E. doi:10.1103 / RevModPhys.57.617. S2CID  18330392.
  15. ^ Bryant, P.; Braun, R .; Abarbanel, H. (1990). "Lyapunovning kuzatilgan vaqt seriyasidagi ko'rsatkichlari". Jismoniy tekshiruv xatlari. 65 (13): 1523–1526. Bibcode:1990PhRvL..65.1523B. doi:10.1103 / PhysRevLett.65.1523. PMID  10042292.
  16. ^ Braun, R .; Bryant, P.; Abarbanel, H. (1991). "Dinamik tizimning Lyapunov spektrini kuzatilgan vaqt qatoridan hisoblash". Jismoniy sharh A. 43 (6): 2787–2806. Bibcode:1991PhRvA..43.2787B. doi:10.1103 / PhysRevA.43.2787. PMID  9905344.
  17. ^ Bryant, P. H. (1993). "Ajablanarlisi uchun o'ziga xoslik o'lchovlari". Fizika xatlari. 179 (3): 186–190. Bibcode:1993 yil PHLA..179..186B. doi:10.1016 / 0375-9601 (93) 91136-S.
  18. ^ Abarbanel, XDI; Braun, R .; Kennel, M.B. (1992). "Mahalliy Lyapunov eksponentlari kuzatilgan ma'lumotlardan hisoblangan". Lineer bo'lmagan fan jurnali. 2 (3): 343–365. Bibcode:1992JNS ..... 2..343A. doi:10.1007 / BF01208929. S2CID  122542761.
  19. ^ Qarang, masalan, Pekora, L. M .; Kerrol, T. L .; Jonson, G. A .; Mar, D. J .; Heagy, J. F. (1997). "Xaotik tizimlarda sinxronizatsiya asoslari, tushunchalari va ilovalari". Xaos: fanlararo jurnal. 7 (4): 520–543. Bibcode:1997 yil Xaos ... 7..520P. doi:10.1063/1.166278. PMID  12779679.

Qo'shimcha o'qish

  • M.-F. Danca va N.V.Kuznetsov (2018). "Fraksiyonel-tartibli tizimlarning Lyapunov eksponentlari uchun Matlab kodi". Xalqaro bifurkatsiya va betartiblik jurnali. 25 (5): san'at. num. 1850067. doi:10.1142 / S0218127418500670.

Dasturiy ta'minot

  • [1] R. Hegger, H. Kantz va T. Shrayber, Lineer bo'lmagan vaqt seriyasini tahlil qilish, TISEAN 3.0.1 (2007 yil mart).
  • [2] Scientio's ChaosKit mahsuloti boshqa xaotik tadbirlar qatorida Lyapunov eksponentlarini hisoblab chiqadi. Kirish onlayn ravishda veb-xizmat va Silverlight demo orqali taqdim etiladi.
  • [3][doimiy o'lik havola ] Doktor Ronald Djo Rekordning matematik rekreatsiya dasturiy laboratoriyasida majburiy logistika xaritasi va birlik oralig'idagi boshqa xaritalarning Lyapunov ko'rsatkichlarini grafik ravishda o'rganish uchun X11 grafik mijozi lyap mavjud. The mazmuni va qo'llanma sahifalari[doimiy o'lik havola ] mathrec dasturiy ta'minot laboratoriyasi ham mavjud.
  • [4] Ushbu sahifadagi dastur eksponentlarning to'liq spektrini samarali va aniq hisoblash uchun maxsus ishlab chiqilgan. Bunga harakat tenglamalari ma'lum bo'lgan holatlar uchun LyapOde va eksperimental vaqt qatorlari bilan bog'liq holatlar uchun Lyap kiradi. LyapOde, "C" da yozilgan manba kodini o'z ichiga oladi, shuningdek, birlashtirilgan tizimlar uchun shartli Lyapunov ko'rsatkichlarini hisoblashi mumkin. Bu foydalanuvchiga o'z model tenglamalarini taqdim etish yoki unga kiritilganlardan birini ishlatishga imkon berish uchun mo'ljallangan. Fortranda yozilgan manba kodini o'z ichiga olgan Lyap, o'zgarmaydiganlar, parametrlar soniga xos cheklovlar yo'q, shuningdek Lyapunov yo'nalishi vektorlarini hisoblab chiqishi va o'ziga jalb etuvchi o'ziga xosligini tavsiflashi mumkin, bu ko'proq hisoblashda qiyinchiliklarning asosiy sababi hisoblanadi. vaqt qatorlari ma'lumotlarining salbiy ko'rsatkichlari. Ikkala holatda ham keng hujjatlar va namunaviy kirish fayllari mavjud. Dastur Windows, Mac yoki Linux / Unix tizimlarida ishlash uchun tuzilishi mumkin. Dastur matnli oynada ishlaydi va grafik imkoniyatlariga ega emas, lekin Excel kabi dastur yordamida osongina tuzilishi mumkin bo'lgan chiqish fayllarini yaratishi mumkin.

Tashqi havolalar