如何将第二个子批次的y轴范围设置为例如[01000]? 我的数据(文本文件中的一列)的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()
其他的改进也值得赞赏!
使用axes objects是一个很好的方法。如果要与多个图形和子图交互,它会有帮助。要直接添加和操作轴对象,请执行以下操作:
有时,在绘制数据之前,您确实需要设置轴限制。在这种情况下,可以设置
Axes
或AxesSubplot
对象的“自动缩放”功能。感兴趣的功能是set_autoscale_on
、set_autoscalex_on
和set_autoscaley_on
。在您的情况下,您希望冻结y轴的限制,但允许x轴展开以容纳您的数据。因此,您需要将
autoscaley_on
属性更改为False
。下面是代码中FFT子块片段的修改版本:如http://www.mofeel.net/582-comp-soft-sys-matlab/54166.aspx所示
注意:命令必须在绘图后执行!
相关问题 更多 >
编程相关推荐