imodwpt

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

Описание

пример

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 формы волны ECG и продемонстрировайте идеальную реконструкцию с помощью обратной 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);

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

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] Персиваль, Д. Б., и А. Т. Уолден. Вейвлет для анализа временных рядов. Кембридж, Великобритания: Cambridge University Press, 2000.

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

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

.

См. также

| |

Введенный в R2016a