I'm testing some code to perform deconvolution of two audio signals to recover the impulse response. Presently, as part of testing, I am simply deconvovolving two identical signals.
I have the assumption that I should be able to recover a simple dirac-delta impulse response from this, but for some reason, there is a lot of spurious noise in the impulse response result.
I am using standard 8-byte doubles for data storage and FFTW as my FFT library. I'm also normalising the impulse response.
Is my assumption correct (regarding the unit impulse response)?
I'm still trying to figure out how to correctly compensate for wrap-around (circular convolution) by correctly padding the input buffers. I am using a similar strategy to my working convolution code [ M+L-1 ]