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

Бабалао Ф.Р.

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

01.03.01 Математика

Математика

Дипломы

Вуз: Белгородский государственный университет - национальный исследовательский университет (НИУ «БелГУ»)

ID: 5c1a4fea7966e104f6f851a8

UUID: eb35d130-e5c4-0136-fbf0-525400005860

Язык: Русский

Опубликовано: больше 5 лет назад

Просмотры: 17

Бабалао Ф.Р.

Источник: Федеральное государственное автономное образовательное учреждение высшего образования «Белгородский государственный национальный исследовательский университет»

Комментировать 0

Рецензировать 0

Скачать - 778,6 КБ

Поделиться работой

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ «БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ» (НИУ «БелГУ») ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК КАФЕДРА ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ДИСКРЕТНОЕ ПРЕОБРАЗОВАНИЕ ФУРЬЕ И ЕГО ПРИМЕНЕНИЯ Выпускная квалификационная работа обучающегося по направлению подготовки 01.03.01 Математика очной формы обучения, группы 07001413 Бабалао Фолакэ Рашель Научный руководитель д. ф. - м. н., профессор, Васильев В.Б. Белгород 2018

Оглавление Введение ................................................................................................................... 2 Глава 1 ...................................................................................................................... 5 1.1Определение и основные свойства интегрального преобразования Фурье ................................................................................................................................... 5 1.2 Свертка функций и преобразование Фурье…………………………………13 1.3 Примеры преобразований Фурье……………………………………………14 Глава 2. ................................................................................................................... 23 2.1. Применение преобразования Фурье для решения задач математической физики .................................................................................................................... 24 2.2. Применение преобразования Фурье для решения интегральных уравнений ............................................................................................................... 29 2.3. Преобразование Фурье и ряды Фурье, дискретное преобразование Фурье ................................................................................................................................. 30 Заключение ............................................................................................................ 45 Список литературы ............................................................................................... 46 Приложение ........................................................................................................... 49

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

интерпретировать экспериментальные наблюдения в астрономии, медицине и химии. Французский математик Жан Батист Джозеф Фурье, чье имя называлось преобразованием, первым человеком сказал об этом методе в мире. В 1789 году он получил уравнение, описывающее распространение тепла в твердом теле. В 1807 году Фурье изобрёл и метод решения этого уравнения: преобразование Фурье. Преобразование Фурье используется во многих областях науки - в области физики, теории чисел, комбинаторики, теории вероятностей, статистики, акустики, океанологии, оптики, геометрии и многих других. Благодаря широкому применению метода Фурье и аналогичных аналитических методов мы можем теперь правильно повторить то, что сказал Лорд Кельвин в 1867 году: «Теорема Фурье - это не только один из самых элегантных результатов современного анализа, но он также дает нам незаменимый инструмент в изучении самых сложных проблем современной физики.» Цель исследования является изучением преобразования Фурье, его свойства и теорем с целю определения применения этого преобразование в области математики и физики. Задачи нашего исследования есть: 1. Изучить научную литературу, касающуюся темы исследования. 2. Определить преобразование Фурье и ряды Фурье. 3. Определить что такой дискретное преобразование Фурье, его основные свойства а также проводить примеры преобразований Фурье 4. Составить примеры применения преобразования Фурье в разных областях математики и физики

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

Глава 1 1.1 Определение и основные свойства интегрального преобразования Фурье Преобразование Фурье (символ ℱ) является операцией, которая связывает функцию действительной переменной с другой функцией действительной переменной. Эта новая функция описывает коэффициенты («амплитуды») при разложении исходной функции на элементарные компоненты - гармонические колебания с разными частотами (так же, как музыкальный аккорд может быть выражен как сумма музыкальных звуков, которые его составляют). Преобразование Фурье функции f вещественной переменной является интегральным и задаётся следующей формулой: 𝑓̂(𝜔) = 1 √2𝜋 ∞ ∫ 𝑓(𝑥)𝑒 −𝑖𝑥𝜔 𝑑𝑥 −∞ Определение 1. Функция ∞ 1 ℱ[𝑓](𝜉) = ∫ 𝑓(𝑥)𝑒 −𝑖𝜉𝑥 𝑑𝑥 2𝜋 (1.1.1) −∞ Называется преобразованием Фурье функции 𝑓: ℝ → ℂ. Интеграл здесь понимается в смысле главного значения ∞ 𝐴 ∫ 𝑓(𝑥)𝑒 −𝑖𝜉𝑥 𝑑𝑥 = lim ∫ 𝑓(𝑥)𝑒 −𝑖𝜉𝑥 𝑑𝑥 𝐴→+∞ −∞ −𝐴 и считается что он существует. Если 𝑓: ℝ → ℂ – абсолютно интегрируемая на ℝ функция, то, поскольку |𝑓(𝑥)𝑒 −𝑖𝜉𝑥 | = |𝑓(𝑥)| при 𝑥, 𝜉 ∈ ℝ, для любой такой функции имеет смысл преобразование Фурье (1.1.1), причем интеграл (1.1.1) сходится абсолютно и равномерно по 𝜉 на всей прямой ℝ.

Определение 2. Если 𝑐(𝜉) = ℱ[𝑓](𝜉) – преобразование Фурье функции 𝑓: ℝ → ℂ, то сопоставляемый 𝑓 интеграл ∞ 𝑓(𝑥)~ ∫ 𝑐(𝜉)𝑒 𝑖𝑥𝜉 𝑑𝜉 (1.1.2) −∞ понимаемый в смысле главного значения, называется интегралом Фурье функции 𝑓. Определение 3. Понимаемые в смысле главного значения интегралы ∞ 1 ℱ𝑐 [𝑓](𝜉) = ∫ 𝑓(𝑥) cos 𝜉𝑥 𝑑𝑥, 𝜋 (1.1.3) −∞ ∞ 1 ℱ𝑠 [𝑓](𝜉) = ∫ 𝑓(𝑥) sin 𝜉𝑥 𝑑𝑥, 𝜋 (1.1.4) −∞ Называются соответственно косинус- и синус-преобразованиями Фурье функции 𝑓. Полагая 𝑐(𝜉) = ℱ[𝑓](𝜉), 𝑎(𝜉) = ℱ𝑐 [𝑓](𝜉), 𝑏(𝜉) = ℱ𝑠 [𝑓](𝜉), получаем отчасти уже знакомое нам по рядам Фурье соотношение 1 𝑐(𝜉) = (𝑎(𝜉) − 𝑙𝑏(𝜉)). 2 (1.1.5) Как видно из соотношений (3), (4), 𝑎(−𝜉) = 𝑎(𝜉), 𝑏(−𝜉) = −𝑏(𝜉). (1.1.6) Формулы (1.1.5), (1.1.6) показывают, что преобразования Фурье вполне определяются на всей прямой ℝ, если они известны лишь для неотрицательных значений аргумента. Различные источники могут давать определения, которые отличаются от вышеприведенного выбором коэффициента перед интегралом, а также знак

«-» в экспоненты. Но все свойства будут одинаковыми, хотя форма некоторых формул может измениться. Прежде чем развивать теорию дальше, давайте вычислим три основных примера преобразований Фурье. Рис. a. График функции Xa Пример1. Пусть Xa-функция, изображенная на рисунке а: 1, если |𝑥| < 𝑎 𝜒𝑎 (𝑥) = { 0, в противном случае Тогда 𝑎 𝜒̂ 𝑎 (𝜉) = ∫ 𝑒 −𝑎 −𝑖𝜉𝑥 𝑒 −𝑖𝑎𝜉 − 𝑒 𝑖𝑎𝜉 𝑠𝑖𝑛𝑎𝜉 𝑑𝑥 = =2 −𝑖𝜉 𝜉 Пример 2. Пусть функция 𝑓(𝑥) = 𝑒 −𝑎𝑥 2 /2 где a>0 Мы наблюдаем, что f удовлетворяет дифференциальное уравнение 𝑓 ′ (𝑥) + 𝑎𝑥𝑓(𝑥) = 0 Если применить преобразование Фурье к этому уравнению, получаем

′ 𝑖𝜉𝑓̂(𝜉) + 𝑖𝑎[𝑓̂] (𝜉) = 0 или, ′ [𝑓̂] (𝜉) + 𝑎−1 𝜉𝑓̂(𝜉) = 0 Это дифференциальное уравнение для 𝑓̂ легко решается: ′ 𝜉 𝜉2 [𝑓̂] (𝜉) 2 ̂ = − ⇒ 𝑙𝑜𝑔𝑓 (𝜉) = − + 𝑙𝑜𝑔𝐶 ⇒ 𝑓̂(𝜉) = 𝐶𝑒 −𝜉 /2𝑎 𝑎 2𝑎 𝑓̂(𝜉) Для оценки константы C мы устанавливаем ξ = 0 𝐶 = 𝑓(0) = ∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑒 −𝑎𝑥 2 /2 2 2𝜋 2 𝑑𝑥 = √ ∫ 𝑒 −𝑦 𝑑𝑦 = √ 𝑎 𝑎 Поэтому, 𝑎𝑥 2 − ℱ [𝑒 2 ] (Аккуратный вывод этого 2𝜋 2 = √ 𝑒 −𝜉 /2𝑎 𝑎 результата с использованием контурных интегралов набросан в примере 1) Пример 3. Пусть функция 𝑓(𝑥) = (𝑥 2 + 𝑎2 )−1 Где a > 0 Мы должны посчитать 𝑓̂ здесь по контурной интеграции, еще один вывод. Если 𝜉 < 0, 𝑒 −𝑖𝜉𝑧 -ограниченная аналитическая функция z в верхней полуплоскости. Поэтому, применяя теорема вычетов на контуре на рис.b и пусть N →∞, получим 𝑓̂(𝜉) = ∫ 𝑒 −𝑖𝜉𝑥 𝑒 −𝑖𝜉𝑧 𝑒 𝑎𝜉 𝜋 𝑎𝜉 𝑑𝑥 = 2𝜋𝑖𝑅𝑒𝑠𝑧=𝑖 2 = 2𝜋𝑖 = 𝑒 (𝜉 < 0) 𝑥 2 + 𝑎2 (𝑧 + 𝑎2 ) 2𝑖𝑎 𝑎

Аналогично, если 𝜉 > 0, мы можем интегрировать вокруг нижней полуплоскости, чтобы получить 𝑓̂(𝜉) = 𝜋 −𝑎𝜉 𝑒 𝑎 Конечно, для 𝜉 = 0 𝑓(0) = ∫ 𝑓(𝑥)𝑑𝑥 = 𝜋/𝑎 Заключения: ℱ[(𝑥 2 + 𝑎2 )−1 ] = 𝜋 −𝑎|𝜉| 𝑒 𝑎 Рис.b Контур примера 3 Мораль здесь, как и в теории рядов Фурье, что чем более гладкий f, чем быстрее 𝑓̂ распадается на безграничности, и наоборот. Приведенные выше примеры иллюстрируют этот принцип. Функция примера 1 исчезает из конечного интервала, но не является непрерывной: ее преобразование Фурье аналитическое, но уменьшается только в виде 1/𝜉 до бесконечности. Функция примера 3 гладкая, но медленно убывается: ее преобразование Фурье быстро убывает, но не дифференцируется при ξ = 0. Функция примера 2 есть одновременно гладкость и распад, она по существу, его собственное преобразование Фурье. Кроме того, существуют разнообразные обобщения данного понятия

Свойства Хотя формула преобразования Фурье имеет четкое значение только для функций класса 𝐿1 (ℝ), преобразование Фурье также может быть определено для более широкого класса функций и даже обобщенных функций. Это возможно благодаря определенному числу свойств преобразования Фурье: • Преобразование Фурье является линейным оператором: + 𝛽𝑔) = 𝛼𝑓̂ + 𝛽𝑔̂ (𝛼𝑓̂ • Справедливо равенство Парсеваля: если 𝑓 ∈ 𝐿1 (ℝ) ∩ 𝐿2 (ℝ), то преобразование Фурье сохраняет L2-норму: 2 ∞ ∞ ∫−∞|𝑓(𝑥)|2 𝑑𝑥 = ∫−∞|𝑓̂(𝜔)| 𝑑𝜔 Это свойство позволяет по непрерывности распространить определение преобразования Фурье на всё пространство 𝐿2 (ℝ). Равенство Парсеваля будет при этом справедливо для всех 𝑓 ∈ 𝐿2 (ℝ). • Формула обращения: 𝑓(𝑥) = ∞ 𝑓̂(𝜔)𝑒 𝑖𝑥𝜔 𝑑𝜔 ∫ √2𝜋 −∞ 1 справедлива, если интеграл в правой части имеет смысл. В частности, это верно, если функция f достаточно гладкая. Если 𝑓 ∈ 𝐿2 (ℝ), то формула также верна, так как равенство Парсеваля позволяет дать смысл правой части интеграла с предельным переходом. Эта формула объясняет физический смысл преобразования Фурье: правая часть — (бесконечная) сумма гармонических колебаний 𝑒 𝑖𝑥𝜔 с частотами 𝜔, амплитудами • 1 √2𝜋 𝑓̂(𝜔) и фазовыми сдвигами 𝑓̂(𝜔) соответственно. Теорема о свёртке: если 𝑓, 𝑔 ∈ 𝐿1 (ℝ), тогда ∗ 𝑔) = √2𝜋𝑓̂𝑔̂, где (𝑓̂

∞ (𝑓 ∗)(𝑡) = ∫−∞ 𝑓(𝑡 − 𝑠)𝑔(𝑠)𝑑𝑠 Эта формула может быть распространена и на случай обобщённых функций. • Преобразование Фурье и дифференцирование. Если 𝑓, 𝑓′ ∈ 𝐿1 (ℝ), то (𝑓̂′ ) = 𝑖𝜔𝑓 Из этой формулы легко выводится формула для n-й производной: (𝑛) ) = (𝑖𝜔)𝑛 𝑓̂ (𝑓̂ Формулы верны и в случае обобщённых функций. • Преобразование Фурье и сдвиг. 𝑓(𝑥 − 𝑥0 ) = 𝑒 −𝑖𝜔𝑥0 𝑓̂(𝜔) Эта и предыдущая формула являются частными случаями теоремы о свёртке, так как сдвиг по аргументу — это свёртка со сдвинутой дельта- функцией 𝛿(𝑥 − 𝑥0 ), а дифференцирование — свёртка с производной дельтафункции. • Преобразование Фурье и растяжение. ̂ = |𝑎|−1 𝑓(𝜔) 𝑓(𝑎𝑥) 𝑎 • Преобразование Фурье обобщённых функций. Преобразование Фурье может быть определено для широкого класса обобщенных функций. Сначала определим пространство гладких функций быстрого распада (пространство Шварца): 𝑆(ℝ) ≔ {𝜑 ∈ 𝐶 ∞ (ℝ): ∀𝑛, 𝑚 ∈ ℕ 𝑥 𝑛 𝜑 (𝑚) (𝑥) → 0 при 𝑥 → +∞}

Ключевое свойство этого пространства состоит в том, что оно является инвариантным подпространством относительно преобразования Фурье. 𝑆 ∗ (ℝ). Это некоторое Теперь определим его двойственное пространство подпространство называемые в пространстве обобщённые функции всех обобщённых медленного функций — роста. Теперь так для функции 𝑓 ∈ 𝑆 ∗ (ℝ). её преобразованием Фурье называется обобщённая функция 𝑓̂ ∈ 𝑆 ∗ (ℝ). , действующая на основные функции по правилу 〈𝑓̂, 𝜑〉 = 〈𝑓, 𝜑̂〉 Например, вычислим преобразование Фурье дельта-функции: 〈𝛿̂ , 𝜑〉 = 〈𝛿, 𝜑̂〉 = 〈𝛿, 1 √2𝜋 ∞ ∫ 𝜑(𝑥)𝑒 −𝑖𝜔𝑥 𝑑𝑥 〉 = −∞ 1 ∞ 1 ∫ 𝜑(𝑥) ∙ 1𝑑𝑥 = 〈 , 𝜑〉 √2𝜋 √2𝜋 −∞ Таким образом, преобразованием Фурье дельта-функции является константа 1 √2𝜋 .

Таблица. Основные преобразования Фурье 1.2 Свертка функций и преобразование Фурье Свертка (взаимная свертка) – интегральная операция получения новой функции по двум исходным. ∞ 𝑆𝑣1,2 (𝑥. 𝑦) = 𝑠1 ⊗ 𝑠2 (𝑥, 𝑦) = ∬ 𝑠1 (𝑢, 𝑣)𝑠2 (𝑥 − 𝑢, 𝑦 − 𝑣)𝑑𝑢𝑑𝑣 −∞ где s1(x,y),s2(x,y)- исходные функции; u,v – переменные интегрирования; x,y – аргументы взаимной свертки, характеризующие некоторый сдвиг. Сделаем замену переменной 𝑥 − 𝑢 = 𝑢′ , 𝑦 − 𝑉 = 𝑉′ , и поменяем порядок интегрирования. Тогда выражение приводится к виду:

∞ ∞ 𝑆𝑣12 (𝑣𝑥 , 𝑣𝑦 ) = ∫−∞ ∫−∞ 𝑠1 (𝑢, 𝑉)𝑒 −𝑖2𝜋(𝑣𝑥 𝑢+𝑣𝑦 𝑉) 𝑑𝑢𝑑𝑉 × ∞ ∞ ∫−∞ ∫−∞ 𝑠2 (𝑢′ , 𝑉 ′ )𝑒 −𝑖2𝜋(𝑣𝑥 𝑢′+𝑣𝑦 𝑉′) 𝑑𝑢′ 𝑑𝑉 ′ = 𝑠̃1 (𝑣𝑥 , 𝑣𝑦 )𝑠̃2 (𝑣𝑥 , 𝑣𝑦 ) , т.е. Фурье-образ взаимной свертки равен произведению Фурье-образов свертываемых функций. Отсюда для свертки двух одинаковых функций (автосвертки) имеем: 𝑆𝑣𝑠 (𝑣𝑥 , 𝑣𝑦 ) = 𝐴𝑆𝑣 (𝑣𝑥 , 𝑣𝑦 ) = [𝑠̃ (𝑣𝑥 , 𝑣𝑦 )] 2 1.3 Примеры преобразований Фурье Пример 1. Найти преобразование Фурье функции 0, 𝑓(𝑥) = [ 1, 0, −∞ < 𝑥 < 1 1<𝑥<2 2<𝑥 Заданная функция абсолютно интегрируема на ℝ, действительно, ∞ 2 ∫ |𝑓(𝑥)| 𝑑𝑥 = ∫ 𝑑𝑥 = 1 < ∞ −∞ 1 Так как ∞ ∫ 𝑓(𝑥)𝑒 −∞ 2 −𝑖𝜉𝑥 𝑑𝑥 = ∫ 𝑒 1 −𝑖𝜉𝑥 𝑖 ∙ (𝑒 −2𝑖𝜉 − 𝑒 −𝑖𝜉 ) 𝑑𝑥 = , 𝜉 То 𝑖 ∙ (𝑒 −2𝑖𝜉 − 𝑒 −𝑖𝜉 ) ℱ[𝑓](𝜉) = . 2𝜋𝜉 Пример 2. Найти косинус - и синус - преобразования Фурье функции

0, 𝑓(𝑥) = [ 1, 0, −∞ < 𝑥 < 1 1<𝑥<2 2<𝑥 Как показано в примере 1, заданная функция абсолютно интегрируема на ℝ. Найдем ее косинус - преобразование Фурье по формуле (3): ∞ 2 1 1 sin(2 ∙ 𝜉) − sin(𝜉) ℱ𝑐 [𝑓](𝜉) = ∫ 𝑓(𝑥) ∙ cos(𝜉 ∙ 𝑥) 𝑑𝑥 = ∫ cos(𝜉 ∙ 𝑥) 𝑑𝑥 = . 𝜋 𝜋 𝜋∙𝜉 −∞ 1 Аналогично, нетрудно найти синус – преобразование Фурье функции f(x) по формуле (4): ∞ 2 1 1 cos(𝜉) − cos(2 ∙ 𝜉) ℱ𝑐 [𝑓](𝜉) = ∫ 𝑓(𝑥) ∙ sin(𝜉 ∙ 𝑥) 𝑑𝑥 = ∫ sin(𝜉 ∙ 𝑥) 𝑑𝑥 = . 𝜋 𝜋 𝜋∙𝜉 −∞ 1 Используя примеры 1 и 2, нетрудно непосредственной подстановкой убедиться, что для f(x) выполняется соотношение (5). Если функция 𝑓 вещественнозначна, то из формул (5), (6) в этом случае следует ̅̅̅̅̅̅, 𝑐(−𝜉) = 𝑐(𝜉) (7) Поскольку в этом случае 𝑎(𝜉) и 𝑏(𝜉) – вещественные функции на R, что видно из их определений (3), (4). Впрочем, равенство (7) при условии ̅̅̅̅̅̅ 𝑓(𝑥) = 𝑓(𝑥) получается и непосредственно из определения (1) преобразования Фурье, если учесть, что знак сопряжения можно вносить под знак интеграла. Последние наблюдение позволяет заключить, что для любой функции 𝑓: 𝑅 → ℂ справедливо равенство ̅̅̅̅̅̅̅̅̅̅, ℱ[𝑓](−𝜉) = ℱ[𝑓](𝜉) (8)

Полезно также заметить, что если 𝑓 – вещественная и четная функция, т.е. ̅̅̅̅̅̅ 𝑓(𝑥) = 𝑓(𝑥) = 𝑓(−𝑥), то ̅̅̅̅̅̅̅̅̅̅̅ ℱ𝑐 [𝑓](𝜉) = ℱ𝑐 [𝑓](𝜉), ℱ𝑠 [𝑓](𝜉) ≡ 0, ̅̅̅̅̅̅̅̅̅̅ = ℱ[𝑓](𝜉) = ℱ[𝑓](−𝜉); ℱ[𝑓](𝜉) (9) если 𝑓 – вещественная и нечетная функция, т.е. ̅̅̅̅̅̅ 𝑓(𝑥) = 𝑓(𝑥) = −𝑓(−𝑥), то ̅̅̅̅̅̅̅̅̅̅̅ ℱ𝑐 [𝑓](𝜉) ≡ 0, ℱ 𝑠 [𝑓](𝜉) = ℱ[𝑓](𝜉) ̅̅̅̅̅̅̅̅̅̅ = −ℱ[𝑓](𝜉) = ℱ[𝑓](−𝜉) ℱ[𝑓](𝜉) (10) а если 𝑓 – чисто мнимая функция, т.е. . ̅̅̅̅̅̅ 𝑓(𝑥) = −𝑓(𝑥), то ̅̅̅̅̅̅̅̅̅̅. ℱ[𝑓](−𝜉) = −ℱ[𝑓](𝜉) (11) Заметим что если 𝑓 – вещественнозначная функция, то интеграл Фурье можно записать также в виде ∞ ∞ ∫ √𝑎2 (𝜉) + 𝑏 2 (𝜉)cos(𝑥𝜉 + 𝜑(𝜉)) 𝑑𝜉 = 2 ∫ |𝑐(𝜉)|cos(𝑥𝜉 + 𝜑(𝜉)) 𝑑𝜉, 0 0 где 𝜑(𝜉) = − arctg Пример 3. Найдем 𝑏(𝜉) = arg 𝑐 (𝜉). 𝑎(𝜉) преобразование Фурье функции 𝑓(𝑡) = (считая 𝑓(0) = 𝑎 ∈ ℝ) 𝐴 𝐴 −𝐴 −𝐴 1 sin 𝑎𝑡 −𝑖𝛼𝑡 1 sin 𝑎𝑡 cos 𝛼𝑡 ℱ[𝑓](𝛼) = lim ∫ 𝑒 𝑑𝑡 = lim ∫ 𝑑𝑡 = 𝐴→+∞ 2𝜋 𝐴→+∞ 2𝜋 𝑡 𝑡 +∞ +∞ 2 sin 𝑎𝑡 cos 𝛼𝑡 1 sin(𝑎 + 𝛼)𝑡 sin(𝑎 − 𝛼)𝑡 = ∫ 𝑑𝑡 = ∫ ( + ) 𝑑𝑡 = 2𝜋 𝑡 2𝜋 𝑡 𝑡 −∞ 0 sin 𝑎𝑡 𝑡

∞ 1 1 sin 𝑢 , если |𝛼| ≤ |𝑎|, (𝑠𝑔𝑛(𝑎 + 𝛼) + 𝑠𝑔𝑛(𝑎 − 𝛼)] ∫ 𝑑𝑢 = {2 2𝜋 𝑢 0, если |𝛼| > |𝑎|, 0 поскольку нам известно значение интеграла Дирихле ∞ ∫ 0 sin 𝑢 𝜋 𝑑𝑢 = . 𝑢 2 (12) Рассмотренная в примере функция 𝑓 не является абсолютно интегрируемой на ℝ и её преобразование Фурье имеет разрывы. О том, что преобразование Фурье абсолютно интегрируемых функций не имеет разрывов, говорит следующая Пример 4. Найдем преобразование Фурье функции 𝑓(𝑡) = 𝑒 −𝑡 +∞ 2 ⁄2 : 𝐴 ℱ[𝑓](𝛼) = ∫ 𝑒 −𝑡 2 ⁄2 𝑒 −𝑖𝛼𝑡 𝑑𝑡 = ∫ 𝑒 −𝑡 −∞ 2 ⁄2 cos 𝛼𝑡 𝑑𝑡. −𝐴 Дифференцируя последний интеграл по параметру 𝛼 и интегрируя затем по частям, находим, что 𝑑ℱ[𝑓] 𝛼 (𝛼) + ℱ[𝑓](𝛼) = 0, 𝑑𝛼 2 или 𝑑 𝛼 ln ℱ[𝑓](𝛼) = − . 𝑑𝛼 2 Значит, ℱ[𝑓](𝛼) = 𝑐𝑒 −𝛼 2 /2 , где 𝑐 – постоянная, которую, пользуясь интегралом Эйлера-Пуассона находим из соотношения ∞ 𝑐 = ℱ[𝑓](0) = ∫ 𝑒 −𝑖 2 ⁄2 𝑑𝑡 = √2𝜋. −∞ Итак, мы нашли, что ℱ𝑐 [𝑓](𝛼) = √2𝜋𝑒 −𝛼 ℱ𝑐 [𝑓](𝛼) = √2𝜋𝑒 −𝛼 2 ⁄2 , а ℱ𝑠 [𝑓](𝛼) ≡ 0. 2 ⁄2 , и одновременно показали, что

Пример 5. Представить функцию 𝑓(𝑥) интегралом Фурье, если −1, 1, 𝑓(𝑥) = [ 0, −1 < 𝑥 < 0 0<𝑥<1 |𝑥| > 1 Данная функция 𝑓(𝑥) является нечетной и непрерывной на ℝ, кроме точек 𝑥 = −1, 𝑥 = 0, 𝑥 = 1. В силу нечетности и вещественности функции 𝑓(𝑥) имеем: ℱ𝑐 [𝑓](𝜉) = 0 , и из равенств (5) и (10) следует, что +∞ 𝑖 𝑖 𝑐(𝜉) = − ∙ 𝑏(𝜉) = − ∙ ∫ 𝑓(𝑦) ∙ sin(𝜉 ∙ 𝑦) 𝑑𝑦 = 2 2𝜋 −∞ 1 𝑖 =− ∙ 2 ∙ ∫ sin(𝜉 ∙ 𝑦) 𝑑𝑦 2𝜋 = 0 𝑖 ∙ cos(𝜉 ∙ 𝑦) 1 𝑖 ∙ (cos(𝜉) − 1) |0 = . 𝜋∙𝜉 𝜋∙𝜉 В точках непрерывности функции 𝑓 имеем: +∞ 𝑓(𝑥) = ∫ 𝑐(𝜉) ∙ 𝑒 𝑖𝜉𝑥 𝑑𝜉 = −∞ +∞ 1 𝑖 ∙ (cos(𝜉) − 1) = ∫ ∙ (cos(𝜉 ∙ 𝑥)) + 𝑖 ∙ sin(𝜉 ∙ 𝑥))𝑑𝜉 = 𝜋 𝜉 −∞ = 1 ∙ 𝜋 +∞ +∞ (1 − cos(𝜉)) (cos(𝜉) − 1) ∙{ ∫ ∙ sin(𝜉 ∙ 𝑥) 𝑑𝜉 + 𝑖 ∙ ∫ ∙ 𝜉 𝜉 −∞ −∞ ∙ cos(𝜉 ∙ 𝑥) 𝑑𝜉} . Но функция 𝑔(𝜉) = (cos(𝜉)−1) 𝜉 ∙ cos(𝜉 ∙ 𝑥) − нечетная, поэтому

+∞ ∫ −∞ (cos(𝜉) − 1) ∙ cos(𝜉 ∙ 𝑥) 𝑑𝜉 = 0 , 𝜉 так как интеграл вычисляется в смысле главного значения. Функция ℎ(𝜉) = (1−cos(𝜉)) 𝜉 ∙ sin(𝜉 ∙ 𝑥) − четная, поэтому +∞ (1 − cos(𝜉)) 2 𝑓(𝑥) = ∙ ∫ ∙ sin(𝜉 ∙ 𝑥) 𝑑𝜉 , 𝜋 𝜉 0 если 𝑥 ≠ 0, 𝑥 ≠ ±1. При 𝑥 = 1 должно выполняться равенство +∞ (1 − cos(𝜉)) 𝑓(1 − 0) + 𝑓(1 + 0) 1 2 = = ∙∫ ∙ sin(𝜉) 𝑑𝜉. 2 2 𝜋 𝜉 0 Полагая 𝜉 = 2𝑡, отсюда находим +∞ sin3 𝑡 cos 𝑡 𝜋 ∫ 𝑑𝑡 = . 𝑡 16 0 Если в последнем выражении для 𝑓(𝑥) положить 𝑥 = ½, то +∞ (1 − cos(𝜉)) 1 2 𝜉 𝑓( ) = 1 = ∙∫ ∙ sin ( ) 𝑑𝜉 . 2 𝜋 𝜉 2 0 Полагая здесь 𝜉 = 2𝑡, найдем +∞ ∫ 0 sin3 𝑡 𝜋 𝑑𝑡 = . 𝑡 4 Если вещественнозначная функция 𝑓(𝑥) кусочно непрерывна на любом отрезке действительной прямой, абсолютно интегрируема на 𝑅 и имеет конечные односторонние производные в каждой точке 𝑥 ∈ 𝑅, то в точках непрерывности 𝑓 она представляется как интеграл от Фурье

+∞ +∞ 1 𝑓(𝑥) = ∫ 𝑑𝑦 ∫ 𝑓(𝑡) cos 𝑦(𝑥 − 𝑡) 𝑑𝑡, 𝜋 0 (13.1) −∞ а в точках разрыва функции 𝑓 левую часть равенства (13.1) следует заменить на 𝑓(𝑥 + 0) + 𝑓(𝑥 − 0) . 2 Если непрерывная, абсолютно интегрируемая на 𝑅 функция 𝑓 имеет в каждой точке 𝑥 ∈ 𝑅 конечные односторонние производные, то в случае, когда это функция является четной, справедливо равенство +∞ 𝑓(𝑥) = ∫ 𝑎(𝑦) cos 𝑥𝑦 𝑑𝑦, (13.2) 0 где +∞ 2 𝑎(𝑦) = ∫ 𝑓(𝑡) cos 𝑦𝑡 𝑑𝑡, 𝜋 (13.2′) 0 а в случае, когда 𝑓 - нечетная функция, выполняется равенство +∞ 𝑓(𝑥) = ∫ 𝑏(𝑦) sin 𝑥𝑦 𝑑𝑦, (13.3) 0 где +∞ 2 𝑏(𝑦) = ∫ 𝑓(𝑡) sin 𝑦𝑡 𝑑𝑡, 𝜋 (13.3′) 0 Пример 5’. Представить функцию 𝑓(𝑥) интегралом Фурье, если:

𝑓(𝑥) = 𝑒 −𝑎|𝑥| , 𝑎 > 0; Так как 𝑓 - непрерывная на 𝑅 четная функция, то, используя формулы (13.2), (13.2’), имеем +∞ 2 2𝛼 𝑎(𝑦) = ∫ 𝑒 −𝑎𝑡 cos 𝑦𝑡 𝑑𝑡 = , 𝜋 𝜋(𝛼 2 + 𝑦 2 ) 0 +∞ 𝑒 −𝛼|𝑥| 2𝛼 cos 𝑦𝑥 = ∫ 2 𝑑𝑦, 𝜋 𝛼 + 𝑦2 𝑥 ∈ 𝑅. 0 Обозначим символом 𝐹[𝜑] понимаемый в смысле главного значения интеграл ∞ 𝐹[𝜑](𝑥) = ∫ 𝜑(𝜉)𝑒 𝑖𝑥𝜉 𝑑𝜉. (14) −∞ Пример 6. Пусть 𝑎 > 0 и 𝑒 −𝑎𝑥 при 𝑥 > 0, 𝑓(𝑥) = { 0 при 𝑥 ≤ 0; тогда +∞ 1 1 1 ℱ[𝑓](𝜉) = ∫ 𝑒 −𝑎𝑥 𝑒 −𝑖𝑥𝜉 𝑑𝜉 = ∙ . 2𝜋 2𝜋 𝑎 + 𝑖𝜉 0 Заметим, что если 𝑓− (𝑥) = 𝑓(−𝑥), то при любой функции 𝑓 ∞ ∞ −∞ −∞ 1 1 ℱ[𝑓− ](𝜉) = ∫ 𝑓(−𝑥)𝑒 −𝑖𝑥𝜉 𝑑𝑥 = ∫ 𝑓(𝑥)𝑒 𝑖𝑥𝜉 𝑑𝑥 = ℱ[𝑓](−𝜉). 2𝜋 2𝜋 Возьмем теперь функцию 𝑒 −𝑎|𝑥| = 𝜑(𝑥) = 𝑓(𝑥) + 𝑓(−𝑥). Тогда ℱ[𝜑](𝜉) = ℱ[𝑓](𝜉) + ℱ[𝑓](−𝜉) = 1 𝑎 . 𝜋 𝑎2 + 𝜉 2

Если же взять функцию 𝜓(𝑥) = 𝑓(𝑥) − 𝑓(−𝑥), являющуюся нечетным продолжением функции 𝑒 −𝑎𝑥 , 𝑥 > 0, на всю числовую ось, то ℱ[𝜓](𝜉) = ℱ[𝑓](𝜉) − ℱ[𝑓](−𝜉) = − 𝑖 𝜉 𝜋 𝑎2 + 𝜉 2 Используя теорему 1, получаем, что 𝑒 −𝑎𝑥 , если 𝑥 > 0, 1 𝑒 𝑖𝑥𝜉 1 ∫ 𝑑𝜉 = { , если 𝑥 = 0, 2𝜋 𝑎 + 𝑖𝜉 2 −∞ 0, если 𝑥 < 0, +∞ +∞ 1 𝑎𝑒 𝑖𝑥𝜉 ∫ 2 𝑑𝜉 = 𝑒 −𝑎|𝑥| , 2 𝜋 𝑎 +𝜉 −∞ +∞ 𝑒 −𝑎𝑥 , если 𝑥 > 0, 𝑖 𝜉𝑒 𝑖𝑥𝜉 ∫ 2 𝑑𝜉 = { 0, если 𝑥 = 0, 𝜋 𝑎 + 𝜉2 −𝑒 𝑎𝑥 , если 𝑥 < 0. −∞ Все интегралы здесь понимаются в смысле главного значения Отделяя в двух последних интегралах действительные и мнимые части, находим интегралы Лапласа +∞ ∫ 0 +∞ ∫ 0 cos 𝑥𝜉 𝜋 −𝑎|𝑥| 𝑑𝜉 = 𝑒 , 𝑎2 + 𝜉 2 2𝑎 𝜉 sin 𝑥𝜉 𝜋 −𝑎|𝑥| 𝑑𝜉 = 𝑒 sgn 𝑥. 𝑎2 + 𝜉 2 2𝑎

Глава 2. Преобразование Фурье - полезный инструмент для анализа различных проблем математики и физических наук. Основная часть этих приложений является следующим основным фактом. Предположим, что L является линейным оператором для функций на R, коммутирующих с переводами: то есть, если 𝐿[𝑓(𝑥)] = 𝑔(𝑥), то 𝐿[𝑓(𝑥 + 𝑠)] = 𝑔(𝑥 + 𝑠) для любого s принадлежащей R. Тогда любая экспоненциальная функция eax (𝑎 ∈ ∁), которая принадлежит к области L есть собственной функцией L. Доказательство этого очень простое: Пусть 𝑓(𝑥) = 𝑒 𝑎𝑥 и 𝑔 = 𝐿[𝑓]. Тогда для любого 𝑠 ∈ 𝑅 𝑔(𝑥 + 𝑠) = 𝐿[𝑒 𝑎(𝑥+𝑠) ] = 𝐿[𝑒 𝑎𝑠 𝑒 𝑎𝑥 ] = 𝑒 𝑎𝑠 𝐿[𝑒 𝑎𝑥 ] = 𝑒 𝑎𝑠 𝑔(𝑥) Установляясь x=0 мы находим ,что 𝑔(𝑠) = 𝑔(0)𝑒 𝑎𝑠 для всех 𝑠 ∈ 𝑅. Другими словами, g = Cf где C = g(0) . Таким образом 𝐿[𝑓] = 𝐶𝑓 . Предположим, в частности, что домен L включает в себя все мнимые экспоненциальные 𝑒 𝑖𝜉𝑥 и пусть h() - собственное значение для 𝑒 𝑖𝜉𝑥 ; таким образом, 𝐿[𝑒 𝑖𝜉𝑥 ] = ℎ(𝜉)𝑒 𝑖𝜉𝑥 . Если L удовлетворяет некоторым очень мягким условиям непрерывности, можно прочитать действие L на более или менее произвольная функция f из формулы обращения Фурье. Действительно, эта формула выражает f как непрерывную суперпозицию экспонент 𝑒 𝑖𝜉𝑥 𝑓(𝑥) = и тогда, 1 ∫ 𝑓(𝜉)𝑒 𝑖𝜉𝑥 𝑑𝜉 2𝜋

𝐿[𝑓](𝑥) = 1 1 ∫ 𝑓̂(𝜉)𝐿[𝑒 𝑖𝜉𝑥 ]𝑑𝜉 = ∫ 𝑓̂(𝜉)ℎ(𝜉)𝑒 𝑖𝜉𝑥 𝑑𝜉 2𝜋 2𝜋 (условия непрерывности на L используются здесь, чтобы оправдать рассмотрение интеграла как если бы это была конечная сумма). Таким образом, в терминах преобразования Фурье 𝑓̂ , действия L сводится к простой алгебраической операции умножения на функцию h, ̂ = ℎ𝑓̂, поэтому проход от f к 𝑓̂ может упростить анализ L больш. (𝐿[𝑓]) Альтернативно, если h -преобразование Фурье функции H, мы можем выразить L как свертку: Lf = f*H. Для того, чтобы эти расчеты работали в соответствии с теорией, f и H должно быть L1 или L2 функции. Однако, можно продлить домен Фурье включить гораздо более общие виды функций, и настоящее обсуждение затем расширяется к более общей ситуации. 2.1. Применение преобразования Фурье для решения задач математической физики Приведем пример решения задачи с помощью преобразования Фурье. Задача. Построить функцию Грина для уравнения, описывающего распространение тепла в бесконечном однородном стержне, и с ее помощью найти температурный режим стержня при отсутствии внешних источников. Решение. Выбираем координатную систему так, чтобы ось x совпала с положением стержня. Пусть u(x, t) — температура стержня в точке x в момент t. Тогда функция u(x, t) удовлетворяет уравнению теплопроводности 𝜕𝑢 𝜕𝑡 = 𝑎2 𝜕2 𝑢 𝜕𝑥 2 (2.1.1)

Как известно, для стержня конечной длины задача отыскания функции Грина сводится к решению уравнения (2.1.1) с заданными однородными граничными условиями и начальным условием в виде δ-функции. Попробуем и для бесконечного стержня найти решение u0(x,t) уравнения (2.1.1), для которого u0(x,0)=δ(ξ−x) и покажем, что оно обладает свойствами, аналогичными свойствам функции Грина в задаче о тепловом режиме стержня конечной длины. Однако, чтобы оставаться в рамках математической строгости, не прибегая к теории обобщенных функций, мы возьмем начальное условие в виде δ-последовательности δn(ξ, x)∗, найдем отвечающее δn(ξ, x) решение un(x,t) и потом перейдем к пределу при n → ∞. Итак, пусть 𝜕𝑢𝑛 𝜕𝑡 = 𝑎2 𝜕 2 𝑢𝑛 𝜕𝑥 2 , (2.1.2) 𝑢𝑛 (𝑥, 0) = 𝛿𝑛 (𝜉, 𝑥), (2.1.3) где ξ — произвольное число, δn(ξ,x) ∈ C(−∞, +∞), δn(ξ,x)= δn(x, ξ), 𝛿𝑛 (𝜉, 𝑥) = { 0, |𝑥 − 𝜉| ≥ 1 𝑛 ≥ 0, |𝑥 − 𝜉| < 1 , ∫ 𝛿𝑛 (𝜉, 𝑥) 𝑑𝑥 = 1. (2.1.4) 𝑛 Пример. Найдём преобразование Фурье прямоугольного импульса 2𝑡 1, |𝑥| ≤ 𝑎 𝑓(𝑥) = 𝜂(𝑎 − |𝑥|) = { ≡ 𝑟𝑒𝑐𝑡 ( ) , 𝑎 > 0 ∶ 0, |𝑥| > 𝑎 𝑎 +∞ +𝑎 𝑓̂(𝑦) = ∫−∞ 𝑓(𝑥)𝑒 −𝑖𝑥𝑦 𝑑𝑥 = ∫−𝑎 𝑒 −𝑖𝑥𝑦𝑑𝑥 = Функция 𝑠𝑖𝑛 𝑦 𝑦 𝑒 𝑖𝑦𝑎 −𝑒 −𝑖𝑦𝑎 𝑖𝑦 = 2 𝑠𝑖𝑛 𝑦𝑎 𝑦 . ̂ = 𝑠𝑖𝑛𝑐 𝑦. обозначается как sinc y, то есть 𝑟𝑒𝑐𝑡𝑥 2 Пример. Найдём преобразование Фурье функции 𝑓(𝑥) = 𝑒 −𝑎𝑥 : 2 +∞ 𝑓̂(𝑦) = ∫−∞ 𝑒 −𝑎𝑥 𝑒 −𝑖𝑥𝑦 𝑑𝑥.

Преобразуем степень экспоненты, выделив полный квадрат: 𝑖𝑦 −𝑎𝑥 2 − 𝑖𝑥𝑦 = −𝑎(𝑥 + Тогда 𝑓(𝑦) = 𝑒 −𝑦2 4𝑎 2𝑎 )2 − 𝑦2 4𝑎 . 𝑖𝑦 2 +∞ −𝑎(𝑥+ ) 2𝑎 𝑑𝑥 . ∫−∞ 𝑒 Под интегралом аналитическая в C функция переменной x ∈ C, стремящаяся к нулю вдоль каждой прямой, параллельной вещественной оси (при |ℝ x| → ∞). Поэтому в силу теоремы Коши и леммы Жордана интеграл не изменит своего значения, если его взять не по действительной оси, а по любой другой прямой, параллельной этой оси. Возьмём его по прямой 𝑥 = 𝑧 − +∞). +∞ Тогда 𝑓(𝑦) = 𝑒 − 𝑦2 4𝑎 𝑑𝑧 = 𝑒 − 𝑦2 4𝑎 𝑖𝑦 2𝑎 , z ∈ (−∞, 𝜋 √𝑎 так , как интеграл Пуассона 2 ∫−∞ 𝑒 −𝑧 𝑑𝑧 = √𝜋 . 𝜋 Итак, 𝑓(𝑦) = √ 𝑒 𝑎 − 𝑦2 4𝑎 . • Преобразования Фурье обобщенных функций Пример. Дельта-функция Дирака может быть определена по формуле (𝛿, 𝜑) = 𝜑(0) для любой непрерывной на R функции ϕ, в частности, для любой 𝜑 ∈ 𝑆(ℝ). Таким образом, 𝛿 ∈ 𝑆′(ℝ). Теперь определим сходимость в пространстве 𝑆′(ℝ). Определяется она аналогично сходимости в 𝐷′(ℝ). Пример. Найдем преобразование Фурье δ-функции Дирака: +∞ (𝛿̂ , 𝜑) = (𝛿, 𝜑̂) = 𝜑̂(0) = ∫ −∞ +∞ 𝜑(𝑦)𝑒 𝑖𝑥𝑦 𝑑𝑦| 𝑥=0 =∫ 𝜑(𝑦). 1𝑑𝑦 = (1, 𝜑). −∞ Таким образом, 𝛿̂ = 1. Пример. Найдем преобразование Фурье обобщенной функции f(x) ≡ 1. Заметим, что обычное преобразование Фурье не существует, т.к. константа 1

̂ ̃ = не является абсолютно интегрируемой на ℝ. Итак, 1̂ = 𝛿̂ = 2𝜋𝛿(−𝑥) 2𝜋𝛿(−𝑥) = 2𝜋𝛿. • Двумерное преобразование Фурье. Задача о распространении волн на плоскости Двумерное преобразование Фурье функции g(x, y) есть по определению функция 𝑔̅ (𝑝, 𝑞) = ∬ 𝑔(𝑥, 𝑦)𝑒 −𝑖(𝑝𝑥+𝑞𝑦) 𝑑𝑥𝑑𝑦 (2.1.1’) Для применимости этого преобразования достаточно, чтобы функция |g(x, y)| была суммируема на всей плоскости XY = {(x, y) | −∞ < x, y < +∞}. Формула обращения для двумерного ИП Фурье дается равенством: 𝑔(𝑥, 𝑦) = Условие её справедливости 1 4𝜋2 — ∬ 𝑔̅ (𝑝, 𝑞)𝑒 𝑖(𝑝𝑥+𝑞𝑦) 𝑑𝑝𝑑𝑞 это условие (2.1.2’) справедливости двух последовательно примененных одномерных формул обращения. Поскольку двумерное ИП Фурье можно рассматривать как произведение двух одномерных (по разным переменным), то для образов вторых производных 𝜕2 𝑔 𝜕𝑥 2 и 𝜕2 𝑔 𝜕𝑦 2 мы получим равенства: 𝜕 2 𝑔 −𝑖𝑝𝑥 −𝑖𝑞𝑦 𝜕 2 𝑔 −𝑖𝑝𝑥 ∬ 2𝑒 𝑒 𝑑𝑥𝑑𝑦 = ∫(∫ 2 𝑒 𝑑𝑥)𝑒 −𝑖𝑞𝑦 𝑑𝑦 = 𝜕𝑥 𝜕𝑥 = −𝑝2 ∫ 𝑔̅1 (𝑝, 𝑦)𝑒 −𝑖𝑞𝑦 𝑑𝑦 = −𝑝2 𝑔̅ (𝑝, 𝑞) (2.1.3a) 𝜕 2 𝑔 −𝑖𝑝𝑥 −𝑖𝑞𝑦 𝜕 2 𝑔 −𝑖𝑞𝑦 ∬ 2𝑒 𝑒 𝑑𝑥𝑑𝑦 = ∫(∫ 2 𝑒 𝑑𝑦)𝑒 −𝑖𝑝𝑥 𝑑𝑥 = 𝜕𝑦 𝜕𝑦 = −𝑞 2 ∫ 𝑔̅2 (𝑥, 𝑞)𝑒 −𝑖𝑝𝑥 𝑑𝑥 = −𝑞 2 𝑔̅ (𝑝, 𝑞) где 𝑔̅1 (𝑝, 𝑦) = ∫ 𝑔(𝑥, 𝑦)𝑒 −𝑖𝑝𝑥 𝑑𝑥 ,𝑔̅2 (𝑥, 𝑞) = ∫ 𝑔(𝑥, 𝑦)𝑒 −𝑖𝑞𝑦 𝑑𝑦, (2.1.3b)

и мы предполагаем, что gxx, gyy ∈ L1(XY ), g, g'x , g''xx ⊂L1(−∞, +∞) при любом фиксированном y и g, g'y , g''yy ∈ L1(−∞, +∞) при любом фиксированном x. Задача. Рассмотрим задачу о распространении волн на плоскости. Обозначим через u(x, y, t) возмущение в точке (x, y) в момент t. Тогда, как известно, функция u(x, y, t) должна удовлетворять уравнению 𝜕2 𝑢 𝜕𝑡 2 = 𝑎2 ( 𝜕2 𝑢 𝜕2 𝑢 𝜕𝑥 𝜕𝑦 2 + 2 ) (2.1.4’) Пусть начальное возмущение и начальная скорость распространения возмущений суть функции ϕ(x, y) и ψ(x, y), т. е. a) u(x, y, 0) = ϕ(x, y), b) ut(x, y, 0) = ψ(x, y). (2.1.5’) Делаем двойное преобразование Фурье, умножая (2.1.4’), (2.1.5’) на e−ipx−iqy и интегрируя по x и по y от −∞ до +∞. Тогда в силу (2.1.1’), (2.1.3a), (2.1.3b) получим ̅ 𝑑2𝑢 𝑑𝑥 2 = −𝑎2 (𝑝2 + 𝑞 2 )𝑢̅, (2.1.6) a) 𝑢̅(𝑝, 𝑞, 0) = 𝜙̅(𝑝, 𝑞), b) 𝑢̅𝑡 (𝑝, 𝑞, 0) = 𝜓̅(𝑝, 𝑞), (2.1.7) где 𝜙̅(𝑝, 𝑞) = ∬ 𝜙(𝑥, 𝑦) 𝑒 −𝑖(𝑝𝑥+𝑞𝑦) 𝑑𝑥 𝑑𝑦, 𝜓̅(𝑝, 𝑞) = ∬ 𝜓(𝑥, 𝑦) 𝑒 −𝑖(𝑝𝑥+𝑞𝑦) 𝑑𝑥 𝑑𝑦. Из (8.6) получаем, что 𝑢̅(𝑝, 𝑞, 𝑡) = 𝑑1 𝑠𝑖𝑛𝑎𝜌𝑡 + 𝑑2 𝑐𝑜𝑠𝑎𝜌𝑡, где 𝜌 = √𝑝2 + 𝑞 2 , откуда в силу начальных условий (8.7) имеем 𝑎𝜌𝑡 𝑢̅(𝑝, 𝑞, 𝑡) = 𝜙̅(𝑝, 𝑞)𝑐𝑜𝑠𝑎𝜌𝑡 + 𝜓̅(𝑝, 𝑞)𝑠𝑖𝑛 . 𝑎𝜌 (2.1.8) Далее в конкретной задаче, подставив в (2.1.8) значения 𝜙̅(𝑝, 𝑞) и 𝜓̅(𝑝, 𝑞) и воспользовавшись формулой обращения (2.1.2), получим

𝑢(𝑥, 𝑦, 𝑡) = 1 4𝜋2 ∬[𝜙̅(𝑝, 𝑞) cos 𝑎𝜌𝑡 + 𝜓̅(𝑝, 𝑞) sin 𝑎𝜌𝑡 . 𝜌−1 𝑎−1 ] × exp[𝑖(𝑝𝑥 + 𝑞𝑦)]𝑑𝑝 𝑑𝑞. (2.1.9) Однако зависимость решения u(x, y, t) от начальных условий можно получить и не находя 𝜙̅(𝑝, 𝑞), 𝜓̅(𝑝, 𝑞). Для этого надо воспользоваться формулой свертки для двумерного ИП Фурье и рассуждениями. 2.2. Применение преобразования Фурье для решения интегральных уравнений Преобразование Фурье может быть использовано для решений некоторых интегральных уравнений, то есть уравнений, для которых неизвестная функция входит под знак интеграла. Рассмотрим, например, уравнение: +∞ ∫−∞ 𝜑(𝑥)𝑒 −𝑖𝜉𝑥 𝑑𝑥 = 2𝜋𝑒 −|𝜉| , (2.2.1) где 𝜑(𝑥)- искомая функция. Записав (2.2.1) в виде 1 +∞ ∫ 𝜑(𝑥)𝑒 −𝑖𝜉𝑥 𝑑𝑥 = √2𝜋𝑒 −|𝜉| √2𝜋 −∞ (2.2.2) Замечаем, что левую часть (2.2.2) можно рассматривать как преобразование Фурье функции 𝜑(𝑥),так что (2.2.2) равносильно следующему равенству: ℱ[𝜑] = √2𝜋𝑒 −|𝜉| Тогда, по формуле обращения +∞ 𝜑(𝑥) = ∫ 0 +∞ 𝑒 −|𝜉| 𝑒 𝑖𝜉𝑥 𝑑𝜉 = ∫ 𝑒 𝜉(1+𝑖𝑥) 𝑑𝜉 + ∫ −∞ −∞ 𝜉=0 𝑒 −𝜉(1+𝑖𝑥) 𝑑𝜉 = 0 𝜉=+∞ 𝑒 𝜉(1+𝑖𝑥) 𝑒 −𝜉(1+𝑖𝑥) = − | | 1 + 𝑖𝑥 𝜉=−∞ 1 − 𝑖𝑥 𝜉=0 = 2 1 + 𝑥2

Функция 𝜑(𝑥) = 2 1+𝑥 2 есть решение уравнения (2.2.1) 2.3. Преобразование Фурье и ряды Фурье, дискретное преобразование Фурье Ряд Фурье. Ряд Фурье - это представление произвольной функции с некоторым периодом в виде ряда. В общем случае такое решение называется разложением элемента на ортогональном базисе. Разложение функций в ряд Фурье является достаточно мощным инструментом для решения множества задач, связанных с свойствами этого преобразования в интегрировании, дифференцировании, а также путем перемещения выражения в аргумент и свертку. Рядом Фурье функции f(x) на интервале (-π;π) называется тригонометрический ряд вида: 𝑓(𝑥) ≅ 𝑎0 2 + ∑∞ 𝑛=1 𝑎𝑛 𝑐𝑜𝑠𝑛𝑥 + 𝑏𝑛 𝑠𝑖𝑛𝑛𝑥 Где, 1 𝜋 1 𝜋 𝜋 1 𝑎0 = ∫−𝜋 𝑓(𝑥)𝑑𝑥 , 𝑎𝑛 = ∫−𝜋 𝑓(𝑥)𝑐𝑜𝑠𝑛𝑥𝑑𝑥 , 𝑏𝑛 = ∫−𝜋 𝑓(𝑥)𝑠𝑖𝑛𝑛𝑥𝑑𝑥. 𝜋 𝜋 𝜋 Рядом Фурье функции f(x) на интервале (-l;l) называется тригонометрический ряд вида: , где .

Вопрос о сходимости срабатывает всякий раз, когда необходимо суммирование бесконечных чисел. Чтобы понять это явление, рассмотрим классический пример. Можете ли вы достигнуть стены, если каждый шаг половину предыдущего? Предположим, вы находитесь в двух метрах от цели, первый шаг ближе к середине, следующий знак в три четверти, а пятый, вы преодолеете почти 97 процентов пути. Но независимо от того, сколько шагов вы делаете, цель состоит в том, чтобы достичь строгого математического значения. Используя численные вычисления, мы можем доказать, что в конечном итоге вы можете получить как можно ближе к заданному короткому расстоянию. Это доказательство эквивалентно доказательству того, что совокупность второго, четвертого значения и так далее стремится к одному. Известно, что периодическая с периодом X ограниченная кусочно-гладкая функция f(x) может быть представлена своим рядом Фурье: 2𝑛𝑠𝑖𝑛𝑥/𝑋 𝑓(𝑥) = ∑𝑛=+∞ , 𝑛=−∞ 𝐹𝑛 𝑒 где 1 𝑋/2 𝐹𝑛 = ∫ 𝑓(𝑥) 𝑒 −2𝑛𝑠𝑖𝑛𝑥 𝑑𝑥 𝑋 −𝑋/2 Рассмотрим снова пространство L2[-;] функций с суммируемым квадратом на отрезке [-;]. Это -полное бесконечномерное евклидово пространство, те. гильбертово пространство. Функции 1, cos nx, sin nx (n = 1, 2, …) образуют в нем полную ортогональную систему, поэтому для каждой функции 𝑓𝜖𝐿2 [−𝜋; 𝜋] ряд Фурье 𝑎0 2 Где + ∑∞ 𝑛=1 𝑎𝑛 𝑐𝑜𝑠𝑛𝑥 + 𝑏𝑛 𝑠𝑖𝑛𝑛𝑥 (2)

1 𝜋 𝑎𝑛 = ∫−𝜋 𝑓(𝑥)𝑐𝑜𝑠𝑛𝑥𝑑𝑥, 𝜋 Сходится к f 1 𝜋 𝑏𝑛 = ∫−𝜋 𝑓(𝑥)𝑠𝑖𝑛𝑛𝑥𝑑𝑥 𝜋 (3) в среднем квадратичном, T. e. в метрике пространства 𝐿2 [−𝜋; 𝜋]. Однако, в связи с применением рядов Фурье к задачам математической физики и другим вопросам будет существенно установить условия, гарантирующие сходимость ряда Фурье к не только в среднем, но и в данной точке, всюду, или даже равномерно. Мы установим сейчас условия, достаточные для сходимости тригонометрического ряда в данной точке. Сделаем некоторые предварительные замечания. Вместо функций, заданных на отрезке [−𝜋; 𝜋], мы можем говорить о периодических функциях с периодом 2 на всей прямой, поскольку каждую функцию, заданную на отрезке, можно периодически продолжить. Далее, функции, образующие тригонометрическую систему, ограничены, поэтому формулы (3), определяющие коэффициенты Фурье по этой системе, имеют смысл для любой суммируемой функции) (а не только для функций с суммируемым квадратом) Таким образом, каждой функции 𝑓𝜖𝐿1 [−𝜋; 𝜋] отвечают совокупность ее коэффициентов Фурье и ее ряд Фурье Преобразование Фурье. Для непериодических функций ряд Фурье заменяется интегралом Фурье: 𝑓(𝑥) = 1 √2𝜋 +∞ ∫ 𝐹(𝜒) 𝑒 𝑖𝜒𝑥 𝑑𝑥, −∞ где 𝐹(𝜒) = 1 √2𝜋 +∞ ∫ 𝑓(𝑥) 𝑒 −𝑖𝜒𝑥 𝑑𝑥 −∞

Связь между рядом Фурье и преобразованием Фурье. Рассмотрим функцию f(x), равную нулю вне интервала [-X/2;X/2]. С одной стороны, для нее можно определить преобразование Фурье, а с другой стороны, ее можно считать периодически продолженной и определить для нее коэффициенты ряда Фурье. Тогда имеем + 𝐹(𝜒) = 1 √2𝜋 𝑋 2 ∫ 𝑓(𝑥) 𝑒 −𝑖𝜒𝑥 𝑑𝑥, − 𝑋 2 +𝑋/2 𝐹𝑛 = 1 ∫ 𝑓(𝑥)𝑒 −2𝑛𝑠𝑖𝑛𝑥/𝑋 𝑑𝑥. 𝑋 −𝑋/2 Сравнивая полученные выражения, можно сделать вывод, что 𝐹 (𝜒𝑛 = 2𝜋𝑛 𝑋 𝐹𝑛 , )= 𝑋 √2𝜋 т.е. коэффициенты ряда Фурье для периодического продолжения функции f(x) определяют значение преобразования Фурье в дискретных точках 𝜒𝑛 = 2𝜋𝑛 𝑋 Дискретное преобразование Фурье. Для последовательности {fn, действительных или n=0,1,2,…,N-1}, комплексных чисел состоящей определяется преобразование Фурье (ДПФ){Dn, n=0,1,2,…,N-1}: 𝑁−1 1 𝐷𝑛 = ∑ 𝑓𝑚 𝑊 𝑚𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1 𝑁 𝑚=0 из N дискретное

𝑊 = 𝑒 −2𝜋𝑖/𝑁 — где Так как дискретные дискретные экспоненциальные экспоненциальные функции Wmn функции. являются ортогональными, т.е. удовлетворяют условию 𝑁−1 𝑁, 𝑚1 ≡ 𝑚2 (𝑚𝑜𝑑 𝑁) ∑ 𝑊 𝑚1𝑛 𝑊 −𝑚2𝑛 = { 0, 𝑚1 ≠ 𝑚2 (𝑚𝑜𝑑 𝑁) 𝑚=0 то справедливо обратное дискретное преобразование Фурье (ОДПФ): 𝑁−1 𝑓𝑛 = ∑ 𝐷𝑚 𝑤 −𝑚𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1 𝑚=0 Приведем некоторые свойства ДПФ. 1) Теорема линейности. Если даны последовательности такие, что , и , , — ДПФ , , ,соответственно, то 2) Теорема комплексной сопряженности. Если -четное последовательность а число и - действительных -ДПФ последовательности чисел, , то

, где - величина, комплексно-сопряженная величине частности, коэффициенты и . В всегда действительны. 3) Теорема циклического сдвига. Если даны последовательности , , такие, что , и — , ДПФ соответственно последовательностей , то , . 4) Теорема циклической свертки. Если даны последовательности , такие, что , , и , — ДПФ то , , , ,соответственно,

, причем последовательность последовательностей и называют циклической сверткой . 5) Теорема циклической корреляции. Если три последовательности {𝑓𝑛 } = {𝑓𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1} , {𝑓𝑛′ } = {𝑓𝑛′ , 𝑛 = 0,1,2, … , 𝑁 − 1}, {𝑓𝑛" } = {𝑓𝑛" , 𝑛 = 0,1,2, … , 𝑁 − 1} такие, что 1 ′ " 𝑓𝑛 = ∑𝑁−1 𝑚=0 𝑓𝑚 𝑓𝑚−𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1, 𝑁 и {𝐷𝑛 } = {𝐷𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1}, {𝐷𝑛′ } = {𝐷𝑛′ , 𝑛 = 0,1,2, … , 𝑁 − 1}, {𝐷𝑛" } = {𝐷𝑛" , 𝑛 = 0,1,2, … , 𝑁 − 1}— ДПФ {fn}, {fn’}, {fn”},соответственно, то 𝐷𝑛 = ̅̅̅̅ 𝐷𝑛′ 𝐷𝑛" , 𝑛 = 0,1,2, … , 𝑁 − 1, При этом последовательность {fn} называют циклической функцией корреляции последовательностей {fn’} и {fn”}. Если {fn’} и {fn”} — идентичные последовательности, то 𝐷𝑛 = |𝐷𝑛′ |2 , 𝑛 = 0,1,2, … , 𝑁 − 1. Находя ОДПФ последовательности {Dn}, получим 𝑁−1 𝑁−1 𝑓𝑛 = ∑ 𝐷𝑚 𝑊 −𝑚𝑛 = ∑ |𝐷𝑛′ |2 𝑊 −𝑚𝑛 𝑚=0 𝑚=0 𝑁−1 1 ′ = ∑ 𝑓𝑚′ 𝑓𝑚+𝑛 , 𝑛 = 0,1,2, … , 𝑁 − 1 𝑁 𝑚=0 В частном случае при n=0 это выражение сводится к 1 𝑁−1 ′ 2 ′ 2 ∑𝑁−1 𝑚=0|𝐷𝑛 | = ∑𝑚=0(𝑓𝑚 ) , 𝑁 что является теоремой Парсеваля последовательности {f’n}={f’n,n=0,1,2,…,N-1}. для дискретной

Связь между ДПФ, рядом Фурье и преобразованием Фурье. Для того чтобы установить связь между ДПФ, рядом Фурье и преобразованием Фурье рассмотрим дискретизацию непрерывной функции f(x), равной нулю вне интервала [-X/2, +X/2]. Заменим на ее дискретное представление f*(x) в виде суммы -функций в точках дискретизации xn = nx (n =…,-1,0,1,…) 𝑁 𝑛= −1 2 𝑛=+∞ 𝑓 ∗ (𝑥) = ∑ 𝑓(𝑛∆𝑥)∆𝑥𝛿(𝑥 − 𝑛∆𝑥) = 𝑛=−∞ ∑ 𝑓𝑛 ∆𝑥𝛿(𝑥 − 𝑛∆𝑥) 𝑛=−𝑁/2 где N=X/x. Таким образом, для функции f*(x) можно определить как преобразование Фурье, так и, считая ее периодически продолженной, коэффициенты ряда Фурье. Тогда для коэффициентов ряда Фурье получим + 𝑋 2 2𝜋𝑖𝑛𝑥 1 ∫ 𝑓 ∗ (𝑥)𝑒 − 𝑋 𝑑𝑥 = 𝑋 𝐹𝑛 = − 𝑋 2 𝑋 + 𝑚=𝑁−1 2 2 1 = ∫ 𝑋 − = ∆𝑥 𝑋 𝑋 2 ∑ 𝑓𝑚 ∆𝑥𝛿(𝑥 − 𝑚=− 𝑁 𝑚= −1 2 1 𝑁 𝑁 2 2𝜋𝑖𝑛𝑚∆𝑥 𝑋 𝑑𝑥 ∑ 𝑓𝑚 𝑒 − 𝑚=− = 2𝜋𝑖𝑛𝑥 − 𝑚∆𝑥)𝑒 𝑋 𝑑𝑥 𝑁 2 𝑁 𝑚= −1 2 ∑ 𝑚=−𝑁/2 𝑓𝑚 𝑒 −2𝜋𝑖𝑛𝑚/𝑁 = = =

= = 1 𝑁 1 𝑁 𝑁 𝑚= −1 2 𝑚=−1 𝑓𝑚 𝑒 −2𝜋𝑖𝑛𝑚/𝑁 + ∑ 𝑚=−𝑁/2 1 ∑ 𝑓𝑚 𝑒 −2𝜋𝑖𝑛𝑚/𝑁 = 𝑁 𝑚=0 𝑁 𝑚= −1 2 𝑚=𝑁−1 𝑓𝑚−𝑁 𝑒 −2𝜋𝑖(𝑚−𝑛)/𝑁 + ∑ 𝑚=−𝑁/2 1 ∑ 𝑓𝑚 𝑒 −2𝜋𝑖𝑛𝑚/𝑁 = 𝑁 𝑚=0 𝑚=𝑁−1 𝑁−1 𝑚=0 𝑚=0 1 1 = ∑ 𝑓𝑚 𝑒 −2𝜋𝑖𝑛𝑚/𝑁 = ∑ 𝑓𝑚 𝑊 𝑛𝑚 = 𝐷𝑛 𝑁 𝑁 При этом, в соответствии с циклическим продолжением, считая, что fm=fm-N , если m>(N/2)-1. Аналогично, для преобразования Фурье 𝐹(𝜒) в дискретных точках 𝜒𝑛 = 2𝜋𝑛/𝑋 получим 𝐹(𝜒𝑛 ) = = = √2𝜋 √2𝜋 √2𝜋 ∫ 𝑓 ∗ (𝑥)𝑒 −𝑖𝜒𝑛𝑥 𝑑𝑥 = −∞ 𝑁 +∞ 𝑚= 2 −1 1 1 +∞ 1 ∫ ∑ 𝑓𝑚 ∆𝑥𝛿(𝑥 − 𝑚∆𝑥)𝑒 −𝑖𝜒𝑛𝑥 𝑑𝑥 −∞ 𝑚=−𝑁/2 𝑁 2 𝑚= −1 ∑𝑚=−𝑁/2 𝑓𝑚 Δ𝑥𝑒 −𝑖𝜒𝑛 𝑚Δ𝑥 = Δ𝑥 1 √2𝜋 1 𝑁 2 𝑚= −1 𝑖2𝜋𝑛𝑚 𝑁 𝑁 ∑𝑚=−𝑁/2 𝑓𝑚 𝑒 − 𝑁 = 𝑋 √2𝜋 𝐷𝑛 . Таким образом, ДПФ является дискретным аналогом коэффициентов ряда Фурье и преобразования Фурье. Быстрое преобразование Фурье. Быстрое преобразование Фурье (БПФ, FFT) - это алгоритм быстрого вычисления дискретного преобразования Фурье (ДПФ). То есть, алгоритм вычисления количества действий меньше О(N2), необходимый для прямого вычисления (по формуле) ДПФ. Иногда FFT является одним из самых

быстрых алгоритмов, называемым алгоритмом прореживания частоты/ времени, который имеет сложность O(Nlog(N)).. Быстрое преобразование Фурье (БПФ) является таким способом вычисления дискретного преобразования Фурье для некоторых значений N, что дает резкое сокращение числа операций, необходимых для вычисления ДПФ. Проще запрограммировать мощности числа БПФ для N, равных мощности 2. Поэтому при использовании БПФ обычно предполагается, что N=2k. Для других значений N коэффициент усиления также может быть достигнут выигрыш, тем больше, чем больше количество множителей в разложении числа N на простые множители. Однако для определения реального выигрыша необходимо учитывать сложность вычисления элементарного Pi — точечного ДПФ, где Pi – все различные простые факторы, участвующие в разложении N на простые множители: 𝒏 𝑵 = ∏𝒊 𝒑𝒊 𝒊 . Если для прямого вычисления ДПФ необходимого порядка N требуются операции умножения на комплексные числа, тогда БПФ требует только порядка (N/2)log2N операций умножения для N=2k. Для более точного определения коэффициента усиления при использовании БПФ необходимо также учитывать операции сложения, а также различные времена выполнения операций сложения и умножения. Существует большое количество алгоритмов БПФ. Однако все они являются частными случаями одного алгоритма, основанного на проблеме деления одного массива чисел на два. Тот факт, что это может быть сделано более чем одним способом, определяет разнообразие алгоритмов БПФ. Рассмотрим два из них. Первый алгоритм называется алгоритмом БПФ с прореживанием по времени.

Пусть N-отсчетный дискретный сигнал образца x(n). Предположим, что N равно степени 2. Если это не так, всегда можно легко дополнить сигнал, нулевыми выборками, до количества ударов, равных мощности, ближайшей к двум. Разобьём исходный сигнал x(n) на два N/2-отсчетных сигнала x1(n) и x2(n), составленных соответственно из четных и нечетных отсчетов исходного сигнала x(n) (а ). N-точечное ДПФ сигнала x(n) можно записать как (b ). С учетом того, что (c) можно записать (d) или (e), где X1(k) и X2(k) соответственно. – N/2-отсчетные ДПФ сигналов x1(n) и x2(n)

Таким образом N-точечное ДПФ X(k) может быть разложено на два N/2точечных ДПФ, результаты которых объединяются согласно (e) Если бы N/2-точечные ДПФ вычислялись прямым способом, то для вычисления N-точечного ДПФ потребовалось бы (N2/2+N) комплексных умножений. При больших N (когда N2/2>>N) это позволяет сократить время вычислений на 50%. Поскольку X(k) определено при при , необходимо , а X1(k) и X2(k) определены доопределить формулу (e) для . Учитывая, что X1(k) и X2(k) – периодические функции с периодомN/2, можно записать (f), поскольку . Поэтому окончательно для N-точечного ДПФ можно записать (g). На рис.A представлена последовательность операций при выполнении восьмиточечного ДПФ с использованием двух четырехточечных ДПФ. Во-первых, входной сигнал x(n) делится на два сигнала x1(n) и x2(n), составленных соответственно из нечетных и четных выборок x(n). После этого вычисляются ДПФ X1(k) и X2(k). Тогда, согласно (g), получим X(k). Рассматриваемая схема расчета также может быть использована для вычисления N/2-точечных ДПФ. Соответственно, каждый из сигналов x1(n) и x2(n) разделен на последовательности, состоящие из нечетных и четных

отсчетов исходных сигналов. Аналогично, N/2-точечные ДПФ могут быть записаны как комбинации двух N/4-точечных ДПФ (h). С учетом того, что можно записать (i). На рис.B представлена последовательность операций при выполнении восьмиточечного ДПФ с использованием двух четырехточечных ДПФ и четырех двухточечных ДПФ. Таким образом, процесс уменьшения размера ДПФ может быть продолжен до тех пор, пока не останутся только двухточечные ДПФ, которые могут быть рассчитаны без операции умножения. (j). Поскольку , то окончательно получим (k). На рис.C представлена порядок операций при последовательном вычислении восьмиточечного ДПФ в соответствии с описанным алгоритмом. Анализ рисунки C показывает, что на каждом этапе БПФ необходимо выполнить комплексные умножения N/2 раз. Поскольку общее число шагов равно , то количество комплексных умножений,

необходимое для нахождения N-точечного ДПФ, приблизительно равно . Грубость оценки означает, что умножения фактически сводятся просто к сложениям и вычитаниям комплексных чисел. Таким образом, на первом этапе алгоритма, показанного на рисунке С, содержатся только сложения и вычитания комплексных чисел поскольку . Даже на втором шаге используются только сложения и вычитания комплексных чисел поскольку . На самом деле, вместо ожидаемых 12 ( ), достаточно выполнить только два нетривиальных умножения. Однако при больших значениях N фактическое число нетривиальных умножений хорошо аппроксимируется выражением . Базовая операция алгоритма называемая «бабочка») состоит с прореживанием в том, по что времени два (так входных числа A и B объединяются для получения двух выходных чисел X и Y по правилу (l). На рис.D изображен направленный граф базовой операции. Каждый из шагов в ДПФ содержит N/2 основных операций. В случае, когда является нетривиальным фактором, для каждой базовой операции необходимо выполнить лишь одно умножение, так как значение можно вычислить и сохранить. Таким образом, структура основных операций такова, что для выполнения БПФ N-отсчетного сигнала, чьи отсчеты расположены в памяти, достаточно иметь только одну дополнительную ячейку памяти. Результаты всех промежуточных шагов БПФ могут быть помещены в те же ячейки памяти, где исходные данные. Алгоритм БПФ, в

котором тот же массив памяти используется для размещения входной и выходной последовательности, называется алгоритмом замещением. Чтобы реализовать алгоритм, необходимо переставить отсчетов входного сигнала, так что выходная последовательность имела естественный порядок расположения отсчетов, то есть k = 0,1 , ..., N-1. В примере, приведенный для 8-точечного БПФ, для этого требовался следующий порядок организации выборок входного сигнала: x(0), x(4), x(2), x(6), x(1), x(5), x(3), x(7). Регулярность перестановки заключается в том, что образцы входного сигнала должны быть помещены в память в двоично-инверсном порядке. Это означает, что номер ячейки памяти, необходимый для размещения следующего отсчета входного сигнала, определяется обратной перестановкой двоичных разрядов в двоичном представлении номера отсчета. Для случая N=8 соответствие номеров отсчетов входного сигнала и их номеров ячеек памяти представлено в таблице 1. Н омер Двоичное представление Двоичная инверсия Двоичноинверсный порядок 0 000 000 0 1 001 100 4 2 010 010 2 3 011 110 6 4 100 001 1 5 101 101 5 6 110 011 3 7 111 111 7

На рис.E показан алгоритм определения двоично-инверсного номера, предложенный Рейдером. Начиная с первого числа с прямым номером 0, этот алгоритм позволяет формировать последовательно все остальные двоично-инверсные номера. Заключение В первой главе мы рассмотрели понятие и основные свойства интегрального преобразования Фурье, а также свертка функций и примеры преобразований Фурье. Во второй главе мы определили применения дискретного преобразования Фурье к решению задач математической физики и к решению интегральных уравнений. Окончательно, мы рассмотрели ряды Фурье и дискретное преобразование Фурье. В этой главе мы рассмотрели представления по Фурье последовательностей конечной длины, названные дискретным преобразованием Фурье. Это представление было основано на связи между последовательностями конечной длины и периодическими последовательностями. А именно, если периодическая последовательность сформирована так, что каждый ее период идентичен последовательности конечной длины, то дискретное преобразование Фурье последовательности конечной длины соответствует дискретному ряду Фурье периодической последовательности. Поэтому мы сначала рассмотрели представление периодических последовательностей рядами Фурье, включая свойства рядов Фурье и интерпретацию их коэффициентов, как равноудаленных выборок на единичной окружности из 2-преобразований одного периода периодических последовательностей. Затем дискретные ряды Фурье были применены к представлению последовательностей конечной длительности. Были рассмотрены свойства

дискретного преобразования Фурье, а также вопросы вычисления линейной свертки с использованием дискретного преобразования Фурье. Список литературы 1) Оппенгейм А. В., Р.В. Шафер «Цифровая обработка сигналов» М.: Связь, 1979, с 74-96 2) А. Н. Колмогоров, С. В. Формин «Элементы функционального анализа» -М.:Наука. 1976, с 423-435 3) Brian Davis “Integral transforms and their applications”, third edition, Spinger Australia, 2001, с 111-123 4) Gerald B. Folland “Fourier analysis and its applications”, Wadsworth & Brooks/Cole Advanced Books and Software Pacific Groove. California 1992, с 18-57, 204-225 5) Борзенков А. В. «Дифференциальные уравнения в частных производных» -Минск БГУИР 2009, с 69 6) Богданов, Ю. С. «Дифференциальные уравнения» Ю. С. Богданов, Ю.Б.Сыроид. –М., 1996

7) Владимиров, В. С. «Уравнения математической физики»/ В. С. Владимиров. – М., 1981. 8) Васильева А.Б., Тихонов Н.А. «Интегральные уравнения.» М.: Физматлит, 2004. 9) Васильева А.Б., Медведев Г.Н., Тихонов Н.А., Уразгильдина Т.А. «Дифференциальные и интегральные уравнения, вариационное исчисление в примерах и задачах»/ М.: Физматлит, 2003. 10) Волков В.Т., Ягола А.Г. «Интегральные уравнения. Вариационное исчисление: методы решения задач.» М.: КДУ, 2007. 11) Краснов М.Л. «Интегральные уравнения. Введение в теорию.» М.: КомКнига, 2006. 12) Краснов М.Л., Киселев А.И., Макаренко Г.И. «Интегральные уравнения: Задачи и примеры с подробными решениями.» М.: Едиториал УРСС, 2003. 13) Под ред. А.В. Ефимова и А.С. Поспелова. «Сборник задач для втузов.» Ч. 3./ М.: Физматлит, 2003. 14) Жислин Г. М. «Интегральные преобразования в задачах математической физики» Электронное учебно-методическое пособиеуравнениям. М.: Физматлит, Нижний Новгород: Нижегородский госуниверситет, 2012. – 80 с 15) Тихонов А. Н., Самарский А. А. «Уравнения математической физики.» — М.: Наука, 1977. 16) Трантер К. Д. «Интегральные преобразования в математической физике.» — М.: Гос. изд-во технико-теоретической лит- ры, 1956. 17) Диткин В. А., Прудников А. П. «Интегральные преобразования и операционное исчисление.» — М.: Физматгиз, 1961.

18) Привалов И. И. «Введение в теорию функций комплексного переменного.» — М.: Наука, 1984 19) Лебедев Н. Н. «Специальные функции и их приложения.» — М.: Физматгиз, 1963. 20) Прудников А. П., Брычков Ю. А., Маричев О. И. «Интегралы и ряды. Элементарные функции.» — М.: Наука, физматгиз, 1984. 21) Фихтенгольц Г. М. «Курс дифференциального и интегрального исчисления.» Т. 2. — М.: Гос. изд-во техникотеоретич Градштейн И. С., Рыжик Н. М. «Таблицы интегралов, сумм, рядов и произведений.» — М.: Физматгиз, 1962. 9. Снеддон И. Преобразование Фурье. — М.: Изд-во ин. лит-ры, 1955.еской литры, 1948. 22) Fb.ru/Статья «Ряды Фурье: история и влияние математического механизма на развитие науки» 23) Wikipedia «Быстое преобразование Фурье», 10/04/2018 24) Wikipedia «Ряд Фурье», 10/04/2018 25) Ets.imo.ru «Быстрое преобразование Фурье» 20/03/2018 26) Youtube «But what is the Fourier transform? A visual introduction» 11/02/2018 27) Youtube «La transformation de Fourier» 11/02/2018

Приложение

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

Рис. А Рис. В Рис. С Рис. D

Рис. E

Рецензии:

Авторизуйтесь, чтобы добавить рецензию

- у работы пока нет рецензий -

Отзывы:

Авторизуйтесь, чтобы оставить отзыв