sigprocess.m

一种快速的独立主元算法

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 size

%v=[0:N-1];% sig=[];
load('f:\gzga09');
load('f:\gzgs11');
load('f:\gzpl16');
load('f:\gzta07');
sig(1,:)=gzga09(20001:30000,2); %sinusoid
sig(2,:)=gzgs11(20001:30000,2); %funny curve
sig(3,:)=gzpl16(1:10000,2); %saw-tooth
sig(4,:)=gzta07(10001:20000,2); %impulsive noise
figure(1);
subplot(411):plot(sig(1,:));
subplot(412):plot(sig(2,:));
subplot(413):plot(sig(3,:));
subplot(414):plot(sig(4,:));
for t=1:4
sig(t,:)=sig(t,:)/std(sig(t,:));
end

%remove mean (not really necessary)

[sig mean]=remmean(sig);

%create mixtures

Aorig=rand(size(sig,1));
mixedsig=(Aorig*sig);
figure(2);
subplot(411):plot(mixedsig(1,:));
subplot(412):plot(mixedsig(2,:));
subplot(413):plot(mixedsig(3,:));
subplot(414):plot(mixedsig(4,:));
gza=FASTICA(mixedsig);
figure(3);
subplot(411):plot(gza(1,:));
subplot(412):plot(gza(2,:));
subplot(413):plot(gza(3,:));
subplot(414):plot(gza(4,:));