wpspectrum

Пакетный спектр вейвлета

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

Синтаксис

[SPEC,TIMES,FREQ] = wpspectrum(WPT,Fs)
[...] = wpspectrum(WPT,Fs,'plot')
[...,TNFO] = wpspectrum(...)

Описание

[SPEC,TIMES,FREQ] = wpspectrum(WPT,Fs) возвращает матрицу пакетных оценок спектра вейвлета, SPEC, для бинарного пакетного объекта дерева вейвлета, WPT. Fs является частотой дискретизации в Герц. SPEC является 2J-by-N матрица, где J является уровнем пакета вейвлета, преобразовывают, и N является длиной временных рядов. TIMES является 1 N вектором времен, и FREQ является 1 2J вектором частот.

[...] = wpspectrum(WPT,Fs,'plot') отображает пакетный спектр вейвлета.

[...,TNFO] = wpspectrum(...) возвращает терминальные узлы пакетного дерева вейвлета в порядке частоты.

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

WPT

WPT является бинарным пакетным деревом вейвлета класса wptree.

Fs

Частота дискретизации в Герц как скаляр класса дважды.

Значение по умолчанию: 1

plot

Вектор символов 'plot' отображает пакетный спектр вейвлета. Введите 'plot' после Fs, чтобы произвести график пакетного спектра вейвлета.

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

SPEC

Пакетный спектр вейвлета. SPEC является матрицей 2J-by-N, где J является уровнем пакета вейвлета, преобразовывают, и N является длиной узла 0 в пакетном объекте дерева вейвлета.

Частотным интервалом между строками SPEC является Fs/2J+1.

TIMES

Временной вектор. TIMES является вектором времен в секундах, равных в длине к узлу 0 из пакетного объекта дерева вейвлета. Время, располагая с интервалами между элементами является 1/Fs.

FREQ

Вектор частоты. FREQ является вектором частот длины 2J, где J является уровнем пакетного объекта дерева вейвлета. Частотным интервалом в FREQ является Fs/2J+1.

TNFO

Терминальные узлы. TNFO является вектором терминальных узлов пакетного объекта дерева вейвлета в порядке частоты.

Примеры

свернуть все

Скрипт Open Live Script

Этот пример показывает пакетный спектр вейвлета для сигнала, состоящего из двух синусоид с непересекающейся поддержкой.

Задайте вейвлет.

fs = 500;
t = 0:1/fs:4;
y = sin(32*pi*t).*(t<2) + sin(128*pi*t).*(t>=2);
plot(t,y); 
axis tight
title('Analyzed Signal');

Задайте пакетный спектр вейвлета.

level = 6;
wpt = wpdec(y,level,'sym6');
figure;
[S,T,F] = wpspectrum(wpt,fs,'plot');

Скрипт Open Live Script

Создайте сигнал щебета.

fs = 1000;
t = 0:1/fs:2; 
% create chirp signal        
y = sin(256*pi*t.^2);

Постройте анализируемый сигнал.

plot(t,y); 
axis tight
title('Analyzed Signal');

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

level = 6;
wpt = wpdec(y,level,'sym8');
figure;
[S,T,F] = wpspectrum(wpt,fs,'plot');

Больше о

свернуть все

Пакетный спектр вейвлета

Пакетный спектр вейвлета содержит абсолютные значения коэффициентов от упорядоченных частоте терминальных узлов входного двоичного пакетного дерева вейвлета. Терминальные узлы обеспечивают, самый прекрасный уровень разрешения частоты в пакете вейвлета преобразовывают. Если J обозначает, что уровень пакета вейвлета преобразовывает, и Fs является частотой дискретизации, терминальные узлы аппроксимированные полосовые фильтры формы:

[nFs2J+1,(n+1)Fs2J+1) n=0,1,2,3,…2J−1

На терминальном уровне пакетного дерева вейвлета преобразование делит интервал от 0 до частоты Найквиста в полосы аппроксимированной ширины Fs/2J+1.

Алгоритмы

wpspectrum вычисляет пакетный спектр вейвлета можно следующим образом:

  • Извлеките пакетные коэффициенты вейвлета, соответствующие терминальным узлам. Примите абсолютное значение коэффициентов.

  • Закажите пакетные коэффициенты вейвлета упорядоченным расположением частоты.

  • Определите степень времени на исходной оси времени, соответствующей каждому пакетному коэффициенту вейвлета. Повторите каждый пакетный коэффициент вейвлета, чтобы заполнить разрывы времени между соседними пакетными коэффициентами вейвлета и создать вектор, равный в длине к узлу 0 из пакетного объекта дерева вейвлета.

Ссылки

Wickerhauser, лекции т-х по пакетным алгоритмам вейвлета, техническому отчету, Вашингтонскому университету, отделу математики, 1992.

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

|

Темы

Представленный в R2010b

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