Stokastisk genomsnittning är en kraftfull metod som används för att förenkla och analysera komplexa system där både långsamma och snabba förändringar förekommer. Metoden har sin grund i att göra systemets dynamik mer hanterbar genom att approximera snabba variabler och reducera systemets dimensioner. Stokastisk genomsnittning används för att behandla system där externa excitationer, som ofta är bredbandiga, påverkar systemet. Det här tillvägagångssättet gör det möjligt att få fram en förenklad modell som kan ge insikter i systemets långsiktiga beteende.

I allmänhet kan stokastisk genomsnittning genomföras i två huvudsakliga procedurer. Den första är att approximera excitationerna som Gaussiska vita brusprocesser och att modellera systemets respons som en Markov-diffusionsprocess. Den andra proceduren involverar att genomföra tidsgenomsnittning för att eliminera de snabbare variablerna och reducera systemets dimensioner, vilket resulterar i en förenklad modell som kallas den "utjämnade versionen". Denna typ av förenkling gör det möjligt att beräkna drift- och diffusionskoefficienter som inte längre explicit beror på tid.

I många fall, efter tidsgenomsnittningen, kan det visa sig att endast en del av systemets respons (eller en del av responsvektorn) består av en diffusionsprocess. Detta innebär en ytterligare reduktion av systemets dimensioner. Ett exempel på detta är enkeldimensionella system där amplituden och fasen behandlas som långsamt varierande storheter, trots att de initialt inte är långsamt varierande. Detta gör det möjligt att analysera dessa system med hjälp av de förenklade, stokastiskt genomsnittsberäknade ekvationerna.

För ett specifikt SDOF-system med linjär styvhet och svag icke-linjär dämpning kan vi skriva systemets rörelseekvation som:

mX¨+ϵh(X,X˙)+ω02X=ϵ1/2gl(X,X˙)ξl(t)\sum m \ddot{X} + \epsilon h(X, \dot{X}) + \omega_0^2 X = \epsilon^{1/2} \sum g_l(X, \dot{X}) \xi_l(t)

Där XX och X˙\dot{X} representerar systemets förflyttningar och hastigheter, och ξl(t)\xi_l(t) är externa excitationer. Här kan vi introducera en transformation där systemets förflyttning XX skrivs som en amplitudprocess A(t)A(t) multiplicerad med en trigonometrisk funktion:

X=A(t)cos(θ),X˙=A(t)ω0sin(θ)X = A(t) \cos(\theta), \quad \dot{X} = -A(t)\omega_0 \sin(\theta)

Där θ=ω0t+φ(t)\theta = \omega_0 t + \varphi(t) representerar den fasvinkel som beror på den långsamt varierande fasfunktionen φ(t)\varphi(t). Efter att ha differentierat och substituerat dessa uttryck i rörelseekvationen, kan man hitta uttryck för de långsamt varierande variablerna A(t)A(t) och φ(t)\varphi(t). Eftersom dessa är långsamt varierande, innebär det att systemet kan betraktas som en Markov-diffusionsprocess, där den stokastiska dynamiken endast beror på amplituden A(t)A(t).

Vidare, efter att ha genomfört tidsgenomsnittning på systemet, kan vi uttrycka drift- och diffusionskoefficienterna m(A)m(A) och σ(A)\sigma(A), vilket resulterar i en förenklad differentialekvation:

dA=m(A)dt+σ(A)dB(t)dA = m(A)dt + \sigma(A)dB(t)

Denna differentialekvation är en enklare version av systemets dynamik som gör det lättare att analysera och förstå systemets långsiktiga beteende. Den stationära sannolikhetsdensiteten för amplitudprocessen A(t)A(t) kan sedan härledas och skrivs som en exponentiell fördelning, vilket gör att vi kan studera systemets statistiska egenskaper över tid.

Det är också viktigt att beakta att även om amplitudprocessen A(t)A(t) är en viktig egenskap hos systemets rörelse, kan den gemensamma stationära sannolikhetsdensiteten, såväl som de marginala fördelningarna för X(t)X(t) och X˙(t)\dot{X}(t), vara nödvändiga för en mer detaljerad analys. Genom att beräkna dessa fördelningar kan man få en djupare förståelse för systemets beteende i olika tillstånd och över tid.

För system med icke-linjär dämpning och excitationer som är bredbandiga, där rörelsen inte längre kan beskrivas som en enkel harmonisk oscillator, kan vi istället analysera systemets energibehov. Detta görs genom att använda energinivåerna hos systemet, vilket ger en mer exakt förståelse för systemets dynamik, särskilt i de fall där amplituden beror på den totala energin i systemet. Genom att studera energin som en funktion av amplituden, kan vi hitta samband mellan perioder och energi, vilket ytterligare förenklar analysen av sådana komplexa system.

Stokastisk genomsnittning är därför inte bara ett matematiskt verktyg för att förenkla modeller utan också ett sätt att få insikter i hur ett komplext system beter sig under påverkan av externa, osäkra excitationer. Det gör det möjligt att beskriva systemets dynamik i enklare, mer hanterbara termer och att identifiera nyckelparametrar som styr dess beteende.

Hur kan stokastisk genomsnittsmetod tillämpas på kvasi-Hamiltonianska system under Gaussisk vitbrus-excitation?

Stokastiska genomsnittsmetoder är centrala för att studera dynamiken hos system utsatta för slumpmässiga excitationer, särskilt när dessa system är icke-linjära eller har flera frihetsgrader. Ett av de viktigaste verktygen för att analysera sådana system är att använda en förenklad modell som fångar de grundläggande fenomenen utan att behöva lösa den fullständiga systemdynamiken. I denna kontext är kvasi-Hamiltonianska system en särskild kategori som spelar en viktig roll i många tillämpningar där svaga excitationer och dämpning förekommer.

För att förstå tillämpningen av stokastiska genomsnittsmetoder på kvasi-Hamiltonianska system måste vi först definiera systemets dynamik. Ett n–DOF system (system med n frihetsgrader) som är stochastiskt exciterat och dämpat kan beskrivas genom Hamiltons ekvationer, där Q representerar de generaliserade förskjutningarna och P de generaliserade impulserna. Hamiltonfunktionen H(Q, P) beskriver systemets totala energi, inklusive både kinetisk och potentiell energi, medan dämpningskoefficienterna (cij) och excitationsamplituderna (gil) styr systemets respons på externa krafter, såsom Gaussiskt vitbrus.

När systemet är svagt exciterat och dämpat kan vi modellera det som ett kvasi-Hamiltonianskt system. Detta innebär att de systemiska beteendena domineras av långsam förändring av den totala energin H(t), medan de andra variablerna som beskriver rörelsen i systemet ändras snabbt. I detta fall kan vi använda stokastiska differentialekvationer för att få en realistisk beskrivning av systemets dynamik.

För att lösa sådana system effektivt kan vi tillämpa Itôs stokastiska differentialekvationer, vilket innebär att vi inför korrigeringstermer, kända som Wong-Zakai-termer. Dessa korrigeringar hanterar effekten av den stochastiska bruset på återställande krafter och dämpande krafter. Genom att kombinera dessa termer med systemets grundläggande fysik kan vi få fram en uppdaterad version av Hamiltonfunktionen som är mer exakt för de stokastiska system vi studerar.

En viktig aspekt är att även om systemet kan beskrivas av stokastiska differentialekvationer, är den totala energin H(t) ett långsamt varierande process, medan de övriga variablerna ändras snabbt. Genom att använda Khasminskiis stokastiska genomsnittsmetod, som säger att för små parametrar (ε → 0) konvergerar systemets energi till en Markov-diffusionsprocess, kan vi få en förenklad beskrivning av systemet. Denna process kan representeras av en genomsnitts- Itô-ekvation, vilket gör det möjligt att få fram en effektiverad modell för systemets långsiktiga beteende.

För att få fram drift- och diffusionskoefficienterna för denna genomsnittsprocess genomförs en tidsgenomsnittsoperation över de stokastiska differentialekvationerna. Eftersom alla variabler i dessa ekvationer är stokastiska processer, kan det vara utmanande att utföra en fullständig tidsgenomsnittsoperation. Men genom att utnyttja ergodicitetsegenskaperna hos ett n-DOF icke-integrerbart Hamiltonianskt system, kan vi ersätta tidsgenomsnittsoperationen med en spatial genomsnittsoperation över isoenergetiska ytor. Denna metod gör det möjligt att förenkla den stokastiska analysen och ge en realistisk beskrivning av systemets långsiktiga beteende.

Vid tillämpning av denna metod är det avgörande att förstå att även om systemets totala energi är långsamt varierande, påverkar de stokastiska krafterna systemets dynamik på en mycket finare skala. Därför blir det viktigt att förstå hur systemets olika komponenter samverkar för att ge en helhetsbild av responsen på externa excitationer. Den exakta formen av de stokastiska differentialekvationerna, tillsammans med den specifika strukturen hos de olika koefficienterna (dämpning, excitation etc.), avgör hur effektivt och exakt vi kan beskriva systemets dynamik under slumpmässiga förhållanden.

För att förstå den praktiska tillämpningen av dessa metoder är det också viktigt att överväga systemets resonansbeteenden. I de flesta praktiska tillämpningar, särskilt de som involverar starkt icke-linjära system, är resonansfenomen en central aspekt som måste beaktas. Vid resonans kan systemet uppvisa starka och potentiellt instabila svar, vilket kräver noggrant val av parametrar för att undvika oönskade effekter. I dessa fall blir valet av metod för stokastisk genomsnittsoperation ännu mer kritiskt, eftersom resonans kan påverka hela systemets stabilitet.

Endtext

Hur kan stokastisk medelvärdesbildning tillämpas på kvasi-Hamiltonska system med stokastisk excitation?

Stokastisk medelvärdesbildning för kvasi-Hamiltonska system, särskilt när dessa är exciterade av stokastiska processer med hopp och diffusion, erbjuder en kraftfull metod för att approximera systemets dynamik och sannolikhetsfördelningar i stationärt tillstånd. Genom att utgå från stokastiska differentialekvationer (SDE) och använda transformeringar, som Di Paola och Falsone-regeln, kan man formulera en SDE för Hamiltonian i systemet, vilken sedan kan förenklas och avrundas till en trunkiert genomsnittlig SDE samt en motsvarande Fokker–Planck–Kolmogorov (FPK) ekvation.

Vid denna medelvärdesbildning integreras systemets koefficienter över nivåytor definierade av Hamiltonian, där variabler som generaliserade positioner och rörelsemängder byts ut enligt definierade transformationer. Dessa integraler, ofta av hög dimension, kräver noggranna numeriska metoder, men leder till approximativa uttryck för stationära sannolikhetsdensiteter (PDF) av Hamiltonian och andra systemvariabler. Genom att lösa den genomsnittliga FPK-ekvationen med hjälp av perturbationsmetoder erhålls lösningar som kan jämföras med numeriska simuleringar, som Monte Carlo, vilket visar att metoden ofta är mer precis än simplare Gaussiska approximationer.

För kvasi-integrerbara Hamiltonska system, där systemet är nära ett fullständigt integrerbart system med tillhörande actions- och vinkelvariabler, kan stokastiska differentialekvationssystem (SIDEs) för actions och vinklar härledas via en stokastisk kedjeregel för hopp-diffusionsprocesser. Genom att expandera förändringarna i actions- och vinkelvariabler i Taylorserier kan dessa system förenklas till genomsnittliga SIDEs och FPK-ekvationer.

Det finns en avgörande skillnad i behandlingen beroende på resonansförhållanden mellan systemets frekvenser. Om systemets frekvenser inte uppfyller svaga interna resonansvillkor, dvs. om det inte finns linjära relationer mellan frekvenserna med små multiplikatorer av ordning ε², så konvergerar action-variablerna till en Markovprocess som kan beskrivas med genomsnittade SIDEs. Eftersom det icke-resonanta systemet är ergodiskt på den n-dimensionella tori som definieras av vinklarna, kan tidsmedelvärden ersättas med rumsmedelvärden i de genomsnittade ekvationerna. Detta förenklar beräkningar och ger en konkret metod för att beskriva dynamiken hos det stokastiskt exciterade Hamiltonska systemet på lång tidsskala.

I dessa sammanhang är det viktigt att förstå att stokastisk medelvärdesbildning inte bara reducerar komplexiteten i modellerna utan också gör det möjligt att erhålla meningsfulla approximativa fördelningar för systemets variabler i stationärt tillstånd. Resultaten kan tillämpas i mekaniska system, vibrationsproblem och andra fysiska modeller där stokastiska krafter påverkar systemets utveckling.

Det är också avgörande att beakta att den approximativa naturen i dessa metoder förutsätter små störningsparametrar (ε) och att högre ordningens termer ofta kan försummas. Därmed är noggrannheten beroende av systemets avstånd från ideal integrerbarhet och styrkan i de stokastiska excitationerna. I praktiken ger metoderna dock ett robust ramverk för att analysera komplexa stokastiskt påverkade Hamiltonska system där direkta numeriska simuleringar är resurskrävande eller svåra att tolka.

En annan viktig aspekt är att metoden bygger på antaganden om ergodicitet och resonansförhållanden, vilket kräver noggrann bedömning av systemets parametrar och frekvenser. Bristande uppfyllande av dessa antaganden kan leda till att approximativa modeller avviker från verklig dynamik. Därför bör läsaren ha en förståelse för när dessa metoder är applicerbara och när alternativa analyser krävs.

Slutligen bör läggas märke till att de stokastiska processerna som beskriver excitationerna ofta är komplexa med både kontinuerliga diffusionskomponenter och diskreta hopp, vilket speglar många verkliga system med plötsliga förändringar eller impulsiva krafter. Att behärska transformeringar och medelvärdesbildningsmetoder för sådana processer är avgörande för att korrekt beskriva och förutsäga dynamiken hos dessa avancerade system.