www.pudn.com > icaML.zip > script_demo.m

% Demo script 
% version 1.0 
 
 
% Generate data 
K_true=6 % number of sources 
interval=[2:10]; % Number of IC components interval to search 
 
% Generate mixed signals from artificial super Gaussian sources and mixing matrix 
Atrue=randn(20,K_true); 
Strue=sinh(randn(K_true,1000))/.1; 
X=Atrue*Strue; 
 
 
 
% - 1th Part - Finding best model wrt. number of components using BIC 
 
%BIC 
P=icaML_bic(X,interval,1); 
 
[most_prop,most_prop_k]=max(P); 
most_prop 
most_prop_k=interval(most_prop_k) 
 
% Draw 
figure(1) 
clf 
bar(interval,P); 
drawnow; 
 
 
 
% - 2nd Part - finding ICA components using BIC estimat of K (number of source components) 
 
% ICA 
[S,Aica,U] = icaML( X , most_prop_k , [] , 1); 
A = U(:,1:most_prop_k) * Aica; 
 
% Draw first two IC components in scatterplot 
figure(2) 
clf 
plot(S(1,:),S(2,:),'.') 
drawnow;