Exponenta Banner
exponenta event banner

lwtcoef

Извлеките или восстановите 1-D LWT коэффициенты небольшой волны и ортогональные прогнозы

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

    Синтаксис

    y = lwtcoef (ca, cd)
    y = lwtcoef (ca, cd, имя, значение)

    Описание

    y = lwtcoef(ca,cd) возвращает коэффициенты аппроксимации уровня 1, соответствующие коэффициентам аппроксимации и детализации, ca и cdсоответственно. ca и cd являются выходами lwt.

    пример

    y = lwtcoef(ca,cd,Name,Value) указывает параметры, использующие один или несколько аргументов значения имени. Например, y = lwtcoef(ca,cd,'OutputType','coefficients') задает выходные коэффициенты.

    Примеры

    свернуть все

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

    Загрузите сигнал 1-D длиной 2048. Постройте график сигнала.

    load wecg
plot(wecg)
title('Signal')
ylabel('Amplitude')
axis tight

    Figure contains an axes. The axes with title Signal contains an object of type line.

    Создание схемы подъема, связанной с db4 вейвлет. Используйте схему подъема для получения вейвлет-разложения сигнала до максимального уровня. Подтвердите, что длина массива ячеек коэффициентов детализации равна floor(log2(N)), где N - длина сигнала. 

    wv = 'db4';
lsc = liftingScheme('Wavelet',wv);
[ca,cd] = lwt(wecg,'LiftingScheme',lsc);
[length(cd) floor(log2(length(wecg)))]
    ans = 1×2

    11    11

    Извлеките и постройте график коэффициентов аппроксимации на уровне 3. Подтвердите, что длина извлечения равна одной восьмой длины исходного сигнала.

    approxCf = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','coefficients','Level',3);
[2048/(2^3) length(approxCf)]
    ans = 1×2

   256   256


    plot(approxCf)
title('Level 3 Approximation Coefficients')
ylabel('Amplitude')
axis tight

    Figure contains an axes. The axes with title Level 3 Approximation Coefficients contains an object of type line.

    Получение ортогональной проекции коэффициентов аппроксимации уровня 3. Также получают ортогональные проекции коэффициентов детализации на уровнях 1, 2 и 3. Постройте график результатов.

    approx3 = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','projection','Level',3);
det3 = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','projection','Level',3,'Type','detail');
det2 = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','projection','Level',2,'Type','detail');
det1 = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','projection','Level',1,'Type','detail');
subplot(4,1,1)
plot(approx3)
title('Projection - Approximation')
axis tight
subplot(4,1,2)
plot(det3)
title('Projection - Level 3 Details')
axis tight
subplot(4,1,3)
plot(det2)
title('Projection - Level 2 Details')
axis tight
subplot(4,1,4)
plot(det1)
title('Projection - Level 1 Details')
axis tight

    Figure contains 4 axes. Axes 1 with title Projection - Approximation contains an object of type line. Axes 2 with title Projection - Level 3 Details contains an object of type line. Axes 3 with title Projection - Level 2 Details contains an object of type line. Axes 4 with title Projection - Level 1 Details contains an object of type line.

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

    max(abs(wecg-(approx3+det3+det2+det1)))
    ans = 1.3323e-15

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

    свернуть все

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

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

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

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

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

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

    Пример: y = lwtcoef(ca,cd,'LiftingScheme',lsc,'OutputType','coefficients','Level',2) использует схему подъема lsc для извлечения коэффициентов аппроксимации на уровне 2.

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

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

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

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

    Тип вывода, указанный как один из:

    • 'coefficients' - Извлечение коэффициентов аппроксимации или детализации

    • 'projection' - Возврат проекции (реконструкции) коэффициентов приближения или детализации

    Пример: y = lwtcoef(ca,cd,'OutputType','projection','Type','detail') возвращает проекцию, соответствующую коэффициентам детализации в самой точной шкале.

    Тип коэффициентов для извлечения или восстановления, указанный как 'approximation' или 'detail'.

    Пример: y = lwtcoef(ca,cd,'Type','detail') извлекает коэффициенты детализации в лучшем масштабе.

    Уровень коэффициентов для извлечения или восстановления, указанный как целое число в диапазоне [1,N], где N - длина cd.

    Пример: y = lwtcoef(ca,cd,'LiftingScheme',lsc,'Level',3) использует схему подъема lsc для извлечения коэффициентов аппроксимации на уровне 3.

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

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

    свернуть все

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

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

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

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

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

    См. также

    | |

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

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

    Поддержка