Python, Matplotlib, subplot: 軸の範囲を設定するには?

2つ目のサブプロットのY軸の範囲を[0,1000]などに設定するにはどうしたらいいですか? 私のデータ(テキストファイルの列)のFFTプロットは、実際のデータが見えないほどの(無限の)スパイクになってしまいます。

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])

注意:コマンドはプロットの後に実行しなければなりません。

解説 (4)

axes objects]1を使うのは、このための素晴らしいアプローチです。複数の図形やサブプロットを操作したい場合に役立ちます。軸オブジェクトを直接追加して操作するには

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()
解説 (1)

データをプロットする前に、どうしても軸の限界を設定したい場合があります。 そのような場合には、AxesまたはAxesSubplotオブジェクトの"autoscaling"機能を設定することができます。 対象となる関数は 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))
解説 (0)