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