Основные свойства преобразований Фурье

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

a n s n (t) ? a n S n (щ)

2. Свойства четности

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

3. Изменение аргумента функции (сжатие или расширение сигнала) приводит к обратному изменению аргумента ее фурье-образа и обратно пропорциональному изменению его модуля.
4. Теорема запаздывания. Запаздывание (сдвиг, смещение) сигнала по аргументу функции на интервал t o приводит к изменению фазочастотной функции спектра (фазового угла всех гармоник) на величину -щt o без изменения модуля (амплитудной функции) спектра.

5. Преобразование производной (дифференцирование сигнала):

s(t) = d/dt = d/dt =Y(щ) dщ= = jщ Y(щ) exp(jщt) dщ jщ Y(щ).

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

6. Преобразование интеграла сигнала в частотной области при известном спектре сигнала может быть получено из следующих простых соображений. Если имеет место s(t) = d/dt jщY(щ) = S(щ), то должна выполняться и обратная операция: y(t) =s(t) dt Y(щ) = S(щ)/jщ. Отсюда следует:

s(t)dt ? (1/j щ)S(щ).

Оператор интегрирования в частотной области (1/j щ) при щ >1 ослабляет в амплитудном спектре высокие частоты и при щ <1 усиливает низкие. Фазовый спектр сигнала смещается на -90 0 для положительных частот и на 90 0 для отрицательных.

7. Преобразование свертки сигналов y(t) = s(t) * h(t):

Y(щ) =y(t) exp(-jщt) dt =s(ф) h(t- ф) exp(-jщt) dфdt

Y(щ) =s(ф) d ф h(t- ф) exp(-jщt) dt.

По теореме запаздывания:

h(t- ф) exp(-jщt) dt = H(щ) exp(-jщt).

Y(щ) =H(щ) s(ф) exp(-jщ ф) dф= H(щ)·S(щ).

s(t) * h(t)?S(щ)H(щ).

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

8. Преобразование произведения сигналов y(t) = s(t)·h(t):

Y(?) =s(t) h(t) exp(-j?t) dt =s(t) [(1/2?)H(?") exp(j?"t) d?"] dt = (1/2?)s(t)H(?") exp(-j(?-?")t) d?"dt = (1/2?)H(?") d?"s(t) exp(-j(?-?")t) dt = (1/2?)H(?") S(?-?") d?" = (1/2?) H(?) * S(?).

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

9. Умножение сигнала на гармоническую функцию заполняет сигнал гармонической частотой и формирует радиоимпульс.

10. Спектры мощности. Если функция s(t) имеет фурье-образ S(?), то спектральная плотность мощности данной функции определяется выражениями:

w(t) = s(t) s*(t) = |s(t)| 2 |S(?)| 2 = S(?) S*(?) = W(?).

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

11. Равенство Парсеваля. Полная энергия спектра сигнала:

E s =W(f)df=|S(f)| 2 df.

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

|s(t)| 2 dt =|S(f)| 2 df,

т.е. энергия сигнала равна интегралу модуля его частотного спектра - сумме энергий всех частотных составляющих сигнала.

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

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

a n s n (t) Û a n S n (w). (4.21)

Пример суммирования сигналов и его отображения в спектральной области на рис. 4.18.

Рис. 4.18. Сигналы и их спектры. s0(k)=s1(k)+s2(k) Û S1(w)+S2(w) = S0(w)

2. Свойства симметрии преобразования определяются косинусными (четными, действительными) и синусными (нечетными, мнимыми) частями разложения и подобием прямого и обратного преобразований.

На рис. 4.19. приведены примеры, поясняющие свойства четности преобразования. Сигнал s1(k) является четным, s1(k) = s1(-k), и имеет только вещественный четный спектр (мнимая часть спектральной функции представлена нулевыми значениями). Сигнал s2(k) = -s2(-k) нечетный и имеет мнимый нечетный спектр, а нулевыми значениями представлена его действительная часть. Сигнал s3(k) образован суммой сигналов s1(k) и s2(k). Соответственно, спектральная функция сигнала представлена и действительной четной частью (принадлежащей s1(k)), и мнимой нечетной частью (принадлежащей s2(k)). При обратном преобразовании Фурье раздельно действительной и мнимой части спектра S3(w), равно как и любых других комплексных спектров, будут раздельно восстановлены четная и нечетная части исходного сигнала.

Произвольный исходный сигнал может быть задан в одностороннем варианте (0-Т), но четная и нечетная части этого сигнала занимают интервал от –Т до Т, при этом на левой половине числовой оси (от –Т до 0) эти два сигнала компенсируют друг друга, давая нулевые значения.

Рис. 4.19. Свойства четности преобразования

3. Изменение аргумента функции (сжатие или расширение сигнала) приводит к обратному изменению аргумента ее фурье-образа и обратно пропорциональному изменению его модуля. Так, если s(t) Û S(w), то при изменении длительности сигнала с сохранением его формы (растяжении сигнала по временной оси), т.е. для сигнала с новым аргументом s(x) = s(at) при x=at, получаем:

s(at) Ûs(at)exp(-jwt) dt = (1/a)s(x)exp(-jxw/a) dx

s(at) Û (1/a) S(w/a). (4.22")

Выражение (4.22") действительно при а>0. При а<0 происходит зеркальный поворот сигнала относительно вертикальной оси, а замена переменной t=x/a вызывает перестановку пределов интегрирования и, соответственно, изменение знака спектра:

s(at) Û -(1/a) S(w/a). (4.22"")

Обобщенная формула изменения аргумента:

s(at) Û (1/|a|) S(w/a), a ≠ 0 (4.22)

Если под аргументом функции и ее спектра понимать определенные физические единицы, например, время - частота, то отсюда следует: чем короче по своей длительности сигнал, тем шире по частоте его спектр, и наоборот. Это можно наглядно видеть на рис. 4.18. для сигналов s1(k) и s2(k) и их спектров S1(w) и S2(w).

От изменения аргумента функций следует отличать изменение масштаба представления функций. Изменение масштаба аргументов изменяет оцифровку числовых осей отображения сигналов и их спектров, но не изменяет самих сигналов и спектров. Так, при масштабе оси времен t=1 секунда, масштаб оси частот f=1/t=1 герц, а при t=1 мксек f=1/t=1 МГц (t=at, f=1/at, a=10 -6).

4. Теорема запаздывания. Запаздывание (сдвиг, смещение) сигнала по аргументу функции на интервал t o приводит к изменению фазочастотной функции спектра (фазового угла всех гармоник) на величину -wt o . Применяя замену переменной t-t o = x, получаем:

s(t-t o)Ûs(t-t o)exp(-jwt) dt =

S(x)exp(-jwx)exp(-jwt o) dx = S(w)exp(-jwt o). (4.23)

Совершенно очевидно, что амплитуды гармоник сигнала при его сдвиге изменяться не должны. С учетом того, что |exp(-jwt o)|=1, это следует и из (4.23):

|S(w) exp(-jwt o)| = |S(w)|.

Фазовый спектр сдвигается на -wt o с линейной зависимостью от частоты:

S(w) exp(-jwt o)= R(w) expexp(-jwt o)= R(w) exp. (4.24)

Пример двух одинаковых сигналов, сдвинутых относительно друг друга на t o =1, и соответствующих данным сигналам спектров приведен на рис. 4.20.

Рис. 4.20. Изменение спектра сигнала при его сдвиге

Аналогично нетрудно показать, что сдвиг спектра в частотной области на w 0 вызывает умножение сигнала на exp(jw 0 t):

S(w - w 0) « s(t) exp(jw 0 t),

что эквивалентно модуляции сигнала функцией комплексной экспоненты во временной области.

5. Преобразование производной (дифференцирование сигнала):

s(t) = d/dt = d/dt =Y(w) dw =

Jw Y(w) exp(jwt) dw Û jw Y(w). (4.25)

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

Рис. 4.21. Спектры сигнала и его производной

Пример сигнала, его производной и соответствующих им спектров приведен на рис. 4.21. По изменению аргумента спектра (для четного исходного сигнала он был нулевым) можно видеть, что для всех гармоник спектра появляется сдвиг фаз на p/2 (90 0) для положительных частот, и на -p/2 (-90 0) для отрицательных частот.

В общем случае, для кратных производных:

d n /dt n = (jw) n Y(w). (4.26)

При дифференцировании спектра функции соответственно получаем:

d n /dw n = (-jt) n s(t).

6. Преобразование интеграла сигнала в частотной области при известном спектре сигнала может быть получено из следующих простых соображений. Если имеет место

s(t) = d/dt Û jw Y(w) = S(w),

то должна выполняться и обратная операция: y(t) =s(t) dt Û Y(w) = S(w)/jw.

Отсюда следует:

s(t)dt Û (1/jw)S(w). (4.27)

Оператор интегрирования в частотной области (1/jw) при w>1 ослабляет в амплитудном спектре высокие частоты и при w<1 усиливает низкие. Фазовый спектр сигнала смещается на -90 0 для положительных частот и на 90 0 для отрицательных. Пример модуля спектра сигнала и его интегральной функции приведен на рис. 4.22.

Рис. 4.22. Сигналы и амплитудные спектры сигналов

Формула (4.27) справедлива для сигналов с нулевой постоянной составляющей. При интегрировании сигналов с определенным значением постоянной составляющей С=const в правой части выражения (4.27) появляется дополнительное слагаемое преобразования Фурье постоянной составляющей C, которое представляет собой дельта-функцию на нулевой частоте с весовым коэффициентом, равным значению С:

7. Преобразование свертки сигналов y(t) = s(t) * h(t):

Y(w) =y(t) exp(-jwt) dt =s(t) h(t-t) exp(-jwt) dtdt.

Y(w) =s(t) dth(t-t) exp(-jwt) dt.

По теореме запаздывания (4.23):

h(t-t) exp(-jwt) dt = H(w) exp(-jwt).

Y(w) =H(w) s(t) exp(-jwt) dt = H(w)·S(w).

s(t) * h(t) Û S(w) H(w). (4.28)

Пример выполнения свертки в частотной области приведен на рис. 4.23.

Рис. 4.23. Сигналы и амплитудные спектры сигналов

Отметим, что частотное представление H(w) импульсного отклика h(t) линейной системы (или соответствующей линейной операции) имеет смысл частотной передаточной функции системы и позволяет определить сигнал на выходе системы (в частотной форме представления) при задании произвольного сигнала (в частотной форме) на ее входе. По существу, функция H(w) представляет собой распределение по частоте коэффициента пропускания частотных составляющих сигнала с входа на выход системы.

Это положение имеет фундаментальное значение в практике обработки данных.

Любая линейная система обработки данных (информационных сигналов) реализует определенную операцию трансформации сигнала, т.е. выполняет операцию свертки входного сигнала s(t) с оператором системы h(t). С использованием преобразования свертки эта операция может производиться как с динамической, так и с частотной формой представления сигналов. При этом обработка данных, представленных в цифровой форме, производится, как правило, в частотной области, т.к. может быть на несколько порядков выше по производительности, чем во временной области. Она представляет собой последовательность следующих операций.

1. Перевод сигнала в частотную область: s(t) Û S(w).

2. Умножение спектра сигнала на передаточную функцию системы: Y(w) = H(w)·S(w).

Передаточная функция системы определяется аналогичным преобразованием h(t) Û H(w) или задается непосредственно в частотном представлении, что позволяет задавать передаточные функции сколь угодно сложной формы, в том числе с разрывами и скачками, для которых во временной области потребуются операторы h(t) с бесконечной импульсной характеристикой.

3. Перевод спектра обработанного сигнала во временную область: Y(w) Û y(t).

8. Преобразование произведения сигналов y(t) = s(t)·h(t):

Y(w) =s(t) h(t) exp(-jwt) dt =s(t) [(1/2p)H(w") exp(jw"t) dw"] dt =

= (1/2p)s(t)H(w") exp(-j(w-w")t) dw"dt =

(1/2p)H(w") dw"s(t) exp(-j(w-w")t) dt =

= (1/2p)H(w") S(w-w") dw" = (1/2p) H(w) * S(w). (4.29)

Таким образом, произведение функций в координатной форме отображается в частотном представлении сверткой фурье-образов этих функций, с нормировочным множителем (1/2p), учитывающем несимметричность прямого и обратного преобразования Фурье функций s(t) и h(t) при использовании угловых частот.

9. Производная свертки двух функций s"(t) = d/dt.

С использованием выражений (4.26) и (4.28), получаем:

s"(t) = jw = (jw X(w)) Y(w) = X(w) (jw Y(w).

s"(t) = x"(t) * y(t) = x(t) * y"(t).

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

10. Спектры мощности. Временная функция мощности сигнала в общей форме определяется выражением:

w(t) = s(t) s * (t) = |s(t)| 2 .

Спектральная плотность мощности, соответственно, равна преобразованию Фурье произведения s(t)·s * (t), которое отобразится в спектральном представлении сверткой Фурье-образов этих функций:

W(f) = S(f) * S * (f) =S(f) S * (f-v) dv. (4.30)

Но для всех текущих значений частоты f интеграл в правой части этого выражения равен произведению S(f)·S * (f), так как для всех значений сдвига v ≠ 0 в силу ортогональности гармоник S(f) и S * (f-v) значения их произведения равны нулю. Отсюда:

W(f) = S(f) * S * (f) = |S(f)| 2 . (4.31)

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

Для функций мощности взаимодействия сигналов в частотной области соответственно имеем частотные спектры мощности взаимодействия сигналов:

W xy (f) = X(f) Y*(f),

W yx (f) = Y(f) X*(f),

W xy (f) = W* yx (f).

Функции мощности взаимодействия сигналов комплексные, даже если обе функции x(t) и y(t) вещественны, при этом Re - четная функция, а Im - нечетная. Отсюда полная энергия взаимодействия сигналов при интегрировании функций мощности взаимодействия определяется только реальной частью спектра:

X(f) Y*(f) df.

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

áx(t),y(t)ñ = áX(f),Y(f)ñ, ||x(t)|| 2 = ||X(f)|| 2 .

Не следует забывать, что при представлении спектров в круговых частотах (по w) в правой части приведенных равенств должен стоять множитель 1/2p.

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

Линейность преобразования Фурье.

Это важнейшее свойство формулируется так: если имеется некоторая совокупность сигналов причем то взвешенная сумма сигналов преобразуется по Фурье следующим образом:

Здесь - произвольные числовые коэффициенты.

Для доказательства формулы (2.26) следует подставить сумму сигналов в преобразование Фурье (2.16).

Свойства вещественной и мнимой частей спектральной плотности.

Пусть - сигнал, принимающий вещественные значения. Его спектральная плотность в общем случае является комплексной:

Подставам это выражение в формулу обратного преобразования фурье (2.18):

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

Это возможно лишь в том случае, если вещественная часть спектральной плотности сигнала есть четная, а мнимая часть - нечетная функция частоты:

Спектральная плотность сигнала, смещенного во времени.

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

Доказательство очень простое. Действительно,

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

Зависимость спектральной плотности сигнала от выбора масштаба измерения времени.

Предположим, что исходный сигнал подвергнут изменению масштаба времени. Это означает, что роль времени t играет новая независимая переменная (k - некоторое вещественное число). Если то происходит «сжатие» исходного сигнала; если же то сигнал «растягивается» во времени.

Оказывается, что если то

Действительно,

откуда следует формула (2.29).

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

К рассматриваемому здесь вопросу близко примыкает Следующая задача.

Дан импульс отличный от нуля на отрезке и характеризуемый спектральной плотностью Требуется иайти спектральную плотность «обращенного во времени» сигнала который представляет собой «зеркальную копию» исходного импульсного колебания. Поскольку очевидно, что то

Выполнив замену переменной находим, что

Спектральная плотность производной и неопределенного интеграла.

Пусть сигнал s(t) и его спектральная плотность заданы. Будем изучать новый сигнал и Поставим цель найти его спектральную плотность - .

По определению,

Преобразование Фурье - линейная операция, значит, равенство (2.31) справедливо и по отношению к спектральным плотностям. Учитывая (2.28), получаем

Представляя экспоненциальную функцию рядом Тейлора: подставляя этот ряд в (2.32) и ограничиваясь первыми двумя членами, находим

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

Формула (2.33) обобщается на случай спектра производной порядка. Легко доказать, что если , то

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

Рассмотренная функция является первообразной (неопределенным интегралом) по отношению к функции Из (2.33) формально следует, что спектр первообразной

Таким образом, множитель служит оператором интегрирования в частотной области.

Спектральная плотность сигнала на выходе интегратора.

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

Здесь - фиксированный параметр.

Определенный интеграл, входящий в (2.36), равен, очевидно, разности двух значений первообразной сигнала одно из которых вычисляется при аргументе t, а другое - при аргументе . Используя соотношения (2.28) и (2.35), получаем формулу связи между спектральными плотностями сигналов на входе и выходе:

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

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

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

– Непрерывное преобразование Фурье (в англоязычной литературе Continue Time Fourier Transform – CTFT или, сокращенно, FT );

– Дискретное преобразование Фурье (в англоязычной литературе Discrete Fourier Transform – DFT );

– Быстрое преобразование Фурье (в англоязычной литературе Fast Fourier transform – FFT ).

Непрерывное преобразование Фурье

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

Существует несколько видов записи непрерывного преобразования Фурье, отличающихся друг от друга значением коэффициента перед интегралом (две формы записи):

или

где и - Фурье-образ функцииили частотный спектр функции ;

- круговая частота.

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

Если функция нечетная на интервале (−∞, +∞), то преобразование Фурье может быть представлено через синус-функцию:

Если функция четная на интервале (−∞, +∞), то преобразование Фурье может быть представлено через косинус-функцию:

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

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

или

где - Фурье-образ функцииили частотный спектр функции ;

- круговая частота.

Если функция нечетная на интервале (−∞, +∞), то обратное преобразование Фурье может быть представлено через синус-функцию:

Если функция четная на интервале (−∞, +∞), то обратное преобразование Фурье может быть представлено через косинус-функцию:

В качестве примера, рассмотрим следующую функцию . График исследуемой экспоненциальной функции представлен ниже.

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

В результате получили зависимость изменения исследуемой экспоненциальной функции на частотном интервале (см. ниже).

Непрерывное преобразование Фурье используют, как правило, в теории при рассмотрении сигналов, которые изменяются в соответствии с заданными функциями, но на практике обычно имеют дело с результатами измерений, которые представляют собой дискретные данные. Результаты измерений фиксируются через равные промежутки времени с определённой частотой дискретизации, например, 16000 Гц или 22000 Гц. Однако в общем случае дискретные отсчёты могут идти неравномерно, но это усложняет математический аппарат анализа, поэтому на практике обычно не применяется.

Существует важная теорема Котельникова (в иностранной литературе встречается название «теорема Найквиста-Шеннона», «теорема отсчетов»), которая гласит, что аналоговый периодический сигнал, имеющий конечный (ограниченный по ширине) спектр (0…fmax), может быть однозначно восстановлен без искажений и потерь по своим дискретным отсчётам, взятым с частотой, большей или равной удвоенной верхней частоте спектра - частота дискретизации (fдискр >= 2*fmax). Другими словами, при частоте дискретизации 1000 Гц из аналогового периодического сигнала можно восстановить сигнал с частотой до 500 Гц. Следует отметить, что дискретизация функции по времени приводит к периодизации ее спектра, а дискретизация спектра по частоте приводит к периодизации функции.

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

Прямое дискретное преобразование Фурье ставит в соответствие временной функции , которая определена N-точками измерений на заданном временном интервале, другую функцию , которая определена на частотном интервале. Следует отметить, что функция на временном интервале задается с помощью N-отсчетов, а функция на частотном интервале задается с помощью K-кратного спектра.

k ˗ индекс частоты.

Частота k-го сигнала определяется по выражению

где T - период времени, в течение которого брались входные данные.

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

Из последних выражений видно, что преобразование раскладывает сигнал на синусоидальные составляющие (которые называются гармониками) с частотами от одного колебания за период до N колебаний за период.

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

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

Амплитуда постоянной составляющей является средним значением функции за выбранный промежуток времени и определяется следующим образом:

Амплитуды и фазы частотных составляющих сигнала определяются по следующим соотношениям:

Полученные значения амплитуды и фазы называют полярным представлением (polar notation). Результирующий вектор сигнала будет определяться следующим образом:

Рассмотрим алгоритм преобразования дискретно заданной функции на заданном интервале (на заданном периоде) с количеством исходных точек

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

В результате преобразования получаем вещественное и мнимое значение функции , которая определена на частотном диапазоне.

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

N ˗ количество значений сигнала, измеренных за период, а также кратность частотного спектра;

k ˗ индекс частоты.

Как уже было сказано, дискретное преобразование Фурье N-точкам дискретного сигнала ставит в соответствие N-комплексных спектральных отсчетов сигнала . Для вычисления одного спектрального отсчета требуется N операций комплексного умножения и сложения. Таким образом, вычислительная сложность алгоритма дискретного преобразования Фурье является квадратичной, другими словами требуется операций комплексного умножения и сложения.

Одним из мощных средств исследования задач математической физики является метод интегральных преобразований. Пусть функция f(x) задана на интервале (а, 6), конечном или бесконечном. Интегральным преобразованием функции f(x) называется функция где К(х, ш) - фиксированная для данного преобразования функция, называемая ядром преобразования (предполагается, что интеграл (*) существуете собственном или несобственном смысле). §1. Интеграл Фурье Всякая функция f(x), которая на отрезке [-f, I] удовлетворяет условиям разложимости в ряд Фурье, может быть на этом отрезке представлена тригонометрическим рядом Коэффициенты а*, и 6„ ряда (1) определяются по формулам Эйлера-Фурье: ПРЕОБРАЗОВАНИЕ ФУРЬЕ Интеграл Фурье Комплексная форма интеграла Преобразование Фурье Косинус и синус преобразования Амплитудный и фазовый спектры Свойства Приложения Ряд в правой части равенства (1) можно записать в иной форме. С этой целью внесем в него из формул (2) значения коэффициентов а» и оп, подведем под знаки интегралов cos ^ х и sin х (что возможно, поскольку переменной интегрирования является т) О) и используем формулу для косинуса разности. Будем иметь Если функция/(ж) первоначально была определена на интервале числовой оси, большем, чем отрезок [-1,1] (например, на всей оси), то разложение (3) воспроизведет значения этой функции только на отрезке [-1,1] и продолжит се на всю числовую ось как периодическую функцию с периодом 21 (рис. 1). Поэтому, если функция f(x) (вообще говоря, непериодическая) определена на всей числовой оси, в формуле (3) можно попытаться перейти к пределу при I +оо. При этом естественно потребовать выполнения следующих условий: 1. f(x) удовлетворяет условиям разложимости в ряд Фурье на любом конечном отрезке оси Ох\ 2. функция f(x) абсолютно интегрируема на всей числовой оси, При выполнении условия 2 первое слагаемое правой части равенства (3) при I -* +оо стремится к нулю. В самом деле, Попытаемся установить, во что перейдет в пределе при I +оо сумма в правой, части (3). Положим так, что Тогда сумма в правой части (3) примет вид В силу абсолютной сходимости интеграла эта сумма при больших I мало отличается от выражения которое напоминает интегральную сумму для функции переменного £ составленную для интервала (0, +оо) изменения Поэтому естественно ожидать, что при сумма (5) перейдет в интеграл Сдругой стороны, при фиксировано) из формулы (3) вытекает, что и мы получаем равенство Достаточное условие справедливости формулы (7) выражается следующей теоремой. Теорема 1. Если функция f(x) абсолютно интегрируема на всей числовой оси и имеет вместе со своей производной конечное число точек разрыва первого рода на любом отрезке [а, 6], то справедливо равенство При этом во всякой точке xq, являющейся точкой разрыва 1-го рода функции /(ж), значение интеграла в правой части (7) равно Формулу (7) называют интегральной формулой Фурье, а стоящий в ее правой части интеграл - интегралом Фурье. Если воспользоваться формулой дня косинуса разности, то формулу (7) можно записать в виде Функции а(£), Ь(£) являются аналогами соответствующих коэффициентов Фурье ап и Ьп 2тг-периодической функции, но последние определены для дискретных значений п, вто время как а(0> НО определеныдля непрерывных значений £ G (-оо, +оо). Комплексная форма интеграла Фурье Предполагая /(х) абсолютно интегрируемой на всей оси Ох, рассмотрим интеграл Этот интеграл равномерно сходится для, так как и потому представляет собой непрерывную и, очевидно, нечетную функцию от Но тогда С другой стороны, интеграл есть четная функция переменной так что Поэтому интегральную формулу Фурье можно записать так: Умножим равенство на мнимую единицу i и прибавим к равенству (10). Получим откуда, в силу формулы Эйлера будем иметь Это - комплексная форма интеграла Фурье. Здесь внешнее интегрирование по £ понимается в смысле главного значения по Коши: §2. Преобразование Фурье. Косинус- и синус-преобразования Фурье Пусть функция f(x) является кусочно-гладкой на любом конечном отрезке оси Ох и абсолютно интегрируема на всей оси. Определение. Функция откуда, в силу формулы Эйлера, будем иметь называется преобразованием Фурье функции /(г) (спектральной функцией). Это - интегральное преобразование функции /(г) на интервале (-оо,+оо) с ядром Используя интегральную формулу Фурье получаем Это так называемое обратное преобразование Фурье, дающее переход от F(£) к /(х). Иногда прямое преобразование Фурье задают так: Тогда обратное преобразование Фурье определится формулой Преобразование Фурье функции /(ж) определяют также следующим образом: ПРЕОБРАЗОВАНИЕ ФУРЬЕ Интеграл Фурье Комплексная форма интеграла Преобразование Фурье Косинус и синус преобразования Амплитудный и фазовый спектры Свойства Приложения Тогда, в свою очередь, При этом положение множителя ^ достаточно произвольно: он может входить либо в формулу (1"), либо в формулу (2"). Пример 1. Найти преобразование Фурье функции -4 Имеем Это равенство допуска ет дифференцирование по £ под знаком интеграла (получающийся после дифференцирования интеграл равномерно сходится, когда { принадлежит любому конечному отрезку): Интегрируя по частям, будем иметь Внеинтегральное слагаемое обращается в нуль, и мы получаем откуда (С - постоянная интегрирования). Полагая в (4) £ = 0, найдем С = F(0). В силу (3) имеем Известно, что В частности, для) получаем, что Пример 2 (разред кокдемсетора через сопропиление). Рассмотрим функцию 4 Для спектрам ыюй функции F(£) получаем Отсюда (рис.2). Условие абсолютной интегри-руемости функции f(x) на всей числовой оси является весьма жестким. Оно исключает, например, такие элементарные функции, как) = cos ж, f(x) = е1, для которых преобразования Фурье (в рассматриваемой здесь классической форме) не существует. Фурье-образ имеют только те функции, которые достаточно быстро стремятся к нулю при |х| -+ +оо (как в примерах 1 и 2). 2.1. Косинус- и синус-преобразования Фурье Используя формулу косинуса, разности, перепишем интегральную формулу Фурье в следующем виде: Пусть f(x) - четная функция. Тогда так что изравснства (5) имеем В случае нечетной f(x) аналогично получаем Если f(x) задана лишь на (0, -foo), то формула (6) продолжает f(x) на всю ось Ох четным образом, а формула (7) - нечетным. (7) Определение. Функция называется косинус-преобразованием Фурье функции f(x). Из (6) следует, что для четной функции f(x) Это означает, что f(x), в свою очередь, является косинус-преобразованием для Fc(£). Иными словами, функции / и Fc являются взаимными косинус-преобразованиями. Определение. Функция называется синус-преобразованием Фурье функции f(x). Из (7) получаем, что для нечетной функции f(x) т.е. f и Fs являются взаимными синус-преобразованиями. Пример 3 (прамоугольный импульс}. Пусть f(t) - четная функция, определенная следующим образом: (рис. 3). Воспользуемся полученным результатом для вычисления интеграла В силу формулы (9) имеем Рис.3 0 0 В точке t = 0 функция f(t) непрерывна и равна единице. Поэтому из (12") получим 2.2. Амплитудный и фазовый спектры интеграла Фурье Пусть периодическая с периодом 2т функция /(х) разлагается в ряд Фурье Это равенство можно записать в виде где - амплитуда колебания с частотой п, - фаза. На этом пути мы приходим к понятиям амплитудного и фазового спектров периодической функции. Для непериодической функции f{x), заданной на (-оо, +оо), при определенньк условиях оказывается возможным представить ее интегралом Фурье осуществляющим разложение этой функции по всем частотам (разложение по непрерывному спектру частот). Определение. Спектральной функцией, или спектральной плотностью интеграла Фурье, называется выражение (прямое преобразование Фурье функции f называется амплитудным спектром, а функция Ф«) = -агgSfc) - фазовым спектром функции /(«). Амплитудный спектр.А(£) служит мерой вклада частоты £ в функцию /(ж). Пример 4. Найти амплитудный и фазовый спектры функции 4 Находим спектральную функцию Отсюда Графики этих функций изображены на рис. 4. §3. Свойства преобразования Фурье 1. Линейность. Если и G(0 - преобразования Фурье функций /(х) и д(х) соответственно, то при любых постоянных а и р преобразованием Фурье функции a f{x) + р д(х) будет функция a Пользуясь свойством линейности интеграла, имеем Таким образом, преобразование Фурье есть линейный оператор. Обозначая его через будем писать. Если F(£) есть преобразование Фурье абсолютно интегрируемой на всей числовой оси функции /(ж), то F(() ограничена при всех. Пусть функция f(x) абсолютно интегрируема на всей оси - преобразование Фурье функции f(x). Тогда 3«fltsJ. Пусть f(x) - функция, допуска кнцэя преобразование Фурье, Л - дойств ительяов число. Фуниция fh(x) = f{z-h) называется сдвигом фунждии f{x). Пользуясь определен нем преобразования Фурье, показать, что Задача. Пусть функция f(z) имеет преобразование Фурье F(0> h - действительное число. Показать, что 3. Преобразование Фурье и ооерэции дифференцирования. Пусть абсолютно интегрируемая функция f(x) имеет производную f"(x), также абсолютно интегрируемую на всей оси Ох, так что /(я) стремится к нулю при |ж| -» +оо. Считая f"(x) гладкой функцией, запишем Интегрируя по частям, будем иметь Внеинтегральноеслагаемое обращается в нуль (так как, и мы получаем Таким образом, дифференцированию функции /(х) отвечает умножение ее образа Фурье ^П/] на множитель Если функция f(x) имеет глад*«е абсолютно интефируемые производные до порядка m включительно и все они, как и сама функция f(x), стремятся к нулю при то, интегрируя по частям нужное число раз, получим Преобразование Фурье очень полезно именно потому, что оно заменяет операцию дифференцирования операцией умножения на величину и тем самым упрощает задачуинтегрирования некоторых видов дифференциальных уравнений. Так как преобразование Фурье абсолютно интегрируемой функции f^k\x) есть ограниченная функция от (свойство 2), то из соотношения (2) получаем для следующую оценку: ПРЕОБРАЗОВАНИЕ ФУРЬЕ Интеграл Фурье Комплексная форма интеграла Преобразование Фурье Косинус и синус преобразования Амплитудный и фазовый спектры Свойства Приложения Из этой оценки следует: чем больше функция f(x) имеет абсолютно интегрируемых производных, тем быстрее ее преобразование Фурье стремится к нулю при. Замечание. Условие является достаточно естественным, поскольку обычная 1еория интегралов Фурье имеет дело с процессами, которые в том или ином смысле имеют начало и коней, но не продолжаются неограниченно с примерно одинаковой интенсивностью. 4. Связь между скоростью убывания функции f(x) при |z| -» -f оо и гладкостью ее преобразования Фурм. Предположим, что не только /(х), но и ее произведение xf(x) является абсолютно интегрируемой функцйей на всей оси Ох. Тогда преобразование Фурье) будет дифференцируемой функцией. Действительно, формальное дифференцирование по параметру £ подынтегральной функции приводит к интегралу который является абсолютно и равномерно сходящимся относительно параметра Следовательно, дифференцирование возможно, и Таким образом, т. е. операция умножения f(x) на аргумент х переходит после преобразования Фурье в операцию t щ. Если вместе с функцией f(x) абсолютно интегрируемыми на всей оси Ох являются функции, то процесс дифференцирования можно продолжить. Получим, что функция имеет производные до порядка m включительно, причем Таким образом, чем быстрее функция f(x) убывает при тем более гладкой получается функция Теорема 2 (о сверле). Пусть- преобразования Фурье функций /,(ж) и f2(x) соответственно. Тогда причем двойной интеграл в правой части сходится абсолютно. Положим - х. Тогда будем иметь или, меняя порядок интегрирования, Функция называется сверткой функций и обозначается символом Формула (1) может быть теперь записана так: Отсюда видно, что преобразование Фурье свертки функций f\(x) и f2(x) равно умноженному на у/2ж произведению преобразований Фурье свертываемых функций, Замечание. Нетрудно установить следующие свойства свертки: 1) линейность: 2) коммутативность: §4. Приложения преобразования Фурье 1. Пусть Р(^) - линейный дифференциальный оператор порядка m с постоянными коэффициентами, Используя формулу для преобразования Фурье производных функции у(х), находим " Рассмотрим дифференциальное уравнение где Р - введенный выше дифференциальный оператор. Предположим, что искомое решение у(х) имеет преобразование Фурье у (О. а функция f(x) имеет преобразование /(£) Применяя преобразование Фурье к уравнению (1), получим вместо дифференциального алгебраическое уравнение на оси относительно откуда так что формально где символ обозначает обратное преобразование Фурье. Основное ограничение применимости этого метода связано со следующим фактом. Решение обыкновенного дифференциального уравнения с постоянными коэффициентами содержит функции вида еЛ*, eaz cos fix, еах sin рх. Они не являются абсолютно интегрируемыми на оси -оо < х < 4-оо, и преобразование Фурье для них не определено, так что, строго говоря, применятьданный метод нельзя. Это ограничение можно обойти, если ввести в рассмотрение так называемые обобщенные функции. Однако в ряде случаев преобразование Фурье все же применимо в своей классической форме. Пример. Найти решение а = а(х, t) уравнения (а = const), при начальных условиях Это - задача о свободных колебаниях бесконечной однородной струны, когда задано начальное отклонение <р(х) точек сгруны, а начальные скорости отсутствуют. 4 Поскольку пространственная переменная х изменяется в пределах от -оо до +оо, подвергнем уравнение и начальные условия преобразованию Фурье по переменной х. Будем предполагать, что 1) функции и(х, t) и