Определите глобальный порог шумоподавления по умолчанию для $$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')