Når vi vurderer en stokastisk prosess X(t), kan denne reduseres til en ren harmonisk prosess med en tilfeldig innledende fase. Dette betyr at prosessen kan beskrives som et system hvor frekvensen er fast, men fasen ved startpunktet er tilfeldig fordelt. Når variansen σ øker, blir båndbredden til prosessen bredere, noe som indikerer en økende grad av tilfeldighet i systemets oppførsel.

Som et spesialtilfelle, når ω₀ = 0, kan korrelasjonsfunksjonen og spektral tettheten uttrykkes ved formlene:

RXX(τ)=1A2exp(σ2τ2)R_{XX}(\tau) = \frac{1}{A^2} \exp\left(-\frac{\sigma^2|\tau|}{2}\right)
SXX(ω)=σ24π(ω2+11)S_{XX}(\omega) = \frac{\sigma^2}{4 \pi \left(\omega^2 + \frac{1}{1}\right)}

Disse uttrykkene viser at når prosessen nærmer seg et lavpass-filter, får vi en korrelasjon som avtar eksponentielt med tiden, og spektral tettheten avtar med frekvensen. Dette er typisk for lavpassprosesser, som ofte brukes til å modellere systemer hvor høyfrekvente signaler dempes, for eksempel i mekaniske systemer som er utsatt for støy.

En annen interessant egenskap ved denne typen stokastiske prosesser er at sannsynlighetsfordelingen (PDF) for den randomiserte harmoniske prosessen X(t) ikke er typisk som de som genereres av lineære eller ikke-lineære filtre. Den generelle formen på PDFen for X(t) er gitt ved:

pX(x)=1πA21(xA)2p_X(x) = \frac{1}{\pi A^2} \sqrt{1 - \left(\frac{x}{A}\right)^2}

Her ser vi at sannsynlighetsfordelingen er mer konsentrert rundt de to grensene A-A og AA, noe som indikerer at prosessen er begrenset. Denne fordelingen avhenger utelukkende av amplitude AA, som må bestemmes basert på den fysiske grensen til fenomenet som modelleres. Parameterne ω₀ og σ påvirker ikke selve sannsynlighetsfordelingen, men de kan justeres for å matche spektral tetthet for prosessen som ønskes modellert.

En av de viktigste fordelene ved å bruke en randomisert harmonisk prosess som denne er at den gir en realistisk fremstilling av systemer som har et naturlig begrenset område, noe som ikke alltid er tilfelle i mer ideelle matematiske modeller. Ved å justere parameterne ω₀ og σ, kan man tilpasse spektral tetthet til den ønskede verdien, og på den måten oppnå en bedre representasjon av det fysiske systemet.

En annen viktig fordel er at modellen for en randomisert harmonisk prosess ikke kun er teoretisk interessant, men kan brukes i praktiske applikasjoner. For eksempel, i ingeniørfag og fysikk, kan slike prosesser brukes til å simulere og analysere systemer som er utsatt for støy og tilfeldige variasjoner, som kan være viktig i alt fra seismiske målinger til maskinlæring og prediksjon.

En avgjørende del ved å forstå denne prosessen er å erkjenne at spektral tetthet og sannsynlighetsfordeling er nært knyttet til de fysiske grensene av systemet som modelleres. Mens parametrene ω₀ og σ har liten effekt på sannsynlighetsfordelingen, er det nødvendig å justere disse for å tilpasse spektral tettheten. Dette gjør modellen fleksibel og egnet for mange typer praktiske anvendelser hvor man trenger å modellere prosesser med varierende grad av tilfeldighet og båndbredde.

Denne typen stokastiske prosesser er spesielt viktig i systemer som ikke bare er utsatt for tilfeldig støy, men også systemer der det er komplekse interaksjoner mellom ulike typer krefter, for eksempel viskoelastiske krefter eller krefter som har fraksjonelle deriverte. Dette gjør at forståelsen og bruken av slike prosesser ikke bare har applikasjoner i ren matematikk, men også i mange praktiske felt, som maskinteknikk, biomedisinsk ingeniørfag, og andre tekniske disipliner hvor kompleksitet og usikkerhet er til stede i modellene.

Hvordan Stasjonære Løsninger Oppstår i Quasi-Hamiltonske Systemer og Deres Stokastiske Besvarelser

Stasjonære løsninger i quasi-Hamiltonske systemer kan uttrykkes som p(I1,I2,h3)=Cexp[λ(I1,I2,h3)]p(I_1, I_2, h_3) = C \exp[-\lambda(I_1, I_2, h_3)], der CC er en konstant og λ(I1,I2,h3)\lambda(I_1, I_2, h_3) er en funksjon som tilfredsstiller spesifikke betingelser gitt i systemets beskrivelse. Denne løsningen, som stammer fra en reduksjon av den gjennomsnittlige FPK-ligningen, innebærer at de stasjonære sannsynlighetsfordelingene for integraler som I1I_1, I2I_2 og h3h_3 kan finnes ved å løse de partielle differensialligningene som er derivert fra systemets dynamikk.

Når man substituerer inn p(I1,I2,h3)p(I_1, I_2, h_3) i den gjennomsnittlige FPK-ligningen og setter den første tidsderivert av pp lik null, kan man bekrive hvordan den stasjonære løsningen utvikler seg under spesifikke forutsetninger. Spesielt når de relevante parametrene i systemet som dempningskoeffisientene og eksitasjonsintensiteten oppfyller kompatibilitetsbetingelsene, kan den eksakte stasjonære løsningen for p(I1,I2,h3)p(I_1, I_2, h_3) uttrykkes som en kompleks eksponentialfunksjon av de relevante parametrene. Dette gir oss en beskrivelse av hvordan systemet oppfører seg i et stasjonært regime, og det er derfor av avgjørende betydning å kunne identifisere og bruke de nødvendige kompatibilitetsbetingelsene for å finne nøyaktige løsninger.

For systemer som involverer resonnanser, spesielt når ω1=ω2\omega_1 = \omega_2, kan de gjennomsnittlige Itô-ligningene for systemet uttrykkes som et sett av stokastiske differensialligninger. Disse ligningene gir en beskrivelse av hvordan systemet utvikler seg over tid under innvirkning av stokastisk støy og vil være avgjørende for å finne stasjonære løsninger. For resonante systemer, der parametrene ω1\omega_1 og ω2\omega_2 er lik, blir de gjennomsnittlige driften og diffusjonskoeffisientene til systemet relatert til dets grunnleggende dynamikk, og løsningen kan uttrykkes i form av en kompleks eksponensialfunksjon som inkluderer alle de relevante systemparametrene.

Når systemet er underlagt støy eller tilfeldige eksitasjoner, er det ofte nyttig å bruke stokastiske gjennomsnittsmetoder for å forenkle analysen av systemets stasjonære tilstand. I mange tilfeller, spesielt i quasi-Hamiltonske systemer som ikke er fullt integrerbare, vil eksakte stasjonære løsninger være vanskelig å finne uten å gjøre forenklinger. Ved å bruke stokastiske gjennomsnittsmetoder kan vi imidlertid finne omtrentlige løsninger som gir en god tilnærming til systemets reelle oppførsel.

Videre er det viktig å merke seg at i systemer der forskjellige elementer, som for eksempel masser eller fjærer, har forskjellige parametre (som dempningskoeffisienter eller stivhetsverdier), vil oppførselen i det stasjonære tilfellet variere. Når systemet er underlagt stokastisk eksitasjon, vil det også være nødvendig å vurdere hvordan disse parameterne påvirker de stasjonære sannsynlighetsfordelingene. For eksempel vil samsvar mellom visse forhold som α12/K1=α21/K2\alpha_{12}/K_1 = \alpha_{21}/K_2 være avgjørende for at den stasjonære løsningen skal være eksakt.

En annen viktig dimensjon er hvordan resonansforholdene, som ω1=ω2\omega_1 = \omega_2, påvirker systemets dynamikk. Ved resonans kan systemet vise sterkere eller mer vedvarende utsving, som kan ha betydelige implikasjoner for stabiliteten og de stasjonære løsningene. Derfor er det viktig å forstå at resonansfenomener ikke bare påvirker systemets responstid, men også dens langsiktige statistiske egenskaper.

I flere praktiske tilfeller, som ved vibrasjonssystemer med støyinnflytelse, kan man bruke stasjonære sannsynlighetsfordelinger for å estimere hvordan systemet vil oppføre seg over lang tid. Dette er spesielt nyttig når man trenger å forstå hvordan systemet stabiliserer seg etter en periode med tilfeldig eksitasjon.

Når systemet er bygget opp av flere komponenter med samspill, som et 2-DOF (grad av frihet) vibrasjon-impakt system, kan det være behov for ytterligere metoder for å bestemme hvilken type stasjonær løsning som er mest passende. Dette kan innebære enten eksakte løsninger i spesifikke tilfeller eller anvendelsen av stokastiske gjennomsnittsmetoder for å få tilnærmede løsninger når det er vanskelig å finne en eksakt løsning.

En viktig utfordring som ofte oppstår når man jobber med stokastisk dynamikk, er hvordan man kan modellere interaksjoner mellom ulike deler av systemet som kan ha ikke-lineære effekter, for eksempel ved kollisjoner eller kontakt mellom deler av systemet. Her vil bruk av kontaktlover, som for eksempel Hertz kontaktlov i et system med to masser som påvirkes av elastiske vegger, være avgjørende for å forstå systemets respons.

Stasjonære løsninger i quasi-Hamiltonske systemer er derfor ikke bare et spørsmål om matematisk beregning, men også om å forstå hvordan ulike fysikalske interaksjoner, resonansforhold og stokastisk støy kombineres for å forme systemets dynamikk i det stasjonære regimet.

Hvordan stokastisk gjennomsnittsmetode kan anvendes på quasi-Hamiltonianske systemer eksitert av fGn

Quasi-Hamiltonianske systemer som er eksitert av fraksjonell Gaussisk støy (fGn) er matematiske modeller som kan brukes til å beskrive et bredt spekter av tilfeldige fenomener i økonomi, finans, vitenskap og ingeniørfag. Et av de mest karakteristiske trekkene ved fGn er dens langdistanseavhengighet eller "langt minne", noe som gjør at den skiller seg betydelig fra standard hvit støy. Denne typen støy har et korrelasjonsmønster som gjør at tidligere hendelser påvirker fremtidige utfall i systemet over lang tid.

Stokastiske prosesser som involverer fGn er ofte ikke Markov-prosesser, noe som betyr at de ikke kan beskrives ved kun å bruke den nåværende tilstanden i systemet. Imidlertid finnes det metoder som kan forenkle analysen av slike systemer. En slik metode er den stokastiske gjennomsnittsmetoden, som kan brukes til å redusere dimensjonen til systemet, og dermed spare betydelig tid når man simulerer systemets respons via Monte Carlo-simuleringer.

I dette kapittelet introduseres den stokastiske gjennomsnittsmetoden for quasi-Hamiltonianske systemer eksitert av fGn. For slike systemer, der den effektive spektrale tettheten for fGn varierer sakte i mellom- til høyfrekvensområdet, kan fGn i dette frekvensområdet behandles som en bredbåndsprosess. Derfor er metoden som er utviklet for quasi-integrerbare Hamiltonianske systemer eksitert av bredbånds støy også anvendelig i dette tilfellet.

Når vi ser på et n-graders frihet (DOF) quasi-Hamiltoniansk system som er eksitert av fGn, beskrives systemet ved et sett differensialligninger. Disse ligningene involverer generaliserte forskyvninger og momenta som er relatert til Hamilton-funksjonen, som er to ganger deriverbar med hensyn til disse variablene. Den uavhengige enheten av fGn som eksiterer systemet, representeres ved flere fGn-er med forskjellige Hurst-indekser, som kontrollerer lengden på minnet i prosessen.

Ved å bruke stokastisk integrasjonsteknikk for fraksjonell Brownsk bevegelse (fBm) kan bevegelsen i systemet omformes til fraksjonelle stokastiske differensialligninger (SDE-er), som beskriver den dynamiske utviklingen av systemets tilstand. Disse SDE-ene kan videre forenkles ved å bruke gjennomsnittsmetoden, som gjør det mulig å redusere systemets kompleksitet. Det er viktig å merke seg at denne metoden gjelder spesifikke tilfeller, der systemet er eksitert av en kombinasjon av fGn med ulike Hurst-indekser, og hvor de stokastiske integrasjonene er definert på en måte som sikrer at den ikke blir påvirket av Wong-Zakai-korreksjonene.

Videre kan man anvende en spesiell form for den stokastiske gjennomsnittsmetoden for ikke-integrerbare Hamiltonianske systemer, som kan beskrives ved fraksjonelle SDE-er. Her kan prosessen til systemet beskrives som en langsom variabel prosess, som konvergerer til en gjennomsnittlig fraksjonell stokastisk prosess etter at man har utført tidsgjennomsnittlig beregning. Ved å bruke denne tilnærmingen kan man erstatte de kompliserte dynamiske beregningene med en enklere modell som fortsatt fanger opp systemets essens.

For eksempel, i et 2-DOF quasi-non-integrerbart Hamiltoniansk system, kan man bruke den stokastiske gjennomsnittsmetoden til å derivere en en-dimensjonal gjennomsnittlig fraksjonell SDE. Her får man frem de nødvendige koeffisientene som er bestemt av de originale ligningene for systemets Hamiltonian. Ved å bruke Monte Carlo-simuleringer kan man videre hente ut den stasjonære sannsynlighetsfordelingen og andre statistiske egenskaper for systemet.

Når man bruker denne metoden, kan man forvente betydelig tidseffektivisering sammenlignet med mer detaljerte simuleringer som involverer fullstendig dynamisk modellering av systemet. Dette er spesielt nyttig i praktiske anvendelser hvor det er behov for å modellere komplekse systemer med støy, som i økonomiske og ingeniørmessige simuleringer, der de involverte prosessene har lang rekkeviddeavhengighet.

Et viktig aspekt å forstå når man arbeider med slike systemer er at fGn kan føre til dynamikk som er langt mer kompleks enn hva man ville forvente med hvit støy. Langt minne gjør at systemet kan vise seg å være sensitivt overfor tidligere tilstander, noe som skaper en form for avhengighet mellom tilstandene over tid. Dette må tas i betraktning når man analyserer systemets respons eller når man skal bruke modellene i praktiske situasjoner.

Når man gjennomfører slike analyser, bør man også være oppmerksom på at den langsomme variabiliteten som er karakteristisk for fGn, kan føre til utfordringer med konvergens i numeriske metoder. Det er derfor viktig å bruke teknikker som kan håndtere slike utfordringer, som for eksempel Monte Carlo-simuleringer, som gir statistisk robuste resultater til tross for de komplekse dynamiske prosessene som er involvert.