Welch's overlapped method for amplitude by frequency estimation pwelch  is a useful matlab function for computing real/magnitudes for power using WOSA method, but I'd like to apply the same thinking in order to output complex pairs for amplitude (not power) spectra.

Does anyone know if this is possible? Or if there is an existing script that I might be able to modify?

Thanks!
I've found this:
Implementing Welch&#39;s method for Power Spectral Density

That seems to answer the question. However, the author encountered a problem, which I believe may be addressed by modifying the code that was identified as problematic in the thread by:

dataA(istart:iend,j) = data1(istart:iend,j).*window;
dataB(istart:iend,j) = data2(istart:iend,j).*window;