В этом примере показано, как выполнять частотно-временной анализ с использованием непрерывного вейвлет-преобразования (CWT). Непрерывный вейвлет-анализ обеспечивает временной/частотный анализ сигналов и изображений. Программное обеспечение Vavelet Toolbox™ имеет как функции командной строки, так и интерактивные функции для поддержки непрерывного вейвлет-анализа 1-D сигналов.
Создайте сигнал, состоящий из двух синусоид с частотами 100 и 50 Гц и белого шума. Опора двух синусоид разобщена. Волна синуса на 100 Гц начинается в 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)
