Kaiser window

The Kaiser window for several values of its parameter

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<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> but which is difficult to compute.<ref>Template:Cite book</ref>

DefinitionEdit

The Kaiser window and its Fourier transform are given by:

<math>

w_0(x) \triangleq \left\{ \begin{array}{ccl} \tfrac{1}{L}\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{\sin\bigg(\sqrt{(\pi L f)^2-(\pi \alpha)^2}\bigg)} {I_0(\pi \alpha)\cdot \sqrt{(\pi L f)^2-(\pi \alpha)^2}}, </math> <ref>Template:Cite journal</ref>Template:Efn-ua.</math> }}

File:Kaiser-Window-Spectra.svg

Fourier transforms of two Kaiser windows

where:

Template:Math is the zeroth-order modified Bessel function of the first kind,
Template:Mvar is the window duration, and
Template:Math 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 Template:Math, where Template:Math.

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

<math>w[n] = L\cdot 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,</math>

where the length of the window is <math>N+1,</math> and N can be even or odd. (see A list of window functions)

In the Fourier transform, the first null after the main lobe occurs at <math>f = \tfrac{\sqrt{1+\alpha^2}}{L},</math> which is just <math>\sqrt{1+\alpha^2}</math> in units of N (DFT "bins"). As α increases, the main lobe increases in width, and the side lobes decrease in amplitude. Template:Math = 0 corresponds to a rectangular window. For large Template:Math 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 <math>0.</math><ref name=Oppenheim> Template:Cite book </ref>

Kaiser–Bessel-derived (KBD) windowEdit

File:Kbd-window.svg

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:

<math>

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} </math>

This defines a window of length 2N, where by construction d_n satisfies the Princen-Bradley condition for the MDCT (using the fact that Template:Math): Template:Math (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.

ApplicationsEdit

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

NotesEdit

Template:Notelist-ua

ReferencesEdit