Непрерывный анализ вейвлета
В этом примере показано, как выполнить частотно-временной анализ с помощью непрерывного вейвлета преобразовывает (CWT). Непрерывный анализ вейвлета обеспечивает масштаб времени/частотно-временной анализ сигналов и изображений. Программное обеспечение Wavelet Toolbox™ имеет и командную строку и интерактивную функциональность, чтобы поддержать непрерывный анализ вейвлета 1D сигналов.
Создайте сигнал, состоящий из двух синусоид с частотами 100 и 50 Гц и белым шумом. Поддержка этих двух синусоид является непересекающейся. Синусоида на 100 Гц начинается в t = 0 и имеет длительность 1 секунды. Синусоида на 100 Гц имеет амплитуду 2. Синусоида на 50 Гц начинается в три секунды и имеет длительность двух секунд. Синусоида на 50 Гц имеет амплитуду 1. Частота дискретизации составляет 1 кГц. Длина сигнала является 5 000 выборок.
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 сигнала и построить его scalogram. Величины компонентов синусоиды в шкале палитры являются по существу своими амплитудами даже при том, что они в различных шкалах.
cwt(x,Fs)
