clear all;
close all;
clc;

s1 = wavread('H:\Double-JJ\毕设\答辩\警笛.wav');
s2 = wavread('H:\Double-JJ\毕设\答辩\乐器.wav');
s3 = wavread('H:\Double-JJ\毕设\答辩\男语音.wav');
s4 = wavread('H:\Double-JJ\毕设\答辩\女语音.wav');

x1 = wavread('H:\Double-JJ\毕设\答辩\混合1.wav');
x2 = wavread('H:\Double-JJ\毕设\答辩\混合2.wav');
x3 = wavread('H:\Double-JJ\毕设\答辩\混合3.wav');
x4 = wavread('H:\Double-JJ\毕设\答辩\混合4.wav');

x(1,:)=x1;
x(2,:)=x2;
x(3,:)=x3;
x(4,:)=x4;

icasig = fastica(x);

figure;
subplot(4,1,1);
plot(s1);
xlabel('n');
ylabel('幅度');
title('警笛');
subplot(4,1,2);
plot(s2);
xlabel('n');
ylabel('幅度');
title('乐器');
subplot(4,1,3);
plot(s3);
xlabel('n');
ylabel('幅度');
title('男语音');
subplot(4,1,4);
plot(s4);
xlabel('n');
ylabel('幅度');
title('女语音');

figure;
subplot(4,1,1);
plot(x1);
xlabel('n');
ylabel('幅度');
title('混合信号1');
subplot(4,1,2);
plot(x2);
xlabel('n');
ylabel('幅度');
title('混合信号2');
subplot(4,1,3);
plot(x3);
xlabel('n');
ylabel('幅度');
title('混合信号3');
subplot(4,1,4);
plot(x4);
xlabel('n');
ylabel('幅度');
title('混合信号4');

figure;
subplot(4,1,1);
plot(icasig(1,:));
xlabel('n');
ylabel('幅度');
title('分离信号1');
subplot(4,1,2);
plot(icasig(2,:));
xlabel('n');
ylabel('幅度');
title('分离信号2');
subplot(4,1,3);
plot(icasig(3,:));
xlabel('n');
ylabel('幅度');
title('分离信号3');
subplot(4,1,4);
plot(icasig(4,:));
xlabel('n');
ylabel('幅度');
title('分离信号4');

wavwrite(icasig(1,:),'H:\Double-JJ\毕设\答辩\分离1.wav');
wavwrite(icasig(2,:),'H:\Double-JJ\毕设\答辩\分离2.wav');
wavwrite(icasig(3,:),'H:\Double-JJ\毕设\答辩\分离3.wav');
wavwrite(icasig(4,:),'H:\Double-JJ\毕设\答辩\分离4.wav');