function?W?=?csp（XY）

%?Common?spatial?patterns?for?spatial?filtering
%
%?X:?EEG?data?of?class?1?（channel?x?sample?point?x?trial）
%?Y:?EEG?data?of?class?2?（channel?x?sample?point?x?trial）
%?W:?the?colunms?of?projection?matrix?W?are?the?spatial?filters
%?d:?eigenvalues
%
%?yuzhang?2013.12.16?ECUST?China
%

%%?Compute?covariance?matrix?for?two?classes
Cx?=?zeros（size（X1）size（X1））;
Cy?=?zeros（size（Y1）size（Y1））;
for?i?=?1:size（X3）
????Cx?=?Cx?+?cov（X（::i）‘）;
end
for?i?=?1:size（Y3）
????Cy?=?Cy?+?cov（Y（::i）‘）;
end
Cx?=?Cx/size（X3）;?????%?covariance?matrix?averaged?on?all?trials
Cy?=?Cy/size（Y3）;
C=Cx+Cy;


%%?Solve?CSP+FC?joint?spatial?filtering
%[eigvec_Ceigval_C]=eig（C）;
%P=eigval_C^（-0.5）*eigvec_C‘;
%S1=P*Cx*P‘;
%[Veigval_S1]=eig（S1）;
%[sort_v

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1019??2020-09-20?05:42??csp.m