Exponenta Banner
exponenta event banner

littlewoodPaleySum

Сумма Литтлвуд-Пейли

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

Синтаксис

lpsum = littleweyPaleySum (sf)
lpsum = littleweyPaleySum (sf, fb)
[lpsum, f] = littleweyPaleySum (___)

Описание

lpsum = littlewoodPaleySum(sf) возвращает сумму Литлвуда-Пэли для 2-х банков фильтра в 2-й сети рассеивания небольшой волны sf. lpsum является матрицей M-by-N-by-Nfb, где M-by-N - размер матрицы заполненных фильтров, а Nfb - количество блоков фильтров.

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

lpsum = littlewoodPaleySum(sf,fb) возвращает сумму Littlewood-Paley для указанных банков фильтров fb. fb является положительным целым числом или вектором положительных целых чисел между 1 и numfilterbanks(sf) включительно. lpsum - матрица M-by-N-by-L, где L - количество уникальных элементов в fb.

пример

[lpsum,f] = littlewoodPaleySum(___) возвращает пространственные частоты для суммы Литлвуда-Пейли. f - матрица из двух столбцов с первым столбцом, содержащим пространственные частоты в направлении x, и вторым столбцом, содержащим пространственные частоты в направлении y.

Примеры

свернуть все

Открыть сценарий в реальном времени

В этом примере показано, как получить и отобразить сумму Littlewood-Paley сети рассеяния изображения.

Создайте сеть рассеяния с двумя банками фильтров и коэффициентами качества 2 и 1 соответственно.

sf = waveletScattering2('QualityFactors',[2 1]);

Получите суммы Литлвуда-Пейли и пространственные частоты двух банков фильтров. Отображение максимального значения сумм. Поскольку преобразование рассеяния является сжатым, суммы не превышают 1.

[lpsum,f] = littlewoodPaleySum(sf);
max(max(lpsum(:,:,1)))
ans = 1.0000

max(max(lpsum(:,:,2)))
ans = 1.0000

Отобразите сумму Литлвуда-Пейли для второго набора фильтров с нулевой частотой. Обратите внимание, что 2-й банк фильтра Morlet, используемый в преобразовании рассеивания, не разработан, чтобы захватить самые высокие пространственные частоты совместно в x-и y-направлениях.

f(f>1/2) = f(f>1/2)-1;
surf(fftshift(f(:,1)),fftshift(f(:,2)),fftshift(lpsum(:,:,2)))
shading interp
view(0,90)
xlabel('f_x')
ylabel('f_y')
colorbar
title('Q=1')

Figure contains an axes. The axes with title Q=1 contains an object of type surface.

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

свернуть все

Сеть рассеяния вейвлетного изображения, заданная как waveletScattering2 объект.

Фильтрация индекса банка в сети рассеяния изображения, заданного как положительное целое число или вектор положительных целых чисел от 1 до numfilterbanks(sf) включительно. Количество банков фильтров в sf равно числу указанных QualityFactors в sf.

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

свернуть все

Сумма Littlewood-Paley для банков фильтров в сети рассеяния изображений sf, возвращается в виде вещественной матрицы 3-D. lpsum является матрицей M-by-N-by-L, где M-by-N - размер матрицы заполненных фильтров, и L не превышает число блоков фильтров в sf.

Частоты для суммы Литлвуда-Пейли, возвращаемые в виде вещественно-значимой двухколонной матрицы. Частоты находятся в циклах на пиксель. Первый столбец f содержит пространственные частоты в направлении x, а второй столбец содержит пространственные частоты в направлении y. В этом соглашении преобразование Фурье является 1-периодическим в обеих переменных Фурье.

См. также

Представлен в R2019a

Документация инструментария вейвлет

Поддержка