t60.m

Returns an estimate of t60, a measure of reverberation time.

Pass a recorded noise burst or impulse, and this script returns an estimate of t60. Pass only the part of the recording after the offset of the noise. It uses Schroeder's integrated impulse decay method, and then performs a cumulative r-squared to find the straight part of the slope. eg, if you recorded a noise segment that is 8192 samples long, syntax would be:

rt = t60(noise(8193:length(noise)),srate);

The script will also plot the decay function if desired.