小波分析(时频域分析)——python
本实验基于deap数据对eeg信号某个人的fp1通道信号进行不同频段分解
小波分解
(1)导入pywt以及numpy库
import numpy as np
import matplotlib.pyplot as plt
import pywt
(2)对EEg信号进行小波分解
channel=[1] #选择14个通道工情绪分类
subjectList=['01']
(3)定义小波分解函数
def wave_Processing(sub, channel):
meta = []
with open('data\s' + sub + '.dat', 'rb') as file:
subject = pickle.load(file, encoding='latin1')
for i in range(0, 40)://40个材料进行遍历 loop over 0-39 trails
data = subject["data"][i]
labels = subject["labels"][i]
for j in channel:
sampleData = data[j][:1024] #获取FP1通道的前8秒的数据
wavelet='db4' #选取的小波基函数
X = range(len(sampleData))
wave =pywt.wavedec(sampleData, wavelet, level=4)
#小波重构