Ряд Фурье

редактировать
Разложение периодических функций на суммы более простых синусоидальных форм
Преобразование Фурье
Непрерывное преобразование Фурье
Ряд Фурье
Дискретное преобразование Фурье
Дискретное преобразование Фурье
Дискретное преобразование Фурье по кольцу
Анализ Фурье
Связанные преобразования

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

Функция s (x) {\ displaystyle s (x)}s(x)(красным цветом) представляет собой сумму шести синусоидальных функций разных амплитуд и гармонически связанных частот. Их суммирование называется рядом Фурье. Преобразование Фурье, S (f) {\ displaystyle S (f)}S(f)(синим цветом), которое отображает зависимость амплитуды от частоты, показывает 6 частот (в нечетных гармониках) и их амплитуды (1 / нечетное число).
Содержание
  • 1 История
  • 2 Определение
    • 2.1 Комплексные функции
    • 2.2 Другие общепринятые обозначения
  • 3 Сходимость
  • 4 Примеры
    • 4.1 Пример 1 : простой ряд Фурье
    • 4.2 Пример 2: Мотивация Фурье
    • 4.3 Другие приложения
  • 5 Начало
    • 5.1 Рождение гармонического анализа
  • 6 Расширения
    • 6.1 Ряд Фурье на квадрате
    • 6.2 Ряд Фурье решеточно-периодической функции Браве
    • 6.3 Интерпретация гильбертова пространства
  • 7 Свойства
    • 7.1 Таблица основных свойств
    • 7.2 Свойства симметрии
    • 7.3 Лемма Римана – Лебега
    • 7.4 Производная свойство
    • 7.5 Теорема Парсеваля
    • 7.6 Теорема Планшереля
    • 7.7 Теоремы о свертке
    • 7.8 Компактные группы
    • 7.9 Римановы многообразия
    • 7.10 Локально компактные абелевы группы
  • 8 Таблица общий ряд Фурье
  • 9 Аппроксимация и сходимость ряда Фурье
    • 9.1 Свойство наименьших квадратов
    • 9.2 Сходимость
    • 9.3 Дивергенция
  • 10 См. также
  • 11 Примечания
  • 12 Ссылки
    • 12.1 Дополнительная литература
  • 13 Внешние ссылки
История

Ряд Фурье назван в честь Жана-Батиста Жозефа Фурье (1768–1830), внесшего важный вклад в исследование из тригонометрического ряда, после предварительных исследований Леонарда Эйлера, Жана ле Ронд д'Аламбера и Даниэля Бернулли. Фурье ввел ряд с целью решения уравнения теплопроводности в металлической пластине, опубликовав свои первые результаты в своей книге 1807 Mémoire sur la distribution de la chaleur dans les corps solides (Трактат о распространение тепла в твердых телах) и опубликовал его Théorie analytique de la chaleur (Аналитическая теория тепла) в 1822 году. Mémoire ввел анализ Фурье, в частности, ряды Фурье. Исследованиями Фурье был установлен факт, что произвольная (непрерывная) функция может быть представлена ​​тригонометрическим рядом. Первое сообщение об этом великом открытии было сделано Фурье в 1807 году перед Французской академией. Ранние идеи разложения периодической функции на сумму простых осциллирующих функций восходят к 3 веку до нашей эры, когда древние астрономы предложили эмпирическую модель движения планет, основанную на деферентах и ​​эпициклах.

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

С современной точки зрения результаты Фурье несколько неформальны из-за отсутствия точного понятия функции и интеграла в начале девятнадцатого века. Позже Петер Густав Лежен Дирихле и Бернхард Риман выразили результаты Фурье с большей точностью и формальностью.

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

Определение

Рассмотрим функцию с действительным знаком, s (x) {\ displaystyle s (x)}s(x), то есть интегрируемую на интервале длины P {\ displaystyle P}P, который будет периодом ряда Фурье. Типичные примеры интервалов анализа:

x ∈ [0, 1], {\ displaystyle x \ in [0,1],}{\ displaystyle x \ in [0,1],} и P = 1. {\ displaystyle P = 1.}{\displaystyle P=1.}
x ∈ [- π, π], {\ displaystyle x \ in [- \ pi, \ pi],}{\displaystyle x\in [-\pi,\pi ],}и P = 2 π. {\ displaystyle P = 2 \ pi.}{\ displaystyle P = 2 \ pi.}

Процесс анализа определяет веса, индексированные целым числом n {\ displaystyle n}n, которое также является числом циклов гармоники n-й {\ displaystyle n ^ {\ text {th}}}n^{{\text{th}}}в интервале анализа. Следовательно, длина цикла в единицах x {\ displaystyle x}xравна P / n {\ displaystyle P / n}{\displaystyle P/n}. И соответствующая частота гармоники равна n / P {\ displaystyle n / P}n/P. n-я {\ displaystyle n ^ {th}}n^{th}гармоники: sin ⁡ (2 π xn P) {\ displaystyle \ sin \ left (2 \ pi x {\ tfrac {n } {P}} \ right)}{\displaystyle \sin \left(2\pi x{\tfrac {n}{P}}\right)}и cos ⁡ (2 π xn P) {\ displaystyle \ cos \ left (2 \ pi x {\ tfrac {n} {P}} \ right)}{\displaystyle \cos \left(2\pi x{\tfrac {n}{P}}\right)}, а их амплитуды (веса) находятся путем интегрирования в интервале длины P {\ displaystyle P}P:

коэффициенты Фурье
an = 2 P ∫ P s ( x) cos ⁡ (2 π xn P) dxbn = 2 P ∫ P s (x) ⋅ sin ⁡ (2 π xn P) dx. {\ displaystyle {\ begin {align} a_ {n} = {\ frac {2} {P}} \ int _ {P} s (x) \ cdot \ cos \ left (2 \ pi x {\ tfrac { n} {P}} \ right) \, dx \\ b_ {n} = {\ frac {2} {P}} \ int _ {P} s (x) \ cdot \ sin \ left (2 \ pi x {\ tfrac {n} {P}} \ right) \, dx. \ end {align}}}{\displaystyle {\begin{aligned}a_{n}={\frac {2}{P}}\int _{P}s(x)\cdot \cos \left(2\pi x{\tfrac {n}{P}}\right)\,dx\\b_{n}={\frac {2}{P}}\int _{P}s(x)\cdot \sin \left(2\pi x{\tfrac {n}{P}}\right)\,dx.\end{aligned}}}

(Eq.1)

  • Если s (x) {\ displaystyle s (x)}s(x)является P {\ displaystyle P}P-периодическим, тогда достаточно любого интервала такой длины.
  • a 0 {\ displaystyle a_ {0}}a_{0}и b 0 {\ displaystyle b_ {0}}b_ {0} можно уменьшить до a 0 = 2 P ∫ P s (x) dx {\ displaystyle a_ {0 } = {\ frac {2} {P}} \ int _ {P} s (x) \, dx}{\displaystyle a_{0}={\frac {2}{P}}\int _{P}s(x)\,dx}и b 0 = 0 {\ displaystyle b_ {0} = 0}b_ {0} = 0 .
  • Многие тексты выбирают P = 2 π {\ displaystyle P = 2 \ pi}{\displaystyle P=2\pi }, чтобы упростить аргумент синусоидальных функций.

Процесс синтеза ( фактический ряд Фурье):

ряд Фурье, синус-косинусная форма
s N (x) = a 0 2 + ∑ n = 1 N (an cos ⁡ (2 π nx P) + bn sin ⁡ (2 π nx P)). {\ displaystyle {\ begin {align} s_ {N} (x) = {\ frac {a_ {0}} {2}} + \ sum _ {n = 1} ^ {N} \ left (a_ {n} \ cos \ left ({\ tfrac {2 \ pi nx} {P}} \ right) + b_ {n} \ sin \ left ({\ tfrac {2 \ pi nx} {P}} \ right) \ right). \ end {align}}}{\displaystyle {\begin{aligned}s_{N}(x)={\frac {a_{0}}{2}}+\sum _{n=1}^{N}\left(a_{n}\cos \left({\tfrac {2\pi nx}{P}}\right)+b_{n}\sin \left({\tfrac {2\pi nx}{P}}\right)\right).\end{aligned}}}

(Eq.2)

В общем, целое число N {\ displaystyle N}Nтеоретически бесконечно. Даже в этом случае ряд может не сходиться или точно соответствовать s (x) {\ displaystyle s (x)}s(x)при всех значениях x {\ displaystyle x}x(например, одноточечный разрыв) в интервале анализа. Для «хороших» функций, типичных для физических процессов, обычно предполагается равенство.

Если s (t) {\ displaystyle s (t)}s(t)- функция, содержащаяся в интервале длины P {\ displaystyle P}P( и ноль в другом месте), верхний правый квадрант является примером того, как могут выглядеть его коэффициенты ряда Фурье (A n {\ displaystyle A_ {n}}A_{n}), если их сопоставить с соответствующими частотами гармоник.. Левый верхний квадрант - это соответствующее преобразование Фурье s (t). {\ displaystyle s (t).}s(t).Суммирование ряда Фурье (не показано) синтезирует периодическое суммирование s (t), {\ displaystyle s (t),}s(t),, тогда как обратное преобразование Фурье (не показано) синтезирует только s (t). {\ displaystyle s (t).}s(t).

Использование тригонометрического тождества:

A n ⋅ cos ⁡ (2 π nx P - φ n) ≡ A n cos ⁡ (φ n) ⏟ an ⋅ cos ⁡ (2 π Nx P) + A N грех ⁡ (φ N) ⏟ BN ⋅ грех ⁡ (2 π Nx P), {\ displaystyle A_ {n} \ cdot \ cos \ left ({\ tfrac {2 \ pi nx} {P }} - \ varphi _ {n} \ right) \ \ Equiv \ \ underbrace {A_ {n} \ cos (\ varphi _ {n})} _ {a_ {n}} \ cdot \ cos \ left ({\ tfrac {2 \ pi nx} {P}} \ right) + \ underbrace {A_ {n} \ sin (\ varphi _ {n})} _ {b_ {n}} \ cdot \ sin \ left ({\ tfrac {2 \ pi nx} {P}} \ right),}{\displaystyle A_{n}\cdot \cos \left({\tfrac {2\pi nx}{P}}-\varphi _{n}\right)\ \equiv \ \underbrace {A_{n}\cos(\varphi _{n})} _{a_{n}}\cdot \cos \left({\tfrac {2\pi nx}{P}}\right)+\underbrace {A_{n}\sin(\varphi _{n})} _{b_{n}}\cdot \sin \left({\tfrac {2\pi nx}{P}}\right),}

и определения A n ≜ an 2 + bn 2 {\ displaystyle A_ {n} \ triangleq {\ sqrt {a_ {n} ^ { 2} + b_ {n} ^ {2}}}}{\displaystyle A_{n}\triangleq {\sqrt {a_{n}^{2}+b_{n}^{2}}}}и φ n ≜ arctan2 ⁡ (bn, an) {\ displaystyle \ varphi _ {n} \ triangleq \ operatorname {arctan2} (b_ {n}, a_ {n})}{\displaystyle \varphi _{n}\triangleq \operatorname {arctan2} (b_{n},a_{n})}, пары синус и косинус могут быть выражены как одна синусоида со смещением фазы, аналогично преобразованию между ортогональными (декартовыми) и полярными координатами:

Ряд Фурье, амплитуда-фаза
s N (x) = A 0 2 + ∑ n = 1 NA n ⋅ cos ⁡ (2 π nx P - φ n). {\ displaystyle s_ {N} (x) = {\ frac {A_ {0}} {2}} + \ sum _ {n = 1} ^ {N} A_ {n} \ cdot \ cos \ left ({\ tfrac {2 \ pi nx} {P}} - \ varphi _ {n} \ right).}{\displaystyle s_{N}(x)={\frac {A_{0}}{2}}+\sum _{n=1}^{N}A_{n}\cdot \cos \left({\tfrac {2\pi nx}{P}}-\varphi _{n}\right).}

(уравнение 3)

Обычная форма для обобщения на комплексные значения s (x) {\ displaystyle s (x)}s(x)(следующий раздел) получается с использованием формулы Эйлера для разделения функции косинуса на комплексные экспоненты. Здесь комплексное сопряжение обозначено звездочкой:

cos ⁡ (2 π nx P - φ n) ≡ 1 2 ei (2 π nx P - φ n) + 1 2 e - i ( 2 π nx P - φ n) знак равно (1 2 e - i φ n) ⋅ ei 2 π (+ n) x P + (1 2 e - i φ n) ∗ ei 2 π (- n) x P. {\ displaystyle {\ begin {array} {lll} \ cos \ left ({\ tfrac {2 \ pi nx} {P}} - \ varphi _ {n} \ right) {} \ Equiv {\ tfrac {1 } {2}} e ^ {i \ left ({\ tfrac {2 \ pi nx} {P}} - \ varphi _ {n} \ right)} {} + {\ tfrac {1} {2}} e ^ {- i \ left ({\ tfrac {2 \ pi nx} {P}} - \ varphi _ {n} \ right)} \\ = \ left ({\ tfrac {1} {2}} e ^ {- i \ varphi _ {n}} \ right) \ cdot e ^ {i {\ tfrac {2 \ pi (+ n) x} {P}}} {} + \ left ({\ tfrac {1 } {2}} e ^ {- i \ varphi _ {n}} \ right) ^ {*} \ cdot e ^ {i {\ tfrac {2 \ pi (-n) x} {P}}}. \ end {array}}}{\displaystyle {\begin{array}{lll}\cos \left({\tfrac {2\pi nx}{P}}-\varphi _{n}\right){}\equiv {\tfrac {1}{2}}e^{i\left({\tfrac {2\pi nx}{P}}-\varphi _{n}\right)}{}+{\tfrac {1}{2}}e^{-i\left({\tfrac {2\pi nx}{P}}-\varphi _{n}\right)}\\=\left({\tfrac {1}{2}}e^{-i\varphi _{n}}\right)\cdot e^{i{\tfrac {2\pi (+n)x}{P}}}{}+\left({\tfrac {1}{2}}e^{-i\varphi _{n}}\right)^{*}\cdot e^{i{\tfrac {2\pi (-n)x}{P}}}.\end{array}}}

Следовательно, с определениями:

cn ≜ {A 0/2 = a 0/2, n = 0 A n 2 e - i φ n = 1 2 (an - ibn), n>0 c | п | ∗, n < 0 } = 1 P ∫ P s ( x) ⋅ e − i 2 π n x P d x, {\displaystyle c_{n}\triangleq \left\{{\begin{array}{lll}A_{0}/2=a_{0}/2,\quad n=0\\{\tfrac {A_{n}}{2}}e^{-i\varphi _{n}}={\tfrac {1}{2}}(a_{n}-ib_{n}),\quad n>0 \\ c_ {| n |} ^ {*}, \ quad n <0\end{array}}\right\}\quad =\quad {\frac {1}{P}}\int _{P}s(x)\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx,}{\displaystyle c_{n}\triangleq \left\{{\begin{array}{lll}A_{0}/2=a_{0}/2,\quad n=0\\{\tfrac {A_{n}}{2}}e^{-i\varphi _{n}}={\tfrac {1}{2}}(a_{n}-ib_{n}),\quad n>0 \\ c_ {| n |} ^ {*}, \ quad n <0\end{array}}\right\}\quad =\quad {\frac {1}{P}}\int _{P}s(x)\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx,}

окончательный результат :

Ряд Фурье, экспоненциальная форма
s N (x) = ∑ n = - NN cn ⋅ ei 2 π nx P. {\ displaystyle s_ {N} (x) = \ sum _ {n = -N} ^ {N} c_ {n} \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P}}}.}{\ displaystyle s_ {N} (x) = \ sum _ {n = -N} ^ {N} c_ {n} \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P}}}.}

(Eq.4)

Комплексные функции

Если s (x) {\ displaystyle s (x)}s(x)является комплексной функцией от вещественная переменная x, {\ displaystyle x,}x, оба компонента (действительная и мнимая части) являются функциями с действительными значениями, которые могут быть представлены рядом Фурье. Два набора коэффициентов и частичная сумма определяются по формуле:

c R n = 1 P ∫ P Re ⁡ {s (x)} ⋅ e - i 2 π nx P dx {\ displaystyle c _ {_ {Rn} } = {\ frac {1} {P}} \ int _ {P} \ operatorname {Re} \ {s (x) \} \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P} }} \ dx}{\displaystyle c_{_{Rn}}={\frac {1}{P}}\int _{P}\operatorname {Re} \{s(x)\}\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx}и c I n = 1 P ∫ P Im ⁡ {s (x)} ⋅ e - i 2 π nx P dx {\ displaystyle c _ {_ {In}} = {\ frac {1} {P}} \ int _ {P} \ operatorname {Im} \ {s (x) \} \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P}} } \ dx}{\displaystyle c_{_{In}}={\frac {1}{P}}\int _{P}\operatorname {Im} \{s(x)\}\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx}
s N (x) = ∑ n = - NN c R n ⋅ ei 2 π nx P + i ⋅ ∑ n = - NN c I n ⋅ ei 2 π nx P = ∑ n = - NN (c R n + i ⋅ c I n) ⋅ ei 2 π nx P. {\ displaystyle s_ {N} (x) = \ sum _ {n = -N} ^ {N} c _ {_ {Rn}} \ cdot e ^ {я {\ tfrac {2 \ pi nx} {P}} } + i \ cdot \ sum _ {n = -N} ^ {N} c _ {_ {In}} \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P}}} = \ sum _ { n = -N} ^ {N} \ left (c _ {_ {Rn}} + i \ cdot c _ {_ {In}} \ right) \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P }}}.}{\displaystyle s_{N}(x)=\sum _{n=-N}^{N}c_{_{Rn}}\cdot e^{i{\tfrac {2\pi nx}{P}}}+i \cdot \sum _{n=-N}^{N}c_{_{In}}\cdot e^{i{\tfrac {2\pi nx}{P}}}=\sum _{n=- N}^{N}\left(c_{_{Rn}}+i\cdot c_{_{In}}\right)\cdot e^{i{\tfrac {2\pi nx}{P}}}.}

Определение cn ≜ c R n + i ⋅ c I n {\ displaystyle c_ {n} \ triangleq c _ {_ {Rn}} + i \ cdot c _ {_ {In}} }{\displaystyle c_{n}\triangleq c_{_{Rn}}+i\cdot c_{_{In}}}возвращает:

s N (x) = ∑ n = - NN cn ⋅ ei 2 π nx P. {\ displaystyle s_ {N} (x) = \ sum _ {n = -N} ^ {N} c_ {n} \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P}}}.}{\ displaystyle s_ {N} (x) = \ sum _ {n = -N} ^ {N} c_ {n} \ cdot e ^ {i {\ tfrac {2 \ pi nx} {P}}}.}

(уравнение 5)

Это идентично уравнению 4, за исключением cn {\ displaystyle c_ {n}}c_{n}и c - n {\ displaystyle c _ {- n}}{\displaystyle c_{-n}}больше не являются комплексными конъюгатами. Формула для cn {\ displaystyle c_ {n}}c_{n}также не изменилась:

cn = 1 P ∫ P Re ⁡ {s (x)} ⋅ e - i 2 π nx P dx + i ⋅ 1 P ∫ P Im ⁡ {s (x)} ⋅ e - i 2 π nx P dx = 1 P ∫ P (Re ⁡ {s (x)} + i ⋅ Im ⁡ {s (x)}) ⋅ e - я 2 π nx P dx знак равно 1 P ∫ P s (x) ⋅ e - i 2 π nx P dx. {\ displaystyle {\ begin {align} c_ {n} = {\ frac {1} {P}} \ int _ {P} \ operatorname {Re} \ {s (x) \} \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P}}} \ dx + i \ cdot {\ frac {1} {P}} \ int _ {P} \ operatorname {Im} \ {s (x) \} \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P}}} \ dx \\ [4pt] = {\ frac {1} {P}} \ int _ {P} \ left (\ OperatorName {Re} \ {s (x) \} + i \ cdot \ operatorname {Im} \ {s (x) \} \ right) \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P }}} \ dx \ = \ {\ frac {1} {P}} \ int _ {P} s (x) \ cdot e ^ {- i {\ tfrac {2 \ pi nx} {P}}} \ dx. \ end {align}}}{\displaystyle {\begin{aligned}c_{n}={\frac {1}{P}}\int _{P}\operatorname {Re} \{s(x)\}\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx+i\cdot {\frac {1}{P}}\int _{P}\operatorname {Im} \{s(x)\}\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx\\[4pt]={\frac {1}{P}}\int _{P}\left(\operatorname {Re} \{s(x)\}+i\cdot \operatorname {Im} \{s(x)\}\right)\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx\ =\ {\frac {1}{P}}\int _{P}s(x)\cdot e^{-i{\tfrac {2\pi nx}{P}}}\ dx.\end{aligned}}}

Другие распространенные обозначения

Обозначение cn {\ displaystyle c_ {n}}c_{n}неадекватно для обсуждения коэффициентов Фурье нескольких разные функции. Поэтому его обычно заменяют модифицированной формой функции (в данном случае s {\ displaystyle s}s ), например s ^ (n) {\ displaystyle {\ hat {s}} (n)}{\displaystyle {\hat {s}}(n)}или S [n] {\ displaystyle S [n]}{\displaystyle S[n]}, а функциональная нотация часто заменяет индекс:

s ∞ (Икс) знак равно ∑ N = - ∞ ∞ s ^ (n) ⋅ ei 2 π nx / P = ∑ n = - ∞ ∞ S [n] ⋅ ej 2 π nx / P commonengineeringnotation {\ displaystyle {\ begin {выровнено } s _ {\ infty} (x) = \ sum _ {n = - \ infty} ^ {\ infty} {\ hat {s}} (n) \ cdot e ^ {i \, 2 \ pi nx / P } \\ [6pt] = \ sum _ {n = - \ infty} ^ {\ infty} S [n] \ cdot e ^ {j \, 2 \ pi nx / P} \ scriptstyle {\ mathsf {common \ engineering \ notation}} \ end {align}}}{\displaystyle {\begin{aligned}s_{\infty }(x)=\sum _{n=-\infty }^{\infty }{\hat {s}}(n)\cdot e^{i\,2\pi nx/P}\\[6pt]=\sum _{n=-\infty }^{\infty }S[n]\cdot e^{j\,2\pi nx/P}\scriptstyle {\mathsf {common\ engineering\ notation}}\end{aligned}}}

В инженерии, особенно когда переменная x {\ displaystyle x}xпредставляет время, последовательность коэффициентов называется представление в частотной области. Квадратные скобки часто используются, чтобы подчеркнуть, что область действия этой функции представляет собой дискретный набор частот.

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

S (f) ≜ ∑ n = - ∞ ∞ S [n] ⋅ δ (f - n P), {\ Displaystyle S (е) \ \ треугольник \ \ сумма _ {n = - \ infty} ^ {\ infty} S [n] \ cdot \ delta \ left (f - {\ frac {n} {P} } \ right),}{\displaystyle S(f)\ \triangleq \ \sum _{n=-\infty }^{\infty }S[n]\cdot \delta \left(f-{\frac {n}{P}}\right),}

где f {\ displaystyle f}fпредставляет непрерывную частотную область. Когда переменная x {\ displaystyle x}xимеет единицы измерения в секундах, f {\ displaystyle f}fимеет единицы измерения герц. «Зубцы» гребенки разнесены на интервалы, кратные гармоникам ) 1 / P {\ displaystyle 1 / P}1/P, что называется основной частота. s ∞ (x) {\ displaystyle s _ {\ infty} (x)}{\displaystyle s_{\infty }(x)}можно восстановить из этого представления с помощью обратного преобразования Фурье :

F - 1 {S (f)} = ∫ - ∞ ∞ (∑ n = - ∞ ∞ S [n] ⋅ δ (f - n P)) ei 2 π fxdf, = ∑ n = - ∞ ∞ S [n] ⋅ ∫ - ∞ ∞ δ ( f - n P) ei 2 π fxdf, = ∑ n = - ∞ ∞ S [n] ⋅ ei 2 π nx / P s ∞ (x). {\ displaystyle {\ begin {align} {\ mathcal {F}} ^ {- 1} \ {S (f) \} = \ int _ {- \ infty} ^ {\ infty} \ left (\ sum _ {n = - \ infty} ^ {\ infty} S [n] \ cdot \ delta \ left (f - {\ frac {n} {P}} \ right) \ right) e ^ {i2 \ pi fx} \, df, \\ [6pt] = \ sum _ {n = - \ infty} ^ {\ infty} S [n] \ cdot \ int _ {- \ infty} ^ {\ infty} \ delta \ left (f - {\ frac {n} {P}} \ right) e ^ {i2 \ pi fx} \, df, \\ [6pt] = \ sum _ {n = - \ infty} ^ {\ infty} S [ n] \ cdot e ^ {i \, 2 \ pi nx / P} \ \ \ triangleq \ s _ {\ infty} (x). \ end {align}}}{\ displaystyle {\ begin {align} {\ mathcal {F}} ^ {-1} \ {S (f) \} = \ int _ {- \ infty} ^ {\ infty} \ left (\ sum _ {n = - \ infty} ^ {\ infty} S [n] \ cdot \ delta \ left (f - {\ frac {n} {P}} \ right) \ right) e ^ {i2 \ pi fx} \, df, \\ [6pt] = \ sum _ {n = - \ infty} ^ {\ infty} S [n ] \ cdot \ int _ {- \ infty} ^ {\ infty} \ delta \ left (f - {\ frac {n} {P}} \ right) e ^ {i2 \ pi fx} \, df, \\ [6pt] = \ sum _ {n = - \ infty} ^ {\ infty} S [n] \ cdot e ^ {i \, 2 \ pi nx / P} \ \ \ triangleq \ s _ {\ infty} ( x). \ end {align}}}

Построенная функция S (f) {\ displaystyle S (f)}S(f)поэтому обычно называют преобразованием Фурье, хотя интеграл Фурье периодической функции не сходится на частотах гармоник.

Сходимость

В инженерных приложениях обычно предполагается, что ряды Фурье сходятся почти везде (исключения находящиеся на дискретных разрывах), поскольку функции, встречающиеся в инженерии, ведут себя лучше чем функции, которые математики могут предоставить в качестве контрпримеров этому предположению. В частности, если s {\ displaystyle s}s является непрерывным и производная от s (x) {\ displaystyle s (x)}s(x)(что может не существуют везде) интегрируем в квадрате, то ряд Фурье s {\ displaystyle s}s абсолютно и равномерно сходится к s (x) {\ displaystyle s (x)}s(x). Если функция интегрируема с квадратом на интервале [x 0, x 0 + P] {\ displaystyle [x_ {0}, x_ {0} + P]}{\displaystyle [x_{0},x_{0}+P]}, то ряд Фурье сходится к функции почти в каждой точке. Сходимость рядов Фурье также зависит от конечного числа максимумов и минимумов функции, которая широко известна как одно из условий Дирихле для рядов Фурье. См. Сходимость рядов Фурье. Можно определить коэффициенты Фурье для более общих функций или распределений, в таких случаях обычно представляет интерес сходимость по норме или слабая сходимость.

Здесь можно увидеть интерактивную анимацию .

Примеры

Пример 1: простой ряд Фурье

График пилообразной волны, периодического продолжения линейная функция s (x) = x / π {\ displaystyle s (x) = x / \ pi}s(x)=x/\pi на интервале (- π, π] {\ displaystyle (- \ pi, \ pi]}(-\pi,\pi ] Анимированный график первых пяти последовательных частных рядов Фурье

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

s ( x) = x π, для - π < x < π, {\displaystyle s(x)={\frac {x}{\pi }},\quad \mathrm {for} -\pi s(x) = \frac{x}{\pi}, \quad \mathrm{for } -\pi <x <\pi,
s (x + 2 π k) = s (x), для - π < x < π and k ∈ Z. {\displaystyle s(x+2\pi k)=s(x),\quad \mathrm {for} -\pi {\displaystyle s(x+2\pi k)=s(x),\quad \mathrm {for} -\pi <x<\pi {\text{ and }}k\in \mathbb {Z}.}

В этом случае коэффициенты Фурье определяются как

an = 1 π ∫ - π π s (x) cos ⁡ (nx) dx = 0, n ≥ 0. bn = 1 π ∫ - π π s (x) sin ⁡ (nx) dx = - 2 π n cos ⁡ (n π) + 2 π 2 N 2 грех ⁡ (N π) = 2 (- 1) n + 1 π n, n ≥ 1. {\ displaystyle {\ begin {align} a_ {n} = {\ frac {1} {\ pi}} \ int _ {- \ pi} ^ {\ pi} s (x) \ cos (nx) \, dx = 0, \ quad n \ geq 0. \\ [4pt] b_ {n} = { \ frac {1} {\ pi}} \ int _ {- \ pi} ^ {\ pi} s (x) \ si n (nx) \, dx \\ [4pt] = - {\ frac {2} {\ pi n}} \ cos (n \ pi) + {\ frac {2} {\ pi ^ {2} n ^ {2}}} \ sin (n \ pi) \\ [4pt] = {\ frac {2 \, (- 1) ^ {n + 1}} {\ pi n}}, \ quad n \ geq 1. \ end {align}}}{\displaystyle {\begin{aligned}a_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\cos(nx)\,dx=0,\quad n\geq 0.\\[4pt]b_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\sin(nx)\,dx\\[4pt]=-{\frac {2}{\pi n}}\cos(n\pi)+{\frac {2}{\pi ^{2}n^{2}}}\sin(n\pi)\\[4pt]={\frac {2\,(-1)^{n+1}}{\pi n}},\quad n\geq 1.\end{aligned}}}

Можно доказать, что ряд Фурье сходится к s (x) {\ displaystyle s (x)}s(x)в каждой точке x {\ displaystyle x}xгде s {\ displaystyle s}s дифференцируемый, и, следовательно:

s (x) = a 0 2 + ∑ n = 1 ∞ [an cos ⁡ (nx) + bn sin ⁡ (nx)] = 2 π ∑ n = 1 ∞ (- 1) n + 1 n sin ⁡ (nx), для x - π ∉ 2 π Z. {\ displaystyle {\ begin {align} s (x) = {\ frac {a_ {0}} {2}} + \ sum _ {n = 1} ^ {\ infty} \ left [a_ {n} \ cos \ left (nx \ right) + b_ {n} \ sin \ left (nx \ right) \ right] \\ [4pt] = {\ frac {2} {\ pi}} \ sum _ {n = 1 } ^ {\ infty} {\ frac {(-1) ^ {n + 1}} {n}} \ sin (nx), \ quad \ mathrm {for} \ quad x- \ pi \ notin 2 \ pi \ mathbb {Z}. \ end {align}}}{\displaystyle {\begin{aligned}s(x)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[a_{n}\cos \left(nx\right)+b_{n}\sin \left(nx\right)\right]\\[4pt]={\frac {2}{\pi }}\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx),\quad \mathrm {for} \quad x-\pi \notin 2\pi \mathbb {Z}.\end{aligned}}}

(Eq.7)

Когда x = π {\ displaystyle x = \ pi}x=\pi, ряд Фурье сходится к 0, который представляет собой полусумму левого и правого пределов s в x = π {\ displaystyle x = \ pi}x=\pi. Это частный пример теоремы Дирихле для рядов Фурье.

Распределение тепла в металлической пластине с использованием метода Фурье

Этот пример приводит нас к решению проблемы Базеля.

Пример 2: Мотивация Фурье

Разложение в ряд Фурье нашего функция в примере 1 выглядит более сложной, чем простая формула s (x) = x / π {\ displaystyle s (x) = x / \ pi}{\displaystyle s(x)=x/\pi }, поэтому не сразу понятно, почему потребуется ряд Фурье. Хотя существует множество приложений, мотивация Фурье заключалась в решении уравнения теплопроводности. Например, рассмотрим металлическую пластину в форме квадрата со стороной π {\ displaystyle \ pi}\ pi метров с координатами (x, y) ∈ [0, π] × [0, π] {\ displaystyle (x, y) \ in [0, \ pi] \ times [0, \ pi]}{\displaystyle (x,y)\in [0,\pi ]\times [0,\pi ]}. Если внутри пластины нет источника тепла и если температура на трех из четырех сторон составляет 0 градусов Цельсия, а четвертая сторона, заданная как y = π {\ displaystyle y = \ pi}{\displaystyle y=\pi }поддерживается на уровне температурного градиента T (x, π) = x {\ displaystyle T (x, \ pi) = x}{\displaystyle T(x,\pi)=x}градусов Цельсия для x {\ displaystyle x }xв (0, π) {\ displaystyle (0, \ pi)}{\displaystyle (0,\pi)}, тогда можно показать, что стационарное распределение тепла (или распределение тепла после длительного время истекло) определяется как

T (x, y) = 2 ∑ n = 1 ∞ (- 1) n + 1 n sin ⁡ (nx) sinh ⁡ (ny) sinh ⁡ (n π). {\ Displaystyle Т (х, у) = 2 \ сумма _ {п = 1} ^ {\ infty} {\ гидроразрыва {(-1) ^ {п + 1}} {п}} \ грех (пх) {\ sinh (ny) \ over \ sinh (n \ pi)}.}T ( x, y) = 2 \ sum_ {n = 1} ^ \ infty \ frac {(- 1) ^ {n + 1}} {n} \ sin (nx) {\ sinh (ny) \ over \ sinh (n \ pi)}.

Здесь sinh - это функция гиперболического синуса. Это решение уравнения теплопроводности получается путем умножения каждого члена Eq.7на sinh ⁡ (ny) / sinh ⁡ (n π) {\ displaystyle \ sinh ( ny) / \ sinh (n \ pi)}{\displaystyle \sinh(ny)/\sinh(n\pi)}. Хотя в нашем примере функция s (x) {\ displaystyle s (x)}s(x), кажется, имеет излишне сложный ряд Фурье, распределение тепла T (x, y) {\ displaystyle T (x, y)}{\ displaystyle T (x, y)} нетривиально. Функция T {\ displaystyle T}T не может быть записана как выражение в закрытой форме. Этот метод решения тепловой проблемы стал возможным благодаря работе Фурье.

Другие приложения

Еще одно применение этого ряда Фурье - решение проблемы Базеля с использованием теоремы Парсеваля. Пример обобщает, и можно вычислить ζ (2n) для любого положительного целого числа n.

Начало

Жозеф Фурье писал:

φ (y) = a 0 cos ⁡ π y 2 + a 1 cos ⁡ 3 π y 2 + a 2 cos ⁡ 5 π y 2 + ⋯. {\ displaystyle \ varphi (y) = a_ {0} \ cos {\ frac {\ pi y} {2}} + a_ {1} \ cos 3 {\ frac {\ pi y} {2}} + a_ { 2} \ cos 5 {\ frac {\ pi y} {2}} + \ cdots.}\varphi(y)=a_0\cos\frac{\pi y}{2}+a_1\cos 3\frac{\pi y}{2}+a_2\cos5\frac{\pi y}{2}+\cdots.

Умножение обеих сторон на cos ⁡ (2 k + 1) π y 2 {\ displaystyle \ cos (2k +1) {\ frac {\ pi y} {2}}}\cos(2k+1)\frac{\pi y}{2}, а затем интегрирование от y = - 1 {\ displaystyle y = -1}y=-1до y = + 1 {\ displaystyle y = + 1}y=+1дает:

ak = ∫ - 1 1 φ (y) cos ⁡ (2 k + 1) π y 2 dy. {\ displaystyle a_ {k} = \ int _ {- 1} ^ {1} \ varphi (y) \ cos (2k + 1) {\ frac {\ pi y} {2}} \, dy.}a_k=\int_{-1}^1\varphi(y)\cos(2k+1)\frac{\pi y}{2}\,dy.

— Жозеф Фурье, Память о пропаганде шаллера в солидном корпусе. (1807)

Это немедленно дает любой коэффициент a k из тригонометрического ряда для φ (y) для любой функции, которая имеет такое расширение. Это работает, потому что если φ имеет такое разложение, то (при подходящих условиях сходимости) интеграл

ak = ∫ - 1 1 φ (y) cos ⁡ (2 k + 1) π y 2 dy = ∫ - 1 1 ( a соз ⁡ π y 2 соз ⁡ (2 k + 1) π y 2 + a ′ cos ⁡ 3 π y 2 cos ⁡ (2 k + 1) π y 2 + ⋯) dy {\ displaystyle {\ begin {align} a_ {k} = \ int _ {- 1} ^ {1} \ varphi (y) \ cos (2k + 1) {\ frac {\ pi y} {2}} \, dy \\ = \ int _ {- 1} ^ {1} \ left (a \ cos {\ frac {\ pi y} {2}} \ cos (2k + 1) {\ frac {\ pi y} {2}} + a '\ cos 3 {\ frac {\ pi y} {2}} \ cos (2k + 1) {\ frac {\ pi y} {2}} + \ cdots \ right) \, dy \ end {align}}}\begin{align} a_k=\int_{-1}^1\varphi(y)\cos(2k+1)\frac{\pi y}{2}\,dy \\ = \int_{-1}^1\left(a\cos\frac{\pi y}{2}\cos(2k+1)\frac{\pi y}{2}+a'\cos 3\frac{\pi y}{2}\cos(2k+1)\frac{\pi y}{2}+\cdots\right)\,dy \end{align}

можно проводить посменно. Но все члены, содержащие cos ⁡ (2 j + 1) π y 2 cos ⁡ (2 k + 1) π y 2 {\ displaystyle \ cos (2j + 1) {\ frac {\ pi y} {2} } \ cos (2k + 1) {\ frac {\ pi y} {2}}}\cos(2j+1)\frac{\pi y}{2} \cos(2k+1)\frac{\pi y}{2}для j ≠ k исчезают при интегрировании от -1 до 1, оставляя только k-й член.

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

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

Рождение гармонического анализа

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

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

Расширения

Ряд Фурье на квадрате

Мы также можем определить ряд Фурье для функций двух переменных x {\ displaystyle x}xи y {\ displaystyle y}yв квадрате [- π, π] × [- π, π] {\ displaystyle [- \ pi, \ pi] \ times [ - \ pi, \ pi]}{\displaystyle [-\pi,\pi ]\times [-\pi,\pi ]}:

f (x, y) = ∑ j, k ∈ Z (целые числа) cj, keijxeiky, cj, k = 1 4 π 2 ∫ - π π ∫ - π π f (x, y) e - ijxe - ikydxdy. {\ Displaystyle {\ begin {выровнено} е (х, у) = \ сумма _ {j, к \, \ in \, \ mathbb {Z} {\ text {(целые числа)}}} c_ {j, k } e ^ {ijx} e ^ {iky}, \\ [5pt] c_ {j, k} = {1 \ over 4 \ pi ^ {2}} \ int _ {- \ pi} ^ {\ pi} \ int _ {- \ pi} ^ {\ pi} f (x, y) e ^ {- ijx} e ^ {- iky} \, dx \, dy. \ end {align}}}{\displaystyle {\begin{aligned}f(x,y)=\sum _{j,k\,\in \,\mathbb {Z} {\text{ (integers)}}}c_{j,k}e^{ijx}e^{iky},\\[5pt]c_{j,k}={1 \over 4\pi ^{2}}\int _{-\pi }^{\pi }\int _{-\pi }^{\pi }f(x,y)e^{-ijx}e^{-iky}\,dx\,dy.\end{aligned}}}

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

Ряд Фурье периодической-решеточной-функции Браве

Трехмерная решетка Браве определяется как набор векторов вида:

R = n 1 a 1 + n 2 a 2 + n 3 a 3 {\ displaystyle \ mathbf {R} = n_ {1} \ mathbf {a} _ {1} + n_ {2} \ mathbf {a} _ {2} + n_ {3} \ mathbf {a} _ {3}}{\displaystyle \mathbf {R} =n_{1}\mathbf {a} _{1}+n_{2}\mathbf {a} _{2}+n_{3}\mathbf {a} _{3}}

где ni {\ displaystyle n_ {i}}n_{i}- целые числа, а ai {\ displaystyle \ mathbf { a} _ {i}}\mathbf{a}_i- три линейно независимых вектора. Предположим, у нас есть некоторая функция, f (r) {\ displaystyle f (\ mathbf {r})}f(\mathbf{r}), такая, что она подчиняется следующему условию для любого вектора решетки Браве R: f (г) знак равно е (р + р) {\ Displaystyle \ mathbf {R}: е (\ mathbf {r}) = е (\ mathbf {R} + \ mathbf {r})}{\displaystyle \mathbf {R} :f(\mathbf {r})=f(\mathbf {R} +\mathbf {r})}, мы могли бы составить из него ряд Фурье. Такого рода функцией может быть, например, эффективный потенциал, который один электрон «ощущает» внутри периодического кристалла. Тогда полезно составить ряд Фурье потенциала, применяя теорему Блоха. Во-первых, мы можем записать любой произвольный вектор r {\ displaystyle \ mathbf {r}}\mathbf {r} в системе координат решетки:

r = x 1 a 1 a 1 + x 2 a 2 a 2 + x 3 a 3 a 3, {\ displaystyle \ mathbf {r} = x_ {1} {\ frac {\ mathbf {a} _ {1}} {a_ {1}}} + x_ {2 } {\ frac {\ mathbf {a} _ {2}} {a_ {2}}} + x_ {3} {\ frac {\ mathbf {a} _ {3}} {a_ {3}}},}\mathbf{r} = x_1\frac{\mathbf{a}_{1}}{a_1}+ x_2\frac{\mathbf{a}_{2}}{a_2}+ x_3\frac{\mathbf{a}_{3}}{a_3},

где ai ≜ | а я |. {\ displaystyle a_ {i} \ треугольникq | \ mathbf {a} _ {i} |.}{\displaystyle a_{i}\triangleq |\mathbf {a} _{i}|.}

Таким образом, мы можем определить новую функцию

g (x 1, x 2, x 3) ≜ f ( r) = f (x 1 a 1 a 1 + x 2 a 2 a 2 + x 3 a 3 a 3). {\ displaystyle g (x_ {1}, x_ {2}, x_ {3}) \ треугольник f (\ mathbf {r}) = f \ left (x_ {1} {\ frac {\ mathbf {a} _ { 1}} {a_ {1}}} + x_ {2} {\ frac {\ mathbf {a} _ {2}} {a_ {2}}} + x_ {3} {\ frac {\ mathbf {a} _ {3}} {a_ {3}}} \ right).}{\displaystyle g(x_{1},x_{2},x_{3})\triangleq f(\mathbf {r})=f\left(x_{1}{\frac {\mathbf {a} _{1}}{a_{1}}}+x_{2}{\frac {\mathbf {a} _{2}}{a_{2}}}+x_{3}{\frac {\mathbf {a} _{3}}{a_{3}}}\right).}

Эта новая функция, g (x 1, x 2, x 3) {\ displaystyle g (x_ {1}, x_ { 2}, x_ {3})}g(x_1,x_2,x_3), теперь является функцией трех переменных, каждая из которых имеет периодичность a 1, a 2, a 3 соответственно:

g (x 1, x 2, x 3) = g (x 1 + a 1, x 2, x 3) = g (x 1, x 2 + a 2, x 3) = g (x 1, x 2, x 3 + a 3). {\ displaystyle g (x_ {1}, x_ {2}, x_ {3}) = g (x_ {1} + a_ {1}, x_ {2}, x_ {3}) = g (x_ {1}, x_ {2} + a_ {2}, x_ {3}) = g (x_ {1}, x_ {2}, x_ {3} + a_ {3}).}{\ displaystyle g (x_ {1}, x_ {2}, x_ {3}) = g (x_ {1} + a_ {1}, x_ {2}, x_ {3}) = g (x_ {1}, x_ {2} + a_ {2}, x_ {3 }) = g (x_ {1}, x_ {2}, x_ {3} + a_ {3}).}

Если мы напишем серию для g на интервале [0, a 1 ] для x 1, мы можем определить следующее:

hone (m 1, x 2, x 3) ≜ 1 a 1 ∫ 0 a 1 г (Икс 1, Икс 2, Икс 3) ⋅ е - я 2 π м 1 a 1 Икс 1 dx 1 {\ displaystyle h ^ {\ mathrm {one}} (m_ {1}, x_ {2 }, x_ {3}) \ треугольникq {\ frac {1} {a_ {1}}} \ int _ {0} ^ {a_ {1}} g (x_ {1}, x_ {2}, x_ {3 }) \ cdot e ^ {- i2 \ pi {\ frac {m_ {1}} {a_ {1}}} x_ {1}} \, dx_ {1}}{\displaystyle h^{\mathrm {one} }(m_{1},x_{2},x_{3})\triangleq {\frac {1}{a_{1}}}\int _{0}^{a_{1}}g(x_{1},x_{2},x_{3})\cdot e^{-i2\pi {\frac {m_{1}}{a_{1}}}x_{1}}\,dx_{1}}

И тогда мы можем написать:

g (x 1, x 2, x 3) = ∑ m 1 = - ∞ ∞ точить (m 1, x 2, x 3) ⋅ ei 2 π m 1 a 1 x 1 {\ displaystyle g (x_ {1}, x_ {2}, x_ {3}) = \ sum _ {m_ {1} = - \ infty} ^ {\ infty} h ^ {\ mathrm {one}} (m_ {1}, x_ {2}, x_ {3}) \ cdot e ^ {i2 \ pi {\ frac {m_ {1}} {a_ {1}}} x_ {1}}}g(x_1, x_2, x_3)=\sum_{m_1=-\infty}^\infty h^\mathrm{one}(m_1, x_2, x_3) \cdot e^{i 2\pi \frac{m_1}{a_1} x_1}

Дальнейшее определение:

htwo (m 1, m 2, x 3) ≜ 1 a 2 ∫ 0 a 2 hone (m 1, x 2, x 3) ⋅ e - i 2 π m 2 a 2 x 2 dx 2 = 1 a 2 ∫ 0 a 2 dx 2 1 a 1 ∫ 0 a 1 dx 1 g (x 1, x 2, x 3) ⋅ e - i 2 π (m 1 a 1 x 1 + m 2 a 2 x 2) {\ displaystyle {\ begin {align} h ^ {\ mathrm {two}} (m_ {1}, m_ {2}, x_ {3}) и \ треугольникq {\ frac {1 } {a_ {2}}} \ int _ {0} ^ {a_ {2}} h ^ {\ mathrm {one}} (m_ {1}, x_ {2}, x_ {3}) \ cdot e ^ {-i2 \ pi {\ frac {m_ {2}} {a_ {2}}} x_ {2}} \, dx_ {2} \\ [12pt] = {\ frac {1} {a_ {2} }} \ int _ {0} ^ {a_ {2}} dx_ {2} {\ frac {1} {a_ {1}}} \ int _ {0} ^ {a_ {1}} dx_ {1} g (x_ {1}, x_ {2}, x_ {3}) \ cdot e ^ {- i2 \ pi \ left ({\ frac {m_ {1}} {a_ {1}}} x_ {1} + { \ frac {m_ {2}} {a_ {2}}} x_ {2} \ right)} \ end {align}}}{\ displaystyle {\ begin {align} h ^ {\ mathrm {two}} (m_ {1}, m_ {2}, x_ { 3}) \ треугольникq {\ frac {1} {a_ {2}}} \ int _ {0} ^ {a_ {2}} h ^ {\ mathrm {one}} (m_ {1}, x_ {2 }, x_ {3}) \ cdot e ^ {- i2 \ pi {\ frac {m_ {2}} {a_ {2}}} x_ {2}} \, dx_ {2} \\ [12pt] = {\ frac {1} {a_ {2}}} \ int _ {0} ^ {a_ {2}} dx_ {2} {\ frac {1} {a_ {1}}} \ int _ {0} ^ {a_{1}}dx_{1}g(x_{1},x_{2},x_{3})\cdot e^{-i2\pi \left({\frac {m_{1}}{a_ {1}}}x_{1}+{\frac {m_{2}}{a_{2}}}x_{2}\right)}\end{aligned}}}

Мы можем написать g {\ displaystyle g}gеще раз как:

g (x 1, x 2, x 3) = ∑ m 1 = - ∞ ∞ ∑ m 2 = - ∞ ∞ htwo (m 1, m 2, x 3) ⋅ ei 2 π m 1 a 1 x 1 ⋅ ei 2 π m 2 a 2 x 2 {\ displaystyle g (x_ {1}, x_ {2}, x_ {3}) = \ sum _ {m_ {1} = - \ infty} ^ {\infty }\sum _{m_{2}=-\infty }^{\infty }h^{\mathrm {two} }(m_{1},m_{2},x_{3})\cdot e ^{i2\pi {\frac {m_{1}}{a_{1}}}x_{1}}\cdot e^{i2\pi {\frac {m_{2}}{a_{2}}} x_{2}}}g (x_1, x_2, x_3) = \ sum_ {m_1 = - \ infty} ^ \ infty \ sum_ {m_2 = - \ infty} ^ \ infty h ^ \ mathrm {two} (m_1, m_2, x_3) \ cdot e ^ {i 2 \ pi \ гидроразрыв {m_1} {a_1} x_1} \ cdot e ^ {i 2 \ pi \ frac {m_2} {a_2} x_2}

Finally applying the same for the third coordinate, we define:

hthree ( m 1, m 2, m 3) ≜ 1 a 3 ∫ 0 a 3 htwo ( m 1, m 2, x 3) ⋅ e − i 2 π m 3 a 3 x 3 d x 3 = 1 a 3 ∫ 0 a 3 d x 3 1 a 2 ∫ 0 a 2 d x 2 1 a 1 ∫ 0 a 1 d x 1 g ( x 1, x 2, x 3) ⋅ e − i 2 π ( m 1 a 1 x 1 + m 2 a 2 x 2 + m 3 a 3 x 3) {\displaystyle {\begin{aligned}h^{\mathrm {three} }(m_{1},m_{2},m_{3})\triangleq {\frac {1}{a_{3}}}\int _{0}^{a_{3}}h^{\mathrm {two} }(m_{1},m_{2},x_{3})\cdot e^{-i2\pi {\frac {m_{3}}{a_{3}}}x_{3}}\,dx_{3}\\[12pt]={\frac {1}{a_{3}}}\int _{0}^{a_{3}}dx_{3}{\frac {1}{a_{2}}}\int _{0}^{a_{2}}dx_{2}{\frac {1}{a_{1}}}\int _{0}^{a_{1}}dx_{1}g(x_{1},x_{2},x_{3})\cdot e^{-i2\pi \left({\frac {m_{1}}{a_{1}}}x_{1}+{\frac {m_{2}}{a_{2}}}x_{2}+{\frac {m_{3}}{a_{3}}}x_{3}\right)}\end{aligned}}}{\displaystyle {\begin{aligned}h^{\mathrm {three} }(m_{1},m_{2},m_{3})\triangleq {\frac {1}{a_{3}}}\int _{0}^{a_{3}}h^{\mathrm {two} }(m_{1},m_{2},x_{3})\cdot e^{-i2\pi {\frac {m_{3}}{a_{3}}}x_{3}}\,dx_{3}\\[12pt]={\frac {1}{a_{3}}}\int _{0}^{a_{3}}dx_{3}{\frac {1}{a_{2}}}\int _{0}^{a_{2}}dx_{2}{\frac {1}{a_{1}}}\int _{0}^{a_{1}}dx_{1}g(x_{1},x_{2},x_{3})\cdot e^{-i2\pi \left({\frac {m_{1}}{a_{1}}}x_{1}+{\frac {m_{2}}{a_{2}}}x_{2}+{\frac {m_{3}}{a_{3}}}x_{3}\right)}\end{aligned}}}

We write g {\displaystyle g}gas:

g ( x 1, x 2, x 3) = ∑ m 1 = − ∞ ∞ ∑ m 2 = − ∞ ∞ ∑ m 3 = − ∞ ∞ h t h r e e ( m 1, m 2, m 3) ⋅ e i 2 π m 1 a 1 x 1 ⋅ e i 2 π m 2 a 2 x 2 ⋅ e i 2 π m 3 a 3 x 3 {\displaystyle g(x_{1},x_{2},x_{3})=\sum _{m_{1}=-\infty }^{\infty }\sum _{m_{2}=-\infty }^{\infty }\sum _{m_{3} =-\infty }^{\infty }h^{\mathrm {three} }(m_{1},m_{2},m_{3})\cdot e^{i2\pi {\frac {m_{1}}{a_{1}}}x_{1}}\cdot e^{i2\pi {\frac {m_{2}}{a_{2}}}x_{2}}\cdot e^{i2\pi {\frac {m_{3}}{a_{3}}}x_{3}}}g(x_1, x_2, x_3)=\sum_{m_1=-\infty}^\infty \sum_{m_2=-\infty}^\infty \sum_{m_3=-\infty}^\infty h^\mathrm{three}(m_1, m_2, m_3) \cdot e^{i 2\pi \frac{m_1}{a_1} x_1} \cdot e^{i 2\pi \frac{m_2}{a_2} x_2}\cdot e^{i 2\pi \frac{m_3}{a_3} x_3}

Re-arranging:

g ( x 1, x 2, x 3) = ∑ m 1, m 2, m 3 ∈ Z h t h r e e ( m 1, m 2, m 3) ⋅ e i 2 π ( m 1 a 1 x 1 + m 2 a 2 x 2 + m 3 a 3 x 3). {\displaystyle g(x_{1},x_{2},x_{3})=\sum _{m_{1},m_{2},m_{3}\in \mathbb {Z} }h^{\mathrm {three} }(m_{1},m_{2},m_{3})\cdot e^{i2\pi \left({\frac {m_{1}}{a_{1}}}x_{1}+{\frac {m_{2}}{a_{2}}}x_{2}+{\frac {m_{3}}{a_{3}}}x_{3}\right)}.}g(x_1, x_2, x_3)=\sum_{m_1, m_2, m_3 \in \Z } h^\mathrm{three}(m_1, m_2, m_3) \cdot e^{i 2\pi \left( \frac{m_1}{a_1} x_1+ \frac{m_2}{a_2} x_2 + \frac{m_3}{a_3} x_3 \right)}.

Now, every reciprocal lattice vector can be written as G = ℓ 1 g 1 + ℓ 2 g 2 + ℓ 3 g 3 {\displaystyle \mathbf {G} =\ell _{1}\mathbf {g} _{1}+\ell _{2}\mathbf {g} _{2}+\ell _{3}\mathbf {g} _{3}}{\displaystyle \mathbf {G} =\ell _{1}\mathbf {g} _{1}+\ell _{2}\mathbf {g} _{2}+\ell _{3}\mathbf {g} _{3}}, where l i {\displaystyle l_{i}}l_{i}are integers and g i {\displaystyle \mathbf {g} _{i}}\mathbf{g}_iare the reciprocal lattice vectors, we can use the fact that g i ⋅ a j = 2 π δ i j {\displaystyle \mathbf {g_{i}} \cdot \mathbf {a_{j}} =2\pi \delta _{ij}}\mathbf{g_i} \cdot \mathbf{a_j}=2\pi\delta_{ij}to calculate that for any arbitrary reciprocal lattice vector G {\displaystyle \mathbf {G} }\mathbf {G} and arbitrary vector in space r {\displaystyle \mathbf {r} }\mathbf {r} , their scalar product is:

G ⋅ r = ( ℓ 1 g 1 + ℓ 2 g 2 + ℓ 3 g 3) ⋅ ( x 1 a 1 a 1 + x 2 a 2 a 2 + x 3 a 3 a 3) = 2 π ( x 1 ℓ 1 a 1 + x 2 ℓ 2 a 2 + x 3 ℓ 3 a 3). {\displaystyle \mathbf {G} \cdot \mathbf {r} =\left(\ell _{1}\mathbf {g} _{1}+\ell _{2}\mathbf {g} _{2}+\ell _{3}\mathbf {g} _{3}\right)\cdot \left(x_{1}{\frac {\mathbf {a} _{1}}{a_{1}}}+x_{2}{\frac {\mathbf {a} _{2}}{a_{2}}}+x_{3}{\frac {\mathbf {a} _{3}}{a_{3}}}\right)=2\pi \left(x_{1}{\frac {\ell _{1}}{a_{1}}}+x_{2}{\frac {\ell _{2}}{a_{2}}}+x_{3}{\frac {\ell _{3}}{a_{3}}}\right).}{\displaystyle \mathbf {G} \cdot \mathbf {r} =\left(\ell _{1}\mathbf {g} _{1}+\ell _{2}\mathbf {g} _{2}+\ell _{3}\mathbf {g} _{3}\right)\cdot \left(x_{1}{\frac {\mathbf {a} _{1}}{a_{1}}}+x_{2}{\frac {\mathbf {a} _{2}}{a_{2}}}+x_{3}{\frac {\mathbf {a} _{3}}{a_{3}}}\right)=2\pi \left(x_{1}{\frac {\ell _{1}}{a_{1}}}+x_{2}{\frac {\ell _{2}}{a_{2}}}+x_{3}{\frac {\ell _{3}}{a_{3}}}\right).}

And so it is clear that in our expansion, the sum is actually over reciprocal lattice vectors:

f ( r) = ∑ G h ( G) ⋅ e i G ⋅ r, {\displaystyle f (\ mathbf {r}) = \ sum _ {\ mathbf {G}} h (\ mathbf {G}) \ cdot e ^ {i \ mathbf {G} \ cdot \ mathbf {r}},}{\displaystyle f(\mathbf {r})=\sum _{\mathbf {G} }h(\mathbf {G})\cdot e^{i\mathbf {G} \cdot \mathbf {r} },}

где

h (G) = 1 a 3 ∫ 0 a 3 dx 3 1 a 2 ∫ 0 a 2 dx 2 1 a 1 ∫ 0 a 1 dx 1 f (x 1 a 1 a 1 + x 2 a 2 a 2 + x 3 a 3 a 3) ⋅ e - i G ⋅ r. {\ displaystyle h (\ mathbf {G}) = {\ frac {1} {a_ {3}}} \ int _ {0} ^ {a_ {3}} dx_ {3} \, {\ frac {1} {a_ {2}}} \ int _ {0} ^ {a_ {2}} dx_ {2} \, {\ frac {1} {a_ {1}}} \ int _ {0} ^ {a_ {1 }} dx_ {1} \, f \ left (x_ {1} {\ frac {\ mathbf {a} _ {1}} {a_ {1}}} + x_ {2} {\ frac {\ mathbf {a } _ {2}} {a_ {2}}} + x_ {3} {\ frac {\ mathbf {a} _ {3}} {a_ {3}}} \ right) \ cdot e ^ {- i \ mathbf {G} \ cdot \ mathbf {r}}.}{\displaystyle h(\mathbf {G})={\frac {1}{a_{3}}}\int _{0}^{a_{3}}dx_{3}\,{\frac {1}{a_{2}}}\int _{0}^{a_{2}}dx_{2}\,{\frac {1}{a_{1}}}\int _{0}^{a_{1}}dx_{1}\,f\left(x_{1}{\frac {\mathbf {a} _{1}}{a_{1}}}+x_{2}{\frac {\mathbf {a} _{2}}{a_{2}}}+x_{3}{\frac {\mathbf {a} _{3}}{a_{3}}}\right)\cdot e^{-i\mathbf {G} \cdot \mathbf {r} }.}

Предполагая, что

r = (x, y, z) = x 1 a 1 a 1 + x 2 a 2 a 2 + x 3 a 3 a 3, {\ displaystyle \ mathbf {r} = (x, y, z) = x_ {1} {\ frac {\ mathbf {a} _ {1}} {a_ {1}}} + x_ {2} { \ frac {\ mathbf {a} _ {2}} {a_ {2}}} + x_ {3} {\ frac {\ mathbf {a} _ {3}} {a_ {3}}},}{\displaystyle \mathbf {r} =(x,y,z)=x_{1}{\frac {\mathbf {a} _{1}}{a_{1}}}+x_{2}{\frac {\mathbf {a} _{2}}{a_{2}}}+x_{3}{\frac {\mathbf {a} _{3}}{a_{3}}},}

мы можем решить эту систему трех линейных уравнений относительно x {\ displaystyle x}x, y {\ displaystyle y}yи z {\ displaystyle z}zв виде x 1 {\ displaystyle x_ {1}}x_{1}, x 2 {\ displaystyle x_ {2}}x_{2}и x 3 {\ displaystyle x_ {3} }x_{3}для вычисления элемента объема в исходной декартовой системе координат. Если у нас есть x {\ displaystyle x}x, y {\ displaystyle y}yи z {\ displaystyle z}zв терминах x 1 {\ displaystyle x_ {1}}x_{1}, x 2 {\ displaystyle x_ {2}}x_{2}и x 3 {\ displaystyle x_ {3}}x_{3}, мы можем вычислить определитель Якоби :

| ∂ x 1 ∂ x ∂ x 1 ∂ y ∂ x 1 ∂ z ∂ x 2 ∂ x ∂ x 2 ∂ y ∂ x 2 ∂ z ∂ x 3 ∂ x ∂ x 3 ∂ y ∂ x 3 ∂ z | {\ displaystyle {\ begin {vmatrix} {\ dfrac {\ partial x_ {1}} {\ partial x}} {\ dfrac {\ partial x_ {1}} {\ partial y}} {\ dfrac {\ частичный x_ {1}} {\ partial z}} \\ [12pt] {\ dfrac {\ partial x_ {2}} {\ partial x}} {\ dfrac {\ partial x_ {2}} {\ partial y }} {\ dfrac {\ partial x_ {2}} {\ partial z}} \\ [12pt] {\ dfrac {\ partial x_ {3}} {\ partial x}} {\ dfrac {\ partial x_ {3}} {\ partial y}} {\ dfrac {\ partial x_ {3}} {\ partial z}} \ end {vmatrix}}}{\displaystyle {\begin{vmatrix}{\dfrac {\partial x_{1}}{\partial x}}{\dfrac {\partial x_{1}}{\partial y}}{\dfrac {\partial x_{1}}{\partial z}}\\[12pt]{\dfrac {\partial x_{2}}{\partial x}}{\dfrac {\partial x_{2}}{\partial y}}{\dfrac {\partial x_{2}}{\partial z}}\\[12pt]{\dfrac {\partial x_{3}}{\partial x}}{\dfrac {\partial x_{3}}{\partial y}}{\dfrac {\partial x_{3}}{\partial z}}\end{vmatrix}}}

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

a 1 a 2 a 3 a 1 ⋅ (a 2 × a 3) {\ displaystyle {\ frac {a_ {1} a_ {2} a_ {3}} { \ mathbf {a_ {1}} \ cdot (\ mathbf {a_ {2}} \ times \ mathbf {a_ {3}})}}}\frac{a_1 a_2 a_3}{\mathbf{a_1}\cdot(\mathbf{a_2} \times \mathbf{a_3})}

(может быть полезно для упрощения вычислений, чтобы работать в такой декартовой системе координат, в которой так уж получилось, что a 1 {\ displaystyle \ mathbf {a_ {1}}}{\displaystyle \mathbf {a_{1}} }параллельно оси x, a 2 { \ displaystyle \ mathbf {a_ {2}}}{\ displaystyle \ mathbf {a_ {2}}} лежит в плоскости xy, а a 3 {\ displaystyle \ m athbf {a_ {3}}}{\displaystyle \mathbf {a_{3}} }имеет компоненты всех трех осей). Знаменатель - это в точности объем примитивной элементарной ячейки, которая окружена тремя примитивными векторами a 1 {\ displaystyle \ mathbf {a_ {1}}}{\displaystyle \mathbf {a_{1}} }, a 2 {\ displaystyle \ mathbf {a_ { 2}}}{\ displaystyle \ mathbf {a_ {2}}} и a 3 {\ displaystyle \ mathbf {a_ {3}}}{\displaystyle \mathbf {a_{3}} }. В частности, теперь мы знаем, что

d x 1 d x 2 d x 3 = a 1 a 2 a 3 a 1 ⋅ (a 2 × a 3) ⋅ d x d y d z. {\ displaystyle dx_ {1} \, dx_ {2} \, dx_ {3} = {\ frac {a_ {1} a_ {2} a_ {3}} {\ mathbf {a_ {1}} \ cdot (\ mathbf {a_ {2}} \ times \ mathbf {a_ {3}})}} \ cdot dx \, dy \, dz.}dx_1 \, dx_2 \, dx_3 = \frac{a_1 a_2 a_3}{\mathbf{a_1}\cdot(\mathbf{a_2} \times \mathbf{a_3})} \cdot dx \, dy \, dz.

Теперь мы можем написать h (K) {\ displaystyle h ( \ mathbf {K})}{\displaystyle h(\mathbf {K})}как интеграл с традиционной системой координат по объему примитивной ячейки, вместо x 1 {\ displaystyle x_ {1}}x_{1}, x 2 {\ displaystyle x_ {2}}x_{2}и x 3 {\ displaystyle x_ {3}}x_{3}переменные:

h (G) = 1 a 1 ⋅ ( a 2 × a 3) ∫ C drf (r) ⋅ е - я К ⋅ r {\ displaystyle h (\ mathbf {G}) = {\ frac {1} {\ mathbf {a_ {1}} \ cdot (\ mathbf {a_ {2}} \ times \ mathbf {a_ {3}})}} \ int _ {C} d \ mathbf {r} f (\ mathbf {r}) \ cdot e ^ {- i \ mathbf { K} \ cdot \ mathbf {r}}}{\displaystyle h(\mathbf {G})={\frac {1}{\mathbf {a_{1}} \cdot (\mathbf {a_{2}} \times \mathbf {a_{3}})}}\int _{C}d\mathbf {r} f(\mathbf {r})\cdot e^{-i\mathbf {K} \cdot \mathbf {r} }}

И C {\ displaystyle C}C- примитивная элементарная ячейка, таким образом, a 1 ⋅ (a 2 × a 3) {\ displaystyle \ mathbf {a_ {1}} \ cdot (\ mathbf {a_ {2}} \ times \ mathbf {a_ {3}})}\mathbf{a_1}\cdot(\mathbf{a_2} \times \mathbf{a_3})- объем примитивной элементарной ячейки.

Интерпретация гильбертова пространства

На языке гильбертовых пространств набор функций {en = einx: n ∈ Z} {\ displaystyle \ {e_ {n} = e ^ {inx}: n \ in \ mathbb {Z} \}}{\displaystyle \{e_{n}=e^{inx}:n\in \mathbb {Z} \}}- это ортонормированный базис для пространства L 2 ([- π, π]) {\ displaystyle L ^ {2} ([- \ pi, \ pi])}{\displaystyle L^{2}([-\pi,\pi ])}квадратично интегрируемых функций на [- π, π] {\ displaystyle [- \ pi, \ pi]}[-\pi,\pi ]. Это пространство на самом деле является гильбертовым пространством с внутренним продуктом, заданным для любых двух элементов f {\ displaystyle f}fи g {\ displaystyle g}gпо

⟨f, g⟩ ≜ 1 2 π ∫ - π π f (x) g (x) ¯ dx. {\ Displaystyle \ langle f, \, g \ rangle \; \ треугольник \; {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} f (x) {\ overline {g (x)}} \, dx.}{\displaystyle \langle f,\,g\rangle \;\triangleq \;{\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x){\overline {g(x)}}\,dx.}

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

f = ∑ n = - ∞ ∞ ⟨f, en⟩ en. {\ displaystyle f = \ sum _ {n = - \ infty} ^ {\ infty} \ langle f, e_ {n} \ rangle \, e_ {n}.}f=\sum_{n=-\infty}^\infty \langle f,e_n \rangle \, e_n.
Синусы и косинусы образуют ортонормированный набор, как показано выше. Интеграл синуса, косинуса и их произведения равен нулю (зеленая и красная области равны и сокращаются), когда m {\ displaystyle m}m, n {\ displaystyle n}nили функции различны, а число "пи" - только в том случае, если m {\ displaystyle m}mи n {\ displaystyle n}nравны, а используемая функция одинакова.

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

∫ - π π cos ⁡ (mx) cos ⁡ (nx) dx = 1 2 ∫ - π π cos ⁡ ((n - m) x) + cos ⁡ ((N + м) Икс) dx знак равно π δ мn, м, N ≥ 1, {\ Displaystyle \ int _ {- \ pi} ^ {\ pi} \ cos (mx) \, \ cos (nx) \, dx = {\ frac {1} {2}} \ int _ {- \ pi} ^ {\ pi} \ cos ((nm) x) + \ cos ((n + m) x) \, dx = \ pi \ delta _ {mn}, \ quad m, n \ geq 1, \,}{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\cos(nx)\,dx={\frac {1}{2}}\int _{-\pi }^{\pi }\cos((n-m)x)+\cos((n+m)x)\,dx=\pi \delta _{mn},\quad m,n\geq 1,\,}
∫ - π π sin ⁡ (mx) sin ⁡ (nx) dx = 1 2 ∫ - π π cos ⁡ ((n - м) Икс) - соз ⁡ ((N + м) Икс) dx знак равно π δ мn, м, N ≥ 1 {\ Displaystyle \ int _ {- \ pi} ^ {\ pi} \ sin (mx) \, \ sin (nx) \, dx = {\ frac {1} {2}} \ int _ {- \ pi} ^ {\ pi} \ cos ((nm) x) - \ cos ((n + m) x) \, dx = \ pi \ delta _ {mn}, \ quad m, n \ geq 1}{\displaystyle \int _{-\pi }^{\pi }\sin(mx)\,\sin(nx)\,dx={\frac {1}{2}}\int _{-\pi }^{\pi }\cos((n-m)x)-\cos((n+m)x)\,dx=\pi \delta _{mn},\quad m,n\geq 1}

(где δ mn - дельта Кронекера ) и

∫ - π π cos ⁡ (mx) sin ⁡ (nx) dx = 1 2 ∫ - π π sin ⁡ ((n + m) x) + sin ⁡ ((n - m) x) dx = 0; {\ Displaystyle \ int _ {- \ pi} ^ {\ pi} \ cos (mx) \, \ sin (nx) \, dx = {\ frac {1} {2}} \ int _ {- \ pi} ^ {\ pi} \ sin ((n + m) x) + \ sin ((nm) x) \, dx = 0; \,}{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\sin(nx)\,dx={\frac {1}{2}}\int _{-\pi }^{\pi }\sin((n+m)x)+\sin((n-m)x)\,dx=0;\,}

кроме того, синусы и косинусы ортогональны постоянной функции 1 {\ displaystyle 1}1. Ортонормированный базис для L 2 ([- π, π]) {\ displaystyle L ^ {2} ([- \ pi, \ pi])}{\displaystyle L^{2}([-\pi,\pi ])}, состоящий из действительных функций, образован функции 1 {\ displaystyle 1}1и 2 cos ⁡ (nx) {\ displaystyle {\ sqrt {2}} \ cos (nx)}{\displaystyle {\sqrt {2}}\cos(nx)}, 2 sin ⁡ (nx) {\ displaystyle {\ sqrt {2}} \ sin (nx)}{\displaystyle {\sqrt {2}}\sin(nx)}с n = 1, 2,... Плотность их диапазона является следствием теоремы Стоуна – Вейерштрасса, но также следует из свойств классических ядер, таких как ядро ​​Фейера.

Свойства

Таблица основных свойств

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

  • z ∗ {\ displaystyle z ^ {*}}z^{*}- комплексное сопряжение числа z {\ displaystyle z}z.
  • f (x), g (x) {\ displaystyle f (x), g (x)}{\displaystyle f(x),g(x)}обозначает P {\ displaystyle P}P-периодические функции, определенные на 0 < x ≤ P {\displaystyle 0{\displaystyle 0<x\leq P}.
  • F [n ], G [n] {\ displaystyle F [n], G [n]}{\displaystyle F[n],G[n]}обозначают коэффициенты ряда Фурье (экспоненциальная форма) для f {\ displaystyle f}fи g {\ displaystyle g}g, как определено в уравнении уравнение 5 .
СвойствоВременная областьЧастотная область ( экспоненциальная форма)ПримечанияСсылка
Линейностьa ⋅ f (x) + b ⋅ g (x) {\ displaystyle a \ cdot f (x) + b \ cdot g (x)}{\displaystyle a\cdot f(x)+b\cdot g(x)}a ⋅ F [n] + b ⋅ G [n] {\ displaystyle a \ cdot F [n] + b \ cdot G [n]}{\ displaystyle a \ cdot F [n] + b \ cdot G [n]} комплексные числа a, b {\ displaystyle a, b}a,b
Поворот времени / Изменение частотыf (- x) {\ displaystyle f (-x)}f(-x)F [- n] {\ displaystyle F [-n] }{\displaystyle F[-n]}
Временное сопряжениеf (x) ∗ {\ displaystyle f (x) ^ {*}}{\displaystyle f(x)^{*}}F [- n] ∗ {\ displaystyle F [-n] ^ {*}}{\ displaystyle F [-n] ^ {*}}
Обращение времени и спряжениеf (- x) ∗ {\ displaystyle f (-x) ^ {*}}{\displaystyle f(-x)^{*}}F [n] ∗ {\ displaystyle F [n] ^ {*}}{\displaystyle F[n]^{*}}
Реальная часть времениRe ⁡ (f (x)) {\ displaystyle \ operatorname {Re} {(f (x)) }}{\displaystyle \operatorname {Re} {(f(x))}}1 2 (F [n] + F [- n] ∗) {\ displaystyle {\ frac {1} {2}} (F [n] + F [-n] ^ {*})}{\displaystyle {\frac {1}{2}}(F[n]+F[-n]^{*})}
Мнимая часть времениIm ⁡ (f (x)) {\ displaystyle \ operatorname {Im} {(f (x))}}{\displaystyle \operatorname {Im} {(f(x))}}1 2 i (F [n] - F [- n] ∗) {\ displaystyle {\ frac {1} {2i}} (F [n] -F [-n] ^ {*})}{\displaystyle {\frac {1}{2i}}(F[n]-F[-n]^{*})}
Действительная часть частоты1 2 (f (x) + е (- х) *) {\ Displaystyle {\ гидроразрыва {1} {2}} (е (х) + f (-x) ^ {*})}{\ displaystyle {\ frac {1} {2}} (е (х) + f (-x) ^ {*})} Re ⁡ (F [n]) {\ displaystyle \ operatorname {Re} {(F [n])}}{\displaystyle \operatorname {Re} {(F[n])}}
Мнимая часть по частоте1 2 i (f (x) - f (- x) ∗) {\ displaystyle {\ frac { 1} {2i}} (f (x) -f (-x) ^ {*})}{\displaystyle {\frac {1}{2i}}(f(x)-f(-x)^{*})}Im ⁡ (F [n]) {\ displaystyle \ operatorname {Im} {(F [n])} }{\displaystyle \operatorname {Im} {(F[n])}}
Сдвиг во времени / Модуляция по частотеf (x - x 0) {\ displaystyle f (x-x_ {0})}{\displaystyle f(x-x_{0})}F [n] ⋅ e - i 2 π x 0 T п {\ disp Laystyle F [n] \ cdot e ^ {- i {\ frac {2 \ pi x_ {0}} {T}} n}}{\displaystyle F[n]\cdot e^{-i{\frac {2\pi x_{0}}{T}}n}}вещественное число x 0 {\ displaystyle x_ {0}}x_{0}
Сдвиг частоты / Модуляция во времениf (x) ⋅ ei 2 π n 0 T x {\ displaystyle f (x) \ cdot e ^ {i {\ frac {2 \ pi n_ {0}} { T}} x}}{\displaystyle f(x)\cdot e^{i{\frac {2\pi n_{0}}{T}}x}}F [n - n 0] {\ displaystyle F [n-n_ {0}] \!}{\displaystyle F[n-n_{0}]\!}целое число n 0 {\ displaystyle n_ {0}}n_{0}

Свойства симметрии

Когда действительная и мнимая части сложной функции разлагаются на их четную и нечетную части, есть четыре компонента, обозначенных ниже индексами RE, RO, IE, и IO. И существует взаимно-однозначное соответствие между четырьмя компонентами комплексной функции времени и четырьмя компонентами ее комплексного частотного преобразования:

Временная область f = f RE + f RO + if IE + if IO ⏟ ⇕ F ⇕ F ⇕ F ⇕ F ⇕ F Частотная область F = FRE + i FIO ⏞ + i FIE + FRO {\ displaystyle {\ begin {array} {rccccccccc} {\ text {Time domain}} f = f _ {_ {\ text {RE}}} + f _ {_ {\ text {RO}}} + if _ {_ {\ text {IE}}} + \ underbrace {i \ f _ {_ {\ text {IO}}} } \\ {\ Bigg \ Updownarrow} {\ mathcal {F}} {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \\ {\ text {Частотный домен}} F = F_ {RE} + \ overbrace {i \ F_ {IO}} + i \ F_ {IE} + F_ {RO} \ end {array}}}{\ displaystyle {\ begin {array} {rccccccccc} {\ text {Time domain}} f = f _ {_ {\ text {RE}}} + f _ {_ {\ text {RO}}} + if _ {_ {\ text {IE}}} + \ underbrace {i \ f _ {_ {\ text {IO}} }} \\ {\ Bigg \ Updownarrow} {\ mathcal {F}} {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \ \ {\ Bigg \ Updownarrow} {\ mathcal {F}} \\ {\ text {Frequency domain}} F = F_ {RE} + \ overbrace {я \ F_ {IO}} + я \ F_ {IE} + F_ {RO} \ end {array}}}

Отсюда очевидны различные отношения, например:

  • Преобразование функция с действительным знаком (f RE+ f RO) - это даже симметричная функция F RE+ i F IO. И наоборот, четно-симметричное преобразование подразумевает действительную временную область.
  • Преобразование мнимозначной функции (если IE+, если IO) является нечетно-симметричным функция F RO+ i F IE, и верно обратное.
  • Преобразование четно-симметричной функции (f RE+ if IO) является вещественнозначной функцией F RE+ F RO, и верно обратное.
  • Преобразование нечетно-симметричной функции (f RO+ if IE) является мнимозначной функцией i F IE+ i F IO, и верно обратное.

Лемма Римана – Лебега

Если f {\ displaystyle f}fинтегрируемо, лим | п | → ∞ е ^ (n) знак равно 0 {\ displaystyle \ lim _ {| n | \ rightarrow \ infty} {\ hat {f}} (n) = 0}\ lim_ {| n | \ rightarrow \ infty} \ hat {f} (n) = 0 , lim n → + ∞ an = 0 {\ displaystyle \ lim _ {n \ rightarrow + \ infty} a_ {n} = 0}\lim_{n\rightarrow +\infty}a_n=0и lim n → + ∞ bn = 0. {\ displaystyle \ lim _ {n \ rightarrow + \ infty} b_ {n} = 0.}\lim_{n\rightarrow +\infty}b_n=0.Этот результат известен как лемма Римана – Лебега.

Свойство производной

Мы говорим, что f {\ displaystyle f}fпринадлежит C k (T) {\ displaystyle C ^ {k} (\ mathbb {T})}C^k(\mathbb{T})если f {\ displaystyle f}f- 2π-периодическая функция на R {\ displaystyle \ mathbb {R}}\mathbb {R} , которая k {\ displaystyle k}k раз дифференцируема, а ее k-я производная непрерывна.

  • Если f ∈ C 1 (T) {\ displaystyle f \ in C ^ {1} (\ mathbb {T})}f \in C^1(\mathbb{T}), то коэффициенты Фурье f ′ ^ (n) {\ displaystyle {\ widehat {f '}} (n)}\widehat{f'}(n)производной f ′ {\ displaystyle f'}f'может быть выражено через коэффициенты Фурье f ^ (n) {\ displaystyle {\ widehat {f}} (n)}\widehat{f}(n)функции f {\ displaystyle f}f, по формуле f '^ (n) = inf ^ (n) {\ displaystyle {\ widehat {f'}} (n) = in {\ widehat {f}} (n)}{\displaystyle {\widehat {f'}}(n)=in{\widehat {f}}(n)}.
  • Если f ∈ C k (T) {\ displaystyle f \ in C ^ {k} (\ mathbb {T})}f \in C^k(\mathbb{T}), затем f (k) ^ (n) = (in) kf ^ (n) {\ displaystyle {\ widehat {f ^ {(k)}}} (n) = (дюйм) ^ {k} {\ widehat {f}} (n)}{\displaystyle {\widehat {f^{(k)}}}(n)=(in)^{k}{\widehat {f}}(n)}. В частности, поскольку для фиксированного k ≥ 1 {\ displaystyle k \ geq 1}k\geq 1мы имеем f (k) ^ (n) → 0 {\ displaystyle {\ widehat {f ^ {(k)}}} (n) \ to 0}{\displaystyle {\widehat {f^{(k)}}}(n)\to 0}как n → ∞ {\ displaystyle n \ to \ infty}n\to \infty , следует, что | п | kf ^ (n) {\ displaystyle | n | ^ {k} {\ widehat {f}} (n)}{\ displaystyle | n | ^ {k} {\ widehat {f}} (n)} стремится к нулю, что означает, что коэффициенты Фурье сходятся к нулю быстрее, чем k-я степень из n для любого k ≥ 1 {\ displaystyle k \ geq 1}k\geq 1.

теорема Парсеваля

Если f {\ displaystyle f}fпринадлежит L 2 ([- π, π]) {\ displaystyle L ^ {2} ([- \ pi, \ pi])}{\displaystyle L^{2}([-\pi,\pi ])}, затем ∑ n = - ∞ ∞ | f ^ (n) | 2 = 1 2 π ∫ - π π | f (x) | 2 dx {\ displaystyle \ sum _ {n = - \ infty} ^ {\ infty} | {\ hat {f}} (n) | ^ {2} = {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} | f (x) | ^ {2} \, dx}\ sum_ {n = - \ infty} ^ \ infty | \ hat {f} (n) | ^ 2 = \ frac {1} {2 \ pi} \ int _ {- \ pi} ^ {\ pi} | f (x) | ^ 2 \, dx .

Теорема Планшереля

Если c 0, c ± 1, c ± 2,… {\ displaystyle c_ {0}, \, c _ {\ pm 1}, \, c _ {\ pm 2}, \ ldots}c_0,\, c_{\pm 1},\, c_{\pm 2},\ldots- коэффициенты, а ∑ n = - ∞ ∞ | c n | 2 < ∞ {\displaystyle \sum _{n=-\infty }^{\infty }|c_{n}|^{2}<\infty }\sum_{n=-\infty}^\infty |c_n|^2 <\inftyтогда существует уникальная функция f ∈ L 2 ([- π, π]) {\ displaystyle f \ in L ^ {2} ([- \ pi, \ pi])}f\in L^2([-\pi,\pi])так, что f ^ (n) = cn {\ displaystyle {\ hat {f}} (n) = c_ {n}}\hat{f}(n) = c_nдля каждого n {\ displaystyle n }n.

Теоремы о свертке

  • Первая теорема о свертке утверждает, что если f {\ displaystyle f}fи g {\ displaystyle g}gнаходятся в L 1 ([- π, π]) {\ displaystyle L ^ {1} ([- \ pi, \ pi])}{\displaystyle L^{1}([-\pi,\pi ])}, коэффициенты ряда Фурье 2π-периодичности свертка из f {\ displaystyle f}fи g {\ displaystyle g}gзадаются следующим образом:
[f ∗ 2 π g ^ ] (N) знак равно 2 π ⋅ е ^ (N) ⋅ g ^ (n), {\ displaystyle [{\ widehat {f * _ {2 \ pi} g}}] (n) = 2 \ pi \ cdot { \ hat {f}} (n) \ cdot {\ hat {g}} (n),}[\widehat{f*_{2\pi}g}](n) = 2\pi\cdot \hat{f}(n)\cdot\hat{g}(n),
где:
[f ∗ 2 π g] (x) ≜ ∫ - π π f (u) ⋅ g [pv ⁡ (x - u)] du, (и pv ⁡ (x) ≜ Arg ⁡ (eix) ⏟ главное значение) = ∫ - π π f (u) ⋅ g (x - u) du, когда g (x) 2 π -периодична. = ∫ 2 π f (u) ⋅ g (x - u) d u, когда обе функции 2 π -периодичны, а интеграл ведется по любому интервалу 2 π. {\ Displaystyle {\ begin {align} \ left [е * _ {2 \ pi} g \ right] (x) \ \ треугольникq \ int _ {- \ pi} ^ {\ pi} f (u) \ cdot g [\ operatorname {pv} (xu)] \, du, {\ big (} {\ text {and}} \ underbrace {\ operatorname {pv} (x) \ \ triangleq \ operatorname {Arg} (e ^ {ix})} _ {\ text {главное значение}} \, {\ big)} \\ = \ int _ {- \ pi} ^ {\ pi} f (u) \ cdot g (xu) \, du, {\ text {when}} g (x) {\ text {is}} 2 \ pi {\ text {-periodic.}} \\ = \ int _ {2 \ pi} f (u) \ cdot g (xu) \, du, {\ text {когда обе функции являются}} 2 \ pi {\ text {-периодическими, а интеграл берется по любому}} 2 \ pi {\ text {interval.}} \ end {align}}}{\displaystyle {\begin{aligned}\left[f*_{2\pi }g\right](x)\ \triangleq \int _{-\pi }^{\pi }f(u)\cdot g[\operatorname {pv} (x-u)]\,du,{\big (}{\text{and }}\underbrace {\operatorname {pv} (x)\ \triangleq \operatorname {Arg} (e^{ix})} _{\text{principal value}}\,{\big)}\\=\int _{-\pi }^{\pi }f(u)\cdot g(x-u)\,du,{\text{when }}g(x){\text{ is }}2\pi {\text{-periodic.}}\\=\int _{2\pi }f(u)\cdot g(x-u)\,du,{\text{when both functions are }}2\pi {\text{-periodic, and the integral is over any }}2\pi {\text{ interval.}}\end{aligned}}}
  • Вторая теорема о свертке утверждает, что коэффициенты ряда Фурье произведения f {\ displaystyle f}fи g {\ displaystyle g}gзадаются дискретной сверткой элементов f ^ {\ displaystyle {\ hat {f}}}{\hat {f}}и g ^ {\ displaystyle { \ hat {g}}}\hat gпоследовательности:
[f ⋅ g ^] (n) = [f ^ ∗ g ^] (n). {\ displaystyle [{\ widehat {f \ cdot g}}] (n) = [{\ hat {f}} * {\ hat {g}}] (n).}[\widehat{f\cdot g}](n) = [\hat{f}*\hat{g}](n).
  • A дважды бесконечная последовательность {cn} n ∈ Z {\ displaystyle \ left \ {c_ {n} \ right \} _ {n \ in Z}}{\displaystyle \left\{c_{n}\right\}_{n\in Z}}in c 0 (Z) {\ displaystyle c_ {0} (\ mathbb {Z})}{\displaystyle c_{0}(\mathbb {Z})}- последовательность коэффициентов Фурье функции из L 1 ([0, 2 π]) {\ displaystyle L ^ {1} ( [0,2 \ pi])}{\displaystyle L^{1}([0,2\pi ])}тогда и только тогда, когда это свертка двух последовательностей в ℓ 2 (Z) {\ displaystyle \ ell ^ {2} (\ mathbb {Z})}\ell ^{2}({\mathbb {Z}}). См.

Компактные группы

Одно из интересных свойств преобразования Фурье, о котором мы упоминали, заключается в том, что оно переносит свертки в точечные произведения. Если это свойство, которое мы стремимся сохранить, можно построить ряд Фурье на любой компактной группе. Типичные примеры включают те классические группы, которые являются компактными. Это обобщает преобразование Фурье на все пространства вида L (G), где G - компактная группа, таким образом, что преобразование Фурье переводит свертки в точечные произведения. Ряд Фурье существует и сходится аналогично случаю [−π, π].

Альтернативным расширением компактных групп является теорема Питера – Вейля, которая доказывает результаты о представлениях компактных групп, аналогичные представлениям о конечных группах.

Римановы многообразия

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

Если домен не является группой, то внутренне определенной свертки нет. Однако, если X {\ displaystyle X}X- компакт риманово многообразие, оно имеет оператор Лапласа – Бельтрами. Оператор Лапласа – Бельтрами - это дифференциальный оператор, который соответствует оператору Лапласа для риманова многообразия X {\ displaystyle X}X. Затем по аналогии можно рассмотреть уравнения теплопроводности на X {\ displaystyle X}X. Поскольку Фурье пришел к своей основе, пытаясь решить уравнение теплопроводности, естественным обобщением является использование в качестве основы собственных решений оператора Лапласа – Бельтрами. Это обобщает ряд Фурье на пространства типа L 2 (X) {\ displaystyle L ^ {2} (X)}{\displaystyle L^{2}(X)}, где X {\ displaystyle X}X- риманово многообразие. Ряд Фурье сходится способами, аналогичными случаю [- π, π] {\ displaystyle [- \ pi, \ pi]}[-\pi,\pi ]. Типичный пример - взять X {\ displaystyle X}Xкак сферу с обычной метрикой, и в этом случае базис Фурье состоит из сферических гармоник.

Локально компактные абелевы группы

Обобщение на компактные группы, рассмотренное выше, не распространяется на некомпактные, неабелевы группы. Однако есть прямое обобщение на локально компактные абелевы (LCA) группы.

Это обобщает преобразование Фурье на L 1 (G) {\ displaystyle L ^ {1} (G)}{\displaystyle L^{1}(G)}или L 2 (G) {\ displaystyle L ^ {2} (G)}{\displaystyle L^{2}(G)}, где G {\ displaystyle G}G- группа LCA. Если G {\ displaystyle G}Gкомпактно, также получается ряд Фурье, сходящийся аналогично [- π, π] {\ displaystyle [- \ pi, \ pi ]}[-\pi,\pi ]случай, но если G {\ displaystyle G}Gнекомпактный, вместо него получается интеграл Фурье. Это обобщение приводит к обычному преобразованию Фурье, когда базовая локально компактная абелева группа имеет вид R {\ displaystyle \ mathbb {R}}\mathbb {R} .

Таблица общих рядов Фурье

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

  • f (x) {\ displaystyle f (x)}f(x)обозначает периодическую функцию, определенную на 0 < x ≤ T {\displaystyle 0{\displaystyle 0<x\leq T}.
  • a 0, an, bn {\ displaystyle a_ {0}, a_ {n}, b_ {n}}{\displaystyle a_{0},a_{n},b_{n}}обозначают коэффициенты ряда Фурье (синус-косинусная форма) периодической функции f {\ displaystyle f}f, как определено в Уравнение 4 .
Временная область. f (x) {\ displaystyle f (x)}f(x)ГрафикЧастотная область (синус-косинусная форма). a 0 an для n ≥ 1 bn для n ≥ 1 {\ displaystyle {\ begin {align} a_ {0} \\ a_ {n} \ quad {\ text {for}} n \ geq 1 \\ b_ {n} \ quad {\ text {for}} n \ geq 1 \ end {align}}}{\displaystyle {\begin{aligned}a_{0}\\a_{n}\quad {\t ext{for }}n\geq 1\\b_{n}\quad {\text{for }}n\geq 1\end{aligned}}}ПримечанияСсылка
f (x) = A | sin ⁡ (2 π T x) | для 0 ≤ x < T {\displaystyle f(x)=A\left|\sin \left({\frac {2\pi }{T}}x\right)\right|\quad {\text{for }}0\leq x{\displaystyle f(x)=A\left|\sin \left({\frac {2\pi }{T}}x\right)\right|\quad {\text{for }}0\leq x<T}PlotRectifiedSineSignal.svg a 0 = 4 A π an = {- 4 A π 1 n 2-1 n четный 0 n нечетный bn = 0 {\ displaystyle {\ begin {align} a_ {0} = {\ frac {4A} {\ pi}} \\ a_ {n} = {\ begin {cases} {\ frac {-4A} {\ pi}} {\ frac {1} {n ^ {2} -1} } \ quad n {\ text {even}} \\ 0 \ quad n {\ text {odd}} \ end {cases}} \\ b_ {n} = 0 \\\ end {align}}}{\displaystyle {\begin{aligned}a_{0}={\frac {4A}{\pi }}\\a_{n}={\begin{cases}{\frac {-4A}{\pi }}{\frac {1}{n^{2}-1}}\quad n{\text{ even}}\\0\quad n{\text{ odd}}\end{cases}}\\b_{n}=0\\\end{aligned}}}Двухполупериодный выпрямленный синус
f (x) = {A sin ⁡ (2 π T x) для 0 ≤ x < T / 2 0 for T / 2 ≤ x < T {\displaystyle f(x)={\begin{cases}A\sin \left({\frac {2\pi }{T}}x\right)\quad {\text{for }}0\leq x{\displaystyle f(x)={\begin{cases}A\sin \left({\frac {2\pi }{T}}x\right)\quad {\text{for }}0\le q x<T/2\\0\quad {\text{for }}T/2\leq x<T\\\end{cases}}}PlotHalfRectifiedSineSignal.svg a 0 = 2 A π an = {- 2 A π 1 1 - n 2 n четный 0 n нечетный bn = {A 2 n = 1 0 n>1 {\ displaystyle {\ begin {align} a_ {0} = {\ frac {2A} {\ pi}} \\ a_ {n} = {\ begin {case} {\ frac {-2A} {\ pi}} {\ frac {1} {1-n ^ {2}}} \ quad n {\ text {even}} \\ 0 \ quad n {\ text {odd}} \ end {cases}} \\ b_ {n} = {\ begin {cases} {\ frac {A} {2}} \ quad n = 1 \\ 0 \ quad n>1 \ end {cases}} \\\ end {align}}}{\displaystyle {\begin{aligned}a_{0}={\frac {2A}{\pi }}\\a_{n}={\begin{cases}{\frac {-2A}{\pi }}{\frac {1}{1-n^{2}}}\quad n{\text{ even}}\\0\quad n{\text{ odd}}\end{cases}}\\b_{n}={\begin{cases}{\frac {A}{2}}\quad n=1\\0\quad n>1 \ end {cases}} \\\ end {align}}} Полуволновой выпрямленный синус
{x) A для 0 ≤ x < D ⋅ T 0 for D ⋅ T ≤ x < T {\displaystyle f(x)={\begin{cases}A\quad {\text{for }}0\leq x{\displaystyle f(x)={\begin{cases}A\quad {\text{for }}0\leq x<D\cdot T\\0\quad {\text{for }}D\cdot T\leq x<T\\\end{cases}}}PlotRectangleSignal.svg a 0 = 2 AD an = A n π sin ⁡ (2 π n D) bn = 2 A N π (грех ⁡ (π N D)) 2 {\ displaystyle {\ begin {align} a_ {0} = 2AD \\ a_ {n} = {\ frac {A} {n \ pi }} \ sin \ left (2 \ pi nD \ right) \\ b_ {n} = {\ frac {2A} {n \ pi}} \ left (\ sin \ left (\ pi nD \ right) \ right) ^ {2} \\\ конец {выровнено}}}{\displaystyle {\begin{aligned}a_{0}=2AD\\a_{n}={\frac {A}{n\pi }}\sin \left(2\pi nD\right)\\b_{n}={\frac {2A}{n\pi }}\left(\sin \left(\pi nD\right)\right)^{2}\\\end{aligned}}}0 ≤ D ≤ 1 {\ displaystyle 0 \ leq D \ leq 1}{\displaystyle 0\leq D\leq 1}
f (x) = A x T для 0 ≤ x < T {\displaystyle f(x)={\frac {Ax}{T}}\quad {\text{for }}0\leq x{\displaystyle f(x)={\frac {Ax}{T}}\quad {\text{for }}0\leq x<T}PlotSawtooth1Signal.svg a 0 = A an = 0 bn = - A n π {\ displaystyle {\ begin {align} a_ {0} = A \\ a_ {n} = 0 \\ b_ {n} = {\ frac {- A} {n \ pi}} \\\ конец {выровнено}}}{\displaystyle {\begin{aligned}a_{0}=A\\a_{n}=0\\b_{n}={\frac {-A}{n\pi }}\\\end{aligned}}}
f (x) = A - A x T для 0 ≤ x < T {\displaystyle f(x)=A-{\frac {Ax}{T}}\quad {\text{for }}0\leq x{\displaystyle f(x)=A-{\frac {Ax}{T}}\quad {\text{for }}0\leq x<T}PlotSawtootdiv class= a 0 = A an = 0 bn = A n π { \ displaystyle {\ begin {align} a_ {0} = A \\ a_ {n} = 0 \\ b_ {n} = {\ frac {A} {n \ pi}} \\\ end {align}} }{\displaystyle {\begin{aligned}a_{0}=A\\a_{n}=0\\b_{n}={\frac {A}{n\pi }}\\\end{aligned}}}
f (x) = 4 AT 2 (x - T 2) 2 для 0 ≤ x < T {\displaystyle f(x)={\frac {4A}{T^{2}}}\left(x-{\frac {T}{2}}\right)^{2}\quad {\text{for }}0\leq x{\displaystyle f(x)={\frac {4A}{T^{2}}}\left(x-{\frac {T}{2}}\right)^{2}\quad {\text{for }}0\leq x<T}PlotParabolaSignal.svg a 0 = 2 A 3 an = 4 A π 2 n 2 bn = 0 {\ displaystyle {\ begin { выровнено} a_ {0} = {\ frac {2A} {3}} \\ a_ {n} = {\ frac {4A} {\ pi ^ {2} n ^ {2}}} \\ b_ { n} = 0 \\\ end {align}}}{\displaystyle {\begin{aligned}a_{0}={\frac {2A}{3}}\\a_{n}={\frac {4A}{\pi ^{2}n^{2}}}\\b_{n}=0\\\end{aligned}}}
Аппроксимация и сходимость рядов Фурье

Важным вопросом как для теории, так и для приложений является сходимость. В частности, в приложениях часто бывает необходимо заменить бесконечный ряд ∑ - ∞ ∞ {\ displaystyle \ sum _ {- \ infty} ^ {\ infty}}\sum_{-\infty}^\inftyконечным,

f N (x) знак равно ∑ n = - NN f ^ (n) einx. {\ displaystyle f_ {N} (x) = \ sum _ {n = -N} ^ {N} {\ hat {f}} (n) e ^ {inx}.}f_N(x) = \sum_{n=-N}^N \hat{f}(n) e^{inx}.

Это называется частичной суммой. Мы хотели бы знать, в каком смысле f N (x) {\ displaystyle f_ {N} (x)}{\displaystyle f_{N}(x)}сходится к f (x) {\ displaystyle f (x)}f(x)as N → ∞ {\ displaystyle N \ rightarrow \ infty}N\rightarrow \infty .

Свойство наименьших квадратов

Мы говорим, что p {\ displaystyle p}p - это тригонометрический полином степени N {\ displaystyle N}N, когда он имеет форму

p (x) = ∑ n = - NN pneinx. {\ displaystyle p (x) = \ sum _ {n = -N} ^ {N} p_ {n} e ^ {inx}.}p(x)=\sum_{n=-N}^N p_n e^{inx}.

Обратите внимание, что f N {\ displaystyle f_ {N}}f_{N}- тригонометрический полином степени N {\ displaystyle N}N. Теорема Парсеваля означает, что

теорема. Тригонометрический полином f N {\ displaystyle f_ {N}}f_{N}- единственный лучший тригонометрический полином степени N {\ displaystyle N}N, аппроксимирующий f (x) {\ displaystyle f (x)}f(x)в том смысле, что для любого тригонометрического полинома p ≠ f N {\ displaystyle p \ neq f_ {N}}{\displaystyle p\neq f_{N}}степени N {\ displaystyle N }N, имеем

‖ f N - f ‖ 2 < ‖ p − f ‖ 2, {\displaystyle \|f_{N}-f\|_{2}<\|p-f\|_{2},}\|f_N - f\|_2 <\|p - f\|_2,

где норма гильбертова пространства определяется как:

‖ g ‖ 2 = 1 2 π ∫ - π π | g (x) | 2 д х. {\ displaystyle \ | g \ | _ {2} = {\ sqrt {{1 \ over 2 \ pi} \ int _ {- \ pi} ^ {\ pi} | g (x) | ^ {2} \, dx}}.}\ | g \ | _2 = \ sqrt {{1 \ over 2 \ pi} \ int _ {- \ pi} ^ \ pi | g (x) | ^ 2 \, dx}.

Сходимость

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

Теорема. Если f {\ displaystyle f}fпринадлежит L 2 ([- π, π]) {\ displaystyle L ^ {2} (\ left [- \ pi, \ pi \ right])}{\displaystyle L^{2}(\left[-\pi,\pi \right])}, тогда f ∞ {\ displaystyle f _ {\ infty}}{\displaystyle f_{\infty }}сходится к f {\ displaystyle f }fв L 2 ([- π, π]) {\ displaystyle L ^ {2} (\ left [- \ pi, \ pi \ right])}{\displaystyle L^{2}(\left[-\pi,\pi \right])}, то есть ‖ е N - f ‖ 2 {\ displaystyle \ | f_ {N} -f \ | _ {2}}\|f_N - f\|_2сходится к 0 при N → ∞ {\ displaystyle N \ rightarrow \ infty}N\rightarrow \infty .

Мы уже упоминали, что если f {\ displaystyle f}fнепрерывно дифференцируемо, то (i ⋅ n) f ^ (n) { \ displaystyle (i \ cdot n) {\ hat {f}} (n)}( я \ cdot n) \ hat {f} (n) - n-й коэффициент Фурье производной f ′ {\ displaystyle f '}f'. По сути, из неравенства Коши – Шварца следует, что f ∞ {\ displaystyle f _ {\ infty}}{\displaystyle f_{\infty }}абсолютно суммируем. Сумма этого ряда является непрерывной функцией, равной f {\ displaystyle f}f, поскольку ряд Фурье сходится в среднем к f {\ displaystyle f}f:

Теорема. Если f ∈ C 1 (T) {\ displaystyle f \ in C ^ {1} (\ mathbb {T})}f \in C^1(\mathbb{T}), то f ∞ {\ displaystyle f _ {\ infty}}{\displaystyle f_{\infty }}сходится к f {\ displaystyle f}fравномерно (и, следовательно, также точечно.)

Это результат может быть легко доказан, если f {\ displaystyle f}fдалее предполагается равным C 2 {\ displaystyle C ^ {2}}C^{2}, поскольку в этом case n 2 f ^ (n) {\ displaystyle n ^ {2} {\ hat {f}} (n)}n^2\hat{f}(n)стремится к нулю при n → ∞ {\ displaystyle n \ rightarrow \ infty}n \ rightarrow \ infty . В более общем смысле, ряд Фурье абсолютно суммируем, поэтому сходится равномерно к f {\ displaystyle f}f, при условии, что f {\ displaystyle f}fудовлетворяет условие Гельдера порядка α>1/2 {\ displaystyle \ alpha>1/2}{\displaystyle \alpha>1/2} . В абсолютно суммируемом случае неравенство sup - f | f (x (х) | ≤ ∑ | N |>N | f ^ (n) | {\ displaystyle \ sup _ {x} | f (x) -f_ {N} (x) | \ leq \ sum _ {| n |>N} | {\ hat {f}} (n) |}\sup_x |f(x) - f_N(x)| \le \sum_{|n|>N} | \ hat {f} (n) | доказывает равномерную конвергенцию.

Известно много других результатов, касающихся сходимости рядов Фурье, начиная от умеренно простого результата, что ряд сходится при x {\ displaystyle x}xесли f {\ displaystyle f}fдифференцируем в x {\ displaystyle x}x, в гораздо более сложный результат Леннарта Карлесона что ряд Фурье функции L 2 {\ displaystyle L ^ {2}}L ^ {2} на самом деле сходится почти везде.

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

Дивергенция

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

В 1922 году Андрей Колмогоров опубликовал статью «Une série de Fourier-Lebesgue divergente presque partout», в которой он привел пример интегрируемой по Лебегу функции, ряд Фурье которой расходится почти всюду. Позже он построил пример интегрируемой функции, ряд Фурье которой расходится всюду (Katznelson 1976).

См. Также
Примечания
Ссылки

Дополнительная литература

Внешние ссылки

. Эта статья включает материал из примера серии Фурье на PlanetMath, который находится под лицензией Creative Commons Attribution / Share-Alike License.

Последняя правка сделана 2021-05-20 12:53:34
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте