두 번째 서브 플롯의 y 축 범위를 어떻게 [0,1000]으로 설정합니까? 내 데이터 (텍스트 파일의 열)의 FFT 플롯은 실제 데이터가 표시되지 않도록 (inf ??) 스파이크를 발생시킵니다.
pylab.ylim([0,1000])
불행히도 효과가 없습니다. 이것은 전체 스크립트입니다.
# based on http://www.swharden.com/blog/2009-01-21-signal-filtering-with-python/
import numpy, scipy, pylab, random
xs = []
rawsignal = []
with open("test.dat", 'r') as f:
for line in f:
if line[0] != '#' and len(line) > 0:
xs.append( int( line.split()[0] ) )
rawsignal.append( int( line.split()[1] ) )
h, w = 3, 1
pylab.figure(figsize=(12,9))
pylab.subplots_adjust(hspace=.7)
pylab.subplot(h,w,1)
pylab.title("Signal")
pylab.plot(xs,rawsignal)
pylab.subplot(h,w,2)
pylab.title("FFT")
fft = scipy.fft(rawsignal)
#~ pylab.axis([None,None,0,1000])
pylab.ylim([0,1000])
pylab.plot(abs(fft))
pylab.savefig("SIG.png",dpi=200)
pylab.show()
다른 개선 사항도 감사합니다!