Преобразование ряда фурье. Простыми словами о преобразовании фурье

Одним из мощных средств исследования задач математической физики является метод интегральных преобразований. Пусть функция 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) и

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

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

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

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

Историческая справка

Первым человеком, применившим данный метод, стал французский математик Жан Батист Фурье. Преобразование, названное впоследствии его именем, изначально использовалось для описания механизма теплопроводности. Фурье всю свою сознательную жизнь занимался изучением свойств тепла. Он внес огромный вклад в математическую теорию определения корней алгебраических уравнений. Фурье являлся профессором анализа в Политехнической школе, секретарем Института египтологии, состоял на императорской службе, на которой отличился во время строительства дороги на Турин (под его руководством было осушено более 80 тысяч квадратных километров малярийных болот). Однако вся эта активная деятельность не помешала ученому заниматься математическим анализом. В 1802 году им было выведено уравнение, которое описывает распространение тепла в твердых телах. В 1807 году ученый открыл метод решения данного уравнения, которое и получило название "преобразование Фурье".

Анализ теплопроводности

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

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

Суть анализа

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

Вызов современникам

Алгоритм преобразования Фурье стал вызовом теоретическим основам математики того времени. В начале девятнадцатого века большинство выдающихся ученых, в том числе и Лагранж, Лаплас, Пуассон, Лежандр и Био, не приняли его утверждение о том, что начальное распределение температуры раскладывается на составляющие в виде основной гармоники и более высокочастотные. Однако академия наук не могла проигнорировать результаты, полученные математиком, и удостоила его премии за теорию законов теплопроводности, а также проведение сравнения ее с физическими экспериментами. В подходе Фурье главное возражение вызывал тот факт, что разрывная функция представлена суммой нескольких синусоидальных функций, которые являются непрерывными. Ведь они описывают разрывающиеся прямые и кривые линии. Современники ученого никогда не сталкивались с подобной ситуацией, когда разрывные функции описывались комбинацией непрерывных, таких как квадратичная, линейная, синусоида либо экспонента. В том случае, если математик был прав в своих утверждениях, то сумма бесконечного ряда тригонометрической функции должна сводиться к точной ступенчатой. В то время подобное утверждение казалось абсурдным. Однако, несмотря на сомнения, некоторые исследователи (например Клод Навье, Софи Жермен) расширили сферу исследований и вывели их за пределы анализа распределения тепловой энергии. А математики тем временем продолжали мучиться вопросом о том, может ли сумма нескольких синусоидальных функций сводиться к точному представлению разрывной.

200-летняя история

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

Преобразование Фурье сегодня

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

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

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

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

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

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

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

4. Согласно теореме «свертки», данный процесс превращает сложную операцию в элементарное умножение.

5. Дискретное преобразование Фурье может быть быстро рассчитано на компьютере с использованием «быстрого» метода.

Разновидности преобразования Фурье

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

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

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

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

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

Заключение

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

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

Совокупность операций, позволяющих по заданной функции f(t) находить соответствующую ей спектральную характеристику F(iω ) называется преобразованием Фурье:

Символически формулу (1)будем записывать в виде

Интеграл в правой части (1) как и ранее понимается в смысле главного значения, т.е.

Равенство (1) устанавливает связь между функцией f(t) , аргументом которой служит t , и ей соответствующей комплексной функцией F(iω ), имеющей в качестве аргумента частоту ω .

Формула интеграла Фурье

позволяет от известной функции F(iω ) определить соответствующую ей функцию f(t). На этом основании формулу (3) называют обратным преобразованием Фурье. Символически будем записывать

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

В этом случае f(t) ≡ 0 при t < 0 (1) принимает вид

Преобразование (5) называется прямым односторонним преобразованием Фурье .Обратное преобразование Фурье, соответствующее прямому одностороннему преобразованию, остается двусторонним по переменной ω и дается равенством

При t= 0, значение правой части (6) равно
;

при t < 0 , f(t) ≡ 0

Связь преобразований фурье и лапласа Формула

прямого преобразования Лапласа может рассматриваться как результат определенным образом построенного обобщения одностороннего преобразования Фурье.

Пусть, например, f(t) удовлетворяет условиям Дирихле в интервале 0 ≤ t < ∞ , причем f(t) ≡ 0 при t< 0.

Как известно, преобразование Фурье может быть применено к функциям f(t) , для которых интеграл
существует (условие абсолютной интегрируемости). Этому условию не удовлетворяют многие функции, используемые при анализе процессов в автоматических системах, например 1(t ), Asin(ωt) , Acos(ωt), e αt при α >0, t и др.

Для того чтобы иметь возможность подобную функцию f(t) преобразовать по Фурье, предварительно ее надо умножить на e -ct где вещественное число С>C 0 выбрано таким образом, чтобы интеграл
был бы сходящимся.

Значение С 0 для каждой функции f(t) является вполне определенным. Используя формулу прямого одностороннего преобразования Фурье, будем преобразовывать по Фурье не f(t) , а f(t)e -ct , удовлетворяющую условиям применения этого преобразования.

Введя новую комплексную переменную S=c+jω, получим
.

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

Если F(jω) спектральная х – тика f(t), то функция F(S) комплексной переменой S является спектральной характеристикой затухающей функции времени f(t)e -ct .

Рассмотрим формулу обратного преобразования Фурье:

Заменим в правой и левой частях этого равенства f(t) на f(t)e -ct , получим:

Учитывая, что S=e + jω, dω=dS/j, найдём

Это равенство является формулой обратного преобразования Лапласа, т.е. обратное преобразование Лапласа может рассматриваться как развитие обратного преобразования Фурье.

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

Свойства преобразования Фурье аналогичны свойствам преобразования Лапласа.

Спектральные характеристики некоторых функций

1.Единичная ступенчатая функция. Дельта – функция.

Функция 1(t) вида

называется единичной ступенчатой функцией. Из (1) следует, 1(t) при t=0 имеет разрыв неопределенности первого рода, причем значение функции в точке разрыва не определено. Однако 1(t) при t=0 приписывают вполне определённые значения. Наиболее часто встречаются функции следующего вида:

Выбор того или иного значения единичной функции t=0 связан особенностями решаемой задачи. Например, первое представление удобно в том случае, когда рассматривают функцию 1(t) как предел при λ→∞ последовательности непрерывных функций:

f(t,λ)=1/2+(1/π)arctg λt (3) ,

где λ – параметр и

Последовательность непрерывных функций

при λ→ ∞ также имеет своим пределом первое представление 1(t).

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

Линейность . Рассмотрим функции и
, имеющие спектры
и
:

(12)

Тогда спектр их линейной комбинации будет:

Задержка во времени . Считаем, что известен спектр
сигнала

(14)

Рассчитаем спектр сигнала, сдвинутого во времени:
. Обозначим аргумент функции новой переменной
, тогда
и

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

Изменение масштаба. Считаем, что известен спектр
сигнала
, как через
выражается спектр сигнала
. Вводим новую переменную
, делаем замену переменной интегрирования
.

(16)

Умножение на
. Как и в предыдущем случае, считаем, что известен спектр
сигнала
. Найдем спектр этого сигнала, умноженного на
.

Таким образом, умножение сигнала на
приводит к смещению спектра на.

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

(18)

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

(19)

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

(20)

Преобразование Фурье свёртки двух сигналов даёт произведение спектров этих сигналов.

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

(21)

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

Спектр дискретного сигнала

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

Рис. 1. Дискретизация сигнала.

(22)

Прямоугольные импульсы имеют длительность , период повторения:

(23)

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

(24)

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

(25)

Точно такой же вид имеет разложение в ряд Фурье функции:

(26)

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

(27)

Тогда дискретный сигнал будет иметь вид:

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

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

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

(30)

Вычислим
, используя (25). Так как
периодическая функция, её спектр дискретный.

Таким образом, свёртка (30)

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

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

Если это условие нарушено, то после оцифровки сигнала произойдёт наложение периодически повторяющегося спектра (рис. 2). Получившийся в результате наложения спектр будет соответствовать другому сигналу.

Рис. 2. Перекрывание спектров.

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

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

– Непрерывное преобразование Фурье (в англоязычной литературе 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 операций комплексного умножения и сложения. Таким образом, вычислительная сложность алгоритма дискретного преобразования Фурье является квадратичной, другими словами требуется операций комплексного умножения и сложения.