可以将文章内容翻译成中文,广告屏蔽插件可能会导致该功能失效(如失效，请关闭广告屏蔽插件后再试):
问题:
How can I set the y axis range of the second subplot to e.g. [0,1000] ? The FFT plot of my data (a column in a text file) results in a (inf.?) spike so that the actual data is not visible.
pylab.ylim([0,1000])
has no effect, unfortunately. This is the whole script:
# 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()
Other improvements are also appreciated!
回答1:
As found in http://www.mofeel.net/582-comp-soft-sys-matlab/54166.aspx
pylab.ylim([0,1000])
Note: The command has to be executed after the plot!
回答2:
Using axes objects is a great approach for this. It helps if you want to interact with multiple figures and sub-plots. To add and manipulate the axes objects directly:
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()
回答3:
Sometimes you really want to set the axes limits before you plot the data. In that case, you can set the "autoscaling" feature of the Axes or AxesSubplot object. The functions of interest are set_autoscale_on, set_autoscalex_on, and set_autoscaley_on.
In your case, you want to freeze the y axis' limits, but allow the x axis to expand to accommodate your data. Therefore, you want to change the autoscaley_on property to False. Here is a modified version of the FFT subplot snippet from your code:
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))
