Currently, I find my impulse response by performing an FFT of size 32k
on my entire measurement signal, and my reference signal. I divide
each by the scaling factor 1/sqrt(M)
Scaling here is unnecessary since you are going to divide the two spectra anyway. Any scale factor will cancel out.
I then divide the two signals, and perform an inverse FFT, ...