alpha = 1/2 //Scaling
x <- 1024 samples
w <- hann(1024)
xwin <- x.w
X <- DFT(xwin)
for i = 0 to 1024
	Y[i] <- X[round(i*alpha)]
	YPhase[i] <- Y* exp(i*phi[i])
endfor
y = IDFT(YPhase)
out <- y[n].(w[n]/w[α^-1*n])