MATLAB.Exponenta
MATLAB и Simulink на русском
Технологии разработки и отладки
		сложных технических систем

MATLAB\MATLAB

В.Г.Потемкин "Справочник по MATLAB"
Анализ и обработка данных

  В оглавление книги \ К следующему разделу \ К предыдущему разделу

Преобразование Фурье

FFTSHIFT
Перегруппировка выходных массивов преобразований Фурье

Синтаксис:

            Y = fftshift(X)

Описание:

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

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

v fftshift(v)
1   2   3   4   5 4   5   1   2   3

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

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

X fftshift(X)
I II
III IV
IV III
II I

Пример:

Рассмотрим тот же пример, который рассматривался и для функции 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

58.4112 136.6337 140.1869 62.3762
58.4112 68.9109 140.1869 96.8317

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

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

image840.gif (2432 bytes) image841.gif (2256 bytes)
а) б)

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

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

Ссылки:

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

  В оглавление книги \ К следующему разделу \ К предыдущему разделу

 


Поиск по сайту:

Система Orphus

Яндекс.Метрика