www.pudn.com > cubica.rar > ly.m

 
 
%%%%%%%%%%%%%%%%%%%%   FastICA    %%%%%%%%%%%%%%%% 
m = 100000;%number of observations 
n=2; 
 
 
    x=[]; 
	y=wavread('EXP3_XN1'); 
    x(1,:)=(y(1:m))'; 
    %x(1,:)=modulate(x(1,:),10000,30000,'am',0.3); 
    z=wavread('EXP3_XN2'); 
    x(2,:)=(z(1:m))'; 
    %x(2,:)=modulate(x(2,:),10000,30000,'am',0.3); 
 
 
 
x = inv(diag(std(x'))) * x;  %¹éÒ»»¯ 
 
 
[R,y]=cubica34a(x); 
 
u=R*x;	 
 
subplot(2,2,1);plot(x(1,:)); 
title('mixedsig(1)'); 
 
subplot(2,2,2);plot(x(2,:)); 
title('mixedsig(2)'); 
 
subplot(2,2,3);plot(u(1,:)); 
title('IC(1)'); 
 
subplot(2,2,4);plot(u(2,:)); 
title('IC(2)');