code: def sinegauss_PSD(f, t0, f0, Q):
      """Compute the PSD of the sine-gaussian function at frequency f
  
      .. math::
  
         PSD(f) = |H(f)|^2 + |H(-f)|^2
      """
      a = (f0 * 1.0 / Q) ** 2
      Pf = np.pi / a * np.exp(-2 * np.pi**2 * (f - f0) ** 2 / a)
      Pmf = np.pi / a * np.exp(-2 * np.pi**2 * (-f - f0) ** 2 / a)
      return Pf + Pmf
  
file_path: astroML/astroML/fourier.py
imports: import numpy
label_desc: Use image node_id 6 for calling a global function with example usage: sinegauss_PSD(f, t0, f0, Q) and returns: unknown