Python: FM演示实现

Question

我正在考虑分析一些特定值的时间序列，就像它是一个调频信号一样。

我正在寻找一个调频解调器的Python实现。我知道Matlab和Octave中有一个解调器函数；对于Python，我找到了这个FreqDemod包，但它似乎没有做我想做的事情。

我们将非常感谢您的帮助。

Answer 1

这是一个对复杂样本进行FM解调的Python函数。

def fm_demod(x, df=1.0, fc=0.0):
    ''' Perform FM demodulation of complex carrier.

    Args:
        x (array):  FM modulated complex carrier.
        df (float): Normalized frequency deviation [Hz/V].
        fc (float): Normalized carrier frequency.

    Returns:
        Array of real modulating signal.
    '''

    # Remove carrier.
    n = sp.arange(len(x))
    rx = x*sp.exp(-1j*2*sp.pi*fc*n)

    # Extract phase of carrier.
    phi = sp.arctan2(sp.imag(rx), sp.real(rx))

    # Calculate frequency from phase.
    y = sp.diff(sp.unwrap(phi)/(2*sp.pi*df))

    return y