Exponenta Banner
exponenta event banner

imodwpt

Обратное максимальное перекрытие дискретного вейвлет-пакетного преобразования

свернуть все на странице

Синтаксис

xrec = imodwpt (coefs)
xrec = imodwpt (coefs, wname)
xrec = imodwpt (coefs, lo, hi)

Описание

пример

xrec = imodwpt(coefs) возвращает обратное максимальное перекрытие дискретного вейвлет-пакетного преобразования (обратное MODWPT), в xrec. Обратное преобразование для матрицы коэффициентов терминального узла (coefs) получены с использованием modwpt с длиной по умолчанию 18 Fejér-Korovkin ('fk18') вейвлет.

пример

xrec = imodwpt(coefs,wname) возвращает обратный MODWPT с помощью ортогонального фильтра, указанного wname. Этот фильтр должен быть тем же фильтром, что и в modwpt.

пример

xrec = imodwpt(coefs,lo,hi) возвращает обратный MODWPT с помощью фильтра ортогонального масштабирования, loи вейвлет-фильтр, hi.

Примеры

свернуть все

Открыть сценарий в реальном времени

Получение MODWPT формы сигнала ЭКГ и демонстрация идеальной реконструкции с использованием обратного MODWPT.

load wecg;
wpt = modwpt(wecg);
xrec = imodwpt(wpt);
subplot(2,1,1)
plot(wecg);
title('Original ECG Waveform');
subplot(2,1,2)
plot(xrec);
title('Reconstructed ECG Waveform');

Figure contains 2 axes. Axes 1 with title Original ECG Waveform contains an object of type line. Axes 2 with title Reconstructed ECG Waveform contains an object of type line.

Найдите наибольшую абсолютную разницу между исходным сигналом и реконструкцией. Разница составляет порядка 10-11, что демонстрирует идеальную реконструкцию.

max(abs(wecg-xrec'))
ans = 1.7903e-11
Открыть сценарий в реальном времени

Получить данные MODWPT индекса южного колебания с помощью экстремального фазового импульса Daubechies с двумя моментами исчезновения ('db2'). Реконструируйте сигнал с помощью обратного MODWPT.

load soi;
wsoi = modwpt(soi,'db2');
xrec = imodwpt(wsoi,'db2');
Открыть сценарий в реальном времени

Получить данные MODWPT индекса южного колебания с помощью указанных фильтров масштабирования и вейвлетов с экстремальным фазовым вейвлетом Daubechies с двумя моментами исчезновения ('db2'). 

load soi;
[lo,hi] = wfilters('db2');
wpt = modwpt(soi,lo,hi);
xrec = imodwpt(wpt,lo,hi);

Постройте график исходной формы сигнала КНИ и восстановленной формы сигнала.

subplot(2,1,1)
plot(soi)
title('Original SOI Waveform');
subplot(2,1,2)
plot(xrec)
title('Reconstructed SOI Waveform')

Figure contains 2 axes. Axes 1 with title Original SOI Waveform contains an object of type line. Axes 2 with title Reconstructed SOI Waveform contains an object of type line.

Входные аргументы

свернуть все

Коэффициенты терминального узла дерева вейвлет-пакетов, определенные как матрица. Необходимо получить матрицу коэффициентов из modwpt с использованием 'FullTree',false вариант. 'FullTree',false является значением по умолчанию modwpt.

Типы данных: double

Синтезирующий вейвлет-фильтр, используемый для инвертирования MODWPT, заданного как вектор символов или строковый скаляр. Указанный вейвлет должен быть тем же вейвлетом, что и в анализе с modwpt.

Фильтр масштабирования, заданный как действительный вектор четной длины. lo должен быть тем же фильтром масштабирования, что и в анализе с modwpt. Нельзя указать как пару масштабирование-вейвлет-фильтр, так и wname фильтр.

Вейвлет-фильтр, заданный как действительный вектор четной длины. hi должен быть тем же самым вейвлет-фильтром, что и при анализе с modwpt. Нельзя указать как пару масштабирование-вейвлет-фильтр, так и wname фильтр.

Выходные аргументы

свернуть все

Обратное максимальное перекрытие дискретного вейвлет-пакетного преобразования, возвращаемого в виде вектора строки. Обратное преобразование представляет собой восстановленную версию исходного сигнала на основе коэффициентов терминального узла MODWPT. xrec имеет то же количество столбцов, что и входной coefs матрица.

Ссылки

[1] Персиваль, D. B. и А. Т. Уолден. Вейвлет-методы для анализа временных рядов. Кембридж, Великобритания: Cambridge University Press, 2000.

[2] Уолден, А.Т. и А. Контрерас Кристан. «Скорректированное по фазе недекимированное дискретное вейвлет-пакетное преобразование и его применение для интерпретации синхронизации событий». Труды Лондонского королевского общества A. Vol. 454, выпуск 1976, 1998, стр. 2243-2266.

Расширенные возможности

.

См. также

| |

Представлен в R2016a

Документация инструментария вейвлет

Поддержка