matlab脑电打码,matlab 脑电信号特征提取程序

编写一程序。。驾驶员清醒和疲劳状态时脑电信号特征提取，只需要把脑电信号转换成功率谱

脑电信号已采集

t=0.001:0.001:1;

x=load('D:\EEGdata_data');     %读取文件

y=load('D:EEGdata_newdata\');

xx={};           %每个导联的数据存储

yy={};

n=1000;           %数据数目

sc=7;             %小波包的分解尺度

for i=1:1:8       %导联的数据分离

xx{i}=x(:,i);

yy{i}=y(:,i);

end

for i=1:1:8

%画出原始信号图像

figure

subplot(2,2,1)

plot(t,xx{i})

axis([min(t) max(t) 1.1*floor(min(xx{i})) 1.1*ceil(max(xx{i}))])

title('rest close 原始信号')

ylabel('幅值')

subplot(2,2,2)

plot(t,yy{i})

axis([min(t) max(t) 1.1*floor(min(yy{i})) 1.1*ceil(max(yy{i}))])

title('audio close 原始信号')

ylabel('幅值')

%fft_原始信号的频谱分析

xx1=fft(xx{i},n);

pxx