demosig.m

来自「这是一个关于盲源分离独立成分分析方法（ICA）的软件包」· M 代码 · 共 35 行

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.

% 2 Apr 1998 Aapo Hyv鋜inen

%create source signals (independent components)
N=500; %data size

v=[0:N-1];
sig=[];
sig(1,:)=sin(v/2); %sinusoid
sig(2,:)=((rem(v,23)-11)/9).^5; %funny curve
sig(3,:)=((rem(v,27)-13)/9); %saw-tooth
sig(4,:)=((rand(1,N)<.5)*2-1).*log(rand(1,N)); %impulsive noise
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);
subplot(4,1,1);plot(v,mixedsig(1,:));
subplot(4,1,2);plot(v,mixedsig(2,:));
subplot(4,1,3);plot(v,mixedsig(3,:));
subplot(4,1,4);plot(v,mixedsig(4,:));