WAV files (omnidirectional mics?) from 2/12/07 Dinkelspiel measurements:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
MATLAB scripts and info from Jonathan Abel:
Note: the wavread command puts the stereo response in two rows, rather than the two columns expected by the script(s).
acBE676486.wav - single sine sweep.
acBE676486.m - script
for forming the measured impulse response given the system response to the single sine sweep acBE676486.wav.
ssBE676786.wav - repeated sine sweep.
ssBE676786.m - script
for forming the measured impulse response given
the system response to the repeated sine sweep ssBE676786.wav. Note
that this script uses circular convolution, and processes the system
response during the second sweep. Assuming that the system response
is shorter than the sweep, the response during the second sweep is
equal to the _circular_ convolution of the sweep and the system. The
system is then recovered by dividing the transform of "one period" of
the response by that of the sweep.
ModelLate.m - script
for analyzing reverberation impulse response
late-field initial equalization and decay rate. Set the variable
"ir" to a single impulse response (or balloon pop), arranged as a
column. The script should plot the impulse response and spectrogram,
and then find and plot the decay rates. When finding the decay
rates, the band energy as a function of time is plotted along with a
fitted exponential decay (linear on the dB scale) for each
half-octave band. You might have to adjust the starting and ending
times of the fit so that the modeled section doesn't include the
noise floor or an initial decay, which is sometimes present, and
often faster than the eventual T60. (You might want to run this
script on one of the pair of impulse responses in the variable irSU
in the file SU.mat to compare Dinkelspiel to the CCRMA lobby as Patty
and I measured it before the rennovation.)
pinkbands.m - script
which generated the octave-wide pink noise bands. The responses should be analyzed for their power and their decay rate.