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

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)

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

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

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,:))