Поменять местами левую и правую половины вектора строки. Если вектор имеет нечетное число элементов, то средний элемент считается частью левой половины вектора.

Xeven = [1 2 3 4 5 6];
fftshift(Xeven)

ans = 1×6

     4     5     6     1     2     3

Xodd = [1 2 3 4 5 6 7];
fftshift(Xodd)

ans = 1×7

     5     6     7     1     2     3     4

При анализе частотных составляющих сигналов может быть полезно переместить нулевые частотные составляющие в центр.

Создание сигнала Sвычислите его преобразование Фурье и постройте график мощности.

fs = 100;               % sampling frequency
t = 0:(1/fs):(10-1/fs); % time vector
S = cos(2*pi*15*t);
n = length(S);
X = fft(S);
f = (0:n-1)*(fs/n);     %frequency range
power = abs(X).^2/n;    %power
plot(f,power)

Сместите нулевые частотные компоненты и постройте график нулевой центрированной мощности.

Y = fftshift(X);
fshift = (-n/2:n/2-1)*(fs/n); % zero-centered frequency range
powershift = abs(Y).^2/n;     % zero-centered power
plot(fshift,powershift)

Можно обработать несколько 1-D сигналов, представив их в виде строк в матрице. Затем используйте аргумент размерности для вычисления преобразования Фурье и сдвига нулевых частотных компонентов для каждой строки.

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

fs = 100;               % sampling frequency
t = 0:(1/fs):(10-1/fs); % time vector
S1 = cos(2*pi*15*t);
S2 = cos(2*pi*30*t);
n = length(S1);
A = [S1; S2];
X = fft(A,[],2);
f = (0:n-1)*(fs/n);     % frequency range
power = abs(X).^2/n;    % power
plot(f,power(1,:),f,power(2,:))

Сдвиньте нулевые частотные компоненты и постройте график нулевой центрированной мощности каждого сигнала.

Y = fftshift(X,2);
fshift = (-n/2:n/2-1)*(fs/n); % zero-centered frequency range
powershift = abs(Y).^2/n;     % zero-centered power
plot(fshift,powershift(1,:),fshift,powershift(2,:))