a=1;
m=1;
q=10;
s=wavread('speech');
v=wavread('machine');
y=s+a*v;
tempRsy=XCORR(s,y);
tempRy=XCORR(y);
mitt=ceil(length(tempRy)/2);
Ry=zeros(q,q);
Rsy=zeros(q,1);
j=0;
for g=1:q
    for i=1:q
        Ry(i,g)=tempRy(mitt+j);
        j=j+1;
    end
    j=j-q-1;
end
f=0;
for h=1:q
    Rsy(h,1)=tempRsy(mitt+m+f);
    f=f+1;
end
h=Ry\Rsy;
sut=filter(1,h,y);
[g,w]=freqs(1,h);
%-----------------------------------------Fel----------------------------
sum=0;
for i=1:length(s)
    sum=sum+(s(i)-sut(i))^2;
end
fel=sum/length(s);


%------------------------------------------Plottar-----------------------
subplot(3,3,1)
plot(s);
title('s');
xlabel('tid');
ylabel('förstärkning')
subplot(3,3,2)
plot(sut);
title('Sut');
xlabel('tid');
ylabel('förstärkning')
subplot(3,3,3)
plot(y);
title('y');
xlabel('tid');
ylabel('förstärkning')
subplot(3,3,4)
plot(fft(s));
title('s frekvens');
xlabel('frekvens');
ylabel('förstärkning')
subplot(3,3,5)
plot(fft(sut));
title('sut frekvens');
xlabel('frekvens');
ylabel('förstärkning')
subplot(3,3,6)
plot(fft(y));
title('y frekvens');
xlabel('frekvens');
ylabel('förstärkning')
subplot(3,3,7)
plot(w,g);
xlabel('frekvens');
ylabel('förstärkning')