Displaying 1 result from an estimated 1 matches for "mysxy".
Did you mean:
mysql
2011 Sep 23
1
Cross Spectrum : Conversion of 2-D spectrum into a single complex array
...-pi
# since exp(i pi) = exp(-i pi) = -1 + 0i
# Why isn't it? Sampling issues? Nyquist has been satisfied.
plot( Sxy $ freq, Sxy $ phase[ ,1], type = 'l' )
# The real question (limited to a 2-D input):
# How to combine the amplitude/phase into a single
# complex valued cross spectrum?
mySxy = complex( real = Sxy $ spec[,1] * Sxy$ spec[,2],
imaginary = Sxy $ phase[ ,1] )
# This does give something affine to a plot of Sxy
# Compare
plot( Sxy, log="dB" )
# to:
plot( Sxy$freq, 10*log10( Re(mySxy) ), type='l')