|
Anzeige
|
Hanning Window
Gast
#1130068
Hallo @ all. Ich habe eine FFT programmiert. Jedoch noch ohne Fensterung. Hierzu habe ich einige Fragen: 1. Muss ich die Fensterung z.B. Hanning Window mit dem abgetasteten diskreten Signal multiplizieren? 2. Wo bekomme ich die Koeffizienten für das Hamming Window her? 3. Wie rechne ich aus dem Ergebnis der FFT (Betrag aus Real und Imaginärteil )die Fensterung wieder raus.
Gast
#1130082
Hier schon mal ein guter Link wie man die Fensterfunktion berechnet. Jedoch wie man wieder wieder die Fensterfunktion im Frequenzbereich rausrechnet steht dort nicht. http://www.statistics4u.com/fundstat_germ/ee_fft_windowing.html Johann wrote: > 1. Muss ich die Fensterung z.B. Hanning Window mit dem abgetasteten > diskreten Signal multiplizieren? Du meinst das Hann-Fenster. Es gibt "Hamming", "(von) Hann" - aber nicht "Hanning". Ja, du multiplizierst die Abgetasteten Werte (den Signalausschnitt der Länge N) mit einer Gewichtungsfunktion der Länge N. weighted(n) = original(n) * window(n) (für n = 0..N-1) > 2. Wo bekomme ich die Koeffizienten für das Hamming Window her? Aus der Bildungsvorschrift: http://de.wikipedia.org/wiki/Fensterfunktion#von-Hann-Fenster_.28auch_bekannt_als_.E2.80.9ERaised-Cosinus-Fenster.E2.80.9C.29 > 3. Wie rechne ich aus dem Ergebnis der FFT (Betrag aus Real und > Imaginärteil )die Fensterung wieder raus. Wozu? Im Spektrum ist die Transformierte der Gewichtungsfunktion in das "ursprüngliche" Spektrum gefaltet (vgl. "Leck-Effekt"). Das ist genauso gewollt und je nachdem welche Gewichtungsfunktion man wählt hat man auf der einen Seite gewisse Vorteile die man sich durch Abstriche an anderen Stellen erkauft (z.B. amplitudengenaues vs. frequenzgenaues Spektrum).
Gast
#2354500
T. H. schrieb: > z.B. amplitudengenaues vs. frequenzgenaues Spektrum Welche Möglichkeiten gäbe es, diese Aspekte zu kombinieren? Also wie kann ich "fenstern" um z.B. sehr genaue Amplituden und andererseits sehr genau Frequenzen zu bekommen? Und: Was bedeutet "frequenzgenau"?
Gast
#2359975
Johann schrieb: > Wo bekomme ich die Koeffizienten für das Hamming Window her? Hann ist eine einfache Sinus-Funktion, passend verschoben und als Glocke genutzt. Hamming ist ähnlich, leicht gestaucht und etwas angehoben- ein Kompromiss zwischen geringem Leckeffekt und Repräsentanz der Spektrallinien. Schau Dir auch mal die Gauss-Funktion an, sowie das blackman-Fenster.
Gast
#2366988
Wo könnte man mal nachlesen, welches Fesnter welche Vor und Nachteile hat?
Gast
#2371517
http://old.iem.at/~majdak/alg/projectsSS2002/7/ http://www.comsol.ch/content.php?si=317&id=196 http://www.medialab.ch/archiv/pdf_studien_diplomarbeiten/1sa00/1sa2000-109_SpektrumAnalyzerMitDSP%20(Burri&Schoeb).pdf http://www.imw.tu-clausthal.de/fileadmin/Bilder/Forschung/Publikationen/Mitt_1997/97_14.pdf
Gast
#2512867
VHDL-mässig ist Bartlett am einfachsten.
Gast
#2512915
...ich habe das mal eben aus einem uralten Programm meiner selbst
rausge-
copy-pasted, ohne dass ich das Wochenende darauf verwenden will.
Du brauchst irgendwo
#define SQUARE 1
#define Parzen 2 usw...
Die Funktion erzeugt einen Buffer, dessen Inhalt mit dem Eingabevektor
der FFT multipliziert werden muss, Position für Position. Nach der
Vektor-Mul muss der Eingabepuffer für die FFT in der Mitte große
Zahlen haben & an Anfang & Ende kleine.
Der Buffer time_window muss natürlich genauso groß sein wie Deine
FFT-Daten.
Gruß, Gerhard
-----------------------------
/* Berechnen einer Fensterfunktion. Die Werte werden im Gleitkomma-Array
* window[0..npoints-1] abgelegt.
*/
void mk_window(int wtype, int win_size) {
double f;
int j;
switch(wtype){
case SQUARE:
/* statt 0.999 waere eigentlich 1.0 korrekt, so werden aber
* eventuelle Rundungsfehler beim Uebergang zu Integers harmlos
gemacht.
*/
for(j=0; j<win_size; j++) time_window[j] = 0.999;
break;
case PARZEN:
for(j=0; j<win_size; j++){
time_window[j] =
1.0 - ( fabs( (j - 0.5 * (win_size-1)) / ( 0.5 *
(win_size+1))));
}
break;
case WELCH:
for(j=0; j<win_size; j++){
f= j - 0.5 * (win_size-1);
f = f / ( 0.5 * (win_size + 1));
time_window[j] = 1 - f * f;
}
break;
case HANN:
for(j=0; j<win_size; j++){
time_window[j] = 0.5 * ( 1 - cos( 2 pi j / (win_size-1)));
}
break;
case HAMMING:
for(j=0; j<win_size; j++){
time_window[j] = 0.54 -0.46 * cos( 2 pi j / (win_size-1));
}
break;
case BLACKMAN:
for(j=0; j<win_size; j++){
time_window[j] = 0.42
- 0.5 * cos( 2 pi j / (win_size-1))
+ 0.08 * cos( 4 pi j / (win_size-1));
}
break;
case BLACKMAN_HARRIS:
for(j=0; j<win_size; j++){
time_window[j] = 0.35875
- 0.48829 * cos( 2 pi j / (win_size-1))
+ 0.14128 * cos( 4 pi j / (win_size-1))
- 0.01168 * cos( 6 pi j / (win_size-1));
}
break;
}
} /* mk_window() */
Gast
#2512920
Ich hab' keine Ahnung wieso beim copy&paste aus "2 * pi *j" und so die Sternchen entfernt wurden. Dafür wird PI fett angezeigt. Abweichend vom Kommentar heißt der Buffer time_window und nicht window. Gruß, Gerhard http://www.mikrocontroller.net/articles/Formatierung_im_Forum mit c und /c oder pre /pre in eckigen Klammern müsste es gehen Antwort schreiben |
Anzeige
|
