В этом примере вы демонстрируете экземпляр разрывов в шумных данных, представляемых более редко с помощью вейвлета Хаара чем тогда, когда с помощью вейвлета с большей поддержкой. Этот пример требует Signal Processing Toolbox.
Создайте шумную прямоугольную волну с 512 выборками. Постройте прямоугольную волну.
n = 512; t = 0:0.001:(n*0.001)-0.001; yn = square(2*pi*10*t)+0.02*randn(size(t)); plot(yn) grid on title('Noisy Signal')
Получите максимальное перекрытие дискретный вейвлет преобразовывает (MODWT) сигнала с помощью вейвлета haar
. Вейвлет haar
имеет поддержку длины, равной 1
modhaar = modwt(yn,'haar');
Получите анализ мультиразрешения из матрицы MODWT haar
и постройте детали первого уровня.
mrahaar = modwtmra(modhaar,'haar'); figure hs = stem(mrahaar(1,:)); grid on hs.Marker = 'none'; hs.ShowBaseLine = 'off'; title('First-Level MRA Details Using Haar Wavelet')
Получите MODWT сигнала при помощи вейвлета db4
. Вейвлет db4
имеет поддержку длины, равной 7.
moddb4 = modwt(yn,'db4');
Получите анализ мультиразрешения из матрицы MODWT db4
и постройте детали первого уровня. Разрывы представлены менее редко с помощью вейвлета db4
, чем вейвлет haar
.
mradb4 = modwtmra(modhaar,'db4'); figure hs = stem(mradb4(1,:)); grid on hs.Marker = 'none'; hs.ShowBaseLine = 'off'; title('First-Level MRA Details Using db4 Wavelet')