centfrq
Частота центра вейвлета
Синтаксис
FREQ = centfrq('wname')
FREQ = centfrq('wname',ITER)
[FREQ,XVAL,RECFREQ] = centfrq('wname',ITER,'plot')
Описание
FREQ = centfrq(' возвращает центральную частоту в герц функции вейвлета, wname')'wname' (см. wavefun для получения дополнительной информации).
Для FREQ = centfrq(', wname',ITER)ITER является количеством итераций, выполняемых функциональным wavefun, который используется, чтобы вычислить вейвлет.
[FREQ,XVAL,RECFREQ] = centfrq(' возвращает, кроме того, связанное центральное основанное на частоте приближение wname',ITER,'plot')RECFREQ на 2ITER сетка точек XVAL и строит функцию вейвлета и RECFREQ.
Примеры
Определите центральную частоту
Этот пример показывает, как определить центральную частоту в герц для Добечиса меньше всего - асимметричный вейвлет с 4 исчезающими моментами.
cfreq = centfrq('sym4');Получите вейвлет и создайте синусоиду с частотой, равной центральной частоте, cfreq, вейвлета. Используйте стартовую фазу для синусоиды, чтобы визуализировать, как колебание в синусоиде совпадает с колебанием в вейвлете.
[~,psi,xval] = wavefun('sym4'); y = cos(2*pi*cfreq*xval-pi); plot(xval,psi,'linewidth',2); hold on; plot(xval,y,'r');
Преобразуйте шкалы в частоты
Этот пример показывает, чтобы преобразовать шкалы в частоты для вейвлета Morlet. Существует аппроксимированная обратная связь между шкалой и частотой. А именно, шкала обратно пропорциональна частоте с коэффициентом пропорциональности, являющимся центральной частотой вейвлета.
Создайте вектор шкал с 32 речью на октаву более чем 5 октав для данных, выбранных на уровне 1 кГц.
Fs = 1000; numvoices = 32; a0 = 2^(1/numvoices); numoctaves = 5; scales = a0.^(numvoices:numvoices*numoctaves).*1/Fs;
Преобразуйте шкалы, чтобы аппроксимировать частоты в герц для вейвлета Morlet.
Frq = centfrq('morl')./scales;Можно также использовать scal2frq, чтобы преобразовать шкалы, чтобы аппроксимировать частоты в герц.