Short time FFT with octave smooth

Computes segmented FFT with optional fractional octave smoothing

-----------------------------
UPDATE: easy to use GUI added
-----------------------------
[Px,freq]=averfft(sig,Noct,Nfft)
sig = time domain data (e.g. audio file);
Noct = octave resolution for smoothing (e.g. Noct=3 results in a 1/3-oct. smooth),
default value is 0;
Nfft = FFT resolution, default value is 8192;
sampling frequency assumed to be 44.1 kHz.
If no output is defined averfft plots the frequency spectrum.