Ряды Фурье

Идея разложения в тригонометрический ряд

Жан Батист Жозеф Фурье в 1807 году предложил революционную идею: любую «разумную» функцию можно разложить в ряд по тригонометрическим функциям. Это позволяло решать уравнение теплопроводности — задачу, с которой Фурье работал как инженер.

Тригонометрический ряд Фурье: f(x) = a₀/2 + Σₙ₌₁^∞ (aₙ cos(nx) + bₙ sin(nx)), где

aₙ = (1/π)∫₋π^π f(x)cos(nx)dx, bₙ = (1/π)∫₋π^π f(x)sin(nx)dx.

Ортогональность тригонометрической системы

Функции 1, cos x, sin x, cos 2x, sin 2x, ... ортогональны в пространстве L²[-π, π]:

∫₋π^π cos(mx)cos(nx)dx = 0 при m ≠ n; = π при m = n > 0.

Именно ортогональность позволяет вычислить коэффициенты: умножаем f(x) на cos(nx) и интегрируем — все остальные слагаемые исчезают.

Теорема Дирихле

Если f кусочно-непрерывна и кусочно-монотонна на [-π, π], то ряд Фурье сходится к f(x) в точках непрерывности и к (f(x+)+f(x-))/2 в точках разрыва.

Равенство Парсеваля

a₀²/2 + Σ(aₙ² + bₙ²) = (1/π)∫₋π^π [f(x)]²dx.

Это «теорема Пифагора» в бесконечномерном пространстве: сумма квадратов коэффициентов = квадрат нормы функции. Ряд Фурье — наилучшее приближение в смысле среднеквадратичной нормы.

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

В комплексных обозначениях: f(x) = Σcₙ eⁱⁿˣ, cₙ = (1/2π)∫₋π^π f(x)e⁻ⁱⁿˣdx.

Для функций на всей прямой: преобразование Фурье f̂(ω) = ∫₋∞^∞ f(x)e⁻ⁱωˣdx.

Это главный инструмент в обработке сигналов, квантовой механике, решении дифференциальных уравнений. Быстрое преобразование Фурье (FFT) — алгоритм Кули–Тьюки — один из 10 важнейших алгоритмов XX века.

Сходимость ряда Фурье: тонкости

Теорема Дирихле устанавливает поточечную сходимость для «приличных» функций. Но общий ответ сложнее. Теорема Фейера (1904): средние арифметические частичных сумм (суммы Фейера) сходятся равномерно к f для любой непрерывной периодической функции. Это слабее поточечной сходимости, но зато работает для всех непрерывных функций.

Феномен Гиббса: в окрестности точки разрыва частичная сумма ряда Фурье «перерегулирует» — максимальное отклонение ≈ 9% от величины скачка, независимо от числа удержанных гармоник. Это принципиально ограничивает точность Фурье-аппроксимации у разрывов.

Равенство Парсеваля и энергия сигнала

Равенство Парсеваля a₀²/2 + Σ(aₙ² + bₙ²) = (1/π)∫₋π^π f²(x)dx имеет физический смысл: суммарная энергия сигнала = сумма энергий гармоник. В радиотехнике — это закон сохранения энергии в спектральном представлении. Для цифрового аудио: при сжатии MP3 удаляют гармоники с малыми коэффициентами aₙ², bₙ² — потеря «незначительной» энергии даёт значительное сжатие файла без заметного ухудшения качества звука.

Непериодические функции и преобразование Фурье

Для функций, определённых на всей прямой (не периодических), ряд Фурье обобщается до интегрального преобразования Фурье f̂(ω) = ∫₋∞^∞ f(x)e^(-iωx)dx. Это «предел» ряда Фурье при периоде T → ∞: дискретный спектр {cₙ} переходит в непрерывный f̂(ω). Формула обращения: f(x) = (1/2π)∫₋∞^∞ f̂(ω)e^(iωx)dω.

В механике квантового поля и обработке сигналов используется именно непрерывный спектр. Разность дискретного и непрерывного случаев — это разница между периодическими и непериодическими явлениями.

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

В двух измерениях преобразование Фурье f̂(u, v) = ∫∫ f(x, y) e^(−2πi(ux+vy)) dx dy разлагает изображение по плоским гармоникам. Низкие частоты (малые u, v) соответствуют крупным деталям изображения — общему фону и контурам; высокие частоты — мелким деталям, текстурам, краям объектов. Фильтрация в частотной области: умножение f̂(u,v) на маску H(u,v) и обратное преобразование — это быстрый способ размыть (фильтр низких частот) или заострить (фильтр высоких частот) изображение. Алгоритм JPEG использует дискретное косинусное преобразование (ДКП) — вещественный аналог Фурье — к блокам 8×8 пикселей: при сжатии обнуляют коэффициенты высоких частот. Удаление 80–90% коэффициентов при качестве JPEG=75 даёт визуально приемлемый результат при сжатии в 5–10 раз.

Принцип неопределённости в анализе сигналов

Преобразование Фурье реализует математический аналог принципа неопределённости Гейзенберга: сигнал не может быть одновременно сосредоточен и во времени, и в частотной области. Точная формулировка: σₜ · σ_ω ≥ 1/2, где σₜ и σ_ω — среднеквадратические отклонения сигнала и его Фурье-образа. Это объясняет, почему короткий радиоимпульс имеет широкую полосу частот — ключевой факт в теории связи и цифровой обработке сигналов (DSP).

Вопрос для размышления: Почему ряд Фурье «прямоугольного» сигнала (0 и 1 попеременно) имеет только нечётные гармоники? Как связана симметрия функции с видом её Фурье-разложения?

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

Прямое вычисление коэффициентов Фурье по N точкам требует O(N²) операций умножения. Алгоритм БПФ (FFT, Кули–Тьюки, 1965) сокращает это до O(N log N) — для N = 10⁶ это разница между 10¹² и 2·10⁷ операциями. Идея: разбиваем сигнал на чётные и нечётные отсчёты, вычисляем ДПФ рекурсивно. Применения: обработка звука и изображений (JPEG использует дискретное косинусное преобразование, родственное ДПФ), умножение больших чисел, расчёт свёрток в нейронных сетях. БПФ — один из наиболее влиятельных численных алгоритмов XX века по числу реальных применений.