Непрерывный анализ вейвлета

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

Figure contains an axes. The axes contains an object of type line.

Используйте cwt получить CWT сигнала и построить его scalogram. Величины компонентов синусоиды в шкале палитры являются по существу своими амплитудами даже при том, что они в различных шкалах.

cwt(x,Fs)

Figure contains an axes. The axes with title Magnitude Scalogram contains 3 objects of type image, line, area.

Смотрите также

|

Похожие темы

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