В математике анализ Фурье ( / ˈ f ʊr i eɪ , - i ər / ) [ 1] — это изучение того, как общие функции могут быть представлены или аппроксимированы суммами более простых тригонометрических функций . Анализ Фурье вырос из изучения рядов Фурье и назван в честь Жозефа Фурье , который показал, что представление функции в виде суммы тригонометрических функций значительно упрощает изучение теплопередачи .
Предмет анализа Фурье охватывает широкий спектр математики. В науке и технике процесс разложения функции на колебательные компоненты часто называется анализом Фурье, в то время как операция восстановления функции из этих частей известна как синтез Фурье . Например, определение того, какие частоты компонентов присутствуют в музыкальной ноте, будет включать вычисление преобразования Фурье сэмплированной музыкальной ноты. Затем можно было бы повторно синтезировать тот же звук, включив частотные компоненты, как показано в анализе Фурье. В математике термин анализ Фурье часто относится к изучению обеих операций.
Сам процесс разложения называется преобразованием Фурье . Его выход, преобразование Фурье , часто получает более конкретное название, которое зависит от области определения и других свойств преобразуемой функции. Более того, исходная концепция анализа Фурье со временем была расширена для применения ко все более абстрактным и общим ситуациям, а общая область часто известна как гармонический анализ . Каждое преобразование, используемое для анализа (см. список преобразований, связанных с Фурье ), имеет соответствующее обратное преобразование, которое можно использовать для синтеза.
Для использования анализа Фурье данные должны быть равномерно распределены. Были разработаны различные подходы для анализа неравномерно распределенных данных, в частности, методы спектрального анализа наименьших квадратов (LSSA), которые используют метод наименьших квадратов для подгонки синусоид к выборкам данных, аналогичный анализу Фурье. [2] [3] Анализ Фурье, наиболее используемый спектральный метод в науке, обычно усиливает длиннопериодический шум в записях с большими пробелами; LSSA смягчает такие проблемы. [4]
Анализ Фурье имеет множество научных приложений — в физике , уравнениях с частными производными , теории чисел , комбинаторике , обработке сигналов , цифровой обработке изображений , теории вероятностей , статистике , криминалистике , ценообразовании опционов , криптографии , численном анализе , акустике , океанографии , гидролокации , оптике , дифракции , геометрии , анализе структуры белков и других областях.
Такая широкая применимость обусловлена многими полезными свойствами преобразований :
В криминалистике лабораторные инфракрасные спектрофотометры используют анализ преобразования Фурье для измерения длин волн света, на которых материал будет поглощать в инфракрасном спектре. Метод Фурье используется для декодирования измеренных сигналов и записи данных о длинах волн. И с помощью компьютера эти вычисления Фурье выполняются быстро, так что за считанные секунды компьютерный прибор FT-IR может создать картину поглощения инфракрасного излучения, сравнимую с таковой призматического прибора. [9]
Преобразование Фурье также полезно в качестве компактного представления сигнала. Например, сжатие JPEG использует вариант преобразования Фурье ( дискретное косинусное преобразование ) небольших квадратных фрагментов цифрового изображения. Компоненты Фурье каждого квадрата округляются до более низкой арифметической точности , а слабые компоненты устраняются, так что оставшиеся компоненты могут быть сохранены очень компактно. При реконструкции изображения каждый квадрат изображения собирается заново из сохраненных приближенных преобразованных Фурье компонентов, которые затем подвергаются обратному преобразованию для получения приближения исходного изображения.
При обработке сигналов преобразование Фурье часто берет временной ряд или функцию непрерывного времени и отображает его в частотный спектр . То есть, оно берет функцию из временной области в частотную область; это разложение функции на синусоиды разных частот; в случае ряда Фурье или дискретного преобразования Фурье синусоиды являются гармониками основной частоты анализируемой функции.
Когда функция является функцией времени и представляет собой физический сигнал , преобразование имеет стандартную интерпретацию как частотный спектр сигнала. Величина результирующей комплекснозначной функции на частоте представляет собой амплитуду частотного компонента, начальная фаза которого задается углом (полярные координаты).
Преобразования Фурье не ограничиваются функциями времени и временными частотами. Они в равной степени могут применяться для анализа пространственных частот и, по сути, практически для любой области функций. Это оправдывает их использование в таких различных отраслях, как обработка изображений , теплопроводность и автоматическое управление .
При обработке сигналов, таких как аудио , радиоволны , световые волны, сейсмические волны и даже изображения, анализ Фурье может изолировать узкополосные компоненты сложной формы волны, концентрируя их для более легкого обнаружения или удаления. Большое семейство методов обработки сигналов состоит из преобразования Фурье сигнала, манипулирования преобразованными Фурье данными простым способом и обратного преобразования. [10]
Вот некоторые примеры :
Чаще всего неквалифицированный термин преобразование Фурье относится к преобразованию функций непрерывного действительного аргумента, и оно производит непрерывную функцию частоты, известную как распределение частот . Одна функция преобразуется в другую, и операция обратима. Когда областью определения входной (начальной) функции является время ( ), а областью определения выходной (конечной) функции является обычная частота , преобразование функции на частоте задается комплексным числом :
Оценка этой величины для всех значений дает функцию частотной области . Тогда может быть представлена как рекомбинация комплексных экспонент всех возможных частот :
что является формулой обратного преобразования. Комплексное число передает как амплитуду, так и фазу частоты
Более подробную информацию см. в разделе Преобразование Фурье , в том числе :
Преобразование Фурье периодической функции с периодом становится функцией гребенки Дирака , модулированной последовательностью комплексных коэффициентов :
Обратное преобразование, известное как ряд Фурье , представляет собой представление в виде суммы потенциально бесконечного числа гармонически связанных синусоид или сложных экспоненциальных функций, каждая из которых имеет амплитуду и фазу, заданную одним из коэффициентов :
Любую функцию можно выразить как периодическое суммирование другой функции :
и коэффициенты пропорциональны выборкам на дискретных интервалах :
Обратите внимание, что любое преобразование, имеющее те же дискретные значения выборки, может использоваться в периодическом суммировании. Достаточным условием для восстановления (и, следовательно, ) только из этих выборок (т.е. из ряда Фурье) является то, что ненулевая часть должна быть ограничена известным интервалом длительности , который является частотной областью, двойственной теореме выборки Найквиста–Шеннона .
Более подробную информацию, включая историческое развитие, см. в разделе «Ряды Фурье» .
DTFT является математическим дуалом ряда Фурье во временной области. Таким образом, сходящееся периодическое суммирование в частотной области может быть представлено рядом Фурье, коэффициенты которого являются выборками связанной непрерывной временной функции :
которое известно как DTFT. Таким образом, DTFT последовательности также является преобразованием Фурье модулированной функции гребенки Дирака . [B]
Коэффициенты ряда Фурье (и обратное преобразование) определяются следующим образом :
Параметр соответствует интервалу выборки, и этот ряд Фурье теперь можно распознать как форму формулы суммирования Пуассона . Таким образом, мы имеем важный результат, что когда дискретная последовательность данных пропорциональна отсчетам базовой непрерывной функции, можно наблюдать периодическое суммирование непрерывного преобразования Фурье, Обратите внимание, что любой с теми же дискретными значениями выборки производит то же самое DTFT. Но при определенных идеализированных условиях теоретически можно восстановить и точно. Достаточным условием для идеального восстановления является то, что ненулевая часть должна быть ограничена известным частотным интервалом шириной Когда этот интервал является применимой формулой реконструкции является формула интерполяции Уиттекера–Шеннона . Это краеугольный камень в основе цифровой обработки сигналов .
Еще одна причина, по которой стоит этим заинтересоваться, заключается в том, что этот метод часто дает представление о степени искажения, вызванного процессом выборки.
Приложения DTFT не ограничиваются выборочными функциями. См. Дискретное преобразование Фурье для получения дополнительной информации по этой и другим темам, включая :
Подобно ряду Фурье, ДВПФ периодической последовательности с периодом становится функцией гребенки Дирака, модулированной последовательностью комплексных коэффициентов (см. ДВПФ § Периодические данные ) :
Последовательность обычно известна как ДПФ одного цикла Она также является -периодической, поэтому никогда не требуется вычислять больше, чем коэффициенты. Обратное преобразование, также известное как дискретный ряд Фурье , задается как :
Когда выражается как периодическое суммирование другой функции :
коэффициенты представляют собой выборки с дискретными интервалами :
Наоборот, когда требуется вычислить произвольное количество дискретных выборок одного цикла непрерывного ДВПФ, это можно сделать, вычислив относительно простое ДПФ, как определено выше. В большинстве случаев выбирается равным длине ненулевой части Увеличение, известное как заполнение нулями или интерполяция , приводит к более близко расположенным выборкам одного цикла Уменьшение вызывает перекрытие (добавление) во временной области (аналогично наложению спектров ), что соответствует прореживанию в частотной области. (см. Дискретное преобразование Фурье § L=N×I ) В большинстве случаев, представляющих практический интерес, последовательность представляет собой более длинную последовательность, которая была усечена путем применения оконной функции конечной длины или массива КИХ-фильтров .
ДПФ можно вычислить с помощью алгоритма быстрого преобразования Фурье (БПФ), что делает его практичным и важным преобразованием на компьютерах.
Более подробную информацию см. в разделе Дискретное преобразование Фурье , в том числе :
Для периодических функций и преобразование Фурье, и ДВПФ включают только дискретный набор частотных компонентов (ряды Фурье), и преобразования расходятся на этих частотах. Одна из распространенных практик (не обсуждавшаяся выше) заключается в обработке этого расхождения с помощью дельта-функций Дирака и гребенчатых функций Дирака . Но ту же спектральную информацию можно выделить только из одного цикла периодической функции, поскольку все остальные циклы идентичны. Аналогично функции конечной длительности можно представить в виде ряда Фурье без фактической потери информации, за исключением того, что периодичность обратного преобразования является просто артефактом.
На практике обычно длительность s (•) ограничивается периодом P или N. Но эти формулы не требуют этого условия.
Когда действительная и мнимая части комплексной функции разлагаются на четную и нечетную части , есть четыре компонента, обозначенные ниже индексами RE, RO, IE и IO. И есть взаимно-однозначное отображение между четырьмя компонентами комплексной временной функции и четырьмя компонентами ее комплексного частотного преобразования : [11]
Из этого очевидны различные соотношения, например :
Ранняя форма гармонических рядов восходит к древней вавилонской математике , где они использовались для вычисления эфемерид (таблиц астрономических положений). [12] [13] [14] [15]
Классические греческие концепции деферента и эпицикла в птолемеевской системе астрономии были связаны с рядами Фурье (см. Деферент и эпицикл § Математический формализм ).
В наше время варианты дискретного преобразования Фурье использовались Алексисом Клеро в 1754 году для вычисления орбиты [16], которая была описана как первая формула для ДПФ [17], и в 1759 году Жозефом Луи Лагранжем при вычислении коэффициентов тригонометрического ряда для вибрирующей струны. [17] Технически работа Клеро представляла собой ряд, состоящий только из косинусов (форма дискретного косинусного преобразования ), в то время как работа Лагранжа представляла собой ряд, состоящий только из синусов (форма дискретного синусного преобразования ); истинное косинусное + синусное ДПФ использовалось Гауссом в 1805 году для тригонометрической интерполяции орбит астероидов . [18] Эйлер и Лагранж оба дискретизировали задачу о вибрирующей струне, используя то, что сегодня называется выборками. [17]
Ранним современным развитием анализа Фурье стала работа Лагранжа 1770 года « Réflexions sur la résolution algébrique des équations» , в которой в методе резольвент Лагранжа использовалось комплексное разложение Фурье для изучения решения кубического уравнения : [19] Лагранж преобразовал корни в резольвенты :
где ζ — кубический корень из единицы , что является ДПФ порядка 3.
Ряд авторов, в частности Жан Лерон Д'Аламбер и Карл Фридрих Гаусс , использовали тригонометрические ряды для изучения уравнения теплопроводности [20] , но прорывным достижением стала работа 1807 года «Mémoire sur la propagation de la haleur dans les corps solides» Жозефа Фурье , чьей ключевой идеей было моделирование всех функций тригонометрическими рядами, что ввело ряд Фурье.
Историки расходятся во мнениях относительно того, какую долю заслуги Лагранжа и других в развитии теории Фурье следует отдать Даниилу Бернулли и Леонарду Эйлеру : Даниил Бернулли и Леонард Эйлер ввели тригонометрические представления функций, а Лагранж дал решение волнового уравнения в виде ряда Фурье, так что вклад Фурье в основном заключался в смелом утверждении, что произвольная функция может быть представлена рядом Фурье. [17]
Последующее развитие этой области известно как гармонический анализ и также является ранним примером теории представлений .
Первый алгоритм быстрого преобразования Фурье (БПФ) для ДПФ был открыт около 1805 года Карлом Фридрихом Гауссом при интерполяции измерений орбиты астероидов Юноны и Паллады , хотя этот конкретный алгоритм БПФ чаще приписывают его современным первооткрывателям Кули и Тьюки . [18] [16]
В терминах обработки сигналов функция (времени) — это представление сигнала с идеальным временным разрешением , но без информации о частоте, в то время как преобразование Фурье имеет идеальное частотное разрешение , но без информации о времени.
В качестве альтернатив преобразованию Фурье в частотно-временном анализе используются частотно-временные преобразования для представления сигналов в форме, которая имеет некоторую временную информацию и некоторую частотную информацию – по принципу неопределенности , между ними существует компромисс. Это могут быть обобщения преобразования Фурье, такие как кратковременное преобразование Фурье , преобразование Габора или дробное преобразование Фурье (FRFT), или могут использоваться различные функции для представления сигналов, как в вейвлет-преобразованиях и чирплет-преобразованиях , при этом вейвлет-аналогом (непрерывного) преобразования Фурье является непрерывное вейвлет-преобразование .
Варианты Фурье также могут быть обобщены до преобразований Фурье на произвольных локально компактных абелевых топологических группах , которые изучаются в гармоническом анализе ; там преобразование Фурье переводит функции на группе в функции на двойственной группе. Эта трактовка также позволяет дать общую формулировку теоремы о свертке , которая связывает преобразования Фурье и свертки . См. также двойственность Понтрягина для обобщенных основ преобразования Фурье.
Более конкретно, анализ Фурье можно проводить на смежных классах [21], даже на дискретных смежных классах.