2D преобразования Фурье

fft2 функционируйте преобразовывает 2D данные в пространство частоты. Например, можно преобразовать 2D оптическую маску, чтобы показать ее дифракционный шаблон.

Двумерное преобразование Фурье

Следующая формула задает дискретное преобразование Фурье Y m-by-n матричный X.

$Y_{p + 1, q + 1} = \sum_{j = 0}^{m - 1} \sum_{k = 0}^{n - 1} ω_{m}^{j p} ω_{n}^{k q} X_{j + 1, k + 1}$

_ωm и _ωn являются комплексными корнями из единицы, заданными следующими уравнениями.

$\begin{array}{l} ω_{m} = e^{- 2 π i / m} \\ ω_{n} = e^{- 2 π i / n} \end{array}$

i является мнимой единицей, p и j являются индексами, которые запускаются от 0 до m –1, и q и k являются индексами, которые запускаются от 0 до n –1. Индексы для X и Y смещены 1 в этой формуле, чтобы отразить матричные индексы в MATLAB^®.

Вычисление 2D преобразования Фурье X эквивалентно первому вычислению 1D преобразования каждого столбца X и затем взятия 1D преобразования каждой строки результата. Другими словами, команда fft2(X) эквивалентно Y = fft(fft(X).').'.

2D дифракционный шаблон

Скрипт Open Live Script

В оптике преобразование Фурье может использоваться, чтобы описать дифракционный шаблон, произведенный инцидентом плоской волны на оптической маске с маленькой апертурой [1]. Этот пример использует fft2 функция на оптической маске, чтобы вычислить ее дифракционный шаблон.

Создайте логический массив, который задает оптическую маску с маленькой, круговой апертурой.

n = 2^10;                 % size of mask
M = zeros(n);
I = 1:n; 
x = I-n/2;                % mask x-coordinates 
y = n/2-I;                % mask y-coordinates
[X,Y] = meshgrid(x,y);    % create 2-D mask grid
R = 10;                   % aperture radius
A = (X.^2 + Y.^2 <= R^2); % circular aperture of radius R
M(A) = 1;                 % set mask elements inside aperture to 1
imagesc(M)                % plot mask
axis image

Figure contains an axes. The axes contains an object of type image.

Используйте fft2 вычислить 2D преобразование Фурье маски и использовать fftshift функционируйте, чтобы перестроить выход так, чтобы нулевая частотная составляющая находилась в центре. Постройте получившиеся дифракционные частоты шаблона. Синий указывает, что маленькие амплитуды и желтый указывают на большие амплитуды.

DP = fftshift(fft2(M));
imagesc(abs(DP))
axis image

Figure contains an axes. The axes contains an object of type image.

Чтобы улучшить детали областей с маленькими амплитудами, постройте 2D логарифм дифракционного шаблона. Очень маленькие амплитуды затронуты числовой ошибкой округления, и прямоугольная сетка вызывает радиальную асимметрию.

imagesc(abs(log2(DP)))
axis image

Figure contains an axes. The axes contains an object of type image.

Ссылки

[1] Фаулз, G. R. Введение в современную оптику. Нью-Йорк: Дувр, 1989.

Документация MATLAB

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

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

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

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