mdwtdec

Разложение вейвлета мультисигнала 1-D

Синтаксис

DEC = mdwtdec(DIRDEC,X,LEV,WNAME) DEC = mdwtdec(DIRDEC,X,LEV,LoD,HiD,LoR,HiR) DEC = mdwtdec(...,'mode',EXTMODE)

Описание

DEC = mdwtdec(DIRDEC,X,LEV,WNAME) возвращает разложение вейвлета на уровне LEV каждой строки (если DIRDEC = 'r') или каждый столбец (если DIRDEC = 'c') матричного X, с помощью вейвлета WNAME.

Вывод DEC является структурой со следующими полями:

`'dirDec'`	Индикатор Direction: `'r'` (строка) или `'c'` (столбец)
`'level'`	Уровень разложения DWT
`'wname'`	Имя вейвлета
`'dwtFilters'`	Структура с четырьмя полями `LoD`, `HiD`, `LoR` и `HiR`
`'dwtEXTM'`	Режим расширения DWT (см. `dwtmode`),
`'dwtShift'`	DWT переключают параметр (0 или 1)
`'dataSize'`	Размер `X`
`'ca'`	Коэффициенты приближения на уровне `LEV`
`'cd'`	Массив ячеек коэффициентов детали, от уровня 1 до уровня `LEV`

Коэффициенты cA и cD{k} (for k = 1 to LEV) являются матрицами и хранятся в строках если DIRDEC = 'r' или в столбцах если DIRDEC = 'c'.

DEC = mdwtdec(DIRDEC,X,LEV,LoD,HiD,LoR,HiR) использует четыре фильтра вместо имени вейвлета.

DEC = mdwtdec(...,'mode',EXTMODE) вычисляет разложение вейвлета с режимом расширения EXTMODE, который вы задаете (см. dwtmode для допустимых дополнительных режимов).

Примеры

свернуть все

Анализируйте мультисигналы

Скрипт Open Live Script

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

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

load thinker

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

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

dec = struct with fields:
        dirDec: 'r'
         level: 2
         wname: 'db2'
    dwtFilters: [1x1 struct]
       dwtEXTM: 'sym'
      dwtShift: 0
      dataSize: [192 96]
            ca: [192x26 double]
            cd: {[192x49 double]  [192x26 double]}

Вычислите фильтры, сопоставленные с вейвлетом db2.

[LoD,HiD,LoR,HiR] = wfilters('db2');

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

decBIS = mdwtdec('r',X,2,LoD,HiD,LoR,HiR)

decBIS = struct with fields:
        dirDec: 'r'
         level: 2
         wname: ''
    dwtFilters: [1x1 struct]
       dwtEXTM: 'sym'
      dwtShift: 0
      dataSize: [192 96]
            ca: [192x26 double]
            cd: {[192x49 double]  [192x26 double]}

Подтвердите приближение и детализируйте коэффициенты обоих разложений, идентичны.

max(abs(dec.ca(:)-decBIS.ca(:)))

ans = 0

max(abs(dec.cd{1}(:)-decBIS.cd{1}(:)))

ans = 0

max(abs(dec.cd{2}(:)-decBIS.cd{2}(:)))

ans = 0

Ссылки

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

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

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

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

mdwtrec | wavedec

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