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

Материал из testwiki
Перейти к навигации Перейти к поиску

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

Определение

В математике преобразование Фурье функции Шаблон:Math, зависящей от одной вещественной переменной, является интегральным и задаётся следующей формулой[1]:

f^(ω)=12πf(x)eixωdx,

где iмнимая единица.

Тогда «формула обращения» (обратное преобразование Фурье) имеет вид:

f(x)=12πf^(ω)eixωdω

Эта формула объясняет физический смысл преобразования Фурье: правая часть — (бесконечная) сумма гармонических колебаний eiωx с частотами ω, амплитудами 12π|f^(ω)| и фазовыми сдвигами argf^(ω) соответственно.

В радиотехнике (обработке сигналов) преобразование Фурье задаётся без множителя 12π[2]:

f^(ω)=f(x)eixωdx.

Тогда «формула обращения» (обратное преобразование Фурье) имеет вид:

f(x)=12πf^(ω)eixωdω

Свойства

Хотя формула, задающая преобразование Фурье, имеет ясный смысл только для функций класса L1(), преобразование Фурье может быть определено и для более широкого класса функций и даже обобщённых функций. Это возможно благодаря ряду свойств преобразования Фурье:

(αf+βg)^=αf^+βg^.

при наличии множителя 12π в преобразовании Фурье:

|f(x)|2dx=|f^(w)|2dω.

при отсутствии множителя 12π в преобразовании Фурье:

|f(x)|2dx=12π|f^(w)|2dω.

Это свойство позволяет по непрерывности распространить определение преобразования Фурье на всё пространство L2().

Равенство Парсеваля будет при этом справедливо для всех fL2().

Шаблон:Якорь

  • Теорема о свёртке: если f,gL1(), тогда

при наличии множителя 12π в преобразовании Фурье:

(fg)^=2πf^g^,

при отсутствии множителя 12π в преобразовании Фурье:

(fg)^=f^g^,

где

fg=f(s)g(ts)ds — свертка функций f и g.

Эта формула может быть распространена и на случай обобщённых функций.

  • Преобразование Фурье и дифференцирование. Если f,fL1(), то
(f)^=iωf^.

Из этой формулы легко выводится формула для n-й производной:

(f(n))^=(iω)nf^.

Формулы верны и в случае обобщённых функций.

  • Преобразование Фурье и сдвиг.
f(xx0)^=eiωx0f^(ω).

Эта и предыдущая формула являются частными случаями теоремы о свёртке, так как сдвиг по аргументу — это свёртка со сдвинутой дельта-функцией δ(xx0), а дифференцирование — свёртка с производной дельта-функции.

  • Преобразование Фурье и растяжение.
f(ax)^=|a|1f^(ω/a).
  • Формула суммирования Пуассона для принятого в данной статье определения:

при наличии множителя 12π в преобразовании Фурье:

k=+f(k)=2πn=+f^(2πn)
k=+f^(k)=2πn=+f(2πn)

при отсутствии множителя 12π в преобразовании Фурье:

k=+f(k)=n=+f^(2πn)
k=+f^(k)=n=+f(2πn)
Данные формулы могут быть получены из Шаблон:Iw, которая задана для другой формы определения преобразования Фурье.
  • Преобразование Фурье обобщённых функций. Преобразование Фурье можно определить для широкого класса обобщённых функций. Определим вначале пространство гладких быстро убывающих функций (пространство Шварца):
S():={φC():n,mxnφ(m)(x)x0}.

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

Теперь определим его двойственное пространство S*(). Это некоторое подпространство в пространстве всех обобщённых функций — так называемые обобщённые функции медленного роста. Теперь для функции fS*() её преобразованием Фурье называется обобщённая функция f^S*(), действующая на основные функции по правилу

f^,φ=f,φ^.

Например, вычислим преобразование Фурье дельта-функции (при наличии множителя 12π в преобразовании Фурье):

δ^,φ=δ,φ^=δ,12πφ(x)eiωxdx=12πφ(x)1dx=12π,φ.

Таким образом, преобразованием Фурье дельта-функции является константа 12π.

Принцип неопределённости

Рассмотрим сигнал x(t), для которого преобразование Фурье имеет вид: F(ω)=x(t)eiωtdt.

Перейдя от частоты ω=2πf к частоте f получим: F(f)=x(t)ei2πftdt.

Чем больше концентрация сигнала x(t) во временной области, тем более размазанным должен быть модуль его преобразования Фурье |F(f)|. В частности, свойство масштабирования преобразования Фурье можно представить так: если сжать функцию в Шаблон:Mvar раз, то её преобразование Фурье растягивается в Шаблон:Mvar раз. Невозможно произвольно сконцентрировать как функцию, так и её преобразование Фурье.

Предположим, что f(t) — квадратично-интегрируемая функция. Тогда норма или энергия сигнала выражается как[3]:

E=f2(t)dt.

Среднее значение для распределения энергии сигнала по времени имеет вид[3]:

t0=tf2(t)dt.

В качестве меры длительности сигнала можно использовать удвоенную величину среднеквадратичной длительности Δt=2Dt, называемую эффективной длительностью сигнала, где

Dt=(tt0)2f2(t)dt.

В терминах теории вероятности Dt — это центральный второй момент функции f(t).

Среднее значение для распределения энергии сигнала в частотной области имеет вид:

F0=f|F(f)|2df=0,

так как подынтегральная функция нечётна.

В качестве меры локализации сигнала в частотной области можно использовать величину Δf=2Df, называемую эффективной шириной полосы частот сигнала, где

Df=(fF0)2|F(f)|2df=f2|F(f)|2df.

В терминах теории вероятности Df — это центральный второй момент функции |F(f)|[3].

Принцип неопределённости гласит, что для дифференцируемых вещественных сигналов x(t) с энергией E, для которых интеграл t0 сходится (то есть t0<) и limt±tf2(t)=0, произведение эффективной длительности сигнала Δt и эффективной ширины полосы частот сигнала Δf ограничено снизу[3]:

ΔtΔfEπ,

Равенство ΔtΔf=Eπ достигается только в случае гауссова импульса f(t)=Cekt2, где k и C некоторые константы (k>0)[3].

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

Более сильным принципом неопределённости является принцип неопределённости Хиршмана, который выражается как:

H(|x|2)+H(|F|2)ln(e2)

где H(p) — дифференциальная энтропия функции плотности вероятности p(z):

H(p)=p(z)ln(p(z))dz,

Равенство достигается для функции Гаусса, как и в предыдущем случае.

Применения

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

Разновидности

Многомерное преобразование

Преобразование Фурье функций, заданных на пространстве n, определяется формулой

f^(ω)=1(2π)n/2nf(x)eixωdx.

Здесь ω и x — векторы пространства n, xω — их скалярное произведение. Обратное преобразование в этом случае задаётся формулой

f(x)=1(2π)n/2nf^(ω)eixωdω.

Эта формула может быть интерпретирована как разложение функции f в линейную комбинацию (суперпозицию) «плоских волн» вида eixω с амплитудами 1(2π)n/2|f^(ω)|, частотами ω и фазовыми сдвигами argf^(ω) соответственно. Как и прежде, в разных источниках определения многомерного преобразования Фурье могут отличаться выбором константы перед интегралом.

Замечание относительно области задания преобразования Фурье и его основные свойства остаются справедливыми и в многомерном случае, со следующими уточнениями:

  • Взятие частных производных под действием преобразования Фурье превращается в умножение на одноимённую координату:
fxk^=iωkf^(ω).
  • Изменяется константа в теореме о свёртке:
(fg)^=(2π)n/2f^g^.
  • Преобразование Фурье и сжатие координат:
(f(x|a|))^=|a|nf^(ω|a|).
(f(Ax))^=|det(A)|1f^((AT)1ω).

Ряды Фурье

Шаблон:Основная статья Непрерывное преобразование само фактически является обобщением более ранней идеи рядов Фурье, которые определены для 2π-периодических функций и представляют собой разложение таких функций в (бесконечную) линейную комбинацию гармонических колебаний с целыми частотами:

f(x)=n=f^neinx.

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

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

f^(ω)=2πn=f^nδ(ωn).

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

Дискретное преобразование

Шаблон:Основная статья Дискретное преобразование Фурье — преобразование конечных последовательностей (комплексных) чисел, которое, как и в непрерывном случае, превращает свёртку в поточечное умножение. Используется в цифровой обработке сигналов и в других ситуациях, где необходимо быстро выполнять свёртку, например, при умножении больших чисел.

Пусть x0,x1,,xn1 — последовательность комплексных чисел. Рассмотрим многочлен f(t)=x0+x1t+x2t2++xn1tn1. Выберем какие-нибудь n точек на комплексной плоскости z0,z1,,zn1. Теперь многочлену f(t) мы можем сопоставить новый набор из n чисел: f0:=f(z0),f1:=f(z1),,fn1:=f(zn1). Заметим, что это преобразование обратимо: для любого набора чисел f0,f1,,fn1 существует единственный многочлен f(t) степени не выше n1 с такими значениями в z0,,zn1 соответственно (см. Интерполяция).

Набор {fk} и называется дискретным преобразованием Фурье исходного набора {xk}. В качестве точек zk обычно выбирают корни n-й степени из единицы:

zk=e2πikn.

Такой выбор продиктован тем, что в этом случае обратное преобразование принимает простую форму, а также тем, что вычисление преобразования Фурье может быть выполнено особенно быстро. Так, в то время как вычисление свёртки двух последовательностей длины n напрямую требует порядка n2 операций, переход к их преобразованию Фурье и обратно по быстрому алгоритму может быть выполнен за O(nlogn) операций. Для преобразований Фурье свёртке соответствует покомпонентное умножение, которое требует лишь порядка n операций.

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

Шаблон:Основная статья

F(t,ω)=f(τ)W(τt)eiωτdτ,

где F(t,ω) даёт распределение частот (вообще говоря, несколько искажённое) части оригинального сигнала f(t) в окрестности момента времени t.

Классическое преобразование Фурье имеет дело со спектром сигнала, взятым во всём диапазоне существования переменной. Нередко интерес представляет только локальное распределение частот, в то время как требуется сохранить изначальную переменную (обычно время). В этом случае используется обобщение преобразования Фурье — так называемое оконное преобразование Фурье. Для начала необходимо выбрать некоторую оконную функцию W, причём эта функция должна иметь хорошо локализованный спектр.

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

Некоторые анализаторы спектра используют быстрое (или кратковременное) оконное преобразование. При нём сигнал заданной длительности разбивается на ряд интервалов с помощью скользящего окна того или иного типа. Это позволяет получать, исследовать и строить в виде спектрограмм динамические спектры и анализировать их поведение во времени. Спектрограмма строится в трёх координатах — частота, время и амплитуда. При этом амплитуда задаётся цветом или оттенком цвета каждого прямоугольника спектрограммы. Подобные анализаторы спектра называют анализаторами спектра реального времени. Основным их производителем является корпорация Keysight Technologies (США), Rohde & Schwarz (Германия), Tektronix (США). Такие анализаторы появились в конце прошлого века и ныне бурно развиваются. Частотный диапазон исследуемых ими сигналов достигает сотен гигагерц.

Указанные методы спектрального анализа реализуются и в системах компьютерной математики, например, Mathcad, Mathematica, Maple и MATLAB.

Другие варианты

Дискретное преобразование Фурье является частным случаем (и иногда применяется для аппроксимации) дискретного во времени преобразования Фурье (DTFT), в котором xk определены на дискретных, но бесконечных областях, и таким образом спектр является непрерывным и периодическим. Дискретное во времени преобразование Фурье является по существу обратным для рядов Фурье.

Эти разновидности преобразования Фурье могут быть обобщены на преобразования Фурье произвольных локально компактных абелевых топологических групп, которые изучаются в гармоническом анализе; они преобразуют группу в её дуальную группу. Эта трактовка также позволяет сформулировать теорему о свёртке, которая устанавливает связь между преобразованиями Фурье и свёртками. См. также дуализм Понтрягина.

Интерпретация в терминах времени и частоты

В терминах обработки сигналов преобразование берёт представление функции сигнала в виде временны́х рядов и отображает его в частотный спектр, где ω — угловая частота. То есть оно превращает функцию времени в функцию частоты; это разложение функции на гармонические составляющие на различных частотах.

Когда функция f является функцией времени и представляет физический сигнал, преобразование имеет стандартную интерпретацию как спектр сигнала. Абсолютная величина получающейся в результате комплексной функции F(ω) пропорциональна амплитудам соответствующих частот ω, в то время как фазовые сдвиги являются аргументами этой комплексной функции.

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

Примеры формул

Следующая таблица содержит список формул для преобразования Фурье. F(ω) и G(ω) обозначают Фурье компоненты функций f(t) и g(t), соответственно. f и g должны быть интегрируемыми функциями или обобщёнными функциями.

Функция Образ с множителем 12π Образ с множителем 1 Примечания
f(x) F(ω)=12πf(x)eixωdx F(ω)=f(x)eixωdx
1 af(t)+bg(t) aF(ω)+bG(ω) aF(ω)+bG(ω) Линейность
2 f(ta) eiωaF(ω) eiωaF(ω) Запаздывание
3 eiatf(t) F(ωa) F(ωa) Частотный сдвиг
4 f(at) |a|1F(ωa) |a|1F(ωa) Если a большое, то f(at) сосредоточена около нуля, и |a|1F(ωa) становится плоским
5 dnf(t)dtn (iω)nF(ω) (iω)nF(ω) Свойство преобразования Фурье от n-й производной
6 tf(τ)dτ π2F(0)δ(ω)+12πiωF(ω) πF(0)δ(ω)+1iωF(ω) Свойство преобразования Фурье от интеграла
7 tnf(t) indnF(ω)dωn indnF(ω)dωn Это обращение правила 5
8 f(t)*g(t) 2πF(ω)G(ω) F(ω)G(ω) Запись f*g означает свёртку f и g: f(t)*g(t)=f(τ)g(tτ)dτ. Это правило — теорема о свёртке
9 f(t)g(t) F(ω)*G(ω)2π F(ω)*G(ω)2π Это обращение 8
10 δ(t) 12π 1 δ(t) означает дельта-функцию Дирака
11 1 2πδ(ω) 2πδ(ω) Обращение 10.
12 tn in2πδ(n)(ω) in2πδ(n)(ω) Здесь n — натуральное число, δ(n)(ω) — n-я производная дельта-функции Дирака. Следствие правил 7 и 11. Использование его вместе с правилом 1 позволяет делать преобразования любых многочленов
13 eiat 2πδ(ωa) 2πδ(ωa) Следствие 3 и 11
14 cos(at) 2πδ(ωa)+δ(ω+a)2 π(δ(ωa)+δ(ω+a)) Следствие 1 и 12 с использованием формулы Эйлера cos(at)=12(eiat+eiat)
15 sin(at) 2πδ(ωa)δ(ω+a)2i iπ(δ(ωa)δ(ω+a)) Также из 1 и 13
16 exp(at2) 12aexp(ω24a) πaexp(ω24a) Показывает, что функция Гаусса exp(t2/2) совпадает со своим изображением
17 W2πsinc(Wt) rect(ω2W) 2πrect(ω2W) Прямоугольная функция — передаточная характеристика идеального фильтра нижних частот, а функция sinc(x) — его импульсная характеристика
18 1t iπ2sgn(ω) iπsgn(ω) Здесь sgn(ω) — функция sgn. Это правило согласуется с 7 и 11
19 1tn iπ2(iω)n1(n1)!sgn(ω) iπ(iω)n1(n1)!sgn(ω) Обобщение 18
20 sgn(t) 2π1iω 2iω Обращение 17
21 θ(t) 12π(1iω+πδ(ω)) 1iω+πδ(ω) Здесь θ(t) — функция Хевисайда. Следует из правил 1 и 20

См. также

Примечания

Шаблон:Примечания

Литература

Ссылки

Шаблон:Интегральное исчисление Шаблон:Интегральные преобразования Шаблон:Методы сжатия Шаблон:Производные буквы F Шаблон:Внешние ссылки