clear all;
%VARIABLES
%Scale
sc = 120;
hu = 100;
alpha=sc/hu;
%Read audio
[signal,Fs]=audioread('test.mp3', [2 800001]);

x = signal(:,1);
w = hann(800000);
Y = zeros(1,800000);
xwin = x.*w;
X = fft(xwin,800000);
for i = 1:1:800000
    if round(i*alpha) <800000+1
        Y(i)= X(round(i*alpha));
    else
    Y(i) = 0;
    end
     YPhase(i) = Y(i)* exp(1j*1);
end
y = real(ifft(YPhase,800000));
for i = 1:1:800000
    if round(i*alpha^-1) <800000+1
         wfix(i) = (w(i)/w(round(i*alpha^-1)));
    else
         wfix(i) = 0;

    end
end

out = y.*wfix;
plot(out);
sound(out,Fs);