eeg_complexic.m

计算脑电信号非线性参数

%本程序实现脑电信号非线性特征参数--复杂度的计算%
%林宛华  中南大学生物医学工程医学图像与信号处理实验室   %
%E-mail: [e-mail]%
%2008年3月10日%
function Kc=EEG_Complexic(data)
c=1;
EEG=EEG_data0;
BIN_EEG=zeros(1,length(EEG));%初始化粗粒度序列
SUM_EEG=sum(EEG(1:length(EEG)));
MEAN_EEG=SUM_EEG/length(EEG);
for i=1:length(EEG)  %脑电信号粗粒化
      if (EEG(i)>=MEAN_EEG)
          BIN_EEG(i)=1;
      elseif(EEG(i)<MEAN_EEG)
          BIN_EEG(i)=0;
      end
end
S=[BIN_EEG(1)];Q=[BIN_EEG(2)];
%for i=1:length(BIN_EEG)-2
while((length(S)+length(Q))<length(BIN_EEG))
SQ=[S Q];
SQD=SQ(1:length(SQ)-1);
isbelong = 0;
for r=1:length(SQD)-length(Q)+1  %判断Q是否是SQD的子串
   Sr=SQD(r:r+length(Q)-1);
       if(Sr==Q) 
          isbelong = 1; %Q是否是SQD的子串
       end
end
if(isbelong == 1)
   Q=[Q BIN_EEG(length(SQ)+1)];                                
 elseif(isbelong == 0)
   c=c+1;
   S=[S Q];
   Q=BIN_EEG(length(S)+1);  
end
end
 %最后一组Q是否是新子串，不是c+1
SQ=[S Q];
SQD=SQ(1:length(SQ)-1);
isbelong = 0;
for r=1:length(SQD)-length(Q)+1  %判断Q是否是SQD的子串
   Sr=SQD(r:r+length(Q)-1);
       if(Sr==Q) 
          isbelong = 1; %Q是否是SQD的子串
       end
end
if(isbelong == 0)
   c=c+1;
end
b=length(BIN_EEG)/(log2(length(BIN_EEG)));
Kc=c/b;
return