РЕКУРРЕНТНОЕ ПРЕДСТАВЛЕНИЕ БАЗИСА ПЕРЕХОДА ОТ КОЭФФИЦИЕНТОВ РАЗЛОЖЕНИЯ ПО ФУНКЦИЯМ ЧЕБЫШЕВА-ЭРМИТА К ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЮ | Опубликовать статью ВАК, elibrary (НЭБ)

РЕКУРРЕНТНОЕ ПРЕДСТАВЛЕНИЕ БАЗИСА ПЕРЕХОДА ОТ КОЭФФИЦИЕНТОВ РАЗЛОЖЕНИЯ
ПО ФУНКЦИЯМ ЧЕБЫШЕВА-ЭРМИТА К ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЮ

Научная статья

Бочкарев А.В.^1, *, Сайфуллин Р.Т.²

¹ ORCID: 0000-0003-1615-5659;

² ORCID: 0000-0001-7047-3979;

^{1, 2} Самарский Государственный Технический Университет, Самара, Россия

* Корреспондирующий автор (bochkarevpubliciit[at]gmail.com)

Аннотация

Цель работы заключается в формировании рекуррентного представления для базиса, позволяющего по коэффициентам разложения исходного сигнала в базисе функций Чебышева-Эрмита восстановить массив вейвлет-коэффициентов исходного сигнала. В качестве вейвлетов в работе используются производные функции Гаусса m-го порядка. Рекуррентное представление необходимо для сокращения времени вычислений. Рекуррентное представление рассматриваемого базиса формируется на основе свойств полиномов Эрмита, входящих в состав производных функции Гаусса и функций Чебышева-Эрмита. Благодаря использованию полученных формул расчета вейвлет-коэффициентов с учетом описанного алгоритма их применения удается построить быстрые вычислительные алгоритмы обработки сигналов.

Ключевые слова: функции Чебышева-Эрмита, вейвлет-преобразование, вейвлеты Гаусса, разложение сигнала.

RECURRENT REPRESENTATION OF THE TRANSITION BASIS FROM THE EXPANSION COEFFICIENTS BASED ON THE CHEBYSHEV-HERMITE TO THE WAVELET TRANSFORM

Research article

Bochkarev A.V.^1, *, Saifullin R.T.²

¹ ORCID: 0000-0003-1615-5659;

² ORCID: 0000-0001-7047-3979;

^{1, 2} Samara State Technical University, Samara, Russia

* Corresponding author (bochkarevpubliciit[at]gmail.com)

Abstract

The aim of the work is to form a recurrent representation for the basis that allows for reconstructing an array of wavelet coefficients of the original signal from the expansion coefficients of the source signal in the basis of Chebyshev-Hermite functions. The derivatives of the Gauss function of the m order are used as wavelets.  The recurrent representation is necessary to reduce the calculation time. The recurrent representation of the studied basis is formed according to the properties of Hermite polynomials that are part of the derivatives of the Gauss function and the Chebyshev-Hermite functions. Because of the use of the obtained formulas for calculating the wavelet coefficients, taking into account the described algorithm for their application, it is possible to build fast computational algorithms for signal processing.

Keywords: Chebyshev-Hermite functions, wavelet transform, Gaussian wavelets, signal decomposition.

Введение

Одним из подходов к созданию алгоритмов обработки сигналов является кодирование сигнала в базисе функций Чебышева-Эрмита с последующим декодированием по другим, предварительно рассчитанным базисам; причем, в зависимости от выбора базиса возможно получить сам сигнал [1], [2], производную различных порядков [2], вейвлет преобразование [3], [4] и т.п.

Функции Чебышева-Эрмита находят широкое распространение в различных областях науки и техники. Так, их повсеместно применяют в математике, к примеру, в работе [5] изучается вопрос нахождения приближенных решений нелинейных дифференциальных уравнений с использованием спектрального метода, в работах [6] и [7] исследуется тот же вопрос с применением проекционного метода. В работе [8] описано применение данных функций для построения ортогонального банка фильтров, поскольку функции обладают сглаживающим свойством, при их использовании в качестве базиса для кодирования сигналов.

 Вейвлет-анализ является одним из наиболее мощных и гибких средств исследования как одномерных, так и двумерных сигналов. Так, в работе [9] вейвлет-преобразование служит для анализа поверхностей методом проекционных полос с учетом возникающих дифракционных пятен. Также вейвлет-преобразование способствует в работе [10] определению экономического влияния COVID-19 на производство США. В работе [11] рассматривается вопрос существования решения интегрального уравнения Фредгольма с использованием сразу нескольких вейвлетов. Помимо этого, вейвлет-преобразование широко применяется в хроматографии и спектрометрии, к примеру, в работе [12] данное преобразование позволяет авторам дешево, быстро и точно определять содержание дорзоламида и тимолола в человеческой слезе, при этом используется множество различных вейвлетов одновременно.

Вейвлет-преобразование, при использовании множества вейвлетов одновременно, из-за необходимости осуществлять интегрирование зачастую вызывает затруднение при анализе сигналов произвольной формы, что и составляет рассматриваемую в данной работе научную проблему. Следовательно, для разработки современных алгоритмов обработки сигналов и решения уравнений актуальным является вопрос вычисления вейвлет-преобразования сигналов произвольной формы без необходимости интегрирования. В соответствии с этим, задачей настоящего исследования является формирование быстрого вычислительного алгоритма рекуррентного вычисления вейвлет-преобразования сигналов без интегрирования с применением методики кодирования в базисе функций Чебышева-Эрмита.

Базис перехода от коэффициентов разложения по функциям Чебышева-Эрмита к вейвлет-преобразованию

Ранее, в работе [4] была рассмотрена возможность перехода от коэффициентов  кодирования сигнала в базисе функций Чебышева-Эрмита к вейвлет-преобразованию  этого сигнала при использовании в качестве вейвлета производной функции Гаусса m порядка (далее вейвлет-преобразование m порядка). Если сигнал , закодированный в базисе функций  по  функциям (при ) можно восстановить по выражению:

       (1)

где x₀ – величина сдвига функции Чебышева-Эрмита;

γ – коэффициент масштаба функции Чебышева-Эрмита;

то для восстановления вейвлет-преобразования  сигнала  m порядка можно прибегнуть к подобному (1) выражению, в котором  сменяются на функции  базиса перехода к вейвлет-преобразованию:

      (2)

где a – коэффициент масштаба вейвлета;

b – коэффициент сдвига вейвлета;

и сам базис записывается следующим образом:

где .

Базис (3) получен путем интегрирования произведения n-й базисной функции и вейвлета m порядка:

      (5)

Нетрудно заметить, что расчет (4) при высоких порядках m и n затруднителен в силу большого количества математических операций. Следовательно, и сам базис перехода к вейвлет-преобразованию (3) будет требовать все больше времени при увеличении числа используемых базисных функций, либо с увеличением порядка вейвлета. Следовательно, для эффективного применения базиса требуется сократить количество вычислительных процедур.

Формирование рекуррентного представления базиса перехода от коэффициентов разложения по функциям Чебышева-Эрмита к вейвлет-преобразованию

Одним из очевидных способов упростить расчет полинома (4) является его преобразование к рекуррентной форме, что априори возможно, поскольку известны рекуррентные формулы для составляющих данное выражение многочленов Эрмита [9]:

 

– замены, введенные в [4] для упрощения интегрирования.

Для перехода (3) к рекуррентному представлению следует найти произведение (6) и (7). Обозначив данное произведение , а также, для упрощения записи, введя обозначения  и  получим, после раскрытия скобок и приведения подобных членов:

         (8)

Поскольку каждый член (8) представляет собой, равно как и само выражение (8), произведение , допускается ввести также для каждого члена (8) обозначение

Поскольку полином (8) находится под знаком интеграла при расчете вейвлет-преобразования (5), для упрощения интегрирования, за счет дистрибутивности суммы, следует внести множитель t под знак суммы и привести все множители t с различными степенями к одному основанию. Выполнив данные преобразования возможно совершить в выражении выше замену , которой соответствует табличный интеграл [10]:

      (9)

Данная замена необходима [4] для интегрирования в (5).

Чтобы упростить дальнейшую запись выражений, введем следующее обозначение для полинома с внесенным под знак суммы множителем t:

      (10)

Очевидно, правая часть (10) представляет собой  со степенью t под знаком суммы, соответствующей выбранным m, n, к котором добавляется степень множителя . В соответствии с (10) можно выполнить замену в (8), получив при этом:

        (11)

Подставляя  в (3) вместо  с соответствующей заменой  получим:

        (12)

при этом, следуя терминологии, введенной в (11), в дальнейшем можно обозначать , при этом, аналогично (10), в явном виде такое обозначение будет иметь следующий смысл:

      (13)

Сформировать из (12) рекуррентную формулу можно, если подставить в него выражение (11), раскрыть скобки и выделить в виде (12) все возможные отдельные :

 

Полученное выражение (14) хоть и ссылается на собственные члены низших порядков, но содержит также четыре члена , которые невозможно свести к виду (14), поэтому каждый из данных членов должен быть рассчитан отдельно, также рекуррентно, согласно выражению (11). За счет этого теряются преимущества в скорости вычислений, при этом запись базиса в виде (14) куда более громоздкая по сравнению с (12). Поскольку целью перехода к рекуррентной форме являлось упрощение вычислительного процесса, наиболее оптимальным способом формирования базиса перехода к вейвлет-преобразованию с вейвлетами-производными функции Гаусса, является предварительный расчет полиномов (11) рекуррентно с последующей подстановкой результата в (12).

Заключение

Таким образом, в качестве результата работы выступают выражения (14) и (15) задающие искомые выражения рекуррентного вычисления без интегрирования множества вейвлет-преобразований сигнала, закодированного в базисе функций Чебышева-Эрмита. Полученные выражения основываются на рекуррентных свойствах полиномов Эрмита, из-за чего, при предварительном вычислении всех требуемых начальных членов выражения (14), становится возможным без интегрирования сформировать вейвлет-коэффициенты исследуемого сигнала по целому семейству вейвлетов, заданных производной функции Гаусса, к которым, в частности, относится и широко распространенный вейвлет «Мексиканская шляпа».

Финансирование Работа выполнена при финансовой поддержке гранта аспирантов ФГБОУ ВО «СамГТУ».	Funding The paper is carried out with financial support of SamSTU’s postgraduate student grant.
Конфликт интересов Не указан.	Conflict of Interest None declared.