Fftfile-of-EEG 可以用FFT频谱对脑电信号进行提取。我们可以利用提取出的各个波段脑电信号

可以用FFT频谱对脑电信号进行提取。我们可以利用提取出的各个波段脑电信号，来诊断一些脑部疾病或者对大脑组织的电活动及大脑的功能状态进行分析。 1.将实验测得的脑电数据文件转换为文本文件(已经过50Hz陷波)

function [deltathetaalphabeta]=Fftfil（data）
[rc]=size（data）;
fs=256;
dt=1/fs;
N=r;
%n=0:N-1; 
%f=n/（N*dt）;   %频率序列

    x=data;
    y1=fft（x）;      %对原信号作FFT变换
    y2=y1;
    y3=y2;
    y4=y3;
    
    f1=2;f2=4;f3=8;f4=13;f5=30;
    yy1=zeros（1length（y1））; %设置与y相同元素的数组
    yy2=yy1;
    yy3=yy2;
    yy4=yy3;
    
    %delta
    for m=0:N-1   %将频率落在该频率范围及其大于Nyquist频率的波滤去
        if（m/（N*dt）>=f1&&m/（N*dt）=（1/dt-f2）&&m/（N*dt）<（1/dt-f1）） 
            yy1（m+1）=y1（m+1）;       %置在此频率范围内的振动振幅不变
        else
            yy1（m+1）=0;     %其余频率范围的振动振幅为零
        end
    end
    
    %theta
    for m=0:N-1   
        if（m/（N*dt）>=f2&&m/（N*dt）=（1/dt-f3）&&m/（N*dt）<（1/dt-f2）） 
            yy2（m+1）=y2（m+1）;     
        else
            

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

     文件   10085317  2012-03-21 09:33  Fftfile of EEG\0661.txt

     文件       1485  2012-05-25 11:26  Fftfile of EEG\Fftfil.m

     文件       2058  2012-06-02 10:07  Fftfile of EEG\main.m

     目录          0  2013-08-26 21:57  Fftfile of EEG

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

             10088860                    4