Непрерывный вейвлет-анализ модулированных сигналов

В этом примере показано, как использовать непрерывное вейвлет-преобразование (CWT) для анализа модулированных сигналов.

Загрузите сигнал квадратичной частоты. Частота сигнала начинается приблизительно с 500 Гц при t = 0, уменьшается до 100 Гц при t = 2 и увеличивается обратно до 500 Гц при t = 4. Частота дискретизации составляет 1 кГц.

load quadchirp;
fs = 1000;

Получение частотно-временного графика этого сигнала с использованием CWT с импульсным импульсом. Импульсный импульс является хорошим выбором для CWT, когда ваши сигналы колеблются и вы больше заинтересованы во временном частотном анализе, чем в локализации переходных процессов.

[cfs,f] = cwt(quadchirp,'bump',fs);
helperCWTTimeFreqPlot(cfs,tquad,f,'surf','CWT of Quadratic Chirp','Seconds','Hz')

CWT ясно показывает эволюцию частоты квадратичного чирпа во времени. Квадратичная частота является частотно-модулированным сигналом. Хотя этот сигнал является синтетическим, частотная и амплитудная модуляция также часто происходит и в собственных сигналах. Используйте CWT для получения частотно-временного анализа эхолокационного импульса, испускаемого большой коричневой битой (Eptesicus Fuscus). Интервал выборки составляет 7 микросекунд. Используйте импульсный импульс с 32 голосами на октаву. Спасибо Кертису Кондону, Кену Уайту и Элу Фенгу из Центра Бекмана при Иллинойсском университете за данные о летучей мыши и разрешение использовать её в этом примере.

load batsignal
t = 0:DT:(numel(batsignal)*DT)-DT;
[cfs,f] = cwt(batsignal,'bump',1/DT,'VoicesPerOctave',32);
helperCWTTimeFreqPlot(cfs,t.*1e6,f./1e3,'surf','Bat Echolocation (CWT)',...
    'Microseconds','kHz')

В качестве окончательного примера можно получить частотно-временной анализ некоторых сейсмографических данных, зарегистрированных во время землетрясения в Кобе в 1995 году. Данные являются сейсмографическими (вертикальное ускорение, нм/кв.сек) измерениями, зарегистрированными в Университете Тасмании, Хобарт, Австралия, 16 января 1995 начиная с 20:56:51 (GMT) и продолжающимися в течение 51 минут 1 секундными интервалами. Используйте аналитический вейвлет Морса по умолчанию.

load kobe;
dt = 1;
cwt(kobe,1);
title('CWT of 1995 Kobe Earthquake Seismograph Data');