fftshift
Сдвиньте компонент нулевой частоты к центру спектра
Синтаксис
Описание
Y = fftshift( переставляет преобразование Фурье X)X путем сдвига компонента нулевой частоты к центру массива.
Если
Xявляется вектором, тогдаfftshiftменяет левую и правую половиныX.Если
Xявляется матрицей, тогдаfftshiftменяет первый квадрантXс третьим, и вторым квадрантом с четвертым.Если
Xявляется многомерным массивом, затемfftshiftменяет полупространстваXвдоль каждой размерности.
Примеры
Сдвиг Вектора элементы
Замените левую и правую половины вектора-строки. Если вектор имеет нечетное число элементов, то средний элемент считается частью левой половины вектора.
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
1-D сдвига
При анализе частотных составляющих сигналов может быть полезно сместить компоненты нулевой частоты в центр.
Создайте сигнальное 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,:))
Входные параметры
X - Входной массив
вектор | матрица | многомерный массив
Входной массив, заданный как вектор, матрица или многомерный массив.
Типы данных: double | single | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | logical
Поддержка комплексного числа: Да
dim - Размерность для работы
положительный целочисленный скаляр
Размерность для работы, заданная как положительный целочисленный скаляр Если значение не задано, то fftshift свопы по всем размерностям.
Рассмотрим матрицу входа
Xc. Область операцииfftshift(Xc,1)меняет половины каждого столбцаXc.
Рассмотрим матрицу
Xr. Область операцииfftshift(Xr,2)меняет половины каждой строкиXr.
Типы данных: double | single | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | logical