Преобразование Фурье: ключевые свойства

Преобразование Фурье переводит сигнал из временной области в частотную: вместо вопроса «какова амплитуда в момент t?» оно отвечает на вопрос «на каких частотах сосредоточена энергия сигнала?». Именно благодаря этому переходу задачи фильтрации, свёртки и сжатия данных становятся алгебраическими операциями над спектром, а не сложными интегральными уравнениями. Разберём основные свойства преобразования Фурье с доказательствами и поймём, почему каждое из них важно на практике. Чтобы сразу пощупать, как меняется спектр при изменении параметров сигнала, воспользуйтесь калькулятором ниже.

Определение и пара преобразований

Прямое преобразование Фурье (ПФ) определяется как:

$$X(f) = \int_{-\infty}^{+\infty} x(t)\, e^{-j 2\pi f t}\, dt,$$

а обратное - как:

$$x(t) = \int_{-\infty}^{+\infty} X(f)\, e^{+j 2\pi f t}\, df.$$

Здесь $x(t)$ - вещественный или комплексный сигнал, $X(f)$ - его спектральная плотность, $f$ - частота в герцах. Образ $X(f)$ в общем случае комплексный: $X(f) = |X(f)|\, e^{j\,\arg X(f)}$. Модуль $|X(f)|$ называют спектром амплитуды, аргумент $\arg X(f)$ - спектром фазы.

Рассчитать для своей задачи

Свойство линейности

Формулировка. Для любых констант $a, b \in \mathbb{C}$ и сигналов $x_1(t), x_2(t)$ с образами $X_1(f), X_2(f)$:

$$\mathcal{F}\{a\, x_1(t) + b\, x_2(t)\} = a\, X_1(f) + b\, X_2(f).$$

Доказательство следует из линейности интеграла: интеграл от суммы равен сумме интегралов, константу можно выносить за знак интеграла.

Практический смысл. Спектр суперпозиции сигналов - это суперпозиция спектров. Именно поэтому, зная FT-образы базовых функций (прямоугольника, гауссианы, косинуса), можно без вычислений записать образ любой их линейной комбинации. Например, для суммы двух косинусов:

$$x(t) = A_1 \cos(2\pi f_1 t) + A_2 \cos(2\pi f_2 t) \;\xrightarrow{\mathcal{F}}\; \frac{A_1}{2}\bigl[\delta(f-f_1)+\delta(f+f_1)\bigr] + \frac{A_2}{2}\bigl[\delta(f-f_2)+\delta(f+f_2)\bigr].$$

Спектр состоит из четырёх дельта-функций - по одной паре для каждой гармоники. В калькуляторе выше столбики на частотах $f_1$ и $f_2$ имеют высоту $A_1/2$ и $A_2/2$ именно по этой формуле.

Рассчитать для своей задачи

Свойство сдвига по времени

Формулировка. Если $x(t) \xrightarrow{\mathcal{F}} X(f)$, то сдвиг сигнала на $t_0$ секунд:

$$x(t - t_0) \;\xrightarrow{\mathcal{F}}\; X(f)\, e^{-j 2\pi f t_0}.$$

Доказательство. Подставим $\tau = t - t_0$ ($t = \tau + t_0$, $dt = d\tau$):

$$\int_{-\infty}^{+\infty} x(t - t_0)\, e^{-j 2\pi f t}\, dt = \int_{-\infty}^{+\infty} x(\tau)\, e^{-j 2\pi f (\tau + t_0)}\, d\tau = e^{-j 2\pi f t_0} \int_{-\infty}^{+\infty} x(\tau)\, e^{-j 2\pi f \tau}\, d\tau = X(f)\, e^{-j 2\pi f t_0}.$$

Практический смысл. Сдвиг по времени не меняет спектр амплитуды: $|X(f)\, e^{-j 2\pi f t_0}| = |X(f)|$, поскольку $|e^{j\theta}| = 1$ при любом $\theta$. Изменяется только фаза: к фазе каждой спектральной составляющей прибавляется линейный член $-2\pi f t_0$. Это прямой путь к пониманию того, почему задержка в тракте обработки не портит форму спектра амплитуды, но сдвигает фазу - и почему для звукового сигнала важна линейная фазовая характеристика фильтра.

В калькуляторе: передвиньте ползунок $t_0$ - высота столбиков в спектре амплитуды не изменится, а нижний график фазы покажет набег $-2\pi f t_0$ в градусах для каждой частоты.

Свойство масштабирования по времени

Формулировка. При масштабировании аргумента ($a \neq 0$, вещественное):

$$x(a t) \;\xrightarrow{\mathcal{F}}\; \frac{1}{|a|}\, X\!\left(\frac{f}{a}\right).$$

Практический смысл. Сжатие сигнала по времени (увеличение $|a| > 1$) расширяет его спектр и уменьшает амплитуду - и наоборот. Это принцип неопределённости ширины: чем короче сигнал во времени, тем шире его спектр в частотной области. Принцип неопределённости Гейзенберга в квантовой механике - частный случай той же математики.

$\Delta t \cdot \Delta f \geq \frac{1}{4\pi}$

где $\Delta t$ и $\Delta f$ - «ширины» (среднеквадратические отклонения) сигнала и его спектра.

Свойство дифференцирования

По времени:

$$\frac{d^n x(t)}{dt^n} \;\xrightarrow{\mathcal{F}}\; (j 2\pi f)^n\, X(f).$$

По частоте (дуальное):

$$(-j 2\pi t)^n\, x(t) \;\xrightarrow{\mathcal{F}}\; \frac{d^n X(f)}{df^n}.$$

Эти соотношения позволяют заменить дифференциальные уравнения алгебраическими в частотной области. Дифференциальный оператор $d/dt$ превращается в умножение на $j 2\pi f$ - сложная операция становится умножением числа. Именно на этом построено решение линейных ОДУ с постоянными коэффициентами через передаточную функцию.

Теорема о свёртке

Свёртка во временной области эквивалентна умножению в частотной:

$$(x_1 * x_2)(t) = \int_{-\infty}^{+\infty} x_1(\tau)\, x_2(t-\tau)\, d\tau \;\xrightarrow{\mathcal{F}}\; X_1(f) \cdot X_2(f).$$

И дуально: умножение во временной области - свёртка в частотной:

$$x_1(t)\cdot x_2(t) \;\xrightarrow{\mathcal{F}}\; X_1(f) * X_2(f).$$

Практический смысл. Линейный фильтр в системе обработки сигналов описывается импульсной характеристикой $h(t)$. Его реакция на входной сигнал - свёртка: $y(t) = x(t) * h(t)$. В частотной области это превращается в умножение: $Y(f) = X(f) \cdot H(f)$, где $H(f)$ - передаточная функция. Это и есть причина, по которой частотная фильтрация вычислительно эффективна: свёртку длиной $N$ можно заменить двумя БПФ + поэлементным умножением ($O(N \log N)$ вместо $O(N^2)$).

Рассчитать для своей задачи

Теорема Парсеваля

Энергия сигнала одинакова в обоих представлениях:

$$\int_{-\infty}^{+\infty} |x(t)|^2\, dt = \int_{-\infty}^{+\infty} |X(f)|^2\, df.$$

Величину $|X(f)|^2$ называют спектральной плотностью мощности (PSD). Теорема Парсеваля - это инструмент для проверки: если энергия в частотной области не совпадает с временной, где-то ошибка в расчёте.

Частые ошибки

• Путаница с нормировкой. Существуют три стандарта: с $2\pi f$, с $\omega = 2\pi f$, с нормировочным множителем $1/\sqrt{2\pi}$. Перед работой уточняйте конвенцию.
• Знак в экспоненте сдвига. Задержка сигнала на $t_0$ даёт множитель $e^{-j 2\pi f t_0}$ (минус в показателе). Опережение - $e^{+j 2\pi f t_0}$.
• Неверный вывод масштабирования. Не забывайте, что при $a < 0$ в знаменателе стоит $|a|$, иначе амплитуда уйдёт в отрицательное значение.
• Игнорирование единиц. Если частота в рад/с ($\omega$), то $X(\omega) = \int x(t) e^{-j\omega t} dt$ и формулы свойств имеют другой вид (нет множителей $2\pi$).
• Свёртка с $\delta$-функцией. Часто забывают: $x(t) * \delta(t - t_0) = x(t - t_0)$ - свёртка с $\delta$-функцией - это сдвиг. Если это знать, половина задач на свёртку решается мгновенно.

FAQ

Почему сдвиг по времени не меняет модуль спектра? Потому что множитель $e^{-j 2\pi f t_0}$ имеет модуль 1 при любом вещественном $t_0$ и $f$: $|e^{j\theta}| = \sqrt{\cos^2\theta + \sin^2\theta} = 1$. Информация о задержке уходит целиком в фазу, но не в амплитуду спектра.

Что такое спектральная утечка и как она связана со свойствами ПФ? Спектральная утечка возникает при вычислении ДПФ на конечном отрезке - фактически сигнал умножается на прямоугольное окно. По теореме о свёртке это означает, что спектр оконного сигнала - свёртка истинного спектра с преобразованием Фурье прямоугольника (sinc-функцией). «Юбочки» sinc и дают утечку мощности в соседние частоты. Для борьбы с ней используют гладкие оконные функции (Ханн, Хэмминг), у которых sinc-«юбочки» малы.

Верно ли свойство линейности для дискретного преобразования Фурье (ДПФ)? Да, ДПФ линейно: $\text{DFT}\{a x_1 + b x_2\} = a\, X_1 + b\, X_2$. Все четыре основных свойства - линейность, сдвиг, масштабирование, свёртка - имеют точные аналоги для ДПФ, с поправкой на цикличность (сдвиг оборачивается по кругу).

Коротко

Преобразование Фурье обладает четырьмя ключевыми свойствами. Линейность: образ суммы равен сумме образов. Сдвиг по времени: $x(t-t_0) \to X(f)\, e^{-j 2\pi f t_0}$, модуль спектра не меняется, фаза получает линейный набег. Масштабирование: $x(at) \to |a|^{-1} X(f/a)$, сжатие по времени - расширение по частоте. Свёртка: умножение в одной области соответствует свёртке в другой - именно на этом строится быстрая цифровая фильтрация. Знание этих свойств позволяет решать задачи на спектральный анализ без вычисления интегралов, работая только с алгебраическими соотношениями.