Module rpps.utils.psd

Functions

def psd(samps, Fs=1, vbw_hz=None)

Expand source code

def psd(samps, Fs=1, vbw_hz=None):
    y = np.abs(np.fft.fft(samps))**2 / (len(samps)*Fs)
    y = np.fft.fftshift(10.0*np.log10(y))
    if vbw_hz is not None:
        smooth = vbw_calc(len(samps), Fs, vbw_hz)
        y = vbw(y, smooth)
    return y

def vbw(samps, smooth)

Expand source code

def vbw(samps, smooth):
    return gaussian.gaussian(samps, gaussian.fwhm2sigma(smooth))

def vbw_calc(N, Fs, vbw_hz)

Expand source code

def vbw_calc(N, Fs, vbw_hz):
    return vbw_hz/(Fs / N)