idwt3

Одноуровневый обратный дискретный 3-D вейвлет преобразовывает

Синтаксис

X = idwt3(WT) C = idwt3(WT,TYPE)

Описание

idwt3 команда выполняет одноуровневую 3D реконструкцию вейвлета, начинающую с одноуровневого 3D разложения вейвлета.

X = idwt3(WT) вычисляет одноуровневый восстановленный трехмерный массив X, на основе 3D разложения вейвлета, сохраненного в структуре WT. Эта структура содержит следующие поля.

`sizeINI`	Размер 3D массива `X`.
`mode`	Имя вейвлета преобразовывает дополнительный режим.
`filters`	Структура с 4 полями, `LoD`, `HiD`, `LoR`, `HiR`, которые содержат фильтры, используемые для DWT.
`dec`	2 x 2 x 2 массива ячеек, содержащие коэффициенты разложения. `dec{i,j,k}iJ`, `k = 1` или `2` содержит коэффициенты, полученные фильтрацией lowpass (для `i` или `j` или `k = 1`) или фильтрация высокой передачи (для `i` или `j` или `k = 2`).

C = idwt3(WT,TYPE) вычисляет одноуровневое, восстановленное компонентно-ориентированный на 3D разложении вейвлета. Допустимые значения для TYPE :

Группа из трех символов 'xyz', один на направление, с 'x'Y и 'z' выбранный в наборе {'a'DLH} или в соответствующем прописном наборе {'A'DLH}), где 'A' (или 'L') задает фильтр lowpass и 'D' (или 'H') задает фильтр highpass.
Символьный 'd' (или 'h' или 'D' или 'H') который задает сумму всех компонентов, отличающихся от компонента lowpass.

Примеры

свернуть все

Одноуровневая 3D реконструкция вейвлета

Скрипт Open Live Script

Задайте исходные 3-D данные.

X  = reshape(1:64,4,4,4)

X = 
X(:,:,1) =

     1     5     9    13
     2     6    10    14
     3     7    11    15
     4     8    12    16


X(:,:,2) =

    17    21    25    29
    18    22    26    30
    19    23    27    31
    20    24    28    32


X(:,:,3) =

    33    37    41    45
    34    38    42    46
    35    39    43    47
    36    40    44    48


X(:,:,4) =

    49    53    57    61
    50    54    58    62
    51    55    59    63
    52    56    60    64

Анализируйте X использование 'db1'.

wt = dwt3(X,'db1');

Восстановите X от коэффициентов. Проверьте, что восстановленные данные соглашаются с исходными данными обработать точность машинным способом.

XR = idwt3(wt);

dff = max(abs(X-XR))

dff = 
dff(:,:,1) =

   1.0e-13 *

    0.0266    0.0355    0.0888    0.1066


dff(:,:,2) =

   1.0e-13 *

    0.1066    0.1066    0.2132    0.2132


dff(:,:,3) =

   1.0e-13 *

    0.1421    0.1421    0.2132    0.2132


dff(:,:,4) =

   1.0e-13 *

    0.3553    0.3553    0.2842    0.2842

Вычислите восстановленное приближение, которое состоит из компонента lowpass.

A  = idwt3(wt,'aaa');

Вычислите сумму всех компонентов, отличающихся от компонента lowpass.

D  = idwt3(wt,'d');

Восстановите компонент, сопоставленный с lowpass в x и z направлениях и highpass в направлении Y.

ADA  = idwt3(wt,'ada');

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

dwt3 | wavedec3 | waverec3

Введен в R2010a

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