資源簡介
可以用FFT頻譜對腦電信號進行提取。我們可以利用提取出的各個波段腦電信號，來診斷一些腦部疾病或者對大腦組織的電活動及大腦的功能狀態(tài)進行分析。
1.將實驗測得的腦電數(shù)據(jù)文件轉(zhuǎn)換為文本文件(已經(jīng)過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））;?%設(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)的振動振幅不變
????????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
評論
共有 條評論