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??共空間模式\