如何在不使用频谱图的情况下获取频率与时间的关系?

时间:2019-06-07 14:08:57

标签: python signal-processing fft frequency-analysis spectrogram

如何获得频率随时间变化的信号?      我在下面发布我的代码,以计算正弦波的fft,有没有使用频谱图的解决方案?     有什么帮助吗?

import matplotlib.pyplot as plt
import numpy as np
f = 15 #fréquence
f1 = 2
f_s = 2000 # Fréq d'éch
t = np.linspace(0, 10, 10 * f_s, endpoint=False) # axe des temps
fct = np.exp(1/(t+1))
y = np.sin(f1 * 2 * np.pi * t )
x = np.sin(f * 2 * np.pi * t + y) # sinusoide avec une fréquence variable
fig, figure1 = plt.subplots()
figure1.plot(t, x)
figure1.set_xlim(0, 10)
from scipy import fftpack
X = fftpack.rfft(x)
# création de l'axe des fréquences
freqs = fftpack.rfftfreq(len(x)) * f_s
# Génération d'une figure
fig, figure = plt.subplots()
figure.plot(freqs, np.abs(X))
figure.set_xlim(0, 32)

0 个答案:

没有答案