Novel approach

We go now to frequency domain, where convolution becomes simple multiplication at every frequency k, by taking an N-point FFT of all those impulse responses:

We now try to invert this linear equation system at every frequency k, and for every virtual microphone v:

This over-determined system doesn't admit an exact solution, but it is possible to find an approximated solution with the  Least Squares method