Compute normalized PSD in the Bernstein-Dirichlet parametrization.
qpsd(omega, v, w, k, beta_basis_k, epsilon = 1e-20)
See (5) in Choudhuri et al (2004)