ISSNÜ868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2ÜÜ2, том І2, № l, c. 9Ü-95
ОРИГИНАЛЬНЫЕ СТАТЬИ
УДК 621.519
© Д. С. Потехин, Е. П. Тетерин, И. Е. Тарасов
ВЛИЯНИЕ КОЭФФИЦИЕНТОВ И ПРЕДЕЛОВ ИНТЕГРИРОВАНИЯ ВЕЙВЛЕТ-ФУНКЦИИ МОРЛЕ НА ТОЧНОСТЬ РЕЗУЛЬТАТОВ АНАЛИЗА ГАРМОНИЧЕСКИХ СИГНАЛОВ С НЕСТАЦИОНАРНЫМИ ПАРАМЕТРАМИ
В статье рассмотрена зависимость величины интеграла действительной части вейвлет-функции Морле от пределов интегрирования и коэффициента в показателе степени экспоненты. Найдены оптимальные значения этих факторов, при которых исследуемый интеграл стремится к нулю.
ВВЕДЕНИЕ
Создание современных высокоточных приборов предполагает решение целого ряда задач, одна из которых заключается в определении параметров нестационарных сигналов, имеющих различные гармонические составляющие. Для их определения применяют либо фильтрацию сигналов с помощью устройств с избирательной амплитудно-частотной характеристикой, либо спектральное разложение в ряд Фурье, либо методы спектрального анализа, основанные на ином математическом аппарате. И если первый подход часто — аппаратное решение задачи в аналоговой форме со всеми присущими ему достоинствами и недостатками, то последние два предполагают цифровую форму представления сигналов с последующей обработкой числовых массивов средствами вычислительной техники. В настоящее время наибольшее распространение получил Фурье-анализ как хорошо отработанный и испытанный многолетней практикой метод анализа. Однако преобразованию Фурье присущи по крайней мере два существенных недостатка: невозможность развертки во времени амплитудных и фазовых характеристик исследуемых сигналов и наличие эффекта Гиббса, заключающегося в искажении величины спектральной плотности вследствие представления функции конечной длительности бесконечными гармоническим функциями.
От этих недостатков свободен интенсивно развивающийся в последнее время представитель третьего подхода анализа сигналов — вейвлет-анализ [1]. Вейвлет-анализ предусматривает представление исследуемого сигнала в виде функций конечной энергии, быстро убывающей во времени. Спектральная плотность вейвлет-функции находится вычислением интеграла вида
1 x2 ^ f _ t Л
W(t,a) = — Í (p(x)^v\x-|dx ,
aa
где ф(х) — исследуемая функция; ^
вейвлет-функция; х1 и х2 — границы интеграла, приведенные к единичному масштабу; ^ — момент времени, для которого вычисляется значение вейвлет-функции; а — масштаб.
Поскольку величина вычисляемого интеграла существенным образом зависит от вида используемой вейвлет-функции, важной задачей является получение ее аналитического выражения, позволяющего наиболее точно определить параметры исследуемого нестационарного сигнала.
Вейвлет-функцией может быть любая функция конечной энергии, удовлетворяющая условию
x2
JV(x )dx
=о
АНАЛИЗ ПОГРЕШНОСТИ АППРОКСИМАЦИИ СИГНАЛА ВЕЙВЛЕТ-ФУНКЦИЯМИ МОРЛЕ
Одной из широко известных вейвлет-функций является функция Морле, описываемая выражением
w(x )= (<
cos х + i sin X )• e k ,
где k — некоторый постоянный коэффициент.
У вейвлет-функции Морле выделяется реальная (косинусная) составляющая, обозначаемая i//Re, и мнимая (синусная) — 1//1т, а коэффициент к традиционно выбирают равным 50 с пределами интегрирования (1) Xi и х2, равными ±4п Такой вейвлет-
a
2
x
функцией проведены исследования гармонической амплитуду. На рис. 1 показано влияние фазы сиг-функции ф(0=100-со8(ю с целью определить ее нала на погрешность определения амплитуды.
A (V
Рис. 1. Зависимость измеренной амплитуды тестового сигнала от его фазы при коэффициенте k = 50 и границах интегрирования ±4п вейвлет-функции Морле
5 (к)
Рис. 2. Зависимость величины интеграла реальной части вейвлет-функции Морле от коэффициента k (пределы интегрирования ±4п)
Значения интеграла £ от вейвлет-функции Морле, размах изменений измеренной амплитуды и их среднеквадратическое значение при различных пределах интегрирования и значениях параметра k
Пределы интегрирования вейвлет-функции Морле Величина параметра к вейвлет-функции Морле Значение 5 интеграла от вейвлет-функции Морле Размах изменения амплитуды, при А = 100 Среднеквадратическое отклонение амплитуды, при А = 100
±4п 50 5.8 6.0 2.1
±4.04п 50 4.1 — —
±4п 30.011 0.005 0.04 0.01
±4.04п 30.011 0.2 — —
±1.5п 13.284 0.0007 2.1 0.7
±1.515п 13.284 0.06 — —
±3.5п 26.53353 1.3-10-4 0.08 0.03
±3.535п 26.53353 2.010-2 — —
±5.5п 39.485491 2.310-5 0.003 0.001
±5.555п 39.485491 2.210-3 — —
±9.5п 65.020128 6.5-10-6 7.610-6 3.5-10-6
±9.595п 65.020128 5.510-5 — —
Из анализа рис. 1 следует, что полный размах изменения амплитуды составляет 6 единиц, при заданной амплитуде тестового сигнала в 100 единиц. Невысокая точность определения амплитуды тестового сигнала цифровым методом и зависимость этого значения от фазы исследуемого сигнала связаны с неудачным выбором коэффициента k и границ интегрирования, следствием чего является неравенство нулю реальной части щ -функции
х2
JVue (х)dx = S, (3)
xi
где ^Re(x) — реальная или косинусная составляющая вейвлета Морле. При традиционных величинах k, х1 и х2 интеграл (3), взятый численно, колеблется в пределах 4.1—5.8 (см. таблицу), в то время как интеграл мнимой (синусной) части для симметричных пределов всегда равен нулю вследствие нечетности функции.
НАХОЖДЕНИЕ ОПТИМАЛЬНОГО к
Для нахождения оптимальной величины коэффициента k была построена зависимость интеграла реальной составляющей вейвлет-функции Морле от собственно коэффициента k (рис. 2). Коэффициент k дает нулевое значение этой функции в точке 30.011. При этом значение интеграла (3) S находится в пределах 0.005-0.2, что значительно лучше предыдущего случая. Измеренная амплитуда тестового сигнала для уточненного коэффициента k представлена на рис. 3.
Рис. 3 показывает, что полный размах отклонения амплитуды с использованием уточненного коэффициента k составляет уже 0.04. Приведенные зависимости доказывают, что погрешность определения амплитуды гармонического сигнала зависит от значения интеграла (3). Чтобы иметь более полную картину влияния границ интегрирования вейвлет-функции Морле на точность определения амплитуды, необходимо рассмотреть поведение функции cos х вблизи различных границ интегрирования.
НАХОЖДЕНИЕ ОПТИМАЛЬНЫХ ГРАНИЦ ИНТЕГРИРОВАНИЯ
Пределы интегрирования ±4п выбраны не совсем удачно, т.к. при этом косинус "обрывается" на максимуме, а при цифровом представлении исследуемого сигнала пределы интегрирования невозможно выдержать точно (они будут выдержаны с точностью до одного кванта дискретизации по времени). Неточность установки пределов интегрирования и приводит к возникновению гармонических колебаний на графике амплитуды (рис. 1, 3).
A(t)
100.03
99.97 -|_________|____|__________________________________|_____
0 0 001 0.002 0 003 0.004 0 005 ?
Рис. 3. Зависимость измеренной амплитуды тестового сигнала от фазы сигнала при уточненном коэффициенте к
S(k)
Рис. 4. Зависимость величины интегрирования реальной части вейвлет-функции Морле от коэффициента к в пределах интегрирования ±4.5^
Таким образом, можно остановить свой выбор на границах интегрирования, в которых функция у = = cos х имеет значение, близкое к нулю. Такие границы имеют следующие значения: ±1.5я; ±2.5 п; ±3.5п; ±4.5п; ±5.5п и т.д. Зависимость реальной части вейвлет-функции Морле от параметра k при границах интегрирования ±4.5 п видна на рис. 4.
Эта зависимость имеет ярко выраженный минимум, но нулевого значения функция не достигает. Подобное поведение зависимости S(k) наблюдается и для пределов интегрирования ±2.5 п. Для границ ±1.5 я; ±3.5п; ±5.5п зависимость S(k) приведена на рис. 5.
А
Ж)
а
, к
б
т
ад
в
г
Рис. 5. Зависимость Б(А, к) от пределов интегрирования А и коэффициента к (а). Сечения: зависимость Б(к) для пределов интегрирования ±1.5^ (б), ±3.5^ (в), ±5.5^ (г)
Анализируя зависимости, представленные на рис. 5, следует иметь в виду, что конечная точность цифрового представления вейвлет-функции приведет к отклонению от оптимального значения
коэффициента к. В этом случае существенную
дБ ДБ
роль играет отношение — ~-------, которое опреде-
дк Дк
ляет величину погрешности вычисления вейвлет-
функции ЛБ в зависимости от точности представления Лк коэффициента к. Т.к. АБ = — А к, для мини-
Эк
мизации ошибки следует выбирать величину —
дк
возможно меньшей. Если сравнить рис. 4 и 5, г, то можно заключить, что выбор пределов интегрирования ±5.5п обеспечивает лучшую устойчивость к погрешности цифрового представления вейвлет-функции, чем использование пределов ±4.5 п
Величины коэффициента к и интеграла (3) для различных пределов интегрирования приведены в таблице.
Большой интерес представляет зависимость оптимального значения коэффициента к от пределов интегрирования вейвлет-функции Морле (рис. 6).
ЗАВИСИМОСТЬ ОПТИМАЛЬНОГО к ОТ ПРЕДЕЛОВ ИНТЕГРИРОВАНИЯ
На рис. 6 прослеживается прямолинейная зависимость оптимального значения коэффициента к от пределов интегрирования вейвлет-функции Морле, построенная по данным из таблицы. Эта зависимость становится еще более интересной, если заметить, что на данной прямой лежат и уточненное значение коэффициента к для традиционных пределов интегрирования ±4п, и минимумы коэффициента к для границ ±2.5п и ±4.5п Существование такой зависимости позволяет легко переходить от одних границ интегрирования к другим. Такой переход полезен, например, при уточнении коэффициента к при "плавающих" вследствие дискретности цифровых сигналов пределах интегрирования вейвлет-функции Морле.
ЗАКЛЮЧЕНИЕ
В заключение нелишним будет упомянуть о главном недостатке вейвлет-преобразований — большом объеме вычислений. Этот недостаток существенно снижает диапазон их практического применения, а величина границ интегрирования существенным образом сказывается на времени вычислений параметров гармонического сигнала. Поэтому при выборе границ интегрирования необходимо помнить о сохранении оптимального соотношения погрешность/время преобразования. Практический опыт использования подобных методик анализа позволяет рекомендовать пределы
к = 3.789548 + 6.4582414-ж
1 3 5 7 9 (хп) 1
Рис. 6. Зависимость оптимального значения коэффициента к от границ интегрирования вейвлет-функции Морле
интегрирования ±5.5п для вычислений с погрешностью не более 0.1 %. Пределы интегрирования ±1.5п приводят к появлению большей погрешности, а пределы ±7.5 п и ±9.5 п существенно увеличивают вычислительные затраты.
Модифицированная вышеописанным образом вейвлет-функция Морле была использована при обработке сигналов в ряде задач, таких как определение параметров акустической волны в жидких средах, определение периодов сокращения сердечной мышцы, сдвига фаз между двумя гармоническими сигналами и др.
В процессе решения этих задач была отмечена высокая точность определения параметров гармонических сигналов и крайне низкое влияние дестабилизирующих факторов на точность определения этих параметров.
СПИСОК ЛИТЕРАТУРЫ:
1. Петухов А.П. Введение в теорию базисов всплесков. СПб.: Изд-во СПбГТУ, 1999. 132 с.
Ковровская государственная технологическая академия
Материал поступил в редакцию 14.12.2001.
EFFECT OF COEFFICIENTS AND INTEGRATION LIMITS OF THE MORLET WAVELET FUNCTION ON THE ACCURACY OF HARMONIC SIGNAL ANALYSIS WITH NON-STATIONARY PARAMETERS
D. S. Potekhin, E. P. Teterin, I. E. Tarasov
In this paper we discuss the dependence of the real part integral of the Morlet wavelet function on the integration limits and the coefficient in the exponent. Optimal values of these quantities at which the integral in question tends to zero are found.