JADE α波単離不可？

安静閉眼時の自発脳磁図波形にJADEを用い、α波帯域の成分の単離を試みました。
以下のデータを用います。

megをplanar gradiometer 204chの約10秒の安静閉眼時自発脳磁場データ(SSP処理済)、 timeを時間、sfreqを標本化周波数600.615Hzです。
%% spectrum of MEG2112 & MEG2113
%% signal processing toolbox was employed!
N=2^8;
[B1,F,T]=specgram(meg(158,:),N,sfreq,N,N/2); %MEG2112
B2=specgram(meg(157,:),N,sfreq,N,N/2); %MEG2113
f60=find(F<80);f=[min(F),max(F(f60))];t=[min(T),max(T)];
figure;set(gcf,'color',[1,1,1]);
subplot(211);imagesc(t,f,abs(B1(f60,:))/N^0.5,[0,200]);colorbar;
subplot(212);imagesc(t,f,abs(B2(f60,:))/N^0.5,[0,200]);colorbar;

後頭部のセンサーMEG2112とMEG2113の時間・周波数分布を見ています。 横軸が時間、縦軸が周波数です。単位はfT/cm/√Hzです。 10Hzのところにα波が見られます。60Hzの高値は西日本の交流電源の周波数を意味します。
%% preparation for topogrpaphy
load('.../vv_lout.mat');
W=120;x=vv_mag(:,2);y=vv_mag(:,3);
xlin=linspace(min(x),max(x),W);
ylin=linspace(min(y),max(y),W);
[X,Y]=meshgrid(xlin,ylin);
%% apply JADE
tic;B=jadeR(meg,8);toc
0.9840秒で計算が終わりました。
IC=B*meg;AA=pinv(B);
A=sum(AA.^2,1);P=diag(A)*IC;
V=AA./(ones(size(AA,1),1)*A);
figure;set(gcf,'color',[1,1,1]);
set(gcf,'renderer','zbuffer');
nn=8;for n=1:nn;...
subplot('position',[0,1-n/nn,0.2,1/nn]);...
z=(V(1:2:204,n).^2+V(2:2:204,n).^2).^0.5;...
Z=griddata(x,y,z,X,Y,'linear');...
h=surf(X,Y,Z);axis off;set(h,'edgecolor','none');...
axis tight;daspect([1,1,0.01]);view([0,90]);...
subplot('position',[0.25,1-n/nn+0.02,0.7,1/nn-0.06]);...
plot(time,P(n,:));axis tight;grid on;end;

各独立成分の波形です。第１独立成分は交流電源の成分、 第２・第３独立成分がα波の成分ぽく見えます。 検討のため各独立成分波形の周波数変化を見ていきます。
figure;set(gcf,'color',[1,1,1]);
set(gcf,'renderer','zbuffer');
nn=8;for n=1:nn;...
subplot('position',[0,1-n/nn,0.2,1/nn]);...
z=(V(1:2:204,n).^2+V(2:2:204,n).^2).^0.5;...
Z=griddata(x,y,z,X,Y,'linear');...
h=surf(X,Y,Z);axis off;set(h,'edgecolor','none');...
axis tight;daspect([1,1,0.01]);view([0,90]);...
subplot('position',[0.25,1-n/nn+0.02,0.7,1/nn-0.06]);...
BP=specgram(P(n,:),N,sfreq,N,N/2);...
imagesc(t,f,abs(BP(f60,:))/N^0.5,[0,2e+5]);colorbar;end;

交流電源だけの成分・α波だけの成分が抽出されているというわけではないようです。

時間を長く取りすぎたのが原因かもしれません。そこで2秒間のデータで検討することにします。
meg2=meg(:,1:fix(sfreq*2));
time2=time(1:fix(sfreq*2));
figure;set(gcf,'color',[1,1,1]);
subplot(211);plot(time2,meg2(158,:));grid on;axis tight;
subplot(212);plot(time2,meg2(157,:));grid on;axis tight;

B=jadeR(meg2,8);
IC=B*meg2;AA=pinv(B);
A=sum(AA.^2,1);P=diag(A)*IC;
V=AA./(ones(size(AA,1),1)*A);
figure;set(gcf,'color',[1,1,1]);
set(gcf,'renderer','zbuffer');
nn=8;for n=1:nn;...
subplot('position',[0,1-n/nn,0.2,1/nn]);...
z=(V(1:2:204,n).^2+V(2:2:204,n).^2).^0.5;...
Z=griddata(x,y,z,X,Y,'linear');...
h=surf(X,Y,Z);axis off;set(h,'edgecolor','none');...
axis tight;daspect([1,1,0.01]);view([0,90]);...
subplot('position',[0.25,1-n/nn+0.02,0.7,1/nn-0.06]);...
plot(time2,P(n,:));axis tight;grid on;end;

T=0:N/2/sfreq:2;
t=[min(T),max(T(find(T<2-N/sfreq)))];
figure;set(gcf,'color',[1,1,1]);
set(gcf,'renderer','zbuffer');
nn=8;for n=1:nn;...
subplot('position',[0,1-n/nn,0.2,1/nn]);...
z=(V(1:2:204,n).^2+V(2:2:204,n).^2).^0.5;...
Z=griddata(x,y,z,X,Y,'linear');...
h=surf(X,Y,Z);axis off;set(h,'edgecolor','none');...
axis tight;daspect([1,1,0.01]);view([0,90]);...
subplot('position',[0.25,1-n/nn+0.02,0.7,1/nn-0.06]);...
BP=specgram(P(n,:),N,sfreq,N,N/8*7);...
imagesc(t,f,abs(BP(f60,:))/N^0.5,[0,2e+5]);colorbar;end;

独立成分分析で単離した「α波」には低周波の成分が混在すると主張して認知されるだろうか？
・・・今のところα波と主張するには素直に周波数帯域通過フィルターを使った方が無難のような気がします。