fetal_plots_bsica.m

来自「image separation using neural net」· M 代码 · 共 52 行

% Demonstrate ICA by extracting the fetal heart beat from the ECG of a
% pregnant mother.

%
% David Gleich
% CS 152 - Neural Networks
% 12 December 2003
%

% load the data
X = icadata('ECG');

k = 8;

fprintf('Testing BSICA algorithm...\n');


plotrows(Xecg);
xlabel('Sample Number');
subplot(size(Xecg, 1), 1, 4);
ylabel('Source Signals');



fprintf('Press a key for ICA iterations.\n');
pause;

B = bsica(Xecg, k, 'rate', 0.001, 'niter', 1);

for (ii=1:3)
    fprintf('Iteration %i\n', ii);
    
    Y = B*Xecg;
    

    plotrows(Y);
    xlabel('Sample Number');
    subplot(size(Xecg, 1), 1, 4);    
    ylabel(sprintf('ICA signals, iteration %i', ii));
   
    
     B = bsica(Xecg, k, 'rate', 0.001, 'niter', 1, 'B', B);
    
    fprintf('Press a key for the next iteration.\n');
    pause;
end;