コード例 #1
0
ファイル: SilvaJ_HW6_p2.py プロジェクト: jasminks/ASTR260
#importing necessary components to plot from modules

from pylab import plot,show,xlim,ylim
from numpy import loadtxt
from mydft import dft

#importing sunspot file as an array to plot
z=loadtxt("sunspots.txt",float)

#z is now an array

#t is time in months
t=z[:,0]

#n is number of sunspots occuring within that period
n=z[:,1]

#taking the dft of the sunspots
c = dft(n)

#plotting fourier coefficients
plot (abs(c)**2)
xlim(0,3142/2)
ylim(0,3*10**9)
show()
コード例 #2
0
ファイル: myfft.py プロジェクト: jasminks/ASTR260
from numpy import loadtxt
from pylab import plot,xlim,show
from mydft import dft 

y = loadtxt("pitch.txt",float)
c = dft(y)
plot(abs(c))
xlim(0,500)
show()