Тригонометрия
Тригонометрические ряды и Фурье
Почему JPEG-фотография в 20 раз меньше, чем исходный файл? Как шумоподавление в наушниках Sony знает, какой звук заглушить? Как MRI строит изображения внутренних органов без рентгена? За всем стоит разложение Фурье.
- JPEG/MP3 сжатие: отбрасывание незначимых Фурье-коэффициентов без заметной потери качества
- Фильтрация ЭКГ: выделение сердечного ритма из зашумлённого сигнала
- МРТ: изображение строится из обратного преобразования Фурье k-пространства
Предварительные знания
Тригонометрические ряды: базис и ортогональность
MP3-сжатие использует 1152 сэмпла на фрейм и DCT (аналог ряда Фурье), чтобы сжать 44 100 Гц стерео аудио в 10 раз. Функции sin(nx) и cos(nx) образуют ортогональный базис пространства функций на [-π, π]. Любая «разумная» периодическая функция может быть представлена как бесконечная сумма синусов и косинусов - тригонометрический ряд Фурье.
**Тригонометрический ряд Фурье:** f(x) = a₀/2 + Σ_{n=1}^∞ [aₙ·cos(nx) + bₙ·sin(nx)] **Ортогональность** (ключевое свойство): ∫_{-π}^{π} cos(nx)·cos(mx) dx = π·δ_{nm} (0 при n≠m, π при n=m) ∫_{-π}^{π} sin(nx)·sin(mx) dx = π·δ_{nm} ∫_{-π}^{π} sin(nx)·cos(mx) dx = 0 (всегда!) **Комплексная форма:** f(x) = Σ_{n=-∞}^∞ cₙ·e^{inx}, cₙ = (aₙ - ibₙ)/2
Почему ∫_{-π}^{π} sin(nx)·cos(mx) dx = 0 для любых n, m?
Коэффициенты Фурье: проекция функции на базис
Коэффициенты Фурье вычисляются как скалярные произведения функции с базисными функциями. Это прямой аналог нахождения координат вектора: проекция на каждую ось базиса. Ортогональность делает вычисление каждого коэффициента независимым.
**Формулы для коэффициентов Фурье:** a₀ = (1/π) ∫_{-π}^{π} f(x) dx aₙ = (1/π) ∫_{-π}^{π} f(x)·cos(nx) dx bₙ = (1/π) ∫_{-π}^{π} f(x)·sin(nx) dx **Идея:** умножаем f(x) на нужную базисную функцию и интегрируем - все остальные члены «обнуляются» из-за ортогональности. **Равенство Парсеваля** (сохранение энергии): (1/π)∫|f|² dx = |a₀|²/2 + Σ(aₙ² + bₙ²)
Функция f(x) = x² на [-π, π]. Все ли коэффициенты bₙ равны нулю? Почему?
Сходимость: теорема Дирихле и феномен Гиббса
Не любой ряд Фурье сходится к исходной функции в каждой точке. Теорема Дирихле описывает условия сходимости. У точек разрыва функции ряд сходится к среднему - и при этом всегда возникает «выброс» около разрыва, называемый феноменом Гиббса.
**Теорема Дирихле:** Если f - кусочно-монотонная и кусочно-непрерывная на [-π, π], то ряд Фурье сходится: - К f(x) в точках непрерывности - К (f(x⁺) + f(x⁻))/2 в точках разрыва (среднее между левым и правым пределом) **Феномен Гиббса:** Около каждой точки разрыва частичная сумма превышает функцию примерно на ~9% (независимо от числа гармоник!). Конкретно: выброс → (2/π)∫₀^π sin(t)/t dt - 1 ≈ 0.0895 = 8.95% **Равномерная сходимость:** только при непрерывных функциях с кусочно-непрерывной производной.
Феномен Гиббса: при N → ∞ перебег ряда Фурье около точки разрыва стремится к...
Приложения: фильтрация, сжатие, анализ сигналов
Фурье-анализ - фундамент современной цифровой обработки сигналов. JPEG сжимает изображения через DCT (родственник Фурье). MP3 убирает неслышимые частоты. Фильтр нижних частот в частотной области - просто обнуление коэффициентов выше порогового значения.
**Свёртка = умножение в частотной области:** (f * g)(x) ↔ F̂(ω) · Ĝ(ω) **Фильтрация в частотной области:** 1. FFT(сигнал) → частотное представление 2. Умножить на маску фильтра H(ω) 3. IFFT → отфильтрованный сигнал **Фильтр нижних частот (low-pass):** H(ω) = 1 при |ω| ≤ ω_c, иначе 0 **Фильтр верхних частот (high-pass):** H(ω) = 1 при |ω| > ω_c, иначе 0 **JPEG** использует DCT (Discrete Cosine Transform) - вещественная версия DFT.
Почему фильтр нижних частот (low-pass) в частотной области эквивалентен свёртке со сглаживающим ядром в временной области?
Ключевые идеи
- Ортогональность: ∫sin(nx)cos(mx)dx = 0, ∫cos(nx)cos(mx)dx = π·δₙₘ - основа вычисления коэффициентов
- Теорема Дирихле: сходимость к среднему в точках разрыва
- Феномен Гиббса: перебег ~8.95% у разрывов не исчезает при N → ∞
- Фильтрация = умножение спектра × маска (теорема о свёртке)
Связанные темы
Ряды Фурье объединяют тригонометрию, анализ и инженерию:
- Комплексная тригонометрия — Формула Эйлера - основа комплексной формы рядов Фурье
- Тригонометрия в машинном обучении — Позиционные кодировки трансформеров - это Фурье-признаки
Вопросы для размышления
- Докажите равенство Парсеваля: как энергия сигнала во временной и частотной областях связаны? Используйте ортогональность базиса.
- Феномен Гиббса нельзя устранить увеличением числа гармоник. Как с ним борются на практике? Что такое оконные функции (Hann, Hamming)?
- Чему равен ряд Фурье функции f(x) = |sin x|? Вычислите несколько первых коэффициентов.