imodwt

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

Синтаксис

xrec = imodwt(w)
xrec = imodwt(w,wname)
xrec = imodwt(w,Lo,Hi)
xrec = imodwt(___,lev)
xrec = imodwt(___,'reflection')

Описание

пример

xrec = imodwt(w) возвращает в xrec восстановленную версию сигнала. Восстановленный сигнал основан на w, коэффициентах максимального перекрытия дискретного вейвлета преобразовывает (MODWT) и на уровне реконструкции, который значения по умолчанию обнулить.

пример

xrec = imodwt(w,wname) восстанавливает сигнал с помощью wname, ортогонального вейвлета. wname должен быть тем же вейвлетом, используемым, чтобы анализировать вход сигнала к modwt. Реконструкция до уровня 0, который является совершенной реконструкцией исходного сигнала.

пример

xrec = imodwt(w,Lo,Hi) восстанавливает сигнал с помощью ортогонального фильтра масштабирования, Lo и вейвлет фильтруют Hi. Lo и фильтры Hi должны быть теми же фильтрами, используемыми, чтобы анализировать вход сигнала к modwt. Реконструкция до уровня 0, который является совершенной реконструкцией исходного сигнала.

пример

xrec = imodwt(___,lev) восстанавливает сигнал до уровня lev. xrec является проекцией на масштабирующийся пробел на уровне lev.

пример

xrec = imodwt(___,'reflection') использует отражательное граничное условие в реконструкции. Если вы задаете 'reflection', imodwt принимает, что продолжительность исходной длины сигнала является одной половиной количества столбцов во входной матрице коэффициентов. По умолчанию imodwt принимает периодическое расширение сигнала на контуре.

Примеры

свернуть все

Получите MODWT ECG, сигнализируют и демонстрируют совершенную реконструкцию.

Загрузите данные сигнала ECG и получите MODWT.

load wecg;

Получите MODWT и обратный MODWT.

w = modwt(wecg);
xrec = imodwt(w);

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

norm(abs(xrec'-wecg),Inf)
ans = 2.3253e-12

Получите MODWT данных об обменном курсе немецкой марки-доллара США и продемонстрируйте совершенную реконструкцию.

Загрузите данные об обменном курсе немецкой марки-доллара США.

load DM_USD;

Получите MODWT и Обратный MODWT использование вейвлета 'db2'.

wdm = modwt(DM_USD,'db2');
xrec = imodwt(wdm,'db2');

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

norm(abs(xrec'-DM_USD),Inf)
ans = 1.6370e-13

Получите MODWT сигнала ECG использование фильтров Фейера-Коровкина.

Загрузите данные о ECG.

load wecg;

Создайте фильтры Фейера-Коровкина с 8 коэффициентами.

[Lo,Hi] = wfilters('fk8');

Получите MODWT и обратный MODWT.

wtecg = modwt(wecg,Lo,Hi);
xrec = imodwt(wtecg,Lo,Hi);

Отобразите на графике исходные данные и реконструкцию.

subplot(2,1,1)
plot(wecg)
title('ECG Signal');
subplot(2,1,2)
plot(xrec)
title('Reconstruction')

Получите MODWT сигнала ECG вниз к максимальному уровню и получите проекцию сигнала ECG на масштабирующийся пробел на уровне 3.

Загрузите данные о ECG.

load wecg;

Получите MODWT.

wtecg = modwt(wecg);

Получите проекцию сигнала ECG на V3, масштабирующийся пробел на уровне три при помощи функции imodwt.

v3proj = imodwt(wtecg,3);

Постройте исходный сигнал и проекцию.

subplot(2,1,1)
plot(wecg)
title('Original Signal')
subplot(2,1,2)
plot(v3proj)
title('Projection onto V3')

Обратите внимание на то, что характеристика скачков волн R в ECG отсутствует в V3 приближение. Вы видите недостающие детали путем исследования коэффициентов вейвлета на уровне три.

Постройте уровень три коэффициента вейвлета.

figure
plot(wtecg(3,:))
title('Level-Three Wavelet Coefficients')

Получите обратный MODWT использование отражательной обработки контура для южных Индексных данных о Колебании. Период выборки является одним днем. imodwt с опцией 'reflection' принимает, что входная матрица, которая является modwt вывод, является дважды продолжительностью исходной длины сигнала. отражение imodwt граничная обработка сокращает количество вейвлета и масштабных коэффициентов в каждой шкале наполовину.

load soi;
wsoi = modwt(soi,4,'reflection');
xrecsoi = imodwt(wsoi,'reflection');

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

norm(abs(xrecsoi'-soi),Inf)
ans = 1.6419e-11

Входные параметры

свернуть все

MODWT преобразовывают, заданный как матрица размера L +1-by-N. w является вывод modwt, который является MODWT N - входной сигнал точки вниз, чтобы выровнять L. По умолчанию imodwt принимает, что вы получили MODWT использование вейвлета 'sym4' с периодической граничной обработкой.

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

Вейвлет синтеза, заданный как одно из следующего:

  • 'haar' — Вейвлет Хаара

  • Вейвлет Daubechies фазы 'dbN' — Extremal с N, исчезающим моменты, где N является положительным целым числом от 1 до 45.

  • 'symN' — Вейвлет Symlets с N, исчезающим моменты, где N является положительным целым числом от 2 до 45.

  • 'coifN' — Вейвлет Coiflets с N, исчезающим моменты, где N является положительным целым числом от 1 до 5.

  • 'fkN' — Вейвлет Fejér-Korovkin с коэффициентами N, где N = 4, 6, 8, 14, 18 и 22.

Вейвлет синтеза должен быть тем же вейвлетом, используемым в анализе с modwt.

Масштабирование фильтра, заданного как ровная длина вектор с действительным знаком. Можно задать Lo, только если вы не задаете wname. Lo должен быть тем же фильтром масштабирования, используемым в анализе с modwt.

Фильтр вейвлета, заданный как ровная длина вектор с действительным знаком. Можно задать Hi, только если вы не задаете wname. Hi должен быть тем же фильтром вейвлета, используемым в анализе с modwt.

Уровень реконструкции, заданный как неотрицательное целое число между 0 и size(w,1)-2. Уровень должен быть меньше, чем уровень раньше получал w из modwt. Если lev 0, и вы не изменяете коэффициенты, imodwt производит совершенную реконструкцию сигнала.

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

свернуть все

Восстановленная версия исходного сигнала на основе MODWT и уровня реконструкции, возвращенной как вектор - строка.

Ссылки

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

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

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

|

Введенный в R2015b

Для просмотра документации необходимо авторизоваться на сайте