ik zoek een autocorrelatie rekentruuk

Integraalrekening, afgeleiden, rijen, convergentie & divergentie van reeksen, meervoudige integratie.
wimapon
Nieuw lid
Nieuw lid
Berichten: 5
Lid geworden op: 08 jan 2013, 11:08

ik zoek een autocorrelatie rekentruuk

Bericht door wimapon » 08 jan 2013, 16:23

LS,
Ik ben op zoek naar een rekentruuk om een autocorrelatie op een rij van 300000 getallen los te laten.

Ik ontvang gedurende 15 seconden ruis afkomstig van tussen de sterren. ( SETI - ontvangstation)
Ik laat daar een FFT op los en krijg een prachtig frequentie spectrum.
Als een signaal inzit komt er een "paal" in de grafiek te staan.

Ik heb begrepen dat men tegenwoordig een truuk gebruikt om zich herhalende signalen die in die
ruis zitten eruit te vissen.
Men doet dit door een autocorrelatie op die ruis los te laten.

Ik kan geen wiskundige formules lezen... dus mijn vraag is nu: hoe bereken ik numeriek een
autocorrelatie.

UIt Wikipedia komt de volgende ( voor mijn onbegrijpelijke info)

FR(f) = FFT (X(t)) ;; dus een gewone FFT op de rij getallen die der ruis voorstellen... ( lukt me prima)
S(f) = FR(f) * FR*(f) ;; FR*(f) noemt men de complex geconjugeerde (voor mij abacadabra)
R(tau) = IFFT(S(f) ;; de inverse FFT.. dat lukt me ook

Bovenstaande info komt uit: http://en.wikipedia.org/wiki/Autocorrelation ( efficient computation)


Ik doe de eerste retel:
ik doe een FFT de rij getallen en krijg 2 vectoren: REX(i) en IMX(i)

Nu moet ik de tweede regel uitvoeren.......= abacadabra

De derde regel zal wel lukken als ik weer gevulde REX(i) en IMX(i) heb.


Ik programmeer in Pure Basic .. en kan geen beschikbare routine vinden. dus ik moet het
zelf programmeren.. maar dan moet ik eerst snappen wat er gebeurd...

Kan iemand me helpen

Wim Apon

Plaats reactie