В этом примере показано, как выполнить непрерывный анализ вейвлета сигнала острого выступа. Можно использовать cwt
для анализа с помощью аналитического вейвлета и wtmm
изолировать и охарактеризованная сингулярность.
Загрузите и постройте сигнал острого выступа. Отобразите его определение в командной строке.
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 сигнала острого выступа. CWT использует аналитический вейвлет Морзе с гаммой, равной 2 и параметр пропускной способности времени 2,5. Заметьте узкую область в scalogram, сходящемся к самой прекрасной шкале (самая высокая частота). Это указывает на разрыв в сигнале.
cwt(cuspamax,'WaveletParameters',[2 2.5]);
Получите график линий максимумов вейвлета с помощью вейвлета, преобразовывают максимумы модуля. wtmm
возвращает оценки экспонент Держателя, которые характеризуют изолированную сингулярность в сигнале. Заметьте, что острый выступ показывают очень ясно с помощью wtmm
.
wtmm(cuspamax,'ScalingExponent','local');