2番目のサブプロットの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()
他の改良もまた歓迎されています!
http://www.mofeel.net/582-comp-soft-sys-matlab/54166.aspxに見られるように
pylab.ylim([0,1000])
注:コマンドはプロットの後に実行する必要があります。
axesオブジェクト を使うことはこれのための素晴らしいアプローチです。複数の図やサブプロットと対話したい場合に役立ちます。 Axesオブジェクトを直接追加して操作するには:
import matplotlib.pyplot as plt
fig = plt.figure(figsize=(12,9))
signal_axes = fig.add_subplot(211)
signal_axes.plot(xs,rawsignal)
fft_axes = fig.add_subplot(212)
fft_axes.set_title("FFT")
fft_axes.set_autoscaley_on(False)
fft_axes.set_ylim([0,1000])
fft = scipy.fft(rawsignal)
fft_axes.plot(abs(fft))
plt.show()
時々、あなたは本当に軸の範囲を設定したい/ あなたはデータをプロットする 。その場合、Axes
またはAxesSubplot
オブジェクトの "自動スケーリング"機能を設定できます。対象となる関数はset_autoscale_on
、set_autoscalex_on
、およびset_autoscaley_on
です。
あなたのケースでは、あなたはy軸の限界を凍結したいが、あなたのデータを収容するためにx軸が広がるのを許しなさい。したがって、autoscaley_on
プロパティをFalse
に変更します。これはあなたのコードからのFFTサブプロットスニペットの修正版です:
fft_axes = pylab.subplot(h,w,2)
pylab.title("FFT")
fft = scipy.fft(rawsignal)
pylab.ylim([0,1000])
fft_axes.set_autoscaley_on(False)
pylab.plot(abs(fft))