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

В этом примере показано, как использовать непрерывный вейвлет преобразовывает (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')

Figure contains an axes. The axes with title CWT of Quadratic Chirp contains an object of type surface.

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

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')

Figure contains an axes. The axes with title Bat Echolocation (CWT) contains an object of type surface.

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

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

Figure contains an axes. The axes with title CWT of 1995 Kobe Earthquake Seismograph Data contains 3 objects of type image, line, area.

Для просмотра документации необходимо авторизоваться на сайте