脑电图信号的带宽

Question

我试着用Python语言对EEG信号进行FFT，然后根据带宽确定它是alpha信号还是beta信号。它看起来很好，但得到的曲线图与它们应该的完全不同，频率和震级值并不是我所期望的。如果有人愿意帮忙，代码如下：

from scipy.io import loadmat
import scipy
import numpy as np
from pylab import *
import matplotlib.pyplot as plt

eeg = loadmat("eeg_2013.mat");
eeg1=eeg['eeg1'][0]
eeg2=eeg['eeg2'][0]
fs = eeg['fs'][0][0]
fft1 = scipy.fft(eeg1)
f = np.linspace (fs,len(eeg1), len(eeg1), endpoint=False)
plt.figure(1)
plt.subplot(211)
plt.plot (f, abs (fft1))
plt.title ('Magnitude spectrum of the signal')
plt.xlabel ('Frequency (Hz)')
show()
plt.subplot(212)
fft2 = scipy.fft(eeg2)
f = np.linspace (fs,len(eeg2), len(eeg2), endpoint=False)
plt.plot (f, abs (fft2))
plt.title ('Magnitude spectrum of the signal')
plt.xlabel ('Frequency (Hz)')
show()

情节： 在

Answer 1

如果采样频率是fs，并且有N=len(eeg1)个样本，那么fft过程当然会返回一个N值的数组。其中第一个N/2对应于频率范围{}，频率的后半部分对应于镜像频率范围-fs/2..0。对于实际输入信号，镜像半部分只是正半部分的复共轭，因此在进一步的分析中可以忽略它（但在逆fft中则不能）。在

所以本质上，你应该格式化

f=linspace(0,N-1,N)*fs/N

编辑：或者更简单，只需对初始代码进行最小的更改

f = np.linspace (0,fs,len(eeg1), endpoint=False)

因此f的范围从0到fs之前，忽略输出中fft结果的后半部分：

plt.plot( f(0:N/2), abs( fft1(0:N/2) ) )

---

补充：您可以使用fftshift来交换两半，那么正确的频率范围是

f = np.linspace (-fs/2,fs/2,len(eeg1), endpoint=False)