csp 脑机接口中运动想象

脑机接口中运动想象，经典算法共空间模式CSP

clear;clc;
load（‘B0403T_4_40.mat‘）;
D1=size（data1）;D2=size（data2）;D3=size（data3）;L=2;N1=60;N2=0.5*D3-N1; 
%L:W取前L列向量。N1：训练样本对数。N2：测试样本对数。

%数据按左右运动想象分成12；按训练测试分成w，t
j=1;k=1;
for i=1:size（data3）
    if label（i）==1
        c1（j）=i;j=j+1;
    else
        c2（k）=i;k=k+1;
    end
end

for i=1:N1
    w1（::i）=data（::c1（i））;
    w2（::i）=data（::c2（i））;
end
for i=1:N2
    t1（::i）=data（::c1（N1+i））;
    t2（::i）=data（::c2（N1+i））;
end
%CSP学W
W = CSP_func（w1w2）;
W1=W（:[1:0.5*L end-0.5*L+1:end]）;
%特征提取
for i=1:N1
    f1（::i）=W1‘*w1（::i）;
    f2（::i）=W1‘*w2（::i）;
end
for i=1:N2
    g1（::i）=W1‘*t1（::i）;
    g2（::i）=W1‘*t2（::i）;
end
%对投影后的样本求方差
ff1=[];ff2=[];
for i=1:N1
    ff1=[ff1 var（f1（::i）02）];
    ff2=[ff2 var（f2（::i）0

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1565  2015-05-16 12:58  共空间模式\CSP.m
     文件        1020  2015-04-27 22:26  共空间模式\CSP_func.m
     目录           0  2015-09-19 16:06  共空间模式\