ValueError: минимальный предел просмотра -36761.69947916667 меньше 1 и является недопустимым значением даты Matplotlib.

Question

Я использую эту часть кода;

df = pd.read_csv('Vertikale_Netzlast_2013.csv', header=6, sep=';', parse_dates=[[0, 1]], index_col=0, na_values=['n.v.'])
df.columns = ["time", "Load"]
df.Load = df.Load.interpolate()
plt.figure(figsize=(14, 5))
df.Load.plot()
plt.title('Vertical Grid Load Germany 2013')
plt.ylabel('Power [$MW$]')
plt.savefig('VerticalGridLoadGermany2013.png', bbox_inches='tight', dpi=150, transparent=True)
hann = np.hanning(len(df.Load.values))
Y = np.fft.fft(hann * df.Load.values)
N = ceil(len(Y) / 2 + 1)
print(N)
fa = 1.0 / (15.0 * 60.0) # every 15 minutes
print('fa=%.4fHz (Frequency)' % fa)
X = np.linspace(0, fa / 2, N, endpoint=True)
plt.plot(X, 2.0 * np.abs(Y[:N]) / N)
plt.xlabel('Frequency ($Hz$)')
plt.ylabel('vertical powergrid load ($MW$)')

Но я не запускаюсь, и python выдает эту ошибку

/ home / ozlemb / PycharmProjects / work / venv / lib / python3.6 / site-packages / numpy / core/numeric.py:531: ComplexWarning: приведение комплексных значений к реальному отбрасывает возвращаемый массив мнимой части (a, dtype, copy = False, order = order) dt = 0,00629 с (время выборки) fa = 159,00 Гц (частота) / home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/dateutil/parser/_parser.py:1204: UnknownTimezoneWarning: tzname B идентифицировано, но не понято.Передайте аргумент tzinfos, чтобы правильно возвратить дату / время с учетом часового пояса.В будущей версии это вызовет исключение.category = UnknownTimezoneWarning) 17523 fa = 0,0011 Гц (частота) Исключение при обратном вызове Tkinter (последний вызов был последним): файл "/usr/lib/python3.6/tkinter/init.py", строка 1702, в call__return self.func (* args) Файл "/usr/lib/python3.6/tkinter/__init.py", строка 746, в файле callit func (* args) "/ home / ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/backends/_backend_tk.py ", строка 346, в файле idle_draw self.draw ()" / home / ozlemb / PycharmProjects / work / venv / lib/python3.6/site-packages/matplotlib/backends/backend_tkagg.py ", строка 9, в файле draw super (FigureCanvasTkAgg, self) .draw ()" / home / ozlemb / PycharmProjects / work / venv / lib / python3.6 / site-packages / matplotlib / backends / backend_agg.py ", строка 402, в файле draw self.figure.draw (self.renderer)" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages / matplotlib / artist.py ", строка 50, в файле draw_wrapper возвращаемое рисование (художник, рендер, * args, ** kwargs) Файл" / home / ozlemb / PycharmProjects / work / venv /lib / python3.6 / site-packages / matplotlib / figure.py ", строка 1652, в рендере рисования, self, художники, self.suppressComposite) Файл" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/image.py ", строка 138, в файле _draw_list_compositing_images a.draw (renderer)" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/artist.py ", строка 50, в draw_wrapper возвращаемое рисование (художник, визуализатор, * args, ** kwargs) Файл" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/axes/_base.py ", строка 2604, в рисовании mimage._draw_list_compositing_images (средство визуализации, самообслуживание, исполнители)138, в файле _draw_list_compositing_images a.draw (renderer) "/home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/artist.py", строка 50, в draw_wrapper возвращают рисование (исполнитель,рендер, * args, ** kwargs) Файл "/ home / ozlemb / PycharmProjects / work / venv / lib / python3.6 / site-packages / matplotlib / axis.py ", строка 1185, в Draw Нарисуйте файл ticks_to_draw = self._update_ticks (renderer)/matplotlib/axis.py ", строка 1023, в _update_ticks tick_tups = list (self.iter_ticks ()) # iter_ticks вызывает файл локатора" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/axis.py ", строка 967, в файле iter_ticks majorLocs = self.major.locator () Файл" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/dates.py", строка 1229, в вызов файл self.refresh ()" /home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/dates.py ", строка 1249, в dmin обновления, dmax = self.viewlim_to_dt () Файл "/home/ozlemb/PycharmProjects/work/venv/lib/python3.6/site-packages/matplotlib/dates.py", строка 1000, в viewlim_to_dt.format (vmin)) ValueError: минимальный предел просмотра -36761.69947916667 меньше 1 и является недопустимым значением даты Matplotlib.Это часто случается, если вы передаете значение, отличное от datetime, в ось с датами и временем

Процесс завершен с кодом выхода 0

Это вызывает версию python (3) ??Я новичок и хочу обрабатывать данные сигналов с датчиков с помощью fft.

Спасибо

Answer 1

Дает ли эта строка ошибку?

X = np.linspace(0, fa/2, N, endpoint=True)

Если это так, то это должно быть вызвано тем, что N должно быть целым числом.В вашем случае N может быть float

Давайте проверим ссылку на umpy;

numpy.linspace(start, stop, num=50, endpoint=True, retstep=False, dtype=None)

num: int, необязательно

Количество сэмплов для генерации.

По умолчанию 50. Должно быть неотрицательным.

[https://docs.scipy.org/doc/numpy-1.15.0/reference/generated/numpy.linspace.html]

Answer 2

Если вы используете Python3, N = len(Y)/2+1 может вычислять не целое значение.После этого вы будете пытаться нарезать список с помощью Y[:N], но N может быть равен, например, 2.5.

См. Следующий код:

nums = [1, 2, 3, 4, 5, 6, 7]
N = len(nums) / 2 + 1
print("N =", N);
print(nums[:N])

СпособВы вычисляете N, чтобы сделать его равным 4.5, потому что len(nums) / 2 равно 3.5.Вместо этого я бы сделал from math import ceil, а затем

N = ceil(len(nums) / 2);

Это сделает N равным 4.

Это, конечно, зависит от того, что именно вы пытаетесь нарезать.Суть в том, что вы должны убедиться, что N является целым числом, а не значением с плавающей запятой.