Синтаксис:

            Y = fftshift(X)

Описание:

Функция Y = fftshift(X) перегруппировывает выходные массивы функций fft и fft2, размещая нулевую частоту в центре спектра.

Если v - одномерный массив, то выполняется циклическая перестановка правой и левой его половины:

Если X - двумерный массив, то меняются местами квадранты: I « IV и II « III:

          X = [' I '    ' II ';...
                 ' III '   ' IV ';]
         fftshift(X)

Пример:

Рассмотрим тот же пример, который рассматривался и для функции fft, но введем постоянную составляющую с уровнем 0.3:

           t = 0:0.001:0.6;
           x = sin(2 * pi * 50 * t) + sin(2 * pi * 120 * t);
           y = x + 2 * randn(size(t)) + 0.3;
           plot(y(1:50)), grid

Применим одномерное преобразование Фурье для сигнала y и построим график спектральной плотности. Здесь можно выделить 3 частоты, на которых амплитуда спектра максимальна. Это частоты 0, 58.4 и 140.1 Гц (рис. а).

           Y = fft(y);
           Pyy = Y.*conj(Y)/512;
           f = 1000 * (0:255)/512;
           figure(1), plot(f, Pyy(1:256)), grid
           ginput

Выполним операции

           Y = fftshift(Y);
           Pyy = Y.*conj(Y)/512;
           figure(2), plot(Pyy), grid

а)	б)

Из анализа рис. б) следует, что нулевая частота сместилась в середину спектра.

Сопутствующие функции: FFT, FFT2, Signal Processing Toolbox [1].

Ссылки:

1. Signal Processing Toolbox User’s Guide. Natick: The MathWorks, Inc., 1993.