EEG快速傅里叶变换PYTHON代码

以下是一个基于Python的EEG快速傅里叶变换的示例代码：

import numpy as np
import matplotlib.pyplot as plt

# 生成模拟EEG信号
N = 256
T = 1 / 256
t = np.linspace(0, N*T, N)
y = np.sin(5*np.pi*t) + np.sin(10*np.pi*t) + np.sin(20*np.pi*t)

# 进行傅里叶变换
yf = np.fft.fft(y)
xf = np.linspace(0, 1/(2*T), int(N/2))

# 绘制频谱图
plt.plot(xf, 2.0/N * np.abs(yf[:N//2]))
plt.xlabel('Frequency')
plt.ylabel('Amplitude')
plt.show()

这个代码生成一个包含三个正弦波的模拟EEG信号，并使用快速傅里叶变换来计算其频谱。最后，它绘制了频谱图，显示了信号中每个频率的振幅。