Интегральные преобразования и спектральные методы

Зачем переходить в пространство частот?

Дифференциальные уравнения в «x-t» пространстве содержат производные — это непросто. В «k-ω» пространстве (пространстве частот и волновых чисел) производная ∂/∂t заменяется на умножение на −iω, а ∂/∂x на ik. Дифференциальное уравнение превращается в алгебраическое! Это магия интегральных преобразований.

Дополнительный бонус: многие физические явления прозрачнее в пространстве частот. Период осцилляций, полоса пропускания, дисперсия — всё это свойства частотного спектра.

Преобразование Фурье

Прямое: F̂(k) = ∫_{−∞}^{+∞} f(x) e^{−ikx} dx

Обратное: f(x) = (1/2π) ∫_{−∞}^{+∞} F̂(k) e^{ikx} dk

Ключевые свойства: производная → умножение (f'(x) ↔ ikF̂(k)); свёртка f*g → произведение F̂·Ĝ; теорема Парсеваля: ∫|f(x)|²dx = (1/2π)∫|F̂(k)|²dk (сохранение нормы).

Числовой пример: волновое уравнение. ∂²u/∂t² = c² ∂²u/∂x². Фурье по x: ∂²û(k,t)/∂t² = −c²k²û. Это обычное ОДУ по t! Решение: û(k,t) = A(k)e^{ickt} + B(k)e^{−ickt}. Обратное Фурье даёт формулу Даламбера: u(x,t) = f(x+ct) + g(x−ct) — волны, бегущие в обе стороны.

Уравнение теплопроводности. ∂u/∂t = κ ∂²u/∂x². Фурье по x: ∂û/∂t = −κk²û. Решение: û(k,t) = û(k,0)e^{−κk²t}. При t > 0: высокие частоты затухают быстрее → размывание профиля. Обратное Фурье: u(x,t) = (1/√(4πκt)) ∫ u(x',0) e^{−(x−x')²/(4κt)} dx' — функция Грина теплопроводности в действии.

Преобразование Лапласа

Определение: F(s) = L{f(t)} = ∫₀^{+∞} f(t) e^{−st} dt, s = σ + iω ∈ ℂ

Таблица ключевых пар: L{1} = 1/s, L{e^{at}} = 1/(s−a), L{cos(ωt)} = s/(s²+ω²), L{f'(t)} = sF(s) − f(0).

Числовой пример: задача Коши для ОДУ. y'' + 5y' + 6y = 0, y(0) = 1, y'(0) = 0.

Шаг 1. Применяем L: (s²Y − s − 0) + 5(sY − 1) + 6Y = 0

Шаг 2. Y(s²+ 5s + 6) = s + 5 → Y = (s+5)/((s+2)(s+3))

Шаг 3. Разложение: Y = A/(s+2) + B/(s+3). A = (2+5)/(2−3+5)|{s=−2} = 3/(−1+3) → A = 3, B = −2. (Вычисляем: A = lim{s→−2}(s+2)Y = (−2+5)/(−2+3) = 3, B = lim_{s→−3}(s+3)Y = (−3+5)/(−3+2) = −2.)

Шаг 4. y(t) = 3e^{−2t} − 2e^{−3t}. Проверка: y(0) = 3−2 = 1 ✓, y'(0) = −6+6 = 0 ✓.

Спектральные методы для численного решения ДУ

Идея: разложить неизвестную функцию u(x,t) по базису:

u(x,t) ≈ Σₙ aₙ(t) φₙ(x)

Подстановкой в ДУ получаем систему ОДУ для коэффициентов aₙ(t). Базисы: {e^{ikx/L}} — Фурье (для периодических задач); многочлены Чебышёва {Tₙ(x)} — для непериодических на отрезке; полиномы Лежандра {Pₙ(x)} — для физических задач со сферической симметрией.

Преимущество спектральных методов: Для гладких функций сходимость экспоненциальная: ошибка ~ e^{−αN} при N членах ряда. Конечноразностные схемы дают лишь алгебраическую сходимость ~ N^{−p}. Поэтому спектральные методы применяют в атмосферном моделировании, вычислительной гидродинамике (коды ORCL, DEDALE).

Быстрое преобразование Фурье (FFT): алгоритм Кули–Тьюки вычисляет N-точечное ДПФ за O(N log N) операций вместо O(N²). При N = 10⁶: экономия в 50 раз. FFT — основа цифровой обработки сигналов (MP3, JPEG, связь 4G/5G).

Разложение в ряды Фурье: физический пример

Стоячая волна в замкнутом резонаторе: u(x,t) на [0,L], u(0) = u(L) = 0. Разложение: u(x,t) = Σₙ aₙ(t) sin(nπx/L). Каждая мода sin(nπx/L) — «нормальная мода». Её частота: ωₙ = nπc/L. При n=1: основная мода. При n=2: первый обертон с λ/2 = L.

Применение: оптические резонаторы лазеров; акустика музыкальных инструментов (флейта, орган); микроволновые резонаторы в ускорителях.

Реальные приложения

Wi-Fi (OFDM): сигнал разбивается на ~50 поднесущих частот (мод). Каждая поднесущая кодирует данные независимо — через IFFT/FFT передатчик/приёмник переходят между временной и частотной областями. Спектральная эффективность OFDM — прямое следствие ортогональности тригонометрического базиса.

Сейсмология: Фурье-анализ сигналов от землетрясений позволяет разделить P- и S-волны, определить частотный состав и оценить характеристики очага.

Интегральные преобразования в цифровой обработке сигналов

Преобразование Фурье и его обобщения — фундамент цифровых технологий. Быстрое преобразование Фурье (БПФ, алгоритм Кули–Тьюки, 1965) сократило вычислительную сложность ДПФ с O(N²) до O(N log N) и произвело революцию в цифровой обработке сигналов. БПФ используется в: аудиокодеках (MP3, AAC — вычисление модифицированного ДКП для частотного разложения), OFDM-модуляции (LTE, Wi-Fi, 5G — каждый поднесущий канал — это частотная составляющая из БПФ), радаре (дальность цели — через обратное БПФ свёртки отражённого и опорного сигнала), медицинской томографии (реконструкция КТ-изображений — через 2D БПФ теоремы центрального среза). Преобразование Лапласа — инструмент анализа линейных динамических систем в управлении: передаточные функции, критерий Найквиста устойчивости и синтез ПИД-регуляторов формулируются в s-области. BIBO-устойчивость фильтра эквивалентна расположению всех полюсов Н(s) в левой полуплоскости — прямое следствие теории комплексного анализа. В вычислительной химии быстрые алгоритмы для потенциалов взаимодействия (PME — Particle Mesh Ewald) используют 3D БПФ для ускорения молекулярной динамики.

Спектральные методы применяются в вычислительной гидродинамике (псевдоспектральные схемы для уравнений Навье–Стокса), где точность экспоненциально растёт с числом мод для гладких решений. Алгоритм быстрого синуса/косинуса (DCT/DST) ускоряет разложение Фурье на граничных областях и используется в пакетах LAPACK и FFTW, на которых основана вся научная вычислительная экосистема — от климатических моделей до расчёта турбулентности в авиастроении.

Задание: (а) Решите ∂u/∂t = κ ∂²u/∂x² на [0,L], u(0)=u(L)=0, u(x,0)=sin(πx/L). Найдите u(x,t) через разложение Фурье. Как зависит время «остывания» от L и κ? (б) Преобразованием Лапласа: y'' − y = e^{−t}, y(0)=y'(0)=0. (в) Оцените, при каком N = числе мод Фурье-ряд для функции f(x) = x(L−x) аппроксимирует её с точностью 1% на [0,L].