CC BY
256
УДК 004.67:004.9 Юр Т. В.1, Дубровин В. И.2, Харитонов В. Н.3
1 Аспирант Запорожского национального технического университета 2Канд. техн. наук, профессор Запорожского национального технического университета 3Инженер-конструктор ГП «Ивченко-Прогресс» (г. Запорожье)
ИССЛЕДОВАНИЕ ФИЛЬТРУЮЩИХ СВОЙСТВ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
В работе исследованы фильтрующие свойства непрерывного вейвлет-преобразования. Приведен метод уменьшения влияния краевых эффектов на результаты вейвлет-анализа. Выполнено обоснование выбора типа вейвлета для анализа нестационарных сигналов. Предложен метод построения банка вейвлетных фильтров с необходимыми анализирующими свойствами. Разработанный метод применен к задаче диагностирования подшипников качения по вибрационному сигналу.
Ключевые слова: нестационарный сигнал, вейвлет-фильтрация, банк фильтров, вейвлет Морле.
ВВЕДЕНИЕ
Наиболее широко используемым методом обработки цифровых сигналов является преобразование Фурье. Однако оно имеет ряд недостатков, которые привели к появлению и разработке новых усовершенствованных методов цифрового анализа нестационарных сигналов. Среди недостатков преобразования Фурье можно выделить усреднение характерных особенностей по всей длительности сигнала, что делает нецелесообразным применение данного метода при необходимости анализа изменений сигнала во времени.
На сегодняшний день все большее распространение получает вейвлет-преобразование [1-4]. Этот инструмент позволяет получить частотно-временное представление анализируемого сигнала.
Однако зачастую на практике вейвлет-анализ сигналов выполняется без учета особенностей настройки вейвлетов, которые существенно влияют на результаты преобразования и возможность их правильной интерпретации. Применение вейвлет-преобразования без учета частотных характеристик вейвлета может давать искаженную информацию об особенностях анализируемого сигнала. Такой подход может привести к неуниверсальности разрабатываемых методов в силу того, что при изменении характеристик исследуемого сигнала применение одних и тех же настроек метода может давать совершенно непредсказуемые результаты.
ПОСТАНОВКА ЗАДАЧИ
Для успешного использования вейвлет-преобразования необходим тщательный анализ параметров мате© Юр Т. В., Дубровин В. И., Харитонов В. Н., 2010
ринского вейвлета. Это позволит гибко управлять настройкой вейвлета для получения необходимых анализирующих свойств вейвлет-преобразования и правильно интерпретировать получаемые результаты.
Целью данной работы является исследование фильтрующих возможностей непрерывного вейвлет-преобразования, обоснование выбора материнского вейвлета и разработка метода построения анализирующего банка вейвлетных фильтров.
ПРЕДСТАВЛЕНИЕ ВЕЙВЛЕТ-
ПРЕОБРАЗОВАНИЯ В ВИДЕ ОПЕРАЦИИ
ФИЛЬТРАЦИИ
Непрерывное вейвлет-преобразование представляет собой разложение анализируемого сигнала в базисе анализирующей вейвлетной функции.
Базис вейвлет-преобразования строится путем масштабных преобразований и переносов материнского вейвлета у(*) с непрерывными значениями базисных параметров - масштабного коэффициента а и параметра сдвига Ь:
Уа, ь(*) = \а\72у(*-0), а, Ь е Я, (1)
где уа, Ь (*) - дочерний вейвлет базиса.
В интегральной форме вейвлет-преобразование непрерывного сигнала 5(*) может быть выражено:
1 м
Ж(а, Ь ) = Г 5(*)у*ь (*)Л, (2)
л/И
где * - оператор комплексного сопряжения.
Используя равенство Парсеваля, которое описывает связь между функциями и их образами Фурье, выражение (2) может быть записано в альтернативной форме:
Ж(а, Ь) = „ТаЕ-{Б(/)¥*(а/)}, (3)
где Е 1 - оператор обратного преобразования Фурье; ¥(о/) - преобразование Фурье дочернего вейвлета; Б(/) - преобразование Фурье анализируемого сигнала; / - частота; а - масштабный коэффициент.
Такое представление вейлет-преобразования позволяет значительно ускорить вычисление вейвлет-коэффициентов за счет использования быстрого преобразования Фурье (БПФ). К тому же, преобразование Фурье материнского вейвлета обычно известно в аналитическом виде и, следовательно, количество вычислений сокращается.
Выражение (3) также показывает, что вейвлет-преобразование может быть рассмотрено как специальная операция фильтрации с вейвлетом в качестве ядра фильтра.
ЧАСТОТНО-ВРЕМЕННЫЕ
ХАРАКТЕРИСТИКИ
ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
Рассмотрим частотно-временные характеристики непрерывного вейвлет-преобразования [5, 6].
Выражение (2) представляет вейвлет-преобразование как средство анализа во временно-масштабной области. В таком виде интерпретация масштабов и их соответствие частотным диапазонам могут быть затруднены. Если же рассматривать вейвлет-преобразование в частотной области, частотный образ материнского вейвлета соответствует полосовому фильтру, характеристики которого определяются типом и настройками вейвлета.
Пусть вейвлет у(*) и его образ Фурье ¥(/) являются функциями-окнами с центрами *0 и /0, радиусами 8/ и 8* соответственно.
На рис. 1 изображено частотно-временное разрешение непрерывного вейвлет-преобразования. Оно иллюстрирует анализирующие возможности вейвлетов на различных масштабах.
Благодаря свойствам смещения и масштабирования материнского вейвлета каждый дочерний вейвлет на масштабе а будет иметь ширину окна, равную 28// а, и центральную частоту, равную /0/ а. Их отношение не зависит от масштаба преобразования и является постоянной величиной 28///0, значение которой задается параметрами материнского вейвлета.
Рис. 1. Частотно-временное разрешение вейвлет-преобразования
Таким образом, вейлет-преобразование также может быть рассмотрено как операция фильтрации с постоянной добротностью (с постоянной 0, при ко -торой отношение среднеквадратической ширины полосы пропускания к центральной частоте фильтра является постоянной величиной.
Если рассмотреть частотно-временное окно, характеризующее разрешение вейвлет-анализа, можно отметить, что его площадь постоянна и равна 4 5/ 5,. При масштабировании материнского вейвлета радиус окна дочернего вейвлета во временной области увеличивается прямо пропорционально масштабу (5а), а в частотной радиус окна уменьшается обратно пропорционально масштабу (5/ а). Таким образом, вейвлет-преобразование предоставляет переменное разрешение в частотно-временной плоскости, которое дает ему большое преимущество при анализе сигналов, так как быстрые вариации сигналов (высокочастотные характеристики) хорошо локализованы, а для выявления медленно меняющихся характеристик необходимо большое временное окно.
Согласно принципу неопределенности Гейзенбер -га между временной и частотной локализациями площадь частотно-временного окна вейвлет-преобразования ограничена 5/5, >(4п)-1, что определяет ограничение разрешающей способности вейвлет-анализа [1].
ВЫБОР МАТЕРИНСКОГО ВЕЙВЛЕТА
Среди существующих вейвлет-функций, используемых в непрерывном вейвлет-анализе, вейвлет Мор-ле и вейвлет Мексиканская шляпа более всего по форме подобны импульсным составляющим неста-
ционарных сигналов, что делает их наиболее подходящим средством анализа такого вида сигналов.
Оба эти вейвлета являются экспоненциально убывающими как во временной, так и в частотной области. Принципиальное различие между ними заключается в следующем: вейвлет Морле имеет узкий спектральный диапазон и продолжителен во временной области, тогда как вейвлет Мексиканская шляпа более локализован во временной области за счет более широкого частотного представления [3]. В силу принципа неопределенности Гейзенберга невозможно получить идеальную локализацию одновременно в частотной и временной областях.
При выборе анализирующего вейвлета для фильтрации сигналов предпочтение было отдано вейвлету Морле, который имеет узкий образ Фурье. Присутствие доминирующей частоты позволяет варьировать избирательность вейвлета Морле в частотной области.
Кроме того, в качестве анализирующего был выбран комплексный тип вейвлета Морле. Его преобразование Фурье равно нулю для отрицательных частот, что позволяет разделить фазовые и амплитудные компоненты сигнала при выполнении вейвлет-преобразования.
Во временной области комплексный вейвлет Мор-ле представляет собой комплексную экспоненту, модулируемую функцией Гаусса:
.(*) = ^е-0 * е '
л/П
(4)
где * - время; о и /0 - параметры вейвлета, задающие его форму.
В частотной области вейвлет Морле имеет форму Гауссового окна с центральной частотой /0 и шириной о:
¥(/) = ¥ *(/) = е
(5)
где ¥(/) - преобразование Фурье вейвлета Морле; /0 - центральная частота анализирующего вейвлета; о - ширина анализируемого диапазона частот; * означает комплексное сопряжение.
Так как ¥( /) является действительным числом, то ¥( /) = ¥*( /).
Таким образом, частотный диапазон, покрываемый окном вейвлета Морле, практически ограничен интервалом [/0 - о/2, /0 + о/2 ].
На рис. 2 приведен пример представления вейвлета Морле во временной и частотной областях.
Строго говоря, вейвлет Морле не удовлетворяет условию допустимости для вейвлетов:
с. = {' М л/ <•
/
(6)
Это условие также подразумевает, что вейвлет не должен иметь нулевой частотной компоненты ¥( 0) = = 0 или же, другими словами, вейвлет должен иметь нулевое среднее.
Это означает, что вейвлет Морле не обладает компактным носителем. Однако наибольшая часть его энергии сосредоточена на интервале, ширина которого определяется его полосой пропускания, и на практике при /0» 0 вейвлет Морле может быть использован с минимальной погрешностью [3].
КРАЕВЫЕ ЭФФЕКТЫ И ОБЛАСТЬ
ДОСТОВЕРНОСТИ ДЛЯ ВЕЙВЛЕТ-
ПРЕОБРАЗОВАНИЯ
Рассмотрим механизм возникновения граничного эффекта при выполнении вейвлет-преобразования сигнала [7, 8].
Значение анализируемой функции в точке *0 влияет на значения коэффициентов преобразования в растущих с масштабом временных диапазонах, в так называемом угле влияния (рис. 3).
Если анализирующий вейвлет .(*) хорошо локализован в интервале 28* для а =1, то коэффициенты вейвлетного спектра, соответствующие моменту времени *0 будут содержаться в конусе влияния, ограниченном прямыми а = *0 - *8* и а = *0 + *8* (рис. 3, а). Коэффициент Ш(а0, Ь0) в свою очередь зависит от значений ряда из определяемого тем же углом влияния временного диапазона около значения Ь0 (рис. 3, б).
Диапазон этот тем больше, чем больше анализируемый масштаб а. Таким образом, высокочастотная
а)
Рис. 3. Угол влияния вейвлет-преобразования:
а - влияние точки анализируемой функции на значения вейвлет-коэффициентов; б - зависимость коэффициента Ш(ао, Ь0) от значений анализируемой функции (пунктиром обозначена область достоверности)
(мелкомасштабная) информация вычисляется на основе малых по длительности отрезков сигнала, а низкочастотная - на основе больших отрезков сигнала.
Максимальный угол влияния отмечает область достоверности - область, вне которой коэффициенты Ш(а, Ь) вычисляются уже по отрезкам, выходящим за границы ряда, и поэтому амплитуды вейвлет-коэффициентов носят искаженный характер. Заметим, что при нелинейном изменении масштабов на графиках вейвлет-преобразования угол влияния примет нелинейный вид.
Таким образом возникает граничный эффект - область, в которой коэффициенты вейвлет-преобразования искажаются разрывной природой границы сигнала. Пример граничного эффекта, возникающего при вейвлет-преобразовании, приведен на рис. 4.
Как видно из рис. 4, амплитуды вейвлет-преобразования в области краевых эффектов сравнимы с амплитудами коэффициентов, определяющих особую точку сигнала.
При выполнении вейвлет-преобразования конечной временной реализации длины N в частотной области ошибки в расчете коэффициентов Ш(а, Ь) обусловлены природой преобразования Фурье, которое предполагает, что сигнал 5 (*) периодичен и имеет период длиной N.
Для уменьшения ошибки граничных эффектов ряд дополняется с учетом его поведения (средним значением, известным временным ходом и т. д.) [3].
В данной работе для анализа вибрационных сигналов предлагается использовать дополнение временного ряда нулями до длины 2N. Такой подход приводит к появлению сильной неоднородности на границах исходного временного ряда. Однако возмущения, вносимые этой неоднородностью, лежат в области очень больших временных масштабов, в то время как спектр исходной неоднородности (которая фор-
1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
б)
Рис. 4. Иллюстрация граничного эффекта, возникающего при вейлет-анализе: а - пример анализируемого сигнала; б - вейвлет-пребразование сигнала
Рис. 5. Вейвлет-преобразование сигнала 4, а после дополнения нулями
мально существенно меньше внесенной) лежал бы в области масштабов анализируемого сигнала. Так что внесение такой неоднородности приводит к уменьшению амплитуды Ж(а, Ь) около границ исходного ряда. Очевидно, чем больше будет длина дополнения, заполненного нулями, тем меньше влияние краевых эффектов. Использование для анализа ряда, в котором длина последовательности нулей равна длине исходного сигнала, является оптимальной в смысле баланса скорости вычислений, затрат машинной памяти и точности осуществления вейв-лет-преобра-зования вблизи границ исходного временного ряда [8] (рис. 5).
На рис. 5 показано вейвлет-преобразование сигнала, приведенного на рис. 4, а, с использованием дополнения ряда нулями. Как видно, если сравнить рис. 4, б и рис. 5, влияние краевых эффектов значительно уменьшилось. Амплитуды в области краевых эффектов значительно ниже амплитуд коэффициентов, определяющих особую точку.
РАСЧЕТ ОБЛАСТИ ДОСТОВЕРНОСТИ ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ, ВЫПОЛНЯЕМОГО ПРИ ПОМОЩИ ВЕЙВЛЕТА МОРЛЕ
Пусть вейвлет Морле .(*) представляет собой функцию-окно с радиусом 8*. Тогда каждый дочерний вейвлет на масштабе а будет иметь радиус окна, равный 8,а. Как видно из рис. 3, б, если точка Ж(а, Ь) расположена от границы ряда на расстоянии меньшем, чем 8*а, то для ее расчета необходимы значения, выходящие за пределы ряда. Таким образом, количество точек от краев сигнала, которые носят искаженный характер, может быть задано следующим образом:
п = 8,а*/л, (7)
где / - частота дискретизации сигнала.
Так как вейвлет Морле является комплексной экспонентой, модулируемой функцией Гаусса (см. выр. 4), его ширина во временной области может быть определена при помощи правила трех сигм как:
о72'
(В)
Используя выражения (7)-(8), получим 2 области: область достоверности вейвлет-преобразования и область, в которой нельзя пренебречь наличием краевых эффектов. Предлагается отбрасывать данную область краевых эффектов при дальнейшей манипуляции с рассчитанными коэффициентами вейвлет-преобразования.
На рис. 6 приведен пример вейвлет-преобразования 4, б и 5 с рассчитанным максимальным углом влияния.
Заметим, что с уменьшением количества точек после отсекания области, содержащей краевые
1000 2000 3000 4000 5000 Е000 7000 8000 Э000 10000
а)
Рис. 6. Вейвлет-преобразование анализируемого сигнала (рис. 4, а) с обозначенной областью достоверности: а - без дополнения сигнала нулями; б - с дополнением сигнала нулями
i6i
эффекты, падает разрешение по частоте для спектра вейвлет-преобразования:
А/ =
1
(9)
(N - 2*п /
где N - количество отсчетов сигнала; / - частота дискретизации сигнала; п - длина области краевых эффектов.
МЕТОД ПОСТРОЕНИЯ
АНАЛИЗИРУЮЩЕГО БАНКА ФИЛЬТРОВ
Для разделения различных частотных компонент сигналов используются банки фильтров. Как было отмечено, вейвлет-преобразование сигналов можно рассматривать как фильтрацию с постоянной относительной полосой частот. Поэтому по аналогии с обычными банками (гребенками) из 1/п-октавных полосовых фильтров может быть построен вейвлетный банк фильтров.
Для построения анализирующего банка фильтров на основе комплексного вейвлета Морле предлагается использовать следующую стратегию.
Добротность вейвлета Морле в качестве фильтра определяется как
= 0. (1°)
/0
При использовании набора масштабов а1 набор соответствующих центральных частот / и параметров ширины диапазона частот ог- комплексного вейвлета Морле могут быть получены с помощью константы 2Ш:
(11)
Задавая исходные параметры вейвлета /0 и о, а также набор масштабов для расчета непрерывного вейвлет-преобразования, можно построить банк фильтров с необходимыми анализирующими параметрами.
Для /0 максимальным примем значение /0 = 0, 8/уд, где Nуд - частота Найквиста сигнала.
Такое значение не позволит первому фильтру в банке выйти за границу /^
В качестве основания для расчета масштабов вейвлет-преобразования выберем а1 = (21/к)', что позволит строить банк фильтров с к фильтрами на октаву.
Параметр добротности выберем таким образом, чтобы ширина диапазона каждого фильтра равнялась расстоянию между фильтрами
Яш = 21/к -1.
(12)
При помощи данного параметра можно варьировать ширину полосы пропускания фильтров и степень перекрытия фильтров в банке.
После определения обозначенных параметров при помощи выражения (5) может быть построен банк фильтров.
На рис. 7, а показан банк фильтров, полученный при помощи данной процедуры, а на рис. 7, б - банк фильтров с добротностью, увеличенной в два раза.
Используя выражение (7), производится вейвлет-фильтрация анализируемого сигнала.
Так как используемый вейвлет Морле является комплексным, то получаемые вейвлет-коэффициенты для каждого масштаба/фильтра также являются комплексными значениями. Реальная часть вейвлет-коэффициентов Шге(а, Ь) представляет собой отфильтрованный сигнал, мнимая часть Шт(а, Ь) является фазой отфильтрованного сигнала.
Для анализа особенностей сигнала также предлагается использовать спектр огибающей вейвлет-преобразования, который позволяет демодулировать сигнал в необходимом частотном диапазоне:
ЕШ(а, Ь) = [ Шге(а, Ь)2 + Шт(а, Ь)2Р2. (13)
Правильная настройка банка фильтров и учет области достоверности вейвлет-преобразования позволяет проанализировать особенности сигнала наилучшим образом.
а) б)
Рис. 7. Банк вейвлет-фильтров с тремя фильтрами на октаву, покрывающий 4 октавы: а - при добротности = 21/к - 1; б - при добротности, увеличенной в два раза: = 2 • (21/к - 1)
В качестве признаков для выбора оптимального фильтра из построенного банка предлагается использовать различные статистические параметры. Например, эксцесс коэффициентов вейвлет-преобразования может быть использован для получения отфильтрованного сигнала с наибольшей импульсно-стью [9].
При настройке вейвлетного банка фильтров следует также учитывать требования к необходимой точности обработки сигнала. Достижение высокой точности и детальности описания характеристик сигнала приводит к увеличению длительности обработки сигнала. Поэтому одной из главных задач является рациональная настройка метода обработки сигнала, чтобы с приемлемой скоростью получить требуемое описание сигнала.
Для примера в табл. 1 приведено время обработки сигналов при помощи банков вейвлетных фильтров в зависимости от длины сигнала и характеристик банка фильтров. Приведенные значения получены для конкретного программно-аппаратного комплекса.
Таблица 1. Время обработки сигнала при помощи банка вейвлетных фильтров в зависимости от длины сигнала и характеристик банка фильтров
ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ ПРЕДЛОЖЕННОГО МЕТОДА К ЗАДАЧЕ ДИАГНОСТИКИ ПОДШИПНИКОВ КАЧЕНИЯ
В качестве примера для демонстрации фильтрующих возможностей вейвлет-преобразования были использованы вибрационные сигналы тестовой установки с дефектными подшипниками качения [10]. На установку помещали подшипник с точечным дефектом на внешней дорожке качения и подшипник с точечным дефектом на внутренней дорожке. Дефекты размером 0,007 дюйма нанесены при помощи эрозионного станка.
Скорость вращения вала установки примерно равна / = 29 Гц. Частота дискретизации сигнала 48 кГц. Подшипник имеет следующие параметры: количество шариков г = 9, угол контакта шарика с дорожка-
ми ф = 0°, диаметр шарика d = 7,95 мм, диаметр сепаратора D = 38,5 мм.
Таким образом, может быть рассчитана характеристическая частота перекатывания тел качения по внешней дорожке:
Ebpfiq = 0,5zf[ 1~D^) = 104 Гц
и частота перекатывания тел качения по внешней дорожке:
Fbpfi = 0,5zf[ 1 + ■^-С08ф) = 158 Гц.
На рис. 8 приведены результаты обработки сигнала с подшипником с дефектом на внешней дорожке качения, на рис. 9 - с дефектом на внутренней дорожке качения.
Масштабу 7 соответствует центральная частота фильтра 3809,8 Гц с шириной полосы пропускания 990,23 Гц, а масштабу 8 - 3023,8 Гц и 785,95 Гц соответственно.
Рис. 8 и 9 иллюстрируют эффективность фильтрации нестационарных сигналов при помощи комплексного вейвлета Морле.
ВЫВОДЫ
В работе решена актуальная задача исследования фильтрующих свойств непрерывного вейвлет-преобразования. Приведен метод уменьшения влияния краевых эффектов на результаты вейвлет анализа. Выполнено обоснование выбора типа вейвлета для анализа нестационарных сигналов. Предложен метод настройки материнского вейвлета Морле и построения банка фильтров с необходимыми анализирующими свойствами. Предложенный метод эффективно применен к задаче диагностирования подшипников качения по вибрационному сигналу.
В качестве направления дальнейшего исследования может быть выделена задача построения вектора характеристических признаков анализируемого сигнала из получаемых вейвлет-коэффициентов с целью постановки диагноза и принятия решения о состоянии исследуемого объекта.
СПИСОК ЛИТЕРАТУРЫ
1. Чуи К. Введение в вейвлеты / К. Чуи. - М. : Мир, 2001.- 412 с.
2. Малла С. Вэйвлеты в обработке сигналов / С. Малла. -М. : Мир, 2005. - 671 с.
3. Addison P. S. Illustrated wavelet transform handbook. Introductory Theory and Applications in Science, Engineering, Medicine and Finance / Paul S. Addison. - Bristol : Institute of physics publishing, 2002. - 400 p.
4. Дьяконов В. П. Вейвлеты. От теории к практике /
В. П. Дьяконов. - Изд. 2-е, перераб. и доп. - М. : Со-лон-Пресс, 2004. - 400 с.
Длина сигнала, отсчеты Время обработки, с
6 фильтров в банке 12 фильтров в банке 24 фильтра в банке
1024 0,0156 0,0313 0,0469
2048 0,0313 0,0469 0,0781
4096 0,0313 0,0625 0,1406
8192 0,0496 0,0938 0,2188
16384 0,0781 0,1719 0,4219
32768 0,1875 0,4219 1,0469
виброускорение,д
50 100 150 200 250 300 350 400 450 Гц ^
Рис. 8. Обработка сигнала с дефектом на внешней дорожке качения: а - исходный сигнал; б - спектр необработанного сигнала; в - вейвлет-преобразование сигнала; г - спектр огибающей вейвлет-коэффициентов на масштабе номер 8
Рис. 9. Обработка сигнала с дефектом на внутренней дорожке качения: а - исходный сигнал; б - спектр необработанного сигнала; в - вейвлет-преобразование сигнала; г - спектр огибающей вейвлет-коэффициентов на масштабе номер 8
З. Torrence C. A Practical Guide to Wavelet Analysis /
C. Torrence, G. P. Compo // Bulletin of the American Meteorological Society. - 199В. - Vol. 79. - Pp. 61-7В.
6. Najmi A. The Continuous Wavelet Transform and Variable Resolution Time-Frequency Analysis / A. Najmi, J. Sadowsky. // JOHNS HOPKINS APL TECHNICAL DIGEST. - 1997. - Vol. 1В, No. 1. - Pp.m-140.
7. Астафьева Н. Вейвлет анализ: основы теории и примеры применения / Н. Астафьева // Успехи физических наук. - 1996. - Т. 166, № 11. - C. 114З-1170.
В. Короновский А. А. Непрерывный вейвлетный анализ и
его приложения / А. А. Короновский, А. Е. Храмов. -M. : ФИЗMАТЛИТ, 200З. - 176 с.
9. Shir L. Rolling element bearing fault detection using improved envelope analysis / L. Shir, B. Randall, J. Antoni // IMECHE conference transactions. - 2004. -Vol 2. - P. З01-З12.
10. Case Western Reserve University Bearing Data Center. Seeded Fault Test Data [Электронный ресурс]. - Режим доступа: http://www. eecs. case. edu/ laboratory/be-aring/welcome_overview. htm.
Надійшла 19.04.2010
Юр Т. В., Дубровін В. І., Харитонов В. Н.
ДОСЛІДЖЕННЯ ФІЛЬТРУЮЧИХ ВЛАСТИВОСТЕЙ
ВЕЙВЛЕТ-ПЕРЕТВОРЕННЯ
В роботі досліджено фільтруючі властивості безперервного вейвлет-перетворення. Наведено метод зменшення
впливу крайових ефектів на результати вейвлет-аналізу. Виконано обґрунтування вибору типу вейвлету для аналізу нестаціонарних сигналів. Запропоновано метод побудови банку вейвлетних фільтрів з необхідними аналізуючими властивостями. Розроблений метод застосовано до задачі діагностування підшипників кочення за вібраційним сигналом.
Ключові слова: нестаціонарний сигнал, вейвлет-
фільтрація, банк фільтрів, вейвлет Морле.
Yur T. V., Dubrovin V. I., Kharitonov V. N.
ANALYSIS OF WAVELET TRANSFORM FILTERING PROPERTIES
Filtering properties of continuous wavelet transform are investigated in the paper. The method that decreases the influence of edge effects on wavelet analysis results is described. The choice of wavelet for non-stationary signals analysis is grounded. The method of construction of a filter bank with required analyzing properties is proposed. The given method is applied to the problem of rolling bearing diagnostics by vibration signals.
Key words: non-stationary signal, wavelet filtering, filter bank, Morlet wavelet.
