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

Материал из MachineLearning.

Перейти к: навигация, поиск
Статья написана с использованием LLM DeepSeek-V3-0324 и проверена участником Nikita Elкhin 06:36, 18 июля 2026 (MSD)


Содержание

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

Преобразование Фурье (часто обозначается \mathcal{F}) — интегральное преобразование, ставящее в соответствие функции вещественной или комплексной переменной её спектральное представление в базисе гармонических колебаний. В анализе данных и прикладной математике оно является фундаментальным инструментом перехода от временно́го (пространственного) описания сигнала к частотному, раскрывая внутреннюю структуру процессов и позволяя эффективно решать дифференциальные уравнения, выполнять фильтрацию и сжатие информации[1].

Определение и геометрический смысл

Пусть f: \mathbb{R} \to \mathbb{C} — абсолютно интегрируемая функция, f \in L^1(\mathbb{R}). Прямое непрерывное преобразование Фурье определяется выражением \hat{f}(\omega) = \mathcal{F}\{f(t)\} = \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt, \qquad \omega \in \mathbb{R}. Обратное преобразование имеет вид f(t) = \mathcal{F}^{-1}\{\hat{f}(\omega)\} = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i\omega t} d\omega, \qquad t \in \mathbb{R}, где i — мнимая единица, \omega — круговая частота.

Принципиально важна геометрическая интерпретация. В гильбертовом пространстве L^2(\mathbb{R}) квадратично-интегрируемых функций скалярное произведение определяется как \langle f, g \rangle = \int_{-\infty}^{\infty} f(t) \overline{g(t)} dt. Система комплексных экспонент \{e^{i\omega t}\}_{\omega \in \mathbb{R}} образует континуальный ортонормированный базис. В отличие от конечномерного случая условие ортонормированности записывается через дельта-функцию Дирака: \frac{1}{2\pi} \int_{-\infty}^{\infty} e^{i(\omega_1 - \omega_2)t} dt = \delta(\omega_1 - \omega_2). Преобразование Фурье \hat{f}(\omega) — это в точности координаты функции f в этом базисе, то есть её проекция на гармонику частоты \omega: \hat{f}(\omega) = \langle f(t), e^{i\omega t} \rangle = \int_{-\infty}^{\infty} f(t) \overline{e^{i\omega t}} dt = \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt. Таким образом, модуль |\hat{f}(\omega)| показывает амплитуду, а аргумент — фазу соответствующей частотной составляющей. Обратное преобразование восстанавливает функцию как суперпозицию всех таких проекций с весами \frac{d\omega}{2\pi}.

Это геометрическое понимание приводит к сохранению скалярного произведения (равенству Парсеваля–Планшереля): \langle f, g \rangle = \frac{1}{2\pi} \langle \hat{f}, \hat{g} \rangle, а при g = f — к сохранению энергии: \int_{-\infty}^{\infty} |f(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |\hat{f}(\omega)|^2 d\omega.

О нормировках

В литературе встречаются различные способы нормировки преобразования Фурье. Приведённая выше форма (с множителем 1/(2\pi) только у обратного преобразования) удобна в анализе дифференциальных уравнений и теореме о свёртке, так как не порождает дополнительных констант при дифференцировании. Не менее распространена симметричная нормировка \hat{f}(\omega) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt, f(t) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} \hat{f}(\omega) e^{i\omega t} d\omega, при которой оператор \mathcal{F} становится унитарным в L^2(\mathbb{R}), а равенство Парсеваля принимает особенно простой вид \|f\| = \|\hat{f}\|[1]. В данной статье используется несимметричная нормировка, если не оговорено иное.

Условия существования

Классическое интегральное определение требует абсолютной интегрируемости f \in L^1(\mathbb{R}). При этом \hat{f}(\omega) непрерывна и стремится к нулю на бесконечности (лемма Римана–Лебега). Для принадлежности образа тому же пространству накладывают более жёсткие условия, например, принадлежность функциям из класса Шварца \mathcal{S}(\mathbb{R}) (быстро убывающим гладким функциям), где преобразование Фурье является автоморфизмом.

Если же функция принадлежит L^2(\mathbb{R}), но не L^1(\mathbb{R}), преобразование определяют через предельный переход в смысле L^2, а равенство Парсеваля гарантирует его корректность. Для обобщённых функций умеренного роста (например, констант, тригонометрических функций, дельта-функций) используется аппарат теории распределений; в этом случае спектр содержит сингулярные составляющие, такие как дельта-пики[1].

Основные свойства

Все приведённые свойства справедливы в смысле L^1, L^2 или пространства Шварца, если не оговорено иное.

  • Линейность: \mathcal{F}\{\alpha f(t) + \beta g(t)\} = \alpha \hat{f}(\omega) + \beta \hat{g}(\omega).
  • Сдвиг во временно́й области: \mathcal{F}\{f(t - t_0)\} = e^{-i\omega t_0} \hat{f}(\omega).
  • Модуляция (сдвиг по частоте): \mathcal{F}\{e^{i\omega_0 t} f(t)\} = \hat{f}(\omega - \omega_0).
  • Масштабирование: для ненулевого вещественного a \mathcal{F}\{f(at)\} = \frac{1}{|a|}\hat{f}\!\left(\frac{\omega}{a}\right).
  • Сопряжение: \mathcal{F}\{\overline{f(t)}\} = \overline{\hat{f}(-\omega)}. Для вещественной f это даёт эрмитову симметрию: \hat{f}(-\omega) = \overline{\hat{f}(\omega)}.
  • Дифференцирование: \mathcal{F}\{f^{(n)}(t)\} = (i\omega)^n \hat{f}(\omega), где f^{(n)} — n-я производная (при условии её существования и интегрируемости).
  • Интегрирование (в смысле обобщённых функций умеренного роста): \mathcal{F}\left\{\int_{-\infty}^{t} f(\tau) d\tau\right\} = \frac{\hat{f}(\omega)}{i\omega} + \pi \hat{f}(0) \delta(\omega).
  • Теорема о свёртке: если h(t) = (f * g)(t) = \int_{-\infty}^{\infty} f(\tau) g(t-\tau) d\tau, то \hat{h}(\omega) = \hat{f}(\omega) \hat{g}(\omega).
  • Двойственность (симметрия): если \hat{f}(\omega) = \mathcal{F}\{f(t)\}, то \mathcal{F}\{\hat{f}(t)\} = 2\pi f(-\omega). Это свойство позволяет получать новые пары преобразований.
  • Равенство Парсеваля–Планшереля: \int_{-\infty}^{\infty} f(t) \overline{g(t)} dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega) \overline{\hat{g}(\omega)} d\omega.

Таблица характерных преобразований

Приведём несколько классических пар, демонстрирующих геометрический смысл и свойства.

f(t) \hat{f}(\omega) Примечание
e^{-a |t|}, \; a > 0 \frac{2a}{a^2 + \omega^2} Двустороннее экспоненциальное затухание
e^{-a t^2}, \; a > 0 \sqrt{\frac{\pi}{a}} e^{-\frac{\omega^2}{4a}} Гауссиан переходит в гауссиан, сжатие/растяжение
\operatorname{rect}\left(\frac{t}{T}\right) \equiv \begin{cases} 1, & |t| \leq T/2, \\ 0, & |t| > T/2 \end{cases} T \operatorname{sinc}\!\left(\frac{\omega T}{2}\right), \operatorname{sinc} x = \frac{\sin x}{x} Прямоугольный импульс — прототип идеального фильтра
\frac{\sin(\Omega t)}{\pi t} \operatorname{rect}\!\left(\frac{\omega}{2\Omega}\right) Дуальность с предыдущим случаем
\delta(t) (дельта-функция Дирака) 1 Бесконечно узкий импульс содержит все частоты с равной амплитудой
1 2\pi \delta(\omega) Постоянная составляющая — нулевая частота
e^{i\omega_0 t} 2\pi \delta(\omega - \omega_0) Чистый тон — дельта-пик в частотной области
\cos(\omega_0 t) \pi [\delta(\omega - \omega_0) + \delta(\omega + \omega_0)] Вещественная гармоника

Принцип неопределённости Габора–Гейзенберга

Из геометрической интерпретации следует фундаментальное ограничение: функция и её преобразование Фурье не могут быть одновременно хорошо локализованы. В теории сигналов этот факт формализуется теоремой Габора–Гейзенберга (частотно-временной принцип неопределённости). Пусть f \in L^2(\mathbb{R}) и \|f\| = 1. Определим центры «масс» и дисперсии: \mu_t = \int t|f(t)|^2 dt, \quad \sigma_t^2 = \int (t - \mu_t)^2 |f(t)|^2 dt, \mu_\omega = \frac{1}{2\pi} \int \omega |\hat{f}(\omega)|^2 d\omega, \quad \sigma_\omega^2 = \frac{1}{2\pi} \int (\omega - \mu_\omega)^2 |\hat{f}(\omega)|^2 d\omega. Тогда справедливо неравенство \sigma_t \sigma_\omega \ge \frac{1}{2}. Равенство достигается исключительно на функциях вида гауссова импульса f(t) = C e^{-a (t - t_0)^2} e^{i\omega_0 t}, что делает гауссиан одновременно оптимальным по частотно-временнóй локализации[1]. Эта теорема подчёркивает, что сжатие сигнала во временно́й области неизбежно ведёт к расширению его спектра, и наоборот.

Связь с рядом Фурье

Периодический сигнал с периодом T может быть разложен в ряд Фурье: f(t) = \sum_{k=-\infty}^{\infty} c_k e^{i k \omega_0 t}, \quad \omega_0 = \frac{2\pi}{T}. Его преобразование Фурье, понимаемое в обобщённом смысле, является суммой эквидистантных дельта-импульсов: \hat{f}(\omega) = 2\pi \sum_{k=-\infty}^{\infty} c_k \delta(\omega - k\omega_0). Таким образом, непрерывный спектр периодической функции сосредоточен на дискретном множестве частот, кратных основной частоте. При устремлении периода к бесконечности расстояние между спектральными линиями \omega_0 \to 0, и огибающая коэффициентов переходит в непрерывный образ \hat{f}(\omega). Эта предельная связь объясняет, почему преобразование Фурье рассматривают как непрерывный аналог разложения по ортогональному базису.

Преобразование Фурье в многомерном случае

Естественное обобщение на функции f: \mathbb{R}^n \to \mathbb{C} задаётся формулой \hat{f}(\mathbf{\omega}) = \int_{\mathbb{R}^n} f(\mathbf{x}) e^{-i \langle \mathbf{\omega}, \mathbf{x} \rangle} d\mathbf{x}, \qquad \mathbf{\omega} \in \mathbb{R}^n, где \langle \mathbf{\omega}, \mathbf{x} \rangle = \sum_{j=1}^n \omega_j x_j. Все одномерные свойства переносятся на многомерный случай: сдвиг порождает линейный фазовый множитель, дифференцирование — умножение на компоненты i\omega_j, а теорема о свёртке сохраняет вид. Особую ценность многомерное преобразование представляет в задачах обработки изображений, томографии и спектрального анализа полей.

Оконное преобразование Фурье

Для анализа сигналов, спектральные характеристики которых изменяются во времени, применяют оконное преобразование Фурье (STFT). Оно получается умножением исходной функции на скользящее окно w(t-\tau), локализующее сигнал в окрестности момента \tau: \text{STFT}\{f(t)\}(\tau, \omega) = \int_{-\infty}^{\infty} f(t) w(t - \tau) e^{-i\omega t} dt. Таким образом, STFT даёт разложение по частотам, «привязанное» к временно́й координате. Выбор оконной функции (Хэмминга, Ханна, Гаусса и др.) определяет компромисс между разрешением по времени и по частоте в соответствии с принципом неопределённости Габора–Гейзенберга: сужение окна улучшает временну́ю локализацию, но расширяет спектральную полосу. Квадрат модуля STFT называют спектрограммой — это широко используемое в обработке речи, музыки и вибродиагностике частотно-временное представление.

Применения в анализе данных и прикладной математике

Хотя цифровая обработка использует в основном дискретные версии (см. Тригонометрическая интерполяция), само непрерывное преобразование Фурье служит теоретической базой для построения алгоритмов и интерпретации результатов:

  • Спектральный анализ. Непрерывный спектр |\hat{f}(\omega)|^2 (спектральная плотность мощности) количественно описывает распределение энергии сигнала по частотам, что востребовано в вибродиагностике, сейсмологии и анализе временных рядов.
  • Решение дифференциальных уравнений. Свойство \mathcal{F}\{f^{(n)}\} = (i\omega)^n \hat{f} превращает линейные дифференциальные уравнения с постоянными коэффициентами в алгебраические, что является классическим методом математической физики.
  • Теория сигналов и систем. Понятие частотной характеристики линейной стационарной системы есть в точности преобразование Фурье её импульсной характеристики; анализ устойчивости и фильтрации опирается на свойства \hat{f}(\omega).
  • Оптика и волновая физика. Дифракция Фраунгофера описывается преобразованием Фурье от апертуры, что делает спектральный подход основным в задачах формирования изображений и голографии.
  • Теория вероятностей. Характеристическая функция случайной величины определяется как преобразование Фурье её плотности распределения, а мультипликативное свойство свёртки отвечает суммированию независимых величин.
  • Частотно-временной анализ. Сформулированный выше принцип неопределённости устанавливает фундаментальные рамки для компромисса между разрешением по времени и по частоте, что прямо учитывается при выборе окон в кратковременном преобразовании Фурье и вейвлет-анализе.

Смотрите также

Примечания

Источник — «http://www.machinelearning.ru/wiki/index.php?title=%D0%9F%D1%80%D0%B5%D0%BE%D0%B1%D1%80%D0%B0%D0%B7%D0%BE%D0%B2%D0%B0%D0%BD%D0%B8%D0%B5_%D0%A4%D1%83%D1%80%D1%8C%D0%B5»