ICALAB 変数

先ほどsaveしたtest_ica.matには以下の変数が保存されていました。

SはPeigSで計測された信号、 XはSを混合行列HとSの積、 Wは復元行列、 Yは各独立成分の波形（Source）で、Y=W・X, G=W・Hです。
混合行列が単位行列のとき、即ちH=Iのとき、G=Wとなります。 途中でMixing matrix Hのところはいじっていないので、H,G,Xは本来不要です。 このうちW(&S?)以外は、保存しなくてもMATLABのWorkspace上の変数として存在しています。 但しX=x,Y=yですが。

topopgraphy表示します。
load('.../vv_lout.mat');
Wt=120;xt=vv_mag(:,2);yt=vv_mag(:,3);
xlin=linspace(min(xt),max(xt),Wt);
ylin=linspace(min(yt),max(yt),Wt);
[Xt,Yt]=meshgrid(xlin,ylin);
AA=Veig(:,1:size(W,1))*pinv(W);
At=sum(AA.^2,1);
Pt=diag(At)*Y;
Vt=AA./(ones(size(AA,1),1)*At);
figure;set(gcf,'color',[1,1,1]);
set(gcf,'renderer','zbuffer');
nn=20;
for n=1:nn;...
zt=(Vt(1:2:204,n).^2+Vt(2:2:204,n).^2).^0.5;...
Zt=griddata(xt,yt,zt,Xt,Yt,'linear');...
ha=subplot(nn,2,n*2-1);...
position=get(ha,'position');...
position=[0,position(2),1/nn,position(4)];...
set(ha,'position',position);...
hh=surf(Xt,Yt,Zt);axis off;axis tight;...
set(hh,'edgecolor','none');...
daspect([1,1,1]);view([0,90]);...
hb=subplot(nn,2,n*2);...
plot(time,Pt(n,:));axis tight;grid on;...
position=get(hb,'position');...
position=[1/nn+0.05,position(2),0.9-1/nn,position(4)];...
set(hb,'position',position);...
end;

うまく独立成分に分けられているかどうかよくわかりません。