I systemer som involverer fraksjonell derivasjonsdemping, kan den matematiske beskrivelsen bli betydelig mer kompleks enn i tradisjonelle modeller. Dette skyldes at fraksjonelle derivasjoner representerer en generalisering av de vanlige derivasjonene og innfører et ikke-kommuntativt aspekt i de dynamiske ligningene. En slik modell kan gi innsikt i mer realistiske fysiske systemer, der visse effekter som hukommelseseffekter eller langsom variabilitet i responsen er viktige. Dette er spesielt relevant i analysen av quasi-integrable Hamiltoniansystemer under støyforstyrrelser.

I et klassisk Hamiltoniansystem beskrives systemets tilstand med posisjon (Q) og momentum (P) i et fase rom, og systemet kan ofte forenkles til et sett av bevarte energier og harmoniske bevegelser. Når et system inneholder fraksjonelle derivasjoner som demperende krefter, forandres dette bildet vesentlig. Fraksjonell derivasjon innebærer at responsen til systemet ikke nødvendigvis er lokal i tid, men kan være hukommelsesbasert. Dette fører til en endring i både systemets dynamikk og de tilhørende energifunksjonene, som nå blir et resultat av fraksjonelle krefter.

For et system der fraksjonelle derivasjoner er til stede, som i tilfelle av et quasi-integrabelt Hamiltoniansystem, kan man anvende metoder som den generaliserte harmoniske balanse-teknikken for å beskrive de dempende kreftene. I et slikt system kan de fraksjonelle dempende kreftene uttrykkes som:

Dα1(X1X2)=C11(A1)X1˙C12(A2)X2˙+K11(A1)X1K12(A2)X2D_{\alpha 1}(X_1 - X_2) = C_{11}(A_1) \dot{X_1} - C_{12}(A_2) \dot{X_2} + K_{11}(A_1) X_1 - K_{12}(A_2) X_2
Dα2(X2X1)=C22(A2)X2˙C21(A1)X1˙+K22(A2)X2K21(A1)X1D_{\alpha 2}(X_2 - X_1) = C_{22}(A_2) \dot{X_2} - C_{21}(A_1) \dot{X_1} + K_{22}(A_2) X_2 - K_{21}(A_1) X_1

Disse kreftene kan videre bli implementert i de originale systemligningene, og systemet kan omformes til et quasi-Hamiltoniansystem hvor bevegelsen av de to koordinatene og momentene kan modelleres. De resulterende bevegelseslikningene er i stor grad styrt av støyforstyrrelser, som representeres ved hvit støy, og vi får et stasjonært sannsynlighetsfordelingssystem som kan analyseres ved hjelp av stokastiske gjennomsnittsmetoder.

Når vi ser på resonante systemer, må vi vurdere tilfeller hvor de to frekvensene, ω1\omega_1 og ω2\omega_2, er relatert i en intern resonansforhold. I slike tilfeller kan systemet beskrives som en tre-dimensjonal Markov-prosess, som kan analyseres ytterligere ved å bruke stokastiske gjennomsnittsmetoder for resonante Hamiltoniansystemer under Gaussisk hvit støy. I denne tilnærmingen kan vi derivere de gjennomsnittlige Itô stokastiske differensialligningene og deres tilhørende reduserte Fokker-Planck ligning.

For ikke-resonante tilfeller, hvor ω1ω2\omega_1 \neq \omega_2, kan systemet behandles som et to-dimensjonalt Markov-diffusjonsprosesser, og den stasjonære sannsynlighetsfordelingen kan beregnes numerisk gjennom de reduserte Fokker-Planck-ligningene. For systemer med spesifikke parameterverdier kan man bruke nøyaktige analytiske løsninger for den stasjonære sannsynlighetsfordelingen, som kan sammenlignes med resultater fra Monte Carlo-simuleringer for å validere modellen.

En viktig del av denne analysen er å forstå hvordan systemets respons på støy og fraksjonell demping påvirker dets langsiktige statistiske egenskaper. Det kan være nyttig å utføre numeriske simuleringer for å visualisere de stasjonære fordelingene for systemets koordinater og momenter. Ved å analysere disse distribusjonene kan man få innsikt i hvordan systemet oppfører seg under forskjellige støybetingelser og hvilke faktorer som styrer systemets stasjonære tilstand.

I tillegg til de rent matematiske betraktningene, er det viktig å merke seg at stokastiske metoder som de som er beskrevet her, gir en kraftig verktøykasse for å studere komplekse fysiske systemer som er vanskelig å behandle med tradisjonelle analytiske teknikker. Metodene kan også utvides til systemer med flere grad av frihet, og det er viktig å vurdere hvordan disse metodene kan tilpasses for å håndtere ikke-lineære interaksjoner, eller systemer som utsettes for forskjellige former for støy.

Hvordan kan stokastisk gjennomsnittsmetode anvendes for kvasi-integrerbare Hamiltoniske systemer med stasjonær bredbåndsstøy?

Stokastisk gjennomsnittsmetode er et kraftig verktøy for å forenkle analysen av dynamikken i komplekse systemer, spesielt når disse systemene er utsatt for støy. For kvasi-integrerbare Hamiltoniske systemer som påvirkes av stasjonær bredbåndsstøy, kan metoden brukes til å finne en analytisk løsning for sannsynlighetstettheten til systemets tilstand. Denne tilnærmingen er spesielt nyttig når systemet har flere frihetsgrader og støyen kan anses som et stasjonært bredbåndsprosesser.

I et kvasi-integrerbart Hamiltonisk system, kan systemets tilstand beskrives av et sett med variable som representerer aksjoner og faser. La oss anta at den totale sannsynlighetstettheten, p(a,ψ)p(a, \psi), for systemet er gitt som en eksponentiell funksjon av et potensial λ(a,ψ)\lambda(a, \psi), som representerer sannsynligheten for at systemet befinner seg i en gitt tilstand. Denne potensialfunksjonen kan ekspanderes som en trukket Fourier-serie med hensyn til de ulike fasevariablene ψ1,ψ2,ψ3\psi_1, \psi_2, \psi_3. Denne Fourier-ekspansjonen gir en forenklet måte å håndtere de komplekse periodiske funksjonene som beskriver systemets dynamikk under støypåvirkning.

Når den potensielle funksjonen λ(a,ψ)\lambda(a, \psi) er kjent, kan den substitueres inn i den stokastiske Fokker-Planck-ligningen for å finne den stasjonære sannsynlighetstettheten for aksjonsvariablene I1,I2,I3,I4I_1, I_2, I_3, I_4. Løsningen av denne ligningen gir oss informasjon om hvordan sannsynligheten for forskjellige systemtilstander fordeler seg over tid. Denne stasjonære løsningen kan estimeres både analytisk og ved hjelp av Monte Carlo-simuleringer, og det er viktig at de analytiske resultatene stemmer overens med numeriske simuleringer for å validere metoden.

Når det gjelder støy, er det viktig å merke seg at systemene kan være utsatt for ulike typer støy, som for eksempel bredbånds støy eller fraksjonell Gaussisk støy (fGn). Den fraksjonelle Gaussiske støyen har en langsom frekvensavhengighet, og etter et visst punkt kan denne støyen betraktes som en stasjonær bredbåndsstøy. Hvis de naturlige frekvensene til et kvasi-integrerbart system ligger i dette området, kan fGn behandles som en bredbåndsprosess, og man kan bruke metodene som er utviklet for stasjonær bredbåndsstøy for å forenkle beregningene. Dette reduserer dimensjonaliteten til problemet og gjør det mer håndterbart, både analytisk og numerisk.

I slike systemer, hvor støyen ikke er et Markov-prosess, kan de stokastiske gjennomsnittsmetodene bidra til å forenkle dynamikken betydelig. I stedet for å løse de komplekse differensialligningene direkte, kan man anvende en transformasjon som forvandler systemets tilstand fra Qi,PiQ_i, P_i til aksjons- og fasevariablene Ai,ϕiA_i, \phi_i. Denne transformasjonen forenkler systemet ved å bruke tilfeldige prosesser for AiA_i og ϕi\phi_i, som gir en bedre forståelse av systemets dynamikk under støy.

Etter transformasjonen kan systemet beskrives ved hjelp av Itô-ligninger, som gir en enklere måte å beskrive systemets utvikling i tid. For eksempel, i det ikke-resonante tilfellet, vil aksjonsvariablene A=[A1,A2,,An]TA = [A_1, A_2, \dots, A_n]^T konvergere til et svakt Markov-diffusjonsprosess. De resulterende Itô-ligningene gir en forenklet beskrivelse av systemet, der drift og diffusjonskoeffisienter er relatert til systemets dynamikk og støyforhold.

Stokastiske Fokker-Planck-ligninger (FPK) kan også utledes for systemet i denne forenklede formen, og løsningen av disse gir den stasjonære sannsynlighetstettheten for systemets tilstand i form av aksjonsvariablene. Dette gir et nyttig verktøy for å forstå hvordan systemet vil oppføre seg over tid under påvirkning av støy.

For systemer med intern resonans vil dimensjonen og formen på de gjennomsnittede Itô-ligningene variere, og dette må tas i betraktning når man bruker den stokastiske gjennomsnittsmetoden. I slike tilfeller kan det være nødvendig å bruke spesifikke metoder for å håndtere resonansbetingelsene, og det kan være utfordrende å finne en løsning som er nøyaktig uten numeriske simuleringer.

I tillegg til å redusere dimensjonaliteten til problemet, gir stokastiske gjennomsnittsmetoder en mer effektiv måte å håndtere komplekse systemer som involverer både støy og ikke-lineære dynamikker. Når støyen kan behandles som bredbånds, som i tilfelle av fraksjonell Gaussisk støy, kan metoden gi gode resultater på en mye raskere og mer effektiv måte enn ved å bruke fullstendige simuleringer av systemet.

Det er viktig å merke seg at metodene beskrevet her er spesielt nyttige når systemet er kvasi-integrerbart og når støyen er av en type som kan beskrives som bredbånds. For andre typer støy eller for systemer som ikke er kvasi-integrerbare, kan det være nødvendig med andre tilnærminger eller metoder for å oppnå en nøyaktig løsning.

Hvordan beregnes maksimal Lyapunov-eksponent for stasjonære stokastiske systemer med flere frihetsgrader?

I analysen av lineære systemer med flere frihetsgrader (n-DOF) som er utsatt for stasjonær stokastisk parametrisk eksitasjon, spiller maksimal Lyapunov-eksponent en sentral rolle for å bestemme systemets asymptotiske stabilitet. For et ikke-gyroskopisk system, beskrevet ved en Hamiltonsk struktur, kan den stokastiske bevegelseslikningen formuleres som en Itô-stokastisk differensialligning der parametrene påvirkes av en bredbåndet stasjonær prosess med null forventningsverdi og kjent korrelasjonsfunksjon.

Ved å anvende stokastisk averaging i tilfeller uten intern resonans, reduseres det komplekse systemet til en Itô-likning for amplitudene til hvert enkelt oscillerende ledd. Deretter transformeres amplitudene til energilikninger, hvor drift- og diffusjonskoeffisientene kan utledes eksakt gjennom spectral densiteter av eksitasjonen, korrelert med systemets naturlige frekvenser. Disse koeffisientene må tilfredsstille visse matematiske betingelser (som ergodisitet og markov-egenskaper) for at maksimal Lyapunov-eksponent kan beregnes ved hjelp av integraler over den stasjonære fordelingen av systemets tilstandsvariabler.

For en todimensjonal systemvariant (2-DOF) kan uttrykkene for drift og diffusjon forenkles til å omfatte kombinasjoner av spectral densiteter og systemkonstanter, noe som igjen tillater definisjonen av en norm og en logaritmisk transformasjon av energien for å framstille systemets dynamikk i en form som egner seg for analyse av asymptotisk stabilitet. Itô-differensialligninger for den normaliserte andelen av energien, som begrenses mellom 0 og 1, kan analyseres ved hjelp av grensebetingelser som klassifiserer begge endepunktene som «entrance»-grenser. Dette sikrer at en stasjonær sannsynlighetsfordeling eksisterer og kan finnes ved løsning av tilhørende Fokker-Planck-kvasi-likninger.

Denne distribusjonen tillater deretter evaluering av maksimal Lyapunov-eksponent som et forventningsverdi-integral over driftfunksjonen vektet med den stasjonære tettheten, noe som bekrefter sammenfallet med resultater oppnådd ved klassiske metoder for stokastisk averaging.

En praktisk anvendelse av denne metoden er innen stabilitetsvurdering av gyroskopiske systemer, som gyroskoppendler som benyttes i navigasjonssystemer. Disse pendler må opprettholde retningsstabilitet under påvirkning av stokastiske vertikale vibrasjoner fra bærerplattformen. Her kan lignende Hamiltonske metoder anvendes etter en Legendre-transformasjon, og eksitasjonen modelleres som en bredbåndet stasjonær prosess. Systemets stabilitet analyseres ved hjelp av stokastisk averaging av den Lagrangianske bevegelseslikningen, og transformasjonen til Hamiltonsk form muliggjør en presis beregning av stabilitetskriterier gjennom den samme formalismen for maksimal Lyapunov-eksponent.

Den stokastiske differensiallikningen som styrer gyroskoppendelens dynamikk, inkludert gyroskopiske krefter og små dempingsparametere, kan derfor evalueres for stabilitet ved å beregne eksponentens tegn. En positiv maksimal Lyapunov-eksponent indikerer asymptotisk ustabilitet, mens en negativ verdi tilsvarer stabilitet med sannsynlighet 1. Dette gir et kraftfullt verktøy for å vurdere effekten av stokastiske forstyrrelser i mekaniske systemer med flere frihetsgrader.

Det er viktig å forstå at denne metodikken forutsetter ergodisitet og markovegenskaper i de stokastiske prosessene som påvirker systemet. Videre er antakelsen om svake stokastiske forstyrrelser og linearisering rundt stabile løsninger sentral for gyldigheten av analysen. I praktisk bruk må derfor systemparametrene og eksitasjonens statistiske egenskaper evalueres nøye for å sikre at betingelsene for anvendelse av disse metodene er oppfylt.

Endvidere er den teoretiske basen for denne analysen tett knyttet til Hamiltonsk dynamikk og stokastisk kalkulus, og for en dypere forståelse må leseren være fortrolig med Itô-prosesser, Fokker-Planck-ligninger, og spektral analyse av stokastiske prosesser. Denne helhetlige tilnærmingen gir imidlertid et robust rammeverk for å forutsi og kontrollere stabilitet i avanserte tekniske systemer under tilfeldige påvirkninger.