15.4.8.3 Проблемы точности ¶

Наборы инструментов OpenGL для графики ("qt" и "fltk") используют одинарную точность для отрисовки. Это ограничение особенно применимо к графикам временных рядов по отношению к датам в последовательности, как используется функциями datenum, datestr, datestruct, и datetick.

Последовательные даты кодируют отметки времени как количество дней, прошедших с нулевого года, причём часы, минуты и секунды представлены дробной частью. 31 декабря 1999 года последовательная дата составляла 730485. Переменная двойной точности с этой целой частью позволяет разрешение дробной части 1.2e-10, что соответствует примерно 5 микросекундам. Но с одинарной точностью разрешение снижается до примерно 0.06, что соответствует 45 минутам. Любая попытка построить график временных данных с более высокой точностью приведёт к искажению графика.

В качестве обходного решения можно использовать набор инструментов для графики "gnuplot" или вычесть 2000 лет — т.е., datenum (2000, 0, 0) или 730485 — из значений времени. Из-за того, что структура календаря повторяется каждые 2000 лет, связь между годом, месяцем, днём месяца и днём недели останется неизменной, и метки и подписи меток, созданные функцией datetick, останутся корректными. Только у годов не будет указан тысячелетний разряд. Таким образом, «2020» будет напечатано как «20». Например:

# timestamps of 24 hours in one minute steps
t = datenum (2020, 1, 1):(1/1440):datenum (2020, 1, 2);

# some example time series data
x = -cos (2*pi*t) + rand (size (t)) / 10;

subplot (1, 2, 1);
plot (t, x);
datetick ("x");
xlabel ("serial date");
title ("problem");

subplot (1, 2, 2);
plot (t - 730485, x);
datetick ("x");
xlabel ("2000 years off");
title ("workaround");

Результат можно увидеть на рисунке 15.8.

Аналогично, другие данные можно перевести или перемасштабировать, чтобы обойти эту проблему.