Flytte Gjennomsnittet Spektral Estimering

12.1: Estimering av spektral tetthet Vi har tidligere diskutert periodogrammet, en funksjonsgraf som viser informasjon om periodiske komponenter i en tidsserie. Eventuelle tidsserier kan uttrykkes som summen av cosinus og sinusbølger som svinger ved de grunnleggende (harmoniske) frekvensene jn. med j 1, 2,, n 2. Periodogrammet gir informasjon om de relative frekvensene for de ulike frekvensene for å forklare variasjonen i tidsseriene. Periodogrammet er et utvalgsestimat av en populasjonsfunksjon kalt spektral tetthet, som er en frekvensdomæne karakterisering av en befolkningsstasjonær tidsserie. Spektral tettheten er en frekvens domener representasjon av en tidsserie som er direkte relatert til autokovarians tid domene representasjon. I hovedsak inneholder spektral tetthet og autokovariansfunksjonen den samme informasjonen, men uttrykker den på forskjellige måter. Gjennomgå Merknad. Autokovariansen er telleren av autokorrelasjonen. Autokorrelasjonen er autokovariansen delt av variansen. Anta at (h) er autokovariansfunksjonen til en stasjonær prosess, og at f () er spektral tettheten for samme prosess. I notasjonen til foregående setning, h tidsforsinkelse og frekvens. Autokovariansen og spektral tettheten har følgende forhold: På språket av avansert beregning er autokovarians og spektral tetthet Fourier transformasjonspar. Vi vil ikke bekymre oss om beregningen av situasjonen. Godt fokus på estimeringen av spektral tettheten frekvensdomene karakterisering av en serie. Fourier-transformasjonsligningene er bare gitt her for å fastslå at det er en direkte kobling mellom tidsdomenerepresentasjonen og frekvensdomenerrepresentasjonen av en serie. Matematisk defineres spektral tettheten for både negative og positive frekvenser. På grunn av funksjonens symmetri og dens repeterende mønster for frekvenser utenfor området -12 til 12, trenger vi imidlertid bare å være opptatt av frekvenser mellom 0 og 12. Den totale integrerte spektraldensiteten er lik variansen i serien. Dermed kan spektral tettheten innenfor et bestemt frekvensintervall betraktes som mengden av variansen forklart av disse frekvensene. Metoder for estimering av spektral tetthet Råperiodogrammet er et grovt utvalgsestimat av populasjonsspektraldensiteten. Estimatet er grovt, delvis fordi vi bare bruker de diskrete grunnleggende harmoniske frekvensene for periodogrammet, mens spektral tettheten er definert over et kontinuum av frekvenser. En mulig forbedring av periodogrammet estimat av spektral tetthet er å glatte det ved å bruke sentrert glidende gjennomsnitt. En ytterligere utjevning kan opprettes ved hjelp av tapering metoder som vekter endene (i tid) av serien mindre enn sentrum av dataene. Ikke dekk tapende i denne leksjonen. Interesserte kan se Seksjon 4.5 i boken og ulike Internett-kilder. En alternativ tilnærming til å utjevne periodogrammet er en parametrisk estimeringsmetode basert på det faktum at en hvilken som helst stasjonær tidsserie kan tilnærmet ved en AR-modell av noen rekkefølge (selv om det kan være en høy ordre). I denne tilnærmingen er det funnet en egnet AR-modell, og deretter beregnes spektral tettheten som spektral tettheten for den estimerte AR-modellen. Utjevningsmetode (Nonparametric Estimation of Spectral Density) Den vanlige metoden for å utjevne et periodogram har et så fint navn som det høres vanskelig ut. Faktisk er det bare en sentrert glidende gjennomsnittlig prosedyre med noen få mulige modifikasjoner. For en tidsserie er Daniell-kjernen med parameter m et sentrert glidende gjennomsnitt som skaper en glatt verdi ved tid t ved å beregne alle verdier mellom tidene tm og tm (inkluderende). For eksempel er utjevningsformelen for en Daniell-kjerne med m 2 In R, kan vektningskoeffisientene for en Daniell-kjerne med m 2 genereres med kommandokjernen (daniell, 2). Resultatet er coef-2 0,2 ​​coef-1 0,2 coef 0 0,2 coef 1 0,2 coef 2 0,2 ​​Abonnementene for coef refererer til tidsforskjellen fra midten av gjennomsnittet ved tid t. Dermed er utjevningsformelen i dette tilfellet som er den samme som formelen gitt ovenfor. Den modifiserte Daniell-kjernen er slik at de to endepunktene i gjennomsnittet mottar halv vekten som interiørpunktene gjør. For en modifisert Daniell-kjerne med m 2, er utjevningen i R, vil kommandokjernen (modifisert. daniell, 2) angi vektningskoeffisientene som nettopp er brukt. Enten kan Daniell-kjernen eller den modifiserte Daniell-kjernen bli viklet (gjentatt) slik at utjevningen påføres igjen til de jevne verdiene. Dette gir en mer omfattende utjevning ved gjennomsnittsnivå over et bredere tidsintervall. For eksempel å gjenta en Daniell-kjerne med m 2 på de glattede verdiene som resulterte fra en Daniell-kjerne med m 2, ville formelen være Dette er gjennomsnittet av glattede verdier innen to tidsperioder t. i begge retninger. I R vil kommandokjernen (daniell, c (2,2)) levere koeffisientene som vil bli brukt som vekt i gjennomsnitt av de opprinnelige dataværdiene for en innviklet Daniell-kjerne med m 2 i begge smoothings. Resultatet er gt kernel (daniell, c (2,2)) coef-4 0,04 coef-3 0,08 coef-2 0,12 coef-1 0,16 coef 0 0,20 coef 1 0,16 coef 2 0,12 coef 3 0,08 coef 4 0,04 Dette gir utjevning formel En konvolusjon av den modifiserte metoden hvor endepunktene har mindre vekt er også mulig. Kommandokjernen (modifisert. daniell, c (2,2)) gir disse koeffisientene: coef-4 0,01563 coef-3 0,06250 coef-2 0.12500 coef-1 0.18750 coef 0 0.21875 coef 1 0.18750 coef 2 0.12500 coef 3 0.06250 coef 4 0.01563 Dermed er senterverdiene vektet litt tyngre enn i den umodifiserte Daniell-kjernen. Når vi slipper et periodogram, er vi utjevning over et frekvensintervall i stedet for et tidsintervall. Husk at periodogrammet bestemmes ved de grunnleggende frekvensene j jn for j 1, 2, n 2. La jeg (j) angi periodogramverdien ved frekvens j jn. Når vi bruker en Daniell-kjerne med parameter m for å jevne et periodogram, er den glatte verdien (hat (omegaj)) et vektet gjennomsnitt av periodogramverdier for frekvenser i området (j-m) n til (jm) n. Det er L 2 m 1 grunnfrekvensverdier i området (j-m) n til (jm) n. Utvalget av verdier som brukes til utjevning. Båndbredden for det glatte periodogrammet er definert som Båndbredden er et mål på bredden på frekvensintervallet (ene) som brukes til å jevne ut periodogrammet. Når ujevne vekter brukes i utjevning, defineres båndbreddedefinisjonen. Angi den glatte periodogramverdien ved j jn som hatt (omegaj) sum hk Jeg forlot (omegaj frac høyre). Hk er de muligens ulikvektene som brukes i utjevningen. Båndbreddeformelen blir da endret til Egentlig fungerer denne formelen også for likevekter. Båndbredden skal være tilstrekkelig til å jevne ut vårt estimat, men hvis vi bruker en båndbredde som er for stor, glatt ut periodogrammet for mye og savner å se viktige topper. I praksis tar det vanligvis noen eksperimenter for å finne båndbredden som gir en passende utjevning. Båndbredden styres hovedsakelig av antall verdier som er i gjennomsnitt i utjevningen. Med andre ord, m-parameteren for Daniell-kjerne og om kjerne er innviklet (gjentatt), påvirker båndbredden. Merk: Båndbredder R rapporterer med sine tegninger samsvarer ikke med verdiene som skal beregnes ved hjelp av formlene ovenfor. Vennligst se fotnoten på s. 197 av teksten din for en forklaring. Averagingsmoothing periodogrammet med en Daniell-kjernen kan oppnås i R ved å bruke en sekvens av to kommandoer. Den første definerer en Daniell-kjerne, og den andre skaper det glatte periodogrammet. For eksempel, anta at den observerte serien heter x og vi ønsker å jevne periodogrammet ved hjelp av en Daniell-kjerne med m 4. Kommandoene er k-kjerne (daniell, 4) spec. pgram (x, k, taper0, log no) Den første kommandoen skaper vektningskoeffisientene som trengs for utjevning og lagrer dem i en vektor kalt k. (Det er tilfeldig å kalle det k. Det kan kalles noe.) Den andre kommandoen ber om et spektral tetthets estimat basert på periodogrammet for serien x. ved hjelp av vektningskoeffisientene lagret i k, uten tap, og plottet vil være på vanlig skala, ikke en loggskala. Hvis en konvolusjon er ønsket, kan kjernekommandoen bli endret til noe som k kjerne (daniell, c (4,4)). Det er to mulige måter å oppnå en modifisert Daniell-kjernen på. Du kan enten endre kjernekommandoen for å referere til modifisert. daniell i stedet for daniell, eller du kan hoppe over kjernekommandoen og bruke en spans parameter i kommandoen spec. pgram. Spansparameteren gir lengden (2 m 1) av den ønskede modifiserte Daniell-kjernen. For eksempel har en modifisert Daniell-kjerne med m 4 lengde L 2 m 1 9 slik at vi kunne bruke kommandoen spec. pgram (x, spans9, taper 0, logno) To passerer av en modifisert Daniell-kjerne med m 4 på hvert pass kan gjøres ved hjelp av spec. pgram (x, spansc (9,9), konisk 0, logno) Eksempel. Dette eksemplet vil bruke fiskerekruteringsserien som brukes flere steder i teksten, inkludert flere steder i kapittel 4. Serien består av n 453 månedlige verdier av et mål på en fiskpopulasjon på en sørlig halvkuleplassering. Dataene er i filen rekruttering. dat. Råperiodogrammet kan opprettes ved hjelp av kommandoen (eller det kan opprettes ved hjelp av metoden gitt i leksjon 6). spec. pgram (x, taper0, logno) Merk at i kommandoen bare gitt har vi utelatt parameteren som gir vekter for utjevning. Det røde periodogrammet følger: Det neste diagrammet er et jevnt periodogram med en Daniell-kjerne med m 4. Merk at en utjevning av utjevning er at den dominerende toppen i den ujevne versjonen nå er den nest høyeste toppen. Dette skjedde fordi toppen er så skarpt definert i den ujevne versjonen at når vi gjennomsnittlig det med noen omgivende verdier, blir høyden redusert. Det neste plottet er et jevnt periodogram ved bruk av to passerer av en Daniell-kjerne med m 4 på hvert pass. Legg merke til hvordan det er jevnere enn tidligere. For å lære hvor de to dominerende toppene er plassert, tilordne et navn til spec. pgram-utgangen, og så kan du liste den. For eksempel, spesifikasjoner spec. pgram (x, k, taper0, logno) specvalues ​​Du kan sile gjennom utgangen for å finne frekvensene hvor toppene oppstår. Frekvensene og spektral tetthetsestimatene er oppført separat, men i samme rekkefølge. Identifiser maksimale spektraldensiteter og finn deretter de tilsvarende frekvensene. Her er den første toppen med en frekvens .0229. Perioden (antall måneder) assosiert med denne syklusen 1.0229 43.7 måneder, eller ca 44 måneder. Den andre topp forekommer med en frekvens på 0,083333. Den tilknyttede perioden 1.08333 12 måneder. Den første toppen er forbundet med en El Nino vær-effekt. Den andre er den vanlige 12 måneders sesongvirkningen. Disse to kommandoene vil sette vertikale prikkede linjer på den (estimerte) spektraldensitetsplottet ved omtrentlige plasseringer av toppdensiteter. abline (v144, ltydotted) abline (v112, lty prikket) Heres den resulterende plottet: Weve glatt nok, men for demonstrasjonsformål er neste plot resultatet av spec. pgram (x, spansc (13,13), taper0, logno ) Dette bruker to passeringer av en modifisert Daniell-kjerne med lengde L 13 (så m 6) hver gang. Tomten er litt jevnere, men ikke så mye. Toppene, forresten, er på nøyaktig samme steder som i tomten umiddelbart over. Det er definitivt mulig å glatte for mye. Anta at vi skulle bruke en modifisert Daniell-kjerne med total lengde 73 (m 36). Kommandoen er spec. pgram (x, spans73, taper0, logno) Resultatet følger. Toppene er borte Parametrisk estimering av spektral tetthet Utjevningsmetoden for spektral tetthets estimering kalles en nonparametrisk metode fordi den ikke bruker noen parametrisk modell for den underliggende tidsseriprosessen. En alternativ metode er en parametrisk metode som innebærer å finne den beste passende AR-modellen for serien og deretter plotte spektraldensiteten til modellen. Denne metoden støttes av en teorem som sier at spektral tettheten i en hvilken som helst tidsserieprosess kan tilnærmet av spektraldensiteten til en AR-modell (av en eller annen rekkefølge, muligens en høy). I R, er parametrisk estimering av spektral tetthet enkelt gjort med kommandofunksjonen spec. ar. En kommando som spec. ar (x, logno) vil føre til at R gjør alt arbeidet. Igjen, for å identifisere topper kan vi tildele spesifikke resultater ved å gjøre noe som specvaluesspec. ar (x, log no). For eksempel på rekruttering av fisk, er følgende tomt resultatet. Legg merke til at tettheten plottet er den av en AR (13) modell. Vi kan sikkert finne flere parsimoniske ARIMA-modeller for disse dataene. Brukte bare spektraldensiteten til modellen til å omtrentliggjøre spektral tettheten til den observerte serien. Utseendet til estimert spektral tetthet er omtrent det samme som før. Den estimerte El Nino-toppen ligger på et litt annerledes sted, frekvensen er ca. 0,024 for en syklus på ca. 1.024 omtrent 42 måneder. En serie skal de-trender før en spektralanalyse. En trend vil føre til en slik dominerende spektral tetthet i lav frekvens at andre topper ikke vil bli sett. Som standard utfører R-kommandoen spec. pgram en trend med en lineær trendmodell. Det vil si at spektral tettheten estimeres ved bruk av residualene fra en regresjon gjort hvor de y-variable observerte data og x-variabelen t. Hvis en annen type trend er tilstede, kan en kvadratisk for eksempel bruke en polynomisk regresjon til å de-trene dataene før den estimerte spektraldensiteten utforskes. Vær imidlertid oppmerksom på at R-kommandoen spec. ar. men utfører ikke en de-trending som standard. Anvendelse av Smoothers til Raw Data Merk at smoothers beskrevet her også kunne brukes på rå data. Daniell-kjernen og dens modifikasjoner beveger seg ganske enkelt gjennomsnittlig (eller vektet glidende gjennomsnittlig) glatteevne. Navigation16. Spektral estimering Spektral estimeringsproblemet for en diskret tidsserie generert av en lineær, tidsinvariant prosess kan formuleres i form av tre modeller: autoregressive (AR), glidende gjennomsnitt (MA) og autoregressivt bevegelige gjennomsnitt (ARMA). Analyseprosedyrene varierer hver enkelt, og spesifikasjonsfeil oppstår på grunn av bruk av den uhensiktsmessige algoritmen. AR - og MA-modellene fører henholdsvis til maksimal entropi (MEM) og klassisk lag-vindu-tilnærminger. ARMA-modellen har mye seismisk interesse bemuse enhetens impulsrespons av et horisontalt stratifisert medium er uttrykkelig på denne måten. Siden tilbakemeldingskomponenten har den laveste forsinkelsesegenskapen, har en ARMA spektral estimeringsteknikk som tilfredsstiller dette kravet spesiell seismisk relevans. Et slikt spektralestimat resulterer fra anvendelsen av en iterativ minst-kvadrat-algoritme til utvalgte portene til den observerte tidsserien. Et utvalg av syntetiske tidsserier tjener til å illustrere nedbrytningen i spektralestimatet som følge av en feilaktig spesifikasjon av modellen. Mye har blitt skrevet i de senere år om spektralanalyse av diskrete tidsserier. Det eksisterer ingen enkelt korrekt teknikk for å beregne spektret i fravær av kunnskap om hvilken type prosess som har generert dataene. Som vi har sett i kapittel 9, skiller vi mellom tre mulige prosesser: autoregressiv (AR), glidende gjennomsnitt (MA) og autoregressivt bevegelige gjennomsnitt (ARMA). I tekniske termer beskriver disse prosessene allpolig (eller tilbakemelding), all-zero (eller feedforward) og pole-zero (eller feedbaek-feedforward) - systemene. Generelt sett har vi ikke forutgående kunnskap om tidsserienes generasjonsmekanisme, og vi er tvunget til å anta at våre registrerte data faktisk tilfredsstiller en av disse tre representasjonene. Når denne avgjørelsen er gjort, må vi velge en passende algoritme for beregning av det faktiske spektralestimatet. I tilfelle av AR - eller allpoligmodellen er den maksimale entropi-metoden (MEM) som implementert med en teknikk grunnet Burg (1967, 1975) hensiktsmessig. For MA, eller all-zero modellen, benytter vi klassisk lag-vindu-tilnærming (Blackman og Tukey, 1959). I tillegg 16-1 gir vi matematikken til den klassiske lag-vindu-metoden, og i tillegg 16-2, matematikken til maksimal entropi-metoden. ARMA - eller pole-zero-modellen har også fått oppmerksomhet i nyere litteratur: Relevante spektrale estimateknikker er beskrevet av Anderson (1971, kap. 5), av Box og Jenkins (1970, kap. 6 og 7), og ved Alam (1978). Den rasjonelle representasjonen av impulsresponsen til en ARMA-prosess er gitt ved forholdet mellom to polynomene i den komplekse variabelen z. I dette kapitlet skal vi være spesielt interessert i spektralanalyse av seismogrammer. Som vi har sett i kapittel 13, kan enhetens impulsrespons av et perfekt elastisk, horisontalt stratifisert medium uttrykkes som forholdet mellom to slike polynomier i kraften av z. men med den ekstra begrensningen at nevneren polynom har den laveste forsinkelsesegenskapen. Med andre ord tvinger denne tilstanden polene til systemet til å ligge utenfor periferien av enhetens sirkel z 1 i det komplekse planet, og tillater oss å utvide ARMA-polynomforholdet i form av en konvergerende kraftserie i z. Det vil derfor være ønskelig å søke en ARMA spektral estimeringsalgoritme som garanterer en minimumsforsinkelsesnevner. Selv om det ikke er noen grunnleggende matematisk behov for en ARMA spektral estimeringsmetode for å produsere en minimumsforsinkelsesnevner, har vi nettopp uttalt at en slik oppgave har sterk fysisk motivasjon. Følgelig er minimumsforsinkelsesegenskapen til nevnen et sterkt punkt, og en ikke nødvendigvis. delt av andre ARMA spektral estimatorer. Innholdsfortegnelse Er du medlem av SEG eller EEGS Hvis du er SEG-medlem (med tilgang til SEG - og EEGS-tidsskrifter, utvidede abstrakter og prosesser og nedsatte medlemspriser for individuelle kjøp av SEG-bøker) eller du har allerede kjøpt tilgang til dette innhold separat, klikk her for å logge inn og få tilgang til ønsket innhold. Hvis du er et EEGS-medlem (med tilgang til EEGS-publikasjoner og SEG Technical Program Expanded Abstracts), klikk her for å logge inn og få tilgang til ønsket innhold. Alt eBok innhold kan kjøpes separat av enkeltpersoner og institusjoner. Kjøp dette innholdet Velg blant følgende alternativer: SEG Digital Library Institutional Options SEG har utvidet sin samling av online bøker til mer enn 100 titler med en kombinasjon av nye og eldre verker, og har lagt til et evigvarende kjøpskjøp for hele samlingen. Institusjonelle abonnementer til SEG Digital Library (inkludert SEG og EEGS tidsskrifter, utvidede abstrakter og prosesser) er også tilgjengelige for kjøp. Mens SEG-tidsskrift og magasininnhold er tilgjengelig for SEG og EEGS-medlemmer gratis som en del av medlemskapet, må ebøker kjøpes individuelt (med rabatt til medlemmer), eller via abonnement. Utviklingsgruende gjennomsnittlig modellering for spektralparameterestimering fra en multigradient ekko kjemisk skift oppkjøp T1 - Autoregressiv glidende gjennomsnittlig modellering for spektral parameter estimering fra et multigradient ekko kjemisk skift oppkjøp AU - Taylor, Brian A. AU - Hwang, Ken Pin AU - Hazle, John D. AU - Stafford, R. Jason N2 - Forfatterne undersøkte ytelsen til den iterative Steiglitz-McBride (SM) algoritmen på en autoregressiv glidende gjennomsnittlig (ARMA) modell av signaler fra et raskt, tynt samplet, multiecho, kjemisk skiftbilleddannelse (CSI) oppkjøp ved hjelp av simulering, fantom, ex vivo og in vivo eksperimenter med fokus på potensiell bruk i magnetisk resonans (MR) - styrte inngrep. ARMA-signalmodellen muliggjorde en rask beregning av kjemisk skift, tilsynelatende spin-spin-avslapningstid (T2) og komplekse amplituder av et multipeaksystem fra et begrenset antall ekkoer (16). Numeriske simuleringer av en - og to-toppsystemer ble brukt til å vurdere nøyaktigheten og usikkerheten i de beregnede spektrale parametrene som en funksjon av oppkjøps - og vevsparametere. De målte usikkerhetene fra simuleringen ble sammenlignet med den teoretiske Cramer-Rao nedre grensen (CRLB) for oppkjøpet. Målinger laget i fantomer ble brukt til å validere T2 estimatene og å validere usikkerhetsestimater fra CRLB. Vi viste søknad om realtids MR-styrte inngrep ex vivo ved bruk av teknikken til å overvåke en perkutan etanolinjeksjon i en bovint lever og in vivo for å overvåke en laserinducert termisk terapibehandling i en hjørnetannshjerne. Simuleringsresultater viste at kjemisk skift og amplitude usikkerhet nådde deres respektive CRLB ved et signal-støyforhold (SNR) 5 for ekkologlengder (ETLs) 4 ved hjelp av en fast ekkoavstand på 3,3 ms. T2 estimater fra signalmodellen hadde større usikkerheter, men nådde CRLB ved større SNRs og eller ETLs. Meget nøyaktige estimater for kjemisk skift (lt0.01 ppm) og amplitude (lt1.0) ble oppnådd med 4 ekkoer og for T2 (lt1.0) med 7 ekkoer. Vi konkluderer med at SM-algoritmen over et rimelig SNR-område er en robust estimator for spektrale parametere fra raske CSI-oppkjøp som kjøper 16 ekkoer for en - og to-toppsystemer. Preliminære ex vivo og in vivo eksperimenter bekreftet resultatene fra simuleringseksperimenter og angir videre potensialet av denne teknikken for MR-styrte intervensjonsprosedyrer med høy spatiotemporal oppløsning 1,61,64 mm3 på 5 s. 2009 American Association of Physics in Medicine. AB - Forfatterne undersøkte ytelsen til den iterative Steiglitz-McBride (SM) algoritmen på en autoregressiv glidende gjennomsnittlig (ARMA) modell av signaler fra et raskt, sparsomt samplet, multiecho, kjemisk skiftbilleddannelse (CSI) oppkjøp ved hjelp av simulering, fantom, ex vivo og in vivo eksperimenter med fokus på potensiell bruk i magnetisk resonans (MR) - styrte inngrep. ARMA-signalmodellen muliggjorde en rask beregning av kjemisk skift, tilsynelatende spin-spin-avslapningstid (T2) og komplekse amplituder av et multipeaksystem fra et begrenset antall ekkoer (16). Numeriske simuleringer av en - og to-toppsystemer ble brukt til å vurdere nøyaktigheten og usikkerheten i de beregnede spektrale parametrene som en funksjon av oppkjøps - og vevsparametere. De målte usikkerhetene fra simuleringen ble sammenlignet med den teoretiske Cramer-Rao nedre grensen (CRLB) for oppkjøpet. Målinger laget i fantomer ble brukt til å validere T2 estimatene og å validere usikkerhetsestimater fra CRLB. Vi viste søknad om realtids MR-styrte inngrep ex vivo ved bruk av teknikken til å overvåke en perkutan etanolinjeksjon i en bovint lever og in vivo for å overvåke en laserinducert termisk terapibehandling i en hjørnetannshjerne. Simuleringsresultater viste at kjemisk skift og amplitude usikkerhet nådde deres respektive CRLB ved et signal-støyforhold (SNR) 5 for ekkologlengder (ETLs) 4 ved hjelp av en fast ekkoavstand på 3,3 ms. T2 estimater fra signalmodellen hadde større usikkerheter, men nådde CRLB ved større SNRs og eller ETLs. Meget nøyaktige estimater for kjemisk skift (lt0.01 ppm) og amplitude (lt1.0) ble oppnådd med 4 ekkoer og for T2 (lt1.0) med 7 ekkoer. Vi konkluderer med at SM-algoritmen over et rimelig SNR-område er en robust estimator for spektrale parametere fra raske CSI-oppkjøp som kjøper 16 ekkoer for en - og to-toppsystemer. Preliminære ex vivo og in vivo eksperimenter bekreftet resultatene fra simuleringseksperimenter og angir videre potensialet av denne teknikken for MR-styrte intervensjonsprosedyrer med høy spatiotemporal oppløsning 1,61,64 mm3 på 5 s. 2009 American Association of Physics in Medicine. KW - Autoregressive glidende gjennomsnitt (ARMA) KW - Kemisk skift avbildning (CSI) KW - MR-styrte inngrep KW - Multigradient ekkooppkjøpAutoregressiv Moving Average Spectral Estimation: A Model Equation Error ProsedyrePeriodicity og trend for multisite daglige nedbørstidsserier og utviklede simuleringskriterier ble modellert for å oppnå en lengre tidsserie 22. Det er flere prosedyrer for periodisk estimering som spektralanalyse 23 24. Auto-regressiv glidende gjennomsnittlig metode 25 28 og enkel harmonisk analyse 2932. Selv om flere studier på temperaturtendens har blitt gjennomført i Canada , svært få av dem demonstrerer periodicity komponenter med hensyn til Ontario forhold. sitat Fulltekst Artikkel Jan 2014 Syed Imran Ahmed Ramesh Rudra Trevor Dickinson Motahir Ahmed Vis abstrakt Skjul abstrakt ABSTRAKT: En ny adaptiv metode for effektivt å skaffe et ARMA-modellspektralestimat av en bredformet stasjonær tidsserie presenteres. Det er adaptivt i den forstand at som et nytt element i tidsseriene blir observert, kan koeffisientene til en ARMA-modell (p, p), bli algoritmisk oppdatert. Denne algoritmx27s beregningskompleksitet (dvs. antall multiplikasjoner og tillegg som kreves) er av rekkefølgen p log (p) for en bestemt versjon av fremgangsmåten. Videre er den spektrale estimeringsytelsen av denne nye metoden funnet å være langt bedre enn slike moderne tilnærminger som Box-Jenkins, maksimal entropi og Widrowx27s LMS-metoder. Denne ytelsen i forbindelse med dens beregningseffektivitet markerer denne algoritmen som et primært spektral estimeringsverktøy. Konferansepapir Mai 1981 IEEE-transaksjoner på akustikk Tale - og signalbehandling JA Cadzow K. Ogino Vis abstrakt Skjul abstrakt ABSTRAKT: I dette papiret er effektive metoder for å generere todimensjonale kvartplanet kausal autoregressive (AR) og autoregressive glidende gjennomsnittlig (ARMA) spektral estimeringsmodeller er utviklet. Disse prosedyrene er funnet å gi superoppløsningskapasiteter sammenlignet med andre mer klassiske metoder som Fourier-transformasjonen. ARMA-metoden innebærer manipulering av modellekvationen summin max summin maks ax (n - k, n - m) summin maks summin maks b epsilon (n - k, n - m) og benytter det angitte endelige settet av observasjoner x n) for 1 leq n leq N, 1 leq n leq N. I det ovennevnte forhold blir den tilfeldige eksitasjon, n) ansett for å være hvit. Denne ARMA modelx27s autoregressive km-koeffisientene er valgt for å minimere et vektet minste kvadratkriterium som består av feilelementer mens de bevegelige gjennomsnittlige b km-koeffisientene oppnås ved hjelp av en alternativ tilnærming. Den spektrale estimeringsprestasjonen av AR - og ARMA-metodene vil bli demonstrert empirisk ved å vurdere problemet med å løse to sinusoider innebygd i støy. Artikkel Jul 1981 James A. Cadzow Koji Ogino