Документация

Анализ частоты времени модулируемых сигналов

Этот пример показывает, как использовать непрерывный вейвлет преобразовывает (CWT), чтобы анализировать сигналы совместно вовремя и частоту.

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

load quadchirp;
fs = 1000;
[S,F,T] = spectrogram(quadchirp,100,98,128,fs);
helperCWTTimeFreqPlot(S,T,F,'surf','STFT of Quadratic Chirp','Seconds','Hz')

Получите частотный график времени этого сигнала с помощью 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 речью на октаву. Благодаря Кертису Кондону, Кену Вайту и Аль Фэну из Центра Бекмана в Университете Иллинойса для 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')

Получите и постройте STFT bat данных.

[S,F,T] = spectrogram(batsignal,50,48,128,1/DT);
helperCWTTimeFreqPlot(S,T.*1e6,F./1e3,'surf','Bat Echolocation (STFT)',...
    'Microseconds','kHz')

И для моделируемых и для естественных модулируемых сигналов, CWT обеспечивает результаты, подобные STFT.

Для итогового примера получите анализ частоты времени некоторых данных о сейсмографе, зарегистрированных во время 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');