[s1,fs]=wavread('s1.wav');
[s2,fs]=wavread('s2.wav');
[s3,fs]=wavread('colt16.wav');
sz=min(length(s1),size(s3,1));
Source(1,:)=s1(1:sz)';
Source(2,:)=s2(1:sz)';
Source(3,:)=s3(1:sz,1)';
ns=3;
%%%%%%%%%%%%%%%%%%%%%%linear
nx=8;
A=randn(nx,ns)
xlin=A*Source;
%%%%%%%%%%%%%%%%%convolutive
nh=3
h11=randn(1,nh)*1;
h12=randn(1,nh)*1;
h13=randn(1,nh)*1;
h21=randn(1,nh)*1;
h22=randn(1,nh)*1;
h23=randn(1,nh)*1;
h31=randn(1,nh)*1;
h32=randn(1,nh)*1;
h33=randn(1,nh)*1;
h41=randn(1,nh)*1;
h42=randn(1,nh)*1;
h43=randn(1,nh)*1;
h51=randn(1,nh)*1;
h52=randn(1,nh)*1;
h53=randn(1,nh)*1;
h61=randn(1,nh)*1;
h62=randn(1,nh)*1;
h63=randn(1,nh)*1;
h71=randn(1,nh)*1;
h72=randn(1,nh)*1;
h73=randn(1,nh)*1;
h81=randn(1,nh)*1;
h82=randn(1,nh)*1;
h83=randn(1,nh)*1;
h91=randn(1,nh)*1;
h92=randn(1,nh)*1;
h93=randn(1,nh)*1;
S(1,:)=conv(Source(1,:),h11);
S(2,:)=conv(Source(1,:),h21);
S(3,:)=conv(Source(1,:),h31);
S(4,:)=conv(Source(1,:),h41);
S(5,:)=conv(Source(1,:),h51);
S(6,:)=conv(Source(1,:),h61);
S(7,:)=conv(Source(1,:),h71);
S(8,:)=conv(Source(1,:),h81);
S(9,:)=conv(Source(1,:),h91);
S(1,:)=S(1,:)+conv(Source(2,:),h12);
S(2,:)=S(2,:)+conv(Source(2,:),h22);
S(3,:)=S(3,:)+conv(Source(2,:),h32);
S(4,:)=S(4,:)+conv(Source(2,:),h42);
S(5,:)=S(5,:)+conv(Source(2,:),h52);
S(6,:)=S(6,:)+conv(Source(2,:),h62);
S(7,:)=S(7,:)+conv(Source(2,:),h72);
S(8,:)=S(8,:)+conv(Source(2,:),h82);
S(9,:)=S(9,:)+conv(Source(2,:),h92);
S(1,:)=S(1,:)+conv(Source(3,:),h13);
S(2,:)=S(2,:)+conv(Source(3,:),h23);
S(3,:)=S(3,:)+conv(Source(3,:),h33);
S(4,:)=S(4,:)+conv(Source(3,:),h43);
S(5,:)=S(5,:)+conv(Source(3,:),h53);
S(6,:)=S(6,:)+conv(Source(3,:),h63);
S(7,:)=S(7,:)+conv(Source(3,:),h73);
S(8,:)=S(8,:)+conv(Source(3,:),h83);
S(9,:)=S(9,:)+conv(Source(3,:),h93);
xconv=S(1:nx,:);