Дифференциальные уравнения
Преобразование Фурье для ОДУ
Гренобль, 1807 год. Жозеф Фурье изучает остывание раскалённого железного стержня и заявляет дерзкое: любую функцию можно разложить в сумму синусов. Лагранж и Лежандр в Академии говорят 'нет'. Через 15 лет «Théorie analytique de la chaleur» (1822) переворачивает анализ. Сегодня тот же приём превращает ОДУ y'' + ay' + by = f(t) в алгебру: каждая гармоника входа независимо умножается на H(iω). А MP3, JPEG, MRI и Stable Diffusion - дальние потомки той самой идеи про горячий стержень.
- **Фильтрация сигналов:** умножить спектр на H(ω) - нежелательные частоты исчезнут. Так работают Bose QC45 (active noise cancellation) и эквалайзеры Logic Pro.
- **ОДУ с периодическим воздействием:** генератор переменного тока 50 Гц + RLC-цепь → отклик считается через H(iω) ровно за одну строчку.
- **Gaussian blur:** размытие изображения = решение уравнения теплопроводности = свёртка с гауссовым ядром = умножение спектра. Photoshop делает это за O(N log N) через FFT.
Предварительные знания
Ряды Фурье
Периодическая функция f(t) с периодом T разлагается в **ряд Фурье**: f(t) = a₀/2 + Σ[aₙcos(2πnt/T) + bₙsin(2πnt/T)]. Коэффициенты: **aₙ = (2/T)∫₀ᵀ f(t)cos(2πnt/T)dt**, аналогично bₙ. Ключевая идея - синусы и косинусы ортогональны на периоде.
**Ортогональность:** ∫₀ᵀ cos(2πmt/T)cos(2πnt/T)dt = 0 при m≠n. Это позволяет извлечь каждый коэффициент независимо - умножить f(t) на нужную функцию и проинтегрировать.
**Явление Гиббса:** вблизи разрывов ряд Фурье даёт выброс ~9% от величины скачка. Он не исчезает при увеличении числа гармоник - только сужается.
Жан-Батист Жозеф Фурье и горячий стержень (1807-1822)
В 1807 Фурье подал в Парижскую академию мемуар о распространении тепла в твёрдых телах. Лаплас, Лагранж и Лежандр отказали - идея, что любую функцию можно разложить в сумму синусов, казалась дикой (особенно Лагранжу). Фурье потратил 15 лет на полировку и в 1822 выпустил «Théorie analytique de la chaleur». Сегодня преобразование Фурье - инструмент №1 в обработке сигналов, MP3, JPEG, MRI, квантовой механике и решении PDE.
Без Фурье не было бы ни цифрового аудио, ни томографии, ни Stable Diffusion - все они опираются на разложение по гармоникам.
f(t) = t на [-π, π] - нечётная функция. Чему равны коэффициенты aₙ (при косинусах)?
Преобразование Фурье
При T → ∞ ряд Фурье переходит в **интеграл Фурье**: **F{f}(ω) = ∫_{-∞}^{∞} f(t)e^{-iωt}dt**. Обратное: f(t) = (1/2π)∫F(ω)e^{iωt}dω. F(ω) - частотный спектр: насколько сильно присутствует частота ω.
| f(t) | F{f}(ω) | Особенность |
|---|---|---|
| e^(-a|t|), a>0 | 2a/(a²+ω²) | Лоренцева линия |
| rect(t/T) | T sinc(ωT/2) | sinc = sin(x)/x |
| e^(-t²/2σ²) | σ√(2π) e^(-ω²σ²/2) | Гауссиан → гауссиан |
| δ(t) | 1 | Все частоты равны |
| cos(ω₀t) | π[δ(ω-ω₀)+δ(ω+ω₀)] | Две дельта-функции |
**Принцип неопределённости:** Δt · Δω ≥ 1/2. Узкий импульс → широкий спектр. Чистый тон → узкая спектральная линия, но длится бесконечно.
При сжатии гауссиана по времени (σ→0) что происходит с его фурье-спектром?
Теорема о свёртке и ОДУ
**Теорема о свёртке:** F{f*g}(ω) = F(ω)·G(ω). Свёртка во времени = умножение спектров. Для ОДУ y'' + ay' + by = f(t) с периодическим f каждая гармоника Fₙe^(iωₙt) усиливается в **Yₙ = Fₙ · H(iωₙ)**, где **H(s) = 1/(s²+as+b)** - передаточная функция.
**scipy.fft.fft** - быстрое преобразование Фурье (O(N log N)). Для фильтрации: умножить спектр на H(ω) и применить **ifft**. Для проектирования фильтров: **scipy.signal.butter**, **scipy.signal.cheby1**.
y'' + 4y = cos(2t). Что происходит с амплитудой?
Явление Гиббса
В точке разрыва ряд Фурье сходится к среднему значению. Вблизи разрыва - **выброс ~8.9%** от величины скачка независимо от числа гармоник. Выброс не исчезает - лишь сужается.
**Устранение Гиббса:** умножить коэффициенты на ядро Фейера или Ланцоша. Оконные функции (Hann, Blackman) - тот же принцип в обработке сигналов.
| Условие на f | Тип сходимости |
|---|---|
| f непрерывна и периодична | Равномерная |
| f кусочно-непрерывна | Поточечная + Гиббс на разрывах |
| f ∈ L²[0,T] | В норме L² |
| f абсолютно непрерывна | Равномерная, коэффициенты → 0 |
Ряд Фурье f(t) = |t| на [-π, π]. Как он сходится?
Ключевые идеи
- **Ряд Фурье:** f(t) = a₀/2 + Σ[aₙcos + bₙsin]. Коэффициенты через интеграл (ортогональность).
- **Преобразование Фурье:** F{f}(ω) = ∫f(t)e^{-iωt}dt. Спектр - частотный портрет сигнала.
- **Теорема о свёртке:** свёртка = умножение спектров. Фильтрация = умножение на H(ω).
- **Явление Гиббса:** ~9% выброс на разрывах. Не исчезает при добавлении гармоник.
Связанные темы
Фурье-методы связывают ОДУ с обработкой сигналов и уравнениями в частных производных:
- Преобразование Лапласа — Лаплас - обобщение Фурье на комплексные s. При s=iω получаем частотную характеристику H(iω) из этого урока.
- Уравнение теплопроводности — Решение ∂u/∂t = α∇²u через ряд Фурье - это разложение начального условия по собственным гармоникам sin(nπx/L) (de-07).
- Волновое уравнение — Стоячие волны uₙ = sin(nπx/L)·cos(ωₙt) - те же фурье-гармоники, только во времени они колеблются, а не затухают.
- Численные методы для ОДУ — FFT-солверы для PDE на периодической области (спектральные методы) дают экспоненциальную сходимость - быстрее любых FD/FEM.
Вопросы для размышления
- Принцип неопределённости Δt·Δω ≥ 1/2 ограничивает разрешение по времени и частоте. Как это связано с вейвлет-анализом?
- Явление Гиббса не исчезает при N→∞. Как с ним справляются в MP3 и JPEG на практике?
- DFT предполагает периодичность сигнала. Что происходит при анализе непериодических данных и как помогают оконные функции?