频谱分析代码python

可以直接回答，以下是示例代码：

import numpy as np
import matplotlib.pyplot as plt
from scipy.fftpack import fft

# 生成信号序列
t = np.linspace(0, 1, 1000)
y = np.sin(2*np.pi*50*t) + np.sin(2*np.pi*120*t)

# 频率谱分析
N = len(y)
fft_res = fft(y) / N
fft_res_abs = np.abs(fft_res[range(int(N/2))])

# 绘制频率谱图
freq = np.linspace(0, 1/2, int(N/2))
plt.plot(freq, fft_res_abs)
plt.show()

注意，这个代码仅仅是示例，请根据具体需求进行修改和适配。