Обратное непрерывное преобразование вейвлета

icwt функция реализует обратный CWT. Используя icwt требует, чтобы вы получили CWT от cwt.

Поскольку CWT является избыточным преобразованием, нет уникального способа задать обратное. Обратная CWT, реализованная в Wavelet Toolbox™, использует аналитический вейвлет Морса и L1.

Обратная CWT классически представлена в двойно-интегральной форме. Предположим, что у вас есть вейвлет и преобразование Фурье, которое удовлетворяет условию допустимости:

Cψ=|ψ(ω)|2|ω|dω<

Для вейвлетов, удовлетворяющих условию допустимости и функциям конечной энергии, f(t), можно задать обратную CWT как:

f(t)=1Cψab<f(t),ψa,b(t)>ψa,b(t)dbdaa2

где ψa,b(t)=1aψ(tba).

Для анализа вейвлетов и функций, удовлетворяющих следующим условиям, существует единственная интегральная формула для обратной CWT. Этими условиями являются:

  • Анализируемая функция, f(t), является действительной, и анализирующий вейвлет имеет действительное преобразование Фурье.

  • Анализируемая функция, f(t), является вещественной, и преобразование Фурье анализирующего вейвлета имеет поддержку только на множестве неотрицательных частот. Это упоминается как аналитический вейвлет. Функция, преобразование Фурье которой имеет поддержку только на множестве неотрицательных частот, должна быть комплексно оценена.

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

Чтобы мотивировать единственную составную формулу, позвольте ψ1 и ψ2 быть двумя вейвлетами, которые удовлетворяют следующее условие допустимости с двумя вейвлетами:

|ψ1*(ω)||ψ2(ω)||ω|dω<

Задайте константу:

Cψ1,ψ2=ψ1*(ω)ψ2(ω)|ω|dω

Вышеуказанная константа может быть комплексной. Позвольте f(t) и g(t) быть двумя конечными энергетическими функциями. Если условие приемлемости с двумя вейвлетами удовлетворено, выполняется следующее равенство:

Cψ1,ψ2<f,g>=<f,ψ1><g,ψ2>*dbdaa

где <, > обозначает внутреннее произведение, * обозначает комплексный сопряженный, и зависимость

Ключ к одной интегральной формуле для обратной CWT состоит в том, чтобы признать, что условие приемлемости с двумя вейвлетами может быть удовлетворено, даже если один из вейвлетов неприемлем. Иными словами, нет необходимости в том, чтобы обе эти категории, а также, по отдельности, были приемлемы. Можно также ослабить требования дальше, позволив одной из функций и вейвлетов быть распределения. Вначале позволяя g(t) быть функцией дельты Дирака (распределением), а также позволяя

f(t)=2Re{1Cψ1,δ0<f(t),ψ1(t)>daa}

где Re{ } обозначает вещественную часть.

Предыдущее уравнение демонстрирует, что можно восстановить сигнал, суммируя масштабированные коэффициенты CWT по всем шкалам.

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

icwt реализует дискретизированный вариант вышеописанного интеграла.