Journal influence
Bookmark
Next issue
A software package for nonuniform time series analysis based on continuous wavelet transformation
Abstract:Wavelet transformation is one of the methods of data time-frequency analysis. In practice, the researcher often has to analyse nonequidistant (uneven) time series. For this reason, the paper considers algorithms for their continuous wavelet transformation. The authors propose an algorithm for obtaining an array of transformation even shifts, taking into account uneven source date. It means that they choose the interval of involuntary discretization, determine the number of shifts and then calculate their value. The proposed algorithm is the base for an algorithm of continuous wavelet transformation of nonequidistant time series. The process of estimating transformation coefficients uses only samples of time series that are in the width of a wavelet. The advantage of thу proposed algorithm is that the result of the transformation is an even representation. The developed algorithms are used in a complex of programs. The paper shows the results of experiments for models of nonuniform time series with p-transformation and with “jitter” (with uniform distribution and triangular distribution). The velocity of the algorithm can be increased by taking into account the effective radius of the mother wavelet and calculating its width.
Аннотация:На практике исследователю часто приходится иметь дело с неэквидистантными (нерегулярными) временными рядами. В данной статье рассматриваются алгоритмы их непрерывного вейвлет-преобразования, которое является одним из методов частотно-временного анализа данных. Предлагается алгоритм для получения массива регулярных сдвигов преобразования, учитывающий нерегулярность исходных данных: выбирается интервал принудительной дискретизации, определяется число сдвигов и вычисляются их значения. На основе предложенного алгоритма разработан алгоритм непрерывного вейвлет-преобразования неэквидистантных временных рядов. В процессе вычисления оценки коэффициентов преобразования учитываются лишь те отсчеты временного ряда, которые попадают в ширину вейвлета. Преимуществом данного алгоритма является то, что результат преобразования носит регулярный характер. Разработанные алгоритмы реализованы в комплексе программ. Показаны результаты экспериментов для моделей неэквидистантных временных рядов с p-преобразованием и с дрожанием (с распределением случайной величины по равномерному закону и закону Симпсона). Быстродействие алгоритма повышается за счет учета эффективного радиуса базового вейвлета и вычисления его ширины.
Authors: Prokhorov S.A. (sp.prokhorov@gmail.com) - Samara State Aerospace University, Samara, Russia, Ph.D, A.A. Stolbova (anastasiya.stolbova@bk.ru) - Samara National Research University (Assistant), Samara, Russia | |
Keywords: effective radius, nonuniform time series, continuous wavelet transformation |
|
Page views: 7805 |
PDF version article Full issue in PDF (29.80Mb) |
Одним из активно развивающихся методов спектрального анализа данных является вейвлет-преобразование, коэффициенты которого вычисляются следующим образом:
где x(t) – случайный процесс; ψ(t) – выбранный анализирующий вейвлет; a ≠ 0 – параметр масштаба; b ≥ 0 – параметр сдвига [1–3]. Спектральные методы применяются для анализа медицинских сигналов, таких как электрокардиограмма, электроэнцефалограмма, вариабельность сердечных ритмов (ВСР), анализа космофизических явлений и в других областях [4–7]. Такие процессы, как правило, являются нестационарными по частоте, и спектральный анализ не поз- воляет локализовать частоты во времени. В этом случае применимы частотно-временные методы анализа, к которым относится вейвлет-преобразование. Часто исследователю приходится иметь дело с нестационарными и неэквидистантными данными (например ВСР), то есть Δtk = tk+1 – tk = random. Типовые модели процессов с нерегулярной дискретизацией данных и методы их анализа рассмотрены в [8, 9]. В работе [10] исследуется метод адаптивного вейвлет-преобразования для модели с пропусками наблюдений, недостатком которого является то, что необходимо знать, в какой момент времени отсчет был пропущен. При вычислении оценки вейвлет-коэффициен- тов нерегулярных процессов методом прямоуголь- ников выражение (1) преобразуется к следующему виду: где N – число отсчетов реализации неэквидистантного временного ряда. Выражение для оценки коэффициентов методом трапеций имеет вид
Алгоритм получения массива сдвигов для вейвлет-преобразования неэквидистантных временных рядов. При вычислении вейвлет-преобразования регулярных процессов значения сдвигов определяются как значения, кратные интервалу дискретизации, и для выполнения преобразования достаточно знать лишь порядковый номер сдвига. В случае, когда интервал дискретизации не является константой, предлагается следующий алгоритм получения массива сдвигов. 1. Выбрать принудительный интервал дискретизации ∆t0, определяемый как минимум, максимум или среднее значение интервалов дискретизации неэквидистантного временного ряда: 2. Определить число отсчетов предполагаемого регулярного временного ряда: 3. Определить число сдвигов: 4. Получить массив сдвигов: bj = j×K×Dt0, где j = 0, …, Nb – 1. Алгоритм оценки коэффициентов вейвлет-преобразования неэквидистантных временных рядов. С учетом описанных выше особенностей алгоритм оценки коэффициентов вейвлет-преобразования временных рядов с нерегулярной дискретизацией будет следующим. 1. Загрузить исследуемый процесс, представленный неэквидистантным временным рядом
2. Получить массив масштабов:
3. Получить массив сдвигов по алгоритму, предложенному выше:
4. Для текущего значения масштаба ai рассчитать ширину вейвлета wt: wt = 8ai Dt, где ∆t – эффективный радиус базового вейвлета. 5. Вычислить значение коэффициента вейвлет-преобразования, используя выражение (2) или (3). При этом следует учитывать лишь те отсчеты нерегулярного временного ряда, соответствующие временные метки которых попадают в ширину вейвлета относительно текущего сдвига bj: В связи с тем, что на границах сигнала невозможно использовать все значения вейвлета и появляются погрешности, вычисления рекомендуется производить в случае, если выполняется следую- щее условие: 6. Повторить пункты 4 и 5 для всех масштабов ai и сдвигов bj. Таким образом, результатом работы алгоритма является матрица коэффициентов вейвлет-преоб- разования. Отметим, что результат преобразования носит регулярный характер, так как функция масштаба обратно пропорциональна равномерно дискретизированной частоте, а предложенный алгоритм вычисления сдвигов предполагает их регулярность. Программный комплекс. На основе рассмотренных алгоритмов был разработан программный комплекс на языке C# на платформе .NET с применением технологии паттернового проектирования. Он состоит из пяти модулей. 1. Модуль получения исходного процесса, предназначенный для генерирования исходных последовательностей следующих видов: - случайные стационарные с заданным видом корреляционной функции; - случайные нестационарные с заданным видом корреляционной функции; - детерминированные. Данный модуль позволяет получать процессы как с регулярной дискретизацией данных, так и с нерегулярной. 2. Модуль получения спектральных характеристик сигнала при помощи следующих операций: - преобразование Фурье; - оконное преобразование Фурье; - вейвлет-преобразование с регулярной дискретизацией; - вейвлет-преобразование с нерегулярной дискретизацией. 3. Модуль, предназначенный для определения погрешностей вычисленных коэффициентов вейвлет-преобразования. 4. Модуль построения вейвлетов. Данный модуль предназначен для вычисления вейвлет-функций. В системе используются десять основных видов вейвлетов. Для их построения необходимо указать число отсчетов вейвлета N, интервал дискретизации вейвлета Dt, масштаб a и сдвиг b. 5. Модуль имитационного моделирования предназначен для моделирования вейвлет-преобразования процессов с нерегулярной дискретизацией данных, оценки адекватности разработанных алгоритмов. Экспериментальные исследования. Вычислим вейвлет-преобразование для вейвлета Гаусса 8-го порядка и синусоиды: где y(t) = (t8 – 28t6 + 210t4 – 420t2 + 105)exp(–t2/2). В качестве эксперимента генерировалось по 29 реализаций временных рядов для каждого типа нерегулярной дискретизации, так как для доверительной вероятности Pд = 0,95 число испытаний M = 29 при любом законе распределения погрешностей. Для оценки ошибки вычисления каждого вейвлет-коэффициента определим СКО получен- ных результатов: Для сравнения скейлограмм вейвлет-преобразования введем следующую величину:
где В таблице показаны результаты сравнения скейлограмм для временных рядов с пропусками наблюдений (p-преобразование) для различного числа пропусков q и временных рядов с так называемым дрожанием с распределением случайной величины по равномерному закону и закону Симпсона. Результаты сравнения скейлограмм Results of scale-averaged wavelet power comparison
Графики скейлограмм и СКО для смоделированных процессов показаны на рисунках (см. http://www.swsys.ru/uploaded/image/2017_4/2017-4-dop/10.jpg и http://www.swsys.ru/uploaded/image/ 2017_4/2017-4-dop/11.jpg). Отметим, что при q = 0 исходный временной ряд приобретает регулярный характер, то есть Dtk = const, что является частным случаем. Заключение В статье были предложены алгоритм вейвлет-преобразования неэквидистантных временных ря- дов, а также алгоритм выбора сдвигов данного пре- образования, учитывающий случайный характер интервалов дискретизации. С использованием разработанного комплекса программ были приведены примеры результатов преобразования. Показано, что при преобразовании временных рядов с пропусками наблюдений с увеличением числа пропусков q погрешность оценки коэффициентов преобразования в области высоких частот возрастает. При преобразовании временных рядов с дрожанием данная погрешность значительно меньше, при этом закон Симпсона показывает лучшие результаты по сравнению с равномерным законом распределения. Литература 1 Прохоров С.А., Столбова А.А. Вейвлет-преобразование нерегулярных процессов без восстановления пропущенных отсчетов // Перспективные информационные технологии (ПИТ 2017): тр. Междунар. науч.-технич. конф. Самара, 2017. С. 154–156. 2 Чуи Ч. Введение в вейвлеты; [пер. с англ.]. М.: Мир, 2001. 412 с. 3 Малла С. Вейвлеты в обработке сигналов; [пер. с англ.]. M.: Мир, 2005. 671 с. 4 Ziran Peng, Guojun Wang, Novel ECG eigenvalue detection algorithm based on wavelet transform. BioMed Research Intern., 2017, vol. 2017, 12 p. 5 Ke Liao, Ran Xiao, Jania Gonzalez, Lei Ding. Decoding individual finger movements from one hand using human EEG signals. PLOS ONE, 2014. vol. 9, iss. 1, e85192. 6 Баевский Р.М., Иванов Г.Г., Чирейкин Л.В. [и др.]. Анализ вариабельности сердечного ритма при использовании различных электрокардиографических систем // Вестн. аритмологии. 2002. № 24. С. 65–87. 7 Мандрикова О.В., Заляев Т.Л. Моделирование вариаций космических лучей на основе совмещения кратномасштабных вейвлет-разложений и нейронных сетей переменной структуры // Цифровая обработка сигналов. 2015. № 1. С. 11–16. 8 Прохоров С.А. Прикладной анализ неэквидистантных временных рядов. Самара: Изд-во СГАУ, 2001. 375 с. 9 Витязев В.В. Анализ неравномерных временных рядов. СПб: Изд-во СПб ун-та, 2001. 68 с. 10 Галягин Д.К., Фрик П.Г. Адаптивные вейвлеты (Алгоритм спектрального анализа сигналов) // Математическое моделирование систем и процессов. 1996. № 4. С. 20–28. References
|
Permanent link: http://swsys.ru/index.php?id=4365&lang=en&like=1&page=article |
PDF version article Full issue in PDF (29.80Mb) |
The article was published in issue no. № 4, 2017 [ pp. 668-671 ] |
Perhaps, you might be interested in the following articles of similar topics: