ilwt

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

    Описание

    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-канальные данные EEG 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
    

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

    свернуть все

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

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

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

    Коэффициенты детализации, заданные как массив ячеек L -на-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' аргумент имя-значение, только если все элементы входа являются целыми числами.

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

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

    свернуть все

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

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

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

    расширить все

    Поведение изменено в R2021a

    См. также

    | | | |

    Введенный в R2021a