tqwtmra

Настраиваемый анализ мультиразрешения Q-фактора

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

    Синтаксис

    mra = tqwtmra(wt,n)
    mra = tqwtmra(wt,n,QualityFactor=qf)
    tqwtmra(___)

    Описание

    mra = tqwtmra(wt,n) возвращает настраиваемый анализ мультиразрешения вейвлета Q-фактора (MRA) для анализа TQWT, wt, полученный с добротностью по умолчанию 1.

    пример

    mra = tqwtmra(wt,n,QualityFactor=qf) использует добротность qf в получении настраиваемого Q-factor MRA. qf должен совпадать со значением, используемым в получении wt от tqwt.

    tqwtmra(___) без выходных аргументов строит настраиваемый вейвлет Q-фактора MRA на новом рисунке. Для данных с комплексным знаком действительная часть построена в первом, раскрашивают MATLAB® матрица последовательности цветов и мнимая часть построены во втором цвете. Этот синтаксис не поддерживает многомерный MRAs.

    Примеры

    свернуть все

    Скрипт Open Live Script

    Загрузите сигнал ECG. Получите TQWT сигнала вниз к уровню 6 с добротностью 2.

    load wecg
wt = tqwt(wecg,QualityFactor=2,Level=6);

    Получите настраиваемый Q-factor MRA сигнала.

    mra = tqwtmra(wt,length(wecg),QualityFactor=2);

    Постройте исходный сигнал и поддиапазон lowpass.

    plot(wecg)
hold on
plot(mra(end,:),linewidth=2)
hold off
axis tight
legend(["Original","Lowpass"])

    Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original, Lowpass.

    Подтвердите, что сумма вдоль строк MRA равняется исходному сигналу.

    mraSum = sum(mra,1);
max(abs(mraSum(:)-wecg(:)))
    ans = 9.4369e-16
    Скрипт Open Live Script

    Загрузите данные о землетрясении Кобе. Получите настраиваемое преобразование вейвлета Q-фактора данных с помощью добротности 3. 

    load kobe
qf = 3;
wt = tqwt(kobe,QualityFactor=qf);

    Идентифицируйте поддиапазоны, которые содержат по крайней мере 15% полной энергии. Обратите внимание на то, что последний элемент wt содержит коэффициенты поддиапазона lowpass.

    EnergyBySubband = cellfun(@(x)norm(x,2)^2,wt)./norm(kobe,2)^2*100;
idx15 = EnergyBySubband >= 15;
bar(EnergyBySubband)
title("Percent Energy By Subband")
xlabel("Subband")
ylabel("Percent Energy")

    Figure contains an axes object. The axes object with title Percent Energy By Subband contains an object of type bar.

    Получите анализ мультиразрешения и суммируйте те компоненты MRA, соответствующие ранее идентифицированным поддиапазонам.

    mra = tqwtmra(wt,numel(kobe),QualityFactor=qf);
ts = sum(mra(idx15,:));
plot([kobe ts'])
axis tight
legend("Original Data","Large Energy Components",...
        Location="NorthWest")
xlabel("Time (s)")

    Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original Data, Large Energy Components.

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

    свернуть все

    Настраиваемый вейвлет Q-фактора преобразовывает в виде массива ячеек. Элементы wt содержите поддиапазон вейвлета и коэффициенты lowpass. wt как ожидают, будет выходом tqwt.

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

    Исходная длина сигнала в выборках в виде положительного целого числа. Если исходная длина сигнала n является нечетным, n расширен к n+1 получить MRA и итоговую выборку удалено прежде, чем возвратить MRA.

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

    Добротность в виде скаляра с действительным знаком, больше, чем или равный 1. Если незаданный, значения по умолчанию добротности к 1.

    Пример: wt = tqwtmra(qt,2024,QualityFactor=1.5) задает добротность 1,5.

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

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

    свернуть все

    Анализ мультиразрешения, возвращенный как массив. mra Ns-by-N-by-C-by-B массив, где Ns обозначает, что количество поддиапазонов в настраиваемом вейвлете Q-фактора преобразовывает упорядоченный путем уменьшения центральной частоты, N является количеством выборок сигнала вовремя, C является количеством каналов, и B является пакетным размером.

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

    Ссылки

    [1] Selesnick, Иван В. “Преобразование вейвлета С Настраиваемым Q-фактором”. Транзакции IEEE на Обработке сигналов 59, № 8 (август 2011): 3560–75. https://doi.org/10.1109/TSP.2011.2143711.

    Расширенные возможности

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

    |

    Темы

    Введенный в R2021b

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

    Поддержка

    Для просмотра документации необходимо авторизоваться на сайте