pub fn kaiser<T: Float>(n: usize, beta: T) -> Vec<T>
Generate Kaiser window.
w[n] = I₀(β * sqrt(1 - ((n - (N-1)/2) / ((N-1)/2))²)) / I₀(β)
where I₀ is the modified Bessel function of the first kind.