我的代码
import matplotlib.pyplot as plt
plt.style.use("ggplot")
import numpy as np
from mtspec import mtspec
from mtspec.util import _load_mtdata
data = np.loadtxt('262_V01_C00_R000_TEx_BL_4096H.dat')
spec,freq,jackknife,f_statistics,degrees_of_f = mtspec(data=data, delta= 4930.0, time_bandwidth=4 ,number_of_tapers=5, nfft= 4194304, statistics=True)
fig = plt.figure()
ax2 = fig
ax2.plot(freq, spec, color='black')
ax2.fill_between(freq, jackknife[:, 0], jackknife[:, 1],color="red", alpha=0.3)
ax2.set_xlim(freq[0], freq[-1])
ax2.set_ylim(0.1E1, 1E5)
ax2.set_xlabel("Frequency $")
ax2.set_ylabel("Power Spectral Density $)")
plt.tight_layout()
plt.show()
问题出在我的代码的绘图部分。我应该更改什么?我在Ubuntu上使用Python 2.7。
答案 0 :(得分:9)
您将ax2
分配给未定义figure
方法的plot
对象。您想使用plt.axes
代替
ax2 = plt.axes()
# Instead of ax2 = fig
答案 1 :(得分:0)
而不是通过以下方式调用图形:
ax2 = fig #which just references the figure
尝试使用gca()
提取轴:
ax2 = fig.gca() #which is used to extract the axes