Загрузите и постройте график одномерного сигнала.

load sumsin 
plot(sumsin)
title('Signal')

Выполните 3-уровневое вейвлет-разложение сигнала, используя вейвлет Daubechies порядка 2. Извлеките коэффициенты грубой шкалы приближения и коэффициенты детализации из разложения.

[c,l] = wavedec(sumsin,3,'db2');
approx = appcoef(c,l,'db2');
[cd1,cd2,cd3] = detcoef(c,l,[1 2 3]);

Постройте график коэффициентов.

subplot(4,1,1)
plot(approx)
title('Approximation Coefficients')
subplot(4,1,2)
plot(cd3)
title('Level 3 Detail Coefficients')
subplot(4,1,3)
plot(cd2)
title('Level 2 Detail Coefficients')
subplot(4,1,4)
plot(cd1)
title('Level 1 Detail Coefficients')

Обратитесь к разделу Поддержка GPU по версии (Parallel Computing Toolbox), чтобы узнать, какие графические процессоры поддерживаются.

Загрузите шумный сигнал Доплера. Поместите сигнал на графический процессор с помощью gpuArray. Сохраните текущий режим расширения.

load noisdopp
noisdoppg = gpuArray(noisdopp);
origMode = dwtmode('status','nodisp');

Использование dwtmode для изменения режима внутреннего абонента на нулевое заполнение. Получите трехуровневый DWT сигнала на графическом процессоре с помощью db4 вейвлет.

dwtmode('zpd','nodisp')
[c,l] = wavedec(noisdoppg,3,'db4');

Текущий режим расширения zpd не поддерживается для gpuArray вход. Поэтому DWT выполняется вместо этого с помощью sym режим расширения. Чтобы подтвердить это, установите режим внутреннего абонента на sym и возьмите DWT из noisdoppg, затем сравните с предыдущим результатом.

dwtmode('sym','nodisp')
[csym,lsym] = wavedec(noisdoppg,3,'db4');
[max(abs(c-csym)) max(abs(l-lsym))]

ans =

     0     0

Установите текущий режим расширения на per и получите трехуровневый DWT noisdopp. Режим расширения per поддерживается для gpuArray вход. Подтвердите, что результат отличается от sym результаты.

dwtmode('per','nodisp')
[cper,lper] = wavedec(noisdoppg,3,'db4');
[length(csym) ; length(cper)]

ans = 2×1

        1044
        1024

[lsym ; lper]

ans = 2×5

         134         134         261         515        1024
         128         128         256         512        1024

Восстановите режим расширения к исходной настройке.

dwtmode(origMode,'nodisp')