根据列表数据计算时频谱并显示1Hz~100Hz的频谱的python代码

以下是计算时频谱并显示1Hz~100Hz的频谱的Python代码：

import numpy as np
import matplotlib.pyplot as plt

# 生成随机信号
fs = 1000  # 采样率
t = np.arange(0, 1, 1/fs)
x = np.sin(2*np.pi*10*t) + np.sin(2*np.pi*20*t) + np.sin(2*np.pi*50*t)

# 计算时频谱
N = len(x)
f = np.arange(0, fs/2, fs/N)
X = np.fft.fft(x)
X = X[:N//2]
Pxx = np.abs(X)**2 / (fs*N)

# 显示1Hz~100Hz的频谱
plt.plot(f, Pxx)
plt.xlim([1, 100])
plt.xlabel('Frequency (Hz)')
plt.ylabel('Power Spectral Density (V^2/Hz)')
plt.show()

希望能够帮到您！