Exponenta Banner
exponenta event banner

ilwt

Обратное 1-D подъемное вейвлет-преобразование

свернуть все на странице

    Синтаксис

    xr = ilwt (ca, cd)
    xr = ilwt (ca, cd, имя, значение)

    Описание

    xr = ilwt(ca,cd) возвращает 1-D обратное вейвлет-преобразование на основе коэффициентов аппроксимации, caи массив ячеек коэффициентов детализации, cd. По умолчанию ilwt предполагает, что вы использовали схему подъема, связанную с db1 вейвлет для получения ca и cd. Если коэффициенты не изменяются, xr является идеальной реконструкцией сигнала.

    пример

    xr = ilwt(ca,cd,Name,Value) указывает параметры, использующие один или несколько аргументов значения имени. Например, xr = ilwt(ca,cd,'Wavelet','db2') задает ортогональный вейвлет db2.

    Для идеальной реконструкции все аргументы «имя-значение» должны совпадать с аргументами, используемыми в lwt получить ca и cd.

    Примеры

    свернуть все

    Открыть сценарий в реальном времени

    Создание схемы подъема, связанной с db3 вейвлет. Укажите целочисленный сигнал, длина которого равна степени 2.

    lsc = liftingScheme('Wavelet','db3');
n = 8;
sig = 1:2^n;

    Используйте схему подъема для получения целочисленного LWT сигнала до максимального уровня разложения.

    [ca,cd] = lwt(sig,'LiftingScheme',lsc,'Int2Int',true);

    Подтверждение коэффициентов детализации cd - массив ячеек, длина которого равна степени 2.

    length(cd)
    ans = 8

    Получить обратный LWT до уровня 0. Подтвердите идеальную реконструкцию.

    xrec0 = ilwt(ca,cd,'LiftingScheme',lsc,'Int2Int',true,'Level',0);
max(abs(xrec0(:)-sig(:)))
    ans = 0

    Получить обратный LWT до уровня 1.

    xrec1 = ilwt(ca,cd,'LiftingScheme',lsc,'Int2Int',true,'Level',1);

    Получите разложение сигнала на уровне 1. Подтвердить, что коэффициенты аппроксимации равны xrec1.

    [ca,cd] = lwt(sig,'LiftingScheme',lsc,'Int2Int',true,'Level',1);
max(abs(ca(:)-xrec1(:)))
    ans = 0
    Открыть сценарий в реальном времени

    Загрузка 23-канальных данных ЭЭГ Espiga3. Каналы расположены столбчато.

    load Espiga3
size(Espiga3)
    ans = 1×2

   995    23

    Получение LWT многоканального сигнала с помощью db4 вейвлет вниз до максимального уровня разложения по умолчанию.

    wv = 'db4';
[ca,cd] = lwt(Espiga3,'Wavelet',wv);

    Восстановите многоканальный сигнал.

    xrec = ilwt(ca,cd,'Wavelet',wv);

    Поскольку исходный сигнал имеет нечетное количество выборок в каждом канале, подтвердите, что реконструкция имеет на одну строку больше, чем исходный сигнал.

    size(xrec)
    ans = 1×2

   996    23

    Подтвердите, что последняя строка в реконструкции равна предыдущей строке.

    max(abs(xrec(end-1,:)-xrec(end,:)))
    ans = 5.6843e-14

    Удаление последней строки из реконструкции. Подтвердите, что результат равен исходному сигналу.

    xrec(end,:) = [];
max(abs(Espiga3(:)-xrec(:)))
    ans = 4.5475e-13

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

    свернуть все

    Коэффициенты аппроксимации (нижних частот) на самом крупном уровне, задаваемые как скаляр, вектор или матрица. Коэффициенты являются выходом lwt.

    Если ca и элементы cd - матрицы, xr является матрицей, где каждый столбец является обратным вейвлет-преобразованием соответствующих столбцов в ca и cd.

    Типы данных: single | double
    Поддержка комплексного номера: Да

    Коэффициенты детализации, заданные как массив L-by-1 ячеек, где L - уровень преобразования. Элементы cd в порядке уменьшения разрешающей способности. Коэффициенты являются выходом lwt.

    Если ca и элементы cd - матрицы, xr является матрицей, где каждый столбец является обратным вейвлет-преобразованием соответствующих столбцов в ca и cd.

    Типы данных: single | double
    Поддержка комплексного номера: Да

    Аргументы пары «имя-значение»

    Укажите дополнительные пары, разделенные запятыми Name,Value аргументы. Name является именем аргумента и Value - соответствующее значение. Name должен отображаться внутри кавычек. Можно указать несколько аргументов пары имен и значений в любом порядке как Name1,Value1,...,NameN,ValueN.

    Пример: xr = ilwt(ca,cd,'LiftingScheme',lsc,'Level',1) использует lsc схема подъема для выполнения обратного вейвлет-преобразования до уровня 1.

    Ортогональный или биортогональный вейвлет для использования в инверсном LWT, заданном как вектор символов или строковый скаляр. См. свойство Wavelet для liftingScheme для списка поддерживаемых вейвлетов. Для идеальной реконструкции указанный вейвлет должен быть тем же вейвлетом, который использовался для получения коэффициентов. ca и cd.

    Невозможно указать 'Wavelet' и 'LiftingScheme' аргументы «имя-значение» одновременно.

    Пример: xr = ilwt(ca,cd,'Wavelet','bior3.5') использует bior3.5 биорогональный вейвлет.

    Типы данных: char | string

    Схема подъема для использования в инверсном LWT, указанная как liftingScheme объект. Для идеальной реконструкции указанная схема подъема должна быть той же схемой подъема, которая использовалась для получения коэффициентов. ca и cd.

    Невозможно указать 'Wavelet' и 'LiftingScheme' аргументы «имя-значение» одновременно.

    Пример: xr = ilwt(ca,cd,'LiftingScheme',lScheme) использует lScheme схема подъема.

    Уровень реконструкции, заданный как неотрицательное целое число, меньшее или равное length(cd)-1. Если не указан, уровень реконструкции по умолчанию равен 0 и xr является идеальной реконструкцией сигнала.

    Пример: xr = ilwt(ca,cd,'Level',1) восстанавливает сигнал до уровня 1.

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

    Режим расширения для использования в LWT, указанный как 'periodic' (по умолчанию), 'zeropad', или 'symmetric'. Значение 'Extension' определяет способ удлинения сигнала на границах.

    Пример: xr = ilwt(ca,cd,'Extension','symmetric') определяет симметричный режим расширения.

    Обработка целочисленных данных, указанная как числовая или логическая 1 (true) или 0 (false).

    • 1 (true) - Сохранение целочисленных данных

    • 0 (false) - Не сохранять целочисленные данные

    Укажите 'Int2Int' аргумент name-value, только если все элементы ввода являются целыми числами.

    Пример: xr = ilwt(ca,cd,'Int2Int',true) сохраняет целочисленные данные.

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

    свернуть все

    Обратное вейвлет-преобразование ca и cd, возвращается в виде вектора или матрицы. Если ca является скаляром или вектором, и элементы cd являются векторами, xr является вектором. Если ca и элементы cd - матрицы, xr является матрицей, где каждый столбец является обратным вейвлет-преобразованием соответствующих столбцов в ca и cd.

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

    Вопросы совместимости

    развернуть все

    В R2021a изменилось поведение

    См. также

    | | | |

    Представлен в R2021a

    Документация инструментария вейвлет

    Поддержка