Эффект поддержки вейвлета на шумных данных

В этом примере вы демонстрируете экземпляр разрывов в шумных данных, представляемых более редко с помощью вейвлета Хаара чем тогда, когда с помощью вейвлета с большей поддержкой. Этот пример требует 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')

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

| |