Реконструировать частотно-локализованное приближение данных землетрясения Кобе. Извлеките информацию из CWT для частот в диапазоне [0,030, 0,070] Гц.
load kobeПолучение CWT данных.
[wt,f] = cwt(kobe,1);
Восстановите данные землетрясения, добавив среднее значение сигнала обратно в преобразованные данные.
xrec = icwt(wt,f,[0.030 0.070],'SignalMean',mean(kobe));Постройте график и сравните исходные данные и данные для частот в диапазоне [0,030, 0,070] Гц.
subplot(2,1,1) plot(kobe) grid on title('Original Data') subplot(2,1,2) plot(xrec) grid on title('Bandpass Filtered Reconstruction [0.030 0.070] Hz')
![Figure contains 2 axes. Axes 1 with title Original Data contains an object of type line. Axes 2 with title Bandpass Filtered Reconstruction [0.030 0.070] Hz contains an object of type line.](../../examples/wavelet/win64/cwtftdemo_01.png)
С CWT можно также использовать периоды времени вместо частоты. Загрузите данные El Nino и получите их CWT, указав период времени в годах.
load ninoairdata
[cfs,period] = cwt(nino,years(1/12));Получить обратный CWT за годы 2 до 8.
xrec = icwt(cfs,period,[years(2) years(8)]);
Постройте график CWT восстановленных данных. Обратите внимание на отсутствие энергии вне диапазона периодов от 2 до 8 лет.
figure cwt(xrec,years(1/12))
Сравните исходные данные с восстановленными данными за годы 2- 8.
figure subplot(2,1,1) plot(nino) grid on title('Original Data') subplot(2,1,2) plot(xrec) grid on title('El Nino Data - Years 2-8')
