Вейвлет-преобразование (ВП) представляет собой разложение исходного сигнала на сумму взвешенных функций, образующих полный ортогональный базис в некотором конечномерном пространстве [1]. Для традиционного спектрального анализа используется преобразование Фурье – как разложение исходного сигнала по синусоидам разной частоты. Длявейвлет-преобразования вместо синусоид берутся ограниченные по аргументу функции (вейвлеты), а также их сжатые и растянутые копии. Масштаб является параметром ВП.
Дословно «вейвлет» означает “маленькая волна” (от англ. «wavelet»). В отечественной литературе были попытки внедрить термин “всплеск”, но это не увенчалось успехом. В настоящее время ВП широко применяется в задачах распознавания образов при обработке и синтезе различных сигналов, например, речевых, радиолокационных, а также при обработке двумерных сигналов, для упаковки больших объемов информации, сжатия изображений и во многих других информационных технологиях.
Чтобы некоторая функция могла называться вейвлетом, она должна удовлетворять следующим требованиям: должна быть финитна как во временной, так и в частотной области; сдвиги и растяжения данной функции должны составлять полный ортогональный базис; функция должна быстро убывать при t→∞; её среднее значение должно быть равно нулю. ВП описывается интегралом свертки между исходным сигналом S(t) и вейвлетом y(t) для разных значений параметра а, который является масштабом вейвлета:
, (1)
где τ – переменная сдвига.
На практике реальные сигналы, несущие полезную информацию, наблюдаются на фоне различного рода шумов и помех. Существует класс задач, когда имеется априорная информация только о форме принимаемого сигнала. При этом неизвестными являются следующие характеристики: факт наличия сигнала в анализируемой выборке; масштабный параметр и амплитуда сигнала.
Примером предметной области, где существуют подобные задачи, могут служить технические средства охраны (ТСО), использующие различные физические принципы – радиолокационные, сейсмические, телевизионные и т.п. Традиционно в этих ТСО используется пороговое обнаружение с подстраиваемым порогом. Подобный подход прост в реализации, но характеризуется высоким уровнем ошибок: ложных тревог или пропусков сигнала нарушителя.
Если использовать предварительную обработку в виде ВП исходного сигнала, то это позволит существенно улучшить отношение сигнал/шум. Если затем проанализировать результат ВП, то можно не только обнаружить полезный сигнал с низким уровнем ложных тревог, но и классифицировать его по заданным параметрам (рис. 1).
Рис. 1. Структурная схема использования ВП
На рисунке 2 представлены реальные временные диаграммы сигнала с выхода радиоприемника радиолокационного ТСО (РЛТСО) при проходе нарушителя. При движении нарушителя в зоне охраны РЛТСО возникают биения при сложении прямой ЭМ-волны и отраженной от нарушителя, причем в центре зоны охраны частота биений уменьшается.
Рис. 2. Изменение амплитуды ЭМ-волны при движении нарушителя
В статье [2] приведены основные характеристики разработанной модели сигналообразования на выходе радиоприемника РЛТСО, основой которой является косинусоидальная функция с квадратичным фазовым аргументом в виде , зависящим от координаты движения нарушителя x в охраняемой зоне. С учетом колоколообразной огибающей, которая определяется диаграммой направленности антенны радиопередатчика, и пространственно-временного преобразования x = VH t, модель принимаемого сигнала будет описываться функцией:
, (2)
где с [м] – коэффициент масштаба огибающей; VH [м/c] – скорость перемещения нарушителя; α [1/м ] – параметр, имеющий смысл скорости изменения аргумента или пространственной частоты; φ0 [рад] – начальная фаза гармонической функции.
На рис. 3,а показан вид функции (2) с квадратичным изменением фазового аргумента и вид огибающей сигнала на выходе приемника в РЛТСО при использовании двухполупериодного детектирования – (б). Имея модель сигналообразования в РЛТСО можно перейти от порогового обнаружения факта наличия пересечения охраняемой зоны к более тонкому анализу принимаемого сигнала с использованием методов классификации и распознавания образов.
Для реализации ВП в качестве вейвлета необходимо выбрать функцию вида (2). Вейвлет вида (2) называют вейвлетом Френеля. Графически вейвлет Френеля имеет вид, изображенный на рис. 3,а. На стадии анализа общих свойств ВП представляется несущественным учет двухполупериодного детектирования, которое приводит к асимметрии вейвлета Френеля (рис. 3,б), и которое можно учесть в дальнейшем.
Рис. 3
Рассмотрим как изменяется отношение сигнал/шум после ВП. Будем считать, что сигнал S(t) и шум SШ(t) – это электрическое напряжение. Если в выражение (1) для ВП подставить значение относительного сдвига τ = 0, значение коэффициента масштаба a = 1, что означает совпадение масштаба исходного сигнала и вейвлета, а также учесть совпадение формы вейвлета и исходного сигнала с точностью до постоянного множителя q : ψ(t) = qS(t), то получим выражение для ВП:
[В], (3)
где ЕС – энергия сигнала длительностью ТС , выделяемая на резисторе 1 Ом, rE = 1, [1/Ом] – коэффициент, учитывающий размерность при определении энергии сигнала, q = 1, [1/B] – коэффициент пропорциональности, учитывающий, что отклик ВП имеет размерность – В.
Таким образом, максимум ВП пропорционален энергии сигнала. Форма отклика после ВП изменяется. В книге [4] поясняется, что при вычислении интеграла свертки типа (1) в случае равенства исходного сигнала ядру интеграла свертки для относительного сдвига τ = 0 все спектральные составляющие исходного сигнала складываются в фазе, в результате чего отклик достигает своего максимума, определяемого энергией сигнала. Устройство, которое вычисляет интеграл свертки с ядром, равным исходному сигналу, называется оптимальным по критерию максимума отношения сигнал/шум фильтром или согласованным фильтром. Кроме того, в виду симметрии вейвлета, когда ψ(t) = ψ(-t), ВП для a = 1 соответствует автокорреляционной функции.
Принятие решения о наличии на выходе радиоприемного устройства сигнала о проходе нарушителя через охраняемую зону выполняется на фоне шумов и помех (рис. 2). Основным допущением о характере шума является равномерность спектральной плотности мощности шума в пределах полосы пропускания радиоприемника, то есть FШ(ω) = FШ,0, что соответствует модели «белого шума».
Выразим отношение сигнал/шум до ВП через принятые характеристики сигнала и шума. Для сигнала используем максимальное значение Smax(t=0), а для шума – эффективное (среднеквадратическое) значение напряжения на интервале TC ,
, (4)
которое следует выразить через FШ,0. Для этого воспользуемся теоремой Парсеваля и заменим усреднение реализаций шума по времени усреднением по частоте преобразований Фурье ФШ(jω) от реализаций шума.
[В], (5)
где ΔfC = 1/TC – ширина полосы спектра сигнала.
После различного рода преобразований получено выражение для отношения сигнал/шум после ВП:
, (6)
то есть отношение сигнал/шум после ВП определяется энергией сигнала. Далее определим изменение отношения сигнал/шум после ВП:
. (7)
Из (7) следует, что изменение отношения сигнал/шум после ВП не зависит от мощности шума, а зависит от энергии сигнала. С другой стороны, чем меньше амплитуда исходного сигнала по отношению к его энергии, тем эффект от применения ВП больше.
В теории сигналов существует понятие простых и сложных сигналов, которое вполне применимо для ВП. Простыми сигналами называются сигналы, для которых ΔfCTC ≈ 1, а для сложных – ΔfCTC >> 1. Сложные сигналы имеют внутреннюю модуляцию и после оптимальной по критерию максимума отношения сигнал/шум обработки наблюдается сжатие сигнала в m = ΔfCTC раз и возрастание амплитуды в раз [4]. Подобным вейвлетом является вейвлет Френеля (2), так как он имеет внутреннюю квадратичную фазовую модуляцию.
Рис. 4
Рис. 5
На рис. 4 и 5 представлены результаты моделирования ВП в программной среде MathCAD для аддитивной смеси сигнала и шума. При исходном отношении сигнал/шум = 1 после ВП наблюдается существенное улучшение этого параметра [3]. Таким образом, предварительное вейвлет-преобразование в двухпозиционных РЛТСО приводит к увеличению отношения сигнал/шум, что позволяет уменьшить вероятность ошибок или при сохранении характеристик обнаружения существенно увеличить дальность действия РЛТСО.