How can I reduce the computation time in the evaluation of cross-correlation between spectrograms with size 272x18209

The code used for compute the cross-correlation is the following. Unfortunately I cannot use xcorr2 built-in matlab func because it doesn't exactly the same operation.

Thanks!!