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）;??????%對(duì)原信號(hào)作FFT變換
????y2=y1;
????y3=y2;
????y4=y3;
????
????f1=2;f2=4;f3=8;f4=13;f5=30;
????yy1=zeros（1length（y1））;?%設(shè)置與y相同元素的數(shù)組
????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）;???????%置在此頻率范圍內(nèi)的振動(dòng)振幅不變
????????else
????????????yy1（m+1）=0;?????%其余頻率范圍的振動(dòng)振幅為零
????????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
????????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???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