Тригонометрия

Ряды Фурье

Цели урока

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

Предварительные знания

Тригонометрия
Интегрирование
Линейные пространства

Тригонометрические функции

Как сложить бесконечно много синусов и получить прямоугольник - и почему это работает в JPEG?

JPEG: DCT сжимает 10 МП фото в 500 КБ без видимых потерь
Shazam: акустический отпечаток хранится как спектр Фурье
MP3: удаление частот вне диапазона слышимости по маскировке
МРТ: реконструкция изображений из k-пространства через обратное ФП

История: от теории тепла до цифровой эпохи

В 1807 году Фурье представил в Парижской академии трактат о теплопроводности - и его чуть не отклонили из-за идеи разложить функцию в тригонометрический ряд. Лагранж был скептиком. В 1822 году вышла «Аналитическая теория тепла» - и математика навсегда изменилась. Риман в 1854-м формализовал условия сходимости. Лебег в 1902-м решил проблему через интеграл. Кули и Тьюки в 1965-м дали алгоритм FFT за O(N log N) - и ряды Фурье стали практическим инструментом цифровой эры.

Ряды Фурье: разложение и коэффициенты

Жозеф Фурье в 1822 году показал, что любую периодическую функцию можно представить как бесконечную сумму синусов и косинусов. Это не заметим, что: складывать гладкие функции и получать функцию с разрывами. Shazam распознаёт песню за 5 секунд именно потому, что акустический отпечаток хранится как спектр Фурье - набор амплитуд на разных частотах.

Физический смысл: a_0/2 - среднее значение, a_n и b_n - амплитуды n-й гармоники. Чем гладче функция, тем быстрее убывают коэффициенты.

Как вычисляется n-й коэффициент Фурье a_n функции f?

Ортогональность базиса {cos(nx)} позволяет выделить каждый коэффициент независимо. Интеграл ∫cos(nx)cos(mx)dx = π·δ(n,m), что обнуляет все лишние слагаемые.

Сходимость и явление Гиббса

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

Прямоугольная волна и явление Гиббса

Классический пример: f(x) = sign(sin(x))

Ряд Фурье: f(x) = (4/π)(sin x + sin(3x)/3 + sin(5x)/5 + ...). При N = 100 слагаемых пик около разрыва всё равно превышает 1 примерно на 0.09 - независимо от N.

К чему сходится ряд Фурье функции f в точке разрыва x₀?

Теорема Дирихле-Жордана: в точке разрыва ряд Фурье сходится к (f(x₀⁺) + f(x₀⁻))/2. Вблизи разрыва возникает явление Гиббса: переколебание составляет примерно 8.9% от величины скачка.

Тождество Парсеваля и применения

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

Шазам хранит спектрограммы как разреженный набор пиков Фурье-спектра. Поиск совпадения - это сравнение двух разреженных векторов, а не полных аудиопотоков.

Что гарантирует тождество Парсеваля при усечении ряда Фурье до N членов?

Тождество Парсеваля - это изометрия L² ↔ l². Ошибка усечения ||f - S_N||² = Σ_{|n|>N} |c_n|² - точная формула, не оценка.

Связи с другими темами

Ряды Фурье объединяют анализ, линейную алгебру и теорию сигналов

DCT и JPEG — Связанная тема
Вейвлет-анализ — Связанная тема
Теория операторов — Связанная тема
Цифровая обработка сигналов — Связанная тема

Итоги

Ряд Фурье разлагает периодическую функцию в бесконечную сумму гармоник с коэффициентами через ортогональные проекции
В точке разрыва ряд сходится к среднему (f(x₀⁺) + f(x₀⁻))/2; вблизи разрыва возникает переколебание 8.9% - явление Гиббса
Тождество Парсеваля: энергия сохраняется при переходе в частотную область; ошибка усечения = сумма квадратов отброшенных коэффициентов
Гладкость определяет скорость убывания: C^k даёт O(n^{-k}), аналитическая функция - экспоненциальное убывание

Вопросы для размышления

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

Связанные уроки

trig-25-dct — DCT - дискретный аналог рядов Фурье в JPEG
trig-24-wavelets — Вейвлеты обобщают Фурье-анализ на нестационарные сигналы
trig-26-trig-poly — Тригонометрические полиномы - конечные частичные суммы рядов Фурье
trig-27 — Абстрактный Фурье-анализ на группах обобщает классический случай
la-01-vectors-intro

Цели урока

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

Понимать явление Гиббса и его причины

Применять тождество Парсеваля для оценки ошибки сжатия

Связывать гладкость функции со скоростью убывания коэффициентов

Ряды Фурье: разложение и коэффициенты

Как вычисляется n-й коэффициент Фурье a_n функции f?

Сходимость и явление Гиббса

Прямоугольная волна и явление Гиббса

Классический пример: f(x) = sign(sin(x))

К чему сходится ряд Фурье функции f в точке разрыва x₀?

Тождество Парсеваля и применения

Что гарантирует тождество Парсеваля при усечении ряда Фурье до N членов?

Тождество Парсеваля - это изометрия L² ↔ l². Ошибка усечения ||f - S_N||² = Σ_{|n|>N} |c_n|² - точная формула, не оценка.

Итоги

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

В точке разрыва ряд сходится к среднему (f(x₀⁺) + f(x₀⁻))/2; вблизи разрыва возникает переколебание 8.9% - явление Гиббса

Тождество Парсеваля: энергия сохраняется при переходе в частотную область; ошибка усечения = сумма квадратов отброшенных коэффициентов

Гладкость определяет скорость убывания: C^k даёт O(n^{-k}), аналитическая функция - экспоненциальное убывание

Связанные уроки

trig-25-dct — DCT - дискретный аналог рядов Фурье в JPEG

trig-24-wavelets — Вейвлеты обобщают Фурье-анализ на нестационарные сигналы

trig-26-trig-poly — Тригонометрические полиномы - конечные частичные суммы рядов Фурье

trig-27 — Абстрактный Фурье-анализ на группах обобщает классический случай

la-01-vectors-intro