Непрерывный вейвлет сигнала Cusp

Этот пример показывает, как выполнить непрерывный вейвлет cusp-сигнала. Можно использовать cwt для анализа с использованием аналитического вейвлета и wtmm для выделения и характеристики особенностей.

Загрузите и постройте график сигнала cusp. Отобразите его определение в командной строке.

load cuspamax; 
plot(cuspamax); grid on;

Figure contains an axes. The axes contains an object of type line.

disp(caption)
x = linspace(0,1,1024); y = exp(-128*((x-0.3).^2))-3*(abs(x-0.7).^0.4);

Получите и просмотрите CWT сигнала cusp. CWT использует аналитический вейвлет Морса с гаммой, равной 2, и параметром временной полосы 2,5. Заметьте узкую область в скалограмме, сходящуюся к самой мелкой шкале (самая высокая частота). Это указывает на разрыв сигнала.

cwt(cuspamax,'WaveletParameters',[2 2.5]);

Figure contains an axes. The axes with title Magnitude Scalogram contains 3 objects of type image, line, area.

Получите график вейвлета линий maxima с помощью вейвлета transform modulus maxima. wtmm возвращает оценки экспонентов Holder, которые характеризуют изолированные особенности в сигнале. Заметьте, что cusp показан очень четко с помощью wtmm.

wtmm(cuspamax,'ScalingExponent','local');

Figure contains an axes and an object of type uitable. The axes with title Wavelet Transform Maxima Lines contains 7 objects of type image, line.