Используйте анализатор сигналов для отображения гиперболической чирпы и ее скалограммы. В этом примере требуется лицензия Wavelet Toolbox™.
Генерация сигнала, дискретизированного на частоте 1 кГц в течение 1 секунды. Сигнал состоит из двух гиперболических чирпов. Один из чирпов активен между 0,1 и 0,68 секундами. Другой чирп имеет удвоенную амплитуду и одну треть мгновенной частоты первого чирпа. Второй чирп активен в диапазоне от 0,1 секунды до 0,75 секунды.
n = 1000; t = (0:n-1)/n; ff = [1 2]*(sin(pi*[15 5]'./(t-0.8001)).*(t>0.1 & t<[0.68 0.75]'));
Чирп модулируется огибающей, которая начинается с 0,1 секунды и является ненулевой в течение 0,7 секунды.
envl = zeros(size(t)); M = 0.7*n; envl(n/10+(1:M)) = tukeywin(M,0.8); nw = ff.*envl;
Отображение сигнала в анализаторе сигналов. Нажмите кнопку Time-Frequency
▼ на
вкладке Display и выберите Scalogram.
signalAnalyzer(nw,'TimeValues',t)
На вкладке Скалограмма установите максимальное значение параметра Голоса на октаву. Установите нижний предел амплитуды равным 0,5.
pspectrum | cwt (инструментарий вейвлета) | cwtfilterbank (инструментарий вейвлета)