function [sig,mixedsig]=demosig();%% function [sig,mixedsig]=demosig();% % Returns artificially generated test signals, sig, and mixed% signals, mixedsig. Signals are row vectors of% matrices. Input mixedsig to FastICA to see how it works.% @(#)$Id: demosig.m,v 1.2 2003/04/05 14:23:57 jarmo Exp $%create source signals (independent components)N=500; %data sizev=[0:N-1];sig=[];sig(1,:)=sin(v/2); %sinusoidsig(2,:)=((rem(v,23)-11)/9).^5; %funny curvesig(3,:)=((rem(v,27)-13)/9); %saw-toothsig(4,:)=((rand(1,N)<.5)*2-1).*log(rand(1,N)); %impulsive noisefor t=1:4sig(t,:)=sig(t,:)/std(sig(t,:));end%remove mean (not really necessary)[sig mean]=remmean(sig);%create mixturesAorig=rand(size(sig,1));mixedsig=(Aorig*sig);subplot(411):plot(mixedsig(1,:));subplot(412):plot(mixedsig(2,:));subplot(413):plot(mixedsig(3,:));subplot(414):plot(mixedsig(4,:));gza=FASTICA(mixedsig);subplot(411):plot(gza(1,:));subplot(412):plot(gza(2,:));subplot(413):plot(gza(3,:));subplot(414):plot(gza(4,:));