脑电波MATLAB分析代码

大学里做过脑波控制，参考的一个MATLAB代码，觉得挺好的就分享给大家。

%fft变换
close all;
clc;
clear all;

x1=dlmread（‘5.txt‘）;
x=x1;%带入原始信号
fs=200;%设置采样率
N=length（x）;%计算变换点数，直接影响曲线的平滑度
y=fft（xN）;%进行fft变换

mag3=（abs（y））*2/N;                
                           
f=（0:N-1）*fs/N;
figure;
plot（fmag3）;
 xlim（[10 20]）%将x轴范围限定在10-20Hz
xlabel（‘频率/Hz‘）;ylabel（‘幅值/db‘）;title（‘fft变换‘）;


for i=1:N
    if（i-1）*fs/N<10
       y（i）=0.0001*y（i）;
    end
     if（i-1）*fs/N>190
       y（i）=0.0001*y（i）;
     end
    
      if（（i-1）*fs/N<180）&&（（i-1）*fs/N>20）
       y（i）=0.0001*y（i）;
     end
     
     
end
%定位信号频率值
[ab]=max（abs（y））;

%逐次轮选数据进行测试，选择包含信号频率的
data=zeros（1256）;
for i=1000:N-1000;
    x0=x（i:i-1+256）;
    y0=fft（x0N）;
    
    
    for i=1:N
     if（i-1）*fs/N<10
          y0（i）=0.0001*y（i）;
     end
     
     if（i-1）*fs/N>190
       y0（i）=0.0001*y（i）;
     end
    
     if（（i-1）*fs/N<180）&&（（i-1）*fs/N>20）
       y0（i）=0.0001*y（i）;
     end
   
    end
     
    [a0b0]=max（abs（y0））;
    
    
    if b0==b
        if y0（b0-1）<0.05*y0（b0）
            if y0（b0+1）<0.05*y0（b0）
            data=x0;
            break;
            end
        end
    end
    
end
    
%  选择的数据存为txt

fid=fopen（‘data_256.txt‘‘wt‘）;
fprintf（fid‘%g\t‘data）;
fclose（fid）;

%查看选择的数据的FFT图
yn=data;
 N=256;
y=fft（ynN）;



mag3=（abs（y））*2/N;%求幅值
                            
f=（0:N-1）*fs/N;%进行对应的频率转换


figure;
plot（fmag3）;
 xlim（[10 20]）%将x轴范围限定在10-20Hz
xlabel（‘频率/Hz‘）;ylabel（‘幅值/db‘）;title（‘fft变换‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      22694  2013-09-25 16:22  原始信号\1.txt

     文件      21440  2014-10-21 21:49  原始信号\10.txt

     文件        807  2016-01-15 18:59  原始信号\13hz.txt

     文件      14755  2014-10-21 21:44  原始信号\2.txt

     文件      21869  2014-10-21 21:45  原始信号\3.txt

     文件      16548  2014-10-27 11:41  原始信号\4.txt

     文件      18867  2014-10-21 21:46  原始信号\5.txt

     文件       2685  2016-01-15 19:05  原始信号\5_256.txt

     文件      12664  2014-10-27 20:41  原始信号\6.txt

     文件      17005  2014-10-21 21:47  原始信号\7.txt

     文件      17404  2014-10-21 21:48  原始信号\8.txt

     文件       2606  2016-01-15 18:37  原始信号\8_256.txt

     文件      16331  2014-10-27 11:49  原始信号\9.txt

     文件        841  2016-01-15 22:01  原始信号\data_256.txt

     文件       1650  2016-01-15 21:27  原始信号\data_pro.asv

     文件       1651  2016-01-15 22:01  原始信号\data_pro.m

     文件       1148  2016-01-15 18:46  原始信号\g_fft.asv

    I.A....      1213  2016-01-15 19:06  原始信号\g_fft.m

     目录          0  2016-01-15 21:36  原始信号

----------- ---------  ---------- -----  ----

               192178                    19