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

Нажмите на фортепиано сразу три клавиши, и инструмент выпустит в воздух одну спутанную волну давления. До барабанной перепонки она доходит как одна колеблющаяся линия, а не как три. Тем не менее внутреннее ухо без всяких усилий разделяет эту линию обратно на три ноты, а подготовленный музыкант назовет аккорд. Где-то между воздухом и сознанием смешанный сигнал был разложен на чистые частоты.

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

Большинство студентов знакомится с преобразованием Фурье как с формулой для заучивания: интеграл с комплексной экспонентой внутри и без всяких очевидных причин для этого. Это руководство выстраивает картину с самого основания. Недавно мы разбирали близкое ему преобразование Лапласа, и к концу статьи вы увидите, что перед нами одна и та же машина с разным положением регулятора.

Два описания одного и того же

Подумайте о рецепте торта. Готовый торт можно описать, перечислив то, что вы наблюдаете в каждой его точке, слой за слоем, крошка за крошкой. А можно описать его списком ингредиентов: 200 граммов муки, 3 яйца, 100 граммов сахара. Оба описания содержат один и тот же торт. Первое говорит о том, где что находится, второе о том, из чего он сделан.

У сигнала есть ровно такие же два описания:

• Временная область перечисляет значение сигнала в каждое мгновение: давление на барабанной перепонке на 0,001 секунде, на 0,002 секунде и так далее. Именно в таком виде сигналы приходят к нам из реального мира.
• Частотная область перечисляет ингредиенты: сколько в сигнале синусоиды на 440 герц, сколько на 523 герц, сколько на каждой возможной частоте. Именно это и выдает преобразование Фурье.

Открытие Фурье, впервые опубликованное в 1822 году при изучении распространения тепла, состояло в том, что такие списки ингредиентов полны. Любой разумный сигнал, каким бы изломанным он ни был, можно точно восстановить из суммы гладких синусоид. Никакого третьего ингредиента, спрятанного где-то в стороне, не существует. Синусов достаточно.

Почему алфавитом служат именно синусоиды

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

Возьмите производную синусоиды, и вы получите косинус, то есть ту же волну, сдвинутую на четверть периода. Возьмите производную еще раз, и получится синус со знаком минус: снова та же форма, только масштабированная. Как мы разбирали в статье о тригонометрии, синусоида является окружностью, если смотреть на нее сбоку, а поворот окружности никогда не превращает ее во что-то другое.

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

В этом и заключается вся польза частотной области. Если вы знаете, что система делает с каждой чистой частотой, и умеете разложить любой вход на чистые частоты, вы знаете, что система делает с любым входом. Трудная задача «что произойдет с мостом, когда на него въедет грузовик» превращается в список легких задач: «что происходит с мостом на 1 герце, на 2 герцах, на 3 герцах».

Машина намотки

Теперь к самому механизму. Пусть у вас есть сигнал и вы хотите узнать, сколько в нем спрятано конкретной частоты, скажем, 3 циклов в секунду. Вот физическая процедура, которая отвечает на этот вопрос.

Возьмите график вашего сигнала, то есть линию, вытянутую вдоль оси времени. Теперь намотайте эту линию на окружность со скоростью 3 оборота на каждую секунду сигнального времени. Там, где сигнал велик, намотанная линия уходит далеко от центра, там, где он мал, она жмется к центру. Получается закорючка, нарисованная вокруг начала координат.

Теперь найдите центр масс этой закорючки.

Если в сигнале нет ничего на 3 герцах, максимумы и минимумы сигнала попадают на случайные углы окружности. Они гасят друг друга, и центр масс остается рядом с началом координат. Но если составляющая на 3 герцах в сигнале есть, ее пики оборот за оборотом ложатся на одну и ту же сторону окружности. Они накапливаются, и центр масс уходит от начала координат.

Повторите это для каждой возможной скорости намотки, постройте зависимость расстояния до центра масс от скорости намотки, и вы нарисуете преобразование Фурье. Пики на этом графике соответствуют частотам, присутствующим в сигнале. В этом вся идея. Все остальное является обозначениями.

Формула по частям

Машина намотки, записанная в символах, выглядит так:

$$\hat{f}(\omega) = \int_{-\infty}^{\infty} f(t)\, e^{-i\omega t}\, dt$$

Каждая часть соответствует своему шагу процедуры.

Намотка. Множитель $e^{-i\omega t}$ отвечает за вращение. По формуле Эйлера $e^{-i\omega t} = \cos(\omega t) - i\sin(\omega t)$, то есть это точка, движущаяся по единичной окружности по часовой стрелке с угловой частотой ω. С ростом t эта точка вращается. Ее вращение и есть намотка, а ω служит скоростью намотки.

Наматывание сигнала. Произведение $f(t)\, e^{-i\omega t}$ является той самой закорючкой. Умножение вращающейся точки единичной длины на значение сигнала f(t) растягивает ее наружу или притягивает внутрь. Большие значения сигнала отбрасывают точку далеко от центра, малые подтягивают ее обратно. Произведение вычерчивает в точности ту закорючку из предыдущего раздела.

Центр масс задается интегралом. Складывая положение закорючки в каждое мгновение, от начала времен до их конца, и рассматривая результат как одно комплексное число, мы получаем точку, чье расстояние от начала координат измеряет, насколько сильно сигнал резонирует на частоте ω. Строго говоря, это сумма, а не среднее, поэтому величина растет вместе с длительностью сигнала, но форма спектра остается той же.

Результат $\hat{f}(\omega)$ является комплексным числом для каждого ω, и обе его части важны. Модуль задает амплитуду этой частоты: сколько данного ингредиента в рецепте. Аргумент задает фазу: в какой момент цикла этот ингредиент вступает. Два сигнала могут иметь одинаковые амплитуды на всех частотах и звучать совершенно по-разному, потому что их фазы различаются.

Почему экспонента, а не синус

Студенты часто спрашивают, почему ядром служит комплексная экспонента, а не просто синусоида. Ведь можно, казалось бы, умножить сигнал на $\sin(\omega t)$ и проинтегрировать.

Проблема заключается в фазе. Если в сигнале есть косинус на частоте ω, а вы проверяете его синусом на той же частоте, произведение в среднем дает ноль, потому что синус и косинус разведены на четверть периода. Вы сделаете вывод, что частота отсутствует, хотя она присутствует в полной мере. Чтобы поймать любую фазу, понадобились бы две проверки, одна синусом и одна косинусом, а затем объединение результатов.

Комплексная экспонента выполняет обе проверки сразу. Ее действительная часть является проверкой косинусом, а мнимая часть проверкой синусом, и все это упаковано в один вращающийся объект. С какой бы фазой ни пришла частота, какая-то часть вращения с ней совпадет. Модуль получившегося комплексного числа равен тому, что дал бы объединенный результат двух проверок, а его аргумент бесплатно сообщает фазу. В статье о мнимых числах мы описывали i как поворот, и здесь этот поворот выполняет честную инженерную работу.

Отрицательные частоты не являются ошибкой

Интеграл берется по всем ω, включая отрицательные значения, и преобразование действительного сигнала всегда показывает зеркальное отражение в отрицательной области. Это многих смущает. Что вообще должна физически означать отрицательная частота?

Взгляните на формулу Эйлера еще раз. Действительный косинус является средним двух вращений:

$$\cos(\omega t) = \frac{e^{i\omega t} + e^{-i\omega t}}{2}$$

Первое слагаемое вращается против часовой стрелки со скоростью ω, второе вращается по часовой стрелке с той же скоростью. Вращение по часовой стрелке со скоростью ω преобразование записывает как частоту минус ω. Действительный сигнал не способен различить эти два направления, поэтому он всегда содержит оба в равной мере, а спектр получается симметричным.

Для действительных сигналов отрицательная половина не несет новой информации, и на практике инженеры часто рисуют только положительную сторону. Но сохранение обеих половин делает формулы симметричными и позволяет тому же преобразованию работать с комплексными сигналами, например в радиотехнике и томографии, где два направления вращения действительно различаются.

Обратный путь

Превращение рецепта обратно в торт называется обратным преобразованием Фурье:

$$f(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega)\, e^{i\omega t}\, d\omega$$

Читайте это как инструкцию: для каждой частоты ω возьмите количество $\hat{f}(\omega)$ этого ингредиента, раскрутите его против часовой стрелки с этой скоростью и сложите все вместе. Знак в показателе сменился с минуса на плюс, потому что теперь вы разматываете, а не наматываете. Множитель 2π служит бухгалтерией, поправкой на то, что полный оборот составляет 2π радиан. Некоторые учебники делят его симметрично, ставя $1/\sqrt{2\pi}$ в оба преобразования, а некоторые прячут его целиком, измеряя частоту в герцах, а не в радианах в секунду. Физика от этого не меняется, меняется только константа.

Пара формул почти идеально симметрична. Эта симметрия является глубокой подсказкой: время и частота представляют собой два взгляда на один объект, и ни один из них не более фундаментален, чем другой.

Производные снова становятся умножением

Преобразование Лапласа оправдало себя тем, что превращает производные в умножение на s. Преобразование Фурье проделывает тот же трюк с iω.

Продифференцируйте формулу обратного преобразования по t. Единственное место, где встречается t, находится внутри $e^{i\omega t}$, и производная выносит оттуда множитель iω:

$$\frac{d}{dt} f(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \bigl(i\omega\, \hat{f}(\omega)\bigr)\, e^{i\omega t}\, d\omega$$

Значит, преобразование Фурье функции $f'(t)$ равно просто $i\omega\, \hat{f}(\omega)$. Взятие производной по времени умножает каждую частотную составляющую на iω: высокие частоты усиливаются, низкие подавляются, и каждая составляющая поворачивается на четверть оборота. Именно поэтому дифференцирование делает шумные сигналы еще более шумными. Шум живет на высоких частотах, а производная умножает их на большие числа.

Та же логика превращает дифференциальное уравнение с постоянными коэффициентами в алгебраическое уравнение относительно ω, и именно так физики решают волновое уравнение и уравнение теплопроводности. Ровно для этого Фурье и придумал свой инструмент.

Фурье и Лапласа: одна машина, два положения регулятора

Поставьте два преобразования рядом:

$$F(s) = \int_0^{\infty} f(t)\, e^{-st}\, dt, \qquad \hat{f}(\omega) = \int_{-\infty}^{\infty} f(t)\, e^{-i\omega t}\, dt$$

В статье о Лапласе мы описывали s как панель с двумя ручками, $s = \sigma + i\omega$: сигма отвечает за экспоненциальное затухание, омега за колебания. Поставьте ручку затухания на ноль, так что σ = 0 и s = iω, и лапласово ядро $e^{-st}$ превратится в фурье-ядро $e^{-i\omega t}$. Преобразование Фурье является преобразованием Лапласа, вычисленным вдоль мнимой оси s-области.

Второе отличие касается нижнего предела. Лаплас стартует с t = 0, потому что он создан для причинно-следственных систем с моментом включения и начальными условиями. Фурье идет по всему времени, потому что он создан для устойчивых сигналов, которые были всегда и будут всегда: выдержанная нота, несущая волна, ровный гул.

Что выбрать, следует из самого вопроса. Переходные процессы, переключения и начальные условия: Лаплас. Спектры, фильтры и установившиеся колебания: Фурье. Ручка затухания как раз и позволяет Лапласу работать с неограниченно растущими сигналами, которые интеграл Фурье не в состоянии переварить, а двусторонний интеграл как раз и позволяет Фурье описывать сигналы без начала.

Компромисс, от которого не уйти

Последнее свойство объясняет удивительно много физики. Сигнал, узкий во времени, широк по частоте, а сигнал, узкий по частоте, широк во времени. Получить оба качества сразу невозможно.

У чистой синусоиды одна частота, поэтому ее спектр представляет собой единственный пик. Но чистая синусоида длится вечно. Обрежьте ее до одной секунды, и спектр расплывется вокруг исходной частоты, обрежьте до одной миллисекунды, и он расплывется гораздо шире. В пределе мгновенный щелчок содержит все частоты поровну, и поэтому хлопок в ладоши звучит как белый шум, а камертон звучит как тон.

Картина намотки делает это очевидным. Короткая вспышка сигнала при любой скорости намотки обернется вокруг окружности всего несколько раз, поэтому у центра масс нет шанса аккуратно погаситься. Множество близких скоростей намотки дают похожий умеренный результат. Длинный сигнал наматывается тысячи раз, и гашение переживает только точно совпадающая скорость.

В квантовой механике координата и импульс образуют пару Фурье, и этот же самый компромисс представляет собой принцип неопределенности Гейзенберга. Это не ограничение измерений. Это свойство волн.

Откуда берутся ошибки

Работая с преобразованием Фурье, студенты спотыкаются в трех предсказуемых местах.

Первое место связано с потерей фазы. Поскольку на графике обычно рисуют амплитудный спектр, студенты начинают думать о преобразовании как о списке одних только амплитуд. Затем обратное преобразование выходит неверным, потому что аргументы чисел $\hat{f}(\omega)$ были отброшены. Преобразование является списком комплексных чисел, и половина каждого числа заключена в его аргументе.

Второе место связано с путаницей в знаках и константах. Разные учебники ставят 2π в разные места и выбирают разные знаки в показателе, и все они правы, пока прямое и обратное преобразования согласованы между собой. Выберите одно соглашение, запишите обе формулы вверху страницы и проверьте, что прямое и обратное преобразование простой функции возвращает ее без изменений.

Третье затруднение возникает от ожидания, что конечный сигнал даст чистый спектр. Реальные измерения конечны, поэтому их спектры всегда показывают расплывание из предыдущего раздела плюс рябь от резкого обрыва на краях. Это не ошибка вычислений. Это заявляет о себе компромисс, и для управления им существуют оконные функции.

Где помогает Math Zen

Преобразование Фурье собрано из трех навыков, которые должны работать автоматически, прежде чем интуиция сможет вас нести: тригонометрические тождества, комплексные экспоненты и интегралы от произведений.

Если переход между $\cos(\omega t)$ и комплексными экспонентами требует сознательного усилия или если интегрирование $e^{-at} e^{-i\omega t}$ вручную вызывает сомнения, преобразование навсегда останется черным ящиком. Math Zen тренирует эти составляющие в коротких ежедневных сессиях:

• Тригонометрия: Свободное владение формулами суммы и двойного угла и картиной окружности за ними, в модуле практики тригонометрии.
• Комплексные числа: Показательная форма, формула Эйлера и поворот через умножение, в модуле практики комплексных чисел.
• Преобразования: Та же механика интеграла с ядром, включая преобразования производных, в модуле практики преобразования Лапласа.

Как только все это станет рефлексом, преобразование Фурье перестанет быть формулой и станет тем, чем является: машиной намотки, которая считывает рецепт сигнала.

Главный вывод

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

Формула представляет собой машину намотки: множитель $e^{-i\omega t}$ наматывает сигнал на окружность со скоростью ω, интеграл находит центр масс, а удаленный от начала координат центр масс означает, что эта частота присутствует. Комплексная экспонента ловит сразу все фазы, отрицательные частоты являются зеркальным отражением, которое несет любой действительный сигнал, а обратное преобразование разматывает намотку, вращая в другую сторону.

Поставьте ручку затухания преобразования Лапласа на ноль, и вы получите преобразование Фурье. Используйте Лапласа для переключений и переходных процессов, Фурье для спектров и устойчивых сигналов, и помните, что какой бы из них вы ни взяли, короткий во времени сигнал всегда широк по частоте.