Этот пример показывает, как выполнить частотно-временной анализ с использованием непрерывного вейвлет (CWT). Непрерывный вейвлет обеспечивает временной/частотно-временной анализ сигналов и изображений. Программное обеспечение Wavelet Toolbox™ имеет как командную строку, так и интерактивную функциональность для поддержки непрерывного вейвлет 1-D сигналов.
Создайте сигнал, состоящий из двух синусоидов с частотами 100 и 50 Гц и белым шумом. Поддержка двух синусоидов несвязна. 100-Hz синусоида начинается с t = 0 и имеет длительность 1 секунду. Синусоида 100-Hz имеет амплитуду 2. Синусоида 50-Hz начинается с трех секунд и имеет длительность в две секунды. Синусоида 50-Hz имеет амплитуду 1. Частота дискретизации составляет 1 кГц. Длина сигнала составляет 5000 выборки.
Fs = 1000; t = linspace(0,5,5e3); x = 2*cos(2*pi*100*t).*(t<1)+cos(2*pi*50*t).*(3<t)+0.3*randn(size(t));
Постройте график сигнала.
plot(t,x) xlabel('Time (s)') ylabel('Amplitude')
Использование cwt для получения CWT сигнала и построения его скалограммы. Величины синусоидных компонентов в шкале палитры по существу являются их амплитудами, хотя они находятся в разных шкалах.
cwt(x,Fs)
