Exponenta Banner

wthrmngr

Диспетчер параметров порога

Свернуть все на странице

Синтаксис

thr = wthrmngr(opt,method,C,L)
thr = wthrmngr(opt,method,C,L,alpha)
thr = wthrmngr(opt,method,C,L,scale)
thr = wthrmngr(opt,method,swtdec,alpha)
thr = wthrmngr(opt,method,swtdec,scale)
thr = wthrmngr(opt,method,wpt)
thr = wthrmngr(opt,'rem_n0',X)

Описание

wthrmngr возвращает глобальный порог или зависящие от уровня пороги для основанных на вейвлете шумоподавления и сжатия. Функция выводит пороги из коэффициентов вейвлета при вейвлете разложении.

Пороги используются инструментами шумоподавления и сжатия Wavelet Toolbox™, такими как функции командной строки и приложение Wavelet Analyzer.

пример

thr = wthrmngr(opt,method,C,L) возвращает порог для [C, L] вейвлет сигнала или изображения для сжатия или обесценения. Для сигналов, [C, L] - выходы wavedec. Для изображений, [C, L] - выходы wavedec2.

пример

thr = wthrmngr(opt,method,C,L,alpha) возвращает значение [C, L] вейвлет вейвлет-разложения с использованием параметра разреженности alpha. Для сигналов, [C, L] - выход wavedec. Для изображений, [C, L] - выходы wavedec2.

Чтобы узнать больше о alpha, посмотрите wdcbm или wdcbm2 для сжатия, и wbmpen за шумоподавление.

thr = wthrmngr(opt,method,C,L,scale) возвращает значение [C, L] порог разложения вейвлет с использованием типа мультипликативного порогового перемасштабирования, заданного в scale. Для сигналов, [C, L] - выходы wavedec. Для изображений, [C, L] - выходы wavedec2.

The 'rigrsure', 'heursure', и 'minimaxi' методы шумоподавления применимы только к сигналам.

Чтобы узнать больше о мультипликативном перемасштабировании порога, см. wden.

thr = wthrmngr(opt,method,swtdec,alpha) возвращает зависящий от уровня порог для стационарного разложения вейвлет, swtdec, сигнала или изображения для денуилизации. alpha задает параметр разреженности (см. wbmpen). Для сигналов, swtdec - выходы swt. Для изображений, swtdec - выход swt2.

Пороги получают из подмножества коэффициентов в стационарном разложении вейвлета. Для получения дополнительной информации смотрите Выбор коэффициента.

пример

thr = wthrmngr(opt,method,swtdec,scale) возвращает зависящий от уровня порог для стационарного разложения вейвлет с помощью типа мультипликативного перемасштабирования порога, заданного в scale. Для сигналов, swtdec - выходы swt. Для изображений, swtdec - выходы swt2.

Пороги получают из подмножества коэффициентов в стационарном разложении вейвлета. Для получения дополнительной информации смотрите Выбор коэффициента.

The 'rigrsure', 'heursure', и 'minimaxi' шумоподавление применяются только к сигналам.

Чтобы узнать больше о мультипликативном перемасштабировании порога, см. wden.

пример

thr = wthrmngr(opt,method,wpt) возвращает глобальный порог для вейвлета разложения пакетов, wpt, сигнала или изображения для сжатия или денуилизации.

thr = wthrmngr(opt,'rem_n0',X) возвращает глобальный порог для сжатия сигнала или изображения, X, с использованием заданных опции вейвлетов и метода 'rem_n0'.

Если opt является 'dw1dcompGBL' или 'dw2dcompGBL'пороги основаны на наименьших коэффициентах вейвлет шкалы, полученных с использованием вейвлета Haar. Если opt является 'wp1dcompGBL' или 'wp2dcompGBL'пороги основаны на коэффициентах вейвлет-пакета с самой мелкой шкалой, полученных с использованием вейвлета Haar.

Примеры

свернуть все

Открыть Live Script

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

load noisdopp
plot(noisdopp)
grid on
title('Noisy Signal')

Figure contains an axes. The axes with title Noisy Signal contains an object of type line.

Сгенерируйте вейвлет-разложение уровня 5 сигнала с шумом, используя вейвлет Daubechies порядка 4. Постройте график коэффициентов.

[c,l] = wavedec(noisdopp,5,'db4');
plot(c)
grid on
title('Wavelet Coefficients')

Figure contains an axes. The axes with title Wavelet Coefficients contains an object of type line.

Определите глобальный порог для сжатия сигнала.

thr = wthrmngr('dw1dcompGBL','bal_sn',c,l);

Индекс первого коэффициента детализации вейвлета в c является l(1)+1. Примените порог ко всем коэффициентам детализации. Постройте график пороговых коэффициентов. Заметьте, что большинство коэффициентов было установлено на 0.

c(l(1)+1:end) = c(l(1)+1:end).*(c(l(1)+1:end)>thr);
plot(c)
grid on
title('Thresholded Coefficients')

Figure contains an axes. The axes with title Thresholded Coefficients contains an object of type line.

Восстановите сигнал из пороговых коэффициентов. Постройте график реконструкции.

xrec = waverec(c,l,'db4');
plot(xrec)
grid on
title('Compressed Signal')

Figure contains an axes. The axes with title Compressed Signal contains an object of type line.

Открыть Live Script

Сжатие изображения с помощью стратегии Бирже-Массарта.

Загрузите изображение и добавьте белый Гауссов шум. В целях воспроизводимости установите значение по умолчанию для случайного начального числа.

rng default
load sinsin
x = X+18*randn(size(X));

Получите 2-D дискретное вейвлет-преобразование до уровня 3 с помощью наименее асимметричного вейвлета Daubechies с 4 моментами исчезновения. Получите пороги сжатия с помощью стратегии Бирже-Массарта с параметром разреженности, alpha, равный 2.

[C,L] = wavedec2(x,3,'sym4');
alpha = 2;
THR = wthrmngr('dw2dcompLVL','scarcehi',C,L,alpha);

Сжать изображение и отобразить результат. 

xd = wdencmp('lvd',x,'sym4',3,THR,'s');
image(X)
title('Original Image')

Figure contains an axes. The axes with title Original Image contains an object of type image.

figure
image(x)
title('Noisy Image')

Figure contains an axes. The axes with title Noisy Image contains an object of type image.

figure
image(xd)
title('Compressed Image')

Figure contains an axes. The axes with title Compressed Image contains an object of type image.

Открыть Live Script

Этот пример использует зависящий от уровня порог, выведенный из коэффициентов вейвлета в каждой шкале, чтобы реализовать жесткое пороговое значение со стационарным преобразованием вейвлета.

Загрузите сигнал шумных блоков. Получите стационарное вейвлет-преобразование до уровня 5 с помощью вейвлет Haar. 

load noisbloc
L = 5;
swc = swt(noisbloc,L,'haar');

Создайте копию коэффициентов преобразования вейвлета. Определите универсальный порог Донохо-Джонстоуна на основе коэффициентов детализации для каждой шкалы. Использование 'mln' опция, wthrmngr возвращает вектор 1-by-L с каждым элементом, равным универсальному порогу для соответствующей шкалы. 

swcnew = swc;
ThreshML = wthrmngr('sw1ddenoLVL','sqtwolog',swc,'mln');

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

for jj = 1:L
    swcnew(jj,:) = wthresh(swc(jj,:),'h',ThreshML(jj));
end

Инвертируйте стационарное вейвлет на пороговых коэффициентах, swcnew. Постройте график исходного сигнала и деноизированного сигнала для сравнения. 

noisbloc_denoised = iswt(swcnew,'haar');
plot(noisbloc)
hold on
plot(noisbloc_denoised,'r','linewidth',2)
legend('Original','Denoised')

Figure contains an axes. The axes contains 2 objects of type line. These objects represent Original, Denoised.

Открыть Live Script

Обнулите сигнал с шумом путем применения глобального порога к вейвлету структуре разложения пакетов.

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

load noisdopp
plot(noisdopp)
grid on
title('Noisy Signal')

Figure contains an axes. The axes with title Noisy Signal contains an object of type line.

Сгенерируйте вейвлет-пакетное разложение уровня 3 сигнала с шумом, используя вейвлет Daubechies порядка 4.

T = wpdec(noisdopp,3,'db4');

Определите глобальный порог для шумоподавления сигнала.

thr = wthrmngr('wp1ddenoGBL','sqtwologuwn',T);

Получите листья из вейвлета дерева разложения пакетов T и применить порог к листьям. Используйте жесткое пороговое значение.

T1 = T;
sorh = 'h';
cfs = read(T,'data');
cfs = wthresh(cfs,sorh,thr);
T1 = write(T1,'data',cfs);

Восстановите деноцированный сигнал из пороговых коэффициентов. Постройте график реконструкции.

xrec = wprec(T1);
plot(xrec)
grid on
title('Denoised Signal')

Figure contains an axes. The axes with title Denoised Signal contains an object of type line.

Открыть Live Script

Этот пример использует независимый от уровня порог, основанный на наименьших коэффициентах вейвлет шкалы, чтобы реализовать жесткое пороговое значение со стационарным вейвлет-преобразованием.

Загрузите сигнал шумных блоков. Получите стационарное вейвлет-преобразование до уровня 5 с помощью вейвлет Haar. 

load noisbloc
L = 5;
swc = swt(noisbloc,L,'haar');

Создайте копию коэффициентов преобразования вейвлета. Определите универсальный порог Донохо-Джонстоуна на основе коэффициентов детализации первого уровня. Использование 'sln' опция, wthrmngr возвращает вектор 1-by-L с каждым элементом, равным тому же значению. Возьмем среднее значение вектора, чтобы получить скалярный порог. 

swcnew = swc;
ThreshSL = mean(wthrmngr('sw1ddenoLVL','sqtwolog',swc,'sln'));

Используйте универсальный порог для реализации жесткого порога. Тот же порог применяется к коэффициентам вейвлета на каждом уровне. 

for jj = 1:L
    swcnew(jj,:) = wthresh(swc(jj,:),'h',ThreshSL);
end

Инвертируйте стационарное вейвлет на пороговых коэффициентах, swcnew. Постройте график исходного сигнала и деноизированного сигнала для сравнения. 

noisbloc_denoised = iswt(swcnew,'haar');
plot(noisbloc)
hold on
plot(noisbloc_denoised,'r','linewidth',2)
legend('Original','Denoised')

Figure contains an axes. The axes contains 2 objects of type line. These objects represent Original, Denoised.

Входные параметры

свернуть все

Тип и размерности сжатия или шумоподавления, заданные как одно из значений, перечисленных в следующих таблицах. wthrmngr возвращает пороги, соответствующие заданной опции.

С помощью дискретного вейвлет или вейвлет-пакетного разложения данных можно сжать или обесценить эти данные. При стационарном вейвлет данных можно только денурировать данные.

Для объяснения того, какие коэффициенты используются для определения порогов, смотрите Выбор коэффициентов.

1-D дискретных опций разложения вейвлетов

В этих опциях X является сигналом, коэффициенты вейвлета находятся в векторе C, и длины векторов коэффициентов в L. Аргумент alpha является параметром разреженности и scale определяет мультипликативное перемасштабирование порога.

Для получения дополнительной информации о разложении вейвлетов см. wavedec. Чтобы узнать больше о alpha и scale, см. wdcbm и wden соответственно.

opt

Описание

Допустимые синтаксисы
'dw1dcompGBL'

1-D сжатие с помощью глобального порога

  • thr = wthrmngr('dw1dcompGBL','rem_n0',X)

  • thr = wthrmngr('dw1dcompGBL','bal_sn',C,L)

'dw1dcompLVL'

1-D сжатие с помощью зависящих от уровня порогов

  • thr = wthrmngr('dw1dcompLVL','scarcehi',C,L,alpha), где 2,5 < alpha < 10

  • thr = wthrmngr('dw1dcompLVL','scarceme',C,L,alpha), где 1,5 < alpha < 2.5

  • thr = wthrmngr('dw1dcompLVL','scarcelo',C,L,alpha), где 1 < alpha < 2

'dw1ddenoLVL'

1-D шумоподавление с использованием зависящих от уровня порогов

  • thr = wthrmngr('dw1ddenoLVL','sqtwolog',C,L,scale)

  • thr = wthrmngr('dw1ddenoLVL','rigrsure',C,L,scale)

  • thr = wthrmngr('dw1ddenoLVL','heursure',C,L,scale)

  • thr = wthrmngr('dw1ddenoLVL','minimaxi',C,L,scale)

  • thr = wthrmngr('dw1ddenoLVL','penalhi',C,L,alpha), где 2,5 < alpha < 10

  • thr = wthrmngr('dw1ddenoLVL','penalme',C,L,alpha), где 1,5 < alpha < 2.5

  • thr = wthrmngr('dw1ddenoLVL','penallo',C,L,alpha), где 1 < alpha < 2

2-D дискретные опции разложения вейвлетов

В этих опциях X являются данными, коэффициенты вейвлета находятся в векторе C, и размер матриц коэффициентов в L. Аргумент alpha является параметром разреженности и scale определяет мультипликативное перемасштабирование порога.

Для получения дополнительной информации о разложении вейвлетов см. wavedec2. Чтобы узнать больше о alpha и scale, см. wdcbm2 и wden соответственно.

opt

Описание

Допустимые синтаксисы
'dw2dcompGBL'

2-D сжатие с помощью глобального порога

  • thr = wthrmngr('dw2dcompGBL','rem_n0',X)

  • thr = wthrmngr('dw2dcompGBL','bal_sn',C,L)

  • thr = wthrmngr('dw2dcompGBL','sqrtbal_sn',C,L)

'dw2dcompLVL'

2-D сжатие с помощью зависящих от уровня порогов

  • thr = wthrmngr('dw2dcompLVL','scarcehi',C,L,alpha), где 2,5 < alpha < 10

  • thr = wthrmngr('dw2dcompLVL','scarceme',C,L,alpha), где 1,5 < alpha < 2.5

  • thr = wthrmngr('dw2dcompLVL','scarcelo',C,L,alpha), где 1 < alpha < 2

'dw2ddenoLVL'

2-D шумоподавление с использованием зависящих от уровня порогов

  • thr = wthrmngr('dw2ddenoLVL','sqrtbal_sn',C,L)

  • thr = wthrmngr('dw2ddenoLVL','penalhi',C,L,alpha), где 2,5 < alpha < 10

  • thr = wthrmngr('dw2ddenoLVL','penalme,C,L,alpha), где 1,5 < alpha < 2.5

  • thr = wthrmngr('dw2ddenoLVL','penallo,C,L,alpha), где 1 < alpha < 2

  • thr = wthrmngr('dw2ddenoLVL','sqtwolog',C,L,scale)

1-D Вейвлет разложение пакетов Опций

В этих опциях X является сигналом и wpt - вейвлет структура разложения пакета сигнала.

Для получения дополнительной информации о вейвлете разложении пакетов смотрите wpdec.

opt

Описание

Допустимые синтаксисы
'wp1dcompGBL'

1-D сжатие с помощью глобального порога

  • thr = wthrmngr('wp1dcompGBL','rem_n0',X)

  • thr = wthrmngr('wp1dcompGBL','bal_sn',wpt)

'wp1ddenoGBL'

1-D шумоподавление с использованием глобального порога

  • thr = wthrmngr('wp1ddenoGBL','sqtwologuwn',wpt)

  • thr = wthrmngr('wp1ddenoGBL','sqtwologswn',wpt)

  • thr = wthrmngr('wp1ddenoGBL','bal_sn',wpt)

  • thr = wthrmngr('wp1ddenoGBL','penalhi',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 6.25.

  • thr = wthrmngr('wp1ddenoGBL','penalme',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 2.

  • thr = wthrmngr('wp1ddenoGBL','penallo',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 1.5.

2-D Вейвлеты разложение пакетов Опций

В этих опциях X - данные и wpt - вейвлет структура разложения пакетов данных.

Для получения дополнительной информации о вейвлете разложении пакетов смотрите wpdec2.

opt

Описание

Допустимые синтаксисы
'wp2dcompGBL'

2-D сжатие с помощью глобального порога

  • thr = wthrmngr('wp2dcompGBL','rem_n0',X)

  • thr = wthrmngr('wp2dcompGBL','bal_sn',wpt)

  • thr = wthrmngr('wp2dcompGBL','sqrtbal_sn',wpt)

'wp2ddenoGBL'

2-D шумоподавление с использованием глобального порога

  • thr = wthrmngr('wp2ddenoGBL','sqtwologuwn',wpt)

  • thr = wthrmngr('wp2ddenoGBL','sqtwologswn',wpt)

  • thr = wthrmngr('wp2ddenoGBL','sqrtbal_sn',wpt)

  • thr = wthrmngr('wp2ddenoGBL','penalhi',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 6.25.

  • thr = wthrmngr('wp2ddenoGBL','penalme',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 2.

  • thr = wthrmngr('wp2ddenoGBL','penallo',wpt)

    wpbmpen функция используется с параметром настройки ALPHA = 1.5.

1-D стационарного разложения вейвлет- Опции

Шумоподавление с использованием зависящих от уровня порогов является единственной опцией, доступной для 1-D стационарного разложения вейвлет, swtdec. В этой опции alpha является параметром разреженности и scale определяет мультипликативное перемасштабирование порога.

Для получения дополнительной информации о стационарном разложении вейвлетов см. swt. Чтобы узнать больше о alpha и scale, см. wbmpen и wden соответственно.

optДопустимые синтаксисы
'sw1ddenoLVL'

  • thr = wthrmngr('sw1ddenoLVL','sqtwolog',swtdec,scale)

  • thr = wthrmngr('sw1ddenoLVL','rigrsure',swtdec,scale)

  • thr = wthrmngr('sw1ddenoLVL','heursure',swtdec,scale)

  • thr = wthrmngr('sw1ddenoLVL','minimaxi',swtdec,scale)

  • thr = wthrmngr('sw1ddenoLVL','penalhi',swtdec,alpha), где 2,5 < alpha < 10

  • thr = wthrmngr('sw1ddenoLVL','penalme',swtdec,alpha), где 1,5 < alpha < 2.6

  • thr = wthrmngr('sw1ddenoLVL','penallo',swtdec,alpha), где 1 < alpha < 2

Пороги основаны на подмножестве коэффициентов в стационарном разложении вейвлет. Дополнительную информацию см. в разделе Выбор коэффициентов.

2-D стационарного разложения вейвлет- Опции

Шумоподавление с использованием зависящих от уровня порогов является единственной опцией, доступной для 2-D стационарного разложения вейвлет, swtdec. В этой опции alpha является параметром разреженности и scale определяет мультипликативное перемасштабирование порога.

Для получения дополнительной информации о стационарном разложении вейвлетов см. swt2. Чтобы узнать больше о alpha и scale, см. wbmpen и wden соответственно.

optДопустимые синтаксисы
'sw2ddenoLVL'

  • thr = wthrmngr('sw2ddenoLVL','sqrtbal_sn',swtdec)

  • thr = wthrmngr('sw2ddenoLVL','penalhi',swtdec,alpha) где 2,5 < alpha < 10

  • thr = wthrmngr('sw2ddenoLVL','penalme',swtdec,alpha) где 1,5 < alpha < 2.5

  • thr = wthrmngr('sw2ddenoLVL','penallo',swtdec,alpha) где 1 < alpha < 2

  • thr = wthrmngr('sw2ddenoLVL','sqtwolog',swtdec,scale)

Пороги основаны на подмножестве коэффициентов в стационарном разложении вейвлет. Дополнительную информацию см. в разделе Выбор коэффициентов.

Метод порогового значения, заданный как одно из значений, перечисленных здесь.

methodОписание
'scarcehi'Использует стратегию Бирже-Массарта для определения порогов.
'scarceme'Использует стратегию Бирже-Массарта для определения порогов.
'scarcelo'Использует стратегию Бирже-Массарта для определения порогов.
'sqtwolog'Использует универсальный порог фиксированной формы. См. 'sqtwolog' опция в wden.
'sqtwologuwn'Использует универсальный порог фиксированной формы. См. 'sqtwolog' опция в wden при использовании с 'sln'опция.
'sqtwologswn'Использует универсальный порог фиксированной формы. См. 'sqtwolog' опция в wden при использовании с 'mln' опция.
'rigsure'Использует правило мягкой оценки порога основанное на объективной оценке риска Штейна. См. 'SURE' опция в wdenoise.
'heursure'Использует смесь 'rigsure' и 'sqtwolog'. См. 'heursure' опция в wden.
'minimaxi' Использует выбранный фиксированный порог, который приводит к минимаксной эффективности. См. 'Minimax' опция в wdenoise.
'penalhi'Используется для определения стратегии Бирже-Массарта для определения порогов.
'penalme'Используется для определения стратегии Бирже-Массарта для определения порогов.
'penallo'Используется для определения стратегии Бирже-Массарта для определения порогов.
'rem_n0'Возвращает порог, близкий к 0. Типичное THR значение median(abs(coefficients)).
'bal_sn'Возвращает порог, такой что проценты сохраненной энергии и количество нулей одинаковы.
'sqrtbal_sn'Возвращает порог, равный квадратному корню значения, таким образом, что проценты сохраненной энергии и количество нулей одинаковы.

Типы данных: char

Входные данные, заданные как действительный вектор или вещественная матрица.

Типы данных: double

Вейвлеты вейвлет-расширения данных, которые будут сжиматься или деноизменяться, заданные как действительный вектор. Если данные одномерны, C - выходы wavedec. Если данные двумерны, C - выход wavedec2.

Пример: [C,L] = wavedec(randn(1,1024),3,'db4')

Типы данных: double

Размер вейвлетов вейвлет-расширения сжимаемого или деноизменяемого сигнала или изображения, заданный как вектор или матрица положительных целых чисел.

Для сигналов, L - выходы wavedec. Для изображений, L - выходы wavedec2.

Пример: [C,L] = wavedec(randn(1,1024),3,'db4')

Типы данных: double

Параметр разреженности, используемый для сжатия или шумоподавления данных, заданный как положительная скалярная величина больше 1 и меньше 10. Посмотрите wdcbm, wdcbm2, и wbmpen для получения дополнительной информации.

Типы данных: double

Мультипликативное перемасштабирование порога, заданное как одно из следующего:

  • 'one' - Нет перемасштабирования

  • 'sln' - Пересмотр с использованием одной оценки уровня шума на основе коэффициентов первого уровня

  • 'mln' - Пересмотр с использованием зависящей от уровня оценки шума уровня

Для получения дополнительной информации см. wden.

Стационарная структура вейвлета разложения данных, подлежащих сжатию или деноизации, заданная как действительная матрица. Если данные одномерны, swtdec - выходы swt. Если данные двумерны, swtdec - выход swt2.

Пример: swtdec = swt2(randn(256),3,'db1')

Типы данных: double

Вейвлет разложение пакетов данных, которые будут сжиматься или обесцениваться. Если данные одномерны, wpt - выходы wpdec. Если данные двумерны, wpt - выходы wpdec2.

Пример: wpt = wpdec(randn(1,1024),5,'db1')

Выходные аргументы

свернуть все

Порог, возвращенный как действительный скаляр для глобальных порогов или действительный вектор или матрица для зависящих от уровня порогов.

Типы данных: double

Совет

  • Чтобы обесценить 1-D сигналы, рассмотрите использование Wavelet Signal Denoiser. Приложение визуализирует и денофицирует реальные сигналы 1-D с помощью параметров по умолчанию. Можно также сравнить результаты. В сложение можно также воссоздать деноизмененный сигнал в вашей рабочей области, сгенерировав MATLAB® скрипт, который использует wdenoise функция.

Алгоритмы

свернуть все

Выбор коэффициента

Критически дискретизированное разложение вейвлет или вейвлет-пакета включает децимирование коэффициентов в 2 раза на каждой стадии разложения. Десятикратное уменьшение не происходит при неразрешенном стационарном вейвлет.

wthrmngr выводит шумоподавление и пороги сжатия из коэффициентов вейвлета. Для критически дискретизированного разложения вейвлет или вейвлет-пакета опция и метод определяют, используются ли все вейвлет-коэффициенты или только самые мелкие масштабные коэффициенты.

Для стационарного разложения вейвлета, wthrmngr всегда использует подмножество коэффициентов вейвлета. При вычислении деноизирующих порогов N-уровневое стационарное разложение вейвлет, алгоритм сначала подсемеряет вейвлет-коэффициенты на уровне k в 2 разаk, для k = 1,…,N. Алгоритм использует этот подмножество коэффициентов, чтобы определить пороги. Большинство коэффициентов в стационарном разложении вейвлета не рассматриваются.

Стратегия Бирже-Массарта

Стратегия Бирже-Массарта для определения порогов зависит от нескольких различных параметров. Вы задаете разложение вейвлет и метод порогового значения. Можно также задать параметр разреженности, alpha, или определенного мультипликативного порогового перемасштабирования, scale. Основываясь на ваших входах, wthrmngr выводит необходимые параметры Бирже-Массарта. Параметры зависят от размерности сигнала и общего числа, N, коэффициентов в самой грубой шкале вейвлет.

Если метод порогового значения 'scarcehi', 'scarceme', или 'scarcelo', wthrmngr выполняет либо wdcbm или wdcbm2. Если метод порогового значения 'penalhi', 'penalme', или 'penallo', затем wthrmngr выполняет либо wbmpen или wpbmpen.

Метод порогового значения

Описание

'scarcehi'

  • Если сигнал 1-D, то wdcbm используется с входным параметром M = N.

  • Если сигнал 2-D, то wdcbm2 используется с M = 4*N.

'scarceme'

  • Если сигнал 1-D, то wdcbm используется с входным параметром M = 3*N/2.

  • Если сигнал 2-D, то wdcbm2 используется с входным параметром с M = 16*N/3.

'scarcelo'

  • Если сигнал 1-D, то wdcbm используется с входным параметром M = 2 N.

  • Если сигнал 2-D, то wdcbm2 используется с входным параметром M = 32*N/3.

'penalhi'

  • Если вход является вейвлет, то wbmpen используется с ALPHA = 5* (3* alpha+1)/8.

  • Если вход является вейвлетом разложением пакета, то wpbmpen используется ALPHA = 6.25.

'penalme'

  • Если вход является вейвлет, то wbmpen используется с ALPHA = (alpha+5)/8.

  • Если вход является вейвлетом разложением пакета, то wpbmpen используется ALPHA = 2.

'penallo'

  • Если вход является вейвлет, то wbmpen используется с ALPHA = (alpha+3)/4.

  • Если вход является вейвлетом разложением пакета, то wpbmpen используется ALPHA = 1.5.

Ссылки

[1] Бирге, Л. и П. Массарт. От выбора модели до адаптивной оценки. Festschrift for Lucien Le Cam: Research Papers in Probability and Statistics (E. Torgersen, D. Pollard, and G. Yang, eds.). Нью-Йорк: Springer-Verlag, 1997, pp. 55-88.

См. также

Приложения

Функции

Представлено до R2006a

Документация по Wavelet Toolbox

Поддержка