mdwtrec

Реконструкция вейвлета мультисигнала 1-D

Синтаксис

X = mdwtrec(DEC) X = mdwtrec(DEC,IDXSIG) Y = mdwtrec(DEC,TYPE,LEV) A = mdwtrec(DEC,'a') A = mdwtrec(DEC,'a',LEVDEC) D = mdwtrec(DEC,'d') CA = mdwtrec(DEC,'ca') CA = mdwtrec(DEC,'ca',LEVDEC) CD = mdwtrec(DEC,'cd',MODE) CFS = mdwtrec(DEC,'cfs',MODE) Y = mdwtrec(...,IDXSIG)

Описание

X = mdwtrec(DEC) возвращает исходную матрицу сигналов, начинающих с DEC структуры разложения вейвлета (см. mdwtdec).

X = mdwtrec(DEC,IDXSIG) восстанавливает сигналы, индексы которых даны векторным IDXSIG.

Y = mdwtrec(DEC,TYPE,LEV) извлечения или восстанавливают деталь или коэффициенты приближения на уровне LEV в зависимости от значения TYPE. Максимальным значением для LEV является LEVDEC = DEC.level.

То, когда ВВОДЯТ, равно:

'cd' или 'ca', коэффициенты ЛЕВА уровня извлечены.
'd' или 'a', коэффициенты ЛЕВА уровня восстановлены.
'a' или 'ca', ЛЕВ должен быть таков что 0 ≤ LEV ≤ LEVDEC.
'd' или 'cd', ЛЕВ должен быть таков что 1 ≤ LEV ≤ LEVDEC.

A = mdwtrec(DEC,'a') эквивалентно A = mdwtrec(DEC,'a',LEVDEC).

D = mdwtrec(DEC,'d') возвращает матрицу, содержащую сумму всех подробностей, так, чтобы X = A + D.

CA = mdwtrec(DEC,'ca') эквивалентно CA = mdwtrec(DEC,'ca',LEVDEC).

CD = mdwtrec(DEC,'cd',MODE) возвращает матрицу, содержащую коэффициенты всех подробностей.

CFS = mdwtrec(DEC,'cfs',MODE) возвращает матрицу, содержащую все коэффициенты.

Для MODE = 'descend' коэффициенты конкатенированы от уровня, LEVDEC к уровню 1 и MODE = 'descend' конкатенирует от уровня 1 до уровня LEVDEC). Значением по умолчанию является MODE = 'descend'. Конкатенация сделана построчной если DEC.dirDEC = 'r' или по столбцам если DEC.dirDEC = 'c'.

Y = mdwtrec(...,IDXSIG) извлечения или восстанавливают деталь или коэффициенты приближения для сигналов, индексы которых даны векторным IDXSIG.

Примеры

свернуть все

Восстановите мультисигналы

Скрипт Open Live Script

Этот пример показывает, как восстановить и конкретный сигнал мультисигнала в мультисигнале.

Загрузите 1D мультисигнал.

load thinker

Выполните разложение на уровне 2 с помощью вейвлета db2.

dec = mdwtdec('r',X,2,'db2');

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

XR = mdwtrec(dec);

Вычислите ошибку реконструкции.

errREC = max(abs(X(:)-XR(:)))

errREC = 2.1026e-10

Восстановите исходный сигнал в индексе 31, соответствующее приближение на уровне 2 и детали на уровнях 1 и 2.

idx = 31;
Y = mdwtrec(dec,idx);
A2 = mdwtrec(dec,'a',2,idx);
D2 = mdwtrec(dec,'d',2,idx);
D1 = mdwtrec(dec,'d',1,idx);

Вычислите ошибку реконструкции для сигнала 31.

errREC = max(abs(Y-A2-D2-D1))

errREC = 6.7946e-14

Ссылки

Daubechies, я., Десять лекций по вейвлетам, ряду конференции CBMS-NSF в прикладной математике. SIAM Эд., 1992.

Mallat, S., “Теория для мультиразрешения сигнализирует о разложении: представление вейвлета”, Анальный Шаблон IEEE. и Машина Intell., издание 11, № 7, 1989, стр 674–693.

Мейер, Y., Ondelettes и opérateurs, Том 1, Герман Эд. (Английский перевод: Вейвлеты и операторы, Кембриджское Нажатие Унив. 1993.)

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

mdwtdec | waverec

Документация