ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2016. № 2 (50)
УДК 681.391:543/545
НЕПРЕРЫВНОЕ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЕ СИГНАЛОВ АНАЛИТИЧЕСКИХ ПРИБОРОВ С ИСПОЛЬЗОВАНИЕМ РАДИАЛЬНОГО БАЗИСА
Р.Т. Сайфуллин
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Предложен метод вейвлет-обработки данных аналитических приборов с использованием радиальных базисных функций (РБФ). Вейвлет-преобразование (ВП) осуществляется в два этапа. На первом этапе производят интерполяцию отсчетов исходного дискретизированного сигнала с помощью РБФ. На втором этапе, зная весовые множители базисных функций, аналитическое выражение непрерывного ВП для конкретной базисной функции и вейвлета, на основе свойства линейности ВП восстанавливают непрерывное ВП исходного сигнала.
Ключевые слова: вейвлет, вейвлет-преобразование, радиальная базисная функция, обработка сигналов аналитических приборов.
Вейвлет-анализ является одним из наиболее мощных и гибких средств исследования и цифровой обработки сигналов. Наряду с задачами фильтрации и сжатия сигналов анализ в базисе вейвлет-функций позволяет решать задачи идентификации, моделирования, аппроксимации стационарных и нестационарных процессов, исследовать вопросы наличия разрывов в производных и т. д.
Непрерывное вейвлет-преобразование (ВП) сигнала /(I) определяется соотношением
1 ¥ ( X — Ь ^
уа Ь) т = а —2 Г т(х)У х Ь
Жу(а,Ь)/ = а 2 \/(х)у(- «X, (1)
где функция у
а
'х — Ьл
а
называется вейвлетом; а, Ь - параметры соответственно
масштаба и сдвига. Множитель а обеспечивает единичную норму для любой
—1 ( х — ь Л
функции уаЬ (х) = а 2 у - .
V а 0
Коэффициенты Шу(а,Ь)/ содержат комбинированную информацию как
об исследуемом вейвлете, так и об анализируемом сигнале. Каждый вейвлет имеет характерные особенности во временном и частотном пространстве, поэтому с помощью разных вейвлетов можно полнее выявить или подчеркнуть те или иные свойства анализируемого сигнала. Таким образом, вейвлет-анализ сигнала благодаря изменению масштаба вейвлет-функций способен выявить различия в характеристиках сигнала на различных шкалах, а посредством сдвига можно
Раухат Талгатович Сайфуллин (д.т.н., проф.), профессор кафедры «Информационно-измерительная техника».
проанализировать свойства сигнала в различных точках на всем исследуемом интервале.
В качестве вейвлетов очень часто используются производные функции Гаус-
2
ап —
у(х) = (-1Г1-—е 2, п = 1,2,3,.... ах
При этом наибольшее применение находят гауссовы вейвлеты небольших порядков:
2
Х
у (х) = -хе 2 ;
- х!
у2(х) = (1 -х2)е 2 ;
- х!
у3(х) = (х3 - 3х)е 2 ;
- х!
у4 (х) = (-х4 + 6х2 - 3)е 2; (2)
2
х
у5(х) = (-х5 + 10х2 - 15х)е 2;
- х!
уб (х) = (-6х6 + 15х4 - 45х2 +15)е~ 2 ;
- х!
у7 (х) = (-х7 + 21х5 - 105х3 +105х)е
2 •
у8(х) = (-х8 + 28х6 - 210х4 + 420х2 - 105)е 2 .
Присутствие экспоненциального множителя в вейвлетах обеспечивает их локальность. Из определения гауссовых вейвлетов следует, что производная от вейвлета уп (х) совпадает (с точностью до знака) с вейвлетом уп+1 (х) :
аухх) = -Уп+,( х).
ах
Численная реализация непрерывного ВП сигналов аналитических приборов во временной области на основе выражения (1) рассмотрена в работе [2]. В данной работе с целью расширения функциональных возможностей вейвлет-обработки сигналов для реализации непрерывного ВП предлагается использовать радиальный базис. При этом вейвлет-преобразование осуществляется в два этапа.
На первом этапе производится интерполяция отсчетов исходного дискрети-зированного сигнала с помощью выбранных радиально-базисных функций [3]. В результате получают разложение исходного сигнала в виде взвешенной суммы базисных функций.
На втором этапе, зная весовые множители базисных функций, полученные на первом этапе, и аналитическое выражение непрерывного ВП для конкретных базисных функций и вейвлета, используя свойство линейности вейвлет-преобразования, восстанавливают ВП исходного сигнала. Таким образом, по коэффициентам разложения исходного сигнала в выбранной системе базисных функций можно достаточно просто восстановить вейвлет-преобразование этого сигнала. Для этого достаточно заранее получить аналитические выражения непрерывного ВП выбранной базисной функции и интересующих вейвлетов.
х2
Рассмотрим более подробно каждый этап.
Этап I. Интерполяция с помощью радиально-базисных функций. Для
выходного сигнала аналитического прибора дискретизированные отсчеты представляются рядом участков, каждый длиной п. Внутри участков используются локальные координаты. Для каждого из участков находится интерполяция отсчетов исходного сигнала непрерывной функцией вида
(х-х )2
п —-
/(х) = ]Гр,е 0
(3)
,=1
где 0 - параметр ширины РБФ-функций; р, - неизвестные коэффициенты (веса). Элементам суммы соответствуют центры РБФ-функций в точках х ,.
Для каждого участка получают систему из п линейных уравнений (число неизвестных р, также равно п):
п
Е<
,=1
-х-гг 202
или в матричном виде:
Рг = /(х1 X ] = 1 П
Е • Р = ^,
( х!- х )2
(4)
где матрица Е состоит из элементов е 20 со столбцами г = 1, п и строками
] = 1, п. Если А - интервал равномерной дискретизации сигнала, то матрица Е будет иметь вид
Е
2 Л
1
202
_ (пА )2 202
202
- (пА)! , 202
[(п 1) А]2
202
[(п 1) А]2 202
Вектор Р = (р1,р2,...ривектор ^ = (/(хД/(х2),.../(хп)У , Т - знак транспонирования. Длина участка п и параметр ширины 0 РБФ функций подбираются экспериментально.
Полученная система уравнений (4) из-за плохой обусловленности матрицы Е не всегда может иметь решение
Р = Е- • ^.
В этом случае необходимо применять методы регуляризации - например,
использовать расчет псевдообратной матрицы Р = Е + • ^ (Е + - псевдообратная матрица) либо регуляризацию по А.Н. Тихонову [4]. В результате исходный дис-кретизированный сигнал может быть представлен соотношением (3).
2
А
е
2
А
1
е
е
1
е
е
Аналитическое вычисление непрерывного ВП базисной функции. Пусть
_ (х_х, )2
используется гауссова базисная функция £( х) = ре 2Ь с параметрами веса р и ширины Ь • Вейвлет-образ этой базисной функции при использовании гауссовых вейвлетов (2) может быть вычислен аналитически и представлен в виде [5]
ЖУП (а Ь)' = —ш+1— Уп(^,
т
где г:
х _ Ь
т
т = <у] а2 +Ь .
Например, формулы вейвлет-коэффициентов второго и четвертого порядков имеют вид:
у[ярРаъ
т
1 х _ ь)2
т
(х-Ь 2т2
%4(а, Ь> =
4жрЬ
а
т
6
' х _ Ьу
т
'х _ ^4
т
_ 3
(х _ь )
2т2
(5)
(6)
Этап 2. Восстановление непрерывного ВП сигнала. Пусть требуется произвести вейвлет-преобразование сигнала с использованием гауссового вейвлета второго и четвертого порядков. В каждой точке х, дискретизированного сигнала
известен весовой множитель р1 базисной функции (получен на первом этапе).
Для конкретной базисной функции и вейвлета имеем аналитическое выражение для непрерывного ВП. Для гауссовой базисной функции и гауссового вейвлета второго и четвертого порядков это соотношения (5) и (6). Используя свойство линейности ВП, записываем требуемые непрерывные ВП для исходного сигнала, представленного в виде (4):
^(а, Ь)/ = £ рг
4яраъ
Щу4 (а, Ь)/ = £ р1
1=1
т
т
(х_Ь 2т2
а
1=1
т
6
х
Ь
л2
т
х
Ь
л4
т
■3
(х, _Ь )2
2т2
е
е
е
е
Таким образом, благодаря использованию заранее полученных аналитических формул для непрерывного вейвлет-преобразования базисных функций и вейвлетов удается построить быстрые вычислительные алгоритмы вейвлет-обработки, что в сочетании с их высокой эффективностью является значительным достоинством метода.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Малашкевич И.А. Вейвлет-анализ сигналов. Теория и практика. - Йошкар-Ола: МарГТУ, 2008. - 224 с.
2. Сайфуллин Р.Т., Наумов А.А. Численная реализация непрерывного вейвлет-преообразования во временной и частотных областях // Информационно-измерительные и управляющие системы: сб. науч. статей.- Самара: Самар. гос. тех. ун-т, 2016. - Вып. 2(14). - С. 90-95.
3. Акимов А.В., Дрюченко М.А., Сирота А.А. Модели и алгоритмы внесения деформирующих искажений на изображениях с использованием радиально-базисных функций // Вестник Воронежского государственного университета. Сер. Системный анализ и информационные технологии. - 2014. - № 1. - С. 130-137.
4. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. - М.: Наука, 1979. - 285 с.
5. Сайфуллин Р.Т., Наумов А.А. Обработка выходных сигналов аналитических приборов на основе вейвлет-преобразования // Информационно-измерительные и управляющие системы: Сб. науч. статей. - Самара: Самар. гос. тех. ун-т, 2015. - Вып. 1 (11). - С. 108-113.
Статья поступила в редакцию 1 апреля 2016 г.
CONTINUOUS WAVELET SIGNALS TRANSFORMING IN ANALYTICAL INSTRUMENTS WITH RADIAL BASIS USING
R.T. Saifullin
Samara State Technical University
244, Molodogvardeyskaya st., Samara, Russia, 443100, Russian Federation
In this paper we propose a method of wavelet data processing in analytical instruments by radial basis function (RBF) using. The wavelet transform (VT) is carried out in two stages. On the first stage the interpolation of the initial signal samples by RBF using is carried. On the second stage is reduced continuously VT of the original signal with using of the linearity property of the wavelet transform and with knowing of the weight coefficients of the basis functions, the analytical expression for a continual VT and wavelet basis function.
Keywords: wavelet, wavelet transform, radial basis function, data acquisition in the analytical instrumentation.
Rauhat T. Saifullin (Dr. Sci. (Techn.)), Professor.