clc,clear
[x,fs,bits]=wavread('12.wav');
y=x(1:4096,1);
Y=fft(y);
magY=abs(Y);
b=[];
for i=0:126;
n=4096;
x1=x(1+n*i:n+n*i);
X1=fft(x1);
magX=abs(X1);
S=(magX1'.^2-magY.^2);
S1=abs(S).^0.5;
s1=ifft(S1);
a=s1';
b=[b a];
end
x2=b';
plot(x2);
sound(x2,fs,bits);
wavwrite(x2,fs,'13.wav');