В этом примере показано, как определить центральную частоту в герцах для наименее асимметричного вейвлета Daubechies с 4 моментами исчезновения.

cfreq = centfrq('sym4');

Получить вейвлет и создать синусоидальную волну с частотой, равной центральной частоте, cfreq, вейвлета. Используйте начальную фазу $$-\delta$$ для синусоидальной волны, чтобы визуализировать, как колебания синусоидальной волны соответствуют колебаниям вейвлета.

[~,psi,xval] = wavefun('sym4');
y = cos(2*pi*cfreq*xval-pi);
plot(xval,psi,'linewidth',2); 
hold on;
plot(xval,y,'r');

В этом примере показано преобразование масштабов в частоты для вейвлета Морле. Существует приблизительная обратная зависимость между масштабом и частотой. В частности, масштаб обратно пропорционален частоте, а константа пропорциональности является центральной частотой импульса.

Создайте вектор шкал с 32 голосами на октаву в течение 5 октав для данных, отобранных при частоте 1 кГц.

Fs = 1000;
numvoices = 32;
a0 = 2^(1/numvoices);
numoctaves = 5; 
scales = a0.^(0:numvoices*numoctaves-1).*1/Fs;

Преобразуйте шкалы в приблизительные частоты в герцах для импульса Морле.

Frq = centfrq('morl')./scales;

Также можно использовать scal2frq для преобразования масштабов в приблизительные частоты в герцах.