Функции как векторы
Ряд Фурье — это не просто вычислительный инструмент, а геометрическая операция: ортогональная проекция функции на базис.
Пространство функций
Множество квадратично интегрируемых функций на [0,1] образует векторное пространство L²[0,1]. Сложение & скалярное умножение работают поточечно. Скалярное произведение двух функций f, g:
⟨f, g⟩ = ∫₀¹ f(t) · g*(t) dt
где g* — комплексное сопряжение g. Это удовлетворяет всем аксиомам скалярного произведения.
Ортогональность базиса Фурье
Функции φ_k(t) = e^{i2πkt} образуют ортонормированный базис для L²[0,1]:
⟨φ_k, φ_m⟩ = ∫₀¹ e^{i2πkt} · e^{−i2πmt} dt = ∫₀¹ e^{i2π(k−m)t} dt = δ_{km}
(Это равно 1, если k = m, и 0 иначе — в результате интегрирования чистых колебаний за полный период.)
Коэффициент Фурье как скалярное произведение
k-й коэффициент Фурье функции x(t):
c_k = ⟨x, φ_k⟩ = ∫₀¹ x(t) · e^{−i2πkt} dt
Это проекция x на базисный вектор φ_k. Коэффициент измеряет, насколько φ_k присутствует в x.
Проекция на подпространство
Усечение ряда Фурье до 2N+1 членов проецирует x на подпространство, охватываемое {φ_{−N}, …, φ_N}. Усеченный ряд — это ортогональная проекция x на это конечномерное подпространство.
По неравенству Бесселя проекция минимизирует L² ошибку:
‖x − Σ_{k=−N}^{N} c_k φ_k‖² ≤ ‖x − Σ_{k=−N}^{N} a_k φ_k‖² для любого выбора a_k
Усечение Фурье — это наилучшее приближение в L² из этого подпространства. Оно минимизирует среднеквадратическую ошибку (квадрат нормы L² разности).
Прямоугольное окно → ядро синк
Прямоугольное окно во временной области (сохранение только коэффициентов для |k| ≤ N) соответствует умножению на функцию rect в пространстве индексов коэффициентов.
Умножение в одной области соответствует свертке в другой области.
Преобразование Фурье прямоугольного окна (в дискретном пространстве коэффициентов) — это ядро Дирихле (функция, похожая на синк):
D_N(f) = Σ_{k=−N}^{N} e^{i2πfk} = sin(π(2N+1)f) / sin(πf)
Когда мы усекаем ряд Фурье, мы свертываем идеальный отклик H_ideal(f) с D_N(f).
Почему происходит явление Гиббса
Ядро Дирихле имеет большие боковые лепестки, которые затухают медленно. Рядом с разрывом в H_ideal(f) эти боковые лепестки звенят — они складываются когерентно с одной стороны скачка, создавая ≈9% перескок.
Математическая константа: ∫₀^π sin(t)/t dt = Si(π) ≈ 1.8519. Величина перескока Гиббса = (2/π)·Si(π) − 1 ≈ 0.0895 = 8.95%. Это независимо от N.
Геометрия окна
Гладкое окно (Хэмминг, Ханн, Кайзер) имеет преобразование Фурье с более низкими боковыми лепестками. Свертка H_ideal(f) с ядром, имеющим более слабые боковые лепестки, дает меньше звенения. Компромисс: более низкие боковые лепестки всегда сопровождаются более широким главным лепестком, расширяя полосу перехода.
Константа Гиббса
Перескок Гиббса — это определенный интеграл, не зависящий от N.
Первый максимум N-членного частичного ряда Фурье единичного скачка находится примерно на расстоянии f ≈ 1/(2N) от разрыва. При N → ∞ этот максимум приближается к 1/(2)·(2/π)·Si(π) ≈ 1.0895.
Перескок: 0.0895 или примерно 8.95% высоты скачка.
Окна как частотные ядра
Каждая функция окна имеет преобразование Фурье, которое описывает ядро, используемое для сглаживания идеальной частотной характеристики.
Ключевые геометрические параметры ядра:
1. Ширина главного лепестка: определяет ширину переходной полосы (более широкий главный лепесток → более широкая переходная полоса).
2. Уровень первого бокового лепестка: определяет пульсацию полосы пропускания & полосы затухания (более низкие боковые лепестки → меньше пульсаций).
Эти два параметра не независимы. Для данной длины окна 2N+1, снижение высоты бокового лепестка требует расширения главного лепестка — всегда.
Окно Кайзера дает пользователю один параметр (α) для непрерывной компромиссной оптимизации высоты бокового лепестка против ширины главного лепестка, а не прыгание между фиксированными типами окон.
Идея проектирования
Ширина переходной полосы ΔF ≈ ширина главного лепестка / N. Пульсация δ ≈ уровень бокового лепестка. Обе формулы приблизительны; уравнения Кайзера делают их точными.