资源描述：
《基于matlab的脑电波信号处理》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、做脑电波信号处理滴嘿嘿。。MatlabaddictedCodes%FEATUREEXTRACTERfunction[features]=EEGfeaturetrainmod(filename,m)a=4;b=7;d=12;e=30;signals=0;forindex=1:9;%readinthefirsttenEEGdatabecausethefilesarenumberedasha11test01ratherthanha11test1.s=[filename'0'num2str(index)'.dat'];signal=tread_wfdb(s);ifsignals==0;signals
2、=signal;elsesignals=[signalssignal];endendforindex=10:1:m/2;%readintherestoftheEEGtrainingdatas=[filenamenum2str(index)'.dat'];signal=tread_wfdb(s);ifsignals==0;signals=signal;elsesignals=[signalssignal];endend%%%%%modificationjustforvaryingthetrainingtestingratio------forindex=25:1:25+m/2;%readint
3、herestoftheEEGtrainingdatas=[filenamenum2str(index)'.dat'];signal=tread_wfdb(s);ifsignals==0;signals=signal;elsesignals=[signalssignal];endend%%%%%endofmodificationjustforvaryingthetrainingtestingratio-----forl=1:m%exratingfeatures(powerofeachkindofEEGwaveforms)[Pxx,f]=pwelch(signals(:,l)-mean(sign
4、als(:,l)),[],[],[],200);%relativepowerfdelta(l)=sum(Pxx(find(fa)));falpha(l)=sum(Pxx(find(fb)));fbeta(l)=sum(Pxx(find(fd)));fgama(l)=sum(Pxx(find(f>e)));%gamawaveincludedforadditionalworkendfeatures=[fdelta;ftheta;falpha;fbetaa;fgama];features=features';
5、end%CLASSIFIER%(Hasthreesimilarclassificationmodifation:EEGclassification,EEGclassificationmodandEEGclassificationmod1savedandusedintherunningfileforadditionalworks)function[class,err,classall,errall]=EEGclassification(trainfilename,m,testfilename,n,p,q)%p-waveform1,q-waveformtwoo–waveformthree%1-d
6、elta2-theta3-alpha4–beta5-Gamma[featurestrain]=EEGfeature(trainfilename,m);%modificationtoEEGfeaturemodfunctionforvaryingtestingtrainingratio[featurestest]=EEGfeature(testfilename,n);training=[featurestrain(:,p)featurestrain(:,q)];%modifyhowmanyfeaturestoextractheresample=[featurestest(:,p)features
7、test(:,q)];group=[ones(m/2,1);2*ones(m/2,1)];traininga=featurestrain;samplea=featurestest;[class,err,POSTERIOR,logp,coeff]=classify(sample,training,group,'quadratic');%'mahalanobis','quadratic','linear'asde
