EDA365电子论坛网
标题:
CSP 共空间模式的 Matlab代码实现
[打印本页]
作者:
Colbie
时间:
2020-8-21 13:30
标题:
CSP 共空间模式的 Matlab代码实现
CSP 共空间模式的 Matlab代码实现
( I% X7 J% }$ ^- l- u. C3 G
. v, M5 v# b$ ~
function W = csp(X,Y)
% 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(X,1),size(X,1));
Cy = zeros(size(Y,1),size(Y,1));
for i = 1:size(X,3)
Cx = Cx + cov(X(:,:,i)');
end
for i = 1:size(Y,3)
Cy = Cy + cov(Y(:,:,i)');
end
Cx = Cx/size(X,3); % covariance matrix averaged on all trials
Cy = Cy/size(Y,3);
C=Cx+Cy;
%% Solve CSP+FC joint spatial filtering
%[eigvec_C,eigval_C]=eig(C);
%P=eigval_C^(-0.5)*eigvec_C';
%S1=P*Cx*P';
%[V,eigval_S1]=eig(S1);
%[sort_val,idx]=sort(diag(eigval_S1),'descend');
%final_eigvec=V(:,idx);
%W=(final_eigvec'*P)';
[W,D] = eig(Cx,Cx+Cy);
[d,idx] = sort(diag(D),'descend'); %eigvectors sorted via eigvalues descend
W = W(:,idx);
5 \' K7 c9 h0 L8 p% H: C
+ \2 \4 B. f6 c& y# Y/ N
作者:
BUFSFER
时间:
2020-8-21 14:47
CSP 共空间模式的 Matlab代码实现。
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2