Kaiser window

The Kaiser window, also known as the Kaiser–Bessel window, was developed by James Kaiser at Bell Laboratories. It is a one-parameter family of window functions used in finite impulse response filter design and spectral analysis. The Kaiser window approximates the DPSS window which maximizes the energy concentration in the main lobe^[1] but which is difficult to compute.^[2]

Definition

The Kaiser window and its Fourier transform are given by:

w_{0}(x)\triangleq \left\{{\begin{array}{ccl}{\frac {I_{0}\left[\pi \alpha {\sqrt {1-\left(2x/L\right)^{2}}}\right]}{I_{0}[\pi \alpha ]}},\quad &\left|x\right|\leq L/2\\0,\quad &\left|x\right|>L/2\end{array}}\right\}\quad {\stackrel {\mathcal {F}}{\Longleftrightarrow }}\quad {\frac {L\cdot \sinh \left(\pi \alpha {\sqrt {1-\left(Lf\right/\alpha )^{2}}}\right)}{I_{0}(\pi \alpha )\cdot \pi \alpha {\sqrt {1-\left(Lf\right/\alpha )^{2}}}}}.

^[3]

Fourier transforms of two Kaiser windows

where:

$I 0$ is the zeroth-order modified Bessel function of the first kind,
$L$ is the window duration, and
$α$ is a non-negative real number that determines the shape of the window. In the frequency domain, it determines the trade-off between main-lobe width and side lobe level, which is a central decision in window design.
Sometimes the Kaiser window is parametrized by $β$ , where $β = πα$ .

For digital signal processing, the function can be sampled symmetrically as:

w[n]=w_{0}\left({\tfrac {L}{N}}(n-N/2)\right)={\frac {I_{0}\left[\pi \alpha {\sqrt {1-\left({\frac {2n}{N}}-1\right)^{2}}}\right]}{I_{0}[\pi \alpha ]}},\quad 0\leq n\leq N,

where the length of the window is $N+1,$ and N can be even or odd. (see Window_function#A_list_of_window_functions)

In the Fourier transform, the first null after the main lobe occurs at $f={\tfrac {\sqrt {1+\alpha ^{2}}}{L}},$ which is just ${\sqrt {1+\alpha ^{2}}}$ in units of N (DFT "bins"). As α increases, the main lobe increases in width, and the side lobes decrease in amplitude. $α$ = 0 corresponds to a rectangular window. For large $α,$ the shape of the Kaiser window (in both time and frequency domain) tends to a Gaussian curve. The Kaiser window is nearly optimal in the sense of its peak's concentration around frequency 0.^[4]

Kaiser–Bessel-derived (KBD) window

A related window function is the Kaiser–Bessel-derived (KBD) window, which is designed to be suitable for use with the modified discrete cosine transform (MDCT). The KBD window function is defined in terms of the Kaiser window of length N+1, by the formula:

d_{n}={\begin{cases}{\sqrt {\frac {\sum _{i=0}^{n}w[i]}{\sum _{i=0}^{N}w[i]}}}&{\mbox{if }}0\leq n<N\\{\sqrt {\frac {\sum _{i=0}^{2N-1-n}w[i]}{\sum _{i=0}^{N}w[i]}}}&{\mbox{if }}N\leq n\leq 2N-1\\0&{\mbox{otherwise}}.\\\end{cases}}

This defines a window of length 2N, where by construction d_n satisfies the Princen-Bradley condition for the MDCT (using the fact that $w N - n = w n$ ): $d n 2 + (d n + N) 2 = 1$ (interpreting n and n + N modulo 2N). The KBD window is also symmetric in the proper manner for the MDCT: d_n = d_2N−1−n.

Applications

The KBD window is used in the Advanced Audio Coding digital audio format.

References

^ "Slepian or DPSS Window". ccrma.stanford.edu. Retrieved 2016-04-13.
^ Oppenheim, A. V.; Schafer, R. W. (2009). Discrete-time signal processing. Upper Saddle River, N.J.: Prentice Hall. p. 541. ISBN 9780131988422.
^ Harris, Fredric J. (Jan 1978). "On the use of Windows for Harmonic Analysis with the Discrete Fourier Transform" (PDF). Proceedings of the IEEE. 66 (1): 51–83. CiteSeerX 10.1.1.649.9880. doi:10.1109/PROC.1978.10837.
^ Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). "7.2". Discrete-time signal processing (2nd ed.). Upper Saddle River, N.J.: Prentice Hall. p. 474. ISBN 0-13-754920-2. a near-optimal window could be formed using the zeroth-order modified Bessel function of the first kind Also available at https://d1.amobbs.com/bbs_upload782111/files_24/ourdev_523225.pdf

Definition

Kaiser–Bessel-derived (KBD) window

Applications

References

Further reading