Этот пример показывает, как выполнить непрерывный вейвлет cusp-сигнала. Можно использовать cwt
для анализа с использованием аналитического вейвлета и wtmm
для выделения и характеристики особенностей.
Загрузите и постройте график сигнала cusp. Отобразите его определение в командной строке.
load cuspamax; plot(cuspamax); grid on;
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]);
Получите график вейвлета линий maxima с помощью вейвлета transform modulus maxima. wtmm
возвращает оценки экспонентов Holder, которые характеризуют изолированные особенности в сигнале. Заметьте, что cusp показан очень четко с помощью wtmm
.
wtmm(cuspamax,'ScalingExponent','local');