emmanouil@5: %function [y,yh,ys,fr0] = wave2fft2wave(x,fs,w,N,t,nH,minf0,maxf0,f0et,maxhd,stocf)
emmanouil@5: function [y] = wave2fft2wave(x,w,N)
emmanouil@5: %e.g. y = wave2fft2wave(x,hamming(2025),4096);
emmanouil@5: 
emmanouil@5: 
emmanouil@5: M = length(w);   % analysis window size
emmanouil@5: Ns = 1024;                               % FFT size for synthesis
emmanouil@5: H = 256;                                 % hop size for analysis and synthesis
emmanouil@5: soundlength = length(x);                 % length of input sound array
emmanouil@5: hNs = Ns/2;                              % half synthesis window size
emmanouil@5: hM = (M-1)/2;                            % half analysis window size
emmanouil@5: pin = max(hNs+1,1+hM);   % initialize sound pointer to middle of analysis window
emmanouil@5: pend = soundlength-max(hM,hNs);          % last sample to start a frame
emmanouil@5: fftbuffer = zeros(N,1);                  % initialize buffer for FFT
emmanouil@5: y = zeros(soundlength+Ns/2,1);           % output sine component
emmanouil@5: w = w/sum(w);                            % normalize analysis window
emmanouil@5: sw = zeros(Ns,1);
emmanouil@5: ow = triang(2*H-1);                      % overlapping window
emmanouil@5: ovidx = Ns/2+1-H+1:Ns/2+H;               % overlap indexes
emmanouil@5: sw(ovidx) = ow(1:2*H-1);
emmanouil@5: bh = blackmanharris(Ns);                 % synthesis window
emmanouil@5: bh = bh ./ sum(bh);                      % normalize synthesis window
emmanouil@5: sw(ovidx) = sw(ovidx) ./ bh(ovidx);
emmanouil@5: 
emmanouil@5: 
emmanouil@5: while pin<pend
emmanouil@5:     
emmanouil@5:     xw = x(pin-hM:pin+hM).*w(1:M);         % window the input sound
emmanouil@5:     fftbuffer(1:(M+1)/2) = xw((M+1)/2:M);  % zero-phase window in fftbuffer
emmanouil@5:     fftbuffer(N-(M-1)/2+1:N) = xw(1:(M-1)/2);
emmanouil@5:     X = fft(fftbuffer);                    % compute the FFT
emmanouil@5:     ri= pin-hNs;                           % input sound pointer for residual analysis
emmanouil@5:     yw = ifft(X);    
emmanouil@5:     y(ri:ri+Ns-1) = y(ri:ri+Ns-1)+yw(1:Ns).*sw;
emmanouil@5:     pin = pin+H;     
emmanouil@5:     
emmanouil@5: end
emmanouil@5: 
emmanouil@6: y = (max(x)/max(y))*y; % scale y to original amplitude
emmanouil@6: 
emmanouil@5: %wavwrite(y,44100,'test.wav');