; IDL Version 6.3 (linux x86_64 m64)
; Journal File for jmcbride@pluto
; Working directory: /home/jmcbride
; Date: Tue Feb  2 20:33:37 2010
 
t = findgen(100) * 2 * !pi / 100
x = sin(t)
plot, x, charsize = 2, linestyle = 2, thick = 2
ft = fft(x)
pow = ft * conj(ft)
plot, pow, charsize = 2, linestyle = 2, thick = 2
pow = shift(pow, n_elements(pow) / 2.)
plot, pow, charsize = 2, linestyle = 2, thick = 2
plot, pow, charsize = 2, thick = 2
frq = findgen(100) - 50
plot, frq, pow, charsize = 2, thick = 2
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5]
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5], xtitle = textoidl('$\nu$ (MHz)')
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5], xtitle = textoidl('\nu (MHz)')
x2 = sin(2 * t)
ft2 = fft(x2)
pow = ft2 * conj(ft2)
oplot, pow, linestyle = 2, color = !red
pow = shift(pow, n_elements(pow) / 2.)
pow2 = ft2 * conj(ft2)
pow2 = shift(pow2, n_elements(pow2)/2.)
erase
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5], xtitle = textoidl('\nu (MHz)')
oplot, pow2, linestyle = 2, color = !red
oplot, frq, pow2, linestyle = 2, color = !red
legend, ['Power 1', 'Power 2'], linestyle = [0, 2], color = [!white, !red]
legend, ['Power 1', 'Power 2'], linestyle = [0, 2], color = [!white, !red], /right
sampFreq = 1e7
sigFreq = 1e6
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5], xtitle = textoidl('\nu (MHz)'), title = 'Signal Frequency = ' + strcompress(sigFreq, /remove)
plot, frq, pow, charsize = 2, thick = 2, xrange = [-5, 5], xtitle = textoidl('\nu (MHz)'), title = 'Signal Frequency = ' + strcompress(sigFreq / 1e6, /remove) + ' MHz'
print, sigFreq / 1e6, format = '(F3.1)'
;1.0
print, sigFreq / 1e6, format = '(F3.2)'
;***
print, sigFreq / 1e6, format = '(F4.2)'
;1.00