Определите глобальный порог шумоподавления по умолчанию для $$N$$(0,1) вход белого шума. Создайте $$N$$(0,1) вход белого шума. Смените режим расширения DWT на периодический. Установите генератор случайных чисел в начальные настройки по умолчанию для воспроизводимых результатов.

origmode = dwtmode('status','nodisplay');
dwtmode('per','nodisp')
rng default
x = randn(512,1);

Использование ddencmp чтобы получить глобальный порог по умолчанию для шумоподавления вейвлет. Продемонстрировать, что порог равен универсальному порогу Донохо и Джонстоуна, масштабируемому устойчивой оценкой отклонения.

[thr,sorh,keepapp] = ddencmp('den','wv',x);
[A,D] = dwt(x,'db1');
noiselev = median(abs(D))/0.6745;
thresh = sqrt(2*log(length(x)))*noiselev;

Сравните значение переменной thr к значению thresh.

thr

thr = 3.3639

thresh

thresh = 3.3639

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

dwtmode(origmode,'nodisplay')

Определите порог глобального сжатия по умолчанию для $$N$$(0,1) вход белого шума.

Создайте $$N$$(0,1) вход белого шума. Установите периодический режим расширения DWT. Установите генератор случайных чисел в начальные настройки по умолчанию для воспроизводимых результатов.

origmode = dwtmode('status','nodisplay');
dwtmode('per','nodisp')
rng default
x = randn(512,1);

Использование ddencmp с 'cmp' и 'wp' Входные параметры, чтобы вернуть порог глобального сжатия по умолчанию для вейвлета преобразования пакета.

[thr,sorh,keepapp,crit] = ddencmp('cmp','wp',x)

thr = 0.6424

sorh = 
'h'

keepapp = 1

crit = 
'threshold'

Сравните со значениями по умолчанию, возвращенными для шумоподавления.

[thr,sorh,keepapp,crit] = ddencmp('den','wp',x)

thr = 4.1074

sorh = 
'h'

keepapp = 1

crit = 
'sure'

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

dwtmode(origmode,'nodisplay')