Usage:
[Pxx,freq] = cpsd(x,y,Nfft,Fs,window,overlap,range,plot_type,detrend)

Estimate cross power spectrum of data "x" and "y" by the Welch (1967)
periodogram/FFT method.
See "help pwelch" for description of arguments, hints and references