Python:调频解调实现
Python: FM Demod Implementation
我正在考虑分析一些特定值的时间序列,因为它是调频信号。
我正在寻找 Python 调频解调器的实现。
我知道 Matlab 和 Octave 中有解调器功能;对于 Python 我找到了这个 FreqDemod 包,但它似乎没有做我想做的事情。
将不胜感激。
这是一个 Python 函数,可以对复杂样本进行 FM 解调。
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
我正在考虑分析一些特定值的时间序列,因为它是调频信号。
我正在寻找 Python 调频解调器的实现。 我知道 Matlab 和 Octave 中有解调器功能;对于 Python 我找到了这个 FreqDemod 包,但它似乎没有做我想做的事情。
将不胜感激。
这是一个 Python 函数,可以对复杂样本进行 FM 解调。
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