Matlab implementation

[x, fs, nbits] = wavread(`...');
N = length(x);

Nwin = 256;        % window size
Noverlap = Nwin/2; % 50 percent overlap
zpf = 1;  Nfft = Nwin*zpf;            
X = zeros(Nfft, round(N/Noverlap-1));
win = hanning(Nwin); %COLA window

for i=0:N/Noverlap-2
  ix = Noverlap*i+[1:Nwin];
  X(:,i+1) = fft(x(ix).*win, Nfft);
end

%Reconstruct
y = zeros(N,1);
for i=0:N/Noverlap-2
  ix = Noverlap*i+[1:Nwin];
  y(ix) = y(ix) + real(ifft(X(:,i+1)));
end