Thursday, 21 July 2016

Measure scalloping loss of window function


Is there an easy/general way to calculate this? For instance, the rectangular window has a maximum scalloping loss of 0.3634, which can be derived in frequency domain from a normalized sinc(0.5).


The only way, given a precomputed window function, I can think of, is to apply the window function to a signal whose period is 2pi*(k + 0.5)/N and amplitude A, calculate a N FFT and measure the bin k. Bin k should have a magnitude of A * N * f, where f is a factor giving the worst-case scalloping loss, irrespective of N/sampling size, right?



Answer



Harris defines the scalloping loss as the ratio of coherent gain for a tone located half a bin from a DFT sample to the coherent gain at a sample (see the paper for the definitions) SL=|nw(nT)ejπNn|nw(nT)=W(ω1/2)W(0)

Where w is the window and W is the window DFT


ω1/2=πNT


Guess the value you are after is really 1SL



EDIT: I believe the left hand side of the equation is enough to answer @Shaggi's question, but to clarify a bit the notation of the rest, Harris defines DFTs as functions of ωk, so


W(ωk)=W(2πkNT)==nw(nT)ejωknT=nw(nT)ej2πkn/NW(ω1/2)=nw(nT)ejπn/NW(0)=nw(nT)



F. J. Harris - On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform



No comments:

Post a Comment

readings - Appending 内 to a company name is read ない or うち?

For example, if I say マイクロソフト内のパートナーシップは強いです, is the 内 here read as うち or ない? Answer 「内」 in the form: 「Proper Noun + 内」 is always read 「ない...