% A3)

% Einlesen der Bilder und Umwandeln in double (benötigt für fastica)
Bild1 = double(imread('bss_1.png'));
Bild2 = double(imread('bss_2.png'));
Bild3 = double(imread('bss_3.png'));
Bild4 = double(imread('bss_4.png'));

% Jede Matrix in einee Spaltenvektor umwandeln
Bild1v = Bild1(:);
Bild2v = Bild2(:);
Bild3v = Bild3(:);
Bild4v = Bild4(:);

% Daten mittelwertzentrieren
mBild1v = Bild1v - mean(Bild1v);
mBild2v = Bild2v - mean(Bild2v);
mBild3v = Bild3v - mean(Bild3v);
mBild4v = Bild4v - mean(Bild4v);

% Bilder zu einer 4x(64*64)-Matrix, jede Zeile entspricht einem Signal
mixedsig = [mBild1v'; mBild2v'; mBild3v'; mBild4v'];

[icasig, A, W] = fastica( mixedsig );

% Wieder zurück in 64x64-Matrix
Bild1neu = reshape( icasig(1,:), [64, 64] );
Bild2neu = reshape( icasig(2,:), [64, 64] );
Bild3neu = reshape( icasig(3,:), [64, 64] );
Bild4neu = reshape( icasig(4,:), [64, 64] );

subplot(2, 2, 1), imshow( Bild1neu )
subplot(2, 2, 2), imshow( Bild2neu )
subplot(2, 2, 3), imshow( Bild3neu )
subplot(2, 2, 4), imshow( Bild4neu )