[elektro] Goertzel algoritmus

Info info at kiralyelektronika.hu
Sat Feb 14 23:21:08 CET 2009


Sziasztok !

Végre a wikipédiás példából sikerült beizzítanom első FFT-s
algoritmusom :) ennek okán kérdeznék:

Init:
              coef = 2 * omega / samplerate
Mintánként:
              q[0] = coefs * q[-1] - q[-2] + sample
X periódusonként:
              amplitudo = q[0]^2  + q[-1]^2 - (coefs * q[0] * q[-1])

Nos, ez müxik, eljátszottam vele. Az említett forrásban 8k mintavétel
van és kb. 6 periódusonként értékeli ki az eredményt.
Felezve a mintavételt is müxik, és grafikonon kirajzoltam az
amplitudókat, láthatóan csökkent a jel-zaj viszony.
Duplázva a mintavételt kicsit növelni kellett a kiértékeléshez a mért
periódusok mennyiségét, a grafikonon érdekesen javult a jel-zaj
viszony: olyan szinuszosan nézett ki.

Nemtom, hogy jó-e a fenti számítási módszer, azonban kíváncsi lennék,
hogy hogyan javítható eme digitális szűrő fokszáma, és hogyan kapok
vissza skálázott értéket. Tehát mondjuk bemenő +-1024 ami +-10V lenne
AD előtt az a kimeneten milyen mértékegységben lesz ?

Thx !
Béla




More information about the Elektro mailing list